求助：利用fortran写一个矩阵求逆的程序

来源：蜘蛛抓取(WebSpider) 时间：2016-12-01 11:44 标签：
有偿求fortran矩阵求逆代码_fortran吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：4,380贴子：
有偿求fortran矩阵求逆代码收藏
哪位亲爱的大大，有fortran矩阵求逆的代码，能用有偿，原矩阵是32*32，部分数据如下：
苹果壁纸图片软件哪个好用？
我可以帮你写
可以帮忙，需要私信
给多少钱嘛
module pubimplicit double precision(A-H,O-Z)ALLOCATABLE :: A(:,:),B(:,:)end module pubprogram AAuse imsluse pubimplicit double precision(A-H,O-Z)open(12,file='A.txt') !!!A矩阵的TXT文件要复制到程序的目录下open(13,file='B.txt')ALLOCATE(A(32,2),B(32,2))read(12,*) ((A(i,j),i=1,32),j=1,32)B=.i.A !!!!B是A的逆矩阵wRITE (13,&(32(32(f18.12,a),/))&)((B(I,J),char(9),j=1,32),i=1,32) end program
登录百度帐号推荐应用
为兴趣而生，贴吧更懂你。或帐号:密码:下次自动登录{url:/nForum/slist.json?uid=guest&root=list-section}{url:/nForum/nlist.json?uid=guest&root=list-section}
贴数:4&分页:lily发信人: lilypz (lily), 信区: NumComp
标&&题: 请问谁有fortran的矩阵求逆的程序呢
发信站: 水木社区 (Thu Dec&&3 09:16:35 2009), 站内 && 如果有，能不能共享啊，谢谢
春天天气暖洋洋 && ※ 来源:·水木社区 ·[FROM: 218.94.142.*]
不要拿法律当挡箭牌发信人: figgers (Story of Wisteria Lane), 信区: NumComp
标&&题: Re: 请问谁有fortran的矩阵求逆的程序呢
发信站: 水木社区 (Mon Dec&&7 17:05:11 2009), 站内 && 网上一搜一大把
好不好用就难说了，解解超级奇异的矩阵试试。 && 论效率的话还是得用IMSL，MKL这些库
最近新用了MKL库的求逆程序，感觉比IMSL的快了不是一点半点 && 【 在 lilypz (lily) 的大作中提到: 】
: 如果有，能不能共享啊，谢谢
&&&& -- && ※ 来源:·水木社区 newsmth.net·[FROM: 211.99.222.*]
千与千寻发信人: thousand (壹仟), 信区: NumComp
标&&题: Re: 请问谁有fortran的矩阵求逆的程序呢
发信站: 水木社区 (Tue Dec&&8 04:52:47 2009), 站内 && 为啥 MKL不是IMSL的子集吗 && 【 在 figgers (Story of Wisteria Lane) 的大作中提到: 】
网上一搜一大把
好不好用就难说了，解解超级奇异的矩阵试试。 && 论效率的话还是得用IMSL，MKL这些库
最近新用了MKL库的求逆程序，感觉比IMSL的快了不是一点半点 && 【 在 lilypz (lily) 的大作中提到: 】
: 如果有，能不能共享啊，谢谢
&&&& -- &&&&&& -- && ※ 来源:·水木社区 newsmth.net·[FROM: 211.99.222.*]
建议刑法设立卢瑟罪和猥琐罪发信人: ruster (企鹅之王), 信区: NumComp
标&&题: Re: 请问谁有fortran的矩阵求逆的程序呢
发信站: 水木社区 (Tue Dec&&8 09:53:01 2009), 站内 && intel的程序民工写的和业余程序员写的差别。 && 【 在 thousand (壹仟) 的大作中提到: 】
: 为啥 MKL不是IMSL的子集吗
: 网上一搜一大把
: 好不好用就难说了，解解超级奇异的矩阵试试。
: ...................
7.辛亥到本世纪，中国科学界最伟大的工作是什么？
答：航天医学研究证明熟花生米可以发芽，为亡灵复活术奠定了坚实的基础。 &&&& ※ 来源:·水木社区 newsmth.net·[FROM: 202.120.58.*]
文章数:4&分页:匿名用户不能发表回复！|
每天回帖即可获得10分可用分！小技巧：
你还可以输入10000个字符
(Ctrl+Enter)
请遵守CSDN，不得违反国家法律法规。
转载文章请注明出自“CSDN（www.csdn.net）”。如是商业用途请联系原作者。>> Fortran矩阵求逆简单高效程序
Fortran矩阵求逆简单高效程序
所属分类：
下载地址：
矩阵求逆.zip文件大小：821.52 kB
分享有礼！ 》
请点击右侧的分享按钮，把本代码分享到各社交媒体。
通过您的分享链接访问Codeforge，每来2个新的IP，您将获得0.1 积分的奖励。
通过您的分享链接，每成功注册一个用户，该用户在Codeforge上所获得的每1个积分，您都将获得0.2 积分的分成奖励。
! &Description : 计算逆矩阵
!-----------------------------------------------------
! &Contains & &:
! & & &1. & driver &驱动函数
! & & &2. & solve & 方法函数&
已经经过验证，高效简洁
Sponsored links
源码文件列表
温馨提示： 点击源码文件名可预览文件内容哦 ^_^
3.67 kB 21:09
117.00 B 09:54
125.00 B 09:51
177.00 B 09:53
177.00 B 21:09
196.00 B 09:52
INV_MAT.MOD3.13 kB 21:09
SOURCE1.OBJ37.78 kB 21:09
??i???-?.EXE867.00 kB 21:09
??i???-?.EXE.INTERMEDIATE.MANIFEST381.00 B 21:09
??i???-?.PDB2.21 MB 21:09
SOURCE1.F905.71 kB 08:56
??i???-?.U2D320.00 B 21:14
??i???-?.VFPROJ1.93 kB 14:39
??i???-?.SLN902.00 B 13:57
??i???-?.SUO8.50 kB 21:14
（提交有效评论获得积分）
评论内容不能少于15个字，不要超出160个字。
yysdnn可以高效的求解普通矩阵，可是对于条件数比较大病态矩阵求解的还不准确
评价成功，多谢！
下载矩阵求逆.zip
CodeForge积分（原CF币）全新升级，功能更强大，使用更便捷，不仅可以用来下载海量源代码马上还可兑换精美小礼品了
您的积分不足，优惠套餐快速获取 30 积分
10积分 / ￥100
30积分 / ￥200原价 ￥300 元
100积分 / ￥500原价 ￥1000 元
订单支付完成后，积分将自动加入到您的账号。以下是优惠期的人民币价格，优惠期过后将恢复美元价格。
支付宝支付宝付款
微信钱包微信付款
更多付款方式：、
您本次下载所消耗的积分将转交上传作者。
同一源码，30天内重复下载，只扣除一次积分。
鲁ICP备号-3 runtime:Elapsed:88.433ms - init:0.2;find:2.9;t:1.9;tags:0.2;related:46.6;comment:2.1; 27.69
登录 CodeForge
还没有CodeForge账号？
Switch to the English version?
^_^"呃 ...
Sorry!这位大神很神秘，未开通博客呢，请浏览一下其他的吧求助：matlab矩阵求逆问题！ — 编程爱好者社区
主题：求助：matlab矩阵求逆问题！
[专家分：0]
请问有些矩阵为什么用matlab求得的逆矩阵和用vc自己编程求出的逆矩阵不同呢。vc程序分别采用了初等矩阵变换法、LU分解法、全选主元法，这三种方法求出的逆矩阵是完全相同的。对于一般的矩阵求出的逆矩阵和matlab的计算结果是一致的，但是当矩阵病态的时候求出的结果差别比较大。请问matlab中是不是对病态矩阵的求逆作了处理，vc程序的求逆怎么来解决呢？请问大家有没有什么好的解决方法？谢谢！
本帖地址：&
回复列表 （共6个回复）
[专家分：0]
[专家分：0]
inv（）就是求逆阵。是基本的函数。QQ
[专家分：40]
MATLAB对病态矩阵的求逆结果不是很好，但是在工程领域还可以凑合使用。VC、FORTRAN程序的求逆结果比较好的。
[专家分：0]
可以用&inv()哈
[专家分：370]
matlab的建议是不管矩阵是不是ill-conditioned，能避免求逆的尽量避免。比如Ax=B的时候用&A\B而不是x=inv（A)×B；前者即使在矩阵接近singular的时候也能得的较准确结果，而且速度快
[专家分：0]
楼主可以试试这个pinv()这个是pseudoinverse(就是求A+)对于A-存在的时候A+=A-所以可以避免病态的问题了
您尚未登录，请登录后再回复。

我要回帖

 

随机推荐