基于塞尔维斯特准则的二次型标量函数的正定性判断
- 格式:doc
- 大小:118.15 KB
- 文档页数:8
李雅普诺夫稳定性方法李雅普诺夫第一方法又称间接法,它是通过系统状态方程的解来判断系统的稳定性。
如果其解随时间而收敛,则系统稳定;如果其解随时间而发散,则系统不稳定。
李雅普诺夫第二方法又称直接法,它不通过系统状态方程的解来判断系统的稳定性,而是借助李雅普诺夫函数对稳定性作出判断,是从广义能量的观点进行稳定性分析的。
例如有阻尼的振动系统能量连续减小(总能量对时间的导数是负定的),系统会逐渐停止在平衡状态,系统是稳定的。
由于李雅普诺夫第一方法求解通常很烦琐,因此李雅普诺夫第二方法获得更广泛的应用。
李雅普诺夫第二方法的难点在于寻找李雅普诺夫函数。
迄今为止,尚没有通用于一切系统的构造李雅普诺夫函数的方法。
对于系统[]t ,f x x= ,平衡状态为,0e =x 满足()0f e =x 。
如果存在一个标量函数()x V ,它满足()x V 对所有x 都具有连续的一阶偏导数;同时满足()x V 是正定的;则 (1)若()x V 沿状态轨迹方向计算的时间导数()dt /)(dV Vx x = 为半负定,则平衡状态稳定;(2) 若()x V 为负定,或虽然()x V 为半负定,但对任意初始状态不恒为零,则平衡状态渐近稳定。
进而当∞→∞→)(V x x 时,,则系统大范围渐近稳定;(3) 若()x V为正定,则平衡状态不稳定。
判断二次型x x x P )(V τ=的正定性可由赛尔维斯特(Sylvester )准则来确定,即正定(记作V(x)>0)的充要条件为P 的所有主子行列式为正。
如果P 的所有主子行列式为非负,为正半定(记作V(x)≥0);如果-V(x)为正定,则V(x)为负定(记作V(x)<0);如果-V(x)为正半定,则V(x)为负半定(记作V(x)≤0)。
例:[]正定。
则)(V 01121412110,041110,010x x x 1121412110x x x )(V 321321x x >---->>----=⎥⎥⎦⎤⎢⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡ 例:)x x (x x x )x x (x x x 22212122221121+--=+-=(0,0)是唯一的平衡状态。
第四章系统稳定性及其李雅普诺夫稳定4-1 稳定性一般概念对于一个实际的控制系统,其工作的稳定性无疑是一个极其重要的问题,因为一个不稳定的系统在实际应用中是很难有效地发挥作用的。
从直观上看,系统的稳定性就是一个处于稳态的系统,在某一干扰信号的作用下,其状态偏离了原有平衡位置,如果该系统是稳定的,那么当干扰取消后有限的时间内,系统会在自身作用下回到平衡状态;反之若系统不稳定,则系统永远不会回到原来的平衡位置。
系统的稳定一般有外部稳定和内部稳定两种。
外部稳定又称作输出稳定,也就是当系统在干扰取消后,在一定时间内,其输出会恢复到原来的稳态输出。
输出稳定有时描述为系统的BIBO稳定,即有限的系统输入只能产生有限的系统输出。
系统内部稳定主要针对系统内部状态,反映的是系统内部状态受干扰信号的影响。
当扰动信号取消后,系统的内部状态会在一定时间内恢复到原来的平衡状态,则称系统状态稳定。
在经典控制论中,研究对象都是用高阶微分方程或传递函数描述的单输入单输出(SISO)系统,反映的仅是输入输出的关系,不会涉及系统内部的状态。
因此经典控制论中只讨论系统的输出稳定问题。
系统的稳定性是系统本身的特性,与系统的外部输入(控制)无关。
在经典控制论中,我们通过研究线性定常系统的特征根的情况来判断系统的输出稳定性:如果系统的特征根都有负的实部(即都在复平面的左部),则系统输出稳定。
对于n阶线性连续系统,其特征方程为:…………………………(4-1)当n≥4时,要求出其所有特征根是非常困难的,从而要想通过解出高阶系统的特征根来判别系统稳定性也是不现实的。
所以1877年劳斯(Routh)和1895年霍尔维茨(Hurwitz)分别提出了有名的劳斯-霍尔维茨稳定判据,它可以通过线性定常系统特征方程的系数的简单代数运算来判别系统输出稳定性,而不必求出各个特征根。
有关Routh-Hurwitz判据的详细内容请参阅有关经典控制论教材。
当系统不是线性定常系统时,或者对于系统内部状态稳定问题,经典控制论中的方法就不好解决了,这就需要下面介绍的李雅普诺夫(Lyapunov)稳定性的理论。
二次型与正定性二次型是高等数学中的一个重要概念,正定性则是与二次型紧密相关的性质。
本文将介绍二次型及其性质,深入探讨正定性的定义、判别方法以及与正定矩阵的关系。
一、二次型的定义二次型是指形如\[Q(x_1,x_2,\dots,x_n) = \sum_{i=1}^n \sum_{j=1}^n a_{ij}x_ix_j\]的函数,其中\(a_{ij}\)为实数或复数,称为二次型的系数。
\(x_1,x_2,\dots,x_n\)为实数或复数,称为二次型的变量。
二次型可以用矩阵的语言来表示,即\[Q(\mathbf{x}) = \mathbf{x}^T A \mathbf{x}\]其中\(\mathbf{x} = \begin{bmatrix} x_1 \\ x_2 \\ \dots \\ x_n\end{bmatrix}\)为列向量,\(A\)为二次型的系数矩阵,其元素为\(a_{ij}\)。
二、正定性的定义对于任意非零向量 \(\mathbf{x} = \begin{bmatrix} x_1 \\ x_2 \\ \dots \\ x_n \end{bmatrix}\),如果对应的二次型 \(Q(\mathbf{x})\) 满足条件:1. 当 \(\mathbf{x} \neq \mathbf{0}\) 时, \(Q(\mathbf{x}) > 0\);2. 当且仅当 \(\mathbf{x} = \mathbf{0}\) 时, \(Q(\mathbf{x}) = 0\)。
则称二次型 \(Q(\mathbf{x})\) 是正定的。
三、正定性的判别方法判断一个二次型是否正定存在多种方法,下面介绍两种常见的方法:特征值判别法和合同变换法。
1. 特征值判别法设 \(A\) 为二次型的系数矩阵,将 \(A\) 进行对角化得到对角矩阵\(D\),同时得到可逆矩阵 \(P\),使得 \(A = PDP^{-1}\)。
河北师范大学本科生毕业论文(设计)文献综述河北师范大学本科生毕业论文(设计)翻译文章本科生毕业论文设计正定二次型的性质与应用作者姓名:指导教师:所在学院:数学与信息科学学院专业(系):数学与应用数学班级(届):2013届数学B班二〇一三年四月二十八日目录中文摘要、关键字 (2)1 正定二次型与正定矩阵的概念 (3)1.1 二次型的概念 (3)1.2 二次型的矩阵形式 (3)1.3 正定二次型与正定矩阵的概念 (3)2 实正定矩阵的判定方法及证明 (4)2.1 利用定义判定 (4)2.2 利用标准型判定 (4)2.3 利用主子式判定 (8)2.4 其他常用判定 (11)3 实正定矩阵的应用 (15)3.1 用正定矩阵的定义来证明一些结论 (15)3.2 正定矩阵在数学分析上的应用 (17)3.2.1 多元函数的极值问题 (17)3.2.2 正定矩阵在积分中的应用 (19)3.3 正定矩阵在运筹中的应用 (19)3.3.1 具有约束方程的最优化问题 (19)3.4 用正定矩阵来证明不等式 (20)3.5 正定矩阵在几何中的应用 (21)3.5.1二次曲面的标准型 (21)参考文献 (23)英文摘要、关键字 (24)正定二次型的性质及应用数学与信息科学学院数学与应用数学专业指导教师高锁刚作者王敬摘要:本文以矩阵和向量为工具,研究了一种特殊的函数,即二次型。
然而在它的实际应用中许多二次型都是实二次型,其中最重要的一类是正定二次型。
本文主要阐述的是实矩阵的正定性以及应用,文中给出了实对称正定矩阵的多个判定定理和重要结论,从而使人们能够更好地使用正定矩阵这个工具。
全文共分三章,第一章主要叙述二次型及正定二次型、正定矩阵的定义;第二章主要列举说明正定性矩阵的几个判别方法;第三章简单地罗列一些实例来阐述实矩阵正定性的应用。
关键字:正定矩阵正定二次型特征值实对称矩阵1 正定二次型与正定矩阵的概念1.1[1] 二次型的概念设P 是一个数域,ij a ∈P, n 个文字1x ,2x ,…,n x 的二次齐次多项式()n n n x x a x x a x x a x a x x x f 11311321122111212...22,...,,++++=n n x x a x x a x a 22322322222...2++++......+2n nn x a +=∑∑==n i nj jiij xx a 11()n j i a a ji ij ,...2,1,,==称为数域P 上的一个n 元二次型,简称二次型.当ij a 为实数时, f 称为实二次型.当ij a 为复数时,称f 为复二次型.如果二次型中只含有文字的平方项,即),,,(21n x x x f =2221112...n n d x d x d x +++则称f 为标准型. 1.2 二次型的矩阵形式二次型),,,(21n x x x f 可唯一表示成),,,(21n x x x f =T x Ax ,其中12(,,...,)T n x x x x =,()ij n n A a ⨯=为对称矩阵,称上式为二次型的矩阵形式,称A 为二次型的矩阵(A 必是对称矩阵),称A 的秩为二次型f 的秩.1.3 正定二次型与正定矩阵的概念设),,,(21n x x x f =Tx Ax 是n 元实二次型(A 为实对称矩阵),如果对任意不全为零的实数12,,...,n c c c 都有12(,,...)0n f c c c >,则称f 为正定二次型,称A 为正定矩阵;如果12(,,...)0n f c c c ≥,则称f 为半正定二次型,称A 为半正定矩阵;如果12(,,...)0n f c c c <,则称f 为负定二次型,称A 为负定矩阵;如果0),,,(21≤n c c c f ,称f 为半负定二次型,称A 为半负定矩阵;既不是正定又不是负定的实二次型称为不定的二次型,称A 为不定矩阵.2 实正定矩阵的判定方法及证明2.1 利用定义判定定理1 实对称矩阵A ∈n n R ⨯是正定矩阵的充分而且必要条件是对于任意的n 维非零列向量x , 即n R x ∈≠0,使0>Ax x T .定理2[2] 实对称矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡n d d 1是正定矩阵的充分而且必要条件是0>i d , n i ,2,1=.证明:实对称矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡n d d 1是正定的充要条件是对任意的n 维非零列向量x , 即n R x ∈≠0,有T x ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡n d d 10>x , 令T x )0,,0,1( =,则得01>d ,同理,分别令x 为所有的单位列向量,则可得0>i d ,n i ,,2,1 =,所以定理可证.定理3 实对称矩阵n n R A ⨯∈是正定矩阵的充分而且必要条件是对任意的n R x ∈≠0,使二次型Ax x T 的秩和符号差均等于n .证明:因为实对称矩阵A 是正定矩阵,所以存在二次型Ax x T 为正定二次型,其规范形为22221n y y y +++ ,所以正惯性指数为n ,即得二次型Ax x T 的秩和符号差均等于n .所以A 是正定矩阵.2.2 利用标准型判定定理 4 [2] 实对称矩阵n n R A ⨯∈是正定矩阵的充分而且必要条件是A 与单位矩阵E合同,即存在实非奇异矩阵C ,使E AC C T =.证明:必要性,因为实对称矩阵A 是正定矩阵,所以矩阵A 对应的二次型Ax x T为正定二次型,可经过一适当的非退化线性替换TY X =化为规范形22221ny y y +++ ,对应的矩阵为单位矩阵E . 即()()TY A TY T EY Y T =,所以()EY Y Y AT T Y T T T =,故可证得A 合同于单位矩阵E . 充分性, 若A 合同于矩阵E ,则存在可逆矩阵B ,使得A =T B EB .任意取X≠0, BX Y ==()12,,T n y y y ,则有Y ≠0.于是有Y Y EBX B X AX X T T T T ===22212n y y y ++ >0,定理可以得证.定理5 实对称矩阵n n R A ⨯∈是正定矩阵的充分而且必要条件是A 的所有特征根都大于零.证明:必要性, A 为正定矩阵,若A 的全部特征值n λλλ,,,21 不全大于0,不妨设01≤λ. 则存在正交矩阵P 使得有⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=n TAP P λλλ21成立. 令(),,,,21n P ααα = 则有i i i A αλα=()n i ,,2,1 =,即i α为A 的属于特征值i λ的特征向量.特别的,取单位特征向量01≠β,即111βλβ=A .于是11111βλβββT T A =01≤=λ,而这与A 为正定矩阵相矛盾,所以A 的全部特征值n λλλ,,,21 都大于0.充分性,A 的特征值为n λλλ,,,21 ,则存在正交矩阵T ,使得⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡==-n T AT T AT T λλλ 211则有121-⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=T T A n λλλ. 任意取0≠X ,则有Y Y X T TX AX X n T T n T T ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=λλλλλλ2121, 其中T X Y T T =()0,,,21≠=n y y y ,于是得AX X T 02222211>+++=n n y y y λλλ ,即有A 为正定矩.定理6[3] 实对称矩阵n n R A ⨯∈为正定矩阵的充要条件是A 是半正定矩阵且0≠A . 证明:必要性, 因为A 是正定矩阵,则A 一定是半正定矩阵,且0≠A .充分性, 设A 的特征值为n λλλ,,,21 ,由于A 是半正定矩阵可知,i λ()n i ,,2,10 =≥,又021≠⋅⋅⋅=n A λλλ ,故()n i i ,,2,10 =>λ,所以A 是正定矩阵.定理7 实对称矩阵n n R A ⨯∈ 是正定矩阵的充分而且必要条件是存在实可逆矩阵C ,使得C C A T=.证明:必要性,若A 是实对称正定矩阵,则存在实可逆矩阵C 使得EC C A T =C C T =,其中E 为n 阶单位矩阵.充分性,因为存在实可逆矩阵C ,使得C C A T =,并且有C C A T =EC C T=,其中E 为n 阶单位矩阵.即实对称矩阵A 合同于E ,所以可得A 为正定矩阵.定理8 实对称矩阵n n R A ⨯∈为正定矩阵的充分而且必要条件是存在实列满秩矩阵n m P ⨯, 使P P A T =.证明:必要性, 因为A 为正定矩阵, 则存在n 阶实可逆矩阵C , 使得C C A T =()()n m n TnnC -⨯⨯=0()⎪⎪⎭⎫⎝⎛⨯-⨯n n m n n C 0. 令 =P ()⎪⎪⎭⎫⎝⎛⨯-⨯n n m n n C 0,则 P P A T =, 其中P 为n m ⨯列满秩矩阵. 充分性,n m P ⨯为实列满秩矩阵,则P P T 为n 阶可逆矩阵, 故对任意的n R X ∈,0≠X , 则由秩m C =, 知,0≠CX 并且有0)(>==PX PX PX P X AX X T T T T ,即A 为正定矩阵.定理9[4] 对称矩阵A 是正定矩阵的充分而且必要条件是对任意的实n 阶可逆方阵C ,使得AC C T 都是正定的.证明:必要性,首先()TT AC C AC C T =,对任意n R X ∈,0≠X ,由秩n C =,知,0≠CX 由于A为正定矩阵,故()()(),0>=CX A CX X AC C X TT T即AC C T 为正定矩阵.充分性,AC C T 正定,则对任意的n R X ∈,0≠X ,由秩C n =,知,0≠TX 并且()()CX A CX T =()0>X AC C X T T ,即可得A 为正定矩阵.定理10 实对称矩阵A 是正定矩阵的充分而且必要条件是存在实可逆上三角矩阵R ,使R R A T =.证明:必要性,由于A 是实对称正定矩阵,所以存在实可逆矩阵P ,使得P P A T =.且存在矩阵Q 和R 使得QR P =,其中Q 为n 阶正交矩阵,R 为n 阶主对角元素都大于零的上三角矩阵,从而有P P A T =QR Q R T T =R R T =.充分性,因为存在n 阶主对角元素都大于零的上三角矩阵R ,使得R R A T =. 则显然矩阵R 可逆, 即可证得A 是正定矩阵.定理11 实对称矩阵n n R A ⨯∈为正定矩阵的充分而且必要条件是存在n 阶主对角元素都大于零的下三角矩阵U ,使得U U A T =.(证明同上)2.3 利用主子式判定定理12 实对称矩阵nn R A ⨯∈ 是正定矩阵的充分而且必要条件是A 的各阶顺序主子式都大于零.证明:必要性, 因为A 是实对称正定矩阵,所以存在二次型()n x x x f ,,,21 ∑∑===ni nj j i ij x x a 11是正定的.且对于每个k ,n k ≤≤1令()k k x x f ,,1 ∑∑===ki kj j i ij x x a 11.对于任意一组不全为零的实数k b b ,,1 ,有()k k b b f ,,1 ∑∑===k i kj j i ij b b a 11=()0,,0,,,1 k b b f .0>所以()k k x x f ,,1 是正定的. 由正定矩阵的行列式大于零可知,k f 的行列式,01111>kkk k a a a a n k ,,1 =. 所以可证得矩阵A 的一切顺序主子式都大于0.充分性, 对n 作数学归纳法.当1=n 时, ().21111x a x f =由条件中011>a ,显然可得()1x f 是正定的. 假设对于1-n 元二次型成立,现在来证明n 元二次型的情形.令 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=----1,11,11,1111n n n n a a a a A ,=β⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-n n n a a ,11 , 于是矩阵A 可以分块写成A ⎥⎦⎤⎢⎣⎡=nn Ta A ββ1. 由于A 的顺序主子式全大于零,所以1A 的顺序主子式也全大于零. 由归纳法假设可以知道,1A 是正定矩阵,即存在可逆的1-n 阶矩阵P 使得11-=n T E P A P ,此处1-n E 可代表1-n 阶单位矩阵.令⎥⎦⎤⎢⎣⎡=1001P C , 则有⎥⎦⎤⎢⎣⎡100T P ⎥⎦⎤⎢⎣⎡nn T a A αα1⎥⎦⎤⎢⎣⎡100P ⎥⎦⎤⎢⎣⎡=-nn TT n a P P E αα1. 再令⎥⎦⎤⎢⎣⎡-=-1012αT n P E C , 则有2112C AC C C T T ⎥⎦⎤⎢⎣⎡-=-101P E T n α⎥⎦⎤⎢⎣⎡-nn T T n a P P E αα1⎥⎦⎤⎢⎣⎡--101αT n P E ⎥⎦⎤⎢⎣⎡-=-ααT T nn n PP a E 001.最后再令21C C C =, ,ααT T nn PP a a -=则有⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=a AC C T 11 . 两边同时取行列式,可有a A C =2.因为0>A ,所以0>a . 于是可得⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡a 11 =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡a 11 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡111 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡a 11 . 所以矩阵A 可与单位矩阵E 合同,并且可以证得矩阵A 是正定矩阵.定理13 实对称矩阵A 是正定矩阵的充分而且必要条件是A 的一切主子式均大于零.证明:必要性, (利用反证法)设A =()ij n n a ⨯是正定矩阵,假如可存在k 阶主子矩阵111212122212,0k k k k k k k ki i i i i i i i i i i i i i i i i i i i a a a a a a A A a a a =<则可根据k i A 是k 阶实对称矩阵,并由引理知可存在k 阶正交矩阵P ,使得P P A k T i k ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=βββ21 此处k βββ,,,21 为k i A 的特征值.由于k i A <0,且k i A =k βββ 21可知k i A 的特征值k βββ,,,21 中至少有一个小于0.推至一般性,设1β<0,令T Y =()1,0,,0 ,则可有Y ≠0并且k T i Y A Y =1u <0,再令T X =12(,,,)n x x x ,则有当{}12,,,k i i i i ∈ 时,可得i i x y =;当i 为其他时,得0i x =.则有X ≠0,且T X AX =k T i Y A Y =1u <0,而这与A 为正定矩阵的假设相矛盾.充分性, 假设k i A 是A 的一个k 阶主子矩阵, 则由于k i A 任意的一个顺序主子式均是A 的一个主子式,所以可知它们都大于0.所以可得k i A 为正定矩阵.定理可以得证.定理14[5] 实对称矩阵n n R A ⨯∈为正定矩阵的充分而且必要条件是A 的一切主子矩阵均为正定矩阵.证明:必要性,A 正定, 令⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=nn n n a a a a A 1111,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=k k k k i i i i i i i i k a a a a A 1111, 其中 k A 为A 的主子矩阵, n i i k ≤<<≤ 11()n k ,,2,1 =.显然 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=k k k k i i i i i i i i k a a a a A 1111()n k ,,2,1 =也是实对称矩阵.由于⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=k k k k i i i i i i i i k a a a a A 1111的k 个顺序主子式均为A 的k 个主子式,所以k 个主子式都大于零, 从而有⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=k k k k i i i i i i i i k a a a a A 1111()n k ,,2,1 =为正定矩阵.充分性,若实对称矩阵A 的一切主子矩阵均是正定矩阵,则矩阵A 的一切主子式全都大于零,即可证得A 是正定矩阵.2.4 其他常用判定定理15 若A 是实对称正定矩阵,则1-A 也是实对称正定矩阵. 证明:由于A 是实对称正定矩阵,则0>A ,所以A 可逆.又因()(),111---==A A A T T所以可得1-A 也是实对称矩阵.设A 的特征值为n λλλ,,,21 ,由A 正定有()n i i ,,2,10 =>λ,1-A 的全部特征值为01>iλ()n i ,,2,1 =,即可得1-A 为正定矩阵.定理16 若A 是实对称正定矩阵,则对于任意的整数m ,m A 都是正定矩阵. 证明:I 当0=m 时,显然是正定矩阵.II 当0<m 时,由于m m -=,则有()mm A A 1-=,且1-A 也是正定矩阵,故只需假定m 为正整数即可.(i )当m 为偶数时,由于A A T =,并且⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫ ⎝⎛=22m Tm m A A A ,所以可得m A 是正定的; (ii )当m 为奇数时,由于A 是正定矩阵,所以存在实可逆矩阵C ,使得C C A T=; 由此可得:2121212122----==⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛=m m m m CA C A AA A A A A A Tm m Tm m ⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛=--2121m m CA CA T从而m A 是正定矩阵.定理17 若A 是n 阶实对称正定矩阵,则有*A 也是正定矩阵(其中*A 表示A 的伴随矩阵).证明:已知*A =,1n n R A A ⨯-∈且()(),***==A A A T T又由于A 是正定矩阵,所以0>A .设A 的特征值为n λλλ,,,21 ,则由A 是正定矩阵有()n i i ,,2,10 =>λ,于是有*A 的n 个特征值11211,,,---n A A A λλλ 也都大于零,即可证得*A 也是正定矩阵.定理18 实对称正定矩阵的合同矩阵一定是正定矩阵. 证明:设实对称矩阵A 是正定矩阵,矩阵B 与矩阵A 合同,即存在可逆矩阵P ,使有AP P B T =成立,由于A 是正定矩阵,可知对于任意的n 维非零列向量X , 即nR X ∈≠0,有0>AX X T ,所以令PX Y =,则有0≠PX ,有0)()(>=CX A CX BY Y T T ,所以矩阵B 是正定矩阵,所以定理可得证.定理19 任意两个同阶实对称正定矩阵的和还是正定矩阵,更一般地,正定矩阵的正线性组合也是正定矩阵.证明:设A 、B n n R ⨯∈ 都是正定矩阵,同时又可设0,>b a , 因而对于任意的n R x ∈≠0, 可有0)(>+=+Bx bx Ax ax x bB aA x T T T .所以对于任意的两个同阶的正定矩阵的和仍是正定矩阵.而多于两个矩阵时,可以按照相同的方式进行处理, 并且可以利用数学归纳法给出具体的证明:(1)当2=n 时,由上可知命题结论成立;(2)假设当1+<k n 时有命题结论成立,以下可以证明1+=k n 时命题结论仍成立. 设121,,,+k k A A A A 是同阶的正定矩阵,并且有0,,,,121>+k k b b b b .下证1111+++++k k k k A b A b A b 也为正定矩阵.因而可得对于任意的n R x ∈≠0 有0)(11111111>+++=+++++++x A x b x A x b x A x b x A b A b A b x k T k k T k T k k k k T ,此式中的每一项均为正.所以可以得到当1+=k n 时, 结论成立.综合以上的(1)、(2)可知,对于一切的自然数n ,正定矩阵的正线性组合也仍为正定矩阵.定理20 对于任何的实对称矩阵A ,必存在实数0,0>>βα,使得A E α+与A E +β都是正定矩阵.证明:实对称矩阵A 的所有的特征根都是实数,所以不妨记其中一个绝对值最大的特征根为ολ,只要取οβλ>,则可有A E +β是正定矩阵.假设Q 是正交矩阵,使得⎪⎪⎪⎭⎫⎝⎛=n TAQ Q λλ 1则AQ Q EQ Q Q A E Q T T T +=+ββ)(=ββ⎛⎫ ⎪ ⎪ ⎪⎝⎭ +1n λλ⎛⎫ ⎪⎪ ⎪⎝⎭=1n βλβλ+⎛⎫⎪ ⎪ ⎪+⎝⎭由于0i βλ+>()1,2,,i n = ,可得A E +β也是正定矩阵.而当取1αβ=时,则有0α>,()1E A E A αββ+=+也是正定矩阵,于是定理可以得证.定理21 若A 、B 都是实对称矩阵,并且BA AB =,则AB 也必为正定矩阵. 证明:易知AB 的特征根均大于零,且有当AB BA =时,可有AB BA A B AB T T T ===)(,所以AB 又是对称矩阵,从而可得AB 是正定的.定理22 实对称矩阵=A ⎥⎦⎤⎢⎣⎡3221A A A A T为正定矩阵的充分而且必要条件是1A 和21123A A A A T --都是正定矩阵.证明:当1A 可逆时,有⎥⎦⎤⎢⎣⎡--E A A ET 1120⎥⎦⎤⎢⎣⎡3221A A A A T ⎥⎦⎤⎢⎣⎡--E A A E 0211⎥⎦⎤⎢⎣⎡-=-21123100A A A A A T 必要性, 若A 正定,那么1A 也正定,11-A 存在. 令⎥⎦⎤⎢⎣⎡-=-E A A E P 0211, 则P 可逆,所以AP P T 也正定.从而⎥⎦⎤⎢⎣⎡--21123100A A A A A T 为正定矩阵,因此它的主子矩阵1A 和21123A A A A T--都为正定矩阵.充分性,由于1A 和21123A A A A T --都是正定矩阵,且两个正定矩阵的和也是正定矩阵,可知 ⎥⎦⎤⎢⎣⎡--211231A A A A A T 为正定矩阵. 又可得⎥⎦⎤⎢⎣⎡=3221A A A A A T=()TP 1-⎥⎦⎤⎢⎣⎡--2112300A A A A A T 1-P ,即可证得A 为正定矩阵.定理23 实对称矩阵n n R A ⨯∈为正定矩阵的充分而且必要条件是存在正交的向量组n ααα,,,21 使得.2211Tn n T T A αααααα+++=证明:必要性,因为A 是正定矩阵,所以存在正交矩阵Q ,使得Q Q A n T ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=λλλ 21,T n Q ),,(21βββ =, 令 i i i βλα=()n i ,,2,1 =为正交向量组, 则可得.2211Tn n T T A αααααα+++=充分性,Tn n T T A αααααα+++= 2211= )(21T n TT ααα ⎪⎪⎪⎪⎪⎭⎫ ⎝⎛n ααα 21 T T T = (T 为正交矩阵), 显然可证得A 是正定矩阵.3 正定矩阵的应用3.1 用正定矩阵的定义来证明一些结论例 3.1 设A ,B 是n n ⨯实对称矩阵,A 是正定阵,证明:存在实可逆阵T ,使T B A T )(+'为对角阵.证 由于A 是正定阵,从而合同于单位阵E ,即可知存在实可逆阵Q ,使E AQ Q ='. 而BQ Q '仍是实对称矩阵,从而存在正交阵P ,使⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=''n P BQ Q P λλ 1)(,其中n λλ,,1 是BQ Q '的特征值,若令QP T =,则可有⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡++=+'n T B A T λλ11)(1 . 例 3.2 设B 为n 阶实对称矩阵,且正定. A 为m n ⨯实矩阵, T A 为A 的转置矩阵.试证:BA A T 为正定矩阵的充分而且必要条件是秩m A =)(.证 充分性 因为BA A BA A T T T =)(.0,1≠∈∀⨯x R x n ,由秩m A =,知()n j i a a ji ij ,...2,1,,==,而A 为正定阵,故0)()()(>=Ax B Ax x BA A x T T T ,此即BA A T 为正定阵.必要性 利用反证法.若秩m A <,则有0=Ax 有非零实数解0x 存在,即00=Ax ,但00≠x ,并且由BA A T 为正定矩阵,可知)()()(00000Ax B Ax x BA A x T T T=< ①另一方面,因为00=Ax ,所以m A =.0)()(00=Ax B Ax T ②由于①、②矛盾,故秩m A =)(.例 3.3 设A 是n 阶正定矩阵,B 是n 阶半正定矩阵,求证: A B A B +≥+,当且仅当0B =或n 1=时等号成立.证 由A 0>可知,存在n 阶的可逆矩阵P ,使得T P BP n E =成立,所以有()T T n P A B P E P BP +=+,且T T n P A B P E P BP +=+又因为T P BP 是半正定矩阵,设T P BP C ==()ij C ,则可有Tn E P BP +=11121212221211nnn n nnc c c c c c c c c ++=12121111n n n n n c c c c ---+++++其中i c 是C 的所有i 阶主子式之和,1,2,,i n = .而又因为0T C P BP =≥,并且它的所有主子式都是非负的,因此可得T n E P BP +≥1n +n c =n E +T P BP =T P AP +T P BP所以T P A B P +≥()TP A B P +由此可得A B A B +≥+当0B =或1n =时,显然有A B A B +≥+成立;当0B ≠且1n >时,易知T P BP C =0n n ⨯≠,于是可得至少有一个ij c ≠0,此时C 的一阶主子式ii c ,jj c 均不能为零,否则00ijijc c =2ij c -0<,这与C 是半正定矩阵矛盾.于是1c 0>,进一步可有T n E P BP +1>n c +,从而得A B A B +≥+成立.3.2 正定矩阵在数学分析上的应用3.2.1 多元函数的极值问题例3.4 求函数321232221321212),,(x x x x x x x x x f ++++=的极值.解 因为2211123x x x f +=∂∂,212212x x x f +=∂∂,2233+=∂∂x x f,令01=∂∂x f ,02=∂∂x f,03=∂∂x f ,得驻点T x )1,0,0(0-=,T x )1,144,24(1--=.又)(x f 的各二阶偏导数为12126x xf =∂∂,12212=∂∂∂x x f ,2312=∂∂∂x x f ,2222=∂∂xf ,0322=∂∂∂x x f ,2232=∂∂xf ,得(黑塞)矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=20202122126)(1x x H .在点0x 处,有⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=20202122120)(0x H ,而)(0x H 的顺序主子式:0det 1=H ,0144212120det 2<-==H ,0296)(det det 03<-==x H H ,因此)(0x H 不定,0x 不是极值点.在点1x 处,有⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=2020212212144)(1x H ,而)(1x H 的顺序主子式:0144det 1>=H ,014421212144det 2>==H , 0280220212212144det 3>==H , 故)(1x H 为正定矩阵,T x )1,144,24(1--=为极小值点,极小值为6913)1,144,24()(1-=--=f x f .3.2.2 正定矩阵在积分中的应用例3.5 证明:椭球体331j 11ij i j i a x x ==Ω=∑∑:的体积等于1/24/3,Aπ-其中()33ijA a ⨯=是正定矩阵.证明 A 是正定矩阵,∴∃正交矩阵T ,使得⎪⎪⎪⎭⎫⎝⎛=321λλλAT T T,0>i λ,)3,2,1(=i 为A 的特征值 令⎪⎪⎪⎪⎭⎫⎝⎛=---131211λλλB 作变换TBY y y y TB x x x X =⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=321321,则此变换的Jacobi 行列式为2121321)(--=====AB B T TB J λλλ13312321j 13()ij iji x a x xx x x A x x ==⎛⎫ ⎪= ⎪ ⎪⎝⎭∑∑=Y Y BY B Y ATBY T B Y AX X TTT T T T T =⎪⎪⎪⎭⎫⎝⎛==321λλλ 1/212312312311T T X AX Y Y dx dx dx dx dx dx Ady dy dy -Ω≤≤∴===⎰⎰⎰⎰⎰⎰⎰⎰⎰1/24/3Aπ-3.3 正定矩阵在运筹中的应用3.3.1 具有约束方程的最优化问题例 3.6 某地区计划明年修建公路x 百公里和创建工业园区y 百公顷,假设收益函数为xy y x f =),(,受所能提供的资源(包括资金、设备、劳动力等)的限制,x 和y 需要满足约束条件369422≤+y x ,求使),(y x f 达到最大值的计划数x 和y .解 由于约束方程369422=+y x 刻画的不是坐标平面上单位向量的集合,我们需要做变量变换.将这个约束方程写成1)2()3(22=+yx , 再设31x x =,22yx =,即13x x =,22x y =,则约束方程可以写成 12221=+x x ,而目标函数变成2121216)2)(3()2,3(x x x x x x f ==.现在的问题就成为求216)(x x x F =在1=x x T下的最大值,其中⎪⎪⎭⎫⎝⎛=21x x x .设⎥⎦⎤⎢⎣⎡=0330A ,则 Ax x x F T =)(,A 的特征值是3±.属于31=λ的单位特征向量是⎪⎪⎭⎫⎝⎛2121.由此得,当211=x ,212=x 时,)(x F 取得最大值3,即当12.22331≈==x x 百公里,41.1222≈==x y 百公顷时,收益函数),(y x f 去的最大值3.3.4 用正定矩阵来证明不等式例3.7 证明不等式2224222x y z xy xz ++>-(其中,,x y z 是不全为零的实数)证明 令⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛----=--++=z y x z y x xz xy z y x f 301051111),,(2235222则有⎪⎪⎪⎭⎫⎝⎛----=301051111P 的各阶顺序主子式为 01>,045111>=--,0731051111>=----, 所以P 是正定矩阵00,0x y z ⎛⎫⎛⎫ ⎪ ⎪∴∀≠ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭有0f >故可得原不等式成立.3.5 正定矩阵在几何中的应用3.5.1二次曲面的标准型 例3.8 在3R 中化简二次方程03828322620828102222=-++-+-++-z y x zx yz xy z y x ,并判断其曲面形状.解 二次项相应的对称矩阵为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=10410421410141A .A 的特征多项式为)18)(18)(9(+--=-λλλλI A ,特征值为91=λ,182=λ,183=λ,对应的单位特征向量构成的正交矩阵为⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=12221222131P .令⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡'''=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡z y x P z y x ,方程化为 0938316343222222=-'-'+'-'-'+'z y x z y x , 配方得1)34(2)31(2)31(222=+'-+'+-'z y x .令31-'=x X ,31+'=y Y ,34+'=z Z ,得122222=-+Z Y X ,故原方程表示的曲面为单叶双曲面.参考文献[1] 北京大学数学系几何与代数教研室前代数小组.高等代数(第三版)[M],北京:高等教育出版社,2003.[2] 线性代数/余长安编著.—武汉:武汉大学出版社,2010.1[3] 胡跃进.广义正定矩阵的一个不等式[J],阜阳师范学院学报(自然科学版),2001.18(1):10-11.[4] 张禾瑞,郝丙新. 高等代数(第三版)[M],北京:高等教育出版社,1983.[5] 钱吉林.高等代数解题精粹(修订版)[M],北京:中央民族大学出版社,2002.Properties and Applications of positive definite quadratic form Summary: Based on the matrix and vector tool, we study a kind of special function, quadratic form. However many quadratics in practical application are real quadratic form, with one of the most important class being positive definite quadratic form. This paper focuses on the positive definiteness and application of the real matrix. This paper presents several discrimination methods of the real symmetric positive definite matrix and important conclusions, which allow people to make better use of this tool in the positive definite matrix. The paper is divided into three chapters, the first chapter mainly describes the definition of the quadratic, positive definite quadratic form and the positive definite matrix; the second chapter cited several matrix discrimination method of the description positive definiteness; the third chapter simply list some examples to illustrate the application of the positive definiteness of a real matrix.Keyword: positive definite quadratic form positive definite matrixcharacteristic value necessary and sufficient condition real symmetric matrix。
毕业论文(设计)论文(设计)题目:正定二次型的判定及应用目录摘要 (1)ABSTRACT (2)第1章引言 (3)1.1 研究背景及意义 (3)第2章二次型 (4)2.1 二次型 (4)2.3 正定二次型与正定矩阵 (4)第3章正定二次型的判定及应用 (7)3.1 正定二次型的判别方法 (7)3.2 正定二次型在实际中的应用 (15)第4章结论 (18)参考文献 (19)致谢 (20)摘要在二次型中,正定二次型占有特殊的地位,本文总结了正定二次型的一些判断方法及其在证明不等式与极值问题中的应用。
关键词:正定二次型;正定矩阵;顺序主子式;ABSTRACTIn the quadratic form,the positive definite quadratic form has a special position.This paper has summarized some judjement methods of the positive definite quadratic form and given some applications in inequalities proving and extreme problems.Key words: positive definite quadratic; positive definite matrix; principal minor determinant第1章引言1.1 研究背景及意义在数学中,二次型的理论起源于解析几何中化二次曲线和二次曲面方程为标准形的问题. 二次型的系统研究是从18世纪开始的,柯西在其著作中给出结论:当方程式标准型时,二次曲面用二次型的符号来进行分类.然而,那是并不太清楚,在化简成标准型时,为何总是得到同样数目的正项和负项.西尔维斯特回答了这个问题,他给出了n个变数的二次型的惯性定律,但没有证明.这个定律后被雅克比重新发现和证明.1801年,高斯在《算术研究》中引进了二次型的正定、负定、半正定和半负定等术语.二次型化简的进一步研究涉及二次型或行列式的特征方程的概念。