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