矩阵特征值和特征向量的研究.docx

1、矩阵特征值与特征向量的研究2目录一 矩阵特征值与特征向量研究的背景及意义 .3二、特征值与特征向量的定义及其性质 .42.1 定义 .42.2 性质 4三 特征值及其特征向量的求法及其 MATLAB 的实现 .53.1 QR 方法 .53.1.1 基本原理 53.1.2 具体实例 53.2 用多项式的方法来求解特征值 .10四 特征值与特征向量的简单应用 .12五 小结 163一 矩阵特征值与特征向量研究的背景及意义矩阵的特征值与特征向量是高等代数的重要组成部分，通过对矩阵特征值与特征向量的性质介绍，以及对矩阵特征值与特征向量理论的分析，将特征值与特征向量应用于方程组的求解问题是高等代数中的重

2、要内容。随着社会到的进步，计算机的飞速发展，高等代数这门课程已经渗透到各行各业里面。在许多方面都有着很重要的应用。在多数高等代数教材中,特征值与特征向量描述为线性空间中线性变换 A 的特征值与特征向量。从理论上来讲只要求出线性变换 A 的特征值和特征向量就可以知道矩阵 A 的特征值和特征向量。因此求矩阵的特征值与特征向量就变得尤为重要的引入是为了研究线性空间中线性变换 A 的属性。在物理，力学，工程技术中有很多问题在数学上都归结为求矩阵的特征值和特征向量的问题。现在教材中给出的求解特征值和特征性向量的方法基本上都是通过求解特方程来求解。有时候特征方程会极其的麻烦。有一些文章中虽然给了初等行列变

3、换的方法来较少计算量，但是仍未摆脱参数行列式计算的问题。本文中我们将首先讲解有关特征值和特征向量的相关知识，另外介绍一些简单实用的方法来求解矩阵的特征值与特征向量。4二、特征值与特征向量的定义及其性质2.1 定义设 A 是 n 阶方阵，如果存在数 和 n 维非零向量 x，使得 Ax = x 成立，则称 为 A 的特征值， x 是 A 的对应特征值 的特征向量。2.2 性质（1） 是 A 的特征值0 (0)=|=0（2） 是 A 的属于特征值 的特征向量的重要条件为 为齐次方程组 0 非零解。（ 0） =0（3）n 阶矩阵在复数域上恰好有 n 个特征值（重根按重数计算） 。（4）n 阶矩阵 A

4、为可逆矩阵的重要条件是 A 的特征值全不为 0。（5）A 与 有相同的特征值。（6）设 A 是可逆矩阵，如果 是 A 的一个特征值，对应的特征向量为 ，则0 的一个特征值，对应的特征向量仍然为 。10是 1 5三 特征值及其特征向量的求法及其 MATLAB 的实现3.1 QR 方法3.1.1 基本原理QR 算法是计算矩阵特征值问题最有效的方法之一，也是普遍被用于工程实践中的一种方法。QR 方法的思想是基于对于实的非奇异矩阵都可以分解为正交矩阵 Q 和上三角矩阵 R 的乘积，而且当 R 的对角元素符号取定时，分解是唯一的。QR 算法的基本步骤如下（1）令 ，对 进行正交分解，分解为正交矩阵 和上

5、三角矩阵 的=1 1 1 1乘积：1=11（2）然后将得到的因式矩阵 反序相乘，得到：1和 12=11（3）以 代替 ，重复以上步骤得到 ，所以得到的 QR 算法的计算公式2 1 3为：=+1=性质 1 所有的 都相似，它们具有相同的特征值。性质 2 的 QR 分解式为=其中 ,=12. =123.1.2 具体实例例 1用 QR 算法求矩阵6A=5 2 5 11 0 3 20 2 2 30 0 1 2的特征值。解： 令 ,用施密特正交化过程将 分解为1= 1=111=0.98060.0377 0.19320.10380.1961 0.18870.88040.41920 0.9813 0.176

6、1 0.07040 0 0.3962 0.89895.00921.96125.49120.39220 2.0381 1.5852 2.52880 0 2.5242 3.12940 0 0 0.7822将 逆序相乘，求出1与 1 2= =2114.6517 5.9508 1.5299 0.23900.3977 1.9401 -2.5171 1.53610 2.4770 -0.8525 3.12940 0 0.3099 -0.7031用 代替 A 重复上面过程，计算 11 次得112=4.0000 0 1.8789 3.5901 0 1.3290 0.1211 0 0 0 1.0000 由 不难看

7、出，矩阵 A 的一个特征值是 4，另一个特征值是 -1，其他两个特征12值是方程=0|1.8789 3.59011.32900.1211|的根，求得为 1+2i,1-2i例 2 7已知矩阵 A= ,采用 QR 方法计算 A 的全部特征值。2 1 01 3 10 1 4程序代码如下function namda,time,data_na=tzh(A,tol)if nargin=1;tol=1e-7end%设置初始误差使之能进入循环wucha=1%记录迭代的次数time=0%如果误差没有满足精度，并且迭代次数在 500 次以内，可以循环迭代%否则跳出循环while (wuchatol)A1=r*q;

8、tz0=diag(A1);tz1=diag(A);wucha=norm(tz0-tz1);%迭代赋值A=A1;time=time+1;data_na(time,:)=tz1;endnamda=tz1;%用 QR 方法计算矩阵特征值a=2 1 01 3 180 1 4;%调用方法函数namda,time,data_na=tzh(a);disp(特征值为)namdadisp(迭代次数为)time%用于输出数据n1=length(data_na);%n2 为数组n2=(1:n1);%temp1 为迭代序列与特征值组成的向量temp1=n2,data_na;%第一个特征值subplot(2,2,1:2

9、)plot(data_na(:,1)title(第一个特征值)grid%第二个特征值subplot(2,2,3)plot(data_na(:,2)title(第二个特征值)grid%第三个特征值subplot(2,2,4)plot(data_na(:,3)title(第三个特征值)grid9输出结果为：特征值为namda =4.73213.00001.2679迭代次数为time =2210由图像可以看出在迭代的前几次可能会有一些波动，但是逐渐趋于平稳，总体而言，QR 方法是计算矩阵特征值的一个比较好的方法。3.2 用多项式的方法来求解特征值我们知道，求 n 阶方阵 A 的特征值就是求代数方程(

10、)=|=0的根。 称为 A 的特征多项式。上式展开为()()=+11+其中 ， . 为多项式 的系数。1 2 ()从理论上来讲，求 A 得特征值可分为两步：11第一步：直接展开行列式 求出多项式 ；| ()第二步：求代数方程 =0 的根，即特征值。()对于低阶矩阵，这种方法显然是可行的。但是对于高阶矩阵，计算量则非常的大，这种方法就有其自身的弊端。这里我们将介绍 F-L 方法来求特征方程中的多项式 的系数，也就是求多项式 。由于代数方程求根问题的关键是确() ()定矩阵 A 的特征多项式 ，所以这种方法为多项式方法求特征值问题。()记矩阵 A= 的对角线元素之和为()trA=11+22+利用递

11、归的概念定义以下 n 个矩阵 (K=1、2、 、 、 、n):1=, 1=12=(11), 2=122.=(11), =1.=(11), =1可以证明上式中 ，k=1,2,3.n,即是所求 A 特征多项式 的各项系数，用 ()上式求矩阵特征多项式系数的方法称为 F-L 法，相应的特征方程为(1)(11)=0而且可证矩阵 A 的逆矩阵可表示为1=1(11)特征向量的求法当矩阵 A 的特征值确定以后，将这些特征值逐个代入齐次线性方程组x=0 中，由于系数矩阵 的秩小于矩阵 的阶数 n，因此虽然（ ） 有 n 个方程 n 个未知数，但实际上是解有 n 个未知数的相互独立到 r 个方程（rn）.当矩阵

12、 A 的所有特征值互不相同的时候，这样的问题中要解的其次性12方程组中有 n-1 个独立的方程，其中含有 n 个特征向量分量，因此特征向量分量中至少有一个需要任意假设其值，才能求出其他的特征分量。四 特征值与特征向量的简单应用在经济发展与环境污染的增长模型方面的应用在现在的时代发展中，经济增长飞速发展。但是随着经济增长的同时，环境污染也越发的严重。环境的治理称为当今社会需要注意的有一个关键的问题。所以探讨环境与经济增长之间的关系就变得尤为的重要。在这方面矩阵的特征值与特征向量有着一定程度上的应用，可建立如下数学模型： 设 分别为某地区目前的环境污染水平与经济发展水平， 分别为该地0,xy 1,

13、yx区若干年后的环境污染水平和经济发展水平，且有如下关系： 10432xyy令01,xy4132A则上述关系的矩阵形式为 10aA此式反映了该地区当前和若干年后的环境污染水平和经济发展水平之间的关系.如01xy10 0415132A则由上式得由此可预测该地区若干年后的环境污染水平和经济发展水平. 一般地，若令 分别为该地区 t 年后的环境污染水平与经济发展水平，则tyx,经济发展与环境污染的增长模型为令 14(1,2)32tttxytky ttxy13则上述关系的矩阵形式为 1,2ttAk由此，有102233010,().tttA 由此可预测该地区 t 年后的环境污染水平和经济发展水平.下面作

14、进一步地讨论：由矩阵 A 的特征多项式 41| (4)232EA得 A 的特征值为 124,对 ，解方程 得特征向量 410)(X1对 ，解方程 得特征向量 1=2()EAx2显然, 线性无关12,下面分三种情况分析： 第一种： 01一个性质:若 是矩阵 的属于特征值 的特征向, 也是 的属于特征值aAakA的特征向量度（*）k由(*)及特征值与特征向量的性质知，0114tttttA即 或 4ttxytttxy此式表明：在当前的环境污染水平和经济发展水平的前提下， 年后，当t14经济发展水平达到较高程度时，环境污染也保持着同步恶化趋势. 第二种： 021，所以不讨论此种情况0y第三种： 017

15、不是特征值,所以不能类似分析。但是 可以由 唯一线性表出0021，来： 2103由(*)及特征值与特征向量的性质 0121212 ,(3)3434ttt tttttttAA即 ,ttxy342,tt34tt由此可预测该地区年后的环境污染水平和经济发展水平. 因无实际意义而在第二种情况中未作讨论，但在第三种情况的讨论中2仍起到了重要作用. 由经济发展与环境污染的增长模型易见，特征值和特征向量理论在模型的分析和研究中获得了成功的应用。在其他方面的应用简述在信息处理上的意义由于这些投影的大小代表了 A 在特征空间各个分量的投影，那么我们可以使用最小 2 乘法，求出投影能量最大的那些分量，而把剩下的分

16、量去掉，这样最大限度地保存了矩阵代表的信息，同时可以大大降低矩阵需要存储的维度，简称 PCA 方法。3线性变换 PCA 可以用来处理图像。如 2 维的人像识别：我15们把图像 A 看成矩阵，进一步看成线性变换矩阵，把这个训练图像的特征矩阵求出来(假设取了 n 个能量最大的特征向量)。用 A 乘以这个 n 个特征向量，得到一个 n 维矢量 a，也就是 A 在特征空间的投影。今后在识别的时候同一类的图像(例如，来自同一个人的面部照片)，认为是 A 的线性相关图像，它乘以这个特征向量，得到 n 个数字组成的一个矢量 b，也就是 B 在特征空间的投影。那么 a 和 b 之间的距离就是我们判断 B 是不

17、是 A 的准则。又如 Google 公司的 PageRank，也是通过计算一个用矩阵表示的图。这个图代表了整个 Web 各个网页“节点”之间的关联。用特征向量来对每一个节点打“特征值”分。16五 小结在这个信息飞速发展的时代，我们的科技正在越来越进步。各种先进的产品层出不穷。人们在惊叹于社会科技进步的同时不能忘了一些基础学科在其中起到的重要作用。在本文中我们就简单的介绍了线性代数中特征值与特征向量的一些研究。在这个大数据的时代，许多数据都是以矩阵的形式展现在大家到眼前。在矩阵中标有一个很重要的量就是特征值和特征向量。本文主要分为了四个部分介绍了特征值与特征向量的相关知识。首先在第一部分我们知道

18、了特征值与特征向量的研究背景，我们也明白了其重要性。接下来我们介绍了特征值与特征向量到定义和其一些常用的性质，使我们对其有了更加详细的了解。然后我们介绍了两种常用的求解特征值和特征向量的方法，当然为了区别我们并没有介绍本科期间书本上讲解的一些最基本的例如行列变换的解法。我们介绍了 QR 方法并且用 MATLAB 实践了一下，因为在实际生活中我们需要用计算机来解决这些问题。然后我们还介绍了 F-L 的方法来求解矩阵的特征值和特征向量。最后我们讲解在生活中特征值和特征向量的应用，例如在经济发展与环境污染模型中的应用和在信息处理上的意义。希望通过本文的讲解，大家对特征值和特征向量有一个更加理性更加深刻的理解。