正定矩阵和半正定矩阵的性质及应用
半正定矩阵的性质

半正定矩阵的性质内容摘要矩阵是线性代数的一个重要内容,矩阵这一概念是从其它许多事物中抽象出来的,具有很大的现实意义.矩阵的理论不仅贯穿于线性代数的各个部分,而且在在物理学及其它科学技术领域,在经济及其它社会科学领域都有广泛的应用.本文以半正定矩阵的概念为基本出发点,从特征值、主子式、QR 分解、Gram 矩阵、半正定矩阵的各种运算等等系统研究半正定矩阵的基本性质,尤其是hadamard 积 和kronecker 积 ,更深刻的理解半正定矩阵的内涵和性质.【关键词】 半正定矩阵 hadamard 积 kronecker 积 一、矩阵的相关知识定义1[1]. 矩阵的秩向量组的极大线性无关组所含向量的个数称为向量组的秩.矩阵的行秩就是矩阵行向量组的秩,矩阵的列秩就是矩阵列向量组的秩.矩阵的行秩等于矩阵的列秩,统称为矩阵的秩,记作()R A .定义2[1]. 矩阵的特征值与特征向量设A 是n 阶方阵,如果存在数λ和非零n 维列向量x ,使得x Ax λ=成立,则称 λ 为A 的一个特征值.非零n 维列向量x 称为矩阵A 的属于(对应于)特征值λ的特征向量,简称A 的特征向量.特征向量0≠x .注1. 特征向量不是由特征值唯一确定的,但是特征值都是由特征向量唯一决定的.所以一个特征向量只能属于一个特征值,一个特征值有无穷多个特征向量.注2.对于一个n 阶矩阵A ,λ是矩阵A 的特征值,一般通过求解特征方程A E f -=λλ)(和齐次线性方程组()0E A X λ-=来得到矩阵的特征值和特征向量.定义3[1]. 矩阵的迹设矩阵()ij n n A a ⨯=,那么矩阵A 的迹就是矩阵A 的主对角线元素的之和,记作()tr A .注3.矩阵的迹就是矩阵的所有特征值之和. 定义4. 对角优势矩阵 对于矩阵()ij n n A a ⨯=,如果1nii ij j j ia a =≠≥∑,1,2,,i n =则称矩阵A 为对角优势矩阵.定义5[1].对称矩阵对于矩阵()ij n n A a ⨯=,若元素满足ji ij a a =, n j i 2,1,=或者A A T =, 则称矩阵A 为对称矩阵.定义6[2].酉矩阵对n 阶复矩阵A ,用A -表示以A 的元素的共轭复数作元素的矩阵.如A 满足TTA A A A E --==,则称矩阵A 为酉矩阵.定义7.Gram 矩阵 设12,,,n v v v 是欧氏空间V 的一个向量组,定义矩阵111212122212,,,,,,,,,n n n n n n v v v v v v v v v v v v A v v v v v v <><><>⎛⎫⎪<><><> ⎪= ⎪⎪<><><>⎝⎭A 称为由向量12,,,n v v v 组成的Gram 矩阵,记做()12,,,n Gram v v v . 其中,,<⋅⋅>为欧氏空间V 中定义的内积.定义8. 可对角化如果方阵A 相似于一个对角矩阵,称方阵A 为可对角化,换句话说,即如果存在一个可逆矩阵 P 使得AP P 1-是对角矩阵,那么称矩阵A 可对角化.定义9[2]. 置换矩阵对于矩阵()ij n n P p ⨯=,如果它的每一行和每一列都只有一个元素为1,其它的元素都为零,则称矩阵P 为置换矩阵.定义10[2]. 可约矩阵 对于 矩阵()n n ij a A ⨯=,如果满足 ①1=n 时,0=A ;②2≥n ,存在n 阶置换矩阵P ,使得⎪⎪⎭⎫ ⎝⎛=-D O C B PAP 1,其中B 是k 阶方阵11-≤≤n k ,左下角是()k k n ⨯-阶的零矩阵,则称矩阵A 为可约的.否则,矩阵A为不可约.定义11[1]. 非退化矩阵对于矩阵()n n ij a A ⨯=,如果0≠A ,称矩阵A 为非退化的. 定义12[1]. 矩阵的幂对于矩阵()n n ij a A ⨯=,对任意正整数k ,kA 定义为k kA AAA=,称为矩阵A 的k 次幂.规定E A =0.定义13[4]. 阵的QR 分解实(复)非奇异矩阵A 能够化成正交(酉)矩阵Q 与实(复)非奇异上三角矩阵R 的乘积,即QR A =,称为A 的QR 分解.定义14[5].Kronecker 积设()n m ij R a A ⨯∈=,()n m ij R b B ⨯∈=,A 与B 的Kronecker 积,记作B A ⊗,定义为n m nn n n n n C B a Ba B a B a Ba B a B a Ba B a B A ⨯∈⎪⎪⎪⎪⎪⎭⎫⎝⎛=⊗ 212222111211注4:从定义看矩阵的Kronecker 积被表示成为矩阵的分块运算,即B A ⊗是一个分块矩阵,每一个子块是数乘运算()B a ij . 矩阵的Kronecker 积也称为直积或张量积.定义15[5].Hadamard 积设()n m ij R a A ⨯∈=,()n m ij R b B ⨯∈=,其中A 与B 为同阶矩阵,A 与B 的Hadamard 积,记作A B ,定义为()n m ij ij R b a B A ⨯∈= .注5: 矩阵A 与B 的Hadamard 积即将A 与B 对应元素相乘,矩阵的Hadamard 积也称为Schur 积.注6[5]:矩阵Hadamard 积的性质:① ()()kB A B kA = ② ()C A B A C B A +=+ ③ ()()C B A C B A =④()TT T A B A B =注7:由矩阵的Kronecker 积与Hadamard 积的定义可以看出,B A 是B A ⊗的主子矩阵.二、半正定矩阵的性质(一)半正定矩阵的定义如果矩阵n n R A ⨯∈是实对称矩阵,并且对于一切n R X ∈,有0≥AX X T ,则称矩阵A 为半正定矩阵.记作0A ≥.如果0≥-B A ,记作B A ≥.(二)半正定矩阵的二次型对称矩阵A 的二次型()AX X X f T =,如果对任何非零向量X ,都有0≥AX X T 成立,则称()AX X X f T =为半正定二次型.(三)半正定矩阵的性质性质1:设A 为一个n 阶对角优势对称矩阵,且对角线元素非负,那么矩阵A 为半正定矩阵.证明:设A 是一个n 阶对角优势对称矩阵,且对角线元素非负. 设ij A 矩阵是A 的主子矩阵且i 行j 列如下,且其它对角元素等于0⎪⎪⎭⎫ ⎝⎛ij ij ij ij a a a a那么ij A 是一个半正定矩阵,而∑=-nj i ij A A 1,也为非负对角阵,所以A 为半正定的.注8:半正定矩阵对角化之后,所有元素都大于等于0. 性质2:设A 为一个n 阶对称矩阵,下列命题等价:(a)A 是半正定矩阵;(b)A 的所有特征值为非负; (c)A 的所有的主子式非负;(d)存在一个n 阶矩阵B ,使得T BB A =; (e)存在一个n 阶下三角矩阵L ,使得T LL A =; (f)存在一个n 阶对称矩阵C ,使得2C A =;(g)存在一个k 维欧氏空间V 和向量12,,,n v v v V ∈,使得()12,,,n A Gram v v v =;(h)存在k 个向量12,,,nk b b b R ∈,使得∑==ki T i i b b A 1.证明:(a)⇒(b)设,AX X λ=0X ≠,其中λ为矩阵A 的特征值,由于矩阵A 为半正定矩阵,有0≥=X X AX X T T λ,且0T X X >,则\0T T X AX X X λ=≥, 所以矩阵A 的所有特征值非负.(a)⇒(c) 设[]a A 是A 的主子式,由于A 为半正定的,所以[]a A 也为半正定的,由(a)⇒(b)可知[]a A 的特征值为非负,因此,[]0≥a A .(c)⇒(b) 设A 的特征多项式()()()n nk n k kn n n A P x P x P xP x x 112211-++-+-+-=∆--- 其中k P 为A 的所有k k ⨯阶子矩阵的和,由于 (c),n k P k 2,10=≥,,假设0<x ,如果n 为任意正整数,那么0>n x 并且()0≥∆x A ;如果n 唯一,那么0<n x ,并且()0≥∆x A ,这表明矩阵A 不可能有负特征值且A 为对称矩阵,所以矩阵A 特征值存在且非负.(b)⇒(f) 由于A 为对称矩阵,并且特征值非负,它正交相似与一个非负对 角矩阵D .即T UDU A =,其中U 是正交矩阵,D 是非负对角矩阵()n d d d diag D 21=,但是当T T U D U U D U A =,其中由于U 为正交矩阵,所以有T U U E =,然而),n diagd =.所以T U D U C C A ==,2.(d)⇒(e) 为了证明这个结论,首先利用下列这一点:任何矩阵C 有一个QR 因数i.e.,QR C =,其中Q 的行正交,R 是上三角矩阵.设TA BB =,TB 的QR 分解为T B QR =,有()()TTTT T B B QR R Q ===,那么T T T A R Q QR LL ==,那么就有T T Q R L =,T QR L =,然而TR L =是一个下三角矩阵.所以T LL A =.(d)⇒(g) 设T BB A =,然而B 是n n ⨯阶矩阵,设k R V =,并且设T i V 是 B 的i 行,那么()n v v v Gram A 21=(g)⇒(a)由于()12,,,n A Gram v v v =,且设n R x ∈,那么()()0x v ,211i n 1j 1,1,,1,≥=⎪⎪⎭⎫ ⎝⎛====∑∑∑∑∑∑======ni ii n j j i i nj i j j i i j i nj i j i nj i j i ij Tvx v x v x v x x x v v x x a AX X 0T X AX ≥,所以矩阵A 为半正定矩阵.(b)⇒(h) 设T BB A = ,则有∑==ki T i i b b A 1,其中(1,2)i b i k = 为B 的列向量,由(e)⇒(d),(f)⇒(d),可知结论成立.注9:①性质2 中,证明(b)⇒(f) 中,构造的矩阵C 其实为半正定矩阵,这表明任何一个半正定矩阵A 都有唯一的半正定矩阵C 满足2C A =,那么矩阵C 为A 的平方根,记作C =②中指的是主子矩阵而不是顺序主子式,实际上,只有顺序主子式大于等于零并不能保证A 是半正定的.例1 判定二次型()2123221321245,,x x x x x x x x f -++=的正定性 解 (解法1) 用顺序主矩阵判别 首先,该二次型()123,,f x x x 对应的矩阵为⎪⎪⎪⎭⎫⎝⎛--=100041015A求A 的各阶主子式可得,051>=D 02241152>=--=D 03==A D由于A 的各阶主子式全都大于等于零,所以该二次型()123,,f x x x 为半正定二次型。
正定矩阵的性质和判定方法及应用

内蒙古财经大学本科毕业论文正定矩阵的性质及应用作者郝芸芸系别统计与数学学院专业信息与计算科学年级10级学号102093113指导教师高菲菲导师职称讲师答辩日期成绩内容提要矩阵是数学中的一个重要基本概念,也是一个主要研究对象,同时矩阵论又是研究线性代数的一个有力工具.而矩阵的正定性是矩阵论中的一个重要概念.正定矩阵是一种特殊的矩阵,其等价定理在解题过程中可以灵活使用.且正定矩阵具有一般矩阵不具有的特殊性质,尤其是这些性质广泛地应用于各个领域.本文在第一部分介绍了实矩阵的正定性的相关定义以及其等价条件.在第二部分列举了正定矩阵的一系列性质,主要介绍了正定矩阵的关联矩阵的正定性.本文在第三部分介绍了正定矩阵的相关定理.本文在第四部分介绍了矩阵正定性的判定方法:定义法、主子式法、特征值法、与单位矩阵合同法.且简单地举了一些实例来阐述实矩阵正定性的判定.最后本文分别从不等式的证明和多元函数的极值两个方面介绍了正定矩阵的实际应用.关键词:二次型正定矩阵判定方法应用AbstractMatrix is an important basic concepts in mathematics,but also a main research object,at the same time matrix theory is a powerful tool for the study of linear algebra。
At the same time,the positive definiteness of matrix is an important concept in the matrix theory。
The positive definite matrix is a special matrix, the equivalence theorem in the problem solving process can be used flexibly。
正定矩阵的性质和判定方法及应用

内蒙古财经大学本科毕业论文正定矩阵的性质及应用作者郝芸芸系别统计与数学学院专业信息与计算科学年级10级学号*********指导教师高菲菲导师职称讲师答辩日期成绩内容提要矩阵是数学中的一个重要基本概念,也是一个主要研究对象,同时矩阵论又是研究线性代数的一个有力工具.而矩阵的正定性是矩阵论中的一个重要概念.正定矩阵是一种特殊的矩阵,其等价定理在解题过程中可以灵活使用.且正定矩阵具有一般矩阵不具有的特殊性质,尤其是这些性质广泛地应用于各个领域.本文在第一部分介绍了实矩阵的正定性的相关定义以及其等价条件.在第二部分列举了正定矩阵的一系列性质,主要介绍了正定矩阵的关联矩阵的正定性.本文在第三部分介绍了正定矩阵的相关定理.本文在第四部分介绍了矩阵正定性的判定方法:定义法、主子式法、特征值法、与单位矩阵合同法.且简单地举了一些实例来阐述实矩阵正定性的判定.最后本文分别从不等式的证明和多元函数的极值两个方面介绍了正定矩阵的实际应用.关键词:二次型正定矩阵判定方法应用AbstractMatrix is an important basic concepts in mathematics, but also a main research object, at the same time matrix theory is a powerful tool for the study of linear algebra. At the same time, the positive definiteness of matrix is an important concept in the matrix theory. The positive definite matrix is a special matrix, the equivalence theorem in the problem solving process can be used flexibly. And the positive definite matrix with special properties of general matrix does not have these properties, especially widely used in various fields. In the first part of this thesis introduces the related definition of positive definite real matrix and its equivalent conditions. In the second part are held a series of properties of positive definite matrix, mainly introduced the positive definiteness correlation matrix is positive definite matrix. This paper introduces the related theorem of positive definite matrix in the third part. This paper introduces the method to judge the positive definiteness matrix in fourth parts: the definition, the master method, the eigenvalue method. Determination and simply cited a number of examples of real positive definite matrices. Two aspects of extreme finally this paper from the proof of inequality and multiple function describes the practical application of positive definite matrices.Key words:Quadratic form Positive definite matrix Determination method Application目录引言...................................................... 错误!未定义书签。
正定矩阵和半正定矩阵

正定矩阵和半正定矩阵正定矩阵和半正定矩阵是矩阵理论中非常重要的一部分,它们经常被用于一些科学计算以及应用到数学建模中。
本文将从矩阵的定义出发,讨论正定矩阵和半正定矩阵的基本特性以及它们应用于数学建模中的重要作用。
全文共有三个部分,分别介绍矩阵的定义、正定矩阵和半正定矩阵的特性以及它们的应用。
首先,从最基本的概念开始,矩阵是一种用于表达数学模型的数据结构,由一系列有序的数字组成。
可以看到,它们有着m行n列(m ×n)的结构,即被称为m×n矩阵。
它们可以表示线性空间中的相关概念,如点、线、平面等。
在线性代数中,矩阵的几何意义很重要,它可以用来表示矢量和空间中的变换,并能把相关运算转化为矩阵乘法的关系。
接下来要讨论的是正定矩阵和半正定矩阵,首先,它们都是实矩阵。
正定矩阵是一种特殊的实矩阵,它有两个重要的性质:(1)它的行列式不等于0,(2)对于任意n×n正定矩阵A,有A=A*A,其中A*表示A的共轭转置。
由此可见,当A是一个正定矩阵时,A正定矩阵的转置也是正定矩阵。
半正定矩阵和正定矩阵有一些相似之处:(1)它们也是实矩阵,(2)它们的行列式也不等于0,(3)它们也具有自身的转置矩阵。
然而,半正定矩阵的性质未必和正定矩阵相同,它只是满足A=A*A,其中A*表示A的共轭转置,转置矩阵并不一定是正定矩阵。
最后,让我们来看看正定矩阵和半正定矩阵在数学建模中的重要作用。
正定矩阵和半正定矩阵常被用于统计回归分析,这是一种用于预测事件结果的模型,可以用来预测变量之间的关系。
正定矩阵和半正定矩阵也可以用来解决优化问题,比如最小二乘法的求解。
另外,它们还可以用来提取线性态度的代数体系,会非常有效地提升计算速度、准确率以及可解释性。
综上所述,正定矩阵和半正定矩阵是矩阵理论中重要的一部分,它们具有行列式不为零及正定矩阵的转置也是正定矩阵等重要性质,并且可以用于统计回归分析、解决优化问题以及提取线性态度的代数体系等,在数学建模中发挥着重要作用。
正定矩阵的性质及应用

§1 0 0 ·
例
1
矩阵
A=
¨ ¨
0
2
-1¸¸ 的正惯性指数为 3,所以 A 是正定
¨© 0 -1 3 ¸¹
矩阵.矩阵
§0 1 1 ·
B=
¨ ¨
1
0
-3¸¸
¨© 1 -3 0 ¸¹
的秩为 3,但正惯性指数为 2≠3,即 B 不能合同于 I3,所以 B 不是正定矩阵。
例 2 证明:如果 A 是正定矩阵,那么 A-1 也是正定矩阵.
利用构造法巧解高中数学问题
王运行
(兰州新区舟曲中学 甘肃 兰州 730087)
【摘要】在中学数学中,构造法在技巧与方法中占据着非常重要的地位,它可以起到化繁为简,化难为易的作用,将中学数学中的技巧性展示的 淋漓尽致。下面笔者将从一些常见的数学问题中来阐述构造法的具体应用。
【关键词】构造 ;转化 ;中学数学解题应用 【中图分类号】G633.6 【文献标识码】B 【文章编号】2095-3089(2017)16-0266-02
aij xi x j 是正定二次型。
i1 j1
证明:必要性 设 A 是正定矩阵,则 A 与 In 合同,即存在 n 阶可逆矩阵 P,使得
令 那么
PT AP In X PY
f x1, x2,", xn X T AX PY T A PY
= Y T PT APY Y T InY
正定矩阵和半正定矩阵的性质及应用-毕业论文

---文档均为word文档,下载后可直接编辑使用亦可打印---摘要本文主要针对正定矩阵和半正定矩阵进行讨论,归纳和总结了正定矩阵和半正定矩阵的性质,通过实例介绍了正定矩阵(半正定矩阵)的判别方法诸如:定义法、主子式法、特征值法等,并且给出了它们在不等式的证明问题中以及多元函数极值问题中的一些应用.关键词:正定矩阵;半正定矩阵;二次型;主子式;特征值ABSTRACTThis paper mainly discusses positive definite matrices and positive semi-definite matrix,the properties of positive definite matrix and semi-positive definite matrix are summarized.Through examples, the judgment methods of positive definite matrix and semi-positive definite matrix are introduced, such minor method, master type method, eigenvalue method, etc. Some applications of positive definite matrices and semi-positive definite matrix in the proof of inequality extreme value problems of multivariate functions are given.Keywords:positive definite matrix; positive semi-definite matrix; quadratic form; principal minor determinant;characteristic value第1章 正定矩阵和半正定矩阵的定义及性质1.1. 相关概念定义1[1] 设a ij (i =1,2,⋯,n ,i ≤j)都是实常数,则关于n 个实变量x 1,x 2,⋯,x n 的二次齐次多项式函数f(x 1,x 2,⋯,x n )=a 11x 12+a 22x 22+⋯+a nn x n 2+2a 12x 1x 2+2a 12x 1x 3+⋯+2a n−1,n x n−1x n ,称为n 元实二次型.[9]定义2[1] 实二次型f(x 1,x 2,⋯,x n )为正定的,如果对于一组不全为零的实数c 1,c 2,⋯,c n 都有f(c 1,c 2,⋯,c n )>0,如果都有f(c 1,c 2,⋯,c n )<0,那么称f(x 1,x 2,⋯,x n )为负定的.如果都有f(c 1,c 2,⋯,c n )≥0,那么称f(x 1,x 2,⋯,x n )为半正定的.如果都有f(c 1,c 2,⋯,c n )≤0,那么称f(x 1,x 2,⋯,x n )为半负定的.如果二次型既不是半正定又不是半负定,那么称为不定的.[1]定义3[1] 若实数域上的n 元二次型f(x 1,x 2,⋯,x n )=∑∑a ij nj=1ni=1x i x j =X T AX是正定(半正定)二次型,则A 被称为正定(半正定)矩阵,其中A =(a 11a 12⋯a 1n a 21a 22⋯a 2n ⋮⋮⋱⋮a n1a n2⋯a nn ),X =(x 1x 2⋮x n) 定义4[1] 子式|a 11a 12⋯a 1i a 21a 22⋯a 2i ⋮⋮⋱⋮a i1a i2⋯a ii|,,,,,,,,,0)00()(11121>==∑∑==ki kj k j i ij k k c c f c c a c c c f 称为矩阵A =(a ij )的i 阶顺序主子式i =1,2,⋯,n.1.2. 正定矩阵和半正定矩阵的等价命题1.2.1.正定矩阵的等价命题定理1[9] A 是n 阶实对称矩阵,则下列叙述等价: (1) A 是正定矩阵.(2) A 的所有顺序主子式全大于零. (3) A 的特征值全大于零. (4) 存在正定矩阵B ,使得A =B 2. (5) A 合同于E .(6) A 的一切主子式全大于零. (7) A 的一切主子矩阵都是正定矩阵.(8) 对任意的实列满秩矩阵C n×m ,都有C T AC 为正定矩阵.(9) 任意实可逆矩阵T ,都有T T AT 为正定矩阵. (10) 存在秩为n 的m ×n 实矩阵C 使A =C T C . (11) A =P T P ,P 是n 阶可逆矩阵.(12) A =R T R ,R 是n 阶主对角元素全大于零的上三角形矩阵.A =U T U ,U 是n 阶主对角元素全大于零的下三角形矩阵.[9]证明(1)⇒(2)设二次型f(x 1,x 2,⋯,x n )=∑∑a ij x i nj=1n i=1x j 是正定的.对于每个k ,1≤k ≤n ,设f k =(x 1,x 2,⋯,x n )=∑∑a ij x i n j=1n i=1x j以下证明f k 是一个正定二次型,对于任意一组不全为零的实数c 1,c 2,⋯,c k ,有)13(因此f k(x1,x2,⋯,x n)是正定的.由性质1可得,f k所对应的矩阵行列式|a11⋯a1n ⋮⋱⋮a k1⋯a kk|>0,k=1,2,⋯,n.从而A的各阶顺序主子式都大于零.(1)⇒(3)用反证法,若A的特征根λ1,λ2,⋯,λn不都大于零.不妨设λ1≤0,取A的属于λi的单位特征向量β≠0,就有βT Aβ=λ1≤0,这与A为正定矩阵相矛盾,所以A的特征值全大于零。
正定矩阵在函数极值问题中的应用

正定矩阵在函数极值问题中的应用函数极值问题是数学分析中的重要内容之一,它关注的是函数在其定义域内的最大值和最小值。
在实际生活和科学研究中,很多问题可以通过函数极值问题的求解来得到最优解。
正定矩阵是一个重要的工具,在函数极值问题中有着广泛的应用。
本文将从正定矩阵的定义、性质以及在函数极值问题中的应用三个方面进行阐述。
正定矩阵是一个n阶的对称矩阵,并且对于任意非零向量x,有x^TAx>0,其中^T表示向量的转置运算。
正定矩阵有许多重要的性质,其中最重要的是其特征值的性质。
正定矩阵的特征值都大于0,这使得正定矩阵具有良好的性质,在函数极值问题中具有重要的作用。
首先,正定矩阵可以用于判定函数的极值。
对于一个二次型函数f(x)=x^T Ax,其中x是一个n维向量,A是一个n阶正定矩阵。
如果A 正定,那么f(x)的极值为min或max,且极值点对应的向量是矩阵A的特征向量。
这意味着我们只需要求解正定矩阵A的特征值和特征向量,就能得到函数f(x)的极值点和极值。
其次,正定矩阵在最小二乘法中有着广泛的应用。
最小二乘法是一种常用的数学处理方法,用于求解数据拟合问题。
在最小二乘法中,我们使用一个关于参数β的线性模型y=Xβ+ε来拟合一组观察数据y,其中X 是给定的设计矩阵,ε是误差项。
最小二乘法的目标是寻找一个最佳的参数β使得误差项的平方和最小。
使用正定矩阵可以对最小二乘法进行改进,得到更精确的拟合结果。
具体来说,我们可以将最小二乘法表述为一个二次型的最小化问题。
令误差项ε为残差向量r=y-Xβ,其中β是待求的参数向量。
将残差向量r视为一个n维列向量,X视为一个n×k维矩阵,则残差的平方和可以表示为r^Tr=r^T(y-Xβ)。
(y-Xβ)是残差的向量表示形式,因此r^T(y-Xβ)就是残差的二次型表示形式。
如果我们将正定矩阵A设为X^TX,则最小化残差的平方和的问题可以转化为一个正定矩阵的二次型函数的最小化问题。
关于正定矩阵的性质及应用的研究

,得证。
性质3 、 是正定矩阵,则 证明: 、 是正定矩阵,所以
是实对称矩阵。
对任意的 维列向量 ,
,
,其中 ,因 是任意的,所以
也是正定矩阵。
,
,有
学术研讨 135
,
,所以,
也是正定矩阵,得证。
性质4 是正定矩阵,则
、 、 也是正定矩
阵。
证明: 是正定矩阵,故
,,
, 是实对称矩阵。
若 是 的特征值,则 是 的特征值。由
134
◇朔州师范高等专科学校 董改芳
关于正定矩阵的性质及应用的研究
2019 年 第 6 期
正定矩阵是高等代数矩阵理论中非常重要的内容,本文给出了正定矩阵的一些性 质和判定方法,并在实例中得到了正定矩阵的一些应用。
二次齐次多项式在数学的其它分支、物理以及力学中常常用到,是一类非常重要的多
项式。二次型是数域上的二次齐次多项式,在讨论二次型时,我们把二次型
采用钛酸四丁酯和四氯化钛为钛源可以将TiO2负载在玄武岩纤维 表面,但结合XRD分析,负载型的TiO2可能呈高度分散状态或 者无定形态存在。
4 结论 本文分别以钛酸四丁酯和四氯化钛为钛源,采用湿法化学 发在玄武岩纤维表面负载一层TiO2,采用X射线衍射仪和金相分 析仪对TiO2的负载情况进行初步的探索。实验结果表明,TiO2在 玄武岩纤维表面负载均匀,并且以无定形或者高分散状态存 在。
【参考文献】 [1] 胡显奇, 申屠年. 连续玄武岩纤维在军工及民用领域的应 用[J]. 高科技纤维与应用, 2005, 30(6): 7-13 [2] 曹海琳, 郎海军, 孟松鹤. 连续玄武岩纤维结构与性能试 验研究[J]. 高科技纤维与应用, 2007, 32(5): 8-13 [3] 姚勇, 徐鹏, 刘静, 等. 国内外玄武岩纤维耐腐蚀性能对比 研究[J]. 合成纤维工业, 2015, 38(5): 9 [4] Sim J, Park C. Characteristics of basalt fiber as a strengthening material for concrete structures[J]. Composites Part B: Engineering, 2005, 36(6): 504-512 [5] 王广健, 尚德库, 胡琳娜, 等. 玄武岩纤维的表面修饰及生 态环境复合过滤材料的制备与性能研究[J]. 复合材料学报, 2004, 21(1): 38-44 [6] 董丽茜, 陈进富, 郭春梅,等. 玄武岩纤维在环保领域的应 用研究现状及展望[J]. 当代化工, 2018(2) [7] 余娟, 周蓉, 邢建民. 耐高温针刺毡脱硝催化剂负载预处 理工艺探讨[J]. 山东纺织科技, 2018, 59(2): 1-5 [8] 耐高温玄武岩覆膜滤料的制备与性能的研究[D]. 浙江理 工大学, 2013 [9] 强降解VOC纳米TiO2光催化剂的制备及机理研究[D]. 华 中科技大学, 2015 基 金 项目:1、国家级大学生创新创业训练计划项目 (201810649050);2、乐山师范学院引进教师科研启动项目 (Z16024);3、乐山市科技重点研究项目(17GZD051)。 通讯作者:徐要辉,男,工学博士,乐山师范学院讲师, 主要从事功能材料的研究。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
摘要本文主要针对正定矩阵和半正定矩阵进行讨论,归纳和总结了正定矩阵和半正定矩阵的性质,通过实例介绍了正定矩阵(半正定矩阵)的判别方法诸如:定义法、主子式法、特征值法等,并且给出了它们在不等式的证明问题中以及多元函数极值问题中的一些应用.关键词:正定矩阵;半正定矩阵;二次型;主子式;特征值ABSTRACTThis paper mainly discusses positive definite matrices and positive semi-definite matrix,the properties of positive definite matrix and semi-positive definite matrix are summarized.Through examples, the judgment methods of positive definite matrix and semi-positive definite matrix are introduced, such minor method, master type method, eigenvalue method, etc. Some applications of positive definite matrices and semi-positive definite matrix in the proof of inequality extreme value problems of multivariate functions are given.Keywords:positive definite matrix; positive semi-definite matrix; quadratic form; principal minor determinant;characteristic value目录第1章正定矩阵和半正定矩阵的定义及性质 (1)1.1. 相关概念 (1)1.2. 正定矩阵和半正定矩阵的等价命题 (2)1.2.1. 正定矩阵的等价命题 (2)1.2.2 半正定矩阵的等价命题 (6)1.3. 正定矩阵和半正定矩阵的性质 (8)1.3.1. 正定矩阵的性质 (8)1.3.2. 半正定矩阵的性质 (13)第2章正定矩阵和半正定矩阵的判定方法 (15)2.1. 定义法 (15)2.2. 主子式法 (15)2.4. 与单位矩阵E合同法 (18)第3章定矩阵和半正定矩阵的应用 (20)3.1. 在不等式问题中的应用 (20)3.2. 在多元函数极值问题中的应用 (21)参考文献 (25)致谢 (26)第1章 正定矩阵和半正定矩阵的定义及性质1.1. 相关概念定义1[1] 设a ij (i =1,2,⋯,n ,i ≤j)都是实常数,则关于n 个实变量x 1,x 2,⋯,x n 的二次齐次多项式函数f(x 1,x 2,⋯,x n )=a 11x 12+a 22x 22+⋯+a nn x n 2+2a 12x 1x 2+2a 12x 1x 3+⋯+2a n−1,n x n−1x n ,称为n 元实二次型.[9]定义2[1] 实二次型f(x 1,x 2,⋯,x n )为正定的,如果对于一组不全为零的实数c 1,c 2,⋯,c n 都有f(c 1,c 2,⋯,c n )>0,如果都有f(c 1,c 2,⋯,c n )<0,那么称f(x 1,x 2,⋯,x n )为负定的.如果都有f(c 1,c 2,⋯,c n )≥0,那么称f(x 1,x 2,⋯,x n )为半正定的.如果都有f(c 1,c 2,⋯,c n )≤0,那么称f(x 1,x 2,⋯,x n )为半负定的.如果二次型既不是半正定又不是半负定,那么称为不定的.[1]定义3[1] 若实数域上的n 元二次型f(x 1,x 2,⋯,x n )=∑∑a ij nj=1ni=1x i x j =X T AX是正定(半正定)二次型,则A 被称为正定(半正定)矩阵,其中A =(a 11a 12⋯a 1n a 21a 22⋯a 2n ⋮⋮⋱⋮a n1a n2⋯a nn ),X =(x 1x 2⋮x n) 定义4[1] 子式|a 11a 12⋯a 1i a 21a 22⋯a 2i ⋮⋮⋱⋮a i1a i2⋯a ii|,,,,,,,,,0)00()(11121>==∑∑==ki kj k j i ij k k c c f c c a c c c f 称为矩阵A =(a ij )的i 阶顺序主子式i =1,2,⋯,n.1.2. 正定矩阵和半正定矩阵的等价命题1.2.1.正定矩阵的等价命题定理1[9] A 是n 阶实对称矩阵,则下列叙述等价: (1) A 是正定矩阵.(2) A 的所有顺序主子式全大于零. (3) A 的特征值全大于零. (4) 存在正定矩阵B ,使得A =B 2. (5) A 合同于E .(6) A 的一切主子式全大于零. (7) A 的一切主子矩阵都是正定矩阵.(8) 对任意的实列满秩矩阵C n×m ,都有C T AC 为正定矩阵.(9) 任意实可逆矩阵T ,都有T T AT 为正定矩阵. (10) 存在秩为n 的m ×n 实矩阵C 使A =C T C . (11) A =P T P ,P 是n 阶可逆矩阵.(12) A =R T R ,R 是n 阶主对角元素全大于零的上三角形矩阵.A =U T U ,U 是n 阶主对角元素全大于零的下三角形矩阵.[9]证明(1)⇒(2)设二次型f(x 1,x 2,⋯,x n )=∑∑a ij x i nj=1n i=1x j 是正定的.对于每个k ,1≤k ≤n ,设f k =(x 1,x 2,⋯,x n )=∑∑a ij x i n j=1n i=1x j以下证明f k 是一个正定二次型,对于任意一组不全为零的实数c 1,c 2,⋯,c k ,有)13(因此f k(x1,x2,⋯,x n)是正定的.由性质1可得,f k所对应的矩阵行列式|a11⋯a1n ⋮⋱⋮a k1⋯a kk|>0,k=1,2,⋯,n.从而A的各阶顺序主子式都大于零.(1)⇒(3)用反证法,若A的特征根λ1,λ2,⋯,λn不都大于零.不妨设λ1≤0,取A的属于λi的单位特征向量β≠0,就有βT Aβ=λ1≤0,这与A为正定矩阵相矛盾,所以A的特征值全大于零。
(3)⇒(1)设A的特征值为λ1,λ2,⋯,,,则存在正交矩阵T,使T T AT=diag(λ1,λ2,⋯,λn),A=Tdiag(λ1,λ2,⋯,λn)T T,任意取x≠0,则X T AX=X T Tdiag(λ1,λ2,⋯,λn)T T X=Y T diag(λ1,λ2,⋯,λn)Y,其中Y=X T T=(y1,y2,⋯,y n)≠0,于是X T AX=λ1y12+λ2y22+⋯+λn y n2> 0,即A为正定矩阵.(3)⇒(4)设A的特征值为λ1,λ2,⋯,λn,则存在正交矩阵T,使得T T AT=diag(λ1,λ2,⋯,λn).于是A=Tdiag(λ1,λ2,⋯,λn)T T=[Tdiag(√λ1,√λ2,⋯,√λn)T−1][Tdiag(√λ1,√λ2,⋯,√λn)T−1=B2,其中B=Tdiag(√λ1,√λ2,⋯,√λn)T−1,因为B为实对称矩阵且特征值√λi> 0(i=1,2,⋯,n),,故B为正定矩阵.(4)⇒(5)由A=B2=B T B=B T EB,,且B为正定矩阵,从而为可逆矩阵,也即A≅E.(5)⇒(1) 由A ≅E ,即存在可逆矩阵C ,使得A =C T EC =C T C..任意取X ≠0,令CX =Y =(y 1,y 2,⋯,y n )T ,则Y ≠0,于是X T AX =X T C T ACX =Y T Y =y 12+y 22+⋯+y n 2>0,故A 是正定矩阵.(1)⇒(6) 设A 是正定矩阵,且A k 为其任意k(1≤k ≤n)阶主子式对应的k 阶实对称矩阵(设取自A 的第i 1,i 2,⋯,i k 行和列),从而g =X k TA k X k (X k =(x i 1,⋯,x i k )为k 元实二次型.设c i 1,⋯,c i k 为任意一组不都为零的实数,则0,⋯,0,c i 1,0,⋯,0,c i k ,0,⋯,0.(即c i 1,⋯,c i k 以外的变量全取0)不全为0,又f =X T AX 是正定二次型,故g(c i 1,⋯,c i k )=f(0,⋯,0,c i 1,0,⋯,0,c i k ,0,⋯,0)>0.于是,g 是k 元正定的,A k 正定,从而A 的主子式|A k |>0.(6)⇒(7) 设A (i 1i 2i 3i 4i 1i 2i 3i 4)是n 阶正定矩阵A 的主子阵,则A (i 1i 2i 3i 4i 1i 2i 3i 4) 也是实对称矩阵,又因为A (i1i 2i 3i 4i 1i 2i 3i 4)的k 个顺序主子式均为A 的k 个主子式,由(6)知它们都大于零,从而为正定矩阵.(7)⇒(1) 因为A 本身也是它的一个主子矩阵,由(7)知A 为正定矩阵. (1)⇒(8) 任意取X ≠0,则有CX ≠0.(如果CX =0,则C T CX =0,而C T C 为m阶可逆矩阵,所以X =0与假设X ≠0矛盾)由于A 为正定矩阵,因此X T (C T AC )X =(CX )T A (CX )>0,从而C T AC 为正定矩阵.(8)⇒(9) 因为实可逆矩阵必然是实列满秩矩阵,由(8)知T T AT 为正定矩阵.(9)⇒(1) 设T T AT 为正定矩阵,其中T 为实可逆矩阵,又A =(T −1)T (T T AT)T −1, 由(9)知A 正定的.(10)⇒(1) 设存在秩为n 的m ×n 实矩阵C ,使得A =C T C .则A 为n 阶实对称方阵且对于任意不全为零的实数x 1,x 2,⋯,x n ,即X =(x 1,x 2,⋯,x n )T ≠0都会有CX =(y 1,⋯,y m )T ≠0,因此f(x 1,x 2,⋯,x n )=X T AX =(CX )T (CX )=y 12+⋯+y m 2>0,也就是说二次型f =X T AX 是正定的,于是A 为正定矩阵.(1)⇒(11) A 为正定矩阵A 合同于单位矩阵E 存在可逆矩阵P ,使得A =P T EP =P T P .(1)⇒(12) 设A 为正定矩阵,存在可逆矩阵P ,有A =P T P ,又P =QR (其中Q 为正交矩阵,R 为主对角元素全大于零的上三角形矩阵)从而A =P T P =R T Q T QR =R T R .(12)⇒(1) 若A =R T R ,其中R 为主对角元素全大于零的上三角形矩阵,当然R可逆,因此A 为正定矩阵.(1)⇒(13) 证明过程可类同于(1)⇔(12).1.2.2 半正定矩阵的等价命题定理2[2] 设A 是n 阶实对称矩阵,则下列叙述等价: (1) A 是半正定的.T ⇔⇔(2) A 的所有特征值非负. (3) A 的所有主子式非负. (4) A 的行列式大于等于零.(5) 存在n 阶实方阵B ,使得A =B T B . (6) A 主对角线的元素非负.(7) 存在秩为r 的n ×r 实矩阵B 使得A =B T B . (8)A 合同于主对角元素非负的对角矩阵.[2]证明 (1)⇒(3)设A 的k(1≤k ≤n)阶主子式对应的方阵为A k =(a i 1i 1⋯a i 1i k⋮⋱⋮a i k i 1⋯a i k i k),1≤i 1<⋯<i k ≤n. 又f(x 1,x 2,⋯,x n )=X T AX ,f 1=(x i 1,⋯,x i k )=X k TA k X k .若A 是半正定的,即f 是半正定的,从而f 1是半正定的.于是有k 阶实满秩方阵C k 使C k T A k C k =(λ1⋯⋯⋮⋱⋮⋯⋯λk), 其中λi ≥0.所以|C k T A k C k |=|A k ||C k |2=λ1,⋯,λk ≥0.但是|C |2>0,因此|A i |≥0.(3)⇒(1) 若A 的主子式非负,令B k =(a ij )为A 的k 阶顺序主子式所对应的方阵,则|λE k +B k |=|λ+a 11a 12⋯a 1k a21λ+a 22⋯a 2k⋮⋮⋮a k1a k2⋯λ+a kk|, 其中P i 为B k 中的一切i 阶主子式之和,故P i ≥0.从而当λ>0时,|λE k +B k |>0.即λE +A 的顺序主子式均为正,故λE +A 是正定矩阵.若不是半正定的,则必有X 0=(a 1,⋯,a n )≠0使X 0TAX 0=C <0,取λ=−c X 0T X 0=−ca 12+⋯+a n2>0,则X 0T (λE +A )X 0=λX 0T X 0=λX 0T X 0+X 0TAX 0=c −c =0.这与上面的结论矛盾.故A 必为半正定的.(1)⇒(5) 设A 是半正定的,且秩为r ,则存在实可逆方阵P ,使得A =P (Er00)P T 由上式得A =P (E r 000)(E r 00)P T=BB T , 其中B =P (E r 00). (5)⇒(1) 设A =BB T (B 为实方阵).则对任意x ≠0,有f =X T (BB T )X =(BX )T (BX)≥0.故f 是半正定的.因此,A 是半正定的.(1)⇒(7) 由等价条件(5)的证明过程中得A =P (Er00)P T. 令是的前列构成的秩为r (因为P 满秩)的n ×r 实矩阵,则A =(B ,C)(E r 000)(BC)=BB T . 1.3. 正定矩阵和半正定矩阵的性质1.3.1.正定矩阵的性质性质1[1] 正定矩阵的行列式大于零.证明 设A 是正定矩阵,则A 与I (单位矩阵)合同,于是存在可逆矩阵C 使A B P rA=C T IC=C T C.同时在两边取行列式,有|A|=|C T||C|=|C|2>0.性质2[3]A=(a ij)n×n是n阶正定矩,则A的主对角元全大于零.证明设εi是第i个分量为1而其余分量为零的n维列向量,由于A正定,则a ii=εi T Aεi>0.i=1,2,⋯,n性质3[3]正定矩阵A=(a ij)中绝对值最大的元素必在它的主对角线上.证明鉴于A正定,故A的任一二阶主子式大于零,即|a ii a ija ij a jj|=a ii a jj−a ij2>0,即|a ij|<(a ii a jj)12(i≠j,i,j=1,2,⋯,n).设A的主对角线上最大元素为a kk,则a kk>0.|a ij|<(a ii a jj)12≤(a kk2)12=a kk(i≠j,i,j=1,2,⋯,n).即A中绝对值最大的元素必在它的主对角线上.性质4[5]若A是正定矩阵,则是kA,A+kB正定矩阵,其中k>0.性质5[6]若A是正定矩阵,则A−1,A∗是正定矩阵.证明先证A−1是正定的.因为A是正定矩阵,所以A也是对称矩阵,即A T=A,因为(A−1)T=(A T)−1,从而有(A−1)T=A−1,即A−1是对称矩阵.因为矩阵A是正定的,则对于任何非零的n维实列向量X,恒有X T AX>0,由于X T A−1X=X T A−1AA−1X>0,且A−1X不等于零;故(A−1X)T=AA−1X>0,从而X T A−1X>0,则A−1是正定矩阵.再证明A∗也是正定矩阵.由于A是正定矩阵,且A∗=|A|A−1则(A∗)T=A∗,A∗是实对称矩阵.由于A∗的特征值为|A|λ1,⋯,|A|λn>0.所以A∗是正定矩阵.性质6[1]若A是正定矩阵,则对于任意整数k,A k都是正定矩阵.证明(1)当k=0时,A k=I,A k是正定矩阵.(2)当k<0时,k=−|k|,即A k=(A−1)|k|,根据正定矩阵的性质5可知,A−1也是正定的,故下面的k只要是正整数就行.当k为偶数时,由于A T=A,且A k=(A k2)T(A k2),由正定矩阵的条件(6)可知,是正定矩阵.当k为奇数时,由于...A k=A k−12AAk−12=Ak−12C T CAk−12=(CAk−12)T(CAk−12).因此A k是正定矩阵.性质7[7]A是n阶正定矩阵,则0<|A|<a11a22⋯a nn,其中a ii(i=1,⋯,n)是A的主对角元素.证明由于A正定,则A的特征根均大于零,而|A|=λ1λ2⋯λn,其中λi(i=1,⋯,n)是A的特征根,故|A|>0.设A=(A n−1a a T a nn)式中A n−1=(a11a1,n−1an−1,1an−1,n−1),a=(a1n,a2n,⋯,a n−1,n)T.因A是正定矩阵,故A n−1是正定的,且A n−1是可逆阵,对A进行广义初等变换,得(A n−1aa T a nn)→(A n−1aa a nn−a T A n−1−1a)→(A n−100a T A n−1−1a),于是k A(I n−10−a T A n−1−11)(A n−1a a a nn )(I n−1−A n−1−1a 01)=(A n−1a 0a nn −a T A n−1−1a ), 等式两边同取行列式得,|A |=|A n−1|(a nn −a T A n−1−1a ).由A n−1正定,得A n−1正定,那么a T A n−1−1a ≥0,|A n−1|≥0.因此|A |≤|A n−1|a nn ,同理可得,|A n−1|≤|A n−2|a n−1,n−1其中A n−1是A 的n −2阶顺序主子式,一直这样下去,可得|A |≤|A n−1|a nn ≤|A n−2|a n−1,n−1a nn ≤⋯≤a 11a 22⋯a nn .性质8[5] 若A 是正定矩阵,m 是任意正整数,则存在正定矩阵B ,使得A =B m . 证明 A 是正定矩阵,必存在正交矩阵Q ,使Q T AQ =(λ1⋯0⋮⋱⋮0⋯λn ). 式中λ1,⋯,λn >0,于是A =Q (λ1⋯0⋮⋱⋮0⋯λn )Q T , 令B =Q (√λ1n ⋯0⋮⋱⋮0⋯√λnn )Q T . 则A =B m ,结论得证.性质9[7] 设A ,B 都是n 阶正定矩阵,则AB 是正定矩阵的充要条件是AB =BA.. 证明 ⇒)AB 正定,下证AB =BA...因为AB 正定,则AB 必为实对称的,故由A ,B 正定得AB =(AB)T =B T A T =BA ,⇐)AB=BA,下证AB正定.AB是实对称矩阵.又A,B都正定,必存在实满秩方阵P,Q使得A=P T P,B=Q T Q,这时AB=P T PQ T Q与矩阵QP T PQ T=Q(P T PQ T Q)Q−1=Q(AB)Q−1相似,于是相同的特征根.但QP T PQ T=(PQ T)T PQ T为正定矩阵,所以它的特征根全是正实数,故AB的特征根也都是正实数,从而AB是正定的.性质10 [7]A为n阶实方阵,A正定的充要条件是存在秩为n的m×n实方阵C,使得A=C T C..证明⇒)若A正定,则根据正定矩阵的等价条件,结果显然成立.⇐)设存在秩为n的m×n实方阵C,使得A=C T C,则A为实对称方阵,并且由齐次线性方程组的理论知识可知,对任意不全为零的实数x1,⋯,x n,即X=(x1,⋯,x n)≠0,,必有CX=(y1,⋯,y n)≠0,,且f(x1,⋯,x n)=X T AX=(CX)T(CX)=y12+⋯+y n2>0.即是f正定二次型,因而矩阵A正定的.性质11[7]若A,B都是同阶正定矩阵,则AB的特征值全大于零.证明由于A是正定矩阵,则存在n阶可逆矩阵P,使得A=P T P,则AB=P T PB=P T(PBP T)(P T)T,这表明AB相似于P T BP,又(PBP T)T=PB T P T=PBP T,即PBP T,是实对称矩阵,它的特征根全部为实数,因为相似矩阵具有相同的特征根,则AB的特征根也为实数. 设λ为AB的任一特征根,那么λ必为实数α为AB的属于λ的特征向量,即(AB)α=λα,α≠0.Bα=λA−1α,且αT Bα=λαT A−1α.则因为A正定,所以A−1正定.又B正定,故αT Bα>0,αT A−1α>0从而λ>0.1.3.2.半正定矩阵的性质性质1[1]若矩阵A是半正定矩阵,则对于任何非负实数k,kA也是半正定矩阵.证明设f(x1,⋯,x n)=X T AX,因为A是半正定的,所以对任意的n维实向量X≠0,都有X T AX>0,由已知X T kAX=k(X T AX)≥0,,则有k>0,故kA也是半正定的.性质2[8]可逆半正定矩阵的逆矩阵是半正定矩阵.性质3[7]设A=a ij,B=b ij都是n阶方阵.定义A°B=(a ij b ij)若A,B都是半正定的,则A°B也是半正定的.证明因为B=b ij是半正定矩阵,故存在实方阵P=(p ij)使B=P T P,即有.b ij=p i1p j1+⋯+p in p jn,由于A,B都为实对称的,故A°B也是实对称的,且X T(A°B)X=∑a ij b iji,j x i x j=∑a iji,j(∑p ik p jki,j)x i x j=∑∑a ij(i,ji,jp ik x i)(p jk x j)=∑Y k T AY kk,(1)其中X=(x1,⋯,x n),Y k=(p1k x1⋯p nk x n).因为A是半正定的,故对于任何实向量X,都有Y k T AY k.再由(1)得X T(A°B)X>0,故A°B是半正定矩阵.性质4[7] 若A,B都是半正定矩阵,则A+B也是半正定的,而且|A+B|≥|A|+ |B|.性质5[8] 若A是半正定矩阵,则对任意k=2m,m∈Z∗,有A k也为半正定的.证明因为A为半正定矩阵,且k=2m,m是正整数,那么A k=A2m=A m A m=(A m)T(A m),而A m为实矩阵,则由定理2可知A k为半正定矩阵.第2章 正定矩阵和半正定矩阵的判定方法2.1. 定义法例1 设A 是m ×n 阶正定矩阵,B 是m ×n 实矩阵.证明:B T AB 正定↔r (B )=n..证明 ⇒)反证法.由于A ,B T AB 分别是m 阶,n 阶正定阵.反设r(B)<n ,则齐次线性方程组BX =0有非零解α,即 Bα=0 ,于是.αT (B T AB )α=(Bα)T A (Bα)=0,这与B T AB 为正定矩阵矛盾,故r (B )=n.⇐)由于A 正定,则A 为实对称矩阵,(B T AB)T =B T A T B =B T AB ,即B T AB 是实对称矩阵.下证B T AB 是正定矩阵.由于r (B )=n ,故齐次线性方程组BX =0只有零解,于是对任意n 维实列向量α≠0,有Bα≠0,.又A 为正定矩阵,则(Bα)T A (Bα)>0,即αT (B T AB )α>0,故B T AB 为正定矩阵.2.2. 主子式法例2 判断矩阵A 是否正定,其中A =(303120031)解 由于A 的一阶顺序主子式为|3|=3>0,A 的二阶顺序主子式为|3012|=6−0=6>0, A 的三阶顺序主子式为|303120031|=6+9=15>0, 因此矩阵A 是正定矩阵..例3 t 为何值时矩阵A 为正定矩阵?求出对应的t 值.其中A =|21011t 0t 20|, 解 因为A 正定的充要条件是A 的各阶顺序主子式均大于零,|2|=2>0,|2111|=2−1=1>0,|21011t 20t 21|>0, 解得−√2<t <√2.因此,当−√2<t <√2时矩阵A 为正定矩阵.2.3. 特征值法例4 判断矩阵A 是否正定,其中A =|50003−20−23| 解 矩阵A 的特征多项式为f A (λ)=|λ−5000λ−3202λ−3|=(λ−5)2(λ−1),所以A 的特征根分别为λ1=λ2=5>0,λ3=1>0.因此,矩阵A 是正定矩阵.例5 设A ,C ,为n 阶正定矩阵,B 是矩阵方程AX +XA =C 的唯一解.证明:(1) B 是实对称矩阵;(2) B 是正定矩阵.证明 (1)因为A ,C 为n 阶正定矩阵,A ,C 必为实对称矩阵,故A T =A ,C T =C ,对AB +BA =C 两边取共轭,得AB ̅+B ̅A =C ,.2''''''αλααλααλαααααααA A A BA AB C =+=+=从而B ̅也是AX +XA =C 的解,又由于方程的解是唯一的,故B =B ̅,即B 是实矩阵.对AB +BA =C 两边取转置,得C =C T =(AB +BA)T =B T A +AB T ,这表明B 也是AX +XA =C 的解,由解的唯一性可知B =B T ,即B 是实对称矩阵.设λ是B 的任意特征根,α为B 的属于λ的特征向量,则Bα=λα,于是 由于A ,C 为正定矩阵,故αT Aα>0,αT Cα>0,则B 的特征根λ大于0,则B 是正定矩阵.例6 设A 是n 阶正定矩阵,B 是n 阶半正定矩阵,并且满足A 2=B 2.证明: B 是正定矩阵.证明 因为B 是半正定的,必存在正交矩阵T ,使得B =T −1(λ10000λ20000⋱0000λn), 式中λi ≥0(i =1,2,⋯,n),对照题设条件可知A 2=B 2=T −1( λ120000λ220000⋱0000λn 2)T , 因A 2是正定矩阵,于是λI 2>0,从而λI >0,因此B 为正定矩阵.2.4. 与单位矩阵E 合同法例7 A 是n 阶可逆矩阵,证明:A T A 是正定的.)2(证明 因为(A T A)T =(A T A),则矩阵A T A 对称的.又A T A =A T IA ,同时A 是可逆矩阵,于是A T A ≅I ,则矩阵A T A 是正定的.例8 证明:分块矩阵Q =(A 00B )是正定矩阵,其中A ,B 分别为m ,n 阶正定矩阵.证明 因为矩阵A ,B 为正定矩阵,则必存在可逆矩阵C m×m 和D n×n ,使c T AC =I m ,D T BD =I n ,设P =(C 00D ),那么P T =(C T 00D T),且矩阵P 是m +n 阶可逆的. p T QP =(C T 00D T )(A 00B )(C 00D )=(C T AC 00D T BD )=(I m 00I n ), 所以Q ≅I ,故分块矩阵Q =(A 00B)是正定矩阵.第3章定矩阵和半正定矩阵的应用3.1.在不等式问题中的应用例9设A,B,C是一个三角形的内角.证明:对任意实数x,y,z,都有x2+y2+z2≥2xy cos A+2xz cos B+2yz cos C.证明f(x,y,z)=x2+y2+z2−2xy cos A−2xz cos B−2yz cos C.=(x,y,z)(1−cos A−cos B −cos A1−cos C−cos B−cos C1)(xyz),令上式中的三阶实对称矩阵为T,由于T的一阶顺序主子式为|1|=1>0,T的二阶顺序主子式为T的三阶顺序主子式为|1−cos A−cos B −cos A1−cos C−cos B−cos C1|=1−2cos A cos B cos B−(cos A)2−(cos B)2−(cos C)2=0.故f(x,y,z)半正定,即f(x,y,z)≥0,从而结论成立.例10设∆ABC的三条边分别为a,b,c面积为S.证明:a2+b2+c2≥4√3S.证明由于c2=a2+b2−2ab cos C,s=1ab sin C.所以利用上述公式,可以将问题转化为f(a,b)=a2+b2+a2+b2−2ab cos C−2√3ab sin C=2a2+2b2−2ab(cos C+√3)sin C=2a2+2b2−4ab sin (π6+C).二次型f(a,b)的矩阵为A=(2−2sin (π6+C)−2sin (π6+C)2),由于A的一阶、二阶主子式满足2>0,|A|=4[1−(sin (π6+C)]2≥0,故A半正定,即二次型f(a,b)半正定,从而结论成立.3.2.在多元函数极值问题中的应用定义7[4]设n元实函数f(X1,⋯,X n)在X=(x1,⋯,x n)∈R n的某个领域内存在一阶、二阶连续偏导数.记Vf(x)=(ð(X)ðX1,ð(X)ðX2,⋯,ð(X)ðX n),Vf(x)称为函数f(X)在点X=(x1,⋯,x n)T处的梯度.定义8[4]H(X)=(∂2f(X)ðX i X j)m×n=(∂2f(X)ðX i2∂2f(X)ðX1ðX2⋯∂2f(X)ðX1ðn⋮⋮⋯⋮∂2f(X)ðX nðX1∂2f(X)ðX nðX2⋯∂2f(X)ðX n2)此矩阵称为函数在处的(Hessian)黑塞矩阵,则H(X)是由f(X)的n2个二阶偏导数构成的n阶实对称矩阵.定义9[4](极值充分条件)设f(X)在点X0∈R n的某个领域内存在一阶、二阶连续偏导数时,且Vf(X0)=(ð(X0)ðX1,ð(X0)ðX2,⋯,ð(X0)ðX n(1)当H(X0)是正定矩阵时,f(X)在点X0处取得极小值;(2)当H(X0)是负定矩阵时,f(X)在点X0处取得极大值;(3)当H(X0)是不定矩阵时,f(X)在点X0处不取极值;定义10[4](极值必要条件)设函数f(X)在点X0=(x10,⋯,x n0)T处可微,且X0为该函数的极值点,则1)X0必为f(X)的稳定点,即Vf(X0).2)若f(X)在X0的某领域U(X0)存在连续二阶偏导数,则当f(X0)为极小值时,f(X)在X0的黑塞矩阵为正定或半正定;则当f(X0)为极大值时,f(X)在X0的黑塞矩阵为负定或半负定.例1 求多元函数f(x,y,z)=x2+y2+z2+2x+y−z的极值.解先求f的驻点,由{f x=2x+2=0f y=2y+1=0 f z=2z−1=0,解得x1=−1,x2=−12,x3=12.求得函数驻点为P0=(−1,−12,12).再求(Hessian)黑塞矩阵,因为f xx=2,f xy=0,f xz=0,f yx=0,f yy=2,f yz=0,f z x=0,f z y=0,f z z=2.故H =(200020002).由于H 的各阶顺序主子式全为零,故矩阵H 是正定的.所以P 0=(−1,−12,12)是f(x ,y ,z)的极小值点.且在点P 0=(−1,−12,12) 的极小值为f =(−1,−12,12)=−32.例2 求多元函数f (x )=x 12+2x 1x 2+2x 22+x 32−4x 2x 3+3x 1x 3的黑塞矩阵,并根据结果判断该函数的极值点.解 先求驻点,由{f x 1=2x 1+2x 2+3x 3=0f x 2=4x 2+2x 1−4x 3=0f x 3=2x 3−4x 2+3x 1=0,解得x 1=0,x 2=0,x 3=0.可得驻点为P 0=(0,0,0).再求(Hessian)黑塞矩阵,因为fx 1x 1=2,fx 1x 2=2,fx 1x 3=3,fx 2x 1=2,fx 2x 2=4,fx 1x 3=−4,fx 3x 1=3,fx 3x 2=−4,fx 3x 3=2,故所求得(Hessian)黑塞矩阵为H =(22324−43−42),由于H 11=2>0,H 22=|2224|=4>0,H 33=(22324−43−42)<0,所以黑塞矩阵为不定矩阵,故P 0不是极值点.参考文献[1]王萼芳,石生明.高等代数[M].第四版.北京:高等教育出版社,2013.[2]姜国.正定矩阵的判定及性质[J].湖北师范学院学报(自然科学版),2006,26(1): 97-99.[3]岳贵鑫.正定矩阵及其应用[J].辽宁省交通高等专科学院学报,2008,10(5),31-33.[4]华东师范大学数学系.数学分析(第四版)[M].高等教育出版社.[5]张丹,刘庆平.正定矩阵的性质及相关问题[J].中南大学学报,2011,31(4).[6]黄云美.正定矩阵的性质及其应用[J].烟台职业学院学报,2011,17(3):85~88.[7]陈大新.矩阵理论[M].上海:上海交通大学出版社.[8]刘畅.正定矩阵性质的推广[J].沈阳师范大学学报,2009,27(3),268~271.[9]何亚丽.线性代数[M].科学出版社.[7]杨子胥.高等代数[M].高等教育出版社.[8]张淑娜.半正定二次型及半正定矩阵[J].通化师范学院学报,2003.24(6).[9]史秀英.正定矩阵的等价命题及其应用[J].赤峰教育学院学报,2000(2).致谢光阴似箭,日月如棱。