矩阵函数以与应用毕业设计_说明
数学与应用数学毕业论文——正交矩阵及其应用
本科生毕业设计(论文)正交矩阵及其应用学院:专业:数学与应用数学学号:学生姓名:指导教师:二〇一一年六月摘要如果n阶实矩阵A满足,那么称A为正交矩阵.正交矩阵是由内积引出的.本文例举了正交矩阵在线性代数、化学和物理中的三个应用.在线性代数中,求标准正交基一般用Schimidt正交化方法.本文论证了一种特殊的正交矩阵——初等旋转矩阵——也可以求任一向量空间的标准正交基,并通过实例说明此方法的应用.在化学上,原子轨道的杂化,实际是由一组相互正交的单位基向量,通过线性变换转化为另一组相互正交的单位基向量.而线性代数中由一组标准正交基到另一组标准正交基的过渡矩阵是正交矩阵,因此可以利用正交矩阵的性质求原子轨道的杂化轨道式.在物理上,任一刚体运动都对应一个正交矩阵,本文证明了曲线作刚体运动时曲率和挠率是两个不变量.关键词:正交矩阵;初等旋转矩阵;标准正交基;原子轨道的杂化;曲率;挠率AbstractOrthogonal matrices and its applicationsIf a-dimensional real matrixsatisfies,we call it orthogonal matrix. Orthogonal matrix is extracted by inner product.This paper enumerats the applications of orthogonal matrix inlinear algebra, chemistry, and physics. Schimidt method is always used to find the standard orthogonal basis in linear algebra. A special kind of orthogonal matrix, namely elementary rotational matrix, is established to find the standard orthogonal basis in this paper. The orbital atom heterozygous is actually made by a team of mutually orthogonal unit basis vector, through linear transformation into another group of mutually orthogonal unit basis in linear algebra. Thetransition matrix of a group of standard orgthogonal basis to another group of standard orthogonal basis is an orthogonal matrix. Therefore, properties of orthogonal basis can be used to find the orbital atom heterozygous. In physics, any rigid motion corresponds with anorthogonal matrix. The curvature and torsion rate are proved to be two invariants when a curve is in rigid motion.Keywords:Orthogonal matrix; Elementary rotation matrix; Standard orthogonal basis; The orbital atom heterozygous; Curvature;Torsion rate目录1.引言 12.正交矩阵的基本知识 32.1正交矩阵的定义与判定 32.2 正交矩阵的性质 33.正交矩阵的应用 53.1 正交矩阵在线性代数中的应用 53.2正交矩阵在化学中的应用 113.3正交矩阵在物理学中的应用 14参考文献 18致谢 19正交矩阵及其应用姓名:学号:班级:1.引言因为行列式要求行数等于列数,排成的表总是正方形的,通过对它的研究又发现了矩阵的理论.矩阵也是由数排成行和列的数表,可以行数和烈数相等也可以不等.矩阵和行列式是两个完全不同的概念,行列式代表着一个数,而矩阵仅仅是一些数的有顺序的摆法.利用矩阵这个工具,可以把线性方程组中的系数组成向量空间中的向量;这样对于一个多元线性方程组的解的情况,以及不同解之间的关系等等一系列理论上的问题,就都可以得到彻底的解决.矩阵是数学中的一个重要的基本概念,是代数学的一个主要研究对象,也是数学研究和应用的一个重要工具.“矩阵”这个词是由西尔维斯特首先使用的,他是为了将数字的矩形阵列区别于行列式而发明了这个述语.而实际上,矩阵这个课题在诞生之前就已经发展的很好了.从行列式的大量工作中明显的表现出来,为了很多目的,不管行列式的值是否与问题有关,方阵本身都可以研究和使用,矩阵的许多基本性质也是在行列式的发展中建立起来的.在逻辑上,矩阵的概念应先于行列式的概念,然而在历史上次序正好相反.凯莱先把矩阵作为一个独立的数学概念提出来,同研究线性变换下的不变量相结合,首先引进矩阵以简化记号并发表了关于这个题目的一系列文章.1858年,他发表了关于这一课题的第一篇论文《矩阵论的研究报告》,系统地阐述了关于矩阵的理论.文中他定义了矩阵的相等、矩阵的运算法则、矩阵的转置以及矩阵的逆等一系列基本概念,指出了矩阵加法的可交换性与可结合性.另外,凯莱还给出了方阵的特征方程和特征根(特征值)以及有关矩阵的一些基本结果.凯莱出生于一个古老而有才能的英国家庭,剑桥大学三一学院大学毕业后留校讲授数学,三年后他转从律师职业,工作卓有成效,并利用业余时间研究数学,发表了大量的数学论文.1855年,埃米特(C.Hermite,1822~1901)证明了别的数学家发现的一些矩阵类的特征根的特殊性质,如现在称为埃米特矩阵的特征根性质等.后来,克莱伯施(A.Clebsch,1831.1872)、布克海姆(A.Buchheim)等证明了对称矩阵的特征根性质.泰伯(H.Taber)引入矩阵的迹的概念并给出了一些有关的结论.在矩阵论的发展史上,弗罗伯纽斯(G.Frobenius,1849~1917)的贡献是不可磨灭的.他讨论了最小多项式问题,引进了矩阵的秩、不变因子和初等因子、正交矩阵、矩阵的相似变换、合同矩阵等概念,以合乎逻辑的形式整理了不变因子和初等因子的理论,并讨论了正交矩阵与合同矩阵的一些重要性质.1854年,约当研究了矩阵化为标准型的问题.1892年,梅茨勒(H.Metzler)引进了矩阵的超越函数概念并将其写成矩阵的幂级数的形式.傅立叶、西尔和庞加莱的著作中还讨论了无限阶矩阵问题,这主要是适用方程发展的需要而开始的.矩阵本身所具有的性质依赖于元素的性质,矩阵由最初作为一种工具经过两个多世纪的发展,现在已成为独立的一门数学分支——矩阵论.而矩阵论又可分为矩阵方程论、矩阵分解论和广义逆矩阵论等矩阵的现代理论.矩阵的应用是多方面的,不仅在数学领域里,而且在化学、力学、物理、科技等方面都十分广泛的应用.本文主要介绍正交矩阵与其应用.我们把阶实数矩阵满足,称为正交矩阵.尽管我们在这里只考虑实数矩阵,但这个定义可用于其元素来自任何域的矩阵.正交矩阵是由内积自然引出的,要看出其与内积的联系,考虑在维实数内积空间中的关于正交基写出的向量.的长度的平方是.如果矩阵形式为的线性变换保持了向量长度,所以有限维线性等距同构,比如旋转、反射和它们的组合,都产生正交矩阵.本文例举了正交矩阵在线性代数、化学和物理中的三大应用.其中,在线性代数中,求标准正交基除了用Schimidt正交化方法外,本文论证了正交矩阵的其中一种矩阵...初等旋转矩阵也可以求任一矩阵的标准正交基,此法用实例与Schimidt 正交化方法对比;在化学上,根据原子轨道的杂化理论,杂化的原子都有其轨道杂化式,对于形成对阵的原子轨道杂化,利用正交矩阵的性质可以求解该原子杂化轨道的杂化轨道式;在物理上,任一刚体运动都对应一个正交矩阵, 三维空间一条曲线经过刚体运动,其曲率和挠率是不变的,本文考察了曲线做刚体运动时的不变量——曲率和挠率.2.正交矩阵的基本知识本节中在没有特别说明的情况下,都表示为正交矩阵,记矩阵的秩为,与为矩阵的第列与第列,表示矩阵的第行.表示行列式的值即=.2.1正交矩阵的定义与判定定义2.1.1[3]阶实数矩阵满足(或,或),则称为正交矩阵.判定2.1.2 矩阵是正交矩阵;判定2.1.3 矩阵是正交矩阵;判定2.1.4 矩阵是正交矩阵;备注:判定一个是方阵是否为正交矩阵往往用定义,即(或,或),也可以验证的行向量或列向量是否是两两正交的单位向量.当已知的正交矩阵求证其他的结论时,要用正交矩阵的定义及有关性质2.2 正交矩阵的性质若是正交矩阵,则有以下性质([3]):性质2.2.5,则可逆,且其逆也为正交矩阵.证明显然.所以也是正交矩阵.性质2.2.6,,也是正交矩阵, 即有:(1)当时,, 即;(2)当时,, 即证明若是正交矩阵,, 由性质2.2.5,为正交矩阵.因为,所以,当时,, 即;当时., 即.从而为正交矩阵.性质2.2.7是正交矩阵.证明因为,所以.因此,也是正交矩阵性质2.2.8是正交矩阵的充分必要条件是.证明必要性若是正交矩阵,则另一方面,一方面,于是,,;充分性因为是正交矩阵,若,显然也是正交矩阵.性质2.2.9 若也是正交矩阵, 则,,,都为正交矩阵.证明由可知,故为正交矩阵.同理推知,,,均为正交矩阵.正交矩阵的性质主要有以上几点, 还有例如它的特征值的模为1, 且属于不同特征值的特征向量相互正交; 如果是它的特征值, 那么也是它的特征值, 另外正交矩阵可以对角化, 即存在复可逆矩阵, 使,其中为的全部特征值, 即. 这些性质证明略.3.正交矩阵的应用3.1 正交矩阵在线性代数中的应用在线性代数中我们通常用施密特方法求标准正交基,现在可以用正交矩阵中的一种殊矩阵求标准正交基---初等旋转矩阵即Givens矩阵.定义3.1[1] 设向量则称阶矩阵为向量下的Givens矩阵或初等旋转矩阵,也可记作.下面给出Givens矩阵的三个性质[2],[10]性质3.1.1 Givens矩阵是正交矩阵.证明由,则,故是正交矩阵.性质3.1.2 设,则有.证明由的定义知,,且,即右乘向量,只改变向量第和第个元素,其他元素不变.性质3.1.3 任意矩阵右乘,只改变的第列和列元素; 任意矩阵左乘,只改变的第行和行元素.证明由性质3.1.2和矩阵乘法易得结论.引理3.1.4[2] 任何阶实非奇异矩阵 ,可通过左连乘初等旋转矩阵化为上三角矩阵, 且其对角线元素除最后一个外都是正的.定理3.1.5[10] 设是阶正交矩阵若, 则可表示成若干个初等旋转矩阵的乘积, 即;若, 则可以表示成若干个初等旋转矩阵的乘积再右乘以矩阵, 即, 其中是初等旋转矩阵.().证明由于是阶正交矩阵,根据引理3.1.4知存在初等旋转矩阵,使(这里是阶上三角阵),而且的主对角线上的元素除最后一个外都是正的,于是(3-11)注意到是正交矩阵,由(3-11)式得,,即(3-12)设=,其中,,则=.由上式得所以, (3-13)即,当时,;当时,.记,注意到是初等旋转矩阵,故定理1结论成立.引理3.1.6[1] 设其中是阶正交矩阵,是阶上三角阵,是零矩阵.定理3.1.7[10] 设,则可以通过左连乘初等旋转矩阵,把变为的形式,其中是阶上三角阵,是矩阵.证明由引理3.1.6知,其中是阶正交矩阵,是阶上三角阵.又根据定理1知:,则是初等旋转矩阵.(I)当时,;(II)当时,,则.显然,是阶上三角阵,当时,与除最后一行对应元素绝对相等、符号相反外,其余元素对应相等.当时时,.综上,知本定理的结论成立.设,,,是欧氏空间的子空间的一组基,记是秩为的的矩阵.若满足定理2的条件,则存在初等旋转矩阵,使(3-14)且所以(3-15)由(3-14)(3-15)两式知,对、做同样的旋转变换,在把化为的同时,就将化成了,而的前个列向量属于子空间.综上所述可得化欧氏空间的子空间的一组基为一组标准正交基的方法:(1)由已知基为列向量构成矩阵;(2)对矩阵施行初等旋转变换,化为,同时就被化为正交矩阵,这里是阶上三角阵;(3)取的前个列向量便可得的一组标准正交基.显然,上述方法是求子空间的一组标准正交基的另一种方法.下面,我们通过实例对比Schimidt正交化求标准正交基.例求以向量,,为基的向量空间的一组标准正交基.解方法一用Schimidt正交化把它们正交化:,,再把每个向量单位化,得,,.即,,,就是由,得到的的一组标准正交基.方法二(利用连乘初等旋转矩阵)设矩阵,对分块矩阵依次左乘,,,=,=,=,得=,则,,取,,.那么就是由,得到的的一组标准正交基.对比两者的解法,用Schimidt正交化把它们正交化需要的是记公式,若向量的维数比较多的,计算比较麻烦,而用初等旋转矩阵则可根据向量组成的矩阵的特点来求其标准正交基.3.2正交矩阵在化学中的应用原子轨道的杂化是在一个原子中不同原子轨道的线性组合.在结构化学原子轨道杂化理论中,原子中能级相近的几个原子轨道可以相互混合,从而产生新的原子轨道.杂化过程的数学表达式为,为新的杂化轨道,为参加杂化的旧轨道,为第个杂化轨道中的第个参加杂化轨道的组合系数[4].在杂化过程中,轨道数是守恒的,并且杂化轨道理论有三条基本原则[5]:(1)杂化轨道的归一性.杂化轨道满足;(2)杂化轨道的正交性.;(3)单位轨道贡献.每个参加杂化的单位轨道,在所有的新杂化轨道中该轨道成分之和必须为一个单位,即=1.由杂化轨道原理,原子轨道的杂化,实际是由一组相互正交的单位基向量,通过线性变换转化成为另一组相互正交的单位基向量.在线性代数中由一组标准正交基到另一组标准正交基的过渡矩阵是正交矩阵,那么原子轨道的杂化,就可以转化为求出正交矩阵,作线性变换的过程.(A)杂化轨道.以甲烷分子的结构为例,激发态碳原子的电子组态为,这样在形成分子时,激发态碳原子的一个2原子轨道和3个原子轨道进行杂化形成4个等同的杂化轨道.设在激发态碳原子中四个能量相近的原子轨道,,,是一组相互正交的基向量,再通过线性变换将它们转化成另一组相互正交的基向量,,,,那么线性变换系数矩阵A必为正交矩阵,即=.A为正交矩阵,分别是,,,在四个坐标轴的分量.在等性杂化中,四个基向量,,,在四个坐标轴上的分量是相等的,即由四个能量相近的原子轨道,,,进行杂化时形成四个等同的杂化轨道,在四个杂化轨道上,原子轨道和成份完全相同.根据这些理论,我们来求正交矩阵A.因为A 是正交矩阵,由定义可得,即,所以,得=(取正值).又因为是等性杂化轨道.有,=1,所以=(取正值).即得到.又因,,,取符合条件的,,.同理,,即,,得,,取,.又,,得,,.所以,.可以写出四个杂化轨道的杂化轨道式为,,.(B)杂化轨道一个和一个原子轨道杂化形成两个杂化轨道.同样,线性变换的系数矩阵是正交矩阵.根据等性杂化理论有,,,于是,,(取正值).又,,故,,即,.所以杂化轨道式为.3.3正交矩阵在物理学中的应用任意刚体运动都对应一个正交矩阵, 三维空间一条曲线经过刚体运动, 其曲率和挠率是不变的, 称它们为运动不变量.首先我们来简单认识曲率和挠率.曲线的曲率就是针对曲线上某个点的切线方向角对弧长的转动率,通过微分来定义,表明曲线偏离直线的程度.曲率越大,表示曲线的弯曲程度越大.(为角变量,为弧长)趋向于0的时候,定义就是曲率.即.而挠率,它的绝对值度量了曲线上邻近两点的次法向量之间的夹角对弧长的变化率.平面曲线是挠率恒为零的曲线.空间曲线如不是落在一平面上,则称为挠曲线,又由于挠率体现了密切平面的扭转状况,通常说它表示了曲线的扭曲程度.曲线在某点的挠率记为,=.下面, 我们来考察曲线作刚体运动时的不变量[6],[9].设曲线与曲线只差一个运动, 从曲线到曲线的变换为(3-21)其中,是三阶正交矩阵,是常数.对(3-21)两边求阶导数,得.从而有. (3-22)因为是正交矩阵, 所以也有. (3-23) 另一方面, 由一阶, 二阶, 三阶导数, 可作成矩阵.两边取行列式, 由,得.现在取可类似地讨论.因为, (3-24), (3-25)(3-22)代入(3-24)的右边,得=++. (3-26)因(3-24)与(3-25)右边相等, 有(3-25)右边与(3-26)式右边相等,得,,.由正交矩阵的性质2.2.6知,且由,将上面三式左右分别平方相加,=++=.写成矢量函数, 即得于是我们可推得,.这里的分别是曲线的曲率与挠率.参考文献:[1] 陈景良,陈向晖.《特殊矩阵》.第一版.清华大学出版社,2001:353-360[2] 程云鹏.《矩阵论》.第二版.西北工业大学出版社,1999:94.99,196-215[3] 王萼芳,石生明.《高等代数》.第三版.北京:高等教育出版设,2007:162-392[4] 周公度,段连运.《机构化学基础》.第4版.北京大学出版社,2009:79-187[4] 王立东主编《数学》.第一版.大连理工大学出版社,2008:63-74[5] 赵成大等《物质结构》.人民教育出版社. 1982:219-226[6] 强元棨,程嫁夫.《力学》上册.第一版.中国科学技术大学出版社:2005:332-53[7] 张焕玲等《一种求欧氏空间子空间的标准正交基的新方法》山东大学.1996.3.9卷(1)期:14-16[8] 刘钊南.《正交矩阵的作用》.湘潭师范学院学报.1987.11.16: 3[9] 陈少白.《空间曲线的刚体运动基不变量》. 武汉科技大学学报.2003.12.26卷(4)期:424-426[10] 刘国志.《欧氏空间子空间的标准正交基的全新方法—Givens变换法》.抚顺石油学院学报.1996.3.16卷(1)期:78-81致谢感谢父母,给了我生命,也让我懂得这世上什么是真情!当我们遇到困难的时候,会倾注所有一切来帮助我们的人是父母;当我们受到委屈的时候,能耐心听我们哭诉的人是父母.当我们犯错误时,能够毫不犹豫地原谅我们的人是父母;当我们取得成功的时候,会衷心为我们庆祝与我们分享成功的喜悦的,仍然是父母;而现在我们远在外地学习,依然牵挂着我们还是父母.感谢父母给予我爱,是您们让我感到骄傲与自豪!感谢老师,授予我知识!大学四年,不少老师给予我无微不至的关怀,这将成为我人生中难以忘怀的回忆.我不仅从您们身上学到许多专业知识,更多的是学到了为人处世的道理.在和您们的交流中,我对我的未来有了更好的规划.您们是我人生的航标,让我在迷茫时找到前进的方向;您们是我精神上的支柱,让我在困难时重新振作.大学四年,如果没有您们的博学知识,没有您们的倾注爱心,没有您们的谆谆善诱,我将不可能收获那么多.假如我能搏击蓝天,那是您们给了我腾飞的翅膀;假如我是击浪的勇士,那是您们给了我弄潮的力量;假如我是不灭的火炬,那是您们给了我青春的光亮!感谢帮助过我、教导过我的老师们,是您们,让我懂得给予与付出才是最重要的,是您们,让我明白做人就要不断进取,迎难而上,力争上游!本毕业论文是在我的导师XX的亲切关怀和悉心指导下完成的,她给我的论文提出了不少宝贵的意见;她严肃的科学态度,严谨的治学精神,精益求精的工作作风,深深地感染和激励着我.从课题的选择到项目的最终完成,XX老师都始终给予我细心的指导和不懈的支持,在此谨向XX老师致以诚挚的谢意和崇高的敬意.。
正定矩阵及其应用
本科毕业论文(设计)正定矩阵及其应用学生:学号:专业:指导老师:答辩时间:装订时间:A Graduation Thesis(Project)Submitted to School of Science,Hubei University for Nationalities In Partial Fulfillment of the Requiring for BS DegreeIn the Year of 2016Positive definite matrices and their applicationsStudent Name: Student No.:Specialty:s Supervisor:Date of Thesis Defense: Date of Bookbinding:摘要矩阵是高等代数里的一个基本概念,是代数知识的基础,是矩阵代数的一个主要研究对象. 它不仅是数学的一个重要分支,而且已经成为现在科技领域处理有限维空间形式与数量关系的强有力的工具. 而正定矩阵是从矩阵延伸出来的具有特殊性质的矩阵,是研究二次型的基础,在函数、不等式中都有应用,因此正定矩阵的特殊性质和广泛应用得到了许多学者关注,进而对此进行了大量的研究. 本文从矩阵最基本的概念和性质出发,由浅入深,层层递进. 从矩阵的性质出发,给出了正定矩阵定义及其等价定义,归纳整理了正定矩阵的性质及其部分证明,总结了正定矩阵的判定定理,最后研究正定矩阵在理论证明和在函数极值中的应用.关键词:矩阵正定二次型正定矩阵极值AbstractThe matrix is very important in advanced algebra. It is not only an important branch, but also have become a powerful tool for studying finite dimensional space and quantity r- elationship in the real of modern science and technology. However , extending from the m- atrices, the positive definite matrix is a special matrix, which is a foundation for studying quadratic form and apply properly to both functions and inequality. Thus, its special prop- erty and wide applications have drawn scholars'attention, and a lot of research have been done. This paper begins with the matrix'primary concept and properties, going from the e- asy to the difficult. We define the positive definite matrix and its equivalent one, the sum up its properties and partial evidence, and summarize the determined theorems. At last, we study its application in theory and the solution of the function extremum.Keywords: matrix,positive definite quadratic,positive definite matrix,extremum目录摘要IAbstractII1绪论11.1 课题背景11.2 课题研究的目的和意义11.3 国外研究概况22 预备知识32.1 矩阵32.2二次型53正定矩阵83.1正定二次型83.2正定矩阵的判定定理94正定矩阵的应用134.1正定矩阵的相关命题134.2正定矩阵在函数极值中的应用15总结与展望18致201绪论我们知道矩阵是高等代数中非常重要的容之一. 在学习高等代数时,矩阵方面的知识也经常被用到. 而正定矩阵又是矩阵中的重点,它不单单用来解决数学中的问题,还应用于许多的科学领域. 本课题阐述了正定矩阵研究背景、正定矩阵的研究的目的和意义、正定矩阵的现状以及发展方向,明确指出了研究正定矩阵应用所面临的问题.1.1 课题背景正定矩阵作为一类常用矩阵,对它的研究最早出现在二次型中. 它也是从正定二次型中抽象出来的一个概念,有了正定矩阵的概念后,解决二次型的问题就变得简单方便. 不仅在代数学中应用广泛,在函数学、几何学、图像处理学、概率统计和物理学等学科中都得到了广泛的应用. 因此它的性质、定理以及应用问题一直备受学者关注. 而在实际生活问题中也经常出现一些相关数学问题,而用正定矩阵解决问题可能会更方便简洁一点. 这就需要我们研究正定矩阵的应用,如正定矩阵在四则运算、在函数极值、在不等式中的应用. 因此可以使得我们可以更好地使用正定矩阵这一重要工具. 本文通过对正定矩阵的理解和掌握,查阅各种相关资料,对正定矩阵及其相关知识点进行归纳总结,并且由此给出了正定矩阵在四则运算和函数极值及中的应用.根据课题研究容和手中相关文献资料,了解课题研究现状,学习掌握相关理论基础知识,并进行初步研究,撰写开题报告.1.2 课题研究的目的和意义矩阵是代数中一个非常重要的概念,是研究和解决数学问题的一个重要工具. 而正定矩阵是一类非常重要的矩阵,在矩阵中扮演着重要的角色,因此是我们学习矩阵时不可忽略的重点. 本文对我们对数学感兴趣的学生深入理解和掌握正定矩阵理论有非常重要的意义. 能够加强我们对正定矩阵的掌握,也可以促进正定矩阵理论的进一步完善,丰富正定矩阵的应用,加强我们对正定矩阵的理解,丰富矩阵的理论知识. 有助于我们对整个高等代数知识的一体化的认识. 从而可以培养我们对代数知识的串联思想. 正定矩阵多方面的应用,能够开阔我们的视野,加强我们的联想能力,引起我们对数学的探究欲望,对知识的渴望.研究矩阵的正定性,在代数理论和应用中具有重要意义. 正定矩阵不仅在数学方面,在其他各个领域都具有广泛的应用价值,因此引起了学者们极大的研究兴趣. 这些研究不断丰富了正定矩阵的理论知识,也引起了我们对正定矩阵的兴趣.1.3 国外研究概况随着数学的影响力越来越大,矩阵对数学的研究也显得越来越重要. 在代数方面,正定矩阵也同样占有非常重要的地位. 因此人们对正定矩阵的研究也越来越广泛. 因而对正定矩阵的理解和应用也越来越深入,其应用围也越来越广泛. 在函数学、几何学、经济学、图像处理学、概率统计和物理学等学科中都得到了广泛的应用.在历史上,正定矩阵的相关研究最早出现在二次型和Hermite型中. 但是当时对于的正定矩阵局限于对实对称矩阵或者Hermite矩阵. 1970年,Johnson引入了不再局限于对实对称矩阵或者Hermite矩阵实对称矩阵的概念. 他给出了正定矩阵较为广义的定义. 1985年,炯生也给出了正定矩阵较为广义的定义. 1984年,佟文廷再次将正定矩阵的定义进行了推广. 他给出了推广正定矩阵的各种定义. 1988年,夏长富将实对称矩阵的正定性做了深入推广. 他又进一步极大的丰富了正定矩阵的理论. 1990年,屠伯埙将各类广义正定矩阵进行深度结合. 他重新定义了广义正定矩阵,将它称之为亚正定矩阵. 在研究正定矩阵的过程中,许多学者取得了惊人的理论成果,其成果也得到了广泛的应用. 除了对正定矩阵的研究,许多学者还对正定矩阵相关容进行了研究,同样取得了巨大的成就. 近年来,在完善正定矩阵理论成果的历史中,得出了许多其他的概念和定理,将各类正定阵统一起来. 这些新的研究成果对完善正定矩阵的理论和其应用具有非常大的价值.虽然对正定矩阵的研究这么广泛,但是这些正定矩阵的研究只局限在正定矩阵的理论分析方面. 它的一些实际方面的应用还有待笔者和一些学者去探索挖掘.2 预备知识2.1 矩阵定义2.1.1 由n m ⨯个数),2,1,,2,1(n j m i a ij ==;排成的m 行n 列的数表⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=mn m m n n a a a a a aa a a A 212222111211, 称为n m ⨯矩阵,记作.n m ij a A ⨯=)( 特殊地,当n m =时,矩阵称为方阵.定义2.1.2 把一矩阵A 的行列互换,所得到的矩阵称为A 的转置. 记为T A (或者记为'A ).即, 设⎪⎪⎪⎪⎪⎭⎫⎝⎛=sn s s n n a a a a a a a a a A 212222111211,所谓A 的地转置就是指矩阵.212222111211⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=sn n n s s T a a a a a a a a a A 显然,n s ⨯矩阵的转置是s n ⨯矩阵,即n s ij a A ⨯=)(,则.s n ij a A ⨯=)( 转置矩阵满足以下运算规律()()()().T T TT T T T T TT kA kA A B AB B A B A A A ==+=+=,,,定义2.1.3 数域P 上的n n ⨯矩阵A 称为对称矩阵,如果T A A =.即若⎪⎪⎪⎪⎪⎭⎫⎝⎛=nn n n n n a a a a a a a a a A 212222111211, 且满足=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛nn n n n n a a a a a a a a a 212222111211⎪⎪⎪⎪⎪⎭⎫⎝⎛nn n nn n a a a a a a a a a212221212111,则称A 为对称矩阵.定理2.1.1 任意一个n 阶实对称矩阵A ,都存在一个n 阶正交矩阵T ,使得AT T T 成对角型. 对角线上的元素为矩阵A 的特征根.定义2.1.4 数域P 上的n n ⨯矩阵A 称为非退化的,如果0≠A ;否则称为退化的. 即,若.0212222111211≠=nnn n nna a a a a a a a a A则A 为非退化的.定义2.1.5N 级方阵A 称为可逆的,如果有n 级方阵B ,使得.E BA AB == (1)这里E 是n 级单位阵.如果矩阵B 适合(1),那么B 称为A 的逆矩阵,记为T A . 注1:只有方阵才可能可逆; 注2:非零的矩阵不一定可逆;注3:若A 可逆,则(1)中的B 必唯一; 注4:若AC AB =,且A 可逆,则C B =.设A 是n 阶可逆矩阵,下列结论成立:()()()()()()()()().5);(4;13;2;111111111-*-------=====n TT AA k kA kA AA A A A A 为非零数定理2.1.2 矩阵A 是可逆的充分必要条件是A 是非退化的.定义2.1.6 数域P 上n n ⨯矩阵B A ,称为合同的,如果有数域P 上可逆的n n ⨯矩阵C ,使.AC C B T =合同是矩阵之间的的一个关系,不难看出,合同关系具有: (1) 反身性:;AE E A T =(2) 对称性:由T C B =即得();11--=BC C A T(3) 传递性: 由111AC C A T =和2122C A C A T= 即得()().21212C C A C C A T=定义2.1.7 设n n y y y x x x ,,,,,,2121 ;是两组文字,系数在数域P 中的一组关系式⎪⎪⎩⎪⎪⎨⎧++=+++=+++=.22112222121212121111n nn n n n nn n n y c y c y c x y c y c y c x y c y c y c x ,,(2) 称为由n x x x ,,,21 到n y y y ,,,21 的一个线性替换,或者简称线性替换,如果系数行列式,0≠cij那么线性替换(2)称为非退化的.2.2二次型定义2.2.1设P 是一数域. 一个系数在数域P 中的n x x x ,,,21 的二次齐次多项式() ++=211221113212,,x x a x a x x x f+++++n n n x x a x a x x a 22222212122+2nnn x a + (3) 称为数域P 上的一个n 元二次型,或者,在不致引起混淆时简称二次型. 令ji ij a a =,.j i >由于i j j i x x x x =,所以二次型(3)可以写成()n n x x a x x a x a x x x f 1121122111321,,+++=n n x x a x a x x a 2222221221+++++22211n nn n n n n x a x x a x x a ++++j i n i nj ij x x a ∑∑-==11(5)把(5)的系数排成一个n n ⨯矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛=nn n n n n a a a a a a a a a A 212222111211, 它就称为二次型(5)的矩阵.令.21⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=n x x x X于是,二次型可以用矩阵的乘积表示出来()⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫⎝⎛=n nn n n n n n Tx x x a a a a a a a a a x x x AX X 2121222211121121,,..),,(21AX X x x x f T n =定理2.2.1 在数域P 上,任意一个对称矩阵都合同于一对角阵. 即,对于任意一个对称矩阵A 都可以找到一个可逆矩阵C 使AC C T成对角矩阵.定义2.2.2二次型),,,(21n x x x f 经过非退化的线性替换所变成的平方和形式称为二次型),,,(21n x x x f 的一个标准形. 即222221121),,,(n n n x d x d x d x x x f ++=为二次型),,,(21n x x x f 的标准形.定义2.2.3 任意一个复系数的二次型,经过一适当的非退化线性替换可以变成规形. 且规形是唯一的.即任一复数的对称矩阵合同于⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛0011的对角阵.定义 2.2.4 实二次型),,,(21n x x x f 经过某一个非线性替换,可使),,,(21n x x x f 变成标准形22112211r r p p p p y d y d y d y d --++++ ,再做一次非退化线性替换就变成221221r p p z z z z --+++ ,称为实二次型),,,(21n x x x f 的规形.3正定矩阵在二次型中,正定二次型占有特殊的地位. 作为本章的开始,我们给出了它的定义,引出正定矩阵的定义. 正定矩阵同样占有非常特殊的地位,我们给出了正定矩阵的判定定理.3.1正定二次型定义3.1.1 在实二次型 ),,,(21n x x x f 的标准型形中,正平方项的个数p 称为),,,(21n x x x f 的正惯性指数;负平方项的个数p r -称为),,,(21n x x x f 的负惯性指数;它们的差()r p p r p -=--2称为),,,(21n x x x f 的符号差.定义3.1.2 实二次型),,,(21n x x x f 称为正定的. 如果对于任意一组不全为零的实数 n c c c ,,,21 都有 0),,,(21>n c c c f .定理3.1.1 n 元实二次型 ),,,(21n x x x f 是正定的充分必要条件是它的正惯性指数等于n .推论3.1.1 正定矩阵的行列式大于零.定义3.1.3 在n 阶矩阵中任选k 行,再取相同行号的列,所选取的行列交汇处的2k 个元素组成的新的矩阵称为n 阶矩阵的一个k 阶主子式.定义3.1.4 子式),2,1(212222111211n i a a a a a aa a a P ii i i i i i=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=, 称为矩阵()n n ij a A ⨯=的顺序主子式.定理3.1.2实二次型AX X x x x x x f T ni nj j i ij n a ==∑∑==1121),,,(是正定的充分必要条件为:矩阵A 的顺序主子式全大于零.定义3.1.5 若对于方阵A 存在一个非零向量X 和实数λ,使得X AX λ=成立. 则称λ为矩阵A 的特征值,X 称为A 相对于λ的特征向量.定义3.1.6 设实二次型AX X x x x f T n =),,,(21 (A 为对称矩阵). 如果对于任意的0X 21≠⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=n x x x ,有0),,,(21>=AX X x x x f T n ,则称该二次型为正定二次型. 矩阵A为正定矩阵.注:本文所讨论的都为实正定矩阵.3.2正定矩阵的判定定理定理3.2.1实对称矩阵A 为正定矩阵的充分必要条件是:存在可逆矩阵P ,使得P P A T =.证明 必要性 因为矩阵A 为正定矩阵,所以矩阵A 合同于单位矩阵,即存在可逆矩阵Q ,使得E AQ Q T =, 即()()()1111----==Q Q EQ Q A TT ,若我们记1-=QP ,则有.P P A T =充分性 设存在可逆矩阵P 使得P P A T =,则对任意()0,,,x 21≠=Tn x x x , 有()()PX PX PX P X AX X TT T T ==,若我们记()Tn y y y PX Y ,,,21 ==. 则22221n T T y y y Y Y AX X +++== ,所以矩阵A 为正定矩阵.定理3.2.2实对称矩阵A 为正定矩阵的充分必要条件是:存在可逆矩阵P ,使得n T E AP P =.证明 充分性 因为矩阵A 为正定矩阵,所以矩阵A 对应的是正定二次型. 因此可以经过非退化线性替换PY X =. 其中()Tn y y y Y ,,,21 =. 使得()()()).,,,(),,,(212121n n T T TT n y y y g a a a Y AP P Y PY A PY AX X x x x f=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛====所以有n E AP P ='.必要性 存在可逆矩阵P 使得 n T E AP P =,则其对应的二次型).,,,()()(),,,(2121n T T T T n x x x f PY A PY APY P Y EY Y y y y g ==== 因),,(21n x x x g 为正定二次型,所以),,,(21n x x x f 也为正定二次型. 所以其对应的矩阵A 为正定矩阵.定理3.2.3实对称矩阵A 为正定矩阵的充分必要条件是:矩阵A 的正惯性指数n p =.证明 充分性 因为矩阵A 为正定矩阵. 由定理3.2.2知矩阵A 合同于单位阵E . 所以矩阵A 的正惯性指数为n .必要性 因为矩阵A 的正惯性指数为n ,由定理3.1.1知矩阵A 对应的二次型为正定二次型. 因此矩阵A 为正定矩阵.定理3.2.4实对称矩阵A 为正定矩阵的充分必要条件是:矩阵A 的所有顺序主子式都大于零.证明 充分性 因为矩阵A 为正定矩阵,所以矩阵A 对应的二次型为正定二次型. 则构造函数)00)(,,,(21≤<k x x x f k 也为正定二次型. 所以其对应的矩阵顺序主子式k A 为正定矩阵,即0>k A . 所以正定矩阵A 的所有顺序主子式都大于零. 必要性 因为矩阵A 的所有顺序主子式都大于零,所以矩阵A 的任一顺序主子式k A 对应的二次函数都为正定二次型. 因此当n k =时对应的二次型),,,(21n x x x f 为正定二次型. 即对应的矩阵A 为正定矩阵.例3.2.1 设二次型323121232221214-2224),,,(x x x x x x x x x x x x f n -+++=λ ,求λ的取什么围,使得),,,(21n x x x f 为二次型.解 二次型),,,(21n x x x f 的矩阵为⎪⎪⎪⎭⎫ ⎝⎛=22-1-2-41-1λλA .由定理3.2.4得,011>=A()(),02244122>+-=-==λλλλλA 得.22<<-λ(),02-24222-12-41123>-=+-=--=λλλλλλA 得.20<<λ 综合可知当20<<λ时,),,,(21n x x x f 正定.定理3.2.5实对称矩阵A 为正定矩阵的充分必要条件是:矩阵A 的所有主子式都大于零.证明 设正定矩阵()n n ij a A ⨯=,则它的任一m 阶主子式为()mm m mk k k k k k k k m a a a a A1111=.作二次型AX X T 和().Y A Y m T 对任意()0,10≠=m k k b b Y ,都有().0,,10≠=n c c X 其中⎩⎨⎧==其它时当,0,,,,21m i i k k k i b c ,由于AX X T 正定,所以000>AX X T . 从而().0000Y A Y AX X m TT= 由0Y 的任意性即证()Y A Y m T 是正定二次型,即().0>m A例3.2.2 判断⎪⎪⎪⎭⎫⎝⎛=1-10121011A 是否为正定矩阵.解 我们直接可以看出矩阵A 的主子式不全大于零.定理3.2.6实对称矩阵A 为正定矩阵的充分必要条件是:矩阵A 的所有特征值都大于零.证明 由定理2.1.1知对于对称矩阵A 存在一个n 阶正交矩阵T . 使得AT T T 成对角型. 对角线上的元素为矩阵A 的特征根.充分性 因为矩阵A 为正定矩阵,所以存在正交矩阵P ,满足⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=n T a a a AP P21. 其中n a a a ,,,21 是矩阵A 的全部特征值. 则矩阵A 对应的二次型为AX X x x x f T n =),,,(21 . 令PY X =,则有()()().,,,)(),,,(2121n T T TT n y y y g Y AP P Y PY A PY AX X x x x f ====又因为矩阵A 为正定矩阵,所以二次型为正定二次型. 因此矩阵A 的特征值全部大于零.必要性 因为矩阵A 的特征值),,2,1(n i a i =都是大于零,所以存在正交矩阵P ,满足⎪⎪⎪⎪⎪⎭⎫⎝⎛=n T a a a AP P21. 则矩阵A 所对应的二次型 ()()()),,,(,,,),,,(2121212121n TT T n n n n x x x f PY A PY APY P Y y y y a a a y y y y y y g===⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫⎝⎛=所以二次型()n y y y g ,,,21 是正定二次型. 因此矩阵A 为正定矩阵.4正定矩阵的应用正定矩阵作为本论文的中心容,我们不仅仅只是研究它的定义和性质,它的应用也是我们需要研究的反向. 正定矩阵的应用非常广泛,它在函数学、几何学、图像处理学、概率统计和物理学等学科中都得到了广泛的应用. 本论文主要研究了它在理论证明中和在函数极值中的应用.4.1正定矩阵的相关命题命题4.1.1 若矩阵B A ,是n 阶正定矩阵,则矩阵B A +也是正定矩阵. 证明 因为矩阵B A ,为正定矩阵,所以对所有00,0>>≠BX X AX X X T T ,. 因此0)(>+X B A X T .命题4.1.2 若矩阵A 是n 阶正定矩阵,R k ∈<0,则kA 也为正定矩阵. 证明 因为所有0,0>≠AX X X T ,所以0)()(>=AX X k X kA X T T .命题4.1.3 若矩阵B A ,都是n 阶正定阵,BA AB =,则AB 也是正定阵. 证明 因为BA AB =,所以()AB BA A B AB T T T===. 所以AB 是对称矩阵又因为B A ,为正定矩阵,所以存在可逆矩阵Q P ,,使得.,Q Q B P P A T T == 因此Q PQ P AB T T =.又因为()()TT T PQ PQ PQ QP QABQ ==-1正定, 且与AB 相似,所以AB 正定.命题4.1.4 设矩阵A 是正定阵,则*1A A ,-为正定阵.证明 因为矩阵A 为正定矩阵. 所以存在可逆矩阵C ,使得C C A T =. 因此()()()TTT C C C C CC A 111111------===. 所以1-A 正定.又因为01*>⋅=-A A A A ,此时,所以*A 也是正定阵.命题4.1.5 设矩阵A 为正定阵,则与矩阵A 合同的矩阵也是正定阵. 证明 因为正定矩阵A 合同于单位矩阵E ,又因为合同矩阵具有传递性 所以结论成立.命题4.1.6 若矩阵A 为正定矩阵,那么矩阵A 的绝对值最大的元素一定在矩阵A 的主对角线上.证明 设{}00,max 00j i a a ij j i ==,000000000≤j j i j j i i i a a a a . 这与矩阵A 为正定矩阵矛盾.例4.1.1判断矩阵⎪⎪⎪⎭⎫ ⎝⎛=113121311B 是不是正定矩阵.解 因为绝对值最大的元素不在主对角线上,所以矩阵B 不是正定矩阵.- . -4.2正定矩阵在函数极值中的应用定义4.2.1 设n 元函数()),,,(21n x x x f x f =在n n R x x x x ∈⎪⎪⎪⎪⎪⎭⎫ ⎝⎛= 210的某个邻域存在一阶和二阶连续偏导数.记⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂=∇n x x f x x f x x f x f )(,,)(,)()(02010 .)(x f ∇称为函数)(x f 在点⎪⎪⎪⎪⎪⎭⎫⎝⎛=n x x x x 210处的梯度,或记为).(0x gradf定义4.2.2设n 元函数()),,,(21n x x x f x f =有二阶连续偏导数,并且在⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=n αααα 21处的一阶偏导全部为零. 则称α为()),,,(21n x x x f x f =的一个驻点,则n 阶矩阵⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=n n n n n n x x x x x x x x x x xx x x x x x x f f f f f f f f f x H212221212121)(. 称为()),,,(21n x x x f x f =在α点的黑塞矩阵.定理4.2.1 设函数),,,()(21n x x x f x f =的一阶和二阶连续偏导数存在.并且在),,,(21n αααα =处的一阶偏导为零. 则由函数二阶偏导所确定的n 元黑塞矩阵⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=n n n n n n x x x x x x x x x x xx x x x x x x f f f f f f f f f x H212221212121)(,满足 (1)当)(αH 为正定矩阵时,()),,,(21n x x x f x f =在α处取得极小值; (2)当)(αH 为负定矩阵时,()),,,(21n x x x f x f =在α处取得极大值; (3)当)(αH 为不定矩阵时,()),,,(21n x x x f x f =在α无极值.- . -证明 因为()),,,(21n x x x f x f =在α的所有二阶偏导数都存在,所以由泰勒公式得()n n x x x f ∆+∆+∆+ααα,,,2211()+=n f ααα,,,21 ()n n n n x x x f x x x x x x ∆+∆+∆+⎪⎪⎭⎫⎝⎛∂∂∆++∂∂∆+∂∂∆θαθαθα,,,22112211 ())10(,,,21221122211<<∆+∆+∆+⎪⎪⎭⎫⎝⎛∂∂∆++∂∂∆+∂∂∆+θθαθαθα其中!n n n n x x x f x x x x x x 又因为()),,,(21n x x x f x f =在α处的一阶偏导为零,所以()n n n i x i x x x f x f i∆+∆+∆+∆=∆∑=θαθαθα,,,(!212211122()).,,,2221111n n ni x x jni j ix x x f xx j i ∆+∆+∆+∆∆+∑∑=+=θαθαθα所以我们可以得到()n n x x x x x f j i ∆+∆+∆+θαθαθα,,,2211 ()).,2,1,(,,,21n j i c f ij n x x j i =+=ααα 当()0,,,21→∆∆∆=∆n x x x x 时,.0→ij c 所以()n n i x i if x f ααα,,,(!2121122 ∑=∆=∆()n n i x x j ni j i j i f x x ααα,,,22111∑∑=+=∆∆+)211122∑∑∑=+==∆∆+∆+ni jni j iijn i iix x c x c .因为()0,,,21→∆∆∆=∆n x x x x 时.0→ij c 所以存在x 的一个领域,使得在这个区域f ∆的符号与()n n i x i if x f ααα,,,2112'2 ∑=∆=()n ni x x jni j ij i f xx ααα,,,22111∑∑=+=∆∆+的符号一致.所以由实二次型及正定矩阵的定义可以证明以上定理的正确性.例4.2.1 求三元函数()32123222132162432,,x x x x x x x x x f -+-++=的极值. 解 求驻点⎪⎩⎪⎨⎧=-==⇒⎪⎩⎪⎨⎧=-==+==-=.111.066022044321321321x x x x f x f x f x x x ,,,,所以驻点为()1,1-1,α. 求得二阶偏导分别为.0,6,0,2,0,4313332222111======x x x x x x x x x x x x f f f f f f所以矩阵()⎪⎪⎪⎭⎫ ⎝⎛=600020004αH , 由以上判定定理可知H 为正定矩阵.所以),,(321x x x f 在()1,1-1,α处取得极小值,极小值为()().51,1,1-=-=f f α例4.2.2 求三元函数32123223132126),,(x x x x x x x x x f -+++=的极值. 解 求驻点⎪⎩⎪⎨⎧===⎪⎩⎪⎨⎧===⇒⎪⎩⎪⎨⎧=-==+==+=.1,27,9.1,0,0.022062063321321312221321x x x x x x x f x x f x x f x x x 或,,所以驻点为()1,001,α,()1,2792,α. 求得二阶偏导分别为.0,0,2,0,0,2,6,6,61331332332221221111=========x x x x x x x x x x x x x x x x x x f f f f f f f f x f所以矩阵()⎪⎪⎪⎭⎫ ⎝⎛=2000260661x H α.所以矩阵().2000260601⎪⎪⎪⎭⎫ ⎝⎛=αH在()1,001,α处的顺序主子式为 .7220002606036266000321-==-====H H H ,,由定理3.2.4知矩阵()1αH 不是正定矩阵,所以()1,001,α不是),,(321x x x f 的极值点. ().20002606542⎪⎪⎪⎭⎫⎝⎛=αH在()1,2792,α处的顺序主子式为 .0144200026065407226654054321>==>==>=H H H ,,由定理3.2.4知矩阵()2αH 是正定矩阵,在()1,2792,α处取得极小值,极小值为 ()().29151,27,92==f f α总结与展望正定矩阵在高等代数中有很多重要的应用,其实质就是简化二次型的运算. 本文一共有四章. 第一章主要介绍了本文的研究背景和现状;第二章归纳了部分矩阵知识和二次型知识;第三章通过正定二次型导出正定矩阵的定义,并且整理了正定矩阵的相关知识,着重归纳证明了正定矩阵的六个判定定理及其证明;第四章在前面两部分的知识基础上,给出了正定矩阵的六个命题及其证明,给出了解决了函数极值存在问题的方法,即正定矩阵在函数极值中的应用. 从代数方面解决分析问题,使我意识到数学的跨度非常大,我们应该加强自己的逻辑思维和联想能力并且要学会多方面思考问题.以上这些正定矩阵的研究只局限在正定矩阵的理论分析方面,也可能总结不太完整,归纳的不够完善,这就希望其它研究者完善,还有它的一些实际方面的应用还有待笔者和一些学者去探索挖掘. 本文作者知识和写作水平有限,不足之处请读者和专家批评指正.致在论文完成之际,我首先要向我的指导老师老师表示最真挚的意,本论文是在导师老师的悉心指导下完成的.在论文写作期间,老师一边要兼顾自己的学业一边还耐心认真地指导我的论文,不辞辛苦,花费了许多宝贵时间和心血. 导师渊博的学识,宽厚待人的学者风,严谨求学的治学态度,忘我的敬业精神让我受益匪浅. 能够师从先平老师,是我的幸运,更是我的荣幸.衷心感和我是同一个指导老师的付江林同学. 感他帮助我指正和修改我论文的不足之处. 因为他的帮助我才能顺利完成我的论文.感我的室友们,感他们的督促与各方面的帮助.还有感我的家人们,没有他们的支持,我的论文不可能顺顺利利的完成.最后,向评阅论文和参加论文答辩的老师们表示由衷的感.由于我知识水平的限制,再加上我写此论文的时间仓促. 文中难免有错误和有待改进之处. 真诚欢迎各位老师、同学提出宝贵意见.参考文献[1]慕生. 高等代数[M]. 复旦大学, 2007.9.[2]王蕚芳. 石生明. 高等代数[M]. . 高等教育, 2003.9.[3]岳贵鑫. 正定矩阵及其应用[J]. 省交通高等专科学校学报, 2008, 10(5):31-33.[4]周杰. 矩阵分析及应用[M]. 大学, 2009.7.[5]王松江. 矩阵不等式[M]. 科学. 科学, 2006.5.[6]文杰. 静. 多元函数的极值问题[J]. 工业大学学报:自然科学版, 2004, 24(1):27-30.[7]邵东南. 马鸿. 正定矩阵的性质及应用.大学学报:自然科学版(2), 1999, 59-62.[8]黄云美. 正定矩阵的性质及其应用.职业学院学报, 17.3(2011).[9]王昊. 正定矩阵的性质及应用[J]. 城市建设理论研究:电子版, 2011(20):59-62.[10]路红军. 一类正定矩阵的性质及其应用[J]. 工学院学报, 2003, 12(3):6-7.[11]史秀英. 正定矩阵的等价命题及其应用[J]. 学院学报:自然科学版,2000(2):44-47.[12]Roger A.Horn . Charles R.Johnson. Matrix Analysis (Second Edition), 剑桥大学,20012.11.[13]宋国际. 论正定矩阵在多元函数极值问题中的应用[J]. 旅游职业学院学报,2010,15(1):58-60.[14]丹. 庆平. 正定矩阵的性质及相关问题[J]. 数学理论与应用, 2011(4):124-128.[15]庭骥. 判别正定矩阵的充分必要条件及其等价性[J]. 矿业学院学报, 1992,20(2):107-110.[16]建忠. 关于正定矩阵的一个不等式及应用[J]. Studies in College Mathematics, 2001,4(1):40-40.[17]贤淑. 莉军. 正定矩阵的降阶判别法及其应用[J]. 印刷学院学报, 2000(2):40-43.[18]家昶. 齐远伟. 条件正定矩阵及其在多元插值计算中的应用[J]. 计算数学, 1989,11(4):386-393.[19]景晓. 正定矩阵基与正交矩阵基及其应用[J]. 师大学学报:自然科学版,2003,18(3):15-17.[20]定华. 矩阵的次特征值及其应用[J]. 中国科学技术大学学报, 2012,42(11):920-924.独创性声明本人声明所呈交的论文(设计)是我个人在导师指导下进行的研究工作及取得的研究成果。
矩阵的函数范文范文
矩阵的函数范文范文矩阵是线性代数中非常重要的概念之一,它在许多数学和科学领域中都有广泛的应用。
矩阵的函数则是将矩阵作为输入,输出一个新的矩阵的函数。
在本文中,我将探讨矩阵的函数及其应用。
首先,我们来看一个简单的例子。
考虑一个2x2的矩阵A,其中元素为a,b,c,d。
我们可以定义一个函数f,使得f(A)=B,其中B也是一个2x2的矩阵。
例如,我们可以定义函数f(A)=A^2,即将矩阵A与自身相乘得到一个新的矩阵B。
在这个例子中,函数f(A)可以表示为:f(A)=A^2=[a,b]*[a,b][c,d][c,d]= [a^2 + bc, ab + bd][ca + cd, cb + d^2]这个例子说明了矩阵的函数可以通过简单的数学运算来定义,而且可以应用到许多实际问题中。
矩阵的函数有许多重要的应用。
其中一个应用是在物理学中的量子力学。
在量子力学中,一个物理系统可以用一个向量表示,而一个操作符则可以用一个矩阵表示。
通过将操作符应用于向量,我们可以得到一个新的向量,描述了系统的演化过程。
这里的操作符实际上就是矩阵的函数,它将表示系统状态的向量转变为另一个表示。
另一个应用是在图像处理中。
图像可以用一个矩阵表示,其中每个元素代表一个像素的亮度。
通过应用矩阵的函数,我们可以对图像进行各种处理,例如旋转、缩放、平移等。
这些操作可以通过对矩阵的函数进行定义和计算来实现。
此外,矩阵的函数还在统计学和机器学习中有广泛应用。
在统计学中,我们经常需要处理具有大量样本的数据集。
通过将数据集表示为矩阵,我们可以应用各种函数来计算样本的统计特征,例如均值、方差等。
在机器学习中,矩阵的函数用于计算模型的损失函数,以及模型参数的更新规则。
总结起来,矩阵的函数在数学、物理学、图像处理、统计学和机器学习等领域都有广泛应用。
通过定义和计算矩阵的函数,我们可以对复杂的问题进行建模和求解。
矩阵的函数不仅可以进行简单的数学运算,而且可以描述和分析许多实际问题。
大学矩阵数学论文1200字_大学矩阵数学毕业论文范文模板
大学矩阵数学论文1200字_大学矩阵数学毕业论文范文模板大学矩阵数学论文1200字(一):浅谈矩阵在离散数学中的应用摘要:离散数学是计算机学科的一门重要的专业基础课,扎实的基础是非常重要的。
本文就矩阵在离散数学中的各种应用展开讨论,并实例说明。
关键词:矩阵;离散数学;运用引言:随着计算机科学的发展,重点研究有限系统的离散数学已经成为一门越发重要的科学,数字计算机本质上是一个有限结构,它的许多性质都可以在有限数学系统的框架下得到解释。
矩阵是一种有力的数学工具,本文就矩阵在离散数学中的应用展开讨论,总结了矩阵在离散数学中的应用类型,以期对初学者和数学工作者在学习离散数学时提供学习辅导和参考资料。
定义1给出m×n个数,按一定顺序排成一个m行、n列的矩形数表此数表称为m行n列矩阵。
常记a=,或a=(),或。
有关应用及其举例一、二元关系的表示定义2设a,b为有限集,构造一个矩阵,以a的元素和b的元素分别标注其行与列,对于a∈a和b∈b。
视a,b是否具有关系r,在a行和b列交叉处标上1或0.这样得到的矩阵称为关系矩阵。
例如:a={1,2,3,4},在a上定义二元关系r为大于关系,表示x大于y,采用列举法为r={<2,1>,<3,1>,<4,1>,<3,2>,<4,2>,<4,3>}.则关系矩阵为二、图的表示和邻接矩阵定义3设无向图g=,v={v1,v2,vn},e={e1,e2,,em}。
令为节点vi 与边ej关联的次数,则称矩阵为g的关联矩阵,记为m(g)。
例如:无向图g如下所示,则m(g)为:定义4设图g=为有向图,v={v1,v2,vn},即有n个节点,令是vi邻接到vj的边的数目,则称矩阵为g的邻接矩阵,记为a(g)。
例如:有向图g如下三、用矩阵求关系合成和偏序中的盖住关系(一)关系合成设和分别表示关系r和s的矩阵,令m=,则m中的非零元素表示其对应的元素具有关系。
矩阵分析与应用课程设计
矩阵分析与应用课程设计一、背景介绍在大学数学课程中,矩阵分析成为一个非常重要的内容。
矩阵分析作为现代数学的一个重要分支,被广泛应用于物理、经济、组合优化、图形图像处理以及其他数学领域中。
因此,矩阵分析课程的教学往往也是大学数学课程中不可或缺的一部分。
二、课程设计目标本课程设计旨在通过编写矩阵分析代码实践和应用,帮助学生深入了解矩阵分析的原理和应用。
希望通过本次课程设计,学生能够掌握以下技能:1.熟练掌握Python等语言中进行矩阵计算的基本操作;2.掌握矩阵分析的基本理论和应用;3.熟悉Python等语言中常见的矩阵分析工具,如numpy、scipy等,并能够灵活应用。
三、课程设计内容本课程设计涵盖以下内容:•在Python等语言中利用numpy等工具编写矩阵计算程序,包括矩阵求逆、矩阵乘法、矩阵求秩等操作;•矩阵分析的基本理论和应用,包括线性方程组求解、矩阵特征值和特征向量、最小二乘法等;•利用Python等语言中的matplotlib等工具实现二维、三维图形的矩阵可视化,如矩阵的热度图、散点图等;•矩阵分析的实践应用,如图像处理、信号处理、金融风险评估等。
四、课程设计方案本课程设计采用以下方式进行:第一阶段:矩阵计算程序的编写本阶段主要通过引导学生编写Python等语言中的矩阵计算程序,来帮助学生加深对矩阵计算基本操作的掌握。
此阶段具体内容包括:1.矩阵求逆的实现;2.矩阵乘法的实现;3.矩阵求秩的实现。
第二阶段:矩阵分析理论的学习本阶段将重点介绍矩阵分析的基本理论和应用,并通过具体的例子来加深学生对理论的理解。
此阶段具体内容包括:1.线性方程组求解;2.矩阵特征值和特征向量的求解;3.最小二乘法的应用。
第三阶段:矩阵可视化的实现本阶段将介绍Python等语言中的matplotlib等工具,来帮助学生实现二维、三维图形的矩阵可视化。
此阶段具体内容包括:1.矩阵的热度图;2.矩阵的散点图。
第四阶段:矩阵分析的实践应用本阶段将以图像处理、信号处理和金融风险评估为例,介绍矩阵分析的实践应用。
线性代数论文设计(矩阵在自己专业中地应用及举例)
矩阵在自己专业中的应用及举例摘要:I、矩阵是线性代数的基本概念,它在线性代数与数学的许多分支中都有重要的应用,许多实际问题可以用矩阵表达并用相关的理论得到解决。
II、文中介绍了矩阵的概念、基本运算、可逆矩阵、矩阵的秩等容。
III、矩阵在地理信息系统中也有许多的应用,比如文中重点体现的在计算机图形学中应用。
关键词:矩阵可逆矩阵图形学图形变换正文:第一部分引言在线性代数中,我们主要学习了关于行列式、矩阵、方程、向量等相关性比较强的容,而这些容在我们专业的其他一些学科中应用也是比较广泛的,是其它一些学科的很好的辅助学科之一。
因此,能够将我们所学的东西融会贯通是一件非常有意义的事,而且对我们的学习只会有更好的促进作用。
在计算机图形学中矩阵有一些最基本的应有,但是概念已经与线性代数中的有一些不同的意义。
在计算机图形学中,矩阵可以是一个新的额坐标系,也可以是对一些测量点的坐标变换,例如:平移、错切等等。
在后面的文章中,我通过查询一些相关的资料,对其中一些容作了比较详细的介绍,希望对以后的学习能够有一定的指导作用。
在线性代数中,矩阵也占据着一定的重要地位,与行列式、方程、向量、二次型等容有着密切的联系,在解决一些问题的思想上是相同的。
尤其他们在作为处理一些实际问题的工具上的时候。
图形变换是计算机图形学领域的主要容之一,为方便用户在图形交互式处理过程中度图形进行各种观察,需要对图形实施一系列的变换,计算机图形学主要有以下几种变换:几何变换、坐标变换和观察变换等。
这些变换有着不同的作用,却又紧密联系在一起。
第二部分 研究问题及成果1. 矩阵的概念定义:由n m ⨯个数排列成的m 行n 列的矩阵数表⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡ann an an n a a a n a a a ΛM ΛM M KΛ212222111211 称为一个n m ⨯矩阵,其中an 表示位于数表中第i 行第j 列的数,i=1,2,3,…n ,又称为矩阵的元素。
矩阵函数的性质及其应用word资料9页
§7 矩阵函数的性质及其应用一、矩阵函数的性质: 设n n C B A ⨯∈. 1.A e Ae e dtd At At At⋅== proof : 由 ()∑∑⋅==∞=m m m m AtA t m At m e!1!1对任何t 收敛。
因而可以逐项求导。
()∑∞=--=∴01!11m mm At A t m e dt d ()()⎪⎪⎭⎫ ⎝⎛-⋅=∑∞=-11!11m m At m A ()⎪⎪⎭⎫ ⎝⎛⋅=∑k At k A !1AteA ⋅=()()()A e A At m A A t m Atm m m m m ⋅=⋅⎪⎪⎭⎫ ⎝⎛-=⋅-=∑∑∞=∞=---01111!11!11 可见,A 与At e 使可以交换的,由此可得到如下n 个性质 2.设BA AB =,则 ①.At At Be B e =⋅ ②.B A A B B A e e e e e +=⋅=⋅③.()()AA A AA AB A B A B A BA B A B A BA cos sin 22sin sin cos 2cos sin cos cos sin sin sin sin cos cos cos 22=-=⇒+=+-=+=proof :①,由m m BA B A BA AB =⇒=而∑∑∞=∞==⎪⎭⎫⎝⎛=00!1!1m m m m m m AtB A t m B t A m B e()∑∑∞=∞=⋅==00!1!1m mm m m At m B BA t m At e B ⋅=②令()()A B tAt Bt C t e e e +--=⋅⋅由于()0=t C dtd)(t C ∴为常数矩阵 因而E e e e C C t C =-⋅===000)0()1()(当1=t 时,E e e e B A B A =⋅⋅--+ …………………. (@) 特别地 A B -= 有E e e e A A =⋅⋅-0∴ 有 ()A Ae e --=1∴同理有()B Be e --=1代入(@)式 因而有B A B A e e e ⋅=+ 3.利用绝对收敛级数的性质,可得 ①A i A e iA sin cos +=()()iA iAiA iAe e iA e e A ---=+=⇒21sin 21cos ②()()A A A A sin sin cos cos -=-=-4.E A A =+22cos sin()()A E A AE A cos 2cos sin 2sin ππ+=+A E i A e e =+π2二、矩阵函数在微分方程组中的应用—常用于线性监测系统中 1. 一阶线性常系数齐次方程组的通解 AX dtdX = 其中()Tn n n x x x X C A ,,,21Λ=∈⨯ 则有()K e t X At ⋅=,其中()Tn k k k K ,,,21Λ=1eg解方程:⎪⎪⎪⎩⎪⎪⎪⎨⎧+=+-=+-=313212211234xx dtdx x x dtdxx x dt dx解:原方程变为矩阵形式AX dt dX = ⎪⎪⎪⎭⎫⎝⎛--=201034011A ()Tx x x X 321,,=由()()212--=-λλλA E 得⎪⎪⎪⎭⎫ ⎝⎛=→100110002J A1200000-⎪⎪⎪⎭⎫⎝⎛=∴P e e e e P e t tt tAt⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛=∴-3211200000)(k k k P e e e e P t X t tt t2. 一阶线性常系数微分方程组的定解问题:1Th :一阶线性常数微分方程组的定解问题:()()⎪⎩⎪⎨⎧==Tn x x x X AXdt dX)0(,),0(),0(210Λ 有唯一解)0(X e X At ⋅= proof :实际上,由AX dtdX=的通解为K e t X At ⋅=)( 将初值)0(X 代入,得)0(X k =)0(X e X At =∴由1Th 可的定解问题()⎪⎩⎪⎨⎧==Tn t x t x t x t X AX dt dX)(,),(),()(002010Λ的唯一解为()()00)(t X e t X t t A ⋅=-2eg 求定解问题:()()⎪⎩⎪⎨⎧==Tx Axdt dx1,00,⎪⎪⎭⎫ ⎝⎛--=1221A 的解解:由0=-A E λ 得i x 32,1±=对应的特征向量记为:Ti ⎪⎪⎭⎫ ⎝⎛+=231,1α ⎪⎪⎭⎫⎝⎛-=231,1i β 则,于是矩阵:⎪⎪⎭⎫ ⎝⎛-+=23123111i iP13300--⋅⎪⎪⎭⎫ ⎝⎛⋅=∴P e e P eit itAt⎪⎪⎪⎪⎭⎫⎝⎛+=⎪⎪⎭⎫ ⎝⎛=t t t e t X At 3sin 313cos 3sin 3210)(练习:求微分方程组1132123313383625dx x x dt dx x x x dt dx x x dt ⎧=+⎪⎪⎪=-+⎨⎪⎪=--⎪⎩满足初始条件123(0)1,(0)1,(0)1x x x ===的解。
矩阵的函数范文范文
矩阵的函数范文范文一、矩阵的函数的定义和性质在数学中,一个矩阵的函数将一个矩阵作为输入,经过一系列的运算和变换后,产生一个新的矩阵作为输出。
矩阵的函数可以是线性函数,也可以是非线性函数。
线性函数的定义是:如果f(A + B) = f(A) + f(B)和f(kA) = kf(A)对于任意的矩阵A、B和标量k都成立,则称函数f为线性函数。
非线性函数则没有满足这两个性质。
矩阵的函数可以是一元函数,也可以是多元函数。
一元函数是指只有一个矩阵作为输入,多元函数则有多个矩阵作为输入。
例如,对于方阵A,我们可以定义一个矩阵的一元函数f(A)=A^2,将输入矩阵A乘以它自己得到输出矩阵。
矩阵的函数有一些重要的性质。
首先,函数的定义域是一个矩阵的集合,通常是所有满足一定条件的矩阵组成的集合。
其次,函数的值域也是一个矩阵的集合,通常是满足一定条件的矩阵组成的集合。
最后,函数的图像是一个由输入矩阵和输出矩阵组成的集合。
二、矩阵的函数的运算矩阵的函数可以进行加法、减法、乘法和除法等运算。
这些运算的定义和性质与传统的算术运算类似。
例如,对于两个矩阵函数f(A)和g(A),它们的加法定义为(f+g)(A)=f(A)+g(A),减法定义为(f-g)(A)=f(A)-g(A)。
矩阵函数的乘法定义为(f∘g)(A)=f(g(A)),即先对输入矩阵进行函数g的运算,然后再对结果进行函数f的运算。
矩阵的函数还可以进行逆运算和幂运算。
逆运算的定义是:对于一个矩阵函数f(A),如果存在另一个矩阵函数g(A),使得(f∘g)(A)=(g∘f)(A)=A对于任意的矩阵A都成立,则称函数g为函数f的逆函数,记作g=f^(-1)。
幂运算的定义是:对于一个矩阵函数f(A)和一个非负整数n,我们定义f^n(A)=(f∘f∘...∘f)(A)(共n个f),即对输入矩阵进行n次函数运算。
三、常见的矩阵函数在实际的应用中,矩阵的函数具有广泛的应用。
以下是一些常见的矩阵函数的示例。
数学毕业论文-矩阵函数的性质及其应用
数学毕业论文-矩阵函数的性质及其应用数学毕业论文-矩阵函数的性质及其应用目录中文标题 (1)中文摘要、关键词 (1)英文标题 (1)英文摘要、关键词 (1)正文1引言 (2)2矩阵函数的有关概念及性质 (2)2.1矩阵函数的有关概念 (2)2.2常用矩阵函数的性质 (3)3矩阵函数值的计算 (6)3.1利用Hamilton-Cayley定理 (6)3.2利用相似对角化 (6)3.3利用Jordan标准形 (7)3.4待定系数法 (10)4矩阵函数的应用 (11)4.1矩阵函数在求矩阵的`幂及中的应用 (11)4.2矩阵函数的应用 (13)4.3矩阵函数在求解电路暂态响应中的应用 (15)5结束语 (18)参考文献 (19)致谢 (20)矩阵函数的性质及其应用摘要本文利用收敛的幂级数的和定义了矩阵函数,分析讨论了矩阵函数的有关性质,并介绍了矩阵函数值的几种求解方法:利用Hamilton-Cayley定理、利用相似对角化、利用Jordan标准形、待定系数法,最后探讨了矩阵函数的1些具体应用.关键字:矩阵函数;矩阵函数值;矩阵的幂The properties of Matrix Function and Its ApplicationsABSTRACTIn this paper, matrix function is definded by using convergent power series,and the basic properties of the matrix function is introduced. Several kinds of methods to solve matrix function value,,such as Hamilton Cayley theorem, Similitude diagonal , Jordan canonical form and undetermined method ,have been presented.And finally,some applicaions of matrix function have also been given.Key word:matrix function;matrix function value;powers of matrix【包括:毕业论文、开题报告、任务书】【说明:论文中有些数学符号是编辑器编辑而成,网页上无法显示或者显示格式错误,给您带来不便请谅解。
矩阵的函数范文
矩阵的函数范文矩阵函数是指将一个矩阵作为输入,返回一个新的矩阵作为输出的数学函数。
矩阵函数在许多领域中都有重要的应用,如线性代数、微积分、图论等等。
本文将探讨矩阵函数的定义、性质以及一些常见的矩阵函数的应用。
一、矩阵函数的定义和性质:1.定义:矩阵函数可以定义为一个从矩阵空间到矩阵空间的映射,即对于一个给定的矩阵A,矩阵函数f(A)返回一个新的矩阵B。
一般来说,矩阵函数可以是任意的,它可以是线性的或非线性的,可以是单值的或多值的。
2.线性矩阵函数:线性矩阵函数是指满足以下两个性质的矩阵函数:(1)f(A+B)=f(A)+f(B):对于任意的矩阵A和B,有f(A+B)=f(A)+f(B);(2) f(cA) = cf(A):对于任意的矩阵A和标量c,有f(cA) = cf(A)。
3.非线性矩阵函数:非线性矩阵函数是指不满足线性性质的矩阵函数。
非线性矩阵函数的性质较为复杂,常常需要利用数值方法进行计算。
4.特殊矩阵函数:特殊矩阵函数是指具有一些特定性质的矩阵函数,如对称函数、正定函数等。
特殊矩阵函数在各个领域中都有广泛的应用。
5. 矩阵函数的迹和行列式:对于一个矩阵函数f(A),其迹和行列式可以定义为其矩阵的迹和行列式的函数,即tr(f(A))和det(f(A))。
二、常见的矩阵函数:1.幂函数:幂函数f(A)=A^k将一个矩阵A自乘k次。
2. 指数函数:指数函数f(A) = e^A将一个矩阵A进行Taylor展开,得到一个无限级数。
3. 对数函数:对数函数f(A) = ln(A)将一个矩阵A进行类似于指数函数的Taylor展开,得到一个无限级数。
4. 三角函数:三角函数sin(A)、cos(A)和tan(A)分别将矩阵A中的每个元素作为角度计算其三角函数值。
5. 反三角函数:反三角函数asin(A)、acos(A)和atan(A)分别将矩阵A中的每个元素作为三角函数值计算其对应的角度。
6. 矩阵修正函数:矩阵修正函数f(A) = max(0, A)将矩阵A中的每个元素与0进行比较,将小于0的元素修正为0。
矩阵在数学中的应用本科毕业论文
本科毕业论文(设计)题目矩阵在数学中的应用____________________________________毕业论文(设计)原创性声明本人所呈交的毕业论文(设计)是我在导师的指导下进行的研究工作及取得的研究成果。
据我所知,除文中已经注明引用的内容外,本论文(设计)不包含其他个人已经发表或撰写过的研究成果。
对本论文(设计)的研究做出重要贡献的个人和集体,均已在文中作了明确说明并表示谢意。
作者签名:日期:毕业论文(设计)授权使用说明本论文(设计)作者完全了解**学院有关保留、使用毕业论文(设计)的规定,学校有权保留论文(设计)并向相关部门送交论文(设计)的电子版和纸质版。
有权将论文(设计)用于非赢利目的的少量复制并允许论文(设计)进入学校图书馆被查阅。
学校可以公布论文(设计)的全部或部分内容。
保密的论文(设计)在解密后适用本规定。
作者签名:指导教师签名:日期:日期:注意事项1.设计(论文)的内容包括:1)封面(按教务处制定的标准封面格式制作)2)原创性声明3)中文摘要(300字左右)、关键词4)外文摘要、关键词5)目次页(附件不统一编入)6)论文主体部分:引言(或绪论)、正文、结论7)参考文献8)致谢9)附录(对论文支持必要时)2.论文字数要求:理工类设计(论文)正文字数不少于1万字(不包括图纸、程序清单等),文科类论文正文字数不少于1.2万字。
3.附件包括:任务书、开题报告、外文译文、译文原文(复印件)。
4.文字、图表要求:1)文字通顺,语言流畅,书写字迹工整,打印字体及大小符合要求,无错别字,不准请他人代写2)工程设计类题目的图纸,要求部分用尺规绘制,部分用计算机绘制,所有图纸应符合国家技术标准规范。
图表整洁,布局合理,文字注释必须使用工程字书写,不准用徒手画3)毕业论文须用A4单面打印,论文50页以上的双面打印4)图表应绘制于无格子的页面上5)软件工程类课题应有程序清单,并提供电子文档5.装订顺序1)设计(论文)2)附件:按照任务书、开题报告、外文译文、译文原文(复印件)次序装订3)其它目录摘要 (I)Abstract. (II)1 前言 (1)2 有关概念及重要结论 (1)2.1矩阵的概念 (1)2.2矩阵的秩 (2)2.3矩阵的逆 (3)2.4 用矩阵表示二次型 (3)3 矩阵的应用 (6)3.1矩阵的高次幂 (6)3.1.1 矩阵的幂 (6)3.1.2矩阵高次幂的求法 (7)3.2 解线性方程组 (13)3.2.1线性方程组的有解判定定理 (13)3.2.2 线性方程组一般形式的运用 (14)3.3 解矩阵方程 (16)3.4 矩阵对角化方法 (19)3.4.1 讨论对于有n个特征单根的n阶方阵 (19)3.4.2 讨论对于有特征重根的n阶方阵 (21)结论 (24)致谢 (24)参考文献 (24)矩阵及应用杨灿(重庆三峡学院数学与统计学院数学与应用数学专业2010级重庆万州 404100)摘要:矩阵理论既是学习经典数学的基础,又是一门很有实用价值的数学理论.随着科学技术的发展,这一理论已成为现代各科技领域处理大量数据的有效工具.本文就是利用矩阵的基本理论,把矩阵作为计算工具,对实际问题如方程组的解、矩阵的幂、二次型进行了较为系统的研究并简化了一些计算.关键词: 矩阵;矩阵的幂;线性方程组Matrix and Its ApplicationYANG Can(Grade 2010, Mathematics and Applied Mathematics, College of Mathematics and statistics, Chongqing Three Gorges University, Wan Zhou, Chongqing 404100 )Abstract:Matrix theory is not only the foundation of learning classical mathematics,but also is a very useful mathematical theory.With the development of science and technology,this theory has become the effective tool for modern technology in the field of large amounts of data.This article is on the undamental theory of matrix,the matrix as a calculation tool,the practical problems such as the solution of the equations,the power of matrix,the two type are systematically studied and some simplified calculation.Keywords:Matrix; The power of matrix; Linear equation2014届数学与应用数学专业毕业设计(论文)1 前言矩阵是数学中的一个重要的基本概念,是代数学的主要研究对象之一,也是数学研究和应用的一个重要工具.“矩阵”这个词是由西尔维斯特首先使用的,他是为了将数字的矩形阵列区别于行列式而发明了这个术语.而实际上,矩阵在它的课题诞生之前就已经发展的很好了.18世纪中期,数学家们开始研究二次曲线和二次曲面的方程简化问题,即二次型的化简.在这一问题的研究中,数学家们得到了与后来的矩阵理论密切相关的许多概念和结论.1748年,瑞士数学家欧拉(L .Euler,1707—1783)在将三个变数的二次型化为标准形时,隐含地给出了特征方程的概念.1773年,法国数学家拉格朗日(J .L .Lagrange,1736—1813)在讨论齐次多项式时引入了线性变换.1801年德国数学家高斯(C .F .Gauss,1777一1855)在《算术研究》中,将欧拉与拉格朗日的二次型理论进行了系统的推广,给出了两个线性变换的复合,而这个复合的新变换其系数矩阵是原来两个变换的系数矩阵的乘积.另外,高斯还从拉格朗日的工作中抽象出了型的等价概念,在研究两个互逆变换的过程中孕育了两个矩阵的互逆概念.在线性方程组的讨论中,我们看到,线性方程组的一些重要性质反映在它的系数矩阵和增广矩阵的性质上,并且解线性方程组的过程也表现为变换这些矩阵的过程.除了线性方程组之外,还有大量的各种各样的问题也都提出矩阵的概念,并且这些问题的研究常常反映为有关矩阵的某些方面的研究,甚至于有些性质完全不同的、表面上完全没有联系的问题,归结成矩阵问题以后却是相同的.这使矩阵成为数学中一个极其重要的应用广泛的概念,因而也就使矩阵成为代数特别是线性代数的一个主要研究对象,也是处理高等数学很多问题的有力工具.矩阵的秩是一个基本的概念,也是矩阵最重要的数量特征之一,它在初等变换下是一个不变量.矩阵的秩是反映矩阵固有特性的一个重要概念,无论是在线性代数中,还是在解析几何中,甚至在概率论中,都有不可忽略的作用.矩阵方幂在高等代数题解、矩阵稳定性讨论及预测、控制等方面有广泛的应用,它的求解原理贯穿于代数教学过程的始终,可以用到矩阵各方面的知识.其计算量往往较大,但方法适当,可大大简化其计算难度.本文将给出六种求矩阵方幂地方法.矩阵方程是矩阵运算的一部分,这里我们主要讨论如何求解矩阵方程的问题.掌握简单的矩阵方程的求法,对于求解复杂的矩阵方程有很大帮助.2 有关概念及重要结论2.1矩阵的概念为了便于叙述并考虑以后的应用,我们引进矩阵的概念.由mn 个数排列而成的m 行(横的)n 列(纵的)的表⎪⎪⎪⎪⎪⎭⎫⎝⎛mn m m n n a a a a a a a a a 212222111211称为一个n m ⨯杨灿:矩阵及其应用矩阵.定义 1 把矩阵A 的行换成同序数的列得到的新矩阵, 称为A 的转置矩阵, 记作T A (或A ').即若,212222111211⎪⎪⎪⎪⎪⎭⎫⎝⎛=mn m m n n a a a a a a a a a A 则⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=mn n n m m Ta a a a a a a a a A 212221212111. 2.2矩阵的秩定义2 所谓矩阵的行秩就是指矩阵的行向量组的秩;所谓矩阵的列秩就是指矩阵的列向量组的秩.引理1 如果齐次方程组⎪⎪⎩⎪⎪⎨⎧=+++=+++=+++000221122221211212111n sn s s nn n n x a x a x a x a x a x a x a x a x a 的系数矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛=sn s s n n a a a a a a a a a A 212222111211的行秩n r <,那么它有非零解.定理1 矩阵的行秩与列秩相等.定理 2 n n ⨯矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛=nn n n n n a a a a a a a a a A 212222111211的行列式为零的充分必要条件是A 的秩小于n .推论 1 齐次线性方程组⎪⎪⎩⎪⎪⎨⎧=+++=+++=+++000221122221211212111n nn n n nn n n x a x a x a x a x a x a x a x a x a 有非零解的充分必要条件是它的系数矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛=nn n n n n a a a a a a a a a A 212222111211的行列式等于零.2.3矩阵的逆我们知道,n 阶单位矩阵E 单位性质,即对于任意n 阶方阵A 都有A EA AE ==,是否存在n 阶方阵B 使得E AB =呢?即是否与数域P 中数一样的性质:1)0(1=⋅⇒∈≠∀-a a P a .为此,我们引进逆矩阵的概念.定义1 n 阶方阵A 称为可逆的,如果有n 阶方阵B ,使得E BA AB ==. (2.3.1)这里E 是n 级单位矩阵.并且称B 为A 的一个逆矩阵.定义2 如果矩阵B 适合(2.3.1),那么B 就称为A 的逆矩阵,记为1-A . 定理1 n 阶矩阵A 可逆的充分必要条件是A 非退化,此时,A 的逆矩阵为0,1*1≠==-A d A dA . 定理2 给出了矩阵可逆时逆矩阵的计算公式.下面给出可逆矩阵的一些性质: 性质1 如果n 阶方阵A 可逆,那么0≠=A d ,并且dA 11=-. 性质2 如果矩阵B A ,同级且都可逆,那么T A 与AB 也可逆,且11111)(,)()(-----==A B AB A A T T .性质3 如果n 阶方阵A 可逆,那么kA N k ,∈∀也可逆,并且k k A A )()(11--=. 性质4 如果n 阶方阵A 可逆,那么k A Z k ,∈∀也可逆,并且k k A A )()(11--=.性质5 如果n 阶方阵A 可逆,那么Z l k ∈∀,,有l k l k k l kl l k A A A A A A +===,)()(. 定理3 A 是一个n s ⨯矩阵,如果P 是s s ⨯可逆矩阵,Q 是n n ⨯可逆矩阵,那么)()()(A r AQ r PA r ==.推论1 在定3的假设下有,)()(A r PAQ r =成立.2.4 二次型及矩阵表示定义1 设P 是一个数域,一个系数ij a 在数域P 中的n x x x ,,,21 的二次齐次多项式 jinj i ij i ni ii n xx a x a x x x f ∑∑≤≤≤=+=121212),,,( . (2.4.1)定义2 记ij ji a a =,把n 元二次型(2.4.1),写成对称形式j i ni nj ij n x x a x x x f ∑∑===1121),,,( . (2.4.2)这样,系数ij a 可以构成一个n n ⨯对称矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛==nn n n n n nn ij a a a a a aa a a a A 212222111211)(, (2.4.3) 称(2.4.3)为n 元二次型(1)的矩阵. 令Tn x x x x ),,,(21 =,则有i n i j nj ij j i n i n j ij n x x a x x a x x x f ∑∑∑∑======111121)(),,,( ,=⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛∑∑∑===n j j nj n j j j n j j j n x a xa x a x x x 1121121),,,( ,=⎪⎪⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎪⎪⎭⎫⎝⎛n nn n n n n n x x x a a a a a a a a a x x x 2121222211121121),,,(,=Ax x T, (2.4.4)这就是二次型的矩阵表示.对确定的n 元二次型(2.4.1),就确定唯一的对称矩阵(2.4.3)通过(2.4.4)联系起来,即Ax x xx a x x x f T jin i nj ij n ==∑∑==1121),,,( .因此,一个n 元二次型(2.4.1)对应一个n 阶对称矩阵.每个二次型都有一个对称矩阵与之对应;反之,每个对称矩阵也有一个二次型与之对应.二次型与它的矩阵是相互唯一确定的.一般地,关于二次型的矩阵有下列结果.定理1 设B 是n n ⨯矩阵,则Bx x x x x f Tn =),,,(21 是一个二次型,它的矩阵为2BB T +.2.5 特征值与特征向量n 维线性变换空间V 与矩阵空间nn p ⨯是同构关系,可以通过矩阵来研究线性变换的性质,我们希望找到一组基,,,21n ξξξ 使得线性变换A L 在这组基下的矩阵A 的形式最简单.这个问题的一个简单设想是A 是否可以是对角形式?即),,,(,,,3,2,1,21n j j j A a a a diag A n j a L ===ξξ.这个设想可以归结为:对线性空间V 的线性变换ξξk L A =,P k ∈.这就是线性变换的特征值与特征向量.定义1 设A L 是数域P 上线性空间V 的一个线性变换,如果对于数域P 中一数0λ,存在一个非零向量ξ,使得ξλξ0=A L .那么0λ称为A L 的是一个特征值,而ξ称为A L 的属于特征值0λ的一个特征向量.定义2 设A 是数域P 上一n 级矩阵, λ是一个文字. 矩阵A -E λ的行列式nnn n nn a a a a a a a a a ---------=A -E λλλλ212222111211,称为A 的特征多项式, 这是数域P 上的一个n次多项式.上面的分析说明, 如果0λ是线性变换A L 的特征值, 那么0λ一定是矩阵A 的特征多项式的一个根; 反过来, 如果0λ是矩阵A 的特征多项式在数域P 中的一个根, 即00=-E A λ, 那么齐次线性方程组⎪⎪⎩⎪⎪⎨⎧=-----=---+-=----0)(0)(0)(022111222012112121110n nn n n nn n n x a x a x a x a x a x a x a x a x a λλλ (2.5.1)就有非零解. 这时,如果),,,(00201n x x x 是方程组(2.5.1)的一个非零解, 那么非零解向量.n n x x x ζζζζ0202101+++= .满足(2.5.1)式, 即0λ是线性变换A L 的一个特征值, ζ就是属于特征值0λ的一个特征向量.定理1 设A L 是数域P 上n 维线性空间V 的一个变换,则P ∈0λ是A L 的一个特征值当且仅当0λ是A L 的特征多项式)()(λλA L f f A≡的一个根.定理2 设0λ是线性空间V 的线性变换A L 的一个特征值,则集合{}V L V A ∈==ααλααλ,00 (2.5.2)构成V 的一个子空间.在有限维情形,)(dim 00A E R n V --=λλ,其中,V n dim =,A 是A L 在V 在某个基下的矩阵.定义3 设0λ是线性空间V 的线性变换A L 的一特征值,式(2.5.2)定义的V 的子空间称为A L 的对应特征值0λ的特征子空间0λV因此, 确定一个线性变换A 的特征值与特征向量的方法可以分成一下几步: (1)在线性空间V 中取一组基n ζζζ,,,21 , 写出A L 在这组基下的矩阵A ;(2)求出A 的特征多项式A -E λ在数域P 中全部的根, 它们也就是线性变换A L 的全部特征值;(3)把所得的特征值逐个代入方程组(2.5.1)式, 对于每一个特征值, 解方程组(2.5.1)式,求出一组基础解系, 它们就是属于这个特征值的几个线性无关的特征向量在基n ζζζ,,,21 下的坐标, 这样, 我们也就求出了属于每个特征值的全部线性无关的特征向量.矩阵A 的特征多项式的根有时也称为A 的特征值, 而相应的线性方程组(2.5.1)式的解也就称为A 的属于这个特征值的特征向量.3 矩阵的应用3.1矩阵的高次幂3.1.1 矩阵的幂定义1 设方阵n n ij a A ⨯=)(, 规定.,,0为自然数个k A A A A E A k k⋅⋅⋅==k A 称为A 的k 次幂.方阵的幂满足以下运算规律(假设运算都是可行的): (1) );,(为非负整数n m A A A n m n m +=(2) .)(mn n m A A =注意: 一般地,,)(m m m B A AB ≠ m 为自然数命题1 设B A ,均为n 阶矩阵,,BA AB = 则有,)(m m m B A AB = m 为自然数,反之不成立.3.1.2 矩阵高次幂的求法矩阵方幂在高等代数题解、矩阵稳定性讨论及预测、控制等方面有广泛的应用,它的求解原理贯穿于代数教学过程的始终,可以用到矩阵各方面的知识.其计算量往往较大,但方法适当,可大大简化其计算难度.本文将给出六种求矩阵方幂地方法.3.1.2.1 利用凯莱——哈密尔顿(Cayley —Hamilton )定理求方阵的幂定理1 (Cayley —Hamilton 定理)设A 是n 阶矩阵,)(λf 是A 的特征多项式,则0)(=λf . 设A 是数域P 上n 阶方阵,其特征多项式为)(λf ,为求A n(n 是正整数),令n g λλ=)(,做带余除法,)()()()(λλλλr q f g +=.由定理1知,)()(λλr g =,并且)(λr 的次数小于)(λg 的次数,进而可得n r g A =A =A )()(.利用上定理求幂时在计算过程中可分为两种情形:1、所求矩阵的幂指数相对较低,可直接利用定理1及余式定理求出)(λr .例1 已知 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=101121002A ,求5A .解 令5)(λλ=g 矩 阵A 的 特 征 多 项 式 为)1()2(11121002)det()(2--=-----=A -I =λλλλλλλf 做带余除法,6811649)1750)(()(225+-+++==λλλλλλλf g 于是,由定理1知I +A -A =I +A -A ++A +A A =A =A 68116496811649)1750)(()(2225f g⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=1000100016810112100211610334300449 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=10313132310032 2、所求矩阵的幂指数相对较高,不便用上法直接求出余式.此种情形下矩阵的特征多项式有重根和无重根时分别给出下面的解法.(1)矩阵的特征多项式无重根.对于i ni i c q f r q f g λλλλλλλ∑=+=+=1)()()()()()(,以其n 个不同的特征值分别代入此式即可求出)(λr .例2 已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=A 211110101,求991003A -A .解 令991003)(λλλ-=g .矩阵A 的特征多项式为)3)(1(211110101)det()(--=-------=A -I =λλλλλλλλf .做带余除法,注意到)(λf 的次数是3,即c b a q f g +++=-=λλλλλλλ299100)()(3)(. 以3,1,0=λ分别代入上式得0)0(==c g .2)1(-=++=c b a g .039)3(=++=c b a g . 所以0,3,1=-==c b a .由定理1 ,A -A =I +A +A =A -A =A 33)(2299100c b a g⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0000110112111101013631321312.(3)矩阵的特征多项式有重根.同上法,为获得足够的信息求出)(λr ,可对)()()()(λλλλr q f g +=求导.例3 已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=A 210111111,求100A .解 A 的特征多项式是)2()1()det()(2--=A -E =λλλλf 令100)(λλ=g ,做带余除法0122)()()(b b b q f g +++=λλλλλ以2,1=λ分别代入上式,有⎩⎨⎧=++==++=100012012234)2(1)1(b b b g b b b g 为求)2,1,0(=i b i ,就)(λg 对λ求导得10012'2'1002)()()()]1()2)(1(2[)(λλλλλλλλλ=+++-+--=b b q g q g 以1=λ代入上式,有100212=+b b ,从而求得 1000201110022102,3022,2201-=-=-=b b b , 于是 I +A +A =A0122100b b b .3.1.2.2 对于秩为1的n 阶方阵A 有下面定理定理1 对于n 阶方阵A,若1)(=A rank ,那么A 可分解为一个列向量与一个行向量的乘积'αβ=A ,其中⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=n n b b b b a a a a .,.321321 βα.例4 已知⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=A 1233321231211,求n A . 解 显然1)(=A rank ,并且⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=A 1233321231211⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=3121132`1,而331211321=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡,所以⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=A =⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=A ---123332123121133312113213111n n n n .3.1.2.3 可分解为数量矩阵和零幂矩阵之和的情况要点 观察推敲矩阵A ,看其是否可以分解为一个数量矩阵E λ与一个零幂矩阵P 之和,即P +E =λA ,其中O m ≠P ,但O m =P+1,因为数量矩阵E λ和P 可以交换,于是由二项式定理得m m n kn n k k n nk k k n nk nnm n n k n k n A P ⎪⎪⎭⎫ ⎝⎛++P +=P ⎪⎪⎭⎫ ⎝⎛=P E ⎪⎪⎭⎫ ⎝⎛=P +E =---=-=∑∑λλλλλλ 100)()(.例5 已知矩阵,⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=2000420000210042A ,求n A . 解 观察矩阵A 的特点,可先将其分块写成⎪⎪⎭⎫ ⎝⎛=C O O B A ,其中⎪⎪⎭⎫ ⎝⎛=2142B ,⎪⎪⎭⎫⎝⎛=2042C ,则⎪⎪⎭⎫ ⎝⎛=n nn C OO B A ,下面就先求n B 和nC . 显然1)(=B r ,即pq B =,这里⎪⎪⎭⎫⎝⎛=12p ,⎪⎪⎭⎫⎝⎛=21q ,且4=qp ,所以B B n n 14-=. 至于P +E =⎪⎪⎭⎫⎝⎛+E =⎪⎪⎭⎫ ⎝⎛=2004022042C ,⎪⎪⎭⎫ ⎝⎛=P 0040满足O P =2,代入上述给出的二 次项式公式⎥⎦⎤⎢⎣⎡⋅=P +E =P E +E =+E =---nn nn n n nnnn n n P C 2024222)2()2()2(111. 因此本题得解 ⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡⋅⋅⋅=---n n n n n n n A 2024200004200442111. 3.1.2.4 归纳法例6 已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100101αβαA ,求其n 次幂. 解 先来计算A 的较低次幂2A 和3A ,由矩阵乘法直接计算得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+=10021022122αβααA ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+=100310333123αβααA ,……由此猜想⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+-=100102)1(12αβααn n n n n A n. 以下用数学归纳法加以证明. (1)当1=n 时成立.(2)归纳假设结论对k n =时亦成立,即⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+-=100102)1(12αβααk k k k k A k . 所以当1+=k n 时,A A Ak k =+1,而⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+(++++(=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+-=100)110)1(2)1()11100101100102)1(122αβαααβααβααk k k k k k k k k k A A k , 即当1+=k n 时成立,从而证明结论成立.即⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+-=100102)1(12αβααk k k k k A k. 3.1.2.5 利用相似变换法要点 若已知矩阵可以经过相似变换化为对角阵时,即存在可逆矩阵P ,使Λ=AP P -1,其中Λ为对角阵,其对角线上元素为矩阵A 的特征值.由上可得1-PΛP =A ,1-P PΛ=A n n .于是求A的方幂就转化为求过渡矩阵P 和对角阵nΛ,而对于P 和阵nΛ,我们应用代数知识要好求得多了,具体如下:例7 已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=A 122212221,求其n 次幂. 解 经过计算,矩阵A 的特征值1-=λ和5=λ,对于特征值1-=λ有线性无关特征向量T )101(1-=α和()3011Tα=-()T 1102-=α.对于特征值5=λ有特征向量()T 1113=α.令()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--==P 111110101,,321ααα,即P 可逆,且有,5000)1(000)1(,5000100011⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=Λ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=Λ=AP P -n n n n 于是.,11--P PΛ=A PΛP =A nn计算得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+-+-+-+-+-+-+-+-+-=A ++++++n n nn n n n n n n nn n n n n nn n 52)1(5)1(5)1(5)1(52)1(5)1(5)1(5)1(52)1(31111111.3.1.2.6 利用Jordan 标准形例8 已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----=A 411301621,求k A .解 第一步:首先求矩阵A 的若尔当标准形.由⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--+=A -E 2)1(0001000141131621λλλλλλ.从而初等因子为)1(-λ,2)1(-λ,故A 的若尔当标准形⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100010001J .第二步:求可逆矩阵T 使J AT T =-1,即TJ AT =.设),,(321ααα=T ,所以有332211,,αααααα=A =A =A .由22αα=A 得32)(αα-=A -E ,设()Tx x x 3212,,=α,()Ty y y 3213,,=α,则由⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-----→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=A -E 3221321000311622311311622)(y y y y y y y , 而32)(αα-=A -E 有解,故32y y =,又33αα=A ,从而0)(3=A -E α即0311311622321=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---y y y , 于是有03321=-+y y y ,所以得212y y =.令132==y y ,则21=y .于是T )112(3=α,再解T )001(2-=α.于是求得()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-==101100213,,321αααT . 第三步:由第二步得1-=A TJT .⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+------=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-==A -k k k k k k k kk k TTJ k k 31316221010311110100100011011002131.3.2 解线性方程组3.2.1线性方程组的有解判定定理定理1 (克拉默法则) 如果线性方程组⎪⎪⎩⎪⎪⎨⎧=+++=+++=+++nn nn n n n n n n b x a x a x a b x a x a x a b x a x a x a 22112222212111212111 (4.2.1)的系数矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛=nn n n n n a a a a a a a a a A 212222111211的行列式,0≠=A d 那么线性方程组(4.2.1)有解,并且解是唯一的,解可以通过系数表为,,,,2211dd x d dx d d x n n ===其中j d 是矩阵A 中第j 列换成方程组的常数项n b b b ,,,21 所成的矩阵的行列式,即.,,2,1,1,1,121,221,22111,111,111n j a a b a a a a b a a a a b a a d nnj n nj n n nj j nj j j==+-+-+- 定理(线性方程组的有解判定定理) 线性方程组⎪⎪⎩⎪⎪⎨⎧=+++=+++=+++sn sn s s n n n n b x a x a x a b x a x a x a b x a x a x a 22112222212*********有解的充分必要条件为它的系数矩阵⎪⎪⎪⎪⎪⎭⎫⎝⎛=sn s s n n a a a a a a a a a A 212222111211与增广矩阵⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=s sns s n n b a a a b a a a b a a a A 21222221111211有相同的秩.3.2.2 线性方程组一般形式的运用例9 求下述齐次线性方程组的一个基础解系⎪⎪⎩⎪⎪⎨⎧=+++-=++-+-=---+-=+-+-0931050320117630426354321543215432154321x x x x x x x x x x x x x x x x x x x x把方程组的系数矩阵经过初等行变换化成阶梯形矩阵:⎪⎪⎪⎪⎪⎭⎫⎝⎛---→⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---------000000000078100650219131051312111716341263于是方程组的一般解为:⎩⎨⎧+=--=543542178652x x x x x x x 其中542,,x x x 是自由未知量.令0,0,1542===x x x 得)0,0,0,1,2(1=η 0,1,0542===x x x 得)0,1,8,0,5(2-=η 1,0,0542===x x x 得)1,0,7,0,6(3-=η 这里321,,ηηη就是方程组的一个基础解系.例10 解线性方程组:⎪⎪⎩⎪⎪⎨⎧-=++-+-=++-+-=---+-=++-+2573431272327225354321543215432154321x x x x x x x x x x x x x x x x x x x x解 把此方程组的增广矩阵经过初等行变换化成阶梯形矩阵:⎪⎪⎪⎪⎪⎭⎫⎝⎛------→⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--------→⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----------000000666100121010875001000000666100545110112111257343112111721132712253从而得到此方程组的一般解为:⎪⎩⎪⎨⎧-+=---=-+=66662875543542541x x x x x x x x x 其中54,x x 是自由未知量. 对于方程个数与未知量个数相等的非齐次线性方程组,如果它的系数行列式不为零,我们还可以用克莱姆法则求解.但是这种方法计算量很大,因此我们一般不用它,只是对少数字母系数的方程组采用克莱姆法来进行求解.例11 非齐次线性方程组⎪⎪⎩⎪⎪⎨⎧=+--=++-=+--=+--321934443522134321432143214321x x x x a x x x x x x x x x x x x 求当a 为何值时方程组有解?此时有多少解?解 把方程组的增广矩阵经过初等行变换化成阶梯形矩阵:⎪⎪⎪⎪⎪⎭⎫⎝⎛---→⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-------00000340000211001131132211193444352211311a a ⎪⎪⎪⎪⎪⎭⎫⎝⎛---→⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-------00000340000211001131132211193444352211311a a 显然,当34≠a 时,方程组无解;当34=a 时,方程组无解,此时由于阶梯形矩阵的非零行有2行,而未知量有4个,所以方程组有无穷多个解,易求出一般解为⎩⎨⎧+-=+-=27443421x x x x x 其中42,x x 是自由未知量.3.3 解矩阵方程矩阵方程是矩阵运算的一部分,这里我们主要讨论如何求解矩阵方程的问题.掌握简单的矩阵方程的求法,对于求解复杂的矩阵方程有很大帮助.简单的矩阵方程有三种形式:.,,C AXB C XA C AX ===如果这里的A 、B 都是可逆矩阵,则求解时需要找出矩阵的逆,注意左乘和右乘的区别.它们的解分别为.,,1111----===B A X CA X C A X例如,求解方程C AC =先考察A 是否可逆,如果A 可逆时,方程两边同时左乘1-A ,得,11C A AX A --=即,1C A X -=这里要注意只能左乘不能右乘,因为矩阵的乘法不满足交换律.同样,对于方程,C XA =只能右乘1-A ,得,11--=CA XAA 即.1-=CA X 而对于方程,C AXB =只能是左乘1-A 而右乘1-B ,得,1111----=CB A ACBBA 即.11--=CB A X看下面解矩阵方程例题:例12 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡315432343122321X 解 先求出1-A ,则,111253232313431223211⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-则⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=332123315432111253232313154321343122321X 例13 ⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡212101343122321X解 先求出1-A ,则,111253232313431223211⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-则⎥⎥⎦⎤⎢⎢⎣⎡--=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡-=-27525120111253232312121013431223212121011X 例14 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡3154321325343122321X解 先求出1-A ,则,111253232313431223211⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡-532113251, 则⎥⎦⎤⎢⎣⎡--⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=--532131543211125323231132531543234312232111X ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎦⎤⎢⎣⎡--⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=131148735331332123当矩阵方程C AXB C XA C AX ===,,中的A 、B 不是方阵或者是不可逆的方阵时,前面的方法就不能用了.这时,我们需要用待定元素法来求矩阵方程.设未知矩阵X 的元素为ij x ,即)(ij x X =,然后由所给的矩阵方程列出ij x 所满足的线性方程组,通过解线性方程组求出所有元素ij x ,从而得到所求矩阵)(ij x X =.例15 解矩阵方程⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡-4152102011X解 利用元素法,先确定X 的行数等于左边矩阵的行数3,X 的列数等于积矩阵的列数2,则X 是23⨯的矩阵.设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=2221y y y x x x X ,则⎥⎦⎤⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎦⎤⎢⎣⎡-41521020112121y y y x x x. 即⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡++--4152222111y y x x y y x x ,于是得方程组⎪⎪⎩⎪⎪⎨⎧=+=+=-=-4212522211y y x x y y x x . 解得⎪⎪⎩⎪⎪⎨⎧-=-=-=-=y y x x y y x x 2421522211,所以⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=y y y x x x X 245212,其中y x ,为任意实数.例16 解矩阵方程,C AX =其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=031334213A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=7577111793C . 解 由于0=A ,所以A 是不可逆矩阵,需要用元素法求解.设,222111⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=z y x z y x z yxX 则⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--7577111793031334213222111z y x z y x z y x,即 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++++-+-+-+-+-+-7577111793323334334334232323111212121212121z z y y x x z z z y y y xx x z z z y y y x x x .比较第一列元素得⎪⎩⎪⎨⎧=+=+-=+-73133432312121x x x x x x x x ,解得⎩⎨⎧-=-=9537121x x x x 同样,比较第二、三列元素可得对应方程组,分别解得7537,3535121121-=-=-=-=z z z z y y y y ,所以可得 ⎥⎦⎤⎢⎣⎡------=7573535953711111`1z z y y x x X ,其中111,,z y x 是任意实数. 总之,对于矩阵方程,当系数矩阵是方阵时,先判断是否可逆.如果可逆,则可以利用左乘或右乘逆矩阵的方法求未知矩阵,如果方阵不可逆或是系数矩阵不是方阵,则需要用待定元素法通过解方程确定未知矩阵.3.4 矩阵对角化方法3.4.1 讨论对于有n 个特征单根的n 阶方阵3.4.1.1 基本原理引理1 设A 是秩为r 的n m ⨯阶矩阵,且()n TE A−−−→−行初等变换⎪⎪⎭⎫ ⎝⎛*--n r n mr n rmP D )()(0 其中D 是秩为r 的行满秩矩阵,则齐次线性方程组0=AX 的一个基础解系即为矩阵P 所含的r n -个行向量),,2,1(r n i i -= ξ.引理2 矩阵A 的特征矩阵)(λA 经过一系列行初等变换可化为上三角形的λ-矩阵)(λB ,且)(λB 的主对角线上元素乘积的λ多项式的解为矩阵A 的全部特征根.引理3 对于数域P 上的n 阶方阵A ,若A 的特征多项式在P 内有n 个单根,则由特征向量构成的n 阶可逆矩阵T ,使得⎪⎪⎪⎪⎪⎭⎫⎝⎛=-n AT T λλλ211定理1 若数域P 上的n 阶方阵A 的特征多项式)(λf 在P 内有n 个单根,则A 可通过如下方法对角化:设()())()()(,)(λλλλλQ B E A A E A n TT T −−−→−-=行初等变换且)()1λB 为上三角形矩阵,则有方阵A 的特征根i λ即为)(λB 中主对角线上各个元素乘积的解;)2对于方阵A 的每一个特征根i λ,总有)(i B λ中零行向量所对应的)(i Q λ中的行向量i ξ与之对应.3.4.1.2举例说明例17 设⎪⎪⎪⎭⎫ ⎝⎛=210131012A ,问方阵A 是否可以化为对角形,若可以,求出其对角化后的方阵.解 ()⎪⎪⎪⎭⎫ ⎝⎛-------=100210010131001012)(λλλλE A T−−−−−−→−第一行与第二行互换⎪⎪⎪⎭⎫ ⎝⎛-------100210001012010131λλλ −−−−−−−−−→−-行上乘以第一行再加到第二)2(λ⎪⎪⎪⎭⎫ ⎝⎛---+-+----10021002125500101312λλλλλλ−−−−−−→−第二行与第三行互换⎪⎪⎪⎭⎫ ⎝⎛-+-+------02125501002100101312λλλλλλ −−−−−−−−−−→−+-行上乘以第二行再加到第三)55(2λλ⎪⎪⎪⎭⎫ ⎝⎛+----------5521)4)(2)(1(001002100101312λλλλλλλλ =())()(λλQ B由题意知)4)(2)(1(---λλλ=0⇒11=λ,22=λ,43=λ ,此时方阵A 有3个特征单根,故方阵A 可以化为对角形;将11=λ代入)()(λλQ B 和中知)(λB 的第三行为零,由定理1知)(λQ 的第三行向量)1,1,1(-即为属于1λ的特征向量,同理可知)1,2,1(),1,0,1(-分别为属于32λλ和的特征向量.于是可得⎪⎪⎪⎭⎫ ⎝⎛--=111201111T ,使得⎪⎪⎪⎭⎫ ⎝⎛=-4211AT T .3.4.2 讨论对于有特征重根的n 阶方阵对于有特征重根的方阵,可以通过上述方法将其化为上三角形矩阵,接着再对上三角形矩阵施行一系列初等变换将其化为对角形矩阵,这样就避免了上三角形矩阵中非零行向量可能不构成行满秩的情形. 3.4.2.1基本定理定理2 设TT A E A -=λλ)(,则()())()()(λλλP D E A T −−−→−初等变换且)(λD 为对角形矩阵,则有)1对于A 的每个特征根i λ,)(i P λ中与)(i D λ的零行对应的行向量即为属于i λ的特征向量;)2设s λλλ ,,21为A 的所有不同的特征根,重数分别为s r r r ,,21,则A 可以化成对角形⇔)(i D λ中的零行数目等于i λ的重数),,2,1(s i r i =.由此我们不难得到对于有特征重根的方阵化为对角形方阵的简单步骤如下:)1作()()())()()()()(λλλλλP D Q B E A T −−−→−−−−→−初等变换行初等变换,其中))(),(),(()(21λλλλn d d d diag D =,则A 的特征根恰为0)()()(21=λλλn d d d 的根;)2若A 的特征根全在P 内,且每个i λ有)(i D λ中零行数目等于i λ的重数,则A 可以化为对角形方阵,否则A 不可以化为对角形方阵;)3对于每个特征根i λ,在)(i P λ中取出与)(i D λ中零行对应的行向量),,,(21im i i P P P 得A属于i λ的特征向量且都是线性无关的. 3.4.2.2 举例说明例18 ⎪⎪⎪⎭⎫⎝⎛-=110111110)1A ; ⎪⎪⎪⎭⎫ ⎝⎛--=100112001)2B问方阵A 和B 是否可以化为对角形,若可以,试求出其对角化后的方阵.解 ()⎪⎪⎪⎭⎫ ⎝⎛------=10011101011100101)()1λλλλE A T−−−−−−→−第一行与第三行互换⎪⎪⎪⎭⎫ ⎝⎛------00101010111100111λλλ−−−−−−−−→−-行上乘以第一行再加到第二)1(⎪⎪⎪⎭⎫ ⎝⎛------0010111020100111λλλλ−−−−−−−→−行上乘以第一行再加到第三λ⎪⎪⎪⎭⎫ ⎝⎛-------λλλλλλλ0110110201001112 −−−−−−−−→−-二行上)乘以第三行再加到第(1⎪⎪⎪⎭⎫⎝⎛---------λλλλλλλ011011110111122−−−−−−−−−→−-三行上)乘以第二行再加到第(1λ⎪⎪⎪⎭⎫⎝⎛++------------112)1(001111010*******λλλλλλλλλ−−−−−−−−−→−-列上乘以第二列再加到第三)(2λ⎪⎪⎪⎭⎫ ⎝⎛++----------+--112)1(00111010100111222λλλλλλλλλ−−−−−−−−−−→−-+-列上乘以第一列再加到第三)1(2λλ⎪⎪⎪⎭⎫ ⎝⎛++----------112)1(0011101010001122λλλλλλλ−−−−−−→−第二行加到第一行上⎪⎪⎪⎭⎫⎝⎛++------------112)1(001110101100122λλλλλλλλ())()(λλP D =由题意知0)1(2=-λλ⇒01=λ,)(12二重=λ,因为)(2λD 中零行数目≠1等于2λ的重数,故A 不可以化为对角形方阵.)2 ()⎪⎪⎪⎭⎫ ⎝⎛--+-=100110010010001021)(λλλλE A T2014届数学与应用数学专业毕业(论文)第 23 页 共 24页−−−−−−→−第二行与第三行互换⎪⎪⎪⎭⎫ ⎝⎛+---010*********001021λλλ −−−−−−−−−→−+行上乘以第二行再加到第三)1(λ⎪⎪⎪⎭⎫ ⎝⎛+----1101001001100010212λλλλ −−−−−−−−−→−-列上乘以第二列再加到第三)1(λ⎪⎪⎪⎭⎫ ⎝⎛+----110100100010001)1(2212λλλλ −−−−−−−−→−-列上乘以第一列再加上第三)2(⎪⎪⎪⎭⎫ ⎝⎛+---1101001000100010212λλλ −−−−−−−→−行上乘以第二行再加到第一2⎪⎪⎪⎭⎫ ⎝⎛+---1101001000102010012λλλ())()(λλP D =. 由题意知0)1)(1(2=--λλ⇒)(11二重=λ,12-=λ,此时)(1λD 中零行数等于=21λ的重数,故B 可以化为对角形方阵;将11=λ代人)()(λλP D 和中知)(λD 的第一行和第三行为零,由定理2知)(λP 的第一行向量)2,0,1(和第三行向量)2,1,0(即为属于1λ的特征向量,同理可知)0,1,0(为属于2λ的特征向量.由此可知⎪⎪⎪⎭⎫ ⎝⎛=022110001T 使得⎪⎪⎪⎭⎫ ⎝⎛-=-1111BT T .结 论通过以上对矩阵的学习,我们知道,想要在学习过程中灵活应用矩阵思想,首先要理解矩阵思想,在此基础上,遇到难解的数学问题,能发现矩阵是可以解决此类问题的关键,最后能正确无误的利用矩阵思想把数学问题得以解决.矩阵是代数特别是线性代数的一个主要研究对象,他对于研究矩阵的相关运算、解线性与非线性方程组、特征值和特征向量的求解方法、对角化及二次型矩阵、求解矩阵高次幂等重要问题都有极为广泛的应用.杨灿:矩阵及其应用参考文献[1]李志慧,李永明.高等代数中的典型问题与方法[M].科学出版社,2008.205-211[2]王萼芳,石生明.高等代数(第三版).高等教育出版社[3] 张禾瑞.高等代数(第五版)[M].北京:高等教育出版社,2007[4] 吕林根,许道子.解析几何[M].北京:高等教育出版社,2006[5] 许以超.线性代数与矩阵[M].北京:高等教育出版社,1992[6] 李师正.高等代数解题方法与技巧[M].北京:高等教育出版社,2004[7] 徐仲,张凯院,陆全,冷国伟.矩阵论简明教程[M].北京:科学出版社,2005[8] 贾美娥.矩阵的秩与运算的关系[J].赤峰学院学报,2010,26(9):3-4[9] 钟成义,肖宏儒.方阵秩与零特征值代数重数相关性探讨[J].高等数学研究.2009,12(1):96-97[10] 史明仁. 线性代数600证明题详解[M]. 北京科学技术出版社.1985[11] 徐德余.高等代数(第二版)[M].四川大学出版社.2005:175-178[12] 丘维声. 高等代数[M]. 北京: 高等教育出版社, 1996[13] 赵树嫄. 线性代数(第三版[M]). 北京: 中国人民大学出版社, 2006[14] 程云鹏.矩阵论[M].第二版.西安:西北工业大学出版社,2002[15] 赵树塬.线性代数[M].北京:中国人民大学出版社,1997[16] 李君文.线性代数理论与解题方法[M].长沙 :湖南大学出版社,2002致谢从上学期选题、收集资料到这学期写开题报告,完成初稿,到定稿,期间几个月历经喜悦、聒噪、痛苦、彷徨,在写论文时心情如此复杂,到今天随着论文的完成,都落下了帷幕.在此论文撰写过程中,要特别感谢我的导师向以华老师的指导与督促,同时感谢他的谅解与包容.没有老师的帮助也就没有今天的这篇论文.求学历程是艰苦的,但又是快乐的.感谢我大学所有教过的老师,谢谢他们在这四年中的教诲.在这四年的学期中结识的各位生活和学习上的挚友让我得到了人生最大的一笔财富.在此,也对他们表示衷心感谢.本文参考了大量的文献资料,在此,向各学术界的前辈们致敬!第24页共24 页。
(完整版)幂零矩阵的性质及应用毕业设计
证明:设,
于是 ,
所以是幂零矩阵
性质2.2.3设是阶幂零矩阵,则,均为幂零矩阵。
证明:因为为幂零矩阵,,使得,因为
所以,均为幂零矩阵。
性质2.2.4幂零矩阵的行列式值为零。
证明:设是阶幂零矩阵,则存在一个自然数k,使,由行列式性质得
所以
性质2.2.5与幂零矩阵相似的矩阵是幂零矩阵
2.幂零矩阵的相关概念及简单性质
为了叙述的需要,我们首先引入幂零矩阵的相关概念.
2.1幂零矩阵的相关概念
定义2.1.1令为阶矩阵,若存在正整数,使,则称幂零矩阵。也称为阶幂零矩阵。如为2阶幂零阵,则。
定义2.1.2若为幂零矩阵,满足的最小正整数称为的幂零指数。显然,阶零矩阵是特殊的幂零矩阵且其幂零指数为1。
定义2.1.6形为
,
其中,
由阶数为的若干个若尔当块组成的准对角称为若尔当形矩阵.
定ቤተ መጻሕፍቲ ባይዱ2.1.7设为阶方阵,的首项系数为1的最低次的化零多项式称为的最小多项式。
2.2关于幂零矩阵的一些简单性质
由上述所描述的有关幂零矩阵的定义,可以得出一些幂零矩阵的几条简单性质。
性质2.2.1幂零矩阵都不可逆。
证明:设是任一阶幂零矩阵,则,使,假设可逆,则,于是,故也可逆,这与矛盾。
由幂零矩阵的定义,是幂零矩阵。
借助这个结论,要证明幂零矩阵的伴随矩阵还是幂零矩阵就很方便了,证明如下:
由这个充要条件,可以得出以下的几个推论:
推论3.1.1设是阶幂零矩阵,则为幂零矩阵。
关键词:幂零矩阵;特征值;若尔当形
Abstract
Matrix in important tool to research problem,
矩阵指数函数及其应用
根据代数基本定理,在复数域上可以证明:
性质 2.1.3 设 其最小多项式且有
, λ , λ , … λ 为 A 的 r 个互不相同的特征值,ΨA λ 为
ΨA λ λ λ λ λ … λ λ ,
其中 d 1 i 1,2 … , r , ∑ d m.
如果函数f x 具有足够多阶的导数值,并且下列 m 个值(称 f(x)在影谱上的值)
2
中国矿业大学(北京)07 级本科生设计(论文)
2 矩阵函数定义及矩阵指数函数的性质
2.1 矩阵函数定义
定义 2.1.1 给定矩阵, 如果多项式
pλ α λ α λ
αλ α
满足p
0,则称p λ 是 A 的化零多项式.
定义 2.1.2 在 A 的化零多项式中,次数最低且首项系数为 1 的化零多项式 称为 A 的最小多项式,记为 ( ).
2 矩阵函数的定义及矩阵指数函数的性质………………………………………………2 2.1 矩阵函数定义………………………………………………………………………2 2.2 矩阵指数函数的性质………………………………………………………………4
3 矩阵指数函数的三种计算方法…………………………………………………………8 3.1 第一种方法…………………………………………………………………………8 3.2 第二种方法…………………………………………………………………………12 3.3 第三种方法…………………………………………………………………………16 3.4 三种方法的比较……………………………………………………………………19
院长签字:
指导教师签字:
中国矿业大学(北京)本科生毕业设计(论文)指导教师评阅书
学院 理学院 专业 信息与计算科学 班级 1 学生姓名 郭世钢 题目:矩阵指数函数及其应用 专题题目: 指导教师评语:
正定矩阵的若干应用 毕业设计 精品
本科毕业设计(论文)( 2011届 )题目:正定矩阵的若干应用学院:数理与信息工程学院专业:数学与应用数学学生姓名:学号:指导教师:职称:副教授合作导师:职称:完成时间:20 年月日成绩:本科毕业设计(论文)正文目录摘要 (1)英文摘要 (1)1 引言 (2)1.1 矩阵理论的发展历史 (2)1.2 正定矩阵的发展与地位 (3)2 正定矩阵 (4)2.1 正定矩阵的定义 (4)2.2 正定矩阵的相关理论 (4)2.2.1 正定矩阵的性质 (4)2.2.2 正定矩阵的相关定理 (7)2.2.3 正定矩阵的判别方法 (10)3 正定矩阵应用 (12)3.1 正定矩阵的相关理论推广 (12)3.1.1 广义正定矩阵 (12)3.1.2 准正定矩阵 (144)3.1.3 Schur定理与华罗庚定理的推广 (16)3.1.4 Ky Fan等著名不等式的推广 (16)3.2 正定矩阵在实际问题中的应用 (17)3.2.1 二次型理论的应用 (17)3.2.2 仿射变换 (19)参考文献 (22)正定矩阵的若干应用摘要:正定矩阵是矩阵理论中的一类重要矩阵, 且在多个不同领域内均有重要作用. 本文回顾了正定矩阵的发展历史以及性质, 主要探讨了它的若干应用, 其中包含正定矩阵的理论推广和实际应用等问题.关键词:正定矩阵;性质;理论;推广;应用Several applications of Positive DefiniteMatrixesAbstract:Positive definite matrixes are an important class of matrixes in the matrix theory, which are widely used in different fields. In this paper, we first recall the history of development, then some basic properties of positive definite matrixes, and we mainly discuss several applications of positive definite matrixes, including the extension and practical applications of positive definite matrixes.Key Words:positive definite matrix;properties;theories;extend;application1 引言矩阵是数学中的一个重要的基本概念, 是代数学的一个主要研究对象, 也是数学研究和应用的一个重要工具. 矩阵理论是数学的一个重要的分支, 它不仅是一门基础学科, 也是最具实用价值和广泛应用的数学理论, 现已成为处理有限维空间形式和数量关系的强有力工具.正定矩阵作为一类常用矩阵, 其在数学学科和其他科学技术领域的应用也非常广泛, 因此它的性质定理以及应用问题一直倍受关注, 而在实际生活中也经常出现有关正定矩阵的应用, 如线性规划、二次型理论解决二次曲线问题等. 尽管个别理论已为人们所熟知, 但缺乏系统性的整理.本文对正定矩阵的研究主要集中在对正定矩阵其性质的推广和应用上, 包括理论和实际的应用. 结合当前对正定矩阵已有的成形研究, 从二次型理论入手,弄清其应用及推广, 并研究正定矩阵的仿射变换解决不规则二次曲线问题等实际问题的应用.1.1 矩阵理论的发展历史“矩阵(Matrix)”术语是由西尔维斯特创用并由凯莱首先明确其概念的. 他为了将数字的矩形阵列区别于行列式而发明了这个述语. 而实际上, 矩阵这个课题在诞生之前就已经发展的很好了. 早在公元前1世纪, 矩阵形式解方程组在中国古代数学著作《九章算术》中已相当成熟, 那个时候矩阵只是被看作一种排列形式来解决实际问题, 并没有建立起独立的矩阵理论. 从18世纪末到19世纪中叶, 这种排列形式在求解线性方程组和行列式计算等问题中应用日益广泛, 矩阵思想才得到进一步的发展. 19世纪50年代, 西尔维斯特引入“矩阵”一词来表示“一项由m行n列元素组成的矩形阵列”或“各种行列式组”, 凯莱作为矩阵理论的创立者, 首先为简化记法引进矩阵, 然后系统地阐述了矩阵的理论体系. 随后, 弗罗伯纽斯等人发展完善了矩阵的理论体系形成了矩阵的现代理论.矩阵思想的萌芽由来已久, 早在公元前1世纪中国的《九章算术》就己经用到类似于矩阵的名词. 1748年, 瑞士数学家欧拉在将三个变数的二次型化为标准型时, 隐含地给出了特征方程的概念. 1773年, 法国数学家拉格朗日在讨论齐次多项式时, 引入了线性变换. 1801年德国数学家高斯在《算数研究》中,将欧拉与拉格朗日的二次型理论进行了系统的推广, 给出了两个线性变换的复合, 而这个复合的新变换其系数矩阵是原来两个变换的系数矩阵的乘积. 另外, 高斯还从拉格朗日的工作中抽象出了型的等价概念, 在研究两个互逆变换的过程中孕育了两个矩阵的互逆概念. 1826年, 柯西在《微积分在几何中的应用教程》中讨论了二次型束的特征根使束的行列式为零的情况, 证明了当其中一个二次型对变数的所有非零实数值是正定时, 束的特征根全为实数. 18世纪中期, 数学家们开始研究二次曲线和二次曲面的方程简化问题, 即二次型的化简. 在这一问题的研究中, 数学家们得到了与后来的矩阵理论密切相关的许多概念和结论. 从18世纪中期到19世纪初, 数学家们在研究二次型的过程中涉及到大量的线性变换, 得到了许多重要概念和结论. 由于二次型和线性变换均可以使用矩阵来表示, 所以这些概念和结论也就可以自然而然地移植到矩阵理论之中. 因此二次型理论是矩阵思想得以孕育的重要源泉之一. 与此同时, 这种排列形式在线性方程组和行列式计算中应用日益广泛, 行列式等理论的发展提供了矩阵发展的条件, 矩阵概念由此产生, 矩阵理论得到系统的发展. 20世纪初, 无限矩阵理论得到进一步发展.矩阵及其理论现已广泛地应用于现代科技的各个领域, 矩阵也是高等代数中最重要的内容之一, 矩阵的初等理论现已作为高中数学的选修内容. 由此看出, 矩阵理论在科技发展以及教育教学中的重要地位.1.2 正定矩阵的发展与地位矩阵的正定性源于二次型与Hermite型的研究, 最初只限于在实对称矩阵或Hermite矩阵中讨论. 但它要求矩阵是实对称的或是Hermite的. 刚开始, 人们通过将非对称矩阵加以“对称化”来解决与原来矩阵相关的问题. 1936年, KB用正定阵乘以原矩阵使其乘积成为对称阵的方法, 解决了概率论中的一些问题. 1937年, Johnson在其博士论文中研究了方阵A的对称化AA'+是正定阵的某些不等式. 随着实际应用的需要, 1970年, Johnson给出n阶实矩阵A正定的定义, 把正定性的研究推广到未必对称的矩阵中, 并对这类矩阵的性质、不等式给予研究, 这些结论应用于许多领域. 国内近10年来出现了关于这方面研究的大量文章.将正定矩阵推广到广义正定矩阵之后, 得出了许多正定矩阵的相关理论. 从定义出发, 研究结论包括对称矩阵非对称矩阵、准正定矩阵、次正定矩阵、广义正定矩阵;从正定矩阵的性质出发, 研究结论有:矩阵的三角分解、矩阵跡的问题, 还有相关Hadamard积、Kronecker积、Hermite矩阵方程、Hamilton四元数理论的应用等问题;从矩阵相关理论出发, 得出:惯性定理, Schur定理, 华罗庚定理, Minkowski及Ky Fan不等式, 扩大了Minkowski不等式的指数范围等;同时也包括这些性质的推广与实际应用的探讨.研究矩阵的正定性, 在数学理论或应用中具有重要意义, 是矩阵论中的热门课题之一. 正定矩阵具有广泛的应用价值, 是计算数学、数学物理、控制论等领域中具有广泛应用的重要矩阵类, 其应用引起人们极大的研究兴趣. 在非对称矩阵领域研究正定矩阵突破了其本身定义的限制, 获得了丰富的研究成果, 并且得到了广泛的应用, 如线性规划的最优算法及线性回归模型结构、控制论、矩阵方程论、组合矩阵等. 它的研究成果, 使整个矩阵理论体系得以完善.2 正定矩阵2.1 正定矩阵的定义定义 1 (实正定矩阵)设),,,(21n x x x f 是一个实二次型, 若对任意的一组不全为零的实数n c c c ,,21 都有0),,,(21>n x x x f , 则称),,,(21n x x x f 是实正定二次型, 它所对应的对称矩阵为正定对称矩阵, 简称正定矩阵.定义2 n 阶实对称矩阵A 称为正定矩阵, 如果对于任意的n 维实非零列向量),,,(21n x x x X =都有0>'A X X . 正定的实对称矩阵A 简称为正定矩阵.由此可知, 研究矩阵的正定问题, 可以转化为研究其所对应二次型的正定问题.注记:若未作特别说明, 这里所讨论的矩阵均为实矩阵.2.2 正定矩阵的相关理论2.2.1 正定矩阵的性质对于实方阵来说, 首先具备以下性质:性质1 设矩阵A 为n 阶实方阵, 则下列命题等价:1) A 是正定矩阵;2) 1-A 是正定矩阵;3) A '是正定矩阵;4) A A '+是正定矩阵;5) 对任意n 阶可逆矩阵P , AP P '是正定矩阵;6) A 的各阶主子矩阵是正定矩阵.证2) 若A 是正定的, 则存在实可逆矩阵C 使C C A '=.∴)()(1111'='=----C C C C A ,∵C 可逆, ∴1-C 也是实可逆矩阵. ∴有1-A 也是正定矩阵.充分性:若1-A 是正定矩阵, 则1111)()(----'='=C C C C A .∵C C C C A A '='==----1111))(()(,∴A 是正定的.3) 同2)的证明方法.由性质4)可得如下推论:推论1 若B A ,都是正定矩阵,则B A +也是正定矩阵.证 因为B A ,都是正定矩阵,所以BX X AX X ','都是正定二次型,于是有BX X AX X X B A X '')('+=+也一定是正定二次型,所以B A +是正定矩阵.性质2 设n 阶实对称矩阵A 为正定矩阵, 与下列命题互为充要条件:1) A 的正惯性指数等于A 的维数n ;2) A 合同于单位矩阵E ;3) 存在满秩阵C , 使C C A '=成立;4) A 的n 个特征值全为正值;5) 存在满秩阵P , 使AP P '成对角线元素皆正的对角阵D ;6)存在对称正定阵B , 使2B A =;证3)必要性:若A 是正定矩阵, 则A 合同于E∴存在实可逆矩阵C , 使C C EC C A '='=充分性:若C C A '=, C 是实可逆矩阵, 对0,0≠≠∀CX X , 则0)()(>'=''='CX CX CX C X AX A所以, A 是正定的.4)设λ为矩阵A 的任一特征值, X 为与其相应的特征向量, 则有x Ax λ=, 因而有0),,,(21>'='=X X AX X x x x f n λ , 因而0>λ.推论2 正定矩阵的实特征值都是正的, 而复特征值一定具有正的实数部分. 注记:这个定理的逆命题不一定成立, 即复特征值具有正的实属部分而实特征值全是正的实方阵不一定是正定的, 例如1400010000270012A ⎛⎫ ⎪ ⎪= ⎪ ⎪-⎝⎭, 它的特征值为27i ±, 实特征值是1,1, 而1222221234122344341400010055(,,,)(2)32()0027220012x x x x x x x x x x x x x x ⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪=+-++- ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭, 容易看出, A 不是正定的.性质31)[1] A 的所有顺序主子式大于零;2) 正定矩阵A 的主子式全大于零;3) 正定矩阵的主对角线上元素必全大于零.证2) 设⎪⎪⎪⎭⎫ ⎝⎛=nn n n n n a a a a a a a a a A 212222111211, 而111212122212i i i i i ik i i i i i ik K iki iki ikik a a a a a a A a a a ⎛⎫ ⎪ ⎪= ⎪ ⎪⎝⎭ , n i i i k <<≤≤≤ 211 为A 的任一K 阶主子式, k A 为所对应的K 阶主子式的行列式. 由于A 是正定矩阵, 故二次型AX X x x x f n '=),,,(21 . 对任意不全为零的实数n x x x ,,,21 都有0),,,(21>n x x x f , 从而, 对不全为零的实数0)0,,,,0,,0(21> ik i i x x x f (即在),,,(21n x x x f 中除ik i x x 1外余者取0).对于变量ik i x x 1矩阵k A 的二次型0)000(),,(121>= ik i ik i i x x f x x x g ,故g 是正定二次型. 因而k A 是正定矩阵, 故0>k A3) 当实对称矩阵()ij A a =是正定的时, 它所确定的二次型11'()n nij i j ij ji i j X AX a x x a a ====∑∑必为正定二次型. 假设0ii a ≤,取110,,0,i x x -==11,0,,0i i n x x x +=== , 代入上式得'0ii X AX a =≤, 这即与'X AX 是正定的相矛盾, 所以只有0(1,2,,)ii a i n >= .性质4 正定矩阵的元素的绝对值最大者一定是主对角线上的元素.证 设()ij A a =是正定矩阵, 其中元素()ij a i j ≠的绝对值为最大, 则2ij ii jj a a a ≥, 由性质3可知, A 的一切主子式都大于零, 从而有20iiij ii jj ij ij jj a a a a a a a -=>, 即ij ii jj a a a <, 这与假设矛盾, 故正定矩阵A 中元素的绝对值最大者必定是主对角线上的元素.注记:这个结论常用于判定某些实对称矩阵不是正定的矩阵. 因为只要有一个非主对角线上的元素绝对值不小于主对角线上元素的绝对值的最大者, 那么这个实对称矩阵必不是正定矩阵.性质5[6] 设()n A M R ∈实正定, 则1) 对任意,0,()0n x R x x Ax x R A x ''∈≠=>;2) A 的所有主子阵k A 及/k A A 实正定;3) A 的所有主子阵k A 的特征根()k A λ满足:0min (())Re ()max (()),Im ()max (())k k R A A R A A S A λλλλλ<≤≤≤;4) A 的所有主子式行列式大于0, 特别地0A >;5) 对任意1(,),0n n X x x R X '=∈≠ , 令1(,,)n AX z z '= (或1(,,)n A X z z ''= ), 则存在1k n ≤≤, 使得0k k x z >;6) 存在()n P M C ∈非奇异, 使得**11(),()[,,,n P R A P E P S A P diag bi bi ==-,,0,0]s s b i b i - , 其中(())2,0,k r S A s b => 21,,,1k s i ==- .7) 若()n B M R ∈实正定, 则[0,1],(1)t tA t B ∀∈+-实正定, 即n 阶实正定矩阵集合为一凸集;8) A 关于任一顺序主子阵的Sylvester 矩阵实正定.证1) 只需证'()0,n X S A X X R =∀∈, 而'()('())''()X S A X X S A X X S A X ==- 故'()0X S A X =.先证3) ∵min (())min (()),max (())max (())k k R A R A R A R A λλλλ≥≤, 对()iS A -, 此为Hermite 矩阵, 故只需证k A A =时结论成立即可. 令()A λ为A 的任意特征值, n X C ∈为响应的单位特征向量, 则**()(,)()()A X AX X R A X X S A X λ==+, 易知*()X S A X 为纯虚数, 故***Re ()(),Im ()()((()))A X R A X i A X S A X i X iS A X λλ===-, 由Courant-Fischer 定理直接得证.4) 可由3) 直接得到.下证2)由于我们有*0/k k A A A A ⎛⎫−−−→ ⎪⎝⎭合同变换, 故*0/k k A A A ⎛⎫ ⎪⎝⎭实正定. ,0,k x R x ∀∈≠令*(',0,,0)',''00/k n k k A x x R x x X A X A A ⎛⎫=∈=> ⎪⎝⎭, 故k A 实正定, 同理可证/k A A 实正定.5) 可由'0(''0)x Ax x A x >>或直接得到.6) ()R A 实对称正定, 存在()n Q M R ∈, 使得'(),'(())n Q R A Q E Q iS A Q =-为Hermite 矩阵, 存在酉阵U 使得*1'(())[,,],n i U Q iS A QU diag d d d R -=∈ . 故*11'()[,,,,,0,,0],(())2,0,1,,s s k U Q S A QU diag bi bi b i b i r S A s b k s =--=>= . 令P QU =即可得证.7) 利用((1))(()(1)())R tA t B t R A t R B +-=+-便可证明.8) 记此顺序主子式为k A , 则A 关于k A 的Sylvester 矩阵为/,k k A A A ⋅0k A >,且/k A A 实正定, 故有/k k A A A ⋅实正定.2.2.2 正定矩阵的相关定理定理1 实对称矩阵A 是正定矩阵的充分而且必要条件是对于任意的n 维实非零列向量10n x X x ⎛⎫⎪=≠ ⎪ ⎪⎝⎭, 二次型'X AX 是正定二次型.定理2[3] Hermite 矩阵是半正定的, 当且仅当它的所有特征值都是非负的, 它是正定的, 当且仅当它的所有特征值都是正的.定理3[6] n n C ⨯中所有正定矩阵构成的集合的边界点是半正定矩阵, 但不是正定矩阵.定理 4[14] (Hadamard 不等式)正定矩阵()n n ij A a C ⨯=∈的行列式不超过对角元素之积, 即等式1det nii i A a =≤∏成立当且仅当A 是对角矩阵.定理5[14](Oppenheim 不等式)设(),()n n ij ij A a B b C ⨯==∈是半正定矩阵, 则1(det )det()nii i A b A B =≤∏ .定理 6[14] (Minkowski 不等式)设,n n A B C ⨯∈是半正定矩阵, 则:111[det()](det )(det )n n n A B A B +≥+定理7 设n 阶方阵A S K =+是正定的, 其中S 和K 分别是A 的对称分量和反对称分量, 则det det det A S K ≥+.证 因为A S K =+是正定的, 可知存在可逆方阵P 使得:'S P P =11000'000s s a a a K P P a ⎛⎫ ⎪- ⎪ ⎪ ⎪⎪= ⎪- ⎪⎪⎪ ⎪ ⎪⎝⎭ 11111'111s s a a a A P P a ⎛⎫ ⎪- ⎪ ⎪ ⎪⎪= ⎪- ⎪⎪⎪ ⎪ ⎪⎝⎭, 其中120s a a a ≥≥≥> , 因此2222211det (1)(1)(det )(1)s s A a a P a a =++≥+ 2(det )P , 如果A 的反对称分量K 不可逆, 则det 0K =, 而2det (det )S P =, 因此22221det (1)(det )(det )det det s A a a P P S K ≥+≥=+ . 如果A 的反对称分量K 可逆,则2n s =, 而且1,,s a a , 都不为零, 并且22221det (det ),det (det )s K a a P S P == , 所以有det det det A S K ≥+.2.2.3 正定矩阵的判别方法在研究正定矩阵的时候, 会出现判断被研究矩阵是否正定的问题. 结合性质, 正定矩阵的判定方法有很多, 随着研究的深入, 方法也不断改进, 以下罗列了几个相关判定方法:1) 与正定矩阵合同的矩阵一定是正定矩阵.事实上由合同的传递性及正定矩阵都与单位矩阵合同可知结论明显成立. 2) 正定矩阵的逆矩阵必为正定矩阵.因为正定矩阵与单位矩阵合同, 所以存在可逆矩阵P , 使得''A P EP P P ==, 取逆矩阵11111()()'(')A P E P P E P -----==, 记1()Q P -'=, 即有1'A Q EQ -=, 则1A -与单位矩阵合同, 所以1A -是正定矩阵.3) 正定矩阵的和仍是正定矩阵.事实上若A 与B 是同阶正定矩阵, 则对于任意的非零实列向量12(,,,)0n C c c c =≠ , 必有'0C AC >, 且'0C BC >, 从而'()C AB C+= ''0C AC C BC +>, 所以A B +是正定的.4) 正定矩阵的任何主子式阵必为正定矩阵.假设()ij A a =是一个n 阶正定矩阵, 它的k 阶主子式阵1111(1),k k k kk a a k n A a a ⎛⎫⎪≤≤= ⎪ ⎪⎝⎭, 又由正定矩阵性质可知0k A >, 从而可知A 的任何主子式阵一定是正定的.5) 对于任何的实对称矩阵A , 必有实数0,0αβ>>, 使得E A α+与E A β+是正定矩阵.∵实对称矩阵A 的特征根都是实数, 不妨记其中绝对值最大的一个特征根为0λ, 只要取0βλ>, 即可使E A β+是正定矩阵.这是因为假设Q 是正交矩阵, 可使1'n Q A Q λλ⎛⎫ ⎪= ⎪ ⎪⎝⎭ , 则有11'()''n n Q E A Q Q EQ Q AQ βλβλβββλβλ+⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪+=+=+=⎪ ⎪ ⎪⎪ ⎪ ⎪+⎝⎭⎝⎭⎝⎭,其中由于0(1,2,,)i n βλ+>= , 可知E A β+是正定矩阵.当取1αβ=时, 则10,()E A E A ααββ>+=+是正定矩阵.6) 假设,A B 都是正定矩阵, 并且AB BA =, 则AB 也必为正定矩阵.AB 的特征根都大于零, 当AB BA =时, ()'''AB B A BA AB ===说明AB 又是对称的, 从而可知AB 是正定的.3 正定矩阵应用正定阵具有广泛的应用, 但被局限在对称阵范围内, 随着应用的需要和研究的深入, 加速了突破这一限制的进程. 国内外不少学者研究了多种未必对称的较为广义的正定阵, 获得了丰富的研究成果, 其成果得到了广泛的应用, 但仍不能满足应用上的需要和达到理论上的完善, 正定矩阵的相关性质、定理以及证明相关的著名不等式等问题都将日趋完善.3.1 正定矩阵的相关理论推广 3.1.1 广义正定矩阵按照正定矩阵的严格定义, 其要满足该矩阵是对称的. 而广义正定矩阵突破这一限制, 来进行研究. 对于广义正定矩阵, 有过一系列的推广, 这里给出其部分定义、性质以及相关应用的问题.定义3 设n n A R ⨯∈, 矩阵A 称为正定的, 如果10n x R ⨯∀≠∈, 有'0X AX >. 这种正定矩阵称为J -正定矩阵, n 阶J -正定矩阵的全体记为J P .定义4 设n n A R ⨯∈, 矩阵A 为正定的, 若对任何10n x R ⨯≠∈, 都有正对角阵存在正对角矩阵0x D D =>, 使'0x DAx >, 则称A 为D -正定矩阵, 记为D A P ∈.定义5 设n n A R ⨯∈, 矩阵A 称为正定的, 如果存在S S P +∈, 使得对1120(,,,)n n X x x x R ⨯∀≠=∈ , 有'0X SAX >. 这种正定矩阵称为X -正定矩阵, n阶X -正定矩阵的全体记为X P .定义6 设n n A R ⨯∈, A K H =+, 其中n K S +∈, n H S -∈, 如果S K P +∈, 那么称A 为亚正定矩阵. 这里把n 阶亚正定矩阵全体记为TP (K 为A 的对称分量, H 为A 的反对称分量, 并且这种分解式是唯一的).定义7 设n n A R ⨯∈, 矩阵A 称为正定的, 如果存在J S P ∈使得1n O X R ⨯∀≠∈, 有'0X SAX >, 这里称这种正定矩阵为Y -正定矩阵, n 阶Y -正定矩阵的全体记为YP . 定义8 设n n A R ⨯∈, 矩阵A 称为正定的, 如果存在n n B R ⨯∈, 且det 0B >, 使得1n O X R ⨯∀≠∈, 有'0X SBX >. 这里把这种正定矩阵称为M -正定矩阵, n 阶M -正定矩阵的全体记为M P .定义9 设n n A R ⨯∈, 矩阵A 称为正定的, 如果存在B n n R ⨯∈, 'B B =, 且det 0B >, 使得1n O X R ⨯∀≠∈, 有'0X BAX >, 这里把这种正定矩阵称为N -正定矩阵, n 阶N -正定矩阵的全体记为N P .定义10 设n n A R ⨯∈, 若1n o X R ⨯∀≠∈, 都存在n 阶可逆阵n nX n Q Q R ⨯=∈使T X X Q AX O >, 则称A 为n 阶准正定矩阵. 若X Q Q =与X 无关, 则称A 为n 阶Q -准正定矩阵, 记为{}1|0,0n n n T Q A P A R X R X QAX ⨯⨯∈=∈∀≠∈>, 否则记为Q A P ∈ .这里进一步对定义9所给出的N -正定矩阵的性质定理进行进一步讨论.显然, N A P ∈的充分必要条件是J BA P ∈, 接下来进一步探讨N P 中矩阵的某些性质. 为了方便表达我们记{}|det 0;'n n G A R A A A ⨯=∈>=.定理8[15](1) N A P ∈当且仅当存在Z G ∈, 以及J B P ∈, 使得A ZB =; (2)N A P ∈当且仅当存在Z G ∈, 以及J B P ∈, 使得A BZ =.证 (1)必要性 因为N A P ∈, 1Z G ∃∈, 使得1n O X R ⨯∀≠∈有10XZ AX ->, 则1J B Z A P =∈, 于是得11A Z G -=∈, 则得到A ZB =.充分性 由于J B P ∃∈以及Z G ∈, 使得A ZB =, 则1J Z A B P -=∈, 所以1n O X R ⨯∀≠∈, 有10XZ AX ->, 而11Z G -∈, 于是得N A P ∈, 综上所述, 即得证.(2)必要性 因为N A P ∈, 所以1Z G ∃∈, 使得1n O X R ⨯∀≠∈有10XZ AX ->, 则1(')'0X Z AX >, 于是1'111''0X Z Z A Z X ->, 即1'1'11()''()0Z X Z A Z X ->, 由于1Z G ∈,则11Z G -∈, 由此可得'N A P ∈. 又'Z G ∃∈, 以及J B P ∈, 使得'''A Z B =, 于是得到A BZ =, 其中J B P ∈, Z G ∈,充分性 由于Z G ∃∈以及J B P ∈, 使得A BZ =, 所以对1n O X R ⨯∀≠∈, 有1(')'0X Z AX >, 则11111()''()'(')''0Z X Z A Z X X Z Z AZ X X AZ X -----==>, 由定义可知N A P ∈.定理9[15] 设n n A R ⨯∈, 那么N A P ∈当且仅当Z G ∃∈使得'S ZA A Z P ++∈.证 必要性 由于N A P ∈有定理1(1)可知Z G ∃∈使得J ZA P ∈, 所以有'J A Z P ∈, 则'J A Z ZA P +∈, 而(')''A Z ZA A Z ZA +=+, 于是有'S A Z ZA P ++∈.充分性 由于Z G ∃∈使得'S A Z ZA P ++∈, 于是1n O X R ⨯∀≠∈有'(')0X A Z ZA X +>, 则'(')'0X ZAX X ZAX +>, 于是'0X ZAX >, 故N A P ∈.以上N -正定矩阵的几个充分必要条件还可以推导出几个N -正定矩阵的若干性质, 下面我们给出性质:性质6[15]1) 若N A P ∈, 则(1)det 0A >;(2)'N A P ∈;(3)1N A P -∈. 2)设非奇异矩阵,,,n n N P R B G PB BP A P ⨯∈∀∈=∈, 则'N PAP P ∈ 3)若N A P ∈, 且A 的QR 分解为11A Q R =, 11N R Q P ∈.4)若,,N A P B G BA K H ∈∈=+, 其中A 是实对称正定矩阵, H 是实反对称矩阵, 则11det det det A B K B H --≥+.证 由于,N A P B G ∈∈, 于是得J BA P ∈, 而BA K H =+, 由定理可知:det det det BA K H ≥+, B 可逆(B G ∈), 则det 1det (det det )det det BA A K H B B=≥+, 所以有11det det det A B K B H --≥+.3.1.2 准正定矩阵这里用E 表示n 阶单位阵;J 表次对角线元全为1, 其余元全为0的n 阶方阵;,m n m n nR R ⨯⨯分别表示m n ⨯实矩阵集与n 阶实可逆矩阵集;n S +表n 阶实对称正定阵集;,A B A B ⊗ 分别表示矩阵A 、B 的Hadamard 积与Kronecker 积.对于定义10, 当Q E =时, E P 便是正定阵集或亚正定阵集, 并有n E S P +⊂;当Q D =(n 阶正对角阵)时, D P 便是广义正定阵集;当nQ S S +=∈时, S P 便是由矩阵做进一步推广后的广义正定阵集;当Q H =(n 阶实对称可逆阵)时, H P 便是非对称广义正定阵集;当Q J =时, J P 便是次亚正定阵集;当Q JD =(n 阶次对角阵)时,JD P 便是准次正定阵集;当Q JS =(n 阶实次对称次正定阵)时, JS P 便是广义次正定阵集. 因此准正定矩阵将各类实正定阵与各类实次正定阵统一了起来, 并有,n I D S H Q J JD JS Q S P P P P P P P P P +⊂⊂⊂⊂⊂⊂⊂⊂.定理10[13] 设n nn Q R ⨯∈, 则Q E A P QA P ∈⇔∈.证 由定义1知:10,0n Q E A P X R X QAX QA P ⨯'∈⇔∀≠∈>⇔∈.推论3 设n nn Q R ⨯∈,Q A P ∈,则E QA P∈为完全主正阵, 因而0,Q A A >可逆. 定理11[13] 设n nn Q R ⨯∈, 则Q n A P QA A Q S +''∈⇔+∈.证 由定理10知:[()]/2T Q E n n A P QA P QA QA S QA A Q S ++''∈⇔∈⇔+∈⇔+∈.在定理2中, 取Q S =为实对称正定阵, 即:定理12[13] 设n nn Q R ⨯∈, 则1Q Q A P A P -'∈⇔∈.证 因为Q A P ∈, 所以10n X R⨯∀≠∈, 有110()n Y Q X R -⨯'∀≠=∈且11111()()()0X Q A X X Q A X X A Q X X Q QA Q X Y QAY -----''''''''''====>, 所以、1T Q A P -∈. 反之, 若1T Q A P -∈, 则由必要性知:11()()QQ A A P P --''=∈=. 推论4 设n n n Q R ⨯∈, 则1Q E A P Q A P -'∈⇔∈. 推论5 设n n n Q R ⨯∈, 则1()Q E A P A Q P -'∈⇔∈. 推论6 设n nn Q R ⨯∈, 则Q E A P A Q P ''∈⇔∈. 定理13[13] 设n nn Q R ⨯∈, 则1T Q Q A P A P -∈⇔∈.证 因为Q A P ∈, 所以由定理10知:E QA P ∈, 又1()E QA P -∈, 故10n X R ⨯∀≠∈, 有10n Y QX R ⨯∀≠=∈, 且1111()()()0X Q A X QX A Q QX Y QA Y ----''''==>. 因而1T Q A P -∈. 反之, 若1T Q A P -∈, 则由必要性知:11()()Q Q A A P P --''=∈=.推论7 设n nn Q R ⨯∈, 则11Q E A P A Q P --∈⇔∈.推论8 设n nn Q R ⨯∈, 0,QQ A P >∈, 则A 的伴随矩阵*T Q A P ∈. 定理14[13] 设,n nn Q M R ⨯∈, 则1Q QM A P MA P -∈⇔∈ 证 由定理10知:()Q E E Q A P QA P M QAM Q M AM P M AM P '''∈⇔∈⇔=∈⇔∈.3.1.3 Schur 定理与华罗庚定理的推广定理15[13](广义Schur 乘积定理) 设12,D D 的n 阶可逆对角阵, 12,D D A P B P ∈∈, 且22B D D B '=, 则12D D A B P ∈ .证 因为12,D D A P B P ∈∈, 所以由定理11知:12,E D A D B P ∈, 又'2222()D B B D B D D B '''===, 故2nD B S +∈, 于是由定理7知:1212()()()E D D A B D A D B P =∈ , 故由定理11知:12D D A B P ∈ .定理15是Schur 定理的推广, 其中, 取12,D D E A A '===, 便得著名结论:推论9(Schur 乘积定理) 设,n A B S +∈, 则n A B S +∈ . 推论10 设D 为n 阶可逆对角阵, ,D n A P B S +∈∈, 则D A P ∈.推论11(广义华罗庚定理) 设D 为n 阶可逆对角阵, ()ij D A a P =∈, 且A D DA '=, 则()k kij D M a P =∈(其中k 为正整数).证 因为D A P ∈, 所以E DA P ∈, 且()DA A D A D DA ''''===, 所以n DA S +∈, 所以, 由华罗庚定理知:()()()()k k E D M D A A A DA DA DA P ==∈ , 故k D M P ∈.3.1.4 Ky Fan 等著名不等式的推广定理16[13] 设n nn Q R ⨯∈, 0,,,(),2QQ A B P QB QB t '>∈=为1B A -的非实特征值个数, 且1/()m n t ≥-, 则m m mA B A B +≥+, 特别当0t =时有广义Minkowski 不等式:1/1/1/nnnA BAB+≥+.证 因,Q A B P ∈, 故,E QA QB P∈, 又()QB QB '=, 故n QB S +∈, 再11()()QB QA B A --=的非实特征值个数为2,()1t m n t -≥, 所以由推论可知:mmmQA QB QA QB +≥+, 即1/1/1/nnnA BAB+≥+.定理17[13] 设n nn Q R ⨯∈, 0,,,(),2Q Q A B P QB QB t '>∈=为1B A -的非实特征值个数, 则[0,1]q ∀∈, 有1(1)()/2q qt qA q B A B -+-≥;特别地, 当0t =时, 有广义KyFan 不等式:1(1),([0,1])qqqA q B A Bq -+-≥∈.证 [0,1]q ∀∈, 由,Q A B P ∈, 有,(1),Q qA q B P -∈,11((1))()[/(1)]q QB qQA q q B A ---=-(即1B A -)的非实特征值个数为2t , 于是有:1/1/1/1/1/(1)/1(1)(1)((1))/2((1))/2 /2/2nnnnnn n t tq nq nnq nq nnnttqQA q QBqQAq QBq QAq QBQA QB Q AB--+-≥+-=+-≥=两边消去1nQ, 再n 次方得:1(1)()/2q qt qA q B A B-+-≥. 当10q q ==或时, 不等式显然成立.上述定理、推论、不等式等结论是在对广义正定矩阵的研究的基础上, 获得了许多新的结果, 改进并推广了著名的Schur 定理、华罗庚定理、Minkowski 不等式及Ky Fan 不等式, 并将各类正定矩阵与次正定阵统一了起来, 这对完善正定阵理论和应用很有价值.3.2 正定矩阵在实际问题中的应用正定矩阵在实际生活中有着广泛的应用, 以下给出几个应用举例. 3.2.1 二次型理论的应用二次型理论有着十分广泛的应用, 其在解决二次曲线与二次曲面方面、证明不等式等方面有着显著的实际应用, 下面就这几方面问题举例说明:一、从二次型理论的起源, 即从化二次曲线和二次曲面为标准型的问题入手, 我们发现二次型理论对二次曲线和二次曲面方程的化简有着重要意义.例1 利用直角坐标变换化简如下二次曲面的方程.032682223222=++--+++z y x xy z y x , 其中),1,3,4(),,,(--='='B z y x X⎪⎪⎪⎭⎫ ⎝⎛=200021013A .解:作平移变换:),,(,321ααααα='-=Y X , 则有03)(2)()(=+-'+-'-αααY B Y A Y , 即 0322=+'-'+'+'-'-'αααααB Y B A AY A Y AY Y令32+'-'=αααβB A , 又∵A A AY A Y =''=',αα, ∴0)(2=+'--'βαY B A AY Y适当选取α, 使A B α=, 由A =秩秩A=3知:A B α=(线性方程组)有唯一解:12311,2ααα===, 由,,A B α'可得92β=-.又∵A 是可逆实对称阵,∴存在正交阵T , 使得123T AT λλλ⎛⎫ ⎪'=⎪ ⎪⎝⎭使得12355552,,22λλλ+-===为A 的特征根, 作正交线性替换123,(,,)Y TZ Z Z Z Z ''''==,则2222221122331235555222Y AY Z Z Z Z Z Z λλλ+-'''''''=++=++ 即:原方程可化简为22212355552022Z Z Z +-'''++=二、在不等式的证明中, 恰当运用二次型理论, 将十分有助于问题的解决. 例2 求证:2211()nnii i i n X X ==≥∑∑证 令2211()n nii i i f X X X AX =='=-=∑∑, 则1212111(,,,)111111n n X n X f X X X n n X ⎛⎫---⎛⎫ ⎪ ⎪ ⎪=--- ⎪⎪ ⎪--- ⎪⎝⎭⎝⎭∴111111111n A n n ---⎛⎫⎪=--- ⎪⎪---⎝⎭, ∵A 的顺序主子式大于或等于零 ∴A 是半正定的, ∴二次型f 是半正定的, 即0f ≥ 即2211()nnii i i n X X ==≥∑∑.3.2.2 仿射变换仿射变换是从运动变换到射影变换的桥梁.通过探究和证明我们知道,通过平行射影不改变的性质和数量,称为仿射不变性质和仿射不变量,经过仿射其对应的性质也是不变的. 在初等几何问题中,圆和椭圆都是比较常见的图形,圆比椭圆更特殊,它有很多很好的性质,与圆有关的定理举不胜举,但椭圆则不然.因其本身的定义要比圆复杂,椭圆的性质和定理就很少,解决一个与椭圆有关的问题要比解决一个与圆有关的相应的问题困难得多.因此,我们自然期望有一种方法,使得处理有关椭圆的问题和处理有关圆问题一样容易由仿射变换性质可知:椭圆通过适当的仿射变换可变成圆.所以,只要考虑有关椭圆仿射性质的问题,就可以先转化为有关圆的相应的问题来解决,再把所得的结果推广到椭圆中去, 即可达到我们解题的目的.下面仅就00x a x y b y '⎛⎫⎛⎫⎛⎫= ⎪ ⎪⎪'⎝⎭⎝⎭⎝⎭错误!未找到引用源。
矩阵初等变换及其应用毕业论文
矩阵初等变换及其应用毕业论文矩阵初等变换及其应用毕业论文摘 要:初等变换是高等代数和线性代数学习过程中非常重要的,使用非常广泛的一种工具。
本文列举了矩阵初等变换的几种应用,包括求矩阵的秩、判断矩阵是否可逆及求逆矩阵、判断线性方程组解的状况、求解线性方程组的一般解及基础解系、证向量的线性相关性及求向量的极大无关组、求向量空间两个基的过渡矩阵、化二次型为标准形。
并用具体例子说明矩阵初等变换在以上几种应用中是如何运用的。
关键词:矩阵 初等变换 初等矩阵在代数的学习过程中,我发现矩阵的初等变换有许多应用,几乎贯穿着始终。
本文将对矩阵的初等变换进行介绍并以具体例子说明矩阵初等变换的七种应用。
虽然这些计算格式有不少类似之处,但是也指出由于这些计算格式有不同的原理,所以它们的应用也有一些明显的区别。
定义1:矩阵的行(列)初等变换是指对一个矩阵施行的下列变换: (1)交换矩阵的两行(列)(交换第i ,j 两行(列),记作()ij ij r c );(2)用一个不等于零的数乘矩阵的某一行(列)即用一个不等于零的数乘矩阵的某一行(列)的每一个元素(用数k 乘以第i 行(列),记作()(())i i r k c k ;(3)用某一个数乘矩阵的某一行(列)后加到另一行(列),即用某一数乘矩阵的某一行(列)的每一个元素再加到另一行(列)的对应元素上(第i 行(列)k 倍加到第j 行(列),记作()(())ij ij r k c k 。
初等行、列变换统称为初等变换。
定义2:对单位矩阵I 仅施以一次初等变换后得到的矩阵称为相应的初等矩阵,分别记为第1、2、3类行(列)初等矩阵为()ij ij R C ,()(())i i R k C k ,()(())ij ij R k C k ,有ij R =ij C =10111⎛⎫ ⎪ ⎪ ⎪⎪⎪ ⎪ ⎪ ⎪ ⎪⎝⎭()i R k =()i C k =1k1⎛⎫ ⎪⎪⎪ ⎪ ⎪ ⎪⎝⎭()ij R k =()ij C k =11j 11i k⎛⎫⎪ ⎪ ⎪⎪⎪ ⎪ ⎪ ⎪ ⎪⎝⎭行行 初等变换与初等矩阵之间有下列基本性质。
(完整版)matlab毕业课程设计(关于矩阵操作)
课程设计任务书学生姓名:李圆超专业班级:信息SY1201指导教师:龙毅宏工作单位:信息工程学院题目:MATLAB矩阵操作设计要求完成的任务:1.利用MATLAB-mathematics对矩阵操作进行设计,具体包括创建(普通、单位、零)矩阵、矩阵运算、矩阵变形等。
2.验证如下函数的功能:all、any、find、isempty、isequal、xor。
课程设计的目的:1.理论目的理解掌握所学理论知识,并能用所学理论知识分析矩阵的操作,设计。
2.实践目的熟练MATLAB的使用,验证并掌握MATLAB的一些基本函数,对矩阵进行实际的操作和设计。
时间安排:指导教师签名:年月日系主任(或责任教师)签字:年月日目录摘要............................................................ . (2)Abstract.................................................... .. (3)1引言............................................................ . (4)1.1M A T L A B的介绍 (4)1.2矩阵的介绍 (4)2M A T L A B对矩阵的操作 (4)2.1矩阵的生成 (4)2.1.1直接输入法 (4)2.1.2特殊矩阵的生成 (6)2.2矩阵的运算 (8)2.2.1矩阵的加减 (8)2.2.2矩阵的相乘 (9)2.2.3矩阵的乘方 (9)2.2.4矩阵的除法 (11)2.2.5矩阵的点乘 (11)2.2.6矩阵的数乘 (12)2.2.7矩阵的转置 (12)2.2.8矩阵的逆 (13)2.2.9求矩阵的特征值和特征向量 (13)2.3矩阵的变形 (14)3验证部分函数的功能 (15)3.1a l l和a n y函数 (15)3.2f i n d函数............................................................173.3i s e m p t y函数 (20)3.4i s e q u a l函数 (20)3.5x o r函数............................................................ .2 1 4心得体会............................................................ .....235参考文献............................................................ .....2 4摘要矩阵是高等代数学中的一种运算工具。
矩阵函数的原理与应用
矩阵函数的原理与应用1. 矩阵函数的基本概念矩阵函数是指将一个矩阵作为输入,并输出另一个矩阵的函数。
矩阵函数的输入和输出可以是任意维数的矩阵,且可以进行各种运算。
矩阵函数的原理主要基于线性代数的理论。
2. 矩阵函数的推导与定义矩阵函数的推导过程涉及到矩阵的特征值和特征向量的计算。
通过分析矩阵的特征值和特征向量,可以得到矩阵函数的定义和性质。
常见的矩阵函数包括指数函数、对数函数、幂函数等。
3. 矩阵函数的应用领域矩阵函数在科学研究和工程应用中具有广泛的应用。
以下列举几个典型的应用领域:•线性方程组求解:矩阵函数可以通过求解线性方程组来解决实际问题,如物理模拟、数据拟合等。
•信号处理:矩阵函数可以用于处理信号,如图像处理、音频处理等。
•优化问题:矩阵函数可以用于求解优化问题,如最小二乘法、最大似然估计等。
•自动控制:矩阵函数可以用于设计和分析控制系统,如PID控制、模糊控制等。
•机器学习:矩阵函数在机器学习算法中有着重要的应用,如主成分分析、支持向量机等。
4. 矩阵函数的算法与实现矩阵函数的求解算法有多种,常见的有幂法、矩阵对角化等。
矩阵函数的实现可以通过各种编程语言和数值计算库来完成,如MATLAB、Python的NumPy库和SciPy库等。
5. 矩阵函数的性质与扩展矩阵函数具有一些基本性质,如可逆性、对角化等。
此外,还存在一些特殊的矩阵函数,如矩阵的广义逆、矩阵的广义特征值等。
6. 总结矩阵函数作为线性代数的一个重要分支,在科学研究和工程应用中具有广泛的应用。
通过对矩阵函数的理解和应用,可以更好地解决实际问题,并推动科学技术的发展。
以上是矩阵函数的原理与应用的简要介绍,希望对读者有所帮助。
深入学习和掌握矩阵函数的原理和应用,将有助于扩展自己的专业知识和提升解决实际问题的能力。
矩阵函数以及应用毕业设计
天津科技大学2014届本科生毕业设计
的历史长河中矩阵理论的创立者被一致认为是英国数学家凯莱(Cayley),是他最先将矩阵 作为一个单独的数学上的概念提出来,并且关于矩阵的很多学术论文和著作都是他最早发 表的。事实上最早的矩阵是从对大量行列式的研究中分离出来的,因为和行列式对应的方 阵本身就可以做许多的研究和运用,随着对行列式研究的深入,矩阵的许多知识点也日渐 完善。从逻辑上讲,概念应先于行列式的矩阵的概念和历史上真正的顺序是恰恰相反的。 在19世纪50年代,英国数学家凯莱(Cayley)公开展示了自己关于矩阵的最新研究成果-《矩阵论的研究报告》,这项研究成果使我们对矩阵的认识更深入了一步。本文定义了矩阵 相等、矩阵的算法、矩阵的转置和基本概念,如矩阵的逆矩阵的加法,给出了系列,互换性 和约束力。除此之外,英国数学家凯莱(Cayley)也给出了方阵的特征根(特征值),还有其他 许多结论。矩阵的发展历史,著名的德国数学家弗洛伯纽斯(Frobenius)起着非常重要的作 用,他是第一个对矩阵中最小多项式问题作全面介绍的著名数学家。他还介绍了矩阵的秩 、不变的因素和主要因素、正交矩阵相似变换等知识,矩阵的其他概念如合同,不变的因素 和主要因素理论的逻辑排列的形式等等在他的著作中也有体现。在19世纪50年代,约丹经 过潜心研究首先发表了把一般矩阵化为标准型矩阵的方法。到了19世纪90年代,梅茨勒(M etzler)首先提出了矩阵函数的基本概念,最后找到用幂级数形式将表示矩阵的方法,这些 对矩阵的发展意义重大。此外,傅立叶(Fourier)与庞加莱(Poincare)研究的主要是无穷矩阵 方面。到这时,矩阵已经相当完善了。 矩阵最大的用途就是在实践中解用常规方法难以求解的方程。另外一个在实际操作中 很有意义的作用是代表线性变换,即是像f(x)、4x之类的关于线性函数的推论。矩阵的特征 向量可以揭示一个线性变换的深层次特征。随着两个世纪中无数数学家的无私奉献,矩阵 论已经成为了一门完善的数学分支。矩阵在很多方面都有重要应用,例如数学领域里,力 学、物理学、工程数学、经济管理方面都有矩阵的出现。
《矩阵函数的应用》教案
《矩阵函数的应用》教案教案:矩阵函数的应用一、教学目标1.理解矩阵函数的概念,了解矩阵函数的代数运算规则;2.能够运用矩阵函数解决实际问题;3.培养学生分析问题、运用矩阵函数解决问题的能力。
二、教学内容1.矩阵函数的概念;2.矩阵函数的代数运算规则;3.矩阵函数的应用实例。
三、教学过程1.导入新课教师可用一道题目导入新课:“如何计算矩阵的平方根?”引导学生思考,探讨矩阵的开平方运算是否可行,如果可行该如何计算。
2.提出问题教师通过提出问题引发学生的兴趣:“矩阵函数在实际生活中有哪些应用?”让学生先讨论并列举出一些例子,然后引导学生思考如何利用矩阵函数解决这些实际问题。
3.知识讲解(1)矩阵函数的概念教师向学生讲解矩阵函数的概念:“矩阵函数是指将一个矩阵映射为另一个矩阵的函数。
在矩阵函数中,每一个矩阵元素可以看作是矩阵函数的一个可能值。
”教师通过示例详细阐述矩阵函数的概念。
(2)矩阵函数的代数运算规则教师介绍矩阵函数的代数运算规则:“对于任意一个矩阵函数f(X),满足以下代数运算规则:1) f(X)+f(Y)=f(X+Y),2) f(aX)=af(X),其中a为任意标量。
”4.实例分析教师以实例的形式分析矩阵函数的应用:(1)矩阵函数在图像处理中的应用教师通过图像处理的实例,教学矩阵函数在图像处理中的应用。
“在图像处理中,我们常常需要对图像进行模糊处理,以达到其中一种效果。
通过构造合适的矩阵函数,我们可以对图像进行模糊处理。
”(2)矩阵函数在金融风险控制中的应用教师通过金融风险控制的实例,教学矩阵函数在金融风险控制中的应用。
“在金融风险控制中,我们常常需要对风险指标进行计算,通过构造合适的矩阵函数,我们可以对风险指标进行计算,并做出相应的风险控制策略。
”5.解决问题教师组织学生进行小组讨论,选取一个具体的实际问题,要求学生利用矩阵函数解决该问题。
学生通过小组合作,利用矩阵函数解决问题,并向全班进行汇报。
毕业设计(论文)-用c语言实现矩阵的运算[管理资料]
职称
指导教师评语及对成绩的评定意见
该同学的毕业设计选题来源于对对矩阵运算算法的研究中。该同学按照毕业设计任务书的要求开展了用C语言实现矩阵运算的研究,该研究对矩阵运算的研究有重要意义。
该同学能够按照毕业设计工作计划,在广泛阅读有矩阵运算参考文献的基础上,基于C语言程序,设计了矩阵几种运算的算法;利用C语言编程的特点和矩阵运算的特点开展了算法设计的研究。论文不足之处在于没有通过计算机运行,应在这一方面作进一步的研究。
B=
为下三角形矩阵。
如果方阵的主对角线以外的元全为零,即
A=
则称它为对角矩阵,记作 或 。在对角矩阵 中,未写出的元表示零元。
主对角线上全为1的n阶对角矩阵,即
称为n阶单位阵,记作 或 。
矩阵的应用非常广泛,可以说我们日常生产生活中都会应用到矩阵。
在白酒工业,成品酒的勾兑这一工序尤为重要,勾兑的目的不仅要使成品酒达到规定的酒精度,更重要的是要使成品酒中影响酒体风味的几十种主要微量成分达到预先设计好的平衡比例和具体含量。因此,使用计算机配合相关勾兑软件,就显得非常重要。如基酒及调配液的配比计算模块中,按目标含量、实际含量、调配液密度、成品密度、调配比数N类参数,依照质量守恒原理,建立起质量平衡线性方程组;且编制有行列式计算模块,可对任意阶行列式进行计算,从而改变了传统的逐次迭代算法,可直接对线性方程组进行解的存在性的判定和求解计算;小到针对单一的酒精度勾兑的配比计算,大到同时针对30种目标成分、31种基酒及调配液的调配比进行计算。
我国高速公路网络中环形结构日益增多,发展形成了复杂的网状结构,如何实现准确的收费清分就成为不得不解决的问题。由于高速公路投资费用极高 ,在实际路网中相邻节点间有两条或以上的路段直接相连的情况一般不会出现,故高速路网的结构图一般都是简单图。因此,可以定义若干矩阵,实现对高速公路网(主线、 匝道、 交叉口等) 的参数化描述。如可通过连接矩阵、属性矩阵和标识站矩阵描述高速公路的物理结构,通过规则矩阵描述收费清分规则。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
矩阵函数以及应用毕业设计1 绪论1.1 矩阵(Matrix)的发展与历史人们对矩阵(Matrix)的研究历史非常悠久,在很久以前就已经有人研究过了幻方和拉丁方阵。
在过去的很长时间内,矩阵都是人们解决线性问题的最主要方法。
成书于汉朝前期的《九章算术》,在表示线性方程组的过程中使用了将方程中不同系数分开的方法,这种方法在后来的不断演化下最终得到方程的增广矩阵。
在计算的过程中经常使用矩阵的初等变换进行消元,具体说就是通过一些计算技巧将前面给出的增广矩阵化为行最简型。
但是当时我们能知道的矩阵知识非常的少,虽然过去的标准和现在的矩阵在表示上已经非常的类似了,但这两者都是以线性方程为基本标准。
事实上子宫基质的控制中心和开始生活意义的地方是矩阵最开始的意义,所以说矩阵有生命的意义。
在数学中,开始出现的是对现在数学都有决定性的行列式,但需要行列式的行和列相等,最终的排成的表都是方的,随着研究的深入人们发现行数等于列数的行列式已经无法满足现实生活中的实际需要了。
在这种情况下,矩阵应运而生。
现在对于我们来说非常熟悉的矩阵和行列式,它们的概念是非常的不一样的。
行列式能按照我们的规则计算出它的结果,而矩阵是将数字按一定顺序排列得到的。
在学术研究中恰当地使用矩阵,能用向量空间中的向量表示线性方程组中系数矩阵;因此,一个多元线性方程组的解的情况,以及一系列问题的理论解之间的不同关系,都可以得到彻底解决。
矩阵都有自身的行和列,水平的称之为行,竖直的称之为列。
这些我们现在能看到的关于矩阵的一切都是由无数数学家的摸索得来的。
矩阵(Matrix)在数学发展历史上有着非常重要的位置,它一直是数学研究的一个主要方面,是数学在研究和应用过程中经常用到的知识。
“矩阵”由英国数学家叶(Sylvester)第一次使用,他使用的这个数学术语最后将矩阵的列数和早期的行列式分离开来。
在数学发展的历史长河中矩阵理论的创立者被一致认为是英国数学家凯莱(Cayley),是他最先将矩阵作为一个单独的数学上的概念提出来,并且关于矩阵的很多学术论文和著作都是他最早发表的。
事实上最早的矩阵是从对大量行列式的研究中分离出来的,因为和行列式对应的方阵本身就可以做许多的研究和运用,随着对行列式研究的深入,矩阵的许多知识点也日渐完善。
从逻辑上讲,概念应先于行列式的矩阵的概念和历史上真正的顺序是恰恰相反的。
在19世纪50年代,英国数学家凯莱(Cayley)公开展示了自己关于矩阵的最新研究成果--《矩阵论的研究报告》,这项研究成果使我们对矩阵的认识更深入了一步。
本文定义了矩阵相等、矩阵的算法、矩阵的转置和基本概念,如矩阵的逆矩阵的加法,给出了系列,互换性和约束力。
除此之外,英国数学家凯莱(Cayley)也给出了方阵的特征根(特征值),还有其他许多结论。
矩阵的发展历史,著名的德国数学家弗洛伯纽斯(Frobenius)起着非常重要的作用,他是第一个对矩阵中最小多项式问题作全面介绍的著名数学家。
他还介绍了矩阵的秩、不变的因素和主要因素、正交矩阵相似变换等知识,矩阵的其他概念如合同,不变的因素和主要因素理论的逻辑排列的形式等等在他的著作中也有体现。
在19世纪50年代,约丹经过潜心研究首先发表了把一般矩阵化为标准型矩阵的方法。
到了19世纪90年代,梅茨勒(Metzler)首先提出了矩阵函数的基本概念,最后找到用幂级数形式将表示矩阵的方法,这些对矩阵的发展意义重大。
此外,傅立叶(Fourier)与庞加莱(Poincare)研究的主要是无穷矩阵方面。
到这时,矩阵已经相当完善了。
矩阵最大的用途就是在实践中解用常规方法难以求解的方程。
另外一个在实际操作中很有意义的作用是代表线性变换,即是像f(x)、4x之类的关于线性函数的推论。
矩阵的特征向量可以揭示一个线性变换的深层次特征。
随着两个世纪中无数数学家的无私奉献,矩阵论已经成为了一门完善的数学分支。
矩阵在很多方面都有重要应用,例如数学领域里,力学、物理学、工程数学、经济管理方面都有矩阵的出现。
1.2 本文所做的主要工作矩阵理论包含的内容非常非常多,矩阵函数在矩阵理论中占据非常重要的位置,相比于矩阵函数中的其他知识,矩阵多项式比较容易理解,就是这样容易理解的矩阵多项式是我们对矩阵函数进行研究的理论基础。
矩阵函数的定义方式有多种,本文主要是从多项式和幂级数两个方面进行研究的。
本文主要论述了矩阵函数以及应用。
在文章的第一部分,总结了矩阵函数所必须的基础知识,主要包括代数学多项式理论、行列式与矩阵等方面的一些结论以及数学分析中幂级数的若干法则。
文章的第二部分,总结了矩阵函数的概念、性质、推论,介绍了若干重要的矩阵函数。
文章的第三部分,归纳了矩阵函数的若干计算方法,包括了Hamiltio-Cayley定理、利用相似对角化计算、利用Jordan标准型法进行计算、利用待定系数法求解等四种计算方法。
在这部分的最后对这四种方法进行了比较,在比较中加深对矩阵函数求解的认识。
可以根据计算过程中遇到的实际情形加以选择,将会给计算带来很大方便。
本文的第四部分,通过查阅文献和指导教师交流的方式,在求解线性微分方程过程中有对矩阵函数的应用研究,并介绍了在线性系统的可控性和可观性中矩阵函数的应用。
本文的最后部分,通过Matlab编写能计算常用矩阵函数的程序,将使矩阵函数的计算更方便、迅速。
2 矩阵函数2.1 研究本论文具备的数学基础为了进一步讨论和便于理解,引入以下研究本论文的相关概念:1、线性空间在集合上具有一定的结构或符合一定的要求,那么这个集合就是特定的空间。
如果是非空的集合,是数域。
对里的元素定义代数类运算,叫作加法;就是给出一种规则,使中任意两个元素和,都能在中找到唯一的一个和它匹配,其中是与的和,记为。
在数域与集合中的元素再定义另外一种运算,叫作数量乘法;就是如果数域中任何一数与中的任何一个元素,在中都能找到一个元素和它匹配,是和的数量乘积,记为。
若加法与数乘都同时符合它们的运算法则,那么就叫作数域上的线性空间。
2、级数级数知识是分析科学中一个重要的部分;这个概念经常出现在数学的其他分支。
把数列的项,,…,,…逐项相加得到的函数。
数项级数简称级数。
如:,缩写为,就是级数的通项,记作是级数的部分和。
如果当时,数列极限有,级数就是收敛的,否则就是发散的。
研究函数经常会用到级数,它不管在理论上还是实际中都有很多用途,原因主要有一下两个方面:一、许多经常用到的非初等函数可以用级数表示,级数还可以表示微分方程的解;二、函数可以用来表示级数,也能用级数去探讨函数的性质。
幂级数,是级数中非常重要的一种,被当作基础知识应用在实变型函数、复变型函数和其他许多基本领域中,在这些领域发挥巨大的作用。
幂级数是指每一项均对应着级数项序号的常数倍的的次方(是从0递增的自然数,是常数)。
幂级数与多项式形式非常接近,在许多方面有相似的特征,可以被视为“无限的多项式”。
3、正定矩阵在线性代数里,正定矩阵有时会简称为正定阵。
它的定义有广义和狭义之分。
广义定义:设是阶方阵,如果有任意非零向量,都有,是的转置,称为正定矩阵。
例如:为阶矩阵,为单位矩阵,为正实数。
在充分大时,为正定矩阵。
(必须为对称阵)狭义定义:一个阶的实对称矩阵是正定的当且仅当对于所有的非零实系数向量,都有。
其中表示的转置。
4、线性算子线性算子,有数学运算各领域的线性性质(如线性变换,线性代数理论的微分方程,积分方程理论,微分,积分,积分变换)的抽象概括。
它是研究线性泛函的一个重要目标。
线性算子的用途很广,不但应用在数学的很多分支当中,同时对于量子物理也是重要的数学基础。
5、对称矩阵和反对称矩阵对称矩阵的定义是:(的转置),对称的矩阵元素。
反对称矩阵的定义为:(的转置前加负)它的首行与首列各元素绝对值相等,符号相反。
即, 因此,在对角线上的元素,,有, 在非偶数域中,有,即反对称矩阵对角线元素为零,此性质只在非偶数域中成立。
6、化零(零化)多项式给定矩阵,如果多项式,满足,则称是的化零多项式,(一般取首项系数为1)。
7、矩阵的谱半径 设A 是n n ⨯矩阵,是其特征值, = 1,2,……,。
下面通过数学式子将其表示出来。
假如()A ρ表示A 的谱半径,即(){}max A A ρλλ= 是的特征值。
也就是说矩阵A 的谱半径是矩阵A 的全部特征值求模的最大值;如果特征值是虚数,谱半径就是实部和虚部的平方和求算术平方根。
8、 n n F ⨯表示数域F 上n n ⨯矩阵全体的线性空间;9、 n n C ⨯表示n n ⨯复矩阵集;10、 (),P F λ数域F 上λ的纯量多项式;11、矩阵的谱 矩阵通过数学运算计算出来的特征值的集合就是一个矩阵的谱,通过数学表达式表示出来也就是:()A σ表示A 的谱,即(){}A A σλλ=是的特征值;12、其中次数最低的零化多项式称为矩阵A 的最小多项式,记做()m ψλ;13、文献[1]给出矩阵级数的定义:定义1:设{}k A 是m n C ⨯的矩阵序列,其中()k k m n ij A a C ⨯=∈,无穷和123k A A A A ++++L L 称为矩阵级数,记为1k k A ∞=∑.对正整数1k ≥,记1kk i i S A ==∑称()k S 为矩阵级数1k k A ∞=∑的部分和,如果矩阵序列{}k S 收敛,且有极限S ,即lim k k S S →∞=,则称矩阵级数1k k A ∞=∑收敛,并称S 为矩阵级数1kk A ∞=∑的和,记为1k k A S ∞==∑。
不收敛的矩阵级数称为发散的.定义2:设n n A C ⨯∈,形如20120k k k k k c A c I c A c A c A ∞==+++++∑L L的矩阵级数称为矩阵幂级数.14、相似矩阵 设是阶矩阵,如果存在阶可逆矩阵使,则称矩阵与相似的,记为.相似矩阵代表等价的关系。
15、可对角化矩阵 如果阶方阵能与一个对角矩阵相似,就说可对角化。
阶方阵可对角化的充要条件是它有个线性无关的特征向量。
对角矩阵(diagonal matrix)是一个矩阵主对角线之外的所有元素都是0。
对角线上的元素可以是0或任何其他值。
然后引入线性无关的概念。
对向量组,如果有一组不全为零的数,然后 被称为向量组线性相关.如果没有这样的,换句话就是向量等式当且仅当才成立,就称向量组是线性无关的.16、可逆矩阵 可逆矩阵是线性代数中的一种矩阵,其定义为在线性代数中,给定一个阶方阵,若存在一阶方阵,使得(或、 满足任意一个),其中 为阶单位矩阵,则称 是可逆的,且是的逆矩阵,记作。
2.2 矩阵函数的定义类比于代数中函数的定义,能知道定义域和值域都属于方阵的函数称为矩阵函数。
矩阵函数的定义方式有很多种,为了便于进一步的研究,本文主要从经常使用的多项式和幂级数来定义矩阵函数。