| 查看: 828 | 回复: 4 | |||
[求助]
有没有人用过lapack95求过矩阵的特征向量的。。弄了一整天了，还是算不对 已有2人参与
|
|
好不容易找到到一个geev函数，call geev(a, wr, wi [,vl] [,vr] [,info])，参数是这样设置的，a是原矩阵，vl和vr就是左右特征向量，右特征向量满足A*vr = r*vr。 今天拿了个简单的矩阵a(1,1)=1,a(1,2)=2,a(2,1)=3,a(2,2)=4做测试，我不管是手算还是用Mathematica算，特征向量都应该是{(1/6)* (-3 + Sqrt[33]), 1}, {(1/6)* (-3 - Sqrt[33]), 1}这两个，也满足A*vr = r*vr，但是用geev函数算出来的结果是图上这样的，完全不满足A*vr = r*vr这个等式，不过特征值却是正确的，有没有大神知道是怎么回事，是不是我哪里用错了。。  111.png |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有49人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
三无产品还有机会吗
已经有6人回复
感谢参与，应助指数 +1
|
本帖内容被屏蔽 |
2楼2015-08-18 00:10:15
3楼2015-08-18 00:12:39
4楼2015-08-18 19:45:48
5楼2015-08-19 17:02:43