LA@二次型分类@正定二次型@主子式
文章目录
abstract
- 介绍正定二次型相关概念和性质
- 二次型的分类
- 主子式
引言
- 科学技术上用的较多的二次型是正(或负)惯性指数为 n n n的 n n n元二次型,这类二次型是正定二次型或负定二次型
正定二次型
- 设二次形 f ( x ) = x T A x f(\bold{x})=\bold{x^{T}\bold{A}x} f(x)=xTAx,其中 x = ( x 1 , x 2 , ⋯ , x n ) T \bold{x}=(x_1,x_2,\cdots,x_n)^{T} x=(x1,x2,⋯,xn)T,如果对于任意的 x ≠ 0 \bold{x\neq{0}} x=0都有 f ( x ) > 0 f(\bold{x})>0 f(x)>0,称 f f f为正定二次型
- 其中 A \bold{\bold{A}} A是正定矩阵,显然: f ( 0 ) = 0 T A 0 = 0 f(\bold{0})=\bold{0}^{T}\bold{\bold{A}0}=0 f(0)=0TA0=0当且仅当 x = 0 \bold{x=0} x=0
小结
- 虽然在定义上并没有直接说明正定二次型的标准形项数为 n n n且都同号,但是通过推理可以得出此结论
可逆线性变换不改变二次型的正定性
- 若 f ( x ) = x T A x f(\bold x)=\bold{x^{T}\bold{A}x} f(x)=xTAx是正定的,经过可逆变换 x = C y \bold{x=Cy} x=Cy的 g ( y ) = y T B y g(\bold y)=\bold{y^{T}By} g(y)=yTBy,(其中 B = C T A C \bold{B=C^{T}\bold{A}C} B=CTAC)也是正定的
- 类似的,若 f f f是不正定的,则 g g g也是不正定的
- 用矩阵描述:若 A = C T B C \bold{\bold{A}=C^{T}BC} A=CTBC,且 C \bold{C} C可逆,则 A , B \bold{\bold{A},B} A,B有相同的正定性
- 证明:
-
对于可逆矩阵 C \bold C C,和任意非零向量 z = ( z 1 , ⋯ , z n ) T ≠ 0 \bold z=(z_1,\cdots,z_n)^{T}\neq{0} z=(z1,⋯,zn)T=0,有 C z ≠ 0 \bold{Cz\neq{0}} Cz=0
- 因为 C \bold C C可逆,齐次线性方程 C z = 0 \bold{Cz=0} Cz=0只有零解
- 从而 z ≠ 0 \bold{z\neq{0}} z=0,一定由 C z ≠ 0 \bold{Cz\neq{0}} Cz=0
-
设可逆线性变换 x = C y \bold{x=Cy} x=Cy将 f f f线性变换为 g g g, f ( x ) = x T A x = y T B y = g ( y ) f(\bold x)=\bold{x^{T}\bold{A}x}=\bold{y^{T}By}=g(\bold y) f(x)=xTAx=yTBy=g(y),其中, x , y \bold {x,y} x,y为 n n n维列向量, B = C T A C \bold{B=C^{T}\bold{A}C} B=CTAC
-
为证明 g ( y ) g(\bold y) g(y)是正定的,就是要证明 ∀ z ≠ 0 ⇒ g ( z ) > 0 \forall \bold{z\neq{0}}\Rightarrow{g(\bold z)>0} ∀z=0⇒g(z)>0
- g ( z ) = z T B z = z T C T A C z g(\bold z)=\bold{z^{T}Bz}=\bold{z^{T}C^{T}\bold{A}Cz} g(z)=zTBz=zTCTACz= ( C z ) T A ( C z ) \bold{(Cz)^{T}\bold{A}(Cz)} (Cz)TA(Cz)
- 记 r = C z \bold {r=Cz} r=Cz,前面已经讨论过 C z ≠ 0 \bold {Cz\neq{0}} Cz=0,从而列向量 r ≠ 0 \bold {r\neq{0}} r=0
- g ( z ) = r T A r = f ( r ) > 0 g(\bold z)=\bold {r^{T}\bold{A}r}=f(\bold r)>0 g(z)=rTAr=f(r)>0(由 f f f的正定性)
- 从而 g ( z ) > 0 g(\bold z)>0 g(z)>0,二次型 g g g依然是正定的
-
若 f f f不是正定的,则需证明 ∃ z \exist{\bold{z}} ∃z s.t. g ( z ) ⩽ 0 g(\bold{z})\leqslant{0} g(z)⩽0
- 设非零向量 x 0 \bold x_0 x0设 f ( x 0 ) ⩽ 0 f(\bold{x_0})\leqslant{0} f(x0)⩽0,令 x 0 = C z 0 \bold{x_0=Cz_0} x0=Cz0,则 z 0 = C − 1 x 0 ≠ 0 \bold{z_0=C^{-1}x_0\neq{0}} z0=C−1x0=0,则 f ( C z 0 ) f(\bold{Cz_0}) f(Cz0)= g ( z 0 ) ⩽ 0 g(\bold{z}_0)\leqslant{0} g(z0)⩽0
- 说明 g ( z ) g(\bold{z}) g(z)不满足正定条件, g ( z ) g(\bold{z}) g(z)是非正定的
-
二次型是正定的充要条件
-
n
n
n元实二次型
f
(
x
)
=
x
T
A
x
f(\bold{x})=\bold{x^{T}\bold{A}x}
f(x)=xTAx是正定的当且仅当
f
f
f的正惯性指数为
n
n
n
- 若二次型 f f f的正惯性指数为 n n n,则 f f f是正定二次型
- 若 f f f是正定二次型,则二次型 f f f的正惯性指数为 n n n
- 正惯性指数为
n
n
n的等价描述:
- 规范形系数全为1
- 标准形所有系数为正数
- 证明:
- 设可逆线性变换 x = C y \bold{x=Cy} x=Cy将二次型 f ( x ) = x T A x f(\bold{x})=\bold{x^{T}\bold{A}x} f(x)=xTAx化为标准形 f ( x ) = f ( C y ) f(\bold{x})=f(\bold{Cy}) f(x)=f(Cy)= g ( y ) = y T Λ y g(\bold{y})=\bold{y^{T}\Lambda{y}} g(y)=yTΛy= ∑ i = 1 n k i y i 2 \sum_{i=1}^{n}k_iy_i^2 ∑i=1nkiyi2, Λ = C T A C \bold{\Lambda=C^{T}\bold{A}C} Λ=CTAC= diag ( k 1 , ⋯ , k n ) \text{diag}(k_1,\cdots,k_n) diag(k1,⋯,kn)
- 充分性
- 设标准形系数 k i > 0 , i = 1 , 2 , ⋯ , n k_i>0,i=1,2,\cdots,n ki>0,i=1,2,⋯,n,则任意 x ≠ 0 \bold{x\neq{0}} x=0,满足 f ( x ) > 0 f(\bold{x})>{0} f(x)>0,从而 f f f就是正定的
- 必要性
- 用反正法证明,在 f f f是正定并且标准形存在负系数的情况下,找到一个非零向量 x 0 \bold{x_0} x0能使 f ( x 0 ) ⩽ 0 f(\bold{x_0})\leqslant{0} f(x0)⩽0(或找到一个非零向量 y 0 \bold{y_0} y0使 f ( C y 0 ) ⩽ 0 f(\bold{Cy_0})\leqslant{0} f(Cy0)⩽0即可完成证明
- 假设 f = y T Λ y = ∑ i = 1 n k i y i 2 f=\bold{y^{T}\Lambda{y}}=\sum_{i=1}^{n}k_iy_i^2 f=yTΛy=∑i=1nkiyi2是正定的,且存在 s ∈ { 1 , ⋯ , n } s\in\{1,\cdots,n\} s∈{1,⋯,n},对应 k s ⩽ 0 k_s\leqslant{0} ks⩽0
- 取 y = e s \bold{y=e_s} y=es时(其中 e s \bold{e}_s es是单位坐标向量), x 0 = C e s \bold{x}_0=\bold{Ce}_s x0=Ces, f ( x 0 ) = f ( C e s ) f(\bold{x_0})=f(\bold{Ce_{s}}) f(x0)=f(Ces)= g ( e s ) g(e_s) g(es)= k s ⩽ 0 k_s\leqslant{0} ks⩽0,
- 显然 C e s ≠ 0 \bold{Ce}_s\neq{\bold{0}} Ces=0,说明 f f f不是正定的;这与 f f f正定矛盾,从而 k i > 0 k_i>0 ki>0
- 事实上,这个证明过程可以从标准形开始,而不需要关心线性变换 x = C y \bold{x=Cy} x=Cy,因为任意二次型总是能够标准化,并且(可逆线性变换)标准化前后有相同的正定性
推论:正定矩阵和特征值
- 对称阵
A
\bold{A}
A为正定的充要条件是
A
\bold{A}
A的特征值全为正
- A \bold{A} A的特征值构成的对角阵 Λ = diag ( λ 1 , ⋯ , λ n ) \Lambda=\text{diag}(\lambda_1,\cdots,\lambda_n) Λ=diag(λ1,⋯,λn)是 f = x T A x f=\bold{x^{T}\bold{A}x} f=xTAx的一个标准形二次型 g ( y ) = y T Λ y g(\bold{y})=\bold{y^{T}\Lambda{y}} g(y)=yTΛy的矩阵,其正惯性指数等于 Λ \Lambda Λ对角元素中的正数个数
- 由惯性定理, f f f的正惯性指数和其任意标准形的正惯性指数一致,
- 所以,若 A \bold{A} A的特征值 λ 1 , ⋯ , λ n \lambda_1,\cdots,\lambda_n λ1,⋯,λn全为正,则 f f f的正惯性系数为 n n n,从而 f f f是正定的, A \bold{A} A是正定的
正定二次型(正定矩阵)性质
-
n n n元二次型 f = x T A x f=\bold{x^{T}\bold{A}x} f=xTAx= ∑ i ∑ j a i j x i x j \sum\limits_{i}\sum\limits_{j}a_{ij}x_ix_j i∑j∑aijxixj是正定二次型 ⇔ \Leftrightarrow ⇔ A = ( a i j ) \bold{A}=(a_{ij}) A=(aij)是正定矩阵
-
a i i > 0 , i = 1 , 2 , ⋯ , n a_{ii}>0,i=1,2,\cdots,n aii>0,i=1,2,⋯,n且 ∣ A ∣ > 0 |\bold{A}|>0 ∣A∣>0
-
证明:
-
因为 f f f是正定的,即当 x ≠ 0 \bold{x\neq{0}} x=0时恒满足 f ( x ) > 0 f(\bold{x})>0 f(x)>0,如果能够找到合适的非零向量 ϵ i \epsilon_i ϵi使得 f ( ϵ i ) = a i i f(\epsilon_i)=a_{ii} f(ϵi)=aii,那么自然得证明了 a i i > 0 , i = 1 , 2 , ⋯ , n a_{ii}>0,i=1,2,\cdots,n aii>0,i=1,2,⋯,n
-
令 ϵ i = ( 0 , ⋯ , 1 , ⋯ , 0 ) \epsilon_{i}=(0,\cdots,1,\cdots,0) ϵi=(0,⋯,1,⋯,0)(只有第 i i i个元素是非零元素,而且等于1, ϵ k \epsilon_k ϵk是单位坐标向量), i = 1 , 2 , ⋯ , n i=1,2,\cdots,n i=1,2,⋯,n
-
从而 f ( ϵ k ) = a k k , k = 1 , 2 , ⋯ , n f(\epsilon_k)=a_{kk},k=1,2,\cdots,n f(ϵk)=akk,k=1,2,⋯,n,又 f ( ϵ k ) > 0 f(\epsilon_k)>0 f(ϵk)>0,
-
所以 a i i > 0 a_{ii}>0 aii>0, i = 1 , 2 , ⋯ , n i=1,2,\cdots,n i=1,2,⋯,n
-
∣ A ∣ = ∏ i λ i > 0 |\bold{A}|=\prod_{i}\lambda_i>0 ∣A∣=∏iλi>0
-
-
-
矩阵 A \bold{A} A的特征值均大于0(上一节讨论过)
-
A \bold{A} A和同阶单位阵 E \bold{E} E合同
- 正定二次型 f f f的正惯性指数为 n n n,所以(其标准形 f = y T Λ y f=\bold{y^{T}\Lambda{y}} f=yTΛy矩阵 Λ \Lambda Λ全为正数),其规范形 f = z T E z f=\bold{z^{T}Ez} f=zTEz的矩阵就是单位阵
- 所以矩阵 A \bold{A} A与 n n n阶单位阵 E \bold{E} E合同,即存在可逆矩阵 Q \bold Q Q使得 Q T A Q = E \bold {Q^{T}\bold{A}Q=E} QTAQ=E
-
存在可逆矩阵 P \bold P P,使得 A = P T P \bold {\bold{A}=P^{T}P} A=PTP,即矩阵 A \bold{\bold{A}} A可以表示为两个互为转置矩阵的可逆矩阵的乘积
- 由于 A ≃ E \bold {\bold{A}\simeq{E}} A≃E,即存在可逆矩阵 Q \bold Q Q使得 Q T A Q = E \bold {Q^{T}\bold{A}Q=E} QTAQ=E
- 从而 A = ( Q T ) − 1 E Q − 1 = ( Q − 1 ) T E Q − 1 = ( Q − 1 ) T Q − 1 \bold{\bold{A}=(Q^{T})^{-1}EQ^{-1}=(Q^{-1})^{T}EQ^{-1}=(Q^{-1})^{T}Q^{-1}} A=(QT)−1EQ−1=(Q−1)TEQ−1=(Q−1)TQ−1
- 取 P = Q − 1 \bold {P=Q^{-1}} P=Q−1,即 A = P T P \bold {\bold{A}=P^{T}P} A=PTP
负定二次型
- 二次形 f ( x ) = x T A x f(\bold{x})=\bold{x^{T}\bold{A}x} f(x)=xTAx,如果对于任意的 x ≠ 0 \bold{x\neq{0}} x=0都有 f ( x ) < 0 f(\bold{x})<0 f(x)<0,称 f f f为正定二次型
- 矩阵 A \bold{\bold{A}} A是负定矩阵
负定二次型判定条件
-
和正定二次型的判定条件相仿:
-
f = x T A x f=\bold{x^{T}\bold{A}x} f=xTAx是负定二次型 ⇔ \Leftrightarrow ⇔ A \bold{A} A是负定矩阵
-
矩阵 A \bold{A} A的特征值均为负
-
A \bold{A} A的负惯性指数 n n n
-
k阶顺序主子式
-
设 A = ( a i j ) \bold{A}=(a_{ij}) A=(aij)为 n n n阶矩阵,正整数 k ⩽ n k\leqslant{n} k⩽n,则 A \bold{A} A的 k k k阶顺序主子式定义为 A \bold{A} A的前 k k k行和前 k k k列的交集元素,简称主子式
- Δ k = ∣ a 11 a 12 ⋯ a 1 k a 21 a 22 ⋯ a 2 k ⋮ ⋮ ⋱ ⋮ a k 1 a k 2 ⋯ a k k ∣ \Delta_k= \begin{vmatrix} a_{11}& a_{12}& \cdots & a_{1k} \\ a_{21}& a_{22}& \cdots & a_{2k} \\ \vdots & \vdots & \ddots & \vdots \\ a_{k1}& a_{k2}& \cdots & a_{kk} \end{vmatrix} Δk= a11a21⋮ak1a12a22⋮ak2⋯⋯⋱⋯a1ka2k⋮akk
-
一个 n n n阶方阵只有 n n n个主子式 ( k = 1 , 2 , ⋯ , n ) (k=1,2,\cdots,n) (k=1,2,⋯,n),且 n n n阶主子式是 A \bold{A} A本身
-
主子式是 k k k阶子式中的一种,它们的结果都是一个数
赫尔维茨定理:主子式判定二次型正定性和负定性
-
对称阵 A \bold{A} A是正定的当且仅当 A \bold{A} A的全部主子式均大于0,即
- Δ r > 0 \Delta_r>0 Δr>0, r = 1 , ⋯ , n r=1,\cdots,n r=1,⋯,n
-
对称阵 A \bold{A} A是的负定的充要条件是奇数阶主子式为负,偶数阶主子式为正,即
- Δ 1 , Δ 3 , ⋯ < 0 \Delta_1,\Delta_3,\cdots<0 Δ1,Δ3,⋯<0, Δ 2 , Δ 4 , ⋯ > 0 \Delta_2,\Delta_4,\cdots>0 Δ2,Δ4,⋯>0
- 可以更紧凑地表示为 ( − 1 ) r Δ r > 0 (-1)^{r}\Delta_r>0 (−1)rΔr>0, r = 1 , ⋯ , n r=1,\cdots,n r=1,⋯,n
二次型分类小结
有定二次型
- 设实二次型
f
=
x
T
A
x
f=\bold{x^{T}\bold{A}x}
f=xTAx对于任意
x
=
(
x
1
,
⋯
,
x
n
)
T
≠
0
x=(x_1,\cdots,x_n)^{T}\neq{0}
x=(x1,⋯,xn)T=0,:
- 若恒有 f > 0 f>0 f>0,则 f f f是正定二次型, A \bold{A} A为正定矩阵
- 若恒有 f ⩾ 0 f\geqslant{0} f⩾0,则 f f f是半正定二次型, A \bold{A} A为半正定矩阵
- 若恒有 f < 0 f<0 f<0,则称 f f f为负定二次型, A \bold{A} A为负定矩阵
- 若恒有 f ⩽ 0 f\leqslant{0} f⩽0,则称 f f f是半负定二次型, A \bold{A} A称为半负定矩阵
- 上述4类情况 地二次型是有定的
不定二次型
- 若二次型 f f f不是有定的,称为不定二次型
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
2022-09-08 java_从键盘输⼊⼀串字符串,存⼊ArrayList中,再将其 分别按字典序升序排序、按字典序降序排序、 按字符串长度降序排序、先按字符串长度降序排列 再按字典序排序。
2022-09-08 AM@三角级数三角函数系@傅里叶级数Fourier Series一般形式