大型矩阵的求逆－－5000*5000

penghust · 发表于 2008-10-28 13:37

呵呵，应该是pinv
仁兄好眼力

fourthfour · 发表于 2011-5-5 11:04

回复 10 # ChaChing 的帖子

请问下CFD中是怎么解决的？请指教下，看我能不能得到点想法

ChaChing · 发表于 2011-8-21 00:32

回复 17 # fourthfour 的帖子

对不住, 刚才发现
我知道的是, 猜一组解, 再循环求解直到收歛

fourthfour · 发表于 2011-8-28 21:34

本帖最后由 fourthfour 于 2011-8-28 21:34 编辑

回复 18 # ChaChing 的帖子

你的意思是说用迭代法求解啊，这样速度不知道怎么样，可是迭代法一般比较慢啊

ChaChing · 发表于 2011-8-28 21:57

回复 19 # fourthfour 的帖子

细节非个人专业, 玩CFD的或许清楚些

fourthfour · 发表于 2011-9-10 14:16

回复 20 # ChaChing 的帖子

哦，谢谢了啊

mxlzhenzhu · 发表于 2013-11-5 19:39

本帖最后由 mxlzhenzhu 于 2013-11-5 20:01 编辑

这个是老问题，而且解可能不唯一，
A_inverse=A\eye(size(A,1));

matlab 的pinv不一定适合我们的解。

pinv的算法试了一下，亏秩的时候，好像就是下面的svd代码：
[P,S,Q]=svd(A);
[~,ki]=find(diag(S)>eps);%% ki 就是非零元素对应的序号
d=diag(1./diag(S(ki,ki)));
A_inverse=Q(:,ki)*d*P(:,ki)';

ChaChing · 发表于 2013-11-5 21:37

不知道是否老了, 有些退化又懒了
总以为动不动大矩阵不是件好事!?
现今计算机速度/内存已非昔日可比, 那前人怎处理相关工程问题!?

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 大型矩阵的求逆－－5000*5000

回复 16楼 ChaChing 的帖子

评分