矩阵特征值和特征向量的研究.docx
- 文档编号:24746581
- 上传时间:2023-06-01
- 格式:DOCX
- 页数:13
- 大小:83.59KB
矩阵特征值和特征向量的研究.docx
《矩阵特征值和特征向量的研究.docx》由会员分享,可在线阅读,更多相关《矩阵特征值和特征向量的研究.docx(13页珍藏版)》请在冰豆网上搜索。
矩阵特征值和特征向量的研究
矩阵特征值与特征向
量的研究
一矩阵特征值与特征向量研究的背景及意义3
二、特征值与特征向量的定义及其性质4
2.1定义4
2.2性质4
三特征值及其特征向量的求法及其MATLAB的实现5
3.1QR方法5
3.1.1基本原理5
3.1.2具体实例5
3.2用多项式的方法来求解特征值10
4特征值与特征向量的简单应用12
5小结16
矩阵特征值与特征向量研究的背景及意义
矩阵的特征值与特征向量是高等代数的重要组成部分,通过对矩阵特征值与特征向量的性质介绍,以及对矩阵特征值与特征向量理论的分析,将特征值与特征向量应用于方程组的求解问题是高等代数中的重要内容。
随着社会到的进步,计算机的飞速发展,高等代数这门课程已经渗透到各行各业里面。
在许多方面都有着很重要的应用。
在多数高等代数教材中,特征值与特征向量描述为线性空间中线性变换A的特征值与特征向量。
从理论上来讲只要求出线性变换A的特征值和特征向量就可以知道矩阵A的特征值和特征向量。
因此求矩阵的特征值与特征向量就变得尤为重要的引入是为了研究线性空间中线性变换A的属性。
在物理,力学,工程技术中有很多问题在数学上都归结为求矩阵的特征值和特征向量的问题。
现在教材中给出的求解特征值和特征性向量的方法基本上都是通过求解特方程来求解。
有时候特征方程会极其的麻烦。
有一些文章中虽然给了初等行列变换的方法来较少计算量,但是仍未摆脱参数行列式计算的问题。
本文中我们将首先讲解有关特征值和特征向量的相关知识,另外介绍一些简单实用的方法来求解矩阵的特征值与特征向量。
二、特征值与特征向量的定义及其性质
2.1定义
设A是n阶方阵,如果存在数入和n维非零向量x,使得Ax=入x成立,则称入为A的特征值,x是A的对应特征值入的特征向量。
2.2性质
(1)仏是A的特征值O匚入)=1血E—川=。
(2)口是A的属于特征值血的特征向量的重要条件为口为齐次方程组a声-a)x=o非零解。
(3)n阶矩阵在复数域上恰好有n个特征值(重根按重数计算)。
(4)n阶矩阵A为可逆矩阵的重要条件是A的特征值全不为0。
(5)A与’有相同的特征值。
(6)设A是可逆矩阵,如果兀是A的一个特征值,对应的特征向量为。
,则
c心円的一个特征值,对应的特征向量仍然为口。
3特征值及其特征向量的求法及其MATLAB勺实现
3.1QR方法
3.1.1基本原理
QR算法是计算矩阵特征值问题最有效的方法之一,也是普遍被用于工程实践中的一种方法。
QR方法的思想是基于对于实的非奇异矩阵都可以分解为正交矩阵Q和上三角矩阵R的乘积,而且当R的对角元素符号取定时,分解是唯一的
QR算法的基本步骤如下
(1)令A=对州进行正交分解,分解为正交矩阵Qi和上三角矩阵©的
乘积:
⑵然后将得到的因式矩阵Qi和反序相乘,得到:
⑶以“2代替"d,重复以上步骤得到虫3,所以得到的QR算法的计算公式为:
:
如=QkRk
性质1所有的'都相似,它们具有相同的特征值。
性质2'的QR分解式为
人=Q風
其中
3.1.2具体实例
用QR算法求矩阵
[5
-2
—S—11
1
0
-3
2
0
2
2
-3
0
—1
0
1
-2
的特征值。
解:
令①二①用施密特正交化过程将戶】分解为
r
r0.9806-0.03770.1932-0JO38i
r5.0092
L9612-5.4912
-0.3922]
0.19610.1887-0.8804-0.4192
o
2.0381
L5852
-25288
00.98130.17610.0704
P
0
2.5242
-3.1294
000+3962-0.8989j
0
0
07322
将逆序相乘,求出去
4.65175.95081.52990.23901
0,39771.9401・2^171L5361
0£1770・(L8525:
k12别
A』iQi二
()00.3()95)-0.70311
A
用①代替A重复上面过程,计算11次得
4.0000*4*
01,8789-3,5901*
0L32900.1211*
000-1.0000
由戶蛇不难看出,矩阵A的一个特征值是4,另一个特征值是-1,其他两个特征值是方程
11.8789-X-3.5901
|1.32900J2H-A=0
的根,求得为1+2i,1-2i
2101
131
已知矩阵A」°1引,采用QR方法计算A的全部特征值。
程序代码如下
function[namda,time,data_na]=tzh(A,tol)
ifnargin==1;
tol=1e-7
end
%设置初始误差使之能进入循环
wucha=1
%记录迭代的次数
time=0
%如果误差没有满足精度,并且迭代次数在500次以内,可以循环迭代
%否则跳出循环
while(wucha>tol)&(time<500)
[q,r]=qr(A);
A仁r*q;
tz0=diag(A1);
tz1=diag(A);
wucha=norm(tz0-tz1);
%迭代赋值
A=A1;
time=time+1;
data_na(time,:
)=tz1;
end
namda=tz1;
%用QR方法计算矩阵特征值a=[210
131
014];%调用方法函数[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)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('第三个特征值')
grid
输出结果为:
特征值为
namda=
4.7321
3.0000
1.2679
迭代次数为
time=
由图像可以看出在迭代的前几次可能会有一些波动,但是逐渐趋于平稳,总体而
言,QR方法是计算矩阵特征值的一个比较好的方法。
3.2用多项式的方法来求解特征值
我们知道,求n阶方阵A的特征值就是求代数方程
g)=\A-A1\=0
的根。
称为A的特征多项式。
上式展开为
(p(X)=An+p1An_1+.„+pn
其中卩1,旳…….再「为多项式申㈤的系数。
从理论上来讲,求A得特征值可分为两步:
第一步:
直接展开行列式|求出多项式'■';
第二步:
求代数方程I匕0的根,即特征值。
对于低阶矩阵,这种方法显然是可行的。
但是对于高阶矩阵,计算量则非常的大,
这种方法就有其自身的弊端。
这里我们将介绍F-L方法来求特征方程中的多项式
-的系数,也就是求多项式|■'0由于代数方程求根问题的关键是确定矩阵A
的特征多项式,,所以这种方法为多项式方法求特征值问题。
记矩阵A=的对角线元素之和为
trA=%+旳2+"+%
利用递归的概念定义以下n个矩阵巧(K=1、2、、、、n):
=At
=p2=
%=i-Pk=
可以证明上式中-,k=1,2,3….n,即是所求A特征多项式」的各项系数,用上式求矩阵特征多项式系数的方法称为F-L法,相应的特征方程为
ti11
(一1)(A1_***~Pft)=0
而且可证矩阵A的逆矩阵可表示为
_11
rn
特征向量的求法
当矩阵A的特征值确定以后,将这些特征值逐个代入齐次线性方程组
(A-Ai}x=0中,由于系数矩阵A-嗣的秩小于矩阵A-的阶数n,因此虽然有n个方程n个未知数,但实际上是解有n个未知数的相互独立到r个方程(r 4特征值与特征向量的简单应用 在经济发展与环境污染的增长模型方面的应用 在现在的时代发展中,经济增长飞速发展。 但是随着经济增长的同时,环境污染也越发的严重。 环境的治理称为当今社会需要注意的有一个关键的问题。 所 以探讨环境与经济增长之间的关系就变得尤为的重要。 在这方面矩阵的特征值与特征向量有着一定程度上的应用,可建立如下数学模型: 设xo,yo分别为某地区目前的环境污染水平与经济发展水平,Xi,yi分别为该地区 若干年后的环境污染水平和经济发展水平,且有如下关系: Xi4xoyo yi3xo2yo 令 xox41 o,iA yoyi32 则上述关系的矩阵形式为印Aao此式反映了该地区当前和若干年后的环境污染水平和经济发展水平之间的关系如 Xo o yo 5i 55i5o 则由上式得 由此可预测该地区若干年后的环境污染水平和经济发展水平 般地,若令Xt,yt分别为该地区t年后的环境污染水平与经济发展水平,则经 济发展与环境污染的增长模型为 论: 由矩阵A的特征多项式 得A的特征值为14, F面分三种情况分析: 1 第一种: 01 k的特征向量度(*) 由(*)及特征值与特征向量的性质知, 即xt4t1或xtyt4t yt1 此式表明: 在当前的环境污染水平和经济发展水平的前提下,t年后,当经济发展水平达到较高程度时,环境污染也保持着同步恶化趋势. 第二种: 02 1 2 Qy0 20, 所以不讨论此种情况 第三种: 1 07 Q0不是特征值,所以不能类似分析。 但是0可以由1,2唯一线性表出来: 由(*)及特征值与特征向量的性质 tAt0At(3122)3At12At 由此可预测该地区年后的环境污染水平和经济发展水平. 2因无实际意义而在第二种情况中未作讨论,但在第三种情况的讨论中仍 起到了重要作用. 由经济发展与环境污染的增长模型易见,特征值和特征向量理论在模型的分析和研究中获得了成功的应用。 在其他方面的应用简述 在信息处理上的意义 由于这些投影的大小代表了A在特征空间各个分量的投影,那么我们可以使用最小2乘法,求出投影能量最大的那些分量,而把剩下的分量去掉,这样最大限度地保存了矩阵代表的信息,同时可以大大降低矩阵需要存储的维度,简称PCA方法。 [3]线性变换PCA可以用来处理图像。 如2维的人像识别: 我们把图像A看成矩阵,进一步看成线性变换矩阵,把这个训练图像的特征矩阵求出来(假设取了n个能量最大的特征向量)。 用A乘以这个n个特征向量,得到一个n维矢量a,也就是A在特征空间的投影。 今后在识别的时候同一类的图像(例如,来自同一个人的面部照片),认为是A的线性相关图像,它乘以这个特征向量,得到n个数字组成的一个矢量b,也就是B在特征空间的投影。 那么a和b之间的距离就是我们判断B是不是A的准则。 又如Google公司的PageRank也是通过计算一个用矩阵表示的图。 这个图代表了整个Web各个网 页“节点”之间的关联。 用特征向量来对每一个节点打“特征值”分。 5小结 在这个信息飞速发展的时代,我们的科技正在越来越进步。 各种先进的产品层出不穷。 人们在惊叹于社会科技进步的同时不能忘了一些基础学科在其中起到的重要作用。 在本文中我们就简单的介绍了线性代数中特征值与特征向量的一些研究。 在这个大数据的时代,许多数据都是以矩阵的形式展现在大家到眼前。 在矩阵中标有一个很重要的量就是特征值和特征向量。 本文主要分为了四个部分介绍了特征值与特征向量的相关知识。 首先在第一部分我们知道了特征值与特征向量的研究背景,我们也明白了其重要性。 接下来我们介绍了特征值与特征向量到定义和其一些常用的性质,使我们对其有了更加详细的了解。 然后我们介绍了两种常用的求解特征值和特征向量的方法,当然为了区别我们并没有介绍本科期间书本上讲解的一些最基本的例如行列变换的解法。 我们介绍了QR方法并且用 MATLAB实践了一下,因为在实际生活中我们需要用计算机来解决这些问题。 然后我们还介绍了F-L的方法来求解矩阵的特征值和特征向量。 最后我们讲解在生活中特征值和特征向量的应用,例如在经济发展与环境污染模型中的应用和在信息处理上的意义。 希望通过本文的讲解,大家对特征值和特征向量有一个更加理性更加深刻的理解。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 矩阵 特征值 特征向量 研究