跳转到内容

7.5 正定二次型与正定矩阵

[美] David C. Lay、Steven R. Lay 等著,刘深泉等译机械工业出版社《线性代数及其应用(原书第5版)》大学本科 · 教材双语助读 (Alt+Y)

二次型在工程计算、物理极值问题、优化理论与多元函数极值判别中有广泛的应用。在这些应用中,有一类特别重要的二次型——正定二次型。本节系统讨论正定二次型与正定矩阵的概念、等价判别准则(包括特征值判别法、标准形判别法、Sylvester 顺序主子式判别法),并简要介绍半正定与负定二次型。

定义 14 正定二次型与正定矩阵

设实二次型 f(x₁, x₂, …, xₙ) = x^TAx。如果对任意 x ≠ 0,都有

x^TAx > 0

则称 f 为正定二次型,称 A 为正定矩阵。

例 19 正定矩阵的主对角元均为正数

设 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

下面给出判别二次型是否正定的充分必要条件。

定理 10 正定二次型的充要条件(标准形判别)

实二次型 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。

推论 1 特征值判别法

n 阶实对称矩阵 A 正定的充分必要条件是 A 的 n 个特征值全为正数。

再由规范形的定义可得推论2。

推论 2 合同规范形判别法

n 阶实对称矩阵 A 正定的充分必要条件是 A 的合同规范形为 n 阶单位矩阵 E(即 A 与 E 合同)。

例 20 配方法判定二次型正定性

判别二次型 f(x₁, x₂, x₃) = x₁² + 2x₂² + x₃² + 2x₁x₂ + 4x₂x₃ 是否正定。

解

用配方法将 f 化为

f = (x₁ + x₂)² + (x₂ + 2x₃)² - 3x₃²

f 的标准形中有 2 个系数为正,还有 1 个系数为负(即正惯性指数为 2 ≠ 3),所以 f 不是正定二次型。

接下来再给出一个判别正定矩阵的常用充分必要条件。

定义 15 主子式与顺序主子式

设 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 & 2

A 的三阶主子式有 4 个,如由第 2, 3, 4 行列构成的三阶主子式为

1 & -1 & 5 3 & 0 & 0 4 & 1 & 2 = 21

A 的三阶顺序主子式只有 1 个,为

1 & -1 & 2 2 & 1 & -1 -1 & 3 & 0 = 16

A 的二阶主子式有 6 个,但二阶顺序主子式只有 1 个,为

1 & -1 2 & 1 = 3
定理 11 赫尔维茨定理(Sylvester 准则)

n 阶实对称矩阵 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 是正定的。定理得证。

例 21 顺序主子式判定二次型正定性

判别二次型 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 为正定二次型。

例 22 待定参数使二次型正定的取值范围

求 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。

对于非正定的二次型,还可分为半正定、负定、半负定及不定等几种类型。下面分别给出定义与判别条件。

定义 16 其他有定二次型与不定二次型

设 f(x₁, x₂, …, xₙ) = x^TAx 是实二次型。

  1. 半正定:如果对任意 x ≠ 0,都有 f(x₁, x₂, …, xₙ) = x^TAx ≥ 0,则称 f 为半正定二次型,称 A 为半正定矩阵。
  2. 负定:如果对任意 x ≠ 0,都有 f(x₁, x₂, …, xₙ) = x^TAx < 0,则称 f 为负定二次型,称 A 为负定矩阵。
  3. 半负定:如果对任意 x ≠ 0,都有 f(x₁, x₂, …, xₙ) = x^TAx ≤ 0,则称 f 为半负定二次型,称 A 为半负定矩阵。
  4. 不定:既非正定、半正定,又非负定、半负定的二次型(矩阵)称为不定二次型(不定矩阵)。

类似于定理10,有下面的定理12。

定理 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)。

但我们有下面的定理:

定理 13 半正定矩阵的主子式充要条件

实对称矩阵 A 半正定的充分必要条件是 A 的所有主子式均非负。

定理必要性的证明与本章定理11 中必要性的证明类似,但定理充分性的证明超出了本课程的要求,在此略去。

负定与半负定的判别转化

由定义可知,实对称矩阵 A 负定(半负定)当且仅当 -A 正定(半正定)。因此关于矩阵负定与半负定的判别就不必另做单独讨论了。

例如,A 为负定矩阵 iff 奇数阶顺序主子式全为负,偶数阶顺序主子式全为正,即 (-1)^k Dₖ > 0 (k = 1, 2, …, n)。

共 0 题