[请教]请问关于矩阵特征值和特征向量的计算

不化顽石

请问，当我用
[v,d]=eig(A);
求出特征值及特征向量之后，想对其排个序，如果两边一起进行呢？
因为对特征值的排序相对简单，但如果两者不对应的话，得到的结果就没意义了。
躹躬，请指教。

eight

help sortrows

[ 本帖最后由 ChaChing 于 2010-7-5 20:51 编辑 ]

不化顽石

请问楼上，我用了这个函数，还是那个问题，特征值矩阵是搞定了
之后呢？如何把刚才算好的特征向量与排好了序的特征值矩阵联系起来？

eight

原帖由 不化顽石 于 2007-5-9 11:38 发表
请问楼上，我用了这个函数，还是那个问题，特征值矩阵是搞定了
之后呢？如何把刚才算好的特征向量与排好了序的特征值矩阵联系起来？

1. 
   
2. >> a=[2,5,1,-2,0,7,9,10]
   
3. a =
   
4.      2     5     1    -2     0     7     9    10
   
5. >> [b,c]=sort(a)
   
6. b =
   
7.     -2     0     1     2     5     7     9    10
   
8. 
   
9. c =
   
10.      4     5     3     1     2     6     7     8
    
11. >> sortrows([c',b'])
    
12. ans =
    
13.      1     2
    
14.      2     5
    
15.      3     1
    
16.      4    -2
    
17.      5     0
    
18.      6     7
    
19.      7     9
    
20.      8    10
    

复制代码

我这里的 b，c 就相当于你的特征向量和特征值。

建议下次自己稍微动手help一下

sffei

如求出的特征值里面有实数和复数，有没有什么命令可以选出实数呢

w89986581

index = find(imag(xx)<eps);
判断实数下标

sffei

谢谢w89986581 提醒

eight

原帖由 sffei 于 2007-5-9 17:52 发表
如求出的特征值里面有实数和复数，有没有什么命令可以选出实数呢

还可以用 isreal 命令

sffei

isreal 好像是对矩阵所有元素进行判断，并不能得到实数的索引，除非对矩阵进行循环

rocwoods

的确isreal是对整个矩阵进行判断，只要矩阵含有虚数，就返回false。
但是我们可以这样取出矩阵中的实数：
b=a(a-real(a)==0)

eight

原帖由 sffei 于 2007-5-10 00:59 发表
isreal 好像是对矩阵所有元素进行判断，并不能得到实数的索引，除非对矩阵进行循环

对，谢谢提醒，没想到这个命令是全局的，呵呵