声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: penghust

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

[复制链接]
 楼主| 发表于 2008-10-28 13:37 | 显示全部楼层

回复 16楼 ChaChing 的帖子

呵呵,应该是pinv
仁兄好眼力
回复 支持 反对
分享到:

使用道具 举报

发表于 2011-5-5 11:04 | 显示全部楼层
回复 10 # ChaChing 的帖子

请问下CFD中是怎么解决的?请指教下,看我能不能得到点想法
发表于 2011-8-21 00:32 | 显示全部楼层
回复 17 # fourthfour 的帖子

对不住, 刚才发现
我知道的是, 猜一组解, 再循环求解直到收歛
发表于 2011-8-28 21:34 | 显示全部楼层
本帖最后由 fourthfour 于 2011-8-28 21:34 编辑

回复 18 # ChaChing 的帖子

你的意思是说用迭代法求解啊,这样速度不知道怎么样,可是迭代法一般比较慢啊
发表于 2011-8-28 21:57 | 显示全部楼层
回复 19 # fourthfour 的帖子

细节非个人专业, 玩CFD的或许清楚些
发表于 2011-9-10 14:16 | 显示全部楼层
回复 20 # ChaChing 的帖子

哦,谢谢了啊
发表于 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)';

评分

1

查看全部评分

发表于 2013-11-5 21:37 | 显示全部楼层
不知道是否老了, 有些退化又懒了
总以为动不动大矩阵不是件好事!?
现今计算机速度/内存已非昔日可比, 那前人怎处理相关工程问题!?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-24 07:46 , Processed in 0.137132 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表