矩陣論(五)——矩陣分析

1. 向量範數

向量範數: xV\forall x \in V,若非負實數x||x||滿足
   (1) 正定性:x0x=0    x=0||x|| \ge 0,且 ||x|| = 0 \iff x = 0
   (2) 齊次性:ax=a xaF||ax|| = |a|\ ||x||,a \in F
   (3) 三角不等式:xyVx+yx+y\forall x,y \in V,都有|| x + y|| \leq ||x|| + ||y||
則稱||x||爲向量x的範數[V.][V;||.||]賦範空間

1-範數: x1=Σixi||x||_1 = \Sigma_i |x_i|
2-範數(向量長度,由內積所誘導的範數): x2=(xx)=xHx=Σixi2||x||_2 = \sqrt{(x,x)} = \sqrt{x^H x} = \sqrt{\Sigma_i|x_i|^2}
\infty-範數: x=maxxi||x||_\infty = max|x_i|
p-範數: p(1, +)xp=ΣixippxCn\forall p \in (1,\ +\infty),||x||_p = \sqrt[p]{\Sigma_i |x_i|^p},\forall x \in C^n

例如:
x=(1, i, 1+i)Tx = (1,\ i,\ 1 + i)^T,有
x1=1+i+1+i=1+1+2=2+2||x||_1 = 1 + |i| + |1 + i| = 1 + 1 + \sqrt{2} = 2 + \sqrt{2}
x2=12+i2+1+i2=1+1+22=2||x||_2 = \sqrt{1^2 + |i|^2 + |1 + i|^2} = \sqrt{1 + 1 + \sqrt{2}^2} = 2
x=max{1, 1, 2}=2||x||_\infty = max\{1,\ 1,\ \sqrt{2}\} = \sqrt{2}

在這裏插入圖片描述
向量範數的連續性: α1, , αnCn\alpha_1,\ \cdots,\ \alpha_n爲C^n的任一組基,.Cn||.||爲C^n上任一向量範數,xCnx=Σi xi αixiCnf(x1, , xn)=xx1, , xn\forall x \in C^n,有x= \Sigma_i \ x_i \ \alpha_i,x_i \in C^n,則f(x_1,\ \cdots,\ x_n) = ||x||爲x_1,\ \cdots,\ x_n的連續函數

向量範數的等價性: x(1)x(2)||x||^{(1)}與||x||^{(2)}是線性空間V上定義的兩種向量範數,
c1,c2>0使c1x(2)x(1)c2x(2)xV\exists c_1, c_2 > 0,使c_1 ||x||^{(2)} \leq ||x||^{(1)} \leq c_2 ||x||^{(2)},\forall x \in V,則稱這兩個範數等價
有限維線性空間的任意兩種向量範數都是等價的
在無限維線性空間中,兩個向量範數是可以不等價的

2. 矩陣範數

矩陣範數: AFn×n\forall A \in F^{n \times n},對應一個非負實數||A||滿足
  (1) 正定性:A0A=0    A=0||A|| \ge 0,且||A|| = 0 \iff A = 0
  (2) 齊次性:aA=a AaF||aA|| = |a|\ ||A||,a \in F
  (3) 三角不等式:ABFn×nA+BA+B\forall A,B \in F^{n \times n},都有||A + B|| \leq ||A|| + ||B||
  (4) 相容性:ABFn×nABA B\forall A,B \in F^{n \times n},都有||AB|| \leq ||A|| \ ||B||
則稱||A||爲矩陣A的範數

F(Frobenius)-範數: AF=ΣiΣjaij2)=tr(AHA)=Σiσi2||A||_F = \sqrt{\Sigma_{i} \Sigma_{j} |a_{ij}|^2)} = \sqrt{tr(A^HA)} = \sqrt{\Sigma_i \sigma_i^2}

例如:
在這裏插入圖片描述

A=(aij)Cn×nAF=ΣiΣjaij2A = (a_{ij}) \in C^{n \times n},||A||_F = \sqrt{\Sigma_i \Sigma_j |a_{ij}|^2},則
(1)AF=AHF(1) \quad ||A||_F = ||A^H||_F
(2)UAF=AVF=UAVF=AF(2) \quad ||UA||_F = ||AV||_F = ||UAV||_F = ||A||_F,其中U,V是酉矩陣
(3)tr(AHA)=ΣiΣjaij2(3) \quad tr(A^H A) = \Sigma_i \Sigma_j |a_{ij}|^2

例如:
A=(03i1010112)AHA=(11211123i22+3i5)A = \begin{pmatrix} 0 & 3i & 1 \\ 0 & -1 & 0 \\ -1 & 1 & 2 \end{pmatrix},A^HA = \begin{pmatrix} 1 & -1 & 2 \\ -1 & 11 & 2-3i \\ -2 & 2 + 3i & 5 \end{pmatrix}
AF=9+1+1+1+1+4=tr(AHA)=1+11+5=17||A||_F = \sqrt{9 + 1 + 1 + 1 + 1 + 4} = \sqrt{tr(A^HA)} = \sqrt{1 + 11 + 5} = \sqrt{17}

相容範數: AxA.x||Ax|| \leq ||A|| .||x||,其中||x||是向量範數,||A||是矩陣範數

誘導範數: A=max{Axx}||A|| = max\{\frac{||Ax||}{||x||}\},其中||x||是向量範數且x0x \neq 0,稱||A||爲由向量範數||x||所誘導的誘導範數

矩陣p-範數:xp||x||_p所誘導的矩陣範數。常用的p-範數爲A1A2A||A||_1,||A||_2與||A||_\infty

列和範數: A1=max(Σi=1naij)||A||_1 = max(\Sigma_{i = 1}^n |a_{ij}|)np.max(np.sum(abs(arr), axis=1, keepdims=True), axis=0)

行和範數: A=max(Σj=1naij)||A||_\infty = max(\Sigma_{j = 1}^n |a_{ij}|)np.max(np.sum(abs(arr), axis=0, keepdims=True), axis=1)

譜範數: A2=λ1λ1AHA||A||_2 = \sqrt{\lambda_1},\lambda_1是A^HA的最大特徵值

例如:
在這裏插入圖片描述

3. 向量序列與矩陣序列的極限

3.1 向量序列的極限

x(k)=(x1(k)xn(k))Tk=12Cnx^{(k)} = (x_1^{(k)} \quad \cdots \quad x_n^{(k)})^T,k = 1,2,\cdots是C^n空間的一個向量序列,如果當k+k \rightarrow +\infty時,它的n個分量數列都收斂,即limkxi(k)=aii=12n\lim_{k \to \infty} x_i^{(k)} = a_i,i = 1,2,\cdots,n,則稱向量序列{x(k)}\{x^{(k)}\}是按分量收斂的。向量α=(α1αn)T\alpha = (\alpha_1 \quad \cdots \quad \alpha_n)^T是它的極限,記爲limkx(k)=αx(k)αlim_{k \rightarrow \infty} x^{(k)} = \alpha或 x^{(k)} \rightarrow \alpha
當至少有一個分量數列是發散的,則稱向量序列是發散的。
例如
在這裏插入圖片描述
x(k)CnαCnlimkx(k)=α    limkx(k)α=0x^{(k)} \in C^n,\alpha \in C^n,則\lim_{k \rightarrow \infty}x^{(k)} = \alpha \iff \lim_{k \rightarrow \infty} ||x^{(k)} - \alpha|| = 0,其中.Cn||.||爲C^n中任一範數

3.2 矩陣序列的極限

矩陣序列{A(k)}A(k)=(aij(k))Cn×nlimkaij(k)=aijij=1n\{A^{(k)}\},A^{(k)} = (a_{ij}^{(k)}) \in C^{n \times n},若\lim_{k \rightarrow \infty} a_{ij}^{(k)} = a_{ij},i,j = 1,\cdots,n
則稱矩陣序列{A(k)}\{A^{(k)}\}收斂,A=(aij(k)){A(k)}A = (a_{ij}^{(k)})稱爲\{A^{(k)}\}的極限,記爲limkA(k)=AA(k)=Ak\lim_{k \rightarrow \infty}A^{(k)} = A或A^{(k)} = A,k \rightarrow \infty
例如:
A(k)=(((1+1k)k1+1k1(1)kk)A=(e110)A^{(k)} = \begin{pmatrix} ((1 + \frac{1}{k}) ^ k & 1 + \frac{1}{k} \\ \\ -1 & \frac{(-1)^k}{k} \end{pmatrix} \longrightarrow A = \begin{pmatrix} e & 1 \\ \\ -1 & 0 \end{pmatrix}

{A(k)}Cn×nACn×nlimkA(k)=A    limkA(k)A=0\{A^{(k)}\} \in C^{n \times n},||A|| \in C^{n \times n},則lim_{k \rightarrow \infty} A^{(k)} = ||A|| \iff lim_{k \rightarrow \infty} ||A^{(k)} - A|| = 0,其中.Cn||.||爲C^n中任一範數

4. 矩陣冪級數

譜半徑: ρ(A)=max(λi)λi{λ1, λ2, , λn}{λ1, λ2, , λn}ACn×n\rho(A) = max(|\lambda_i|),\lambda_i \in \{\lambda_1,\ \lambda_2,\ \cdots,\ \lambda_n\},\{\lambda_1,\ \lambda_2,\ \cdots,\ \lambda_n\}矩陣A \in C^{n \times n}的全部特徵值

ρ(Ak)=(ρ(A))k\rho(A^k) = (\rho(A))^k
Ak0(k)    ρ(A)<1A^k \rightarrow 0(k \rightarrow \infty) \iff \rho(A) < 1

ACn×nACn×nρ(A)AA \in C^{n \times n},\forall 矩陣範數||A|| \in C^{n \times n},都有\rho(A) \leq ||A||。即A的譜半徑是A的任意一種矩陣範數的下界

ACn×nϵ>0A \in C^{n \times n},\forall \epsilon > 0,存在某種矩陣範數A使Aρ(A)+ϵ||A||_*,使||A||_* \leq \rho(A) + \epsilon。即A的譜半徑是A的所有矩陣範數的下確界
證明:
在這裏插入圖片描述
矩陣冪級數: Σk=0akAk=a0I+a1A++akAk+ACn×nakC\Sigma_{k = 0}^{\infty} a_k A^k = a_0 I + a_1 A + \cdots + a_k A^k + \cdots,其中A \in C^{n \times n},a_k \in C

矩陣冪級數的部分和: Sn(A)=Σk=0nakAkS_n(A) = \Sigma^n_{k = 0} a_k A^k

{Sn(A)}\{S_n(A)\}收斂,則稱Σk=0akAk\Sigma^\infty_{k = 0}a_kA^k收斂,否則發散

limnSn(A)=S\lim_{n \rightarrow \infty} S_n(A) = S,則稱S爲Σk=0akAk\Sigma_{k = 0}^\infty a_kA^k的和矩陣

收斂性判別: 若復變量z的冪級數Σk=0akzk\Sigma_{k = 0}^{\infty}a_k z^k的收斂半徑爲R,R=limnanan+1R = \lim_{n \to \infty} \frac{a_n}{a_{n + 1}},而方陣ACn×nA \in C^{n \times n}的譜半徑爲ρ(A)\rho(A),則
  (1) ρ(A)<RΣk=0akAk\rho(A) < R,則\Sigma_{k = 0}^\infty a_k A^k收斂
  (2) ρ(A)>RΣk=0akAk\rho(A) > R,則\Sigma_{k = 0}^\infty a_k A^k發散
  (3) ρ(A)=RΣk=0akAk\rho(A) = R,則\Sigma_{k = 0}^\infty a_k A^k收斂性不定

例如:
在這裏插入圖片描述
在這裏插入圖片描述

在這裏插入圖片描述
在這裏插入圖片描述

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章