PTT數位生活區 / Fortran

[問題] 請問矩陣對角化

看板Fortran作者 (新生活!~)時間15年前 (2009/04/24 15:49), 編輯推噓2(2推 0噓 4→)
留言6則, 2人參與, 最新討論串1/1
我要做一個3x3的矩陣對角化 利用"Numerical Recipes in fortran"第二版 第460頁的jacobi這個subroutine下去代 照理說應該可以求出eigenvalue和eigenvactor 可是我現在算出來的結果都只有eigenvalue是對的,eigenvactor是錯的 我直接寫一個小程式就只有run這個subroutine不做多餘的事,答案還是錯 從書上抄上去的地方我也重覆檢查了n次 請問有可能會是書上提供的程式有問題嗎??這本應該是聖經吧 囧 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.109.103.226
推
04/26 22:59, , 1F
矩陣對角化這個我很熟，會不會你給的矩陣不能對角化?
04/26 22:59, 1F
→
04/26 23:00, , 2F
如果用Numerical Recipes的話還可以試試TRED2加上TQLI
04/26 23:00, 2F
→
04/26 23:02, , 3F
或是你輸入的參數不對呢?導致JACOBI只計算本徵值?
04/26 23:02, 3F
→
04/28 02:09, , 4F
我給的是實數對稱矩陣,這是jacobi要求的條件,參數就我的理
04/28 02:09, 4F
→
04/28 02:09, , 5F
解應該是沒有設不對,再不行也只好用TRED2加TQLI去對照了
04/28 02:09, 5F
推
04/28 09:30, , 6F
可以給我那個3x3的矩陣嗎?我也來試試看。
04/28 09:30, 6F
‣ 更多 sinlin 的文章
文章代碼(AID): #19yMxqHQ (Fortran)
Fortran 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 sinlin 的文章
文章代碼(AID): #19yMxqHQ (Fortran)