07学年论文宫琳

《07学年论文宫琳》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、•■••本科学年论文（设计）题目：对于大规模矩阵矩阵特征值解法初步探索学院专业信息与计算科学班级07级信息与计算科学二班学号20070534213姓名官塑指导教师山东财政学院教务处制二O一一年摘要：现实生活中，很多问题不同系数都可以归结为大规模矩阵问题，这为处理解决问题提供了方法，但是人规模矩阵往往过于复杂，难于计算和处理，这就需要寻找一种可行的分类解决方法。在矩阵解法中，矩阵特征值是矩阵内容的核心内容，本文将对于矩阵的特征值解法进行讨论。对于乘幕法，反幕法，子空间代换法的讨论进行初探，并口对于其中的使用解决领域，精度等进行讨论。

2、关键字：矩阵特征值乘幕法反慕法子空间代换法引言：在科学技术的应用领域屮，许多问题都归为求解一个和应的特征系统。特征系统里面的诸多因索通过矩阵的形式表现出来，如动力学系统和结构系统中的振动问题，求系统的频率与振型：物理学中的某些临界值的确定等等。但是i般系统内部因素诸多，数值很多，大规模矩阵情况很复杂，计算繁琐，需要根据不同类型进行矩阵求解。矩阵解法屮利用矩阵特征值进行求解，是一个比较普遍的方法。这在工程计算领域非常常见，如量子物理中的Kohn-Sham方程求解化为哈密顿矩阵某些关键特征值的对■策；财务风险中的决策树法以财务变虽为判

3、别点建立决策树模型，以鼓低误判成木为标准对样木公司进行分类，进行矩阵运算；会计中所有者权益变动表利用矩阵求解等。本论文希望通过对大规模矩阵的主特征值解决的思想和方法进行探讨，例如乘幕法，反幕法和子空间代换法。通过对不同类型的计算方法的结果和适用类型的对比和探讨，以及对于不同类型解决方法的优化等方血•的研究，加深对矩阵实际的丿应用的理解，学以致用。计算矩阵的时候会引入短阵的核心内容一特征值2,下面为相关介绍。设A为川阶方阵，4=(切)若有数几使Ax=Ax则称兄为A的特征值，x为相皿于兄的特征向量。因此，特征问题的求解包括两方面：1.

4、求特征值入满足0(2)=det(A—2/)=02.求特征向量xw/T(xhO),满足齐方程组(A-AI)x=0称卩(久)为A的特征多项式，它是关于A的n次代数方程。关于特征值的求解，一般用乘幕法，反幕法以及子空间代换法，但是不同方法使用的情况不同，楮度不同，过程不同，因而対于不同的需求，计算机上机实现选择方式也不同，本文将对其进行阐述。关于矩阵特征值的解法（一）乘幕法与反幕法在实际工程应用中，如大型结构的振动系统中，往往要计算振动系统的最低频率（或前几个最低频率）及相应的振型，相应的数学问题便为求解矩阵的按模最人或前几个按模最人特

5、征值及相应的特征向量问题，或称为求主特征值问题。1.1乘慕法介绍：乘幕法是川于求大型稀疏矩阵的主特征值的迭代方法，其特点是公式简单，易于上机实现。基木思想：任取一个非零向量，由已知的矩阵A的乘帚构造一个向量序列，通过不断迭代取得主特征值及相关特征向量。乘幕法的公式推算为：设AM笃取初始向量化R",令x⑴=加°）,x（2）=Ax（l…，一般有(1)形成迭代向量序列｛?k）｝o由递推公式，冇(2)“占）=2）这表明兀⑹是用A的k次幕左乘严得到的，因此称此方法为乘幕法，｛0｝称为迭代序列。例1用规范化乘幕法计算矩阵A的主特征值及相应特

6、征向量14130-4-5-1解：首先，通过矩阵的基本运算，不难得到A的特征值21=6,走=3,23=2其次，运用乘幕法运算，对结果进行对比不妨取初始值x（0）=d,1,1）了，用规范化乘帚法公式计算max（jc（o））=1『(0)_x(o)/max(x(()))=(1,1,l)rx(l)=Ay(0)=(10,8,l)r其它结果见表1（表中的向量均为转置向量）。表1kmax(/))兀⑹=产/max（兀伙））严“严01(1,1,1)(10,&1)110(1,0.&0.1)(7.2,5.4,-0.8)27.2(1,0.75,-0」111

7、11)(6.5,4.75,-1.222222)36.57(1,0.730769,-0.203704)(6.230766,4.499997,-1.407408)46.230766(1,0.722222,-0.225880)(6.111108,4.388886,・1.1451767)56.111108(1,0.718182,-0.237561)(6.054548,4.336336,-1.475122)66.054548(1,0.716216,-0.243639)(6.027024,4.310808,-1.487278)76.02702

8、4(1,0.715247,-0.246768)(6.01345&4.298211,-1.483536)86.013458(1,0.714765,0248366)(6.00671,4.291945,-1.496732)96.00671(1,0.714