工程师、经济学家、科学家和数学家常常要寻找在一些特定集合内的 x 值,使得二次型 Q(x) 取最大值或最小值。具有代表性的是,这类问题可化为 x 是在一组单位向量中的变量的优化问题。下面我们将看到,这类条件优化问题有一个有趣且精彩的解。例6及后面的7.5节将讨论如何从实际中引出这类问题。
ℝⁿ 中的一个单位向量 ±bx 可用以下几种等价形式来描述:
| boldsymbol x | = 1, | boldsymbol x | ^ 2 = 1, boldsymbol x ^ T boldsymbol x = 1和
x _ 1 ^ 2 + x _ 2 ^ 2 + … + x _ n ^ 2 = 1 (1)在应用中经常使用 x^(T)x=1 的展开式(1).
当一个二次型没有交叉乘积项时,可以很容易得到在 x^(T)x=1 的条件下 Q(x) 的最大值和最小值.
求 Q(x)=9x₁²+4x₂²+3x₃² 在限制条件 x^(T)x=1 下的最大值和最小值.
解
由于 x₂² 和 x₃² 是非负的,因此
4 x _ 2 ^ 2 ≤slant 9 x _ 2 ^ 2, 3 x _ 3 ^ 2 ≤slant 9 x _ 3 ^ 2所以当 x₁² + x₂² + x₃² = 1 时,
r l Q (boldsymbol x) & = 9 x _ 1 ^ 2 + 4 x _ 2 ^ 2 + 3 x _ 3 ^ 2 & ≤slant 9 x _ 1 ^ 2 + 9 x _ 2 ^ 2 + 9 x _ 3 ^ 2 & = 9 (x _ 1 ^ 2 + x _ 2 ^ 2 + x _ 3 ^ 2) & = 9因此当 x 为单位向量时, Q(x) 的最大值不超过 9. 更进一步,当 x=(1,0,0) 时, Q(x)=9.从而 9 是 Q(x) 在 x^(T)x=1 条件下的最大值.
为求出 Q(±b x) 的最小值,注意到
9 x _ 1 ^ 2 ≥slant 3 x _ 1 ^ 2, 4 x _ 2 ^ 2 ≥slant 3 x _ 2 ^ 2因此当 x₁² + x₂² + x₃² = 1 时,
Q (±b x) ≥slant 3 x _ 1 ^ 2 + 3 x _ 2 ^ 2 + 3 x _ 3 ^ 2 = 3 (x _ 1 ^ 2 + x _ 2 ^ 2 + x _ 3 ^ 2) = 3当 x₁ = 0, x₂ = 0, x₃ = 1 时, Q(x) = 3 ,从而3是 Q(x) 在 ±bx^(T)±bx = 1 条件下的最小值.
从例1可以看到,二次型的矩阵具有特征值9,4和3,且最大和最小特征值分别等于在限制条件下的 Q(x) 的最大值和最小值.我们将会看到,此结论对任何二次型都是成立的.
令 A = 3 & 0 0 & 7 ,当 ±bx 属于 ℝ² 时, Q(±bx) = ±bx^(T)±bA±bx ,图 7 的图形表示。图 7-9表示圆柱体内部的一部分,圆柱与曲面的截面是点集 Q 的图形表示。图7-9表示圆柱体内部的一部分,圆柱与曲面的截面是点集 (x₁, x₂, z) ,表示在 x₁² + x₂² = 1 情况下的 z = Q(x₁, x₂) 。这些点的“高度值”是 Q(±bx) 的约束值,从几何意义上看,条件优化问题确定的是截面曲线上最高点和最低点的位置。
曲线上的两个最高点在 x₁x₂ 平面之上7个单位,出现在点 x₁ = 0 和 x₂ = ± 1 处,这些点对应于 A 的特征值7和特征向量 ±b x = (0,1) 及 -±b x = (0, - 1).类似地,曲线上的两个最低点在 x₁x₂ 平面之上3个单位,它们对应于特征值3和特征向量(1,0)及(-1,0).
图 7-9 中交线上的每一点对应的 z 坐标在 3 和 7 之间,且对任何 3 和 7 之间的数 t,存在一个单位向量 x 使得 Q(x)=t 。换言之,所有 x^(T)Ax 的可能值在 |x|=1 条件下的集合是闭区间 3 ≤slant t ≤slant 7 。
可以证明,对任何对称矩阵 A,在 |x|=1 条件下, x^(T)Ax 所有可能值的集合是实轴上的闭区间(见习题 13)。分别用 m 和 M 表示区间的左端点和右端点,即取
m = min {boldsymbol x ^ T A boldsymbol x: | boldsymbol x | = 1 }, M = max {boldsymbol x ^ T A boldsymbol x: | boldsymbol x | = 1 } (2)习题12要求证明,如果 λ 是 A 的一个特征值,那么 m ≤slant λ ≤slant M 。正像例2一样,下面的定理说明 m 和 M 自身也是 A 的特征值。 ②
设 A 是对称矩阵,且 m 和 M 的定义如(2)式所示,那么 M 是 A 的最大特征值 λ₁ , m 是 A 的最小特征值。如果 ±bx 是对应于 M 的单位特征向量 ±bu₁ ,那么 ±bx^(T)±bA±bx 的值等于 M 。如果 ±bx 是对应于 m 的单位特征向量,那么 ±bx^(T)±bA±bx 的值等于 m 。
证明
A 的正交对角化是 PDP⁻1 ,我们知道
text 当 ±b x = P ±b y text 时 , ±b x ^ T A ±b x = ±b y ^ T D ±b y (3)同样
text 对所有 y , | x | = | P y | = | y |这是因为 P^(T)P = I ,且 | Py| ² = (Py)^(T)Py = y^(T)P^(T)Py = y^(T)y = | y| ².特别地, | y| = 1 充分必要条件是 | x| = 1.这样,如果 ±bx 和 ±by 是任意单位向量,那么 x^(T)Ax 和 y^(T)Dy 可以认为是具有同样值的集合.为简化记号,假设 A 是以 a≥slant b≥slant c 为特征值的 3× 3 矩阵,重排 P 的列(特征向量),使得 P = [u₁,u₂,u₃] 和
D = [ c c c a & 0 & 0 0 & b & 0 0 & 0 & c ]对给定 ℝ³ 中的单位向量 y ,其坐标为 y₁,y₂,y₃ ,注意到
l a y _ 1 ^ 2 = a y _ 1 ^ 2 b y _ 2 ^ 2 ≤slant a y _ 2 ^ 2 c y _ 3 ^ 2 ≤slant a y _ 3 ^ 2将不等式相加,我们有
r l mathbf y ^ T D mathbf y & = a y _ 1 ^ 2 + b y _ 2 ^ 2 + c y _ 3 ^ 2 & ≤slant a y _ 1 ^ 2 + a y _ 2 ^ 2 + a y _ 3 ^ 2 & = a (y _ 1 ^ 2 + y _ 2 ^ 2 + y _ 3 ^ 2) & = a | mathbf y | ^ 2 = a这样由 M 的定义可知, M ≤slant a 。然而,当 y = e₁ = (1,0,0) 时, y^(T)Dy = a ,所以,事实上 M = a 。由(3)可知,对应于 y = e₁ 的向量 x 是矩阵 A 的特征向量 u₁ ,其原因是
±b x = P ±b e _ 1 = [ ±b u _ 1 ±b u _ 2 ±b u _ 3 ] [ l 1 0 0 ] = ±b u _ 1这样 M = a = e₁^(T)De₁ = u₁^(T)Au₁ ,从而证明了关于 M 的论断.同理可证明 m 是最小特征值即 c 值,且 ±bx^(T)A±bx 的这个值当 ±b x = P±be₃ = ±bu₃ 时可以取到.
令 A = 3 & 2 & 1 2 & 3 & 1 1 & 1 & 4 ,求二次型 ±bx^(T)±bA±bx 在限制条件 ±bx^(T)±bx = 1 下的最大值,并求一个可以取到
该最大值的单位向量.
解
由定理6,只需求出 A 的最大特征值,其特征多项式是
0 = - λ^ 3 + 1 0 λ^ 2 - 2 7 λ + 1 8 = - (λ - 6) (λ - 3) (λ - 1)最大特征值为 6.
限制条件下 x^(T)Ax 的最大值可以在特征值 λ=6 对应的单位特征向量 x 处获得.解 (A-6I)x=0 ,可得特征向量 1 1 1 和 u₁=1/√(3) 1/√(3) 1/√(3) .
在定理7和稍后的应用实例中,要计算 x^(T)Ax 定理7设 ±bx 如定理6所示.在如下条件限制下:A, λ₁ 和 ±bu₁ 如定理6所示.在如下条件限制下:
boldsymbol x ^ T boldsymbol x = 1, boldsymbol x ^ T boldsymbol u _ 1 = 0x^(T)Ax 的最大值是第二大特征值 λ₂ ,且这个最大值可以在 x 是对应于 λ₂ 的特征向量 u₂ 处达到.
定理7的证明和上面讨论的类似,即定理化为二次型的矩阵是对角矩阵的情形,下面的例子给出对角矩阵情形下证明的思路.
求 9x₁² + 4x₂² + 3x₃² 的最大值,其限制条件是 ±bx^(T)±bx = 1 和 ±bx^(T)±bu₁ = 0 ,其中 ±bu₁ = (1,0,0) 。注意到 ±bu₁ 是二次型的矩阵的最大特征值 λ = 9 对应的单位特征向量。
解
如果 x 的坐标为 x₁, x₂, x₃ ,那么限制 x^(T)u₁=0 简单意味着 x₁=0 。对这样的一个单位向量, x₂²+x₃²=1 且
r l 9 x _ 1 ^ 2 + 4 x _ 2 ^ 2 + 3 x _ 3 ^ 2 & = 4 x _ 2 ^ 2 + 3 x _ 3 ^ 2 & ≤slant 4 x _ 2 ^ 2 + 4 x _ 3 ^ 2 & = 4 (x _ 2 ^ 2 + x _ 3 ^ 2) & = 4在这样的限制条件下,二次型的最大值不超过4,且这个最大值可在 x=(0,1,0) 处达到,而这就是二次型的矩阵的第二大特征值对应的特征向量.
令 A 表示例3 中的矩阵,且 u₁ 是对应于矩阵 A 的最大特征值的特征向量,求 x^(T)Ax 的最大值,其限制条件是
boldsymbol x ^ T boldsymbol x = 1, boldsymbol x ^ T boldsymbol u _ 1 = 0 (4)解
由例3可知, A 的第二大特征值是 λ = 3 。解方程 (A - 3I)x = 0 求出 λ = 3 对应的特征向量,且单位化得到
boldsymbol u _ 2 = [ c 1 / sqrt 6 1 / sqrt 6 - 2 / sqrt 6 ]由于特征向量对应于不同的特征值,故向量 u₂ 和 u₁ 自然互相垂直。这样,在条件(4)的限制下, x^(T)Ax 的最大值是 3,且在 x = u₂ 处可以达到。
下面的定理是定理7的推广,它和定理6一起给出矩阵 A 所有特征值的有用特性。具体证明从略。
设 A 是一个 n × n 对称矩阵,其正交对角化为 A = PDP⁻1 ,将对角矩阵 D 上的元素重新排列,使得 λ₁ ≥slant λ₂ ≥slant … ≥slant λₙ ,且 P 的列是其对应的单位特征向量 u₁, …, uₙ 。那么对 k = 2, …, n ,在以下限制条件下:
boldsymbol x ^ T boldsymbol x = 1, boldsymbol x ^ T boldsymbol u _ 1 = 0, … , boldsymbol x ^ T boldsymbol u _ k - 1 = 0x^(T)Ax 的最大值是特征值 λₖ ,且这个最大值在 x=uₖ 处可以达到.
定理8 在 7.4 节中和 7.5 节中非常有用,下面的应用仅仅需要定理6.
在下一年度,一个县政府计划修 x 百英里的公路和桥梁,并且修整 y 百英亩的公园和娱乐场所。政府部门必须确定在两个项目上如何分配它的资源(资金、设备和劳动等)。如果更划算的话,可以同时开始两个项目,而不是仅开始一个项目,那么 x 和 y 必须满足下面的限
制条件:
4 x ^ 2 + 9 y ^ 2 ≤slant 3 6见图 7-10. 每个阴影可行集中的点 (x,y) 表示一个可能的该年度的公共工作计划。在限制曲线 4x²+9y²=36 上的点使资源利用达到最大可能。

为选择它的公共工作计划,县政府需要考虑居民的意见。为度量居民分配各类工作计划 (x, y) 的值或效用,经济学家有时利用下面的函数:
q (x, y) = x y其中使 q(x,y) 为常数的 (x,y) 点的集合称为无差别曲线。从图 7-11 中可以看到三条这样的曲线,沿着无差别曲线的点对应的选择表示居民作为一个群体有相同的价值观 ① 。求公共工作计划,使得效用函数 q 最大。
解
限制条件的方程 4x² + 9y² = 36 并没有描述一个单位向量集,但变量代换可以修正这个问题.重写限制条件为如下形式:
(frac x3) ^ 2 + (frac y2) ^ 2 = 1定义
x₁ = x/3, x₂ = y/2,即 x = 3x₁, y = 2x₂
从而限制条件变成
x _ 1 ^ 2 + x _ 2 ^ 2 = 1效用函数变成 q(3x₁,2x₂) = (3x₁)(2x₂) = 6x₁x₂ 。取 ±bx = x₁ x₂ ,那么原问题变为,在限制条件 ±bx^(T)±bx = 1 下 Q(±bx) = 6x₁x₂ 的最大值。注意到 Q(±bx) = ±bx^(T)A±bx ,其中
A = [ l l 0 & 3 3 & 0 ]A 的特征值是 ± 3 ,对应于 λ = 3 的特征向量是 [ l1 / √(2) 1 / √(2) ] ,对应于 λ = -3 的特征向量是 [ l - 1 / √(2) 1 / √(2) ] 这样 Q(x) = q(x₁,x₂) 的最大值是3,在 x₁ = 1 / √(2) 和 x₂ = 1 / √(2) 处可以达到.
根据原来的变量,最优的公共工作计划是修建 x = 3x₁ = 3 / √(2)≈ 2.1 百英里的公路和桥梁以及 y = 2x₂ = √(2)≈ 1.4 百英亩的公园和娱乐场所.最优公共工作计划是限制曲线和无差别曲线 q(x,y) = 3 恰好相交的点,具有更大效用的点 (x,y) 位于和限制曲线不相交的无差别曲线上,见图 7-11.
- 设 Q(x)=3x₁²+3x₂²+2x₁x₂ ,求一个变量代换,将 Q 变换成一个不含交叉乘积项的二次型,且给出新二次型.
- 对练习题 1 中的 Q,求 Q(x) 在限制条件 x^(T)x=1 下的最大值,并求出达到最大值的单位向量.
在习题 1 和习题 2 中,求变量代换 x = Py,将二次型 x^(T)Ax 变换为对应的 y^(T)Dy .
-
5x₁² + 6x₂² + 7x₃² + 4x₁x₂ - 4x₂x₃ = 9y₁² + 6y₂² + 3y₃²
-
3x₁² + 3x₂² + 5x₃² + 6x₁x₂ + 2x₁x₃ + 2x₂x₃ = 7y₁² + 4y₂²
(提示:x 和 y 必须有同样数目的坐标,这样,这里显示的二次型关于 y₃² 的系数一定为零。)
对习题 3~6,求(a)在条件 x^(T)x=1 限制下 Q(x) 的最大值,(b)达到最大值的一个单位向量 u,(c) 在条件 x^(T)x=1 和 x^(T)u=0 限制下 Q(x) 的最大值.
-
Q(x)=5x₁²+6x₂²+7x₃²+4x₁x₂-4x₂x₃ (见习题1)
-
Q(x)=3x₁²+3x₂²+5x₃²+6x₁x₂+2x₁x₃+2x₂x₃ (见习题2)
-
Q(x)=x₁²+x₂²-10x₁x₂
-
Q(x)=3x₁²+9x₂²+8x₁x₂
-
设 Q(x) = -2x₁² - x₂² + 4x₁x₂ + 4x₂x₃ ,在条件 ±bx^(T)±bx = 1 的限制下,求出在 ℝ³ 中的单位向量 ±bx 使 Q(x) 取最大值。(提示:二次型 Q 的矩阵的特征值是 2,-1 和 -4.)
-
设 Q(x) = 7x₁² + x₂² + 7x₃² - 8x₁x₂ - 4x₁x₃ - 8x₂x₃ ,在条件 ±bx^(T)±bx = 1 的限制下,求出 ℝ³ 中的单位向量 ±bx ,使得 Q(±bx) 最大。(提示:二次型 Q 的矩阵的特征值是9和-3。)
-
在条件 x₁² + x₂² = 1 的限制下,求出 Q(x) = 7x₁² + 3x₂² - 2x₁x₂ 的最大值。(不必求出达到最大值的向量。)
-
在条件 x₁² + x₂² = 1 的限制下,求出 Q(x) = -3x₁² + 5x₂² - 2x₁x₂ 的最大值。(不必求出达到最大值的向量。)
-
若 x 是矩阵 A 对应于特征值 3 的一个单位特征向量, x^(T)Ax 的值是什么?
-
设 λ 是对称矩阵 A 的一个特征值,验证本节的一个结论,即 m ≤slant λ ≤slant M ,其中 m 和 M 的定义在(2)式中给出。(提示:求出 x 使得 λ = x^(T)Ax .)
-
设 A 是 n × n 对称矩阵,M 和 m 表示二次型 x^(T)Ax 的最大值和最小值,其中 x^(T)x = 1 ,对应的单位特征向量是 u₁ 和 uₙ 。下面的计算表明,对任意给定的介于 M 和 m 之间的数 t,有一个单位向量 x,使得 t = x^(T)Ax 。验证对 0 和 1 之间的数 α ,有 t = (1 - α)m + α M ,然后令 x = √(1 - α)uₙ + √(α)u₁ ,证明 x^(T)x = 1 且 x^(T)Ax = t 。[M]请依据习题 3~6 给出的说明做习题 14~17。
-
3x₁x₂ + 5x₁x₃ + 7x₁x₄ + 7x₂x₃ + 5x₂x₄ + 3x₃x₄
-
4x₁² - 6x₁x₂ - 10x₁x₃ - 10x₁x₄ - 6x₂x₃ - 6x₂x₄ - 2x₃x₄
-
-6x₁²-10x₂²-13x₃²-13x₄²-4x₁x₂-4x₁x₃-4x₁x₄+6x₃x₄
-
x₁x₂ + 3x₁x₃ + 30x₁x₄ + 30x₂x₃ + 3x₂x₄ + x₃x₄
-
二次型的矩阵是 A = 3 & 1 1 & 3 ,很容易求出特征值为4和2,对应的单位特征向量为 1 / √(2) 1 / √(2) 和 -1 / √(2) 1 / √(2) ,所以期望的变量代换是 x = Py ,其中 P = 1 / √(2) & -1 / √(2) 1 / √(2) & 1 / √(2) 。(一个常见的错误是忘记将向量单位化。)新的二次型是 y^(T)Dy = 4y₁² + 2y₂² 。
-
对单位向量 x, Q(x) 的最大值是 4,且取得最大值时的单位向量是 1/√(2) 1/√(2) (一个常见的错误答案是 1 0 ,这个向量使二次型 y^(T)Dy 而不是 Q(x) 达到最大值)。见图 7-12.