2021-04-30 離散傅里葉變換DFT的純數(shù)學(xué)理解

本文是介紹離散傅里葉變換的,實(shí)際上筆者看過了n多文章或者書籍介紹傅里葉變換

,但是關(guān)于從傅里葉級(jí)數(shù)到連續(xù)傅里葉變換,再到離散時(shí)間傅里葉變換,再到離散傅列葉變換,尚未發(fā)現(xiàn)有文章能把整個(gè)邏輯鏈介紹清楚的。

因此筆者另辟蹊徑,從一個(gè)非常獨(dú)立的視角來

推導(dǎo)離散傅列葉變換公式。

問題:給定一個(gè)離散復(fù)信號(hào)x[t] (0<=t<N ,t\in Z)
以及給定以下N個(gè)離散復(fù)信號(hào)
e_k[t] = e^{2\pi kt\sqrt{-1}/N}/\sqrt{N} (0<=t<N , 0<=k<N ,n,t\in Z)

現(xiàn)在要求復(fù)數(shù):X[0],X[1],...X[N-1]使得
對(duì)任意0<=t<N,t\in Z有:
x[t]=\sum_{k=0}^{N-1}X[k] e_k[t] (1)

解:
那么,我們可以把上面N個(gè)等式,表達(dá)為矩陣形式:
X=(X[0],X[1],...,X[N-1])^T
x=(x[0],x[1],...,x[N-1])^T
a_{i,j} =e_j[i]
A=(a_{i,j})
[N*N矩陣,第i+1行,j+1列為 a_{i,j} ,后文中對(duì)于一般的矩陣T,T_{i,j}表示其i+1行,j+1列元素,不再贅述]
則,方程組(1)可以簡(jiǎn)寫為:
x=AX
=>
X=A^{-1}x

下面證明A是酉矩陣(逆矩陣為共軛轉(zhuǎn)置的矩陣)
設(shè)A的共軛轉(zhuǎn)置矩陣為B
B_{i,j}= \overline e_i[j] =e^{-2\pi ij\sqrt{-1}/N}/\sqrt{N}

AB_{i,j}=\sum_{k=0}^{N-1} A_{i,k}B_{k,j}
= \sum_{k=0}^{N-1} e^{2\pi ik\sqrt{-1}/N} e^{-2\pi kj\sqrt{-1}/N}/N
=\sum_{k=0}^{N-1} e^{2\pi k(i-j)\sqrt{-1}/N}/N (2)
當(dāng)i=j時(shí),上式 = 1
當(dāng) i!=j時(shí), 上式 =0(證明略)

綜上,A確實(shí)是酉矩陣,則
A^{-1}=B
X=Bx 表達(dá)為一般形式,有:
對(duì)0<=k<1,k\in Z
X[k]=\sum_{n=0}^{N-1} e^{-2\pi kn\sqrt{-1}/N}/ \sqrt{N} *x[n]
如果記Y[k]=\sqrt{N}X[k]
我們便得到了標(biāo)準(zhǔn)的離散傅里葉變換表達(dá)式:
Y[k] = \sum_{n=0}^{N-1} e^{-2\pi kn\sqrt{-1}/N}*x[n]

綜上,離散傅里葉變換可以理解為把離散信號(hào):
x[t] (0<=t<N ,t\in Z) 展開為正交基 e_k[t] 的線性組合
之后, 其系數(shù)的求解公式(乘以一個(gè)常量因子)

至此,就得到了離散傅里葉變換DFT的一個(gè)純數(shù)學(xué)的非常簡(jiǎn)明的理解。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

  • 深入理解傅里葉變換Mar 12, 2017 這原本是我在知乎上對(duì)傅立葉變換、拉普拉斯變換、Z變換的聯(lián)系?為什么要進(jìn)...
    價(jià)值趨勢(shì)技術(shù)派閱讀 5,940評(píng)論 2 2
  • 概述 ??希爾伯特空間是一個(gè)完備的內(nèi)積空間,其標(biāo)準(zhǔn)正交函數(shù)系,直觀來看就是向量空間中基的延伸。其為基于任意正交系上...
    殉道者之花火閱讀 2,095評(píng)論 1 5
  • 網(wǎng)上關(guān)于從連續(xù)傅里葉變換推導(dǎo)出離散傅里葉變換公式的資料好像比較少,博主查閱了不少資料,總結(jié)出了一個(gè)推導(dǎo)的思路,現(xiàn)在...
    llooRice閱讀 58,502評(píng)論 1 8
  • 離散傅里葉變換 標(biāo)準(zhǔn)的 FFTs fft(a[, n, axis, norm]) 計(jì)算一維離散傅立葉變換。 iff...
    無賴宵小閱讀 4,897評(píng)論 0 1
  • ??傅里葉變換真是個(gè)磨人的小妖精,只要是搞理工的,就很難躲避它。我有一本買了很久的《信號(hào)與線性系統(tǒng)分析》,主要內(nèi)容...
    lqbii閱讀 4,960評(píng)論 0 1

友情鏈接更多精彩內(nèi)容