作者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