到目前为止,本教材除了第6章计算 x^(T)x 时所遇到的平方和外,我们所关注的主要是线性方程。这类平方和及更一般形式的表达式称为二次型,它常常出现在线性代数在工程(设计标准和优化)和信号处理(例如输出的噪声功率)的应用中。它们也常常出现在物理学(例如势能和动能)、微分几何(例如曲面的法曲率)、经济学(例如效用函数)和统计学(例如置信椭圆体)中,某些这类应用实例的数学背景很容易转化为对对称矩阵的研究。
ℝⁿ 上的一个二次型是一个定义在 ℝⁿ 上的函数,它在向量 x 处的值可由表达式 Q(x) = x^(T)Ax 计算,其中 A 是一个 n× n 对称矩阵.矩阵 A 称为关于二次型的矩阵.
最简单的非零二次型是 Q(x) = x^(T)Ix = | x|².例1和例2说明了对称矩阵 A 和二次型 x^(T)Ax 之间的关系.
令 ±bx = x₁ x₂ ,计算下列矩阵的 x^(T) A x .
a. A = 4 & 0 0 & 3
b. A = 3 & -2 -2 & 7
解
\n b. 在 A 中有两个值为-2的元素,注意观察矩阵 A 中(1,2)元素如何出现在计算中,用粗体字给出.
r l boldsymbol x ^ T A boldsymbol x & = [ x _ 1 - x _ 2 ] [ c c 3 & - 2 - 2 & 7 ] [ l x _ 1 x _ 2 ] = [ x _ 1 - x _ 2 ] [ c 3 x _ 1 - 2 x _ 2 - 2 x _ 1 + 7 x _ 2 ] & = x _ 1 (3 x _ 1 - 2 x _ 2) + x _ 2 (- 2 x _ 1 + 7 x _ 2) & = 3 x _ 1 ^ 2 - 2 x _ 1 x _ 2 - 2 x _ 2 x _ 1 + 7 x _ 2 ^ 2 & = 3 x _ 1 ^ 2 - 4 x _ 1 x _ 2 + 7 x _ 2 ^ 2例1 b 中二次型中出现了项 -4x₁x₂ ,这是因为矩阵 对角线上出现了元素-2。相反,例1( 对角线上出现了元素-2。相反,例 1(a)中涉及的对角矩阵 A 的二次型中没有 x₁x₂ 的交叉乘积项。
对属于 ℝ³ 的 x ,取 Q(x) = 5x₁² + 3x₂² + 2x₃² - x₁x₂ + 8x₂x₃ ,写出 x^(T)Ax 形式的二次型.
解
x₁², x₂², x₃² 的系数仍在对角线上。为使 A 对称,当 i ≠ j 时, xᵢxⱼ 的系数必须平均分配给矩阵 A 中的 (i, j) 元素和 (j, i) 元素。 x₁x₃ 的系数是0,容易验证
Q (boldsymbol x) = boldsymbol x ^ T A boldsymbol x = [ l l l x _ 1 & x _ 2 & x _ 3 ] [ c c c 5 & - 1 / 2 & 0 - 1 / 2 & 3 & 4 0 & 4 & 2 ] [ l x _ 1 x _ 2 x _ 3 ]令 Q(±b x) = x₁² - 8x₁x₂ - 5x₂² ,计算 Q(±b x) 在 ±b x = -3 1 , 2 -2 和 1 -3 处的值.解
l Q (- 3, 1) = (- 3) ^ 2 - 8 (- 3) (1) - 5 (1) ^ 2 = 2 8 Q (2, - 2) = (2) ^ 2 - 8 (2) (- 2) - 5 (- 2) ^ 2 = 1 6 Q (1, - 3) = (1) ^ 2 - 8 (1) (- 3) - 5 (- 3) ^ 2 = - 2 0在某些情况下,没有交叉乘积项的二次型会显得更容易使用,也就是说,二次型对应的矩阵是对角矩阵。幸运的是,交叉乘积项可以通过用适当的变量代换来消去。
二次型的变量代换
Section titled “二次型的变量代换”如果 x 表示 Rⁿ 中的向量变量,那么变量代换是下面形式的等式:
±b x = P ±b y text 或 ±b y = P ^ - 1 ±b x (1)其中 P 是可逆矩阵且 ±by 是 ℝⁿ 中的一个新的向量变量.这里 P 的列可确定 ℝⁿ 的一个基,±by 是相对于该基的向量 ±bx 的坐标向量.(见4.4节.)
如果用变量代换(1)处理二次型 x^(T)Ax ,那么
boldsymbol x ^ T boldsymbol A boldsymbol x = (P boldsymbol y) ^ T A (P boldsymbol y) = boldsymbol y ^ T P ^ T A P boldsymbol y = boldsymbol y ^ T (P ^ T A P) boldsymbol y (2)且新的二次型矩阵是 P^(T)AP 。因为 A 是对称的,故由定理2,存在正交矩阵 P ,使得 P^(T)AP 是对角矩阵 D ,(2)中的二次型变为 y^(T)Dy ,这就是下面例题的解题思路。
求一个变量代换将例3 中的二次型变为一个没有交叉乘积项的二次型.
解
例3中二次型对应的矩阵是
A = [ c c 1 & - 4 - 4 & - 5 ]第一步是将矩阵 A 正交对角化. A 的特征值是 λ = 3 和 λ = -7 ,相应的单位特征向量是
λ = 3: [ c 2 / sqrt 5 - 1 / sqrt 5 ]; λ = - 7: [ c 1 / sqrt 5 2 / sqrt 5 ]这些特征向量自动正交(因为它们属于不同的特征值)且构成 ℝ² 的一个单位正交基.取
P = [ l l 2 / sqrt 5 & 1 / sqrt 5 - 1 / sqrt 5 & 2 / sqrt 5 ], D = [ c c 3 & 0 0 & - 7 ]那么 A = PDP⁻1 ,且 D = P⁻1AP = P^(T)AP ,像前面指出的那样。一个适当的变量代换是
x = Py,其中 x = x₁ x₂ , y = y₁ y₂
那么
r l x _ 1 ^ 2 - 8 x _ 1 x _ 2 - 5 x _ 2 ^ 2 & = mathbf x ^ T A mathbf x & = (P mathbf y) ^ T A (P mathbf y) & = mathbf y ^ T P ^ T A P mathbf y = mathbf y ^ T D mathbf y & = 3 y _ 1 ^ 2 - 7 y _ 2 ^ 2为了说明例4 中二次型相等的意义,我们可以利用新二次型计算 Q(x) 在 x=(2,-2) 时的值。首先,由于 x=Py,故
boldsymbol y = P ^ - 1 boldsymbol x = P ^ T boldsymbol x从而
mathbf y = [ l l 2 / sqrt 5 & - 1 / sqrt 5 1 / sqrt 5 & 2 / sqrt 5 ] [ l 2 - 2 ] = [ l 6 / sqrt 5 - 2 / sqrt 5 ]因此
r l 3 y _ 1 ^ 2 - 7 y _ 2 ^ 2 & = 3 (6 / sqrt 5) ^ 2 - 7 (- 2 / sqrt 5) ^ 2 = 3 (3 6 / 5) - 7 (4 / 5) & = 8 0 / 5 = 1 6这就是例3中 Q(x) 时的值,见图 7-2.x = (2, - 2) 时的值,见图7-2.
例4 说明了以下定理,定理的证明在例4 之前已基本上给出.
椭圆
设 A 是一个 n × n 对称矩阵,那么存在一个正交变量代换 x = Py,它将二次型 x^(T)Ax 变换为不含交叉乘积项的二次型 y^(T)Dy .
定理中矩阵 P 的列称为二次型 tildex^(T)Ax 的主轴,向量 y 是向量 x 在由这些主轴构造的 Rⁿ 空间的单位正交基下的坐标向量.
主轴的几何意义
Section titled “主轴的几何意义”设 Q(±b x) = ±bx^(T)±b A±bx ,其中 A 是一个 2× 2 可逆对称矩阵, c 是一个常数.可以证明 ℝ² 中所有满足
boldsymbol x ^ T boldsymbol A boldsymbol x = c (3)的 x 的集合对应于一个椭圆(或者圆)、双曲线、两条相交直线或单个点,或根本不含任意点。如果 A 是一个对角矩阵,如图 7 不是一个对角矩阵,如图 7-4所示,则(3)的图形是标准位置的旋转。找到主轴(由 A 不是一个对角矩阵,如图7-4所示,则(3)的图形是标准位置的旋转。找到主轴(由 A 的特征向量确定)等同于找到一个新的坐标系统,在该坐标系统下其图形是在标准位置下的图形。
图 7-4b 中的双曲线是方程 x^(T)Ax=16 的图形表示,其中 是例4 轴的正向是例4 中矩阵 是例 4 中的矩阵。图 7-4b 中 y₁ 轴的正向是例 4 中矩阵 P 第一列的方向, y₂ 轴的正向是矩阵 P 第二列的方向。
b) x₁²-8x₁x₂-5x₂²=16

a) 5x₁² - 4x₁x₂ + 5x₂² = 48
图 7-4a中,椭圆的方程为 5x₁² - 4x₁x₂ + 5x₂² = 48 ,求一个变量代换,将方程中的交叉乘积项消去.
解
二次型对应的矩阵为 A = 5 & -2 -2 & 5 , A 的特征值是3和7,对应的单位特征向量为
boldsymbol u _ 1 = [ l 1 / sqrt 2 1 / sqrt 2 ], boldsymbol u _ 2 = [ c - 1 / sqrt 2 1 / sqrt 2 ]令 P = [ u₁ u₂ ] = 1/√(2) & -1/√(2) 1/√(2) & 1/√(2) ,那么 P 可将 A 正交对角化,所以变量代换 x = Py 得到的二次型为 y^(T) D y = 3 y₁² + 7 y₂² ,变量代换的新坐标轴如图 7-4a 所示.
二次型的分类
Section titled “二次型的分类”当 A 是一个 n × n 矩阵时,二次型 Q(x) = x^(T)Ax 是一个定义域为 ℝⁿ 的实值函数。图 7-5显示了定义域为 ℝⁿ 的四个二次型的图形。对二次型 Q(x) 定义域中的每一个点 x = (x₁, x₂) ,可画出点 (x₁, x₂, z) ,其中 z = Q(x) 外,图 7 的值都是正的,图 7 的值都是负的。图 7-5a和图 7-5d中图形的水平截面是椭圆,而图 7-5c的水平截面是双曲线。x = 0 外,图7-5a中所有 Q(x) 的值都是正的,图7-5d中所有 Q(x) 的值都是负的。图7-5a和图7-5d中图形的水平截面是椭圆,而图7-5c的水平截面是双曲线。

a) z = 3x₁² + 7x₂²



图 7-5 中简单的 2 × 2 例子说明下列定义.
定义 一个二次型 Q 是:
a.正定的,如果对所有 x ≠ 0 ,有 Q(x) > 0 .
b.负定的,如果对所有 x ≠ 0 ,有 Q(x) < 0 。
c.不定的,如果 Q(x) 既有正值又有负值.
此外,Q 被称为半正定的,如果对所有 x, Q(x) ≥slant 0 ;Q 被称为半负定的,如果对所有 x, Q(x) ≤slant 0 。图 7-5 中 a 和 b 的二次型都是半正定的。但 a 最好描述为正定的。
根据特征值为二次型分类.见图 7-6.
设 A 是 n × n 对称矩阵,那么一个二次型 x^(T)Ax 是:
a.正定的,当且仅当 A 的所有特征值是正数.
b.负定的,当且仅当 A 的所有特征值是负数.
c.不定的,当且仅当 A 既有正特征值,又有负特征值.
证明
由主轴定理,存在一个正交的变量代换 x = P y ,使得
Q (boldsymbol x) = boldsymbol x ^ T A boldsymbol x = boldsymbol y ^ T D boldsymbol y = λ_ 1 y _ 1 ^ 2 + λ_ 2 y _ 2 ^ 2 + … + λ_ n y _ n ^ 2 (4)
其中 λ₁, λ₂, …, λₙ 是 A 的特征值.由于 P 是可逆的,因此非零向量 x 和非零向量 y 之间存在一个一一对应,这样,x ≠ 0 时,Q(x) 的值与 (4) 右边表达式的值完全一致.显然,像定理所描述的三类方式一样,它由特征值 λ₁, λ₂, …, λₙ 的符号所确定.
Q(±b x) = 3x₁² + 2x₂² + x₃² + 4x₁x₂ + 4x₂x₃ 是正定的吗?
解
由于所有项的系数是正数,因此二次型表面上看是正定的.但二次型对应的矩阵为
正定的

且 A 的特征值是 5, 2 和 -1,所以 Q(x) 是不定二次型,而不是正定二次型.
负定的

利用二次型的分类,相应地得到矩阵的形式分类.一个正定矩阵A是一个对称矩阵,二次型 x^(T)Ax 是正定的.其他形式的矩阵(如半正定矩阵)的概念可类似定义.
不定的
图 7-6
数值计算的注解 确定对称矩阵 A 是正定的最快捷的方式是尝试将矩阵 A 分解为 A = R^(T) R ,其中 R 是具有正对角线元素的上三角阵(可以采用一个略作修改的 LU 分解算
法). 这样的楚列斯基分解算法可行的充分必要条件是 A 是正定的.见第7章补充习题7. 练习题
根据矩阵 A 的特征值,给出一个半正定矩阵 A.
-
计算二次型 x^(T)Ax ,其中 A=5&1/3 1/3&1 a. x=x₁ x₂ b. x=6 1 c. x=1 3
-
计算二次型 x^(T)Ax ,其中 A=3&2&0 2&2&1 0&1&0 . a. x=x₁ x₂ x₃ b. x=-2 -1 5 c. x=1/√(2) 1/√(2) 1/√(2)
-
求二次型的矩阵,假设 x 属于 R² . a. 3x₁²-4x₁x₂+5x₂² b. 3x₁²+2x₁x₂
-
求二次型的矩阵,假设 x 属于 R² . a. 5x₁² + 16x₁x₂ - 5x₂² b. 2x₁x₂
-
求二次型的矩阵,假设 x 属于 R³ . a. 3x₁² + 2x₂² - 5x₃² - 6x₁x₂ + 8x₁x₃ - 4x₂x₃ b. 6x₁x₂ + 4x₁x₃ - 10x₂x₃
-
求二次型的矩阵,假设 x 属于 R³ . a. 3x₁²-2x₂²+5x₃²+4x₁x₂-6x₁x₃ b. 4x₃²-2x₁x₂+4x₂x₃
-
求一个变量代换 x = Py,将二次型 x₁² + 10x₁x₂ + x₂² 变换为没有交叉乘积项的形式,给出 P 和新的二次型.
-
设 A 是下列二次型的矩阵:
可以证明 A 的特征值是 3, 9 和 15,求正交矩阵 P,使得变量代换 x = Py 将 x^(T)Ax 变换为没有交叉乘积项的形式,给出 P 和新的二次型.
给出习题 9~18 的二次型分类,然后求一个变量代换 x = Py,将二次型变换成没有交叉乘积项的形式,写出新的二次型。利用 7.1 节的方法构造 P。
-
4x₁²-4x₁x₂+4x₂² 10. 2x₁²+6x₁x₂-6x₂² 11. 2x₁²-4x₁x₂-x₂² 12. -x₁²-2x₁x₂-x₂² 13. x₁²-6x₁x₂+9x₂² 14. 3x₁²+4x₁x₂
-
[M]-3x₁²-7x₂²-10x₃²-10x₄²+4x₁x₂+4x₁x₃+4x₁x₄+6x₃x₄
-
[M] 4x₁² + 4x₂² + 4x₃² + 4x₄² + 8x₁x₂ + 8x₃x₄ - 6x₁x₄ + 6x₂x₃
-
[M] 11x₁²+11x₂²+11x₃²+11x₄²+16x₁x₂-12x₁x₄+ 12x₂x₃+16x₃x₄
-
[M] 2x₁² + 2x₂² - 6x₁x₂ - 6x₁x₃ - 6x₁x₄ - 6x₂x₃ - 6x₂x₄ - 2x₃x₄
-
如果 x=(x₁,x₂) 和 x^(T)x=1 ,即 x₁²+x₂²=1 ,那么二次型 5x₁²+8x₂² 的最大值是什么?(用 x 的一些例子试试。)
-
如果 x^(T)x=1 ,那么二次型 5x₁²-3x₂² 的最大值是什么?习题 21~22 中,矩阵是 n× n 方阵且向量属于 Rⁿ ,判断每个命题的真假,验证你的结论.
-
a.二次型的矩阵是一个对称矩阵. b. 一个二次型没有交叉乘积项的充分必要条件是二次型的矩阵是对角矩阵. c. 二次型 x^(T)Ax 的主轴是 A 的特征向量. d. 一个正定二次型 Q 满足对所有 x ∈ ℝⁿ, Q(x) > 0 .
e.如果一个对称矩阵 A 的所有特征值是正的,那么二次型 x^(T)Ax 是正定的.
f.一个对称矩阵 A 的楚列斯基分解具有形式 A = R^(T)R ,其中上三角阵 R 具有正的对角线元素.
-
a.表达式 | x|² 不是一个二次型. b. 如果 A 对称且 P 是正交矩阵,那么变量代换 x = Py 将 x^(T)Ax 变换为没有交叉乘积项的二次型. c. 如果 A 是一个 2 × 2 对称矩阵,那么满足 x^(T)Ax = c ( c 是常数)的 x 的集合对应的几何图形是圆、椭圆或双曲线. d. 一个不定二次型既不是半正定形式也不是半负定形式. e. 如果 A 对称且对 x ≠ 0,二次型 x^(T)Ax 仅仅有负值,那么 A 的所有特征值是正的. 习题 23~24 表明,当 A = a & b b & d 且 det A ≠ 0 时,不用求出 A 的特征值就可将二次型 Q(x) = x^(T)Ax 分类.
-
如果 λ₁ 和 λ₂ 是 A 的特征值,那么 A 的特征多项式可以用两种方式写出: det(A - λ I) 和 (λ - λ₁)(λ - λ₂) 。利用这个结论说明 λ₁ + λ₂ = a + d ( a, d 是 A 的对角线上的元素),且 λ₁λ₂ = det A 。
-
验证下列命题: a.如果 det A > 0 且 a > 0,则 Q 是正定的. b. 如果 det A > 0 且 a < 0,则 Q 是负定的. c. 如果 det A < 0 ,则 Q 是不定的.
-
证明:如果 B 是 m × n 矩阵,那么 B^(T)B 是半正定的;如果 B 是 n × n 可逆矩阵,那么 B^(T)B 是正定的.
-
证明:如果 n × n 矩阵 A 是正定的,那么存在一个正定矩阵 B ,使得 A = B^(T)B (提示:写出 A = PDP^(T) 且 P^(T) = P⁻1 ,构造一个对角矩阵 C ,使得 D = C^(T)C ,且令 B = PCP^(T) ,说明 B 满足条件。)
-
令 A 和 B 是 n × n 对称矩阵且所有特征值都为正,证明 A + B 的特征值也是正的。(提示:考虑二次型。)
-
令 A 是 n × n 可逆对称矩阵,证明:如果二次型 x^(T)Ax 是正定的,那么二次型 x^(T)A⁻1x 也是正定的。(提示:考虑特征值。)
构造一个正交变量代换 x = Py ,并且像(4)式一样写出
boldsymbol x ^ T boldsymbol A boldsymbol x = boldsymbol y ^ T D boldsymbol y = λ_ 1 y _ 1 ^ 2 + λ_ 2 y _ 2 ^ 2 + … + λ_ n y _ n ^ 2如果一个特征值(例如 λᵢ )是负数,那么对应于 y = eᵢ ( Iₙ 的第 i 列)的 x 取值使得 x^(T)Ax 是负数,所以,一个半正定的二次型的所有特征值一定非负。反之,如果所有特征值非负,那么上面的展开式说明 x^(T)Ax 一定是半正定的。见图 7-7。