1. 從“周期”二字說起為什么離散周期信號(hào)的頻譜值得單獨(dú)拎出來講很多人學(xué)數(shù)字信號(hào)處理學(xué)到“離散周期信號(hào)的頻譜”這一節(jié)時(shí)第一反應(yīng)是——公式好多但好像跟前面學(xué)的離散時(shí)間傅里葉變換DTFT差不多如果你也這么想過那說明你還沒被實(shí)際工程問題毒打過。我當(dāng)年做音頻降噪和振動(dòng)信號(hào)分析的時(shí)候就是在這兩個(gè)概念上栽了跟頭明明理論上頻譜應(yīng)該是幾根干凈的譜線實(shí)際用計(jì)算機(jī)一算出來的卻是一片糊掉的“山包”。后來才明白問題就出在離散和周期這兩個(gè)詞同時(shí)出現(xiàn)時(shí)頻譜的形態(tài)會(huì)發(fā)生本質(zhì)變化——它不再是連續(xù)譜而是離散譜線而且這些譜線的位置、幅度、相位都有非常嚴(yán)格的數(shù)學(xué)約束。這篇文章要聊的就是離散周期信號(hào)頻譜的完整邏輯鏈它為什么是離散的線譜、跟DFS離散傅里葉級(jí)數(shù)是什么關(guān)系、實(shí)際做頻譜分析時(shí)混疊和泄露這兩個(gè)“老朋友”是怎么毀掉你的結(jié)果的以及怎么用工程手段把它們按在地上。適合已經(jīng)學(xué)過信號(hào)與系統(tǒng)基礎(chǔ)、正在做信號(hào)處理相關(guān)項(xiàng)目、或者被頻譜分析結(jié)果搞得一頭霧水的朋友。我會(huì)盡量少堆公式多用“人話”和實(shí)際踩坑經(jīng)驗(yàn)把這件事講透。先給一個(gè)最核心的結(jié)論后面所有內(nèi)容都圍繞它展開一個(gè)域時(shí)域的離散化和周期化必然導(dǎo)致另一個(gè)域頻域的周期化和離散化。這句話是理解本文所有內(nèi)容的鑰匙。時(shí)域離散頻域就周期時(shí)域周期頻域就離散。當(dāng)信號(hào)同時(shí)滿足離散和周期兩個(gè)條件時(shí)它的頻譜就是周期性的離散譜線——既離散又周期這就是DFS存在的根本原因。2. 離散周期信號(hào)頻譜的核心邏輯拆解2.1 為什么時(shí)域周期會(huì)導(dǎo)致頻域離散先把這個(gè)最反直覺的點(diǎn)講清楚。連續(xù)時(shí)間里的周期信號(hào)比如一個(gè)純正弦波它的頻譜是兩根沖激線位于正負(fù)頻率處。這本身就是“時(shí)域周期→頻域離散”的體現(xiàn)。但到了離散時(shí)間領(lǐng)域事情變得更微妙。你可以這樣理解時(shí)域周期意味著信號(hào)每隔N個(gè)點(diǎn)就重復(fù)一次這種重復(fù)性在頻域里表現(xiàn)為“只有某些特定頻率分量能夠存在”。具體來說只有那些在N點(diǎn)周期內(nèi)完成整數(shù)個(gè)周期的頻率分量才能跟信號(hào)的周期性“合拍”。其他頻率分量在周期延拓時(shí)會(huì)相互抵消。所以頻域被“篩”成了離散的譜線相鄰譜線的間隔是2π/N數(shù)字角頻率或者fs/N模擬頻率fs為采樣率。這里有個(gè)容易混淆的點(diǎn)離散周期信號(hào)的頻譜是離散的但離散非周期信號(hào)的頻譜是連續(xù)的DTFT。很多人把這兩個(gè)搞混結(jié)果在推導(dǎo)時(shí)怎么都對(duì)不上。記住這個(gè)對(duì)應(yīng)關(guān)系時(shí)域特性頻域特性變換名稱連續(xù)、周期離散、非周期傅里葉級(jí)數(shù)FS連續(xù)、非周期連續(xù)、非周期傅里葉變換FT離散、非周期連續(xù)、周期DTFT離散、周期離散、周期DFS這張表建議你直接背下來。我當(dāng)年考試時(shí)就是靠它理清了所有變換的關(guān)系。DFS的“周期”體現(xiàn)在頻域上是因?yàn)闀r(shí)域離散DFS的“離散”體現(xiàn)在頻域上是因?yàn)闀r(shí)域周期。兩個(gè)條件缺一不可。2.2 DFS的數(shù)學(xué)形式與物理意義DFS的正變換和反變換長這樣正變換X[k] Σ(n0到N-1) x[n] * exp(-j2πk*n/N)反變換x[n] (1/N) * Σ(k0到N-1) X[k] * exp(j2πk*n/N)看起來跟DFT一模一樣對(duì)吧確實(shí)DFS和DFT在計(jì)算形式上完全相同。區(qū)別在于解釋DFS是把一個(gè)周期信號(hào)分解成N個(gè)諧波分量X[k]是第k次諧波的復(fù)幅度而DFT是把一個(gè)有限長序列做變換隱含了對(duì)序列做周期延拓的假設(shè)。這個(gè)區(qū)別在實(shí)際工程中極其重要。當(dāng)你用計(jì)算機(jī)對(duì)一個(gè)信號(hào)做FFT時(shí)你實(shí)際上是在計(jì)算DFT而DFT隱含的周期延拓假設(shè)正是頻譜泄露的根源。如果你不理解DFS的周期延拓本質(zhì)你就永遠(yuǎn)搞不懂為什么截?cái)嘈盘?hào)會(huì)導(dǎo)致頻譜糊掉。X[k]的物理意義是在頻率k*fs/N處信號(hào)分量的復(fù)幅度。它的模|X[k]|是幅度譜角度∠X[k]是相位譜。由于DFS的周期性X[k] X[kN]所以通常只看k0到N-1這一個(gè)周期。2.3 離散周期信號(hào)頻譜的三個(gè)關(guān)鍵特征第一個(gè)特征是譜線間隔固定。相鄰譜線的頻率間隔是fs/N這個(gè)間隔由信號(hào)的周期N和采樣率fs共同決定。N越大譜線越密fs越高譜線間隔越大。這個(gè)關(guān)系在做頻率分辨率分析時(shí)是核心。第二個(gè)特征是只有N根獨(dú)立譜線。雖然頻域是無限周期的但一個(gè)周期內(nèi)只有N個(gè)獨(dú)立值。這意味著你無法從DFS中分辨出頻率間隔小于fs/N的兩個(gè)分量。這就是頻率分辨率的物理極限。第三個(gè)特征是幅度與周期延拓的匹配。如果信號(hào)的周期不是整數(shù)個(gè)采樣點(diǎn)或者你截取的N不是信號(hào)周期的整數(shù)倍那么周期延拓時(shí)在邊界處會(huì)出現(xiàn)跳變這個(gè)跳變?cè)陬l域就表現(xiàn)為頻譜泄露。這是后面要重點(diǎn)講的內(nèi)容。3. 混疊采樣定理背后的隱形殺手3.1 混疊的本質(zhì)頻率的“身份錯(cuò)位”混疊Aliasing這個(gè)詞聽起來很玄其實(shí)本質(zhì)很簡單當(dāng)采樣率不夠高時(shí)高頻信號(hào)會(huì)“偽裝”成低頻信號(hào)。這不是比喻是數(shù)學(xué)上的嚴(yán)格等價(jià)。舉個(gè)例子。假設(shè)采樣率fs1000Hz那么奈奎斯特頻率是500Hz?，F(xiàn)在有一個(gè)600Hz的正弦波按照采樣定理它已經(jīng)超過了奈奎斯特頻率應(yīng)該被抗混疊濾波器濾掉。但如果你沒濾采樣后的序列會(huì)是什么樣600Hz的信號(hào)采樣周期Ts1/1000秒。每個(gè)采樣點(diǎn)的相位增量是2π*600/1000 1.2π。由于正弦的周期性1.2π的相位增量等價(jià)于1.2π - 2π -0.8π也就是相當(dāng)于-400Hz的相位增量。所以采樣后的序列看起來就像一個(gè)400Hz的信號(hào)。600Hz“偽裝”成了400Hz這就是混疊。更一般地頻率f的信號(hào)采樣后會(huì)被映射到f - k*fsk為整數(shù)取落在[-fs/2, fs/2]范圍內(nèi)的那個(gè)值。所以混疊的本質(zhì)是頻率的周期性延拓——采樣使得頻譜以fs為周期重復(fù)如果原信號(hào)帶寬超過fs/2相鄰周期就會(huì)重疊重疊部分就產(chǎn)生了混疊。3.2 抗混疊濾波器的設(shè)計(jì)要點(diǎn)避免混疊的唯一工程手段就是在ADC之前加抗混疊濾波器低通濾波器把高于fs/2的頻率分量濾掉。這個(gè)濾波器必須是模擬濾波器因?yàn)閿?shù)字濾波在采樣之后已經(jīng)來不及了。設(shè)計(jì)抗混疊濾波器時(shí)有幾個(gè)關(guān)鍵參數(shù)截止頻率通常設(shè)為fs/2但實(shí)際要留余量一般設(shè)為0.4fs到0.45fs給濾波器的過渡帶留空間。過渡帶寬度從通帶邊緣到阻帶邊緣的頻率范圍。過渡帶越窄濾波器階數(shù)越高成本越高。阻帶衰減通常要求至少-60dB到-80dB確保高頻分量被充分抑制。我做過一個(gè)振動(dòng)監(jiān)測項(xiàng)目最初為了省成本用了截止頻率0.5fs的濾波器結(jié)果頻譜在接近奈奎斯特頻率處總是有莫名其妙的峰。后來把截止頻率降到0.4fs問題消失。這個(gè)教訓(xùn)是永遠(yuǎn)不要貼著奈奎斯特頻率設(shè)計(jì)抗混疊濾波器一定要留保護(hù)帶。注意抗混疊濾波器必須在采樣之前且必須是模擬濾波器。數(shù)字濾波器無法防止混疊因?yàn)榛殳B在采樣瞬間就已經(jīng)發(fā)生了。3.3 過采樣如何幫助對(duì)抗混疊過采樣Oversampling是另一個(gè)實(shí)用手段。如果你用遠(yuǎn)高于信號(hào)最高頻率的采樣率采樣比如信號(hào)最高頻率1kHz你用100kHz采樣那么即使抗混疊濾波器的過渡帶很寬也能保證混疊分量落在有用頻帶之外。過采樣的好處是降低了對(duì)模擬濾波器過渡帶的苛刻要求。你可以用一個(gè)簡單的RC低通濾波器截止頻率設(shè)在10kHz雖然過渡帶很寬但因?yàn)椴蓸勇适?00kHz奈奎斯特頻率是50kHz混疊分量會(huì)落在50kHz附近遠(yuǎn)離你關(guān)心的0-1kHz頻帶。后續(xù)可以用數(shù)字濾波器進(jìn)一步處理。代價(jià)是數(shù)據(jù)量增大后續(xù)處理負(fù)擔(dān)加重。所以實(shí)際工程中需要在濾波器復(fù)雜度和數(shù)據(jù)量之間做權(quán)衡。我的經(jīng)驗(yàn)是如果信號(hào)帶寬遠(yuǎn)小于奈奎斯特頻率過采樣是性價(jià)比很高的方案如果信號(hào)本身就接近奈奎斯特頻率那就老老實(shí)實(shí)設(shè)計(jì)高階模擬濾波器。4. 頻譜泄露截?cái)鄮淼谋厝徽`差4.1 泄露的根源周期延拓不連續(xù)泄露Leakage是DFT分析中最常見的“背鍋俠”。很多人看到頻譜糊了就說“泄露了”但未必清楚泄露到底怎么來的。回到DFS的假設(shè)DFS要求信號(hào)是周期的且你分析的N個(gè)點(diǎn)正好是一個(gè)周期。但實(shí)際中你拿到的信號(hào)往往是無限長非周期信號(hào)的一段截?cái)唷．?dāng)你對(duì)這N個(gè)點(diǎn)做DFT時(shí)隱含的假設(shè)是這N個(gè)點(diǎn)周期延拓后構(gòu)成一個(gè)周期信號(hào)。如果截?cái)嗵幍男盘?hào)值不連續(xù)——也就是x[N-1]和x[0]接不上——那么周期延拓后就會(huì)在邊界處產(chǎn)生跳變。這個(gè)跳變?cè)陬l域表現(xiàn)為什么一個(gè)跳變包含豐富的頻率分量這些分量會(huì)“泄露”到原本不應(yīng)該有能量的頻點(diǎn)上。所以泄露的本質(zhì)是截?cái)鄬?dǎo)致的周期延拓不連續(xù)使得原本集中在某個(gè)頻點(diǎn)的能量擴(kuò)散到了整個(gè)頻帶。用一個(gè)直觀的例子假設(shè)你有一個(gè)10Hz的正弦波采樣率1000Hz你截取100個(gè)點(diǎn)0.1秒。10Hz的信號(hào)在0.1秒內(nèi)正好完成1個(gè)周期所以截?cái)嗵庍B續(xù)DFT結(jié)果是一根干凈的譜線。但如果你截取105個(gè)點(diǎn)10Hz信號(hào)在0.105秒內(nèi)完成1.05個(gè)周期截?cái)嗵幉贿B續(xù)DFT結(jié)果就會(huì)在10Hz附近出現(xiàn)一系列逐漸衰減的譜線——這就是泄露。4.2 窗函數(shù)用“軟截?cái)唷睋Q“低泄露”解決泄露的標(biāo)準(zhǔn)手段是加窗。窗函數(shù)的作用是讓截?cái)嗵幍男盘?hào)值平滑過渡到零從而減小周期延拓時(shí)的跳變。常用的窗函數(shù)有窗函數(shù)主瓣寬度旁瓣衰減適用場景矩形窗最窄最差-13dB信號(hào)周期正好整數(shù)倍截?cái)酀h寧窗中等較好-31dB通用頻譜分析漢明窗中等好-43dB需要低旁瓣布萊克曼窗寬很好-58dB強(qiáng)干擾下的弱信號(hào)檢測凱澤窗可調(diào)可調(diào)需要靈活權(quán)衡選擇窗函數(shù)的核心權(quán)衡是主瓣寬度 vs 旁瓣衰減。主瓣越窄頻率分辨率越高旁瓣越低泄露越小。但兩者不可兼得你只能根據(jù)具體需求選。我個(gè)人的經(jīng)驗(yàn)法則如果做通用頻譜分析漢寧窗是默認(rèn)選擇它的旁瓣衰減和主瓣寬度平衡得最好。如果信號(hào)里有強(qiáng)弱分量并存弱信號(hào)容易被強(qiáng)信號(hào)的旁瓣淹沒這時(shí)候用布萊克曼窗或凱澤窗。如果做精確的頻率測量且能保證整周期采樣矩形窗反而最好因?yàn)樗靼曜钫?。注意加窗?huì)改變信號(hào)的幅度需要在頻域做幅度修正。不同窗函數(shù)的相干增益不同修正系數(shù)也不同。漢寧窗的相干增益是0.5所以幅度要乘以2漢明窗是0.54乘以1/0.54≈1.85。這個(gè)細(xì)節(jié)很多人忽略導(dǎo)致幅度測量偏差。4.3 整周期采樣從源頭消滅泄露比加窗更徹底的方法是整周期采樣。如果你能讓采樣點(diǎn)數(shù)N正好覆蓋信號(hào)周期的整數(shù)倍那么截?cái)嗵幾匀贿B續(xù)泄露為零。這在旋轉(zhuǎn)機(jī)械振動(dòng)分析中特別有用因?yàn)檗D(zhuǎn)速信號(hào)可以提供周期參考。實(shí)現(xiàn)整周期采樣的關(guān)鍵是同步采樣用轉(zhuǎn)速脈沖觸發(fā)ADC確保每次采樣的起始點(diǎn)都在同一個(gè)相位。這樣即使轉(zhuǎn)速有波動(dòng)也能保證整周期截?cái)?。我做過一個(gè)電機(jī)故障診斷項(xiàng)目用了光電編碼器做同步觸發(fā)頻譜干凈得跟仿真一樣旁瓣完全消失。但整周期采樣需要額外的硬件支持不是所有場景都能用。如果做不到退而求其次用窗函數(shù)也能把泄露壓到可接受的水平。5. 實(shí)操從仿真到實(shí)測的完整頻譜分析流程5.1 仿真驗(yàn)證用Python復(fù)現(xiàn)混疊和泄露理論講再多不如自己跑一遍代碼。下面這段Python代碼演示了混疊和泄露的產(chǎn)生過程你可以直接復(fù)制運(yùn)行。import numpy as np import matplotlib.pyplot as plt # 參數(shù)設(shè)置 fs 1000 # 采樣率 T 1.0 # 信號(hào)時(shí)長 t np.arange(0, T, 1/fs) # 場景1正常采樣信號(hào)頻率100Hz f1 100 x1 np.sin(2*np.pi*f1*t) X1 np.fft.fft(x1) freq1 np.fft.fftfreq(len(x1), 1/fs) # 場景2混疊信號(hào)頻率600Hz超過奈奎斯特頻率500Hz f2 600 x2 np.sin(2*np.pi*f2*t) X2 np.fft.fft(x2) freq2 np.fft.fftfreq(len(x2), 1/fs) # 場景3泄露截取非整數(shù)周期 f3 105 # 105Hz在1秒內(nèi)完成105個(gè)周期但截取0.1秒只有10.5個(gè)周期 t3 np.arange(0, 0.1, 1/fs) x3 np.sin(2*np.pi*f3*t3) X3 np.fft.fft(x3) freq3 np.fft.fftfreq(len(x3), 1/fs) # 繪圖 fig, axes plt.subplots(3, 2, figsize(12, 10)) axes[0,0].plot(t[:200], x1[:200]) axes[0,0].set_title(正常采樣: 100Hz正弦波) axes[0,1].plot(freq1[:len(freq1)//2], np.abs(X1[:len(X1)//2])) axes[0,1].set_title(頻譜: 干凈的單峰) axes[1,0].plot(t[:200], x2[:200]) axes[1,0].set_title(混疊: 600Hz正弦波 (fs1000Hz)) axes[1,1].plot(freq2[:len(freq2)//2], np.abs(X2[:len(X2)//2])) axes[1,1].set_title(頻譜: 混疊到400Hz) axes[2,0].plot(t3[:100], x3[:100]) axes[2,0].set_title(泄露: 105Hz正弦波, 截?cái)?.1秒) axes[2,1].plot(freq3[:len(freq3)//2], np.abs(X3[:len(X3)//2])) axes[2,1].set_title(頻譜: 能量泄露到鄰近頻點(diǎn)) plt.tight_layout() plt.show()運(yùn)行這段代碼你會(huì)看到三個(gè)明顯的現(xiàn)象正常采樣的頻譜是一根干凈的譜線600Hz的信號(hào)混疊到了400Hz105Hz信號(hào)截?cái)嗪箢l譜出現(xiàn)了明顯的旁瓣。這三個(gè)現(xiàn)象分別對(duì)應(yīng)了本文的核心內(nèi)容。5.2 實(shí)測中的參數(shù)選擇與陷阱實(shí)際做頻譜分析時(shí)參數(shù)選擇直接決定結(jié)果質(zhì)量。以下是我總結(jié)的參數(shù)選擇流程第一步確定信號(hào)最高頻率fmax。這決定了采樣率的下限。根據(jù)采樣定理fs 2fmax。但實(shí)際中建議fs 2.5fmax到5*fmax給抗混疊濾波器留過渡帶。第二步確定頻率分辨率Δf。Δf fs/N其中N是采樣點(diǎn)數(shù)。如果你需要分辨間隔為1Hz的兩個(gè)頻率分量那么Δf必須小于1Hz即N fs/1 fs。所以采樣點(diǎn)數(shù)N fs/Δf。第三步確定窗函數(shù)。如果信號(hào)周期已知且能整周期采樣用矩形窗否則用漢寧窗或漢明窗。如果強(qiáng)信號(hào)旁邊有弱信號(hào)用布萊克曼窗。第四步確定平均次數(shù)。單次FFT的頻譜波動(dòng)很大方差等于均值的平方需要多次平均來平滑。平均次數(shù)越多頻譜越穩(wěn)定但測量時(shí)間越長。一般至少平均10次要求高的話平均100次。這里有個(gè)常見的陷阱采樣點(diǎn)數(shù)N和FFT點(diǎn)數(shù)不是一回事。你可以采集N個(gè)點(diǎn)然后補(bǔ)零到M點(diǎn)做FFTM N這樣頻譜看起來更“光滑”但頻率分辨率并沒有提高——分辨率仍然由N決定補(bǔ)零只是插值。很多人被這個(gè)假象騙了以為補(bǔ)零能提高分辨率。記住頻率分辨率只取決于實(shí)際采樣時(shí)長T N/fs補(bǔ)零不增加信息量。5.3 從頻譜反推信號(hào)特征的實(shí)戰(zhàn)技巧拿到一個(gè)頻譜怎么判斷它是干凈的還是有問題的我總結(jié)了一個(gè)快速檢查清單看底噪如果底噪平坦且低說明沒有混疊和泄露如果底噪在某些頻段隆起可能是混疊?？醋V線寬度如果譜線很窄接近一根線說明泄露小如果譜線很寬像山包說明泄露嚴(yán)重?？磁园耆绻靼陜蓚?cè)有規(guī)律的衰減旁瓣那是窗函數(shù)的特征如果旁瓣不規(guī)則可能是信號(hào)本身有調(diào)制?？搭l率位置如果峰值頻率跟預(yù)期不符檢查是否混疊?；殳B頻率 |f - k*fs|取最小值。我遇到過一個(gè)典型案例某設(shè)備振動(dòng)信號(hào)在頻譜上顯示主頻為25Hz但設(shè)備轉(zhuǎn)速明明是3000rpm50Hz。后來發(fā)現(xiàn)采樣率設(shè)成了100Hz奈奎斯特頻率50Hz50Hz的信號(hào)正好在奈奎斯特頻率上混疊后變成了0Hz附近的直流分量而25Hz是50Hz的次諧波。把采樣率提高到500Hz后主頻正確顯示為50Hz。這個(gè)坑讓我明白采樣率的選擇必須基于對(duì)信號(hào)頻率范圍的充分了解不能拍腦袋。6. 常見問題與排查技巧實(shí)錄6.1 頻譜分析中的典型問題速查表現(xiàn)象可能原因排查方法解決方案頻譜出現(xiàn)莫名低頻峰混疊檢查信號(hào)最高頻率是否超過fs/2提高采樣率或加抗混疊濾波器譜線周圍有大量旁瓣泄露檢查截?cái)嗵幨欠襁B續(xù)加窗或整周期采樣頻譜整體模糊分辨率不足檢查采樣時(shí)長T增加采樣點(diǎn)數(shù)N幅度測量偏小窗函數(shù)增益未修正檢查窗函數(shù)類型乘以相干增益修正系數(shù)頻譜波動(dòng)大未平均檢查平均次數(shù)增加平均次數(shù)峰值頻率偏移柵欄效應(yīng)檢查峰值是否落在譜線上補(bǔ)零或調(diào)整采樣率6.2 獨(dú)家避坑經(jīng)驗(yàn)?zāi)切┪臋n里不會(huì)寫的細(xì)節(jié)第一個(gè)坑抗混疊濾波器的截止頻率不是越高越好。我見過有人把截止頻率設(shè)在0.49fs覺得這樣能保留最多信號(hào)。但實(shí)際上濾波器的過渡帶需要空間截止頻率太靠近fs/2會(huì)導(dǎo)致過渡帶內(nèi)的分量混疊回來。建議截止頻率不超過0.4fs。第二個(gè)坑加窗之后不要忘記做重疊處理。如果信號(hào)是分段處理的加窗會(huì)導(dǎo)致邊界處的信息丟失。標(biāo)準(zhǔn)做法是50%重疊即每次移動(dòng)N/2個(gè)點(diǎn)。這樣雖然計(jì)算量翻倍但能保證不丟信息。第三個(gè)坑FFT點(diǎn)數(shù)不是越多越好。FFT點(diǎn)數(shù)超過采樣點(diǎn)數(shù)N時(shí)補(bǔ)零雖然讓頻譜更光滑但不會(huì)提高分辨率反而增加計(jì)算量。FFT點(diǎn)數(shù)等于N或略大于N取2的冪次即可。第四個(gè)坑相位譜比幅度譜更容易受泄露影響。如果你需要精確的相位信息泄露會(huì)導(dǎo)致相位劇烈波動(dòng)。這時(shí)候整周期采樣幾乎是唯一的選擇。第五個(gè)坑不要忽略窗函數(shù)的主瓣寬度對(duì)分辨率的影響。矩形窗主瓣寬度是2fs/N漢寧窗是4fs/N。加漢寧窗后頻率分辨率實(shí)際上降低了一半。如果你需要分辨間隔為Δf的兩個(gè)分量用漢寧窗時(shí)采樣時(shí)長需要加倍。6.3 一個(gè)真實(shí)項(xiàng)目的排查記錄某次做電機(jī)軸承故障診斷頻譜上在2000Hz附近出現(xiàn)了一系列等間隔的峰間隔約50Hz。按照軸承故障頻率計(jì)算不應(yīng)該有這個(gè)特征。排查過程如下第一步檢查采樣率。采樣率10kHz奈奎斯特頻率5kHz2000Hz遠(yuǎn)低于奈奎斯特頻率排除混疊。第二步檢查泄露。信號(hào)截取長度0.5秒分辨率2Hz。2000Hz處的峰間隔50Hz遠(yuǎn)大于分辨率不是泄露導(dǎo)致的模糊。第三步檢查信號(hào)本身。用示波器看原始信號(hào)發(fā)現(xiàn)時(shí)域波形上有明顯的周期性沖擊沖擊間隔20ms對(duì)應(yīng)50Hz。原來這是電源工頻干擾耦合到了信號(hào)里。2000Hz處的峰是50Hz的40次諧波。第四步解決方案。在信號(hào)調(diào)理電路里加了50Hz陷波器問題解決。這個(gè)案例的教訓(xùn)是頻譜上的異常特征先回到時(shí)域找原因。很多問題在時(shí)域一目了然在頻域反而繞暈了。7. 離散周期信號(hào)頻譜的工程應(yīng)用場景7.1 旋轉(zhuǎn)機(jī)械振動(dòng)分析旋轉(zhuǎn)機(jī)械的振動(dòng)信號(hào)是典型的離散周期信號(hào)。轉(zhuǎn)速?zèng)Q定了基頻各階諧波對(duì)應(yīng)不同的機(jī)械故障。比如不平衡對(duì)應(yīng)1倍頻不對(duì)中對(duì)應(yīng)2倍頻軸承故障對(duì)應(yīng)高頻沖擊。在這類應(yīng)用中整周期采樣是標(biāo)配。用轉(zhuǎn)速脈沖觸發(fā)采樣保證每次分析的信號(hào)正好是整數(shù)個(gè)旋轉(zhuǎn)周期。這樣頻譜干凈故障特征清晰。我做過的一個(gè)案例中整周期采樣后1倍頻的幅度測量重復(fù)性從±15%提高到±2%故障判斷準(zhǔn)確率大幅提升。7.2 電力系統(tǒng)諧波分析電力系統(tǒng)的電壓電流信號(hào)是50Hz或60Hz的周期信號(hào)但含有大量諧波。用DFS分析各次諧波的幅度和相位是電能質(zhì)量分析的核心。這里的關(guān)鍵是同步采樣。電網(wǎng)頻率會(huì)在50Hz附近波動(dòng)如果采樣率固定會(huì)導(dǎo)致非整周期截?cái)喈a(chǎn)生泄露。解決方案是用鎖相環(huán)跟蹤電網(wǎng)頻率動(dòng)態(tài)調(diào)整采樣率保證每個(gè)周期采樣點(diǎn)數(shù)固定。這樣諧波分析精度可以做到0.1%以內(nèi)。7.3 通信系統(tǒng)中的OFDM信號(hào)分析OFDM信號(hào)在時(shí)域是多個(gè)子載波疊加看起來像噪聲但在頻域是離散的譜線。每個(gè)子載波對(duì)應(yīng)一個(gè)頻點(diǎn)子載波間隔固定。分析OFDM信號(hào)的頻譜需要保證分析帶寬覆蓋所有子載波且頻率分辨率小于子載波間隔。這里混疊和泄露的影響特別大?；殳B會(huì)導(dǎo)致子載波間干擾泄露會(huì)導(dǎo)致誤碼率上升。所以O(shè)FDM系統(tǒng)對(duì)采樣率和同步的要求極高。實(shí)際系統(tǒng)中前導(dǎo)序列就是用來做同步和信道估計(jì)的。8. 寫在最后幾個(gè)讓我少走彎路的習(xí)慣做信號(hào)分析這些年我養(yǎng)成了幾個(gè)習(xí)慣幫我避免了很多低級(jí)錯(cuò)誤。第一個(gè)習(xí)慣是先仿真再實(shí)測。任何新的分析方法先在Python或MATLAB里用合成信號(hào)驗(yàn)證一遍確認(rèn)理解了原理和參數(shù)影響再上實(shí)測數(shù)據(jù)。這樣出了問題能快速定位是方法問題還是數(shù)據(jù)問題。第二個(gè)習(xí)慣是同時(shí)看時(shí)域和頻域。頻譜上看到異常馬上回到時(shí)域找對(duì)應(yīng)特征。很多時(shí)候時(shí)域波形上的一個(gè)小毛刺在頻域就是一大片。反過來時(shí)域看起來正常的信號(hào)頻域可能藏著問題。第三個(gè)習(xí)慣是記錄每次分析的參數(shù)。采樣率、采樣點(diǎn)數(shù)、窗函數(shù)、平均次數(shù)、濾波器設(shè)置這些參數(shù)全部記下來。下次遇到類似信號(hào)直接復(fù)用參數(shù)效率高很多。而且出了問題能回溯是哪一步導(dǎo)致的。第四個(gè)習(xí)慣是對(duì)頻譜保持懷疑?？吹礁蓛舻淖V線先想想是不是整周期采樣了看到模糊的頻譜先想想是不是泄露了看到奇怪的峰先想想是不是混疊了。把混疊和泄露這兩個(gè)可能性過一遍大部分問題都能解釋。離散周期信號(hào)的頻譜分析理論不復(fù)雜但工程細(xì)節(jié)極多?；殳B和泄露這兩個(gè)問題本質(zhì)上都是“域之間的不匹配”造成的混疊是采樣率與信號(hào)帶寬不匹配泄露是截?cái)嚅L度與信號(hào)周期不匹配。理解了這一點(diǎn)你就知道所有解決方案的方向——要么讓條件匹配提高采樣率、整周期采樣要么讓不匹配的后果可控抗混疊濾波器、窗函數(shù)。這個(gè)思路比死記公式有用得多。