• 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=0g(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=C1x0=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. 规范形系数全为1
      2. 标准形所有系数为正数
    • 证明:
      • 设可逆线性变换 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} ks0
        • 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} ks0,
        • 显然 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是正定的

    正定二次型(正定矩阵)性质

    1. 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 ijaijxixj是正定二次型 ⇔ \Leftrightarrow A = ( a i j ) \bold{A}=(a_{ij}) A=(aij)是正定矩阵

    2. 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

    3. 矩阵 A \bold{A} A的特征值均大于0(上一节讨论过)

    4. 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
    5. 存在可逆矩阵 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}} AE,即存在可逆矩阵 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)1EQ1=(Q1)TEQ1=(Q1)TQ1
      • P = Q − 1 \bold {P=Q^{-1}} P=Q1,即 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} kn,则 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=
        |a11a12a1ka21a22a2kak1ak2akk|" role="presentation">|a11a12a1ka21a22a2kak1ak2akk|
        Δk= a11a21ak1a12a22ak2a1ka2kakk
    • 一个 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} f0,则 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} f0,则称 f f f半负定二次型, A \bold{A} A称为半负定矩阵
    • 上述4类情况 地二次型是有定

    不定二次型

    • 若二次型 f f f不是有定的,称为不定二次型
  • 相关阅读:
    EndNote21 | 杂志输出样式(参考文献格式模板)下载及安装
    如何批量图片重命名不同名字?
    Mybatis缓存
    【Android】adb无法连接设备原因及解决办法
    FPGA 学习笔记:Vivado 2019.1 工程创建
    并查集(Disjoint Set)
    UVM中config_db机制的使用方法
    本地部署FastGPT使用在线大语言模型
    56. 合并区间
    探究美颜算法:直播实时美颜SDK的集成和定制
  • 原文地址:https://blog.csdn.net/xuchaoxin1375/article/details/132767228