f*****y
发帖数: 1997 | 1 A是稀疏对称矩阵，diagonal dominant
结果用MATLAB 求 x = A\b
x里面有些元素是复数，这个是啥回事？ |
v*******e
发帖数: 11604 | 2
matlab忘记了，Ab是啥意思？A^(-1)*b看看
【在 f*****y 的大作中提到】
: A是稀疏对称矩阵，diagonal dominant
: 结果用MATLAB 求 x = A\b
: x里面有些元素是复数，这个是啥回事？
|
s****r
发帖数: 31686 | 3 啥会是？就是只有复数解
【在 f*****y 的大作中提到】
: A是稀疏对称矩阵，diagonal dominant
: 结果用MATLAB 求 x = A\b
: x里面有些元素是复数，这个是啥回事？
|
m**********e
发帖数: 12525 | 4 我肏,杀老鼠的?
实矩阵难道不能有复数特征值了??
【在 f*****y 的大作中提到】
: 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 的大作中提到】
: 我肏,杀老鼠的?
: 实矩阵难道不能有复数特征值了??
|
s****r
发帖数: 31686 | 7 属实，可以忽略
【在 f*****y 的大作中提到】
: 这个矩阵是准对称的， 并且是 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 的大作中提到】
: 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 的大作中提到】
: 这个矩阵是准对称的， 并且是 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 的大作中提到】
: 肏,纯数学上讲,矩阵特征值为实数,需要外加强制条件,比如 self-adjoint
|
w****y
发帖数: 2501 | 12 苦主都说对称了，怎么还能nonsuqare?
【在 k******y 的大作中提到】
: Check if A is a nonsquare matrix, probably MATLAB is being smart and
: calculate Moore-Penrose pseudoinverse of A
|
