广义特征值与极大极小原理

第二十一讲 广义特征值与极小极大原理

一、 广义特征值问题

1、定义：设A 、B 为n 阶方阵，若存在数λ，使得方程Ax Bx =λ存在非零解，则称λ为A 相对于B 的广义特征值，x 为A 相对于B 的属于广义特征值λ的特征向量。

● 是标准特征值问题的推广，当B ＝I （单位矩阵）时，广义特征值问题退化为标准特征值问题。 ● 特征向量是非零的 ● 广义特征值的求解

()A B x 0-λ= 或者 ()B

A x 0λ-=

→

特征方程 ()det A B 0-λ=

求得λ后代回原方程Ax Bx =λ可求出x

本课程进一步考虑A 、B 厄米且为正定矩阵的情况。 2、等价表述

（1） B 正定，1B -存在

→1

B A x x

-=λ，广义特征值问题化为了标准

特征值问题，但一般来说，1B A -一般不再是厄米矩阵。 （2） B 正定，存在Cholesky 分解，H B G G =，G 满秩 H A x G G x =λ 令H G x y = 则 ()

1

1

H

G A G

y y --=λ 也成为标准特征值问题。

(

)

1

1

H

G A G

--为厄米矩阵，广义特征值是实数，可以按大小顺序

排列12n λ≤λ≤≤λ ，一定存在一组正交归一的特征向量，即存在

12n y ,y ,y 满足

()

1

1

H

i i G A G

y y --=λ

H i

j ij 1i j y y 0

i j

=⎧=δ=⎨≠⎩

还原为()1

H i i x G y -= (i=1,2, ,n)，则 ()()

H H H

H i

j i

j i

j ij 1

i j y y x G G x x Bx 0

i j

=⎧===δ=⎨

≠⎩ （带权正交）

二、 瑞利商

A 、

B 为n 阶厄米矩阵，且B 正定，称()()H

H

x A x R x x 0x Bx

=≠为A

相对于B 的瑞利商。

12n x ,x ,x 线性无关，所以，n

x C

∀∈，存在12n a ,a ,a C ∈ ，使

得 n

i

i

i 1

x a x ==

∑

H

n

n n

n

2

H

H

i i i j j j i

j i

i 1j 1i ,j 1

i 1

x Bx a x B a x a a x Bx a ====⎛⎫⎛⎫

==

=

⎪ ⎪⎝⎭⎝⎭

∑∑∑

∑

n

n

n

2

H

H H i i j i

j j i

i j i i

i ,j 1

i ,j 1

i 1

x A x a a x A x a a x Bx a ====

=

λ=

λ∑

∑

∑

∴ ()n

2

i i i 1n

2

i

i 1

a R x a ==λ=

∑

∑

●()1x 0

min R x ≠=λ ()n x 0

max R x ≠=λ

证明：()()()()

()

H

H

H

H

kx A kx x A x R x x Bx

kx B kx =

=

k 为非零常数

可取1k x

=，

kx 1=

∴ ()H

H

x 1

x A x R x x B x

==

（闭区域）

当1x

x =或()i a 0i 2,3,,n == 时，()1R x =λ

i 1λ≥λ ()n

2

i i 11

1n

2

i

i 1

a R x a ==≥λ=λ∑

∑

∴

()1x 0

min R x ≠=λ

另一方面，i n λ≤λ ()n

2

i i 1n

n n

2

i

i 1

a R x a ==≤λ=λ∑

∑

∴ ()n x 0

max R x ≠=λ

[证毕] 当B ＝I 时，标准特征值问题 A x x =λ （H A A =）

12n

H

i j ij

x x λ≤λ≤≤λ⎧⎨=δ⎩

则 ()

H

1H

x 0x A x m in

x x

≠=λ

()

H

n H

x 0x Ax max

x x

≠=λ

进一步分析可得

()

12x 0

a 0

m in R x ≠==λ ()

n n 1x 0

a 0

m ax R x -≠==λ

()

12k k 1x 0

a a a 0

m in R x +≠=====λ ()

n n 1n k n k 1x 0

a a a 0

m ax R x ----≠=====λ

定理1.设{}r r 1s L span x ,x ,,x += ()r r 1s +λ≤λ≤≤λ ,则 ()r x 0x L

m in R x ≠∈=λ ()s x 0x L

m ax R x ≠∈=λ