24小时热门版块排行榜    

查看: 7711  |  回复: 42
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

hk_kong36

金虫 (著名写手)

[求助] 大型矩阵求逆该怎么做

如题,
是576*576的矩阵,直接用matlab的函数inv求解好像结果不太对,有下面的显示:
Matrix is close to singular or badly scaled.
         Results may be inaccurate. RCOND = 2.372912e-021.

matlab矩阵求逆还有没有别的尝试方法。。。。。。,还是我需要对矩阵做一些处理呢?

thanks
回复此楼
thanks,聆听花开的声音,需要的是沉淀
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主

【答案】应助回帖

引用回帖:
35楼: Originally posted by hk_kong36 at 2013-07-20 10:40:47
不管用噻。。。。。  正在补充知识ing...

手里没数据。不好弄。写了段代码,将矩阵中的所有元素的数值都归到1~2之间。通过求这个矩阵的逆进而求得原来矩阵的逆。
这样做的目的是为了排除计算过程中出现数值过大,或者过小,导致无法计算下去,但是如果矩阵是奇异的,仍然无能为力。
求逆函数如下
function Ani=qiuni(A)
nx = size(A,1);
A1(1:nx,1:nx) = 1;
I = eye(nx,nx);
amin = min(min(A));
a = amin(1);
amax = max(max(A));
b = amax(1);
delt = b-a;
Abar = (A-a)/delt+1;
Abarni = Abar\I;
temp = I-(1-a/delt)*Abarni*A1;
tempni = temp\I;
Ani = tempni*Abarni/delt;
每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。
38楼2013-07-23 16:09:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 43 个回答

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

感谢参与,应助指数 +1
是稀疏矩阵吗?阶数太高,用直接法有点问题,试着用迭代法。还有就是,非要matlab做吗?

» 本帖已获得的红花(最新10朵)

2楼2013-07-10 11:50:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

刚才没太看清楚,你是要求逆矩阵还是解方程啊?求逆矩阵就很难了。No offense,你要逆矩阵干什么呢?解是不能直接由逆矩阵直接导出的,那样很不准确。
3楼2013-07-10 12:13:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyqiao

金虫 (小有名气)

这是很麻烦的事情。

宋端磊;鹿涛;乔园园;
大型矩阵求逆的实现与应用[A];
中国化学会第26届学术年会化学信息学与化学计量学分会场论文集[C];
2008年
博客:blog.sina.com.cn/nkuqiao;应用化学:blog.sina.com.cn/nkapc;化学工程:blog.sina.com.cn/nkceg
4楼2013-07-10 12:50:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见