| 查看: 1139 | 回复: 5 | |||
fagege2010铁虫 (小有名气)
|
[交流]
【求助】急求！！！
|
|
我在前面的子程序中计算出了一个矩阵AR（N,N） 现在想用ISML求它的本征值和本征矢 请问如何求；如何把子程序中的这个矩阵中的数据传过去啊？具体怎么编？ 问题已经解决谢谢各位 [ Last edited by fagege2010 on 2011-1-8 at 23:54 ] |
回复此楼» 猜你喜欢
什么时候开奖？
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急求RT-PCR电泳图分析
已经有22人回复
- 急求大家的意见！
已经有13人回复
- 急求帮助我翻译成中文
已经有1人回复
- 急求统计分析问题！
已经有10人回复
- 一段翻译 急求
已经有2人回复
- 急求翻译，谢谢
已经有1人回复
- 理解不透！急求帮忙！！！
已经有8人回复
- 【急求】正庚烷-苯二元体系已知折光率求组分
已经有3人回复
- 急求工艺流程图
已经有7人回复
- 急求！！！！！！翻译摘要，谢谢啦
已经有1人回复
- 急求翻译！！！！！汉译英！！！！在线等！！！！
已经有12人回复
- 急求英文翻译
已经有2人回复
- 急求汉译英翻译
已经有3人回复
- 急求
已经有8人回复
heyansheng12345
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 634.1
- 散金: 54
- 红花: 4
- 帖子: 71
- 在线: 4.9小时
- 虫号: 543364
- 注册: 2008-04-10
- 性别: GG
2楼2010-12-14 13:47:04
cyclone985
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 157.8
- 红花: 1
- 帖子: 175
- 在线: 34.9小时
- 虫号: 516204
- 注册: 2008-03-02
- 专业: 凝聚态物性 II ：电子结构
3楼2010-12-15 12:44:13
wangxn06
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 576.7
- 散金: 35
- 帖子: 509
- 在线: 250.5小时
- 虫号: 1156835
- 注册: 2010-11-27
- 专业: 流体力学
4楼2010-12-17 18:41:20
luoxitc
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1103.9
- 红花: 1
- 帖子: 37
- 在线: 7.6小时
- 虫号: 635873
- 注册: 2008-10-24
- 性别: GG
- 专业: 力学中的基本问题和方法
5楼2010-12-27 19:44:26
haiqing1987
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 293.9
- 帖子: 36
- 在线: 39.8小时
- 虫号: 1155042
- 注册: 2010-11-24
- 性别: GG
- 专业: 原子和分子物理
我有例子 刚刚研究过
★ ★
余泽成(金币+2):谢谢参与应助！ 2010-12-31 16:34:25
fagege2010(金币+1): 2011-01-08 23:46:01
fagege2010(金币+1): 2011-01-08 23:46:42
fagege2010(金币+1): 2011-03-07 09:00:45
余泽成(金币+2):谢谢参与应助！ 2010-12-31 16:34:25
fagege2010(金币+1): 2011-01-08 23:46:01
fagege2010(金币+1): 2011-01-08 23:46:42
fagege2010(金币+1): 2011-03-07 09:00:45
|
不废话 直接上代码 ： program main use IMSL ！调用IMSL,一定要用，而且得保证你的编译器有这个库 implicit none real :: A(3,3)=(/-1,-4,1,& 1,3,0,& 0,0,2 /) real :: eigenvalue(3),top real :: eigenvector(3,3) integer :: i eigenvalue=eig(A,v=eigenvector) ！A(n,n) 输入的矩阵，eigenvalue(n) 得到的本征矢，eigenvector(n,n)对应的本 ！证矢量，其中本征矢是降序排列 ！如果想得到复数数组 用eigenvalue=eig(A,w=eigenvector) do i=1,3 write(*,"('eigenvalue=',f5.2)"  eigenvalue(i)write(*,"('eigenvector=['3(f5.2'')']')" eigenvector(:,i)end do ！top 是验证是否归一 top=0.d0 do i=1,3 top=top+eigenvector(i,2)**2 enddo print*,top stop end [ Last edited by haiqing1987 on 2010-12-28 at 19:39 ] |
6楼2010-12-28 19:34:10
