求解矩阵特征值问题地一种新算法——非线性算法.docx

中南大学硬士毕业论文： 求解矩阵特征值闩寇的一种新算法一非线性算法 第一章特征值问题与非线性算法 特征值问题的提出，出现在50年代后期。由于计算机的出现与计 算机技术的发展，使得人们进行繁重的复杂计算成为可能。从此人们 便开始考虑利用计算机计算特征值和特征向量。在这方面，J．H． 砌蝤nson，G．W．Stewem，G．BMoler，B．N．Parlett等为此作出了重要 贡献。 §1．1特征值问题在工程领域中的应用 计算矩阵特征值特征向量的问题涉及许多领域，最早出现于近似 求解各种物理系统的“稳恒运动”或“微小”运动之中。 下面考虑一根长度为f的等截面均匀弦，两端被固定并且沿轴上 绷紧，一端置于原点O，另一端置于横坐标l处。在长度为出的元 素上该弦受到垂直方向的力zm)出的作用(r为弦的张力)。我们要研 究在垂直平面上的“微小”运动，即寻求一个定义在osx≤z，f≥o上 的函数”∽O，它表示在横坐标x处于时刻t时的弦的垂直变形。可 以证明，函数甜应该满足以下的偏微分方程： 专窑㈤一襄c咖mm。川一。 其中c=√形(p是弦的线密度)。它称为一维波动方程。 在没有外力作用时(，卸)，可以求弦的稳恒运动。一般地，一维波 中南大学硕士毕业论文： 求解矩阵特征值问题的一种新算法一非线性算法 动方程可以变为： 。一(对“o：{球(∞，一(o 这化为函数和v的两个微分方程： 一M弋工)=五搿(工)，一v。(f)=知21，(f)， 五是待定系数。由于 碳O，垆坝|lf)=o，对一切f， 故问题归结为求实数A与不恒等于零的函数“，使得 I一越。(x)=．五牡(x)，0≤工≤z IM(o)=甜(D=o 这就构成一个二阶导算子的特征值问题。 在信号处理的许多应用中，有用的许多信息必须从一个在宽带噪声中观测的低秩向量过程抽取。数据矩阵的协方差矩阵的特征值分解可以用来揭开隐藏的低秩信号。特别地，利用特征值的大小，我们可以将协方差矩阵的几何子空间即特征子空间分为信号子空间和噪声子空间两部分。所谓的特征子空间法就是用信号或噪声子空间进行低秩信息的提取。 最近十几年，特征子空间方法广泛应用于高分辨率谱估计和ARMA建模，谐波恢复、传感器阵列处理、系统辨识、信号增强，甚至于滤波器设计和模式识别等许多问题中。 在现代控制论中，系统的稳定性问题也离不开矩阵特征值问题的讨论。 中南大学硬士毕业论文： 求解矩阵特征值闩题的一种新算法一非线性算法 §1．2矩阵特征问题的算法概述 矩阵的特征值问题，不管在理论上还是在工程计算中，都是一个 十分重要的问题。长期以来，人们对这个问题的研究一直没有停止过。 至今为止，已出现了许许多多各有特色的求解矩阵特征值问题的方 法，如经典的LU分解、Ja∞bi方法、QR或双重步QR方法、二 分／多分(BS／MS)法、分治法、同伦连续法、迭代法等等。 二分／多分法是当前最常用的求解实对称矩阵的并行方法、这个 算法也可解He砌血i粕矩阵特征值问题、实对称矩阵和Hc砌施粕矩阵 的广义特征值问题，更适应于解实对称三对角线特征值问题。 设4是实对称三对角线矩阵， q 62 62 63 “一，％。6月 6。 4的五矩阵4一五I的各级前主子式为鼎五)，f=1，2 ．一，则有州工)具 有stll咖序列性质，即序列的相邻项符号不同数等于A的小于五的特 征值个数。在实际计算中，通常将以五)改用序列{献f)}： 口l(互)i口1-五， 似工)==(即五)_6产／雏1(五)，卢2，．．．，珂。 在这个序列{“i)}(f=2 ．，功中，负数项的个数等于4的小于2的 值的个数。二分法的思想就是对于某个区间k6】，计算区间中点 中南大学磺士毕业论文： 求解矩阵特征值问题的一种新算法—章线性算法 卢!芸上的序列似D，由“jf)的负数项个数确定特征值的分布情况。 通过逐步计算可得到只含一个特征值的区间【cI田，再将区间二等分， 逐次计算||}次，就可确定特征值的精度g。其中仁【lo以皇≥)】。 口 实对称矩阵的特征值问题可以先用Household盯变换、Giv曲s交 换把它转化为实对称三对角线矩阵的特征值问题。 分治法是一个较新的适应于并行运算的矩阵特征值问题的算法。 许多经验表明，分治法在共享存储机上的功效比二分／多分法还好。 分治法的思想是把要求的特征值问题的矩阵4分解成两个新的简单 的矩阵，然后分别对这两个简单矩阵求特征值，从而构成一个比原问 题容易解的特征值问题。容易求解的原因是给出了较好的近似值估计 并且矩阵有特殊的形状。 设实对称三对角矩阵 q 62 62 色 63 6¨吒一1 6。 “ 口。 可以把r化为产(互疋]+趔。 一般来说，五、正是把丁对半分得来的。设矗的阶为刀。，疋的阶 为，Ml，6=6^l，z是一个向量扛 ，其中第栉l、地个元素为l。 中南大