| 6 | 1/1 | 返回列表 |
| 查看: 1037 | 回復(fù): 5 | |||
fagege2010鐵蟲 (小有名氣)
|
[交流]
【求助】急求!!
|
|
我在前面的子程序中計算出了一個矩陣AR(N,N) 現(xiàn)在想用ISML求它的本征值和本征矢 請問如何求;如何把子程序中的這個矩陣中的數(shù)據(jù)傳過去?具體怎么編? 問題已經(jīng)解決謝謝各位 [ Last edited by fagege2010 on 2011-1-8 at 23:54 ] |

銀蟲 (小有名氣)
銀蟲 (小有名氣)
金蟲 (正式寫手)
新蟲 (初入文壇)

銅蟲 (初入文壇)
|
不廢話 直接上代碼 : program main use IMSL !調(diào)用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)對應(yīng)的本 !證矢量,其中本征矢是降序排列 !如果想得到復(fù)數(shù)數(shù)組 用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 是驗(yàn)證是否歸一 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 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +7 | 深V宿舍吧 2026-03-29 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 332求調(diào)劑 +14 | 032500 2026-03-25 | 14/700 |
|
|
[考研] 309求調(diào)劑 +15 | 誰不是少年 2026-03-29 | 15/750 |
|
|
[考研] 化工專碩求調(diào)劑 +5 | question挽風(fēng) 2026-03-24 | 5/250 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
[考研] 求調(diào)劑 +7 | 研研,接電話 2026-03-24 | 8/400 |
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10 | 景隨風(fēng) 2026-03-29 | 15/750 |
|
|
[考研] 337求調(diào)劑 +6 | 《樹》 2026-03-29 | 6/300 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
[考研] 考研調(diào)劑 +7 | 小蠟新筆 2026-03-29 | 7/350 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|