【线性代数与矩阵论】矩阵的谱半径与条件数

矩阵的谱半径与条件数

2023年11月18日


文章目录

  • 矩阵的谱半径与条件数
    • 1. 矩阵的谱半径
    • 2. 谱半径与范数的关系
    • 3. 矩阵的条件数
    • 下链


1. 矩阵的谱半径

定义 A ∈ C n × n {A\in \mathbb C^{n \times n} } ACn×n λ 1 , λ 2 , ⋯   , λ n { \lambda_1,\lambda_2,\cdots,\lambda_n} λ1,λ2,,λn 是A的特征值,则称
ρ ( A ) = max ⁡ 1 ≤ i ≤ n ∣ λ i ∣ \rho(A)=\max_{1\le i\le n}| \lambda_i| ρ(A)=1inmaxλi
为矩阵 A {A} A 的谱半径。矩阵的谱指的是一个矩阵的特征值的集合。

定理 A ∈ C n × n {A\in \mathbb C^{n \times n} } ACn×n ,则

  1. ρ ( A H ) = ρ ( A T ) = ρ ( A ) {\rho(A^ \mathrm H)=\rho(A^ \mathrm T)=\rho(A)} ρ(AH)=ρ(AT)=ρ(A)
  2. ρ ( A k ) = [ ρ ( A ) ] k {\rho(A^k)=[\rho(A)]^k} ρ(Ak)=[ρ(A)]k
  3. A {A} A正规矩阵时, ρ ( A ) = ∣ ∣ A ∣ ∣ 2 {\rho(A)= || A ||_{2 }} ρ(A)=∣∣A2

2. 谱半径与范数的关系

A ∈ C n × n {A\in \mathbb C^{n \times n}} ACn×n λ {\lambda} λ A {A} A 的一个特征值, x {x} x A {A} A 属于 λ { \lambda} λ 的特征向量,则
A x = λ x Ax= \lambda x Ax=λx
C n × n { \mathbb C^{n \times n} } Cn×n 中任一矩阵范数 ∣ ∣ ⋅ ∣ ∣ m { || \cdot ||_{m }} ∣∣m 以及与它相容的向量范数 ∣ ∣ ⋅ ∣ ∣ v { || \cdot ||_{v }} ∣∣v ,有
∣ λ ∣ ⋅ ∣ ∣ x ∣ ∣ v = ∣ ∣ λ x ∣ ∣ v = ∣ ∣ A x ∣ ∣ v ≤ ∣ ∣ A ∣ ∣ ⋅ ∣ ∣ x ∣ ∣ v | \lambda| \cdot || x ||_{v }= || \lambda x ||_{ v}= || Ax ||_{ v} \le || A ||_{ } \cdot || x ||_{v } λ∣∣xv=∣∣λxv=∣∣Axv∣∣A∣∣xv
从而
∣ λ ∣ ≤ ∣ ∣ A ∣ ∣ | \lambda| \le || A ||_{ } λ∣∣A
于是有如下定理。

定理 ρ ( A ) ≤ ∣ ∣ A ∣ ∣ { \rho(A)\le || A ||_{ }} ρ(A)∣∣A ,其中 ∣ ∣ A ∣ ∣ { || A ||_{ }} ∣∣A A {A} A 的任一矩阵范数。
定理 A ∈ C n × n {A\in \mathbb C^{n \times n} } ACn×n ,则任意 ϵ > 0 { \epsilon>0} ϵ>0 ,必存在 C n × n { \mathbb C^{n \times n} } Cn×n 上矩阵范数 ∣ ∣ ⋅ ∣ ∣ m { || \cdot ||_{m }} ∣∣m ,使得
∣ ∣ A ∣ ∣ m ≤ ρ ( A ) + ϵ || A ||_{ m} \le \rho(A)+ \epsilon ∣∣Amρ(A)+ϵ
也就是虽然矩阵范数可能大于谱近谱半径的矩阵半径,却又总是存在无限接范数。数值分析中,谱半径可以认为是算子范数。


3. 矩阵的条件数

引理 P ∈ C n × n {P\in \mathbb C^{n \times n} } PCn×n ,若对 C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ ∣ ∣ P ∣ ∣ < 1 { || P ||_{ }<1} ∣∣P<1 ,则 I − P {I-P} IP 可逆。
定理 A ∈ C n n × n {A\in \mathbb C_n^{n \times n} } ACnn×n 可逆, δ A ∈ C n × n { \delta A\in \mathbb C^{n \times n} } δACn×n 。若对 C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ ∣ ∣ A − 1 δ A ∣ ∣ < 1 { || A^{-1} \delta A ||_{ }<1} ∣∣A1δA<1 ,则

  1. A + δ A {A+ \delta A} A+δA 可逆
  2. ∣ ∣ ( A + δ A ) − 1 ∣ ∣ ≤ ∣ ∣ A − 1 ∣ ∣ 1 − ∣ ∣ A − 1 δ A ∣ ∣ {|| (A+ \delta A)^{-1} ||_{ }\le \frac{|| A^{-1} ||_{ }}{1- || A^{-1} \delta A ||_{ }}} ∣∣(A+δA)11∣∣A1δA∣∣A1
  3. ∣ ∣ A − 1 − ( A + δ A ) − 1 ∣ ∣ ∣ ∣ A − 1 ∣ ∣ ≤ ∣ ∣ A − 1 δ A ∣ ∣ 1 − ∣ ∣ A − 1 δ A ∣ ∣ {\frac{||A^{-1}- (A+ \delta A)^{-1} ||_{ }}{|| A^{-1} ||_{ }} \le \frac{|| A^{-1} \delta A ||_{ }}{1- || A^{-1} \delta A ||_{ }}} ∣∣A1∣∣A1(A+δA)11∣∣A1δA∣∣A1δA 矩阵扰动后逆矩阵的相对误差小于右端式子

定义 A ∈ C n n × n {A\in \mathbb C_n^{n \times n} } ACnn×n 可逆, ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数,称
cond ( A ) = ∣ ∣ A ∣ ∣ ⋅ ∣ ∣ A − 1 ∣ ∣ \text{cond}(A)= || A ||_{ } \cdot || A^{-1} ||_{ } cond(A)=∣∣A∣∣A1
为矩阵 A {A} A条件数

推论 A ∈ C n n × n {A\in \mathbb C_n^{n \times n} } ACnn×n 可逆, δ A ∈ C n × n { \delta A\in \mathbb C^{n \times n} } δACn×n 。若对 C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ ∣ ∣ A − 1 ∣ ∣ ⋅ ∣ ∣ δ A ∣ ∣ < 1 { || A^{-1}|| \cdot || \delta A ||_{ }<1} ∣∣A1∣∣∣∣δA<1 ,则
∣ ∣ A − 1 − ( A + δ A ) − 1 ∣ ∣ ∣ ∣ A − 1 ∣ ∣ ≤ ∣ ∣ A ∣ ∣ ∣ ∣ A − 1 ∣ ∣ ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ 1 − ∣ ∣ A ∣ ∣ ∣ ∣ A − 1 ∣ ∣ ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ = cond ( A ) ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ 1 − cond ( A ) ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ \frac{||A^{-1}- (A+ \delta A)^{-1} ||_{ }}{|| A^{-1} ||_{ }} \le \frac{|| A ||_{ }|| A^{-1} ||_{ }\frac{|| \delta A ||_{ }}{|| A ||_{ }}}{1-|| A ||_{ }|| A^{-1} ||_{ }\frac{|| \delta A ||_{ }}{|| A ||_{ }}}= \frac{\text{cond}(A) \frac{|| \delta A ||_{ }}{|| A ||_{ }}}{1-\text{cond}(A) \frac{|| \delta A ||_{ }}{|| A ||_{ }}} ∣∣A1∣∣A1(A+δA)11∣∣A∣∣A1∣∣A∣∣δA∣∣A∣∣A1∣∣A∣∣δA=1cond(A)∣∣A∣∣δAcond(A)∣∣A∣∣δA

定理 A ∈ C n n × n {A\in \mathbb C_n^{n \times n} } ACnn×n 可逆, δ A ∈ C n × n { \delta A\in \mathbb C^{n \times n} } δACn×n b {b} b δ b ∈ C n { \delta b\in \mathbb C^n} δbCn 。若对 C n × n { \mathbb C^{n \times n} } Cn×n 上的某一矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ ∣ ∣ A − 1 ∣ ∣ ⋅ ∣ ∣ δ A ∣ ∣ < 1 { || A^{-1}|| \cdot || \delta A ||_{ }<1} ∣∣A1∣∣∣∣δA<1 ,则非齐次线性方程组
A x = b     与     ( A + δ A ) ( x + δ x ) = b + δ b Ax=b \,\,\, 与 \,\,\, (A+ \delta A)(x+ \delta x)=b+ \delta b Ax=b(A+δA)(x+δx)=b+δb
的解满足
∣ ∣ δ x ∣ ∣ v ∣ ∣ x ∣ ∣ v ≤ cond ( A ) 1 − cond ( A ) ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ ( ∣ ∣ δ A ∣ ∣ ∣ ∣ A ∣ ∣ + ∣ ∣ δ b ∣ ∣ v ∣ ∣ b ∣ ∣ v ) \frac{|| \delta x ||_{ v}}{|| x ||_{v }}\le \frac{\text{cond}(A) }{1-\text{cond}(A) \frac{|| \delta A ||_{ }}{|| A ||_{ }}} \bigg( \frac{|| \delta A ||_{ }}{|| A ||_{ }}+\frac{|| \delta b ||_{ v}}{||b ||_{ v}} \bigg) ∣∣xv∣∣δxv1cond(A)∣∣A∣∣δAcond(A)(∣∣A∣∣δA+∣∣bv∣∣δbv)
其中 ∣ ∣ ⋅ ∣ ∣ v { || \cdot ||_{v }} ∣∣v C n { \mathbb C^n} Cn 上与矩阵范数 ∣ ∣ ⋅ ∣ ∣ { || \cdot ||_{ }} ∣∣ 相容的向量范数。
根据函数
f ( x ) = x 1 − x f(x)= \frac{x}{1-x} f(x)=1xx
的图像,当 cond ( A ) { \text{cond}(A)} cond(A) 很大,则说求逆或求解线性方程组是病态的。


下链


你可能感兴趣的:(线性代数与矩阵论,线性代数,矩阵,概率论,条件数)