matlab 求矩阵奇异值,MATLAB矩阵特征值和奇异值.

4.2  矩阵特征值和奇异值

对于n阶方阵A，求数λ和向量X，使得等式AX=λX成立，满足等式的数λ称为A的特征值，向量X称为A的特征向量。方程AX=λX和(A-λI)X=0是两个等价方程，要使方程(A-λI)X=0有非0解X，则必须使其行列式等于0，即|A-λI |=0。

由线性代数可知，行列式|A-λI |是一个关于λ的n阶多项式，因此方程|A-λI |=0是一个n次方程，有n个根(包含重根)。n个根就是矩阵A的n个特征值，每一个特征值对应无穷多个特征向量。所以，矩阵的特征值问题有确定的解，但特征向量问题没有确定的解。

4.2.1  特征值和特征向量的求取

特征值和特征向量在科学研究和工程计算中的应用非常广泛。在MATLAB中，计算矩阵A的特征值和特征向量的函数是eig(A)，常用的调用格式有以下3种。

(1)E=eig(A)：用于求矩阵A的全部特征值，构成向量E。

(2)[V,D]=eig(A)：用于求矩阵A的全部特征值，构成对角矩阵D，并求A的特征向量，构成V的列向量。

(3)[V,D]=eig(A,’nobalance’)：与上一种格式类似，只是上一种格式中是先对A作相似变换后再求矩阵A的特征值和特征向量，而本格式中则是直接求矩阵A的特征值和特征向量。

一个矩阵A的特征向量有无穷多个，eig函数只找出其中的n个，A的其他特征向量均可由这n个特征向量组合表示。

【例4-11】  简单实数矩阵的特征值示例。

>> A=[1,-3;2,2/3]

1.0