二次型在工程计算、物理极值问题、优化理论与多元函数极值判别中有广泛的应用。在这些应用中,有一类特别重要的二次型——正定二次型。本节系统讨论正定二次型与正定矩阵的概念、等价判别准则(包括特征值判别法、标准形判别法、Sylvester 顺序主子式判别法),并简要介绍半正定与负定二次型。
一、正定二次型与正定矩阵
Section titled “一、正定二次型与正定矩阵”设实二次型 f(x₁, x₂, …, xₙ) = x^TAx。如果对任意 x ≠ 0,都有
x^TAx > 0则称 f 为正定二次型,称 A 为正定矩阵。
设 n 阶矩阵 A = (aᵢj) 为正定矩阵,求证:aᵢi > 0 (i = 1, 2, …, n)。
证
因为 f = x^TAx 是正定二次型,所以对任意 1 ≤ i ≤ n,取 x = εᵢ 为第 i 个分量为 1、其余分量为 0 的列向量,则 x ≠ 0,从而
x^TAx = aᵢi > 0下面给出判别二次型是否正定的充分必要条件。
实二次型 f(x₁, x₂, …, xₙ) = x^TAx 是正定二次型的充分必要条件为 f 的标准形中 n 个系数全为正数(即 f 的正惯性指数为 n)。
证明
设在可逆线性变换 x = Cy 下,f 化为标准形:
f = x^TAx = y^T(C^TAC)y = k₁y₁² + k₂y₂² + … + kₙyₙ²充分性:如果 kᵢ > 0 (i = 1, 2, …, n),那么对于任意的 x ≠ 0,由于 C 可逆,有 y = C⁻1x ≠ 0,于是
x^TAx = k₁y₁² + k₂y₂² + … + kₙyₙ² > 0这说明 f 是正定的。
必要性:对于任意取定的 1 ≤ i ≤ n,取 y = (y₁, y₂, …, yₙ)^T = εᵢ,即 yᵢ = 1,其他 yₖ (k ≠ i) 都取零值。令 x = Cy,则 x ≠ 0。因为 f 正定,所以
x^TAx = kᵢ > 0在正交变换下,二次型 f(x₁, x₂, …, xₙ) = x^TAx 的标准形中 n 个系数是 A 的全体特征值,所以有推论1。
n 阶实对称矩阵 A 正定的充分必要条件是 A 的 n 个特征值全为正数。
再由规范形的定义可得推论2。
n 阶实对称矩阵 A 正定的充分必要条件是 A 的合同规范形为 n 阶单位矩阵 E(即 A 与 E 合同)。
判别二次型 f(x₁, x₂, x₃) = x₁² + 2x₂² + x₃² + 2x₁x₂ + 4x₂x₃ 是否正定。
解
用配方法将 f 化为
f = (x₁ + x₂)² + (x₂ + 2x₃)² - 3x₃²f 的标准形中有 2 个系数为正,还有 1 个系数为负(即正惯性指数为 2 ≠ 3),所以 f 不是正定二次型。
接下来再给出一个判别正定矩阵的常用充分必要条件。
设 A = (aᵢj) 是 n 阶方阵。任取 1 ≤ i₁ < i₂ < … < iₖ ≤ n,由 A 的第 i₁, i₂, …, iₖ 行与第 i₁, i₂, …, iₖ 列交叉处的元素所构成的 k 阶行列式
aᵢ₁ i₁ & aᵢ₁ i₂ & … & aᵢ₁ iₖ aᵢ₂ i₁ & aᵢ₂ i₂ & … & aᵢ₂ iₖ vdots & vdots & & vdots aᵢₖ i₁ & aᵢₖ i₂ & … & aᵢₖ iₖ称为 A 的一个 k 阶主子式。
特别地,当 i₁ = 1, i₂ = 2, …, iₖ = k 时,这个 k 阶主子式称为 A 的 k 阶顺序主子式。
例如,设
A = 1 & -1 & 2 & 3 2 & 1 & -1 & 5 -1 & 3 & 0 & 0 3 & 4 & 1 & 2A 的三阶主子式有 4 个,如由第 2, 3, 4 行列构成的三阶主子式为
1 & -1 & 5 3 & 0 & 0 4 & 1 & 2 = 21A 的三阶顺序主子式只有 1 个,为
1 & -1 & 2 2 & 1 & -1 -1 & 3 & 0 = 16A 的二阶主子式有 6 个,但二阶顺序主子式只有 1 个,为
1 & -1 2 & 1 = 3n 阶实对称矩阵 A = (aᵢj) 是正定矩阵 iff A 的 n 个顺序主子式都为正数,即
a₁1 > 0, a₁1 & a₁2 a₂1 & a₂2 > 0, …, |A| > 0证*(Sylvester 定理证明)
必要性:已知 A 正定,即二次型 f = x^TAx 正定。
对任意 1 ≤ k ≤ n,记
Aₖ = a₁1 & a₁2 & … & a₁k a₂1 & a₂2 & … & a₂k vdots & vdots & & vdots aₖ1 & aₖ2 & … & aₖk将 A 分块表示为 A = Aₖ & B B^T & C,下面证 |Aₖ| > 0。
对任意 k 维非零向量 vₖ = (x₁, x₂, …, xₖ)^T ≠ 0,取 n 维向量 x^* = (x₁, …, xₖ, 0, …, 0)^T,则 x^* ≠ 0。因为 f = x^TAx 正定,所以
(x^*)^TAx^* = (vₖ^T, 0) Aₖ & B B^T & C vₖ 0 = vₖ^TAₖvₖ > 0这说明 k 元二次型 vₖ^TAₖvₖ 也是正定的,所以 Aₖ 的特征值全大于零,于是 |Aₖ| > 0。
充分性:用数学归纳法。
当 n = 1 时,A = (a₁1) 为一阶矩阵,已知 a₁1 > 0。因为
(frac1sqrta₁1)A(frac1sqrta₁1) = (1)即 A 与一阶单位矩阵合同,所以 A 是正定的。
假设当 n = k - 1 (k ≥ 2) 时,定理中的充分性成立。
当 n = k 时,A = (aᵢj) 为 k 阶矩阵,将 A 分块表示为
A = Aₖ-1 & α α^T & aₖk因为 A 的 k 个顺序主子式全大于零,由归纳假设可知 Aₖ-1 是正定的,所以 Aₖ-1 与 k - 1 阶单位矩阵 Eₖ-1 合同,即存在 k - 1 阶可逆矩阵 Cₖ-1,使得
Cₖ-1^TAₖ-1Cₖ-1 = Eₖ-1已知 |Aₖ-1| > 0,所以 Aₖ-1 可逆。
令
P = Eₖ-1 & -Aₖ-1⁻1α 0 & 1 (其中 0 表示零矩阵)则
P^TAP &= Eₖ-1 & 0 -α^TAₖ-1⁻1 & 1 Aₖ-1 & α α^T & aₖk Eₖ-1 & -Aₖ-1⁻1α 0 & 1 &= Aₖ-1 & 0 0 & aₖk - α^TAₖ-1⁻1α注意因为 Aₖ-1 是对称矩阵,所以 (Aₖ-1⁻1)^T = (Aₖ-1^T)⁻1 = Aₖ-1⁻1。
又因 |P| = 1,故
|A| = |P^TAP| = |Aₖ-1|(aₖk - α^TAₖ-1⁻1α)已知 |A| > 0, |Aₖ-1| > 0,所以 aₖk - α^TAₖ-1⁻1α > 0。
记 a = aₖk - α^TAₖ-1⁻1α > 0,再令
C = Cₖ-1 & 0 0 & (1)/(√(a))则
C^TP^TAPC &= Cₖ-1^T & 0 0 & (1)/(√(a)) Aₖ-1 & 0 0 & a Cₖ-1 & 0 0 & (1)/(√(a)) &= Cₖ-1^TAₖ-1Cₖ-1 & 0 0 & 1 = Eₖ这说明 A 与 k 阶单位矩阵 Eₖ 合同,所以 A 是正定的。定理得证。
判别二次型 f(x₁, x₂, x₃) = 3x₁² + 4x₂² + 5x₃² + 4x₁x₂ - 4x₂x₃ 是否正定。
解
写出二次型的矩阵形式 f = x^TAx,其中
x = x₁ x₂ x₃ , A = 3 & 2 & 0 2 & 4 & -2 0 & -2 & 5计算 A 的三个顺序主子式:
D₁ = 3 > 0, D₂ = 3 & 2 2 & 4 = 8 > 0, D₃ = |A| = 3(20 - 4) - 2(10) = 48 - 20 = 28 > 0由定理11 可知,A 为正定矩阵,f 为正定二次型。
求 k 为何值时,以下三元二次型为正定二次型:
(1) f(x₁, x₂, x₃) = (k + 1)x₁² + (k - 1)x₂² + (k - 2)x₃²;
(2) f(x, y, z) = 5x² + 4xy + y² - 2xz + kz² - 2yz。
解
(1) f 为标准形,f 是正定二次型当且仅当它的所有系数都是正数,即
k + 1 > 0 k - 1 > 0 k - 2 > 0 iff k > 2(2) f 对应的对称矩阵为
A = 5 & 2 & -1 2 & 1 & -1 -1 & -1 & k计算顺序主子式:
D₁ = 5 > 0, D₂ = 5 & 2 2 & 1 = 1 > 0, D₃ = 5 & 2 & -1 2 & 1 & -1 -1 & -1 & k = k - 2由定理11 可知,它是正定二次型当且仅当 D₃ > 0,即 k > 2。
二*、其他有定二次型
Section titled “二*、其他有定二次型”对于非正定的二次型,还可分为半正定、负定、半负定及不定等几种类型。下面分别给出定义与判别条件。
设 f(x₁, x₂, …, xₙ) = x^TAx 是实二次型。
- 半正定:如果对任意 x ≠ 0,都有 f(x₁, x₂, …, xₙ) = x^TAx ≥ 0,则称 f 为半正定二次型,称 A 为半正定矩阵。
- 负定:如果对任意 x ≠ 0,都有 f(x₁, x₂, …, xₙ) = x^TAx < 0,则称 f 为负定二次型,称 A 为负定矩阵。
- 半负定:如果对任意 x ≠ 0,都有 f(x₁, x₂, …, xₙ) = x^TAx ≤ 0,则称 f 为半负定二次型,称 A 为半负定矩阵。
- 不定:既非正定、半正定,又非负定、半负定的二次型(矩阵)称为不定二次型(不定矩阵)。
类似于定理10,有下面的定理12。
实二次型 f(x₁, x₂, …, xₙ) = x^TAx 是半正定二次型的充分必要条件是它的标准形中 n 个系数全为非负数。
定理12 的证明与定理10 的证明完全类似,在此不再赘述。
关于矩阵半正定的判别,没有类似于定理11 的关于顺序主子式的性质。这一点不难理解:
例如 A = 1 & 0 & 0 0 & 0 & 0 0 & 0 & -1,A 的三个顺序主子式为 D₁ = 1 ≥ 0, D₂ = 0 ≥ 0, D₃ = 0 ≥ 0,都非负,但二次型 f(x₁, x₂, x₃) = x^TAx = x₁² - x₃² 却不是半正定的(当取 x = (0, 0, 1)^T 时 f = -1 < 0)。
但我们有下面的定理:
实对称矩阵 A 半正定的充分必要条件是 A 的所有主子式均非负。
定理必要性的证明与本章定理11 中必要性的证明类似,但定理充分性的证明超出了本课程的要求,在此略去。
由定义可知,实对称矩阵 A 负定(半负定)当且仅当 -A 正定(半正定)。因此关于矩阵负定与半负定的判别就不必另做单独讨论了。
例如,A 为负定矩阵 iff 奇数阶顺序主子式全为负,偶数阶顺序主子式全为正,即 (-1)^k Dₖ > 0 (k = 1, 2, …, n)。