f*****y 发帖数: 1997 | 1 A是稀疏对称矩阵,diagonal dominant
结果用MATLAB 求 x = A\b
x里面有些元素是复数,这个是啥回事? |
v*******e 发帖数: 11604 | 2
matlab忘记了,Ab是啥意思?A^(-1)*b看看
【在 f*****y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : A是稀疏对称矩阵,diagonal dominant : 结果用MATLAB 求 x = A\b : x里面有些元素是复数,这个是啥回事?
|
s****r 发帖数: 31686 | 3 啥会是?就是只有复数解
【在 f*****y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : A是稀疏对称矩阵,diagonal dominant : 结果用MATLAB 求 x = A\b : x里面有些元素是复数,这个是啥回事?
|
m**********e 发帖数: 12525 | 4 我肏,杀老鼠的?
实矩阵难道不能有复数特征值了??
【在 f*****y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : A是稀疏对称矩阵,diagonal dominant : 结果用MATLAB 求 x = A\b : x里面有些元素是复数,这个是啥回事?
|
m**********e 发帖数: 12525 | 5 肏,纯数学上讲,矩阵特征值为实数,需要外加强制条件,比如 self-adjoint |
f*****y 发帖数: 1997 | 6 这个矩阵是准对称的, 并且是 diagonally dominant
用 x = A \ b
算出的复数解, 虚数部分比实数部分要小的多,
比如:
这是X的前几个数:
1.508492629370379e+005 +1.241971211869262e-011i
1.315757251394335e+005 +5.302048302827487e-011i
1.086773104076401e+005 +2.957613529394406e-010i
8.242399458762619e+004 +1.776110658155558e-009i
5.262795581796929e+004 +1.127793683619662e-008i
3.599937155297921e+004 +4.027286494040506e-008i
2.792014354044897e+004 +9.034197192842009e-008i
这个虚数部分这么小,是不是数值误差引起的?
用实数部分的答案做图,发现还是符合物理意义的,
是不是不用理会虚数的部分了?
【在 m**********e 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 我肏,杀老鼠的? : 实矩阵难道不能有复数特征值了??
|
s****r 发帖数: 31686 | 7 属实,可以忽略
【在 f*****y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 这个矩阵是准对称的, 并且是 diagonally dominant : 用 x = A \ b : 算出的复数解, 虚数部分比实数部分要小的多, : 比如: : 这是X的前几个数: : 1.508492629370379e+005 +1.241971211869262e-011i : 1.315757251394335e+005 +5.302048302827487e-011i : 1.086773104076401e+005 +2.957613529394406e-010i : 8.242399458762619e+004 +1.776110658155558e-009i : 5.262795581796929e+004 +1.127793683619662e-008i
|
D*******o 发帖数: 3229 | 8 Try this:
x=inv(A)*b;
Please let me know the results. |
f*****y 发帖数: 1997 | 9 还是有虚数, 不过虚数部分仍然比实数部分小15个数量级以上。
难道是我的物理系统有问题?
【在 D*******o 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Try this: : x=inv(A)*b; : Please let me know the results.
|
k******y 发帖数: 1407 | 10 Check if A is a nonsquare matrix, probably MATLAB is being smart and
calculate Moore-Penrose pseudoinverse of A
【在 f*****y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 这个矩阵是准对称的, 并且是 diagonally dominant : 用 x = A \ b : 算出的复数解, 虚数部分比实数部分要小的多, : 比如: : 这是X的前几个数: : 1.508492629370379e+005 +1.241971211869262e-011i : 1.315757251394335e+005 +5.302048302827487e-011i : 1.086773104076401e+005 +2.957613529394406e-010i : 8.242399458762619e+004 +1.776110658155558e-009i : 5.262795581796929e+004 +1.127793683619662e-008i
|
w****y 发帖数: 2501 | 11 他的矩阵即便有复数,如果是Hermitian,特征值也是实数。
【在 m**********e 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 肏,纯数学上讲,矩阵特征值为实数,需要外加强制条件,比如 self-adjoint
|
w****y 发帖数: 2501 | 12 苦主都说对称了,怎么还能nonsuqare?
【在 k******y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : Check if A is a nonsquare matrix, probably MATLAB is being smart and : calculate Moore-Penrose pseudoinverse of A
|