[筆記][總結] MIT線性代數 Gilbert Strang 復矩陣

作者水平有限,歡迎大家提出文中錯誤

複線性空間

重新提一下,線性空間的定義

VV是一個非空集合,其元素x,y,zx,y,z被稱爲向量;KK是一個數域,有元素k,l,mk,l,mVV被稱爲一個向量空間或線性空間,當:

  1. VV中定義加法運算,當x,yVx,y\in V,有唯一的和x+yVx+y\in V,且加法滿足
    a) 結合律 x+(y+z)=(z+y)+zx+(y+z)=(z+y)+z
    b) 交換律 x+y=y+xx+y=y+x
    c) 存在零元00,使x+0=xx+0=x
    d) 每個元素都存在負元素,即xV,yV,x+y=0\forall x \in V,\exists y \in V,x+y=0,稱y爲x的負元素,記作x-x
  2. VV中定義數乘運算,當xV,kKx\in V,k\in K,有唯一的乘積kvVkv\in V,且數乘運算滿足:
    e) 數因子分配律 k(x+y)=kx+kyk(x+y)=kx+ky
    f) 分配律 (k+l)x=kx+lx(k+l)x=kx+lx
    g) 結合律 k(lx)=(kl)xk(lx)=(kl)x
    h) 1x=x1x=x

注意線性空間VV是基於數域KK的,一般的默認數域是實數域R\mathbb R,本部分將討論數域爲複數域C\mathbb C的線性空間的中的矩陣運算及性質。

從實數域到複數域

向量模長

實向量vv的模長是
v=vTv \left\|v\right\|=v^Tv
對於復向量zz
zTz=i=1nzizi z^Tz=\sum\limits_{i=1}^n z_i\cdot z_i
關鍵問題在於zizi=Re2{zi}Im2{zi}z_i\cdot z_i=Re^2\{z_i\}-Im^2\{z_i\}

ziziˉ=Re2{zi}+Im2{zi}z_i\cdot \bar{z_i}=Re^2\{z_i\}+Im^2\{z_i\}纔是正確的求和項
所以對於復向量,其模長爲
z=zˉTz \left\|z\right\|=\bar z^Tz
其中算子 ˉT\bar \cdot^T又被記作 H\cdot^H,叫做HermitianHermitian算子
z=zHz \left\|z\right\|=z^Hz

內積

相應的,復向量x,yx,y,內積爲xHyx^Hy

HermitianHermitian矩陣

對於實矩陣而言,要想使得特徵值都是實數,且存在相互正交的特徵向量,需要A=ATA=A^T,即實對稱矩陣

而對於復矩陣,同樣的要求需要AH=AA^H=A,爲HermitianHermitian矩陣

可見HermitianHermitian矩陣的對角元都是實數

酉矩陣

實標準正交矩陣有性質
qiTqj={0, if ij1, if i=j q_i^Tq_j= \left\{ \begin{aligned} 0,\ if\ i\not=j\\ 1,\ if\ i=j \end{aligned} \right.
QTQ=I Q^TQ=I
在復空間,與標準正交陣相對應的矩陣稱爲酉矩陣U=[u1 u2  un]U=[u_1\ u_2\ \cdots \ u_n]
U for unitaryU\ for\ unitary
uiHuj={0, if ij1, if i=j u_i^Hu_j= \left\{ \begin{aligned} 0,\ if\ i\not=j\\ 1,\ if\ i=j \end{aligned} \right.
UHU=I U^HU=I

FourierFourier矩陣

FourierFourier矩陣是十分出名的酉矩陣
Fn=[11111wnwn2wnn11wnn1wn2(n1)wn(n1)(n1)] F_n= \left[ \begin{matrix} 1&1&1&\cdots&1\\ 1&w_n&w^2_n &\cdots&w^{n-1}_n\\ \vdots&\vdots&\vdots&\ddots&\vdots\\ 1&w^{n-1}_n&w^{2(n-1)}_n&\cdots&w^{(n-1)(n-1)}_n\\ \end{matrix} \right]
注意到FnF_n是對稱的,但不是對稱矩陣,因爲FnF_n是復矩陣,另外也不是HermitianHermitian矩陣FourierFourier矩陣的各列是相互正交的,是一個酉矩陣

通項
(Fn)ij=wnij,  i,j=0,1(n1) (F_n)_{ij}=w^{ij}_n,\ \ i,j=0,1\cdots(n-1)
其中wn=ei2πnw_n=e^{i\frac{2\pi}{n}},可見wmnm=wnw^m_{m*n}=w_nwnw_n在複平面單位圓,輻角爲2πn\frac{2\pi}{n}

F4=[11111ii2i31i2i4i61i3i6i9]=[11111i1i11111i1i] F_4= \left[ \begin{matrix} 1&1&1&1\\ 1&i&i^2 &i^3\\ 1&i^2&i^4&i^6\\ 1&i^3&i^6&i^9\\ \end{matrix} \right]= \left[ \begin{matrix} 1&1&1&1\\ 1&i&-1 &-i\\ 1&-1&1&-1\\ 1&-i&-1&i\\ \end{matrix} \right]
F4F_4可以完成對一個四維複列向量的離散傅里葉變換

FourierFourier逆矩陣

Fn1F_n^{-1}同樣是十分有用的復矩陣,由於其列向量(或行向量)並不標準(並非單位向量),但是各列向量模長都是n\sqrt nFnHFn=nIF_n^HF_n=nI
Fn1=1nFnH=1nFˉn F_n^{-1}= \frac{1}{n}F^H_n=\frac{1}{n}\bar F_n
Fn1=1n[11111wˉnwˉn2wˉnn11wˉnn1wˉn2(n1)wˉn(n1)(n1)] F_n^{-1}=\frac{1}{n} \left[ \begin{matrix} 1&1&1&\cdots&1\\ 1&\bar w_n&\bar w^2_n &\cdots&\bar w^{n-1}_n\\ \vdots&\vdots&\vdots&\ddots&\vdots\\ 1&\bar w^{n-1}_n&\bar w^{2(n-1)}_n&\cdots&\bar w^{(n-1)(n-1)}_n\\ \end{matrix} \right]
上面F4F_4的逆矩陣爲
F41=14[11111i1i11111i1i] F_4^{-1}= \frac{1}{4} \left[ \begin{matrix} 1&1&1&1\\ 1&-i&-1 &i\\ 1&-1&1&-1\\ 1&i&-1&-i\\ \end{matrix} \right]

FFTFFT

FourierFourier矩陣可以使得時域序列變換到頻域,十分有用的矩陣,如果直接使用FnF_{n}左乘時域序列,其時間複雜度是O(n2)O(n^2),實際上FourierFourier矩陣可以被分解爲一系列稀疏矩陣,分解後對序列的變換算法的時間複雜度大大降低到O(nlogn)O(n\cdot logn),被稱爲快速傅里葉變換(Fast fourier transformation)(Fast\ fourier\ transformation),在[筆記][總結] MIT線性代數 Gilbert Strang 矩陣的應用中有介紹

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