特征值问题的MATLAB实践.docxVIP

特征值问题的 MATLAB 实践 !鱼:塑 ScienceandTechnologyInnovationHerald 特征值问题的 MATLAB 实践① 赵娜吕剑峰 (内蒙古科技大学数理与生物工程学院数学教研室内蒙古包头 014010) 工程技术 摘要 :特征值理论是线性代数中的一个重要内容 ,当方阵阶数较高时 ,实际计算比 较繁琐?本文从实际案例入手,利~MATIAB软件讨论 了求解特征值问题的全过程 ,并对案例加以分析 . 关键词 :MATLAB 特征值特征向量 中图分类号：0151.26文献标识码:A文章编号：1674— 098x(2010)10(c) — 0112— 01 在计算机广泛应用的今天 ,线性代数 在计算技术中的基础性地位也日益突出 , 计算机图形学 ,计算机辅助设计 ,密码学 , 虚拟现实等技术无不以线性代数为其理论 和算法基础的一部分 ,这时用代数方法解 决实际问题也渗透到了各个领域 .但当前 的线性代数课程教学并不尽如人意 ,传统 的教学方式和内容过分追求逻辑的严密性 和理论体系的完整性 ,重理论而轻实践 ,结 果不但没有提高学生的实践能力和创新能 力,而且也不能让学生有效地衔接后续专 业课程 .为了改变这种现状 ,我们有必要引 入一些计算机软件工具来增强线性代数教 学，减少抽象性,加强实用性.当前,Matlab 软件在线性代数中的应用尤为突出 . Matlab 是由美国 Mathwork 公司于 1984 年正式推出的 ,它是一种科学计算软件 ,主 要适用于矩阵运算及控制和信息处理领域 的分析设计 .它使用方便 ,输入简洁 ,运算 高效,内容丰富 ,并且很容易由用户自行扩 展.Matlab是矩阵实验室的缩写，它是以矩 阵运算为基础的交互式程序语言 ,专门针 对科学 ,工程计算及绘图的需求 .既然 Matlab 是基于矩阵语言发展起来的 计算机 演算纸,其语言表述与数学表述非常类 似.于是,在学习线性代数时 ,可用 Matlab 帮助我们深入理解和掌握基本概念和方 法.因此 ,Mariab 软件已成为国内外高等院 校学生必须掌握的工具了 .下面就特征值 问题 ,试图用 Matlab 来求解一下 . 1 特征值与特征向量 方阵的特征值与特征向量是矩阵问题 中最重要的内容之一 ,此类问题具有重要 的理论价值与广泛的应用背景 .在工程技 术中有许多与振动和稳定性有关的问题 （如:桥梁或建筑物的振动问题 ,机械机件 的振动 , 飞机机翼的颤动 , 电磁振荡中的某 些临界值的确定等 ）,在数学中 ,解微分方 程组及简化矩阵的计算等 ,也会遇到这样 的问题 ,而这些问题的求解往往都可以转 化为求矩阵的特征值和特征向量问题 . 设A是27X力阶矩阵，如果存在数，使 得 A=X, 则称 A 为矩阵 A 的一个特征 值 ,为对应特征值的特征向量 .求解特 征值问题时 ,当方阵阶数较高时 ,手工演算 将变得比较麻烦 ,而 Mat1ab 中的命令计算 特征值和特征向量就很方便 ,并可以得到 不同的子结果和分解 .求解特征值与特征 向量的算法很多 ,你可以先求出特征多项 式再求多项式的根 .在求方阵特征多项式 时调用命令poly，然后用roots函数求得特征 值,但这种算法不能直接得出特征向量 .在 Matlab中还有一个方便而用的命令eig,它 可同时求得特征值与一组特征向量 .具体 调用命令如下 ,我们一般使用下面的第二 条命令: (1)d=eig()%输入矩阵为方阵，仅 输出特征值 ,并以向量形式 d 存放 . ⑵[X,lamda】=e()%参数为一个矩 阵,它的各列为方阵的特征向量，参lamda是 个对角阵 ,其对角元素为的特征值 . (3)[,lamda]=eig(A, nobaJalTCe)% 当矩 阵 A 中有元素与截断误差数量级相差不远 时 ,该命令可能更精确 . 『X,lamda]:e 埋(), no 2 通过一个实例来直观展示 Matlab 在特征 值问题上的应用——生态模型 例:考察栖息在同一地区的兔子和狐 狸的生态模型 ,对两种动物的数量的相互 依存的关系可用以下模型描述 : JXn=IAx.1_Y1,2,…)【=0.1x1+o.85y,l — 其中 Xn,Y 分别表示第 /7 年时免子和狐 狸的数量比例 ,而 ,,表示开始观察时 0)兔子和狐狸的数量比例 ,问 lO 年以 后两种动物的数量比例是多少 ?30 年或 6O 年以后又如何 ? 解 :这个问题可用矩阵乘法来描述 .令 厂, IJ(H=o,1,2,…),若初始状态为： (]:(:], 则一年以后的比例为 al=Aao,从开始观察到第/7年,A始终保持 不变，因此=1=A A, 下面用 MATLAB 程序进行计算 ： gt;gt;A= 【1.1,