我这个感觉编出来不对劲，因为我们已经有图了画出来也不对
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]]