x*********s 发帖数: 5554 | 1 现在有一个30x30的矩阵A, 想用matlab简单求逆
rank(A)=30没问题,但是cond(A)=1.0058e+14,det(A)=-5.5468e-141
是一个病态矩阵,
matlab有没有简单的办法求的一个稳定的逆呀?
或者有什么简单的算法没有?
我数学基础太差,
仅有的几个包子求助,谢谢了 |
w***n 发帖数: 1084 | 2 既然只有30-by-30,你可以先算个SVD。
【在 x*********s 的大作中提到】 : 现在有一个30x30的矩阵A, 想用matlab简单求逆 : rank(A)=30没问题,但是cond(A)=1.0058e+14,det(A)=-5.5468e-141 : 是一个病态矩阵, : matlab有没有简单的办法求的一个稳定的逆呀? : 或者有什么简单的算法没有? : 我数学基础太差, : 仅有的几个包子求助,谢谢了
|
x*********s 发帖数: 5554 | 3 虽然还是每太弄清楚,不过还是谢谢
【在 w***n 的大作中提到】 : 既然只有30-by-30,你可以先算个SVD。
|
w**a 发帖数: 4743 | 4 [u,s,v]=svd(A)
A^{-1} = v*s^(-1)*u'
【在 x*********s 的大作中提到】 : 虽然还是每太弄清楚,不过还是谢谢
|
x*********s 发帖数: 5554 | 5 谢谢
【在 w**a 的大作中提到】 : [u,s,v]=svd(A) : A^{-1} = v*s^(-1)*u'
|
d********f 发帖数: 43471 | 6 svd truncation
【在 x*********s 的大作中提到】 : 现在有一个30x30的矩阵A, 想用matlab简单求逆 : rank(A)=30没问题,但是cond(A)=1.0058e+14,det(A)=-5.5468e-141 : 是一个病态矩阵, : matlab有没有简单的办法求的一个稳定的逆呀? : 或者有什么简单的算法没有? : 我数学基础太差, : 仅有的几个包子求助,谢谢了
|
w**a 发帖数: 4743 | 7 有包子,多谢拉
【在 x*********s 的大作中提到】 : 谢谢
|