请问大佬 解线性方程组无解
mathematica吧
全部回复
仅看楼主
level 5
古乱今◆ 楼主
我用的do循环解线性方程组 循环一次就解一次方程 再把相应的x,y赋值到向量里面 最后连线画图
求求哪位大佬抽空看看小弟这个[小乖][小乖]代码如下
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
level 8
试了一下,你的A矩阵秩均为8,增广矩阵秩均为9,所以确实无解[小乖]
2023年10月16日 08点10分 2
我把第一行只取3个1 算出来了[小乖]我也不知道为啥
2023年10月16日 08点10分
@古乱今◆ 总之系数矩阵和增广矩阵的秩要一样啦[微微一笑]应该是你取1了之后就一样了
2023年10月16日 09点10分
@qflwq 有道理[真棒]
2023年10月16日 10点10分
1