有大佬看看光学双稳态的程序吗
mathematica吧
全部回复
仅看楼主
level 1
二星rial✨ 楼主
我这个感觉编出来不对劲,因为我们已经有图了画出来也不对
Remove["Global`*"]Clear["Global`*"]\[Omega]d = 2*Pi*200*10^12; \[Kappa]1 = 2*Pi*10^9; \[Kappa]2 = 2*Pi*10^9; \[Kappa]1e = 2*Pi*0.7*10^9; \[Kappa]2e = 2*Pi*0.7*10^9; pin = 0.1/10^3; \[Kappa]0 = \[Kappa]1; \[CapitalDelta]0 = m*\[Kappa]0; \[CapitalDelta]2 = \[CapitalDelta]1 = \[CapitalDelta]0; wm = 2*Pi*6*10^9; J = 2*Pi*4*10^9; \[Gamma]m = 2*Pi*5*10^6; g = 1.6*Pi*10^6; \[CapitalDelta] = \[CapitalDelta]1 - (4* J^2*\[CapitalDelta]2)/(\[CapitalDelta]2^2 + 4*\[CapitalDelta]2^2); \[Kappa] = \[Kappa]1 + (4*\[Kappa]2* J^2)/(\[Kappa]2^2 + 4*\[CapitalDelta]2^2); U21 = ((\[Kappa]2^2 + 4*\[CapitalDelta]2^2)/(4*\[Kappa]2e*J^2))* U; \[Epsilon] = ((-2*I*J*Sqrt[\[Kappa]1e*\[Kappa]2e])/(\[Kappa]2 + 2*I*\[CapitalDelta]2))* sin; sin = Sqrt[ pin/((\[Omega]d*1.05457266)/10^34)]; \[Alpha]1n = sin; U = -(g^2/ wm); NSolve[(\[Kappa]^2/4 + \[CapitalDelta]^2)*\[Alpha] + 2*\[CapitalDelta]*U21*\[Alpha]^2 + U21^2*\[Alpha]^3 == ((20*J^2*\[Kappa]1e*\[Kappa]2e)/(\[Kappa]2^2 + 4*\[CapitalDelta]2^2))*sin, \[Alpha]]
2020年11月05日 06点11分 1
吧务
level 10
你是要别人看什么?
觉得不对劲,那你是不对劲在哪?
画的图不对,画图的代码在哪里,对错的标准是什么?
2020年11月06日 04点11分 2
算了打扰了,贵吧我真的惹不起
2020年11月06日 06点11分
@二星rial✨ ……你的问题确实很不清楚,现阶段恐怕只有非常了解“光学双稳态”的人才有可能回答你的问题——如果你真的指望着遇到这样的回答者,那我建议你看看这帖4楼:https://tieba.baidu.com/p/2964416898
2020年11月06日 16点11分
@xzcyr 没事了我自己做出来了,打扰了
2020年11月07日 00点11分
@xzcyr 我本来也想链这个4L的回复来着,不过当时我不记得是在哪了。
2020年11月07日 03点11分
1