二次型 & 有定性
二次型
以三元二次型举例:x2+x1y1+y2−2z2(有几个变量就是几元,所有项必须是二次)
- x2,y2,z2叫平方项
- xy叫交叉项
x2+xy2、x2+y−6、x2+8不是二次型
二次型的矩阵一定是对称的,经常会用到:AT=A
二次型->矩阵表达式
x12+2x1x2+x22−x2x3+2x32−2x1x3
- 平方项的系数直接做成主对角线元素
- 交叉项的系数除以2放到2个对称的相应位置上
系数位置如下:
⎝⎜⎜⎜⎜⎛x1221×2x1x221×−2x1x321×2x1x2x2221×−x2x321×−2x1x321×−x2x32x32⎠⎟⎟⎟⎟⎞
xTAx如下(其中A叫作二次型矩阵):
(x1,x2,x3)×⎝⎜⎜⎛11−111−21−1−212⎠⎟⎟⎞×⎝⎛x1x2x3⎠⎞
矩阵表达式->二次型
注意:矩阵一定是对称的,如果不是对称的,需要转换下(2x1x3+3x3x1=5x1x3):
⎝⎛11310121−1⎠⎞−>⎝⎜⎜⎛1125101251−1⎠⎟⎟⎞
- 主对接线元素直接做为平方项的系数
- 取主对角线右上角元素乘以2作为交叉项的系数
x12−x32+2x1x2+5x1x3+2x2x3
标准形 & 线性替换
定义:只有平方项的二次型叫标准形,d1y12+d2y22+...+dnyn2如:x12−3x32
线性替换:f(X)=XTAX,X=CY−>(CY)TACY=YTCTACY=YT(CTAC)Y
- 如果B=CTAC= ∧,那么引入CY后的二次型是标准形
- 如果|C| != 0,那么这个线性替换叫作:可逆替换(或非退化替换或满秩替换或非奇异替换)
- 如果|C| = 0,那么这个线性替换叫作:退化替换
两次替换(X=C1C2Z):
证明:X=C1Y,Y=C2Z
XTAX=(C1Y)TAC1Y=YT(C1TAC1)Y=(C2Z)TC1TAC1C2Z
=ZT(C2TC1TAC1C2)Z=ZT(C1C2)TA(C1C2)Z
定理:二次型经过线性替换后仍然是一个二次型
证明:f(X)=XTAX=YT(CTAC)Y
验证CTAC是对称矩阵,设B=CTAC
BT=(CTAC)T=CTAT(CT)T=CTAC=B(A是对称矩阵,所以AT=A)
BT=B,所以B=CTAC是对称的
合同
合同:A、B是n阶方阵,如果存在可逆矩阵C,使得CTAC=B,那么A和B合同。记作:A≃B
- 反身性:A和A合同,A≃A(ETAE=A)
- 对称性:如果A和B合同,那么B和A也合同(证明:同时左乘(CT)−1和右乘C−1)
- 传递性:如果A、B合同,B、C合同,那么A、C也合同(证明:PTAP=B,QTBQ=C,QTPTAPQ=C,(PQ)TA(PQ)=C)
性质1:如果A、B合同,那么A、B的秩一定相同
证明:CTAC=B,A左乘一个可逆矩阵,右乘一个可逆矩阵,矩阵的秩不变
性质2:如果A、B合同,那么A对称的充要条件是B也是对称的(AT=A <=> BT=B)=>
证明:CTAC=B,B=BT
CTAC=(CTAC)T=CTATC
同时左乘(CT)−1,右乘C−1
(CT)−1CTACC−1=(CT)−1CTATCC−1
得:A=AT
性质3:A、B合同,如果A、B均可逆,那么A−1≃B−1
证明:CTAC=B
因为两边都可逆,所以:B−1=(CTAC)−1=C−1A−1(CT)−1=C−1A−1(C−1)T
得:[(C−1)T]TA−1(C−1)T
性质4:如果A、B合同,那么AT≃BT
证明:CTAC=B,BT=(CTAC)T=CTATC
化二次型化为标准形
- 配方法
- 初等变换法
- 正交替换法(就是实对称矩阵对角化,因为正交矩阵肯定是对称的)
配方法
x12−3x22+4x32−2x1x2+2x1x3−6x2x3
注意:先x1再x2再x3...,当x1配完后,后面的式子不能再出现x1,以此类推...
=x12−2x1(x2−x3)+(x2−x3)2−(x2−x3)2−3x22+4x32−6x2x3
=(x1−x2+x3)2−4x22−4x2x3+3x32
=(x1−x2+x3)2−(4x22+4x2x3+x32)+x32+3x32
=(x1−x2+x3)2−(2x2+x3)2+4x32
=y12−y22+4y32
⎩⎪⎨⎪⎧y1=x1−x2+x3y2=2x2+x3y3=x3
写出线性替换(X=CY):
⎩⎪⎨⎪⎧x1=22y1+y2+y3x2=2y2−y3x3=y3
配方法(只有交叉项)
2x1x2−4x1x3+10x2x3
直接代入下式:
⎩⎪⎨⎪⎧x1=y1−y2x2=y1+y2x3=y3
注意:如果有x4写成:x4=y4、x5写成x5=y5...
得:2y12−2y22+6y1y3+14y2y3
再用上面的方式正常配方就可以了
初等变换法
f(X)=XTAX,X=CY,CTC= ∧,①C=?,②∧=?
因为C可逆,C=P1P2...Ps
(P1P2...Ps)TAP1P2...Ps= ∧
PsTP2T...P1TAP1P2...Ps= ∧
左乘一个单位阵,EP1P2...Ps=C
矩阵左乘一些初等矩阵,相当于做初等行变换
矩阵右乘一些初等矩阵,相当于做初等列变换
- 对A和E做同样的初等列变换
- 只对A做相应的初等行变换(相应的举例如下)
- 交换1、3列,交换1、3行
- 2乘以第3列,2乘以第3行
- 2乘以第2列加到第3列,2乘以第2行加到第3行
- A化成对角阵之时,E化成的就是C
注意:列变换和行变换配套做,列行、列行、列行...
计算题(初等变换法)
A=⎝⎛111122121⎠⎞
解:(AE)=⎝⎜⎜⎜⎜⎜⎜⎛111100122010121001⎠⎟⎟⎟⎟⎟⎟⎞
第1列乘以-1加到第2列(整体):⎝⎜⎜⎜⎜⎜⎜⎛111100011−110121001⎠⎟⎟⎟⎟⎟⎟⎞
第1行乘以-1加到第2行(只对A):⎝⎜⎜⎜⎜⎜⎜⎛101100011−110111001⎠⎟⎟⎟⎟⎟⎟⎞
...
最终求得:⎝⎜⎜⎜⎜⎜⎜⎛100100010−11000−10−11⎠⎟⎟⎟⎟⎟⎟⎞
∧ = ⎝⎛10001000−1⎠⎞
C = ⎝⎛100−1100−11⎠⎞
步骤省略:多套列行变换可以一起做
检查方法:成套变换后,矩阵还是对称的
规范形
整理对角形元素依次为:1、1、...、1、-1、-1、...、-1、0、0...0
y12+...+yp2−yp+12−...−yr2
- x12+x22+x32 ✅
- x12−x22+x32 ❌
- x12+x22+x42 ❌
任意二次型都能化成规范形
- 正项数:正惯性指数(正项的个数)
- 负项数:负惯性指数(负项的个数)
- 符号差:正项数-负项数
化成规范型后,1和-1的个数,就是二次型的秩,等于矩阵的秩
定理:任意矩阵A都与规范形的∧合同
定理:合同的充要条件,有相同的秩、相同的正项数、相同的负项数
正交替换
参照:对角化 -> 计算步骤总结
有定性
有定性:二次型f(X)=XTAX,X!=0(当X=0,结果肯定等于0,所以没有讨论意义),存在以下4种情况:
- 正定:XTAX>0,举例:f(x1,x2,x3)=x12+2x22+x32,f > 0
- 半正定:XTAX>=0,举例:f(x1,x2,x3)=x12+2x22+0x32,f >= 0(x取(0,0,1)T往里带)
- 负定:XTAX<0,举例:f(x1,x2,x3)=−x12−2x22−x32,f < 0
- 半负定:XTAX<=0,举例:f(x1,x2,x3)=−x12−2x22+0x32,f <= 0(x取(0,0,1)T往里带)
不定:举例:f(x1,x2,x3)=x12−x22+x32,代入(1,0,0)T,结果大于0。代入(0,1,0)T,结果小于0
定理:f(x1,x2,...,xn)=a1x12+a2x22+...+anxn2 正定 <=> ai>0 =>
先证必要性:代入(1,0,0...,0),f(1,0,0,...,0)=a1>0
依次一直取到,f(0,0,0,...,1)=an>0
再证充分性:因为ai>0,a1x12+a2x22+...+anxn2>=0(平方和大于等于0)
又因为,X !=0,所以a1x12+a2x22+...+anxn2>0
定理:正定二次型经过线性替换(非退化)后仍是正定的
证明:设:f(X)=XTAX是正定的,X=CY
因为非退化,所以C可逆
f(X)=XTAX=(CY)TACY=YTCTACY
任给Y != 0,只要证明YTCTACY>0,就是正定
因为X=CY,只要Y不等于0,X就不等于0,所以f(X)>0,所以正定
定理:f(X)=XTAX正定 <=> 标准形是 d1y12+d2y22+...+dnyn2(di>0)
=>
定理:f(X)=XTAX正定 <=> 正惯性指数是n
=>
定理:f(X)=XTAX正定 <=> A≃E(单位阵是特殊的规范形)
=>
定理:f(X)=XTAX正定 <=> n个特征值(所有)均大于0
=>
定理:f(X)=XTAX正定 => |A| > 0
证明:CTAC=E(A和单位阵合同)
同时左乘(CT)−1,右乘C−1
A=(CT)−1C−1=(C−1)T(C−1)
把C−1写成P,PTP
∣A∣=∣PT∣×∣P∣=∣P∣2
因为P可逆,所以P != 0
所以∣P∣2>0
顺序主子式
定理:f(X)=XTAX正定 <=> 各阶顺序主子式 > 0=>
判断 f(x1,x2,x3)=2x12+2x1x2+x22+x32 是正定的?
A=⎝⎛210110001⎠⎞
1阶顺序主子式(a11):|2| > 0
2阶顺序主子式(a11到a22):∣∣∣∣2111∣∣∣∣ > 0
3阶顺序主子式:|A| > 0
所以f(x1,x2,x3)=2x12+2x1x2+x22+x32 是正定的
如果有4阶、5阶,依次类推
已知,A=⎝⎛11112010t⎠⎞是正定的,求t的取值范围?
1阶顺序主子式:|1| = 1 > 0
2阶顺序主子式:∣∣∣∣1112∣∣∣∣=2−1>0
3阶顺序主子式:∣∣∣∣∣∣11112010t∣∣∣∣∣∣=∣∣∣∣∣∣111−10−210t∣∣∣∣∣∣
=(−1)1+2∣∣∣∣−1−21t∣∣∣∣=−(−t+2)>0
所以,t>2
性质1:如果A是正定,那么A−1也正定
证明:A的特征值λi>0,A−1的特征值λi1>0
性质2:如果A是正定,那么A∗也正定
证明:A的特征值λi>0,A−1=∣A∣1A∗,所以A∗=∣A∣A−1
所以A∗的特征值是λi∣A∣,因为|A|大于0,所以A∗的特征值也大于0
性质3:如果A是正定,那么Ak也正定
证明:因为λi>0,所以λik>0
性质4:A正定,B正定(或半正定),那么A+B正定
证明:任给X != 0,XT(A+B)X=XTAX+XTBX,左半部分大于0,右半部分大于等于0,所以A+B>0
性质5:A正定,A的主对角线元素全大于0,即:aii>0
证明:因为f(X) > 0,X != 0
(x1x2...xn)⎝⎜⎜⎛a11a21⋮an1a12a22⋮an2⋯⋯⋱⋯a1na2n⋮ann⎠⎟⎟⎞⎝⎜⎜⎛x1x2⋮xn⎠⎟⎟⎞
分别代入(1,0,0...,0)、(0,1,0...,0)、...、(0,0,0...,1),求得aii>0