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

PAGE \* MERGEFORMAT2 矩阵特征值与特征向量的研究 目录 TOC \o 1-3 \h \z \u 一 矩阵特征值与特征向量研究的背景及意义 h 3 二、特征值与特征向量的定义及其性质 h 4 2.1 定义 h 4 2.2 性质 h 4 三 特征值及其特征向量的求法及其MATLAB的实现 h 5 3.1 QR方法 h 5 3.1.1 基本原理 h 5 3.1.2 具体实例 h 5 3.2 用多项式的方法来求解特征值 h 10 四 特征值与特征向量的简单应用 h 12 五 小结 h 16  一 矩阵特征值与特征向量研究的背景及意义 矩阵的特征值与特征向量是高等代数的重要组成部分，通过对矩阵特征值与特征向量的性质介绍，以及对矩阵特征值与特征向量理论的分析，将特征值与特征向量应用于方程组的求解问题是高等代数中的重要内容。 随着社会到的进步，计算机的飞速发展，高等代数这门课程已经渗透到各行各业里面。在许多方面都有着很重要的应用。在多数高等代数教材中,特征值与特征向量描述为线性空间中线性变换A的特征值与特征向量。从理论上来讲只要求出线性变换A的特征值和特征向量就可以知道矩阵A的特征值和特征向量。因此求矩阵的特征值与特征向量就变得尤为重要的引入是为了研究线性空间中线性变换A的属性。 在物理，力学，工程技术中有很多问题在数学上都归结为求矩阵的特征值和特征向量的问题。现在教材中给出的求解特征值和特征性向量的方法基本上都是通过求解特方程来求解。有时候特征方程会极其的麻烦。有一些文章中虽然给了初等行列变换的方法来较少计算量，但是仍未摆脱参数行列式计算的问题。本文中我们将首先讲解有关特征值和特征向量的相关知识，另外介绍一些简单实用的方法来求解矩阵的特征值与特征向量。 二、特征值与特征向量的定义及其性质 2.1 定义 设 A是n阶方阵，如果存在数λ 和n维非零向量x，使得 Ax =λ x成立，则称λ 为 A的特征值，x是A 的对应特征值λ 的特征向量。 2.2 性质 （1）λ0是A的特征值 （2）α是A的属于特征值λ0的特征向量的重要条件为α为齐次方程组λ0 （3）n阶矩阵在复数域上恰好有n个特征值（重根按重数计算）。 （4）n阶矩阵A为可逆矩阵的重要条件是A的特征值全不为0。 （5）A与AT （6）设A是可逆矩阵，如果λ0是A的一个特征值，对应的特征向量为α，则λ0- 三 特征值及其特征向量的求法及其MATLAB的实现 3.1 QR方法 3.1.1 基本原理 QR算法是计算矩阵特征值问题最有效的方法之一，也是普遍被用于工程实践中的一种方法。QR方法的思想是基于对于实的非奇异矩阵都可以分解为正交矩阵Q和上三角矩阵R的乘积，而且当R的对角元素符号取定时，分解是唯一的。 QR算法的基本步骤如下 （1）令A=A1，对A1进行正交分解，分解为正交矩阵Q A （2）然后将得到的因式矩阵Q1 A （3）以A2代替A1，重复以上步骤得到 A 性质1 所有的Ak 性质2 Ak A 其中Qk= 3.1.2 具体实例 例1 用QR算法求矩阵 A=5 - 的特征值。 解： 令A1=A,用施密特正交化过程将 A1= =0.9806-0.0377 0.1932-0.1038 将R1与 A2=R1 用A1 A 由A12 1.8789-λ 的根，求得为1+2i,1-2i 例2 已知矩阵A=21 程序代码如下 function [namda,time,data_na]=tzh(A,tol) if nargin==1; tol=1e-7 end %设置初始误差使之能进入循环 wucha=1 %记录迭代的次数 time=0 %如果误差没有满足精度，并且迭代次数在500次以内，可以循环迭代 %否则跳出循环 while (wuchatol)(time500) [q,r]=qr(A); A1=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=[2 1 0 1 3 1 0 1 4]; %调用方法函数 [namda,time,data_na]=tzh(a); disp(特征值为) namda disp(迭代次数为) t