1. 從“聽音辨譜”到“萬物皆波”一個(gè)工程師眼中的傅里葉世界如果你玩過音樂軟件一定見過那種隨著節(jié)奏跳動(dòng)的頻譜柱狀圖。一首復(fù)雜的交響樂在軟件里被分解成不同頻率、不同強(qiáng)度的“柱子”低音鼓對應(yīng)低頻小提琴對應(yīng)高頻。這個(gè)將復(fù)雜聲音拆解成簡單正弦波的過程其背后的數(shù)學(xué)靈魂就是傅里葉分析。但它的威力遠(yuǎn)不止于此。在我十多年的信號處理、圖像分析和通信系統(tǒng)開發(fā)經(jīng)歷中傅里葉變換FT及其前身傅里葉級數(shù)FS是如同“螺絲刀”般的基礎(chǔ)且強(qiáng)大的工具。它解決的是一個(gè)根本性問題如何用“頻率”的語言去描述和理解一個(gè)“隨時(shí)間變化”的信號這不僅僅是數(shù)學(xué)家的游戲更是工程師將現(xiàn)實(shí)世界物理量聲音、圖像、電壓、震動(dòng)轉(zhuǎn)化為可計(jì)算、可處理、可設(shè)計(jì)的數(shù)字信息的橋梁。今天我們不堆砌公式而是從一個(gè)實(shí)踐者的角度聊聊傅里葉級數(shù)和傅里葉變換究竟在干什么為什么它們?nèi)绱酥匾约霸趯?shí)際項(xiàng)目中我們是如何“用”起來并避開那些教科書里不會寫的“坑”。簡單來說傅里葉級數(shù)專治“周期性”信號。它告訴我們?nèi)魏我粋€(gè)周期性的、滿足一定條件的復(fù)雜波形都可以看作是一系列頻率成整數(shù)倍關(guān)系的正弦波和余弦波或者復(fù)數(shù)指數(shù)波疊加而成。就像一道復(fù)雜的菜可以被分解為鹽、糖、醋等基本調(diào)料按不同比例的組合。而傅里葉變換則是傅里葉級數(shù)的“威力加強(qiáng)版”它將這個(gè)思想推廣到了非周期信號。它不再要求信號重復(fù)出現(xiàn)而是將任何有限能量的信號都看作是無限多個(gè)連續(xù)頻率的正弦波的疊加。如果說傅里葉級數(shù)給了我們一把分析重復(fù)節(jié)奏如心臟跳動(dòng)、交流電的解剖刀那么傅里葉變換則給了我們一臺可以分析任何一段聲音、任何一幅圖像、任何一段傳感器數(shù)據(jù)的CT掃描儀。理解這兩者是進(jìn)入數(shù)字信號處理、通信、音頻編解碼、圖像壓縮、甚至金融時(shí)間序列分析等眾多領(lǐng)域的門票。2. 傅里葉級數(shù)為周期信號“譜曲”當(dāng)我們面對一個(gè)周期性重復(fù)的信號時(shí)比如一個(gè)方波、一個(gè)三角波或者電網(wǎng)中標(biāo)準(zhǔn)的50Hz正弦交流電傅里葉級數(shù)是我們首選的解析工具。它的核心思想極具美感任何復(fù)雜的周期振動(dòng)都是由一系列最簡單的“純音”正弦/余弦波和諧地組合而成。這些純音的頻率是基礎(chǔ)頻率信號周期的倒數(shù)的整數(shù)倍。2.1 核心公式與物理意義傅里葉級數(shù)有兩種常見的表達(dá)形式。首先是三角形式f(t) a?/2 Σ [a? cos(nω?t) b? sin(nω?t)] 其中求和從 n1 到無窮。 這里ω? 2π/T是基波角頻率T是信號的周期。a?/2 代表信號的直流分量也就是整個(gè)周期內(nèi)的平均值。如果你測一個(gè)交流電壓信號這個(gè)值可能為零但如果你測的是一個(gè)疊加了直流偏置的脈搏信號這個(gè)值就包含了重要的靜態(tài)信息。a?和b? 這就是傅里葉系數(shù)是我們要計(jì)算的核心。a?衡量了信號與cos(nω?t)這個(gè)余弦波的“相似程度”b?則衡量了與正弦波的相似程度。它們的計(jì)算公式是信號在一個(gè)周期內(nèi)與對應(yīng)基函數(shù)的“內(nèi)積”即乘積的積分再除以歸一化因子。這本質(zhì)上是一種“投影”操作把信號投影到一組完備的正交基sin和cos函數(shù)系上看它在每個(gè)基方向上的分量有多大。另一種更簡潔、在工程中更常用的是復(fù)數(shù)形式指數(shù)形式f(t) Σ c? e^(j n ω? t) 求和從 n-∞ 到 ∞。 這里c?是復(fù)數(shù)系數(shù)。這個(gè)形式的優(yōu)勢在于數(shù)學(xué)處理極其方便一個(gè)公式統(tǒng)一了正弦和余弦并且c?直接包含了第 n 次諧波的幅度和相位信息|c?|是幅度∠c?是相位。同時(shí)由于e^(jωt)的微分、積分性質(zhì)極好它在求解微分方程、分析線性系統(tǒng)時(shí)具有無可比擬的優(yōu)勢。注意很多初學(xué)者會困惑于負(fù)頻率的出現(xiàn)。在復(fù)數(shù)形式中n 可以為負(fù)。這并非物理上存在“負(fù)頻率”的振動(dòng)而是一種數(shù)學(xué)上的便利表達(dá)。一個(gè)實(shí)信號我們實(shí)際測量的都是實(shí)信號的頻譜其正負(fù)頻率部分是共軛對稱的兩者共同決定了一個(gè)實(shí)際正弦波的幅度和相位。你可以理解為一對正負(fù)頻率的復(fù)數(shù)指數(shù)分量組合起來才對應(yīng)一個(gè)物理上可觀測的實(shí)正弦波。2.2 實(shí)戰(zhàn)中的“頻譜圖”與吉布斯現(xiàn)象在實(shí)際工程中我們很少手算這些系數(shù)。無論是使用 MATLAB、Python (NumPy/SciPy) 還是 LabVIEW都有現(xiàn)成的函數(shù)如fft 但注意 FFT 是計(jì)算離散傅里葉變換的快速算法對于周期信號其結(jié)果經(jīng)過適當(dāng)縮放和解釋就對應(yīng)傅里葉級數(shù)的系數(shù)。當(dāng)我們計(jì)算出系數(shù)后通常會繪制兩張圖幅度譜和相位譜。幅度譜 以頻率nω?為橫坐標(biāo)以|c?|或sqrt(a?2b?2)為縱坐標(biāo)的柱狀圖。它一目了然地告訴我們這個(gè)復(fù)雜周期信號里包含哪些頻率的成分以及各自的“音量”有多大。例如一個(gè)理想的方波其頻譜只包含奇數(shù)次諧波1倍 3倍 5倍...基頻且幅度與諧波次數(shù)成反比。相位譜 以頻率為橫坐標(biāo)以∠c?為縱坐標(biāo)。它告訴我們每個(gè)頻率成分的“起始時(shí)間”或波形對齊關(guān)系。在有些應(yīng)用中如通信系統(tǒng)的均衡、音頻的相位校正相位信息至關(guān)重要。一個(gè)經(jīng)典的、必須理解的工程現(xiàn)象是吉布斯現(xiàn)象。當(dāng)你用有限項(xiàng)比如前N項(xiàng)的傅里葉級數(shù)去逼近一個(gè)具有間斷點(diǎn)比如方波的跳變沿的信號時(shí)在間斷點(diǎn)附近會出現(xiàn)持續(xù)的過沖和振蕩并且即使N趨于無窮這個(gè)過沖的峰值也不會消失大約為跳變值的9%左右。這不是計(jì)算錯(cuò)誤而是傅里葉級數(shù)在間斷點(diǎn)收斂特性的體現(xiàn)。實(shí)操心得在濾波器設(shè)計(jì)中吉布斯現(xiàn)象直接導(dǎo)致了設(shè)計(jì)出的濾波器在通帶和阻帶邊界會出現(xiàn)紋波。為了抑制它工程師們發(fā)明了各種窗函數(shù)如漢寧窗、漢明窗、布萊克曼窗。加窗的本質(zhì)就是用一種平滑過渡的方式“柔和地”截?cái)酂o限長的理想濾波器脈沖響應(yīng)雖然犧牲了一些過渡帶陡峭度但換來了紋波的減小。這是理論聯(lián)系實(shí)際的一個(gè)絕佳案例。3. 傅里葉變換打開非周期信號的頻率之門傅里葉級數(shù)很美但它有個(gè)硬性要求信號必須周期重復(fù)。現(xiàn)實(shí)世界中大量的信號是非周期的比如一段語音、一個(gè)雷達(dá)脈沖、一幅圖像的一行像素亮度變化。傅里葉變換應(yīng)運(yùn)而生它可以看作是傅里葉級數(shù)在周期T趨于無窮大時(shí)的極限形式。3.1 從離散譜到連續(xù)譜概念的飛躍在傅里葉級數(shù)中頻譜是離散的只有基頻整數(shù)倍的那些頻率點(diǎn)上有值。當(dāng)周期T → ∞時(shí)基頻ω? 2π/T → 0 離散的譜線間隔無限縮小最終連成一片變成了連續(xù)的頻譜。同時(shí)求和Σ也變成了積分∫。傅里葉變換對定義如下 正變換從時(shí)域到頻域F(ω) ∫_{-∞}^{∞} f(t) e^(-jωt) dt逆變換從頻域到時(shí)域f(t) (1/2π) ∫_{-∞}^{∞} F(ω) e^(jωt) dω這里的F(ω)是一個(gè)連續(xù)復(fù)數(shù)函數(shù)稱為頻譜密度函數(shù)。它的模|F(ω)|稱為幅度譜密度表示信號能量在頻率上的分布密度它的輻角∠F(ω)稱為相位譜。為什么是“密度”這是關(guān)鍵。對于能量有限的非周期信號其在任何一個(gè)單一頻率點(diǎn)上的能量是零因?yàn)檫B續(xù)。|F(ω)|2在頻域上對頻率的積分等于信號在時(shí)域上的總能量帕斯瓦爾定理。所以|F(ω)|2描述的是單位頻率帶寬內(nèi)包含的信號能量即能量譜密度。這在分析噪聲如白噪聲具有平坦的功率譜密度、評估信號帶寬時(shí)極其重要。3.2 典型信號的傅里葉變換與物理直覺理解幾個(gè)典型信號的變換能極大增強(qiáng)你對頻域的直覺。單頻正弦波cos(ω? t) 它的傅里葉變換是在ω?和-ω?處的兩個(gè)沖激函數(shù)狄拉克δ函數(shù)。這完美印證了一個(gè)純正余弦波其所有能量都集中在正負(fù)基頻這兩個(gè)點(diǎn)上。這也是頻譜分析儀測量單一頻率信號幅度的原理基礎(chǔ)。矩形脈沖門函數(shù) 一個(gè)持續(xù)時(shí)間為τ、幅度為A的矩形脈沖。它的頻譜是Aτ · Sa(ωτ/2) 其中Sa(x) sin(x)/x是抽樣函數(shù)。這個(gè)頻譜是連續(xù)的、呈sinc函數(shù)形狀。重要結(jié)論時(shí)域上越“窄”的脈沖τ越小其頻譜主瓣越“寬”帶寬越大。時(shí)域和頻域的這種“展寬”關(guān)系是反比的這是不確定性原理在信號處理中的體現(xiàn)。在設(shè)計(jì)雷達(dá)脈沖或通信符號時(shí)這是一個(gè)核心權(quán)衡。高斯脈沖 高斯函數(shù)的傅里葉變換仍然是高斯函數(shù)。這是唯一一個(gè)在時(shí)域和頻域形狀相同的函數(shù)除了尺度因子因此在需要同時(shí)要求時(shí)域和頻域集中性的場合如超寬帶通信高斯脈沖及其衍生波形備受青睞。實(shí)操心得在仿真或處理實(shí)際信號時(shí)我們永遠(yuǎn)無法處理無限長的信號。我們處理的都是一段有限時(shí)間的信號。這相當(dāng)于用一個(gè)矩形窗去截取原始信號。時(shí)域的截?cái)喑司匦未霸陬l域等價(jià)于原始信號的頻譜與一個(gè)sinc函數(shù)矩形窗的頻譜進(jìn)行卷積。這會導(dǎo)致頻譜泄漏——原本集中在某個(gè)頻率的能量會“泄漏”到旁邊的頻率區(qū)間使得頻譜看起來變模糊、出現(xiàn)虛假的旁瓣。這也是為什么在做頻譜分析前常常要對數(shù)據(jù)加窗非矩形窗的原因加窗就是為了抑制頻譜泄漏的副作用雖然代價(jià)是降低了頻率分辨率。4. 離散傅里葉變換與快速算法數(shù)字世界的基石我們生活在數(shù)字時(shí)代計(jì)算機(jī)只能處理離散的、有限長的數(shù)據(jù)。對應(yīng)于連續(xù)傅里葉變換我們有了離散傅里葉變換對應(yīng)于無限長序列我們處理的是有限長序列。而讓這一切變得實(shí)用的是快速傅里葉變換算法。4.1 離散傅里葉變換對連續(xù)世界的采樣與近似DFT 的定義是針對一個(gè)長度為 N 的離散序列x[n](n0,1,...,N-1)X[k] Σ_{n0}^{N-1} x[n] · e^{-j (2π/N) k n} k0,1,...,N-1。 逆變換為x[n] (1/N) Σ_{k0}^{N-1} X[k] · e^{j (2π/N) k n}。你可以這樣理解 DFT時(shí)域采樣 我們對連續(xù)時(shí)間信號f(t)以間隔Ts進(jìn)行采樣得到離散序列x[n] f(n·Ts)。頻域采樣 理論上這個(gè)離散時(shí)間序列的頻譜是周期性的這是采樣定理的結(jié)論。DFT 所做的就是在一個(gè)周期[0, Fs)內(nèi)Fs1/Ts是采樣率等間隔地取 N 個(gè)點(diǎn)計(jì)算這 N 個(gè)頻率點(diǎn)上的頻譜值X[k]。X[k]對應(yīng)的是數(shù)字頻率ω_k (2πk)/(N) 或模擬頻率f_k k·Fs/N。關(guān)鍵參數(shù)關(guān)系頻率分辨率Δf Fs / N。 這意味著為了區(qū)分兩個(gè)靠得很近的頻率成分你需要增加采樣點(diǎn)數(shù) N 或降低采樣率 Fs在滿足奈奎斯特定理的前提下。N 決定了你頻域譜線的精細(xì)程度。奈奎斯特頻率Fs/2。 這是 DFT 能無混疊地表示的最高頻率。任何高于Fs/2的信號頻率成分都會以“混疊”的形式折疊到[0, Fs/2)區(qū)間內(nèi)造成失真。因此采樣前必須用抗混疊濾波器將高于Fs/2的成分濾除。4.2 快速傅里葉變換讓實(shí)時(shí)處理成為可能DFT 的直接計(jì)算復(fù)雜度是O(N2) 當(dāng) N 很大時(shí)比如 1024 2048計(jì)算量是災(zāi)難性的。庫利和圖基發(fā)明的 FFT 算法通過巧妙的分解將大 N 的 DFT 分解為小 N 的 DFT 組合將復(fù)雜度降低到了O(N log? N)。當(dāng) N1024 時(shí)速度提升超過 100 倍。正是 FFT 的出現(xiàn)才使得實(shí)時(shí)頻譜分析、數(shù)字濾波、正交頻分復(fù)用等技術(shù)在嵌入式系統(tǒng)和消費(fèi)電子產(chǎn)品中得以實(shí)現(xiàn)。在使用 FFT例如 Python 的numpy.fft.fft MATLAB 的fft LabVIEW 的 FFT VI時(shí)有幾個(gè)必須注意的細(xì)節(jié)幅度校正 FFT 輸出的原始結(jié)果X[k]通常需要縮放才能得到真實(shí)的物理幅度。對于周期信號且整周期采樣幅度 |X[k]| * 2 / N(對于 k≠0 和 k≠N/2 的頻率線)。直流分量 (k0) 為|X[0]| / N。對于非周期信號/隨機(jī)信號更常用的是計(jì)算功率譜密度例如使用周期圖法PSD |X[k]|2 / (Fs * N)。這能給出每赫茲的功率。頻率向量的生成 橫坐標(biāo)頻率軸必須正確生成。對于單邊譜通常只顯示 0 到 Fs/2頻率向量為f k * Fs / N 其中k 0, 1, ..., N//2。補(bǔ)零的影響 對數(shù)據(jù)序列后面補(bǔ)零再做 FFT可以增加頻譜的顯示點(diǎn)數(shù)插值讓曲線看起來更光滑但并不能提高真正的頻率分辨率。頻率分辨率只由原始數(shù)據(jù)長度T N * Ts決定即Δf 1/T。補(bǔ)零只是對已有的頻譜進(jìn)行了插值沒有增加新的信息。踩坑實(shí)錄在一次電機(jī)振動(dòng)分析項(xiàng)目中我們采集了轉(zhuǎn)速傳感器的脈沖信號。直接做 FFT 后頻譜圖非常雜亂預(yù)期的轉(zhuǎn)頻及其倍頻被淹沒在噪聲里。問題出在傳感器信號不是穩(wěn)態(tài)的電機(jī)在啟動(dòng)階段轉(zhuǎn)速有輕微變化。我們對非平穩(wěn)信號做了全局 FFT相當(dāng)于假設(shè)整個(gè)時(shí)間段內(nèi)頻率成分不變這必然導(dǎo)致頻譜模糊。解決方案是采用短時(shí)傅里葉變換將長信號分幀對每一幀做 FFT從而觀察頻譜隨時(shí)間的變化成功捕捉到了轉(zhuǎn)頻的軌跡。這個(gè)教訓(xùn)是FFT 是強(qiáng)大的但它假設(shè)信號是平穩(wěn)的。對于時(shí)變信號需要更高級的工具如 STFT、小波變換。5. 二維傅里葉變換圖像處理與計(jì)算機(jī)視覺的眼睛傅里葉變換不僅限于一維時(shí)間信號。對于圖像這種二維空間信號二維傅里葉變換是其頻率分析的基礎(chǔ)。一幅圖像可以看作亮度在 x, y 兩個(gè)方向上的變化。5.1 空間頻率與圖像頻譜二維傅里葉變換的公式是雙重積分連續(xù)或雙重求和離散F(u, v) ∫∫ f(x, y) e^{-j2π(uxvy)} dx dy逆變換f(x, y) ∫∫ F(u, v) e^{j2π(uxvy)} du dv這里(u, v)是空間頻率變量單位可以是“周期每毫米”或“周期每像素”。低頻成分 對應(yīng)圖像中變化緩慢的部分如大面積的天空、墻面、背景。在頻譜圖上低頻集中在中心區(qū)域。高頻成分 對應(yīng)圖像中快速變化的部分如邊緣、紋理、細(xì)節(jié)、噪聲。在頻譜圖上高頻分布在四周。對一幅圖像做二維 FFT如numpy.fft.fft2 得到的是一個(gè)復(fù)數(shù)矩陣。通常我們可視化其對數(shù)幅度譜log(1 |F|)來增強(qiáng)對比。你會發(fā)現(xiàn)頻譜圖中心最亮低頻能量最強(qiáng)。如果圖像中有明顯的周期性紋理如布料、柵欄頻譜上會出現(xiàn)對應(yīng)的亮點(diǎn)對由于實(shí)信號頻譜的共軛對稱性。如果圖像有強(qiáng)邊緣如一個(gè)矩形物體其頻譜在垂直于邊緣的方向上會出現(xiàn)明亮的線條因?yàn)檫吘夘愃朴陔A躍函數(shù)其頻譜衰減較慢。5.2 頻域?yàn)V波的實(shí)際應(yīng)用在頻域?qū)D像進(jìn)行操作最經(jīng)典的應(yīng)用就是濾波。其流程是原圖 → 2D FFT → 頻譜中心化fftshift → 與濾波函數(shù)頻域掩模相乘 → 反中心化 → 2D IFFT → 輸出圖像。低通濾波 保留低頻抑制高頻。可以用于圖像去噪和圖像平滑模糊。常用的濾波函數(shù)有理想低通濾波器ILPF、巴特沃斯低通濾波器BLPF、高斯低通濾波器GLPF。理想低通會產(chǎn)生振鈴效應(yīng)類似吉布斯現(xiàn)象高斯低通效果最平滑。實(shí)操步驟生成一個(gè)與頻譜圖同樣大小的掩模H(u, v)。對于高斯低通H(u, v) exp(-D(u,v)2 / (2*D02)) 其中D(u,v)是點(diǎn)到頻譜中心的距離D0是截止頻率。將中心化后的頻譜F_shift與H逐點(diǎn)相乘得到G_shift F_shift * H 然后進(jìn)行反變換。高通濾波 保留高頻抑制低頻。可以用于圖像銳化和邊緣檢測。高通濾波后圖像中平滑區(qū)域變暗邊緣區(qū)域被突出。帶阻/帶通濾波 在頻譜上“挖掉”或“保留”特定區(qū)域。最著名的應(yīng)用是去除周期性噪聲。例如掃描圖像中常見的摩爾紋、傳感器引起的固定模式噪聲在頻譜上會表現(xiàn)為遠(yuǎn)離中心的一對對稱亮點(diǎn)。我們可以設(shè)計(jì)一個(gè)帶阻濾波器比如在兩個(gè)亮點(diǎn)位置置零的小區(qū)域?qū)⑵錇V除再反變換回去噪聲就能被有效抑制。經(jīng)驗(yàn)技巧在 Python 中使用cv2.dft()或numpy.fft進(jìn)行圖像濾波時(shí)一定要注意數(shù)據(jù)類型和值域。傅里葉變換后的頻譜值動(dòng)態(tài)范圍極大直接顯示可能是一片黑。一定要用np.log()進(jìn)行壓縮顯示。進(jìn)行濾波乘法時(shí)確保濾波掩模是浮點(diǎn)類型。逆變換后需要用np.abs()取模或cv2.magnitude()計(jì)算幅度并可能需要進(jìn)行歸一化 (cv2.normalize) 轉(zhuǎn)換回uint8格式才能正確顯示。一個(gè)常見的錯(cuò)誤是忘記對頻譜進(jìn)行中心化 (np.fft.fftshift)導(dǎo)致濾波操作在錯(cuò)誤的位置進(jìn)行。6. 超越基礎(chǔ)傅里葉分析在現(xiàn)代工程中的延伸掌握了傅里葉級數(shù)、傅里葉變換和 DFT/FFT你已經(jīng)擁有了強(qiáng)大的武器。但在更復(fù)雜的場景下我們需要它的“升級版”。6.1 短時(shí)傅里葉變換處理非平穩(wěn)信號傳統(tǒng) FT 和 FFT 給出的是信號在整個(gè)時(shí)間范圍內(nèi)的全局頻率信息它假設(shè)信號的頻率成分是穩(wěn)定不變的。但對于頻率隨時(shí)間變化的信號如語音、音樂、雷達(dá)回波全局頻譜是無效的。STFT 的基本思想很簡單加窗、分幀、逐幀 FFT。用一個(gè)移動(dòng)的窗函數(shù)如漢明窗截取信號的一小段。對這一小段信號做 FFT得到該時(shí)刻附近的“局部頻譜”。將窗函數(shù)沿時(shí)間軸滑動(dòng)重復(fù)上述過程。將所有局部頻譜排列起來就得到了一個(gè)二維矩陣一個(gè)維度是頻率另一個(gè)維度是時(shí)間。將其可視化就是聲譜圖——語音識別和音樂分析的核心工具。STFT 的缺點(diǎn)是時(shí)間-頻率分辨率矛盾窗長決定了分辨率。窗越長頻率分辨率越高能區(qū)分更近的頻率但時(shí)間分辨率越低無法定位頻率變化的精確時(shí)刻窗越短則反之。這是一個(gè)無法兼得的權(quán)衡。6.2 傅里葉變換在通信與調(diào)制中的應(yīng)用這是傅里葉變換的“主場”之一。正交頻分復(fù)用 4G/5G 和 Wi-Fi 的核心技術(shù)。它將高速數(shù)據(jù)流分割成許多低速子流用不同的正交子載波正弦波同時(shí)傳輸。這里的“正交”就是指子載波在符號周期內(nèi)滿足整數(shù)倍周期關(guān)系在接收端通過相關(guān)積分本質(zhì)是傅里葉變換或等效操作可以完美分離它們。OFDM 的發(fā)射和接收機(jī)核心就是 IFFT 和 FFT 運(yùn)算。單邊帶調(diào)制 為了節(jié)省帶寬我們希望只傳輸調(diào)幅信號的一個(gè)邊帶。利用希爾伯特變換與傅里葉變換密切相關(guān)可以構(gòu)造出解析信號從而濾除一個(gè)邊帶。信道估計(jì)與均衡 在頻域多徑信道可以建模為一個(gè)頻域?yàn)V波器H(f)。通過發(fā)送已知的導(dǎo)頻信號在接收端進(jìn)行 FFT對比發(fā)送和接收的頻域數(shù)據(jù)可以估計(jì)出H(f)。隨后在頻域進(jìn)行簡單的除法運(yùn)算Y(f)/H(f)就能實(shí)現(xiàn)信道均衡補(bǔ)償多徑效應(yīng)帶來的失真。這比時(shí)域的均衡算法往往更高效。6.3 拉普拉斯變換與 Z 變換更廣義的視角傅里葉變換要求信號絕對可積能量有限這對于很多工程信號如階躍信號、指數(shù)增長信號是不滿足的。拉普拉斯變換通過引入一個(gè)衰減因子e^{-σt} 將不滿足絕對可積的信號“壓”下去使其變換存在。它將傅里葉變換的虛軸jω?cái)U(kuò)展到了整個(gè)復(fù)平面s σ jω。拉普拉斯變換是分析線性時(shí)不變系統(tǒng)尤其是電路、控制系統(tǒng)穩(wěn)定性和瞬態(tài)響應(yīng)的終極工具。系統(tǒng)函數(shù)H(s)的極點(diǎn)位置直接決定了系統(tǒng)的穩(wěn)定性。Z 變換則是針對離散時(shí)間系統(tǒng)的拉普拉斯變換將復(fù)平面s映射到z平面。它是數(shù)字濾波器設(shè)計(jì)和分析的基石。一個(gè)數(shù)字濾波器的傳遞函數(shù)H(z) 其零極點(diǎn)分布決定了濾波器的頻率響應(yīng)特性低通、高通、帶通等。設(shè)計(jì)一個(gè) IIR 濾波器本質(zhì)上就是在z平面上配置零極點(diǎn)。傅里葉變換可以看作是拉普拉斯變換在虛軸上的特例 (s jω) 也是 Z 變換在單位圓上的特例 (z e^{jω})。它們共同構(gòu)成了信號與系統(tǒng)分析的完整體系。回顧這趟從周期到非周期、從連續(xù)到離散、從一維到二維的旅程傅里葉分析的精髓在于“變換視角”。它教會我們同一個(gè)物理現(xiàn)象在時(shí)域可能雜亂無章在頻域卻可能井然有序。這種切換視角的能力是解決復(fù)雜工程問題的關(guān)鍵。我個(gè)人的體會是不要被公式嚇倒多從物理圖像和實(shí)際應(yīng)用案例去理解它。當(dāng)你第一次用 FFT 從嘈雜的傳感器數(shù)據(jù)中分離出設(shè)備的故障特征頻率時(shí)當(dāng)你第一次用頻域?yàn)V波讓模糊的老照片重現(xiàn)清晰輪廓時(shí)你會真正感受到這個(gè)兩百多年前的數(shù)學(xué)思想所蘊(yùn)含的磅礴力量。它不僅僅是工具箱里的一個(gè)算法更是一種理解世界波動(dòng)本質(zhì)的語言。