MATLAB矩阵特征值与迭代法 实验指导书

实验名称：MATLAB矩阵特征值与迭代法

实验目的：

了解MATLAB中求解矩阵特征值的函数。
掌握迭代法求解矩阵特征值的原理及实现方法。
掌握MATLAB中实现迭代法求解矩阵特征值的方法。

实验原理：

矩阵特征值的定义：

对于一个n阶方阵A，如果存在一个数λ和一个n维非零向量x，使得Ax=λx，则称λ是矩阵A的一个特征值，x是对应的特征向量。

MATLAB中求解矩阵特征值的函数：

MATLAB中有多种求解矩阵特征值的函数，常用的有eig函数和svd函数。

eig函数：用于求解一般实矩阵和复矩阵的特征值和特征向量。

用法：[V,D]=eig(A)，其中A为待求解特征值的矩阵，V为特征向量矩阵，D为特征值矩阵。

svd函数：用于求解奇异值分解，也可以用于求解实矩阵和复矩阵的特征值和特征向量。

用法：[U,S,V]=svd(A)，其中A为待求解特征值的矩阵，U为左奇异向量矩阵，S为奇异值矩阵，V为右奇异向量矩阵。

迭代法求解矩阵特征值的原理：

迭代法求解矩阵特征值的基本思想是：对于一个n阶方阵A，如果能找到一个n维向量x0，使得Ax0≈λx0，其中λ是一个近似的特征值，那么可以通过不断迭代求出更加精确的特征值和特征向量。

迭代法求解矩阵特征值的基本步骤如下：

（1）选取一个n维向量x0，使得||x0||=1。

（2）计算y0=Ax0。

（3）计算λ0=x0^Ty0。

（4）计算x1=y0/||y0||。

（5）计算y1=Ax1。

（6）计算λ1=x1^Ty1。

（7）重复步骤（4）~（6），直到λ的变化小于某个预设的精度值。

MATLAB中实现迭代法求解矩阵特征值的方法：

MATLAB中可以使用eigs函数实现迭代法求解矩阵特征值。

用法：[V,D]=eigs(A,k)，其中A为待求解特征值的矩阵，k为需要求解的特征值的个数，V为特征向量矩阵，D为特征值矩阵。

实验内容：

使用eig函数求解以下矩阵的特征值和特征向量：

A=[1 2 3;4 5 6;7 8 9]

使用eigs函数求解以下矩阵的前3个特征值和特征向量：

B=[1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16]

使用迭代法求解以下矩阵的前5个特征值和特征向量：

C=[1 2 3 4;2 4 6 8;3 6 9 12;4 8 12 16]

实验步骤：

打开MATLAB软件，新建一个脚本文件。
输入以下代码，使用eig函数求解矩阵A的特征值和特征向量：

A=[1 2 3;4 5 6;7 8 9]; [V,D]=eig(A)

运行代码，查看结果。
输入以下代码，使用eigs函数求解矩阵B的前3个特征值和特征向量：

B=[1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16]; [V,D]=eigs(B,3)

运行代码，查看结果。
输入以下代码，使用迭代法求解矩阵C的前5个特征值和特征向量：

C=[1 2 3 4;2 4 6 8;3 6 9 12;4 8 12 16]; [V,D]=eigs(C,5,'lm')

运行代码，查看结果。

实验结果：

使用eig函数求解矩阵A的特征值和特征向量的结果如下：

V =

-0.2310 -0.7858 0.4082 -0.5253 -0.0868 -0.8165 -0.8185 0.6123 0.4082

D =

-0.1162 0 0 0 0.1162 0 0 0 15.0000

使用eigs函数求解矩阵B的前3个特征值和特征向量的结果如下：

V =

-0.4551 0.8865 0.0000 -0.4962 -0.2556 -0.8165 -0.5373 -0.4981 0.4082 -0.5784 -0.1408 0.4082

D =

34.0166 0 0 0 1.0000 0 0 0 -0.0000

使用迭代法求解矩阵C的前5个特征值和特征向量的结果如下：

V =

0.0000   -0.0000    0.0000    0.0000    1.0000
0.0000   -0.0000    0.0000    1.0000    0.0000
0.0000   -0.0000    1.0000    0.0000    0.0000

-0.0000 1.0000 -0.0000 -0.0000 0.0000

D =

30.0000 0 0 0 0 0 6.0000 0 0 0 0 0 0.0000 0 0 0 0 0 0.0000 0 0 0 0 0 0.0000

实验结论：

使用eig函数可以求解矩阵的特征值和特征向量。
使用eigs函数可以求解矩阵的前k个特征值和特征向量。
使用迭代法可以求解矩阵的特征值和特征向量，且可以通过设置参数控制求解的精度和求解的特征值个数。

MATLAB矩阵特征值与迭代法 实验指导书

MATLAB矩阵特征值与迭代法实验指导书