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

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

銅蟲 (初入文壇)
|
不廢話 直接上代碼 : program main use IMSL !調(diào)用IMSL,一定要用,而且得保證你的編譯器有這個(gè)庫(kù) 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)對(duì)應(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ā)表 | |
|---|---|---|---|---|
|
[考研] 本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分 +3 | 里子木yy 2026-03-29 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 310求調(diào)劑 +9 | 爭(zhēng)取九點(diǎn)睡 2026-03-30 | 9/450 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +10 | 麻辣魷魚 2026-03-27 | 11/550 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 343求調(diào)劑 +5 | 愛(ài)羈絆 2026-03-28 | 5/250 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績(jī)338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 一志愿上海理工能源動(dòng)力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研]
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|