有大佬知道奇异的矩阵怎么求逆么
mathematica吧
全部回复
仅看楼主
level 5
古乱今◆ 楼主
矩阵没有错误 还请各位老哥搭把手[小乖][小乖] 感激不尽
2023年10月31日 08点10分 1
level 5
古乱今◆ 楼主
Clear["@"]
Subscript[K, 1] = 1000;
Subscript[K, 2] = 100;
Subscript[\[Gamma], 1] = Subscript[\[Gamma], 2] = 1;
Subscript[\[CapitalDelta], 1] = Subscript[\[CapitalDelta], 2] = 0;
dd11 = {};
y1 = 5;
Do[
e = 10^-5;
x10 = {2 x2, 0.1, y1/2, 0.1, 0.4, 0.4, 0.2, 0.2, 0.1, 0.3, 0.2, 0.2,
0.1};
F1[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := -2 Subscript[\
\[Gamma], 1]*\[Rho]11 + 2 x1r*\[Rho]12i;
F2[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := -(Subscript[\
\[Gamma], 1] + Subscript[\[Gamma], 2])*\[Rho]12r +
x1r*(\[Rho]22 - \[Rho]11) +
Subscript[\[CapitalDelta], 1]*\[Rho]12i + x2*\[Rho]13i;
F3[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := -(Subscript[\
\[Gamma], 1] + Subscript[\[Gamma], 2])*\[Rho]12i +
x1i*(\[Rho]22 - \[Rho]11) - x2*\[Rho]13r -
Subscript[\[CapitalDelta], 1]*\[Rho]12r;
F4[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := (Subscript[\
\[CapitalDelta], 1] + Subscript[\[CapitalDelta], 2]) \[Rho]13i -
Subscript[\[Gamma], 1]*\[Rho]13r - x1r*\[Rho]23i - x1i*\[Rho]23r +
x2*\[Rho]12i;
F5[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := \
-Subscript[\[Gamma],
1]*\[Rho]13i - (Subscript[\[CapitalDelta], 1] +
Subscript[\[CapitalDelta], 2]) \[Rho]13r + x1r*\[Rho]23r -
x1i*\[Rho]23i - x2*\[Rho]12r;
F6[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := \
-Subscript[\[Gamma], 2]*\[Rho]23r +
Subscript[\[CapitalDelta], 2]*\[Rho]23i - x1r*\[Rho]13i -
x1i*\[Rho]13r;
2023年10月31日 08点10分 2
level 5
古乱今◆ 楼主
代码有点长[小乖]
2023年10月31日 08点10分 3
level 8
你的do循环没给全[小乖]另外你要找的是不是伪逆(PseudoInverse)
2023年10月31日 09点10分 4
求伪逆 可以解 但图像很奇怪[泪]
2023年10月31日 09点10分
谢谢大佬[真棒]
2023年10月31日 09点10分
level 5
奇异矩阵不是没有逆只有广义逆吗
2023年10月31日 09点10分 5
是啊 我都不知道是哪里错了 图像是散点
2023年10月31日 09点10分
吧务
level 15
……为什么你能肯定地说“矩阵没有错误”?还有
Clear["@"]
明显错了,这又是从哪种语言里想当然出来的?
2023年11月04日 05点11分 6
"@" 匹配非大写字符,会匹配到系统的保护符号,可以改为 Clear["`@"]
2023年11月04日 08点11分
@隨意超 ……这还真是被将了一军。长知识了。
2023年11月04日 10点11分
[小乖][小乖]有错误
2023年11月04日 11点11分
1