level 5
古乱今◆
楼主
我用的do循环解线性方程组 循环一次就解一次方程 再把相应的x,y赋值到向量里面 最后连线画图
求求哪位大佬抽空看看小弟这个![[小乖]](/static/emoticons/u5c0fu4e56.png)
代码如下
Clear[lbk]"@"[rbk]Subscript[lbk]G, 1[rbk] := 1; Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk] =Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk] = 0; \[lbk]Gamma[rbk]2 := 1;\[lbk]Gamma[rbk]1 := 1; K := 400; X = { };Do[lbk]A := ({ {-2 \[lbk]Gamma[rbk]1, -I*Subscript[lbk]G, 1[rbk], 0, I*Subscript[lbk]G, 1[rbk], 0, 0, 0, 0, 0}, {-I*Subscript[lbk]G, 1[rbk], -\[lbk]Gamma[rbk]1 - \[lbk]Gamma[rbk]2 - I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk], -I*Subscript[lbk]G, 2[rbk], 0, I*Subscript[lbk]G, 1[rbk], 0, 0, 0, 0}, {0, -I*Subscript[lbk]G, 2[rbk], -\[lbk]Gamma[rbk]1 - I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk] - I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk], 0, 0, I*Subscript[lbk]G, 1[rbk], 0, 0, 0}, {I*Subscript[lbk]G, 1[rbk], 0, 0, -\[lbk]Gamma[rbk]1 - \[lbk]Gamma[rbk]2 + I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk], -I* Subscript[lbk]G, 1[rbk], 0, I*Subscript[lbk]G, 2[rbk], 0, 0}, {2 \[lbk]Gamma[rbk]1, I*Subscript[lbk]G, 1[rbk], 0, -I*Subscript[lbk]G, 1[rbk], -2 \[lbk]Gamma[rbk]2, -I*Subscript[lbk]G, 2[rbk], 0, I*Subscript[lbk]G, 2[rbk], 0}, {0, 0, I*Subscript[lbk]G, 1[rbk], 0, -I*Subscript[lbk]G, 2[rbk], -\[lbk]Gamma[rbk]2 - I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk], 0, 0, I*Subscript[lbk]G, 2[rbk]}, {0, 0, 0, I*Subscript[lbk]G, 2[rbk], 0, 0, -\[lbk]Gamma[rbk]1 + I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk] + I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk], -I*Subscript[lbk]G, 2[rbk], 0}, {0, 0, 0, 0, I*Subscript[lbk]G, 2[rbk], 0, -I*Subscript[lbk]G, 1[rbk], -\[lbk]Gamma[rbk]2 + I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk], -I*Subscript[lbk]G, 2[rbk]}, {0, 0, 0, 0, 2 \[lbk]Gamma[rbk]2, I*Subscript[lbk]G, 2[rbk], 0, -I*Subscript[lbk]G, 2[rbk], 0} });B = {1, 0, 0, 0, 0, 0, 0, 0, 0};\[lbk]Rho[rbk] = LinearSolve[lbk]A, B[rbk];Subscript[lbk]\[lbk]Rho[rbk], 23[rbk] = Simplify[lbk]\[lbk]Rho[rbk][lbk][lbk]6[rbk][rbk][rbk];y = Abs[lbk]Subscript[lbk]G, 2[rbk] - I* K *\[lbk]Gamma[rbk]2*Subscript[lbk]\[lbk]Rho[rbk], 23[rbk][rbk];AppendTo[lbk]X, {y, Subscript[lbk]G, 2[rbk]}[rbk],{Subscript[lbk]G, 2[rbk], 0,25, 0.1}[rbk]ListLinePlot[lbk]{X}, PlotRange -> All, PlotLegends -> {"y", "x"}[rbk]

2023年10月16日 03点10分
1
求求哪位大佬抽空看看小弟这个
Clear[lbk]"@"[rbk]Subscript[lbk]G, 1[rbk] := 1; Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk] =Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk] = 0; \[lbk]Gamma[rbk]2 := 1;\[lbk]Gamma[rbk]1 := 1; K := 400; X = { };Do[lbk]A := ({ {-2 \[lbk]Gamma[rbk]1, -I*Subscript[lbk]G, 1[rbk], 0, I*Subscript[lbk]G, 1[rbk], 0, 0, 0, 0, 0}, {-I*Subscript[lbk]G, 1[rbk], -\[lbk]Gamma[rbk]1 - \[lbk]Gamma[rbk]2 - I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk], -I*Subscript[lbk]G, 2[rbk], 0, I*Subscript[lbk]G, 1[rbk], 0, 0, 0, 0}, {0, -I*Subscript[lbk]G, 2[rbk], -\[lbk]Gamma[rbk]1 - I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk] - I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk], 0, 0, I*Subscript[lbk]G, 1[rbk], 0, 0, 0}, {I*Subscript[lbk]G, 1[rbk], 0, 0, -\[lbk]Gamma[rbk]1 - \[lbk]Gamma[rbk]2 + I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk], -I* Subscript[lbk]G, 1[rbk], 0, I*Subscript[lbk]G, 2[rbk], 0, 0}, {2 \[lbk]Gamma[rbk]1, I*Subscript[lbk]G, 1[rbk], 0, -I*Subscript[lbk]G, 1[rbk], -2 \[lbk]Gamma[rbk]2, -I*Subscript[lbk]G, 2[rbk], 0, I*Subscript[lbk]G, 2[rbk], 0}, {0, 0, I*Subscript[lbk]G, 1[rbk], 0, -I*Subscript[lbk]G, 2[rbk], -\[lbk]Gamma[rbk]2 - I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk], 0, 0, I*Subscript[lbk]G, 2[rbk]}, {0, 0, 0, I*Subscript[lbk]G, 2[rbk], 0, 0, -\[lbk]Gamma[rbk]1 + I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 1[rbk] + I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk], -I*Subscript[lbk]G, 2[rbk], 0}, {0, 0, 0, 0, I*Subscript[lbk]G, 2[rbk], 0, -I*Subscript[lbk]G, 1[rbk], -\[lbk]Gamma[rbk]2 + I*Subscript[lbk]\[lbk]CapitalDelta[rbk], 2[rbk], -I*Subscript[lbk]G, 2[rbk]}, {0, 0, 0, 0, 2 \[lbk]Gamma[rbk]2, I*Subscript[lbk]G, 2[rbk], 0, -I*Subscript[lbk]G, 2[rbk], 0} });B = {1, 0, 0, 0, 0, 0, 0, 0, 0};\[lbk]Rho[rbk] = LinearSolve[lbk]A, B[rbk];Subscript[lbk]\[lbk]Rho[rbk], 23[rbk] = Simplify[lbk]\[lbk]Rho[rbk][lbk][lbk]6[rbk][rbk][rbk];y = Abs[lbk]Subscript[lbk]G, 2[rbk] - I* K *\[lbk]Gamma[rbk]2*Subscript[lbk]\[lbk]Rho[rbk], 23[rbk][rbk];AppendTo[lbk]X, {y, Subscript[lbk]G, 2[rbk]}[rbk],{Subscript[lbk]G, 2[rbk], 0,25, 0.1}[rbk]ListLinePlot[lbk]{X}, PlotRange -> All, PlotLegends -> {"y", "x"}[rbk]
