處理實(shí)驗(yàn)從代碼到頻譜:采樣定理與重建誤差的底層邏輯)
簡(jiǎn)介這份資源是南京郵電大學(xué)「信號(hào)分析與處理實(shí)驗(yàn)」課程的完整實(shí)驗(yàn)報(bào)告面向正在修讀數(shù)字信號(hào)處理、信號(hào)與系統(tǒng)相關(guān)課程的高校學(xué)生以及需要借助 MATLAB 完成信號(hào)處理實(shí)驗(yàn)與課程設(shè)計(jì)的學(xué)習(xí)者。報(bào)告圍繞信號(hào)的產(chǎn)生和運(yùn)算、連續(xù)時(shí)間信號(hào)的頻域分析、時(shí)域采樣與重建、離散傅里葉變換及應(yīng)用等核心實(shí)驗(yàn)展開(kāi)并延伸至連續(xù)時(shí)間系統(tǒng)頻域與復(fù)頻域分析、數(shù)字濾波器設(shè)計(jì)及信號(hào)處理綜合應(yīng)用配有 MATLAB 程序示例與圖形結(jié)果便于對(duì)照理解時(shí)域和頻域特性。資源包共 1 個(gè) docx 文件約 1.02MB內(nèi)容為實(shí)驗(yàn)?zāi)康?、?shí)驗(yàn)內(nèi)容、過(guò)程與結(jié)果分析的完整記錄可直接參考實(shí)驗(yàn)思路、程序編寫與數(shù)據(jù)解析方法。目前已有 952 人學(xué)習(xí)適合作為實(shí)驗(yàn)預(yù)習(xí)、報(bào)告撰寫與 DSP 入門復(fù)習(xí)的參考材料。1. 信號(hào)處理實(shí)驗(yàn)的“黑匣子”為什么你照著敲代碼卻看不懂頻譜很多人第一次做信號(hào)分析與處理實(shí)驗(yàn)時(shí)都有過(guò)這種體驗(yàn)MATLAB 代碼一行行敲進(jìn)去stem、plot、fft跑得飛起圖也畫出來(lái)了可一旦被問(wèn)“這個(gè)頻譜為什么長(zhǎng)這樣”“采樣間隔改一下為什么就混疊了”腦子立刻空白。問(wèn)題不在于不會(huì)用 MATLAB而在于把信號(hào)處理當(dāng)成了“調(diào)包畫圖”的玄學(xué)——代碼是抄的結(jié)果是懵的。這份實(shí)驗(yàn)報(bào)告合集實(shí)驗(yàn) 1 到 4真正的價(jià)值不是給你幾段能跑的腳本而是把“信號(hào)從時(shí)域到頻域、從連續(xù)到離散、從采樣到重建”這條鏈路完整地?cái)傞_(kāi)。它逼著你親手產(chǎn)生單位階躍、復(fù)指數(shù)、正弦序列親手做卷積、算奇偶部、畫雙邊譜、觀察吉布斯過(guò)沖、驗(yàn)證奈奎斯特間隔。每一個(gè)現(xiàn)象背后都對(duì)應(yīng)一個(gè)可推導(dǎo)的機(jī)制。搞懂這些機(jī)制你才算真正跨進(jìn)了 DSP 的門而不是停留在“MATLAB 繪圖員”的層面。下面就從信號(hào)怎么被“造”出來(lái)開(kāi)始一層層拆到采樣定理和重建誤差的底層邏輯。2. 信號(hào)的產(chǎn)生與運(yùn)算從序列定義到卷積的底層邏輯2.1 離散序列不是“連續(xù)曲線的采樣點(diǎn)”而是獨(dú)立定義的映射實(shí)驗(yàn)一里最先碰到的坑是把離散序列當(dāng)成連續(xù)信號(hào)來(lái)理解。比如單位階躍序列u[n]代碼里寫的是x[zeros(1,10) 1 ones(1,20)]配合n-10:20。這行代碼的含義是在索引n-10到n-1上取 0在n0上取 1在n1到n20上取 1。它和連續(xù)階躍函數(shù)u(t)最大的區(qū)別在于離散序列只在整數(shù)點(diǎn)上有定義點(diǎn)與點(diǎn)之間沒(méi)有“過(guò)渡”也沒(méi)有“無(wú)窮小的跳變”。這解釋了一個(gè)常見(jiàn)困惑為什么u[n7]的代碼是[zeros(1,10M) 1 ones(1,20-M)]因?yàn)樽笠?7 位意味著原本在n0處的跳變移到了n-7所以前面補(bǔ)的零從 10 個(gè)變成 17 個(gè)后面剩下的 1 從 20 個(gè)變成 13 個(gè)。移位在離散域就是索引的平移不涉及任何插值或重采樣。理解這一點(diǎn)后面做卷積、做 DFT 的循環(huán)移位才不會(huì)翻車。復(fù)指數(shù)序列x[n]1.5e^{(0.2j3)n}更值得細(xì)看。代碼里c0.23*1i然后xK*exp(c*n)最后分別畫real(x)和imag(x)。這里藏著一個(gè)關(guān)鍵機(jī)制復(fù)指數(shù)序列的實(shí)部和虛部都是衰減或增長(zhǎng)振蕩的。實(shí)部對(duì)應(yīng)1.5e^{0.2n}\cos(3n)虛部對(duì)應(yīng)1.5e^{0.2n}\sin(3n)。很多新手看到實(shí)部圖里幅度在漲以為是“信號(hào)不穩(wěn)定”其實(shí)只是e^{0.2n}這個(gè)包絡(luò)在起作用。把復(fù)指數(shù)拆成包絡(luò)和振蕩兩部分是理解 Z 變換和傅里葉變換的起點(diǎn)。2.2 奇偶分解為什么fliplr能拆出對(duì)稱分量實(shí)驗(yàn)一里有一個(gè)看似簡(jiǎn)單但極其重要的操作對(duì)x[n]{1,2,4,-5,-2}n-2到n2做奇偶分解。代碼只有三行核心y fliplr(x); % 時(shí)間反轉(zhuǎn)x[-n] xe 0.5*(xy); % 偶部 xo 0.5*(x-y); % 奇部這里fliplr做的是左右翻轉(zhuǎn)對(duì)應(yīng)時(shí)域反轉(zhuǎn)x[-n]。但要注意這個(gè)操作成立的前提是序列的索引關(guān)于n0對(duì)稱。如果序列從n0開(kāi)始直接fliplr會(huì)得到錯(cuò)誤的“反轉(zhuǎn)”結(jié)果因?yàn)閤[0]翻轉(zhuǎn)后應(yīng)該還在n0而不是跑到末尾。實(shí)驗(yàn)里特意選了n-2:2這個(gè)對(duì)稱區(qū)間就是為了讓fliplr和真正的時(shí)域反轉(zhuǎn)等價(jià)。偶部xe[n](x[n]x[-n])/2滿足xe[n]xe[-n]奇部xo[n](x[n]-x[-n])/2滿足xo[n]-xo[-n]。這個(gè)分解在信號(hào)處理里的意義是任何序列都可以拆成對(duì)稱部分和反對(duì)稱部分而傅里葉變換的實(shí)部對(duì)應(yīng)偶對(duì)稱、虛部對(duì)應(yīng)奇對(duì)稱。后面做 DFT 時(shí)如果輸入是實(shí)序列頻譜的共軛對(duì)稱性就是從這里來(lái)的。很多人背“實(shí)序列頻譜共軛對(duì)稱”這個(gè)結(jié)論卻不知道它和奇偶分解是同一件事。2.3 線性卷積conv背后到底在算什么實(shí)驗(yàn)一里x[n](-1)^n0≤n≤14和h[n]cos(0.25πn)0≤n≤7做卷積代碼直接調(diào)用yconv(x,h)結(jié)果長(zhǎng)度是158-122對(duì)應(yīng)n0:21。這里有兩個(gè)必須搞清楚的機(jī)制。第一卷積的物理含義是“翻轉(zhuǎn)、平移、相乘、求和”。conv函數(shù)內(nèi)部做的就是這件事把h翻轉(zhuǎn)然后沿著x滑動(dòng)每個(gè)位置做點(diǎn)積。對(duì)于長(zhǎng)度分別為L(zhǎng)和M的序列輸出長(zhǎng)度是LM-1。這個(gè)公式不是隨便定的而是因?yàn)榛瑒?dòng)窗口從完全重疊到完全分離一共經(jīng)歷了LM-1個(gè)位置。第二conv默認(rèn)從n0開(kāi)始計(jì)算它不關(guān)心你原來(lái)的索引是0:14還是-5:9。如果你輸入的序列索引不是從 0 開(kāi)始conv的輸出索引需要你自己平移。比如x的索引是-2:2h的索引是0:3那么卷積結(jié)果的起始索引是-20-2終點(diǎn)索引是235。實(shí)驗(yàn)里兩個(gè)序列都從 0 開(kāi)始所以輸出索引直接是0:21不需要額外平移。這個(gè)細(xì)節(jié)在后續(xù)做系統(tǒng)響應(yīng)、做濾波器時(shí)域分析時(shí)經(jīng)常被忽略導(dǎo)致畫出來(lái)的圖“對(duì)不上”。連續(xù)正弦信號(hào)x_a(t)cos(5000πt)4sin(2000πt)sin(3000πt)的繪制則提醒另一件事連續(xù)信號(hào)在 MATLAB 里是用極密的采樣點(diǎn)“假裝”的。代碼里t0:0.5e-6:6e-3采樣間隔 0.5 微秒對(duì)應(yīng)采樣率 2 MHz。這個(gè)采樣率遠(yuǎn)高于信號(hào)最高頻率分量5000π rad/s 對(duì)應(yīng) 2500 Hz2000π 和 3000π 的乘積產(chǎn)生 1000 Hz 和 5000 Hz 分量所以畫出來(lái)的波形是可信的。如果采樣率不夠畫出來(lái)的“連續(xù)波形”本身就是混疊的產(chǎn)物后面做頻域分析會(huì)全盤皆錯(cuò)。3. 連續(xù)時(shí)間信號(hào)的頻域分析從傅里葉級(jí)數(shù)到吉布斯過(guò)沖3.1 周期矩形脈沖的雙邊譜為什么 T 增大譜線變密實(shí)驗(yàn)二的核心是周期矩形脈沖信號(hào)的傅里葉級(jí)數(shù)。代碼P2_1畫的是A1, τ0.1s, T0.5s的雙邊幅度譜和相位譜。這里的關(guān)鍵參數(shù)是占空比τ/T和基頻Ω?2π/T。雙邊譜的譜線出現(xiàn)在kΩ?處幅度正比于sinc(kΩ?τ/2)。當(dāng)T從 0.5s 增大到 1s 和 2s 時(shí)基頻Ω?從4π減小到2π和π譜線間隔變小所以譜變得更“密”。同時(shí)因?yàn)閏?Aτ/T直流分量隨T增大而減小整體幅度譜的值也下降。這個(gè)現(xiàn)象背后是傅里葉級(jí)數(shù)到傅里葉變換的過(guò)渡周期趨于無(wú)窮大時(shí)離散譜線間隔趨于零譜線連成一片就變成了連續(xù)頻譜密度。實(shí)驗(yàn)里讓你比較T1s和T2s的圖其實(shí)是在用肉眼觀察這個(gè)極限過(guò)程。相位譜同樣值得注意。對(duì)于實(shí)偶對(duì)稱的矩形脈沖c_k是實(shí)數(shù)相位只有 0 和 π 兩種取值對(duì)應(yīng)正負(fù)號(hào)。如果脈沖不是關(guān)于原點(diǎn)對(duì)稱的相位就會(huì)變成連續(xù)的斜線對(duì)應(yīng)時(shí)移。這個(gè)“時(shí)移對(duì)應(yīng)相位線性變化”的規(guī)律在后續(xù)做 DFT 的相位分析時(shí)反復(fù)出現(xiàn)。3.2 諧波疊加與吉布斯現(xiàn)象為什么過(guò)沖不會(huì)消失P2_2用有限項(xiàng)傅里葉級(jí)數(shù)合成矩形脈沖N取 15、25、50?,F(xiàn)象是N越大合成波形越接近矩形但在不連續(xù)點(diǎn)附近始終存在過(guò)沖而且過(guò)沖的峰值不隨N增大而減小只是越來(lái)越靠近跳變點(diǎn)。這就是吉布斯現(xiàn)象。實(shí)驗(yàn)里P2_4用代碼算出了過(guò)沖最大值function p2_4(N) t-2:10^(-4):2; A1; tau1; T2; Omega_0 2*pi/T; c0A*tau/T; f_Nc0*ones(1,length(t)); for k1:1:N f_N f_N 2*A*tau/T*sinc(k*Omega_0*tau/2/pi)*cos(k*Omega_0*t); end overshoot max(f_N); fprintf(The maximum overshoot is %.4f\n, overshoot); end運(yùn)行結(jié)果是1.0896也就是約 8.96% 的過(guò)沖。這個(gè)數(shù)字不是偶然的它等于(2/π)∫?^π (sin t / t) dt - 1是一個(gè)只與跳變幅度有關(guān)的常數(shù)。無(wú)論N取多大只要是用有限項(xiàng)傅里葉級(jí)數(shù)去逼近有跳變的函數(shù)這個(gè)過(guò)沖就存在。要消除它只能換用其他基函數(shù)比如小波或者對(duì)級(jí)數(shù)做求和加速比如 Fejér 求和。實(shí)驗(yàn)里讓你親手算出這個(gè)數(shù)比背十遍“吉布斯現(xiàn)象”定義都管用。3.3 男女聲頻譜差異基音頻率和共振峰的物理含義實(shí)驗(yàn)二最后一部分用 FFT 分析男女聲錄音的頻譜。現(xiàn)象是男聲波形基音頻率低、振幅變化大頻譜能量集中在低頻女聲基音頻率高、波形細(xì)致頻譜能量分布更均勻。這里涉及兩個(gè)概念基音頻率pitch和共振峰formant?;纛l率由聲帶振動(dòng)決定男性通常在 100-150 Hz女性在 200-300 Hz。在頻譜上表現(xiàn)為一系列間隔為基音頻率的諧波。共振峰由聲道形狀決定表現(xiàn)為頻譜包絡(luò)上的峰值。男聲聲道更長(zhǎng)共振峰頻率更低女聲聲道更短共振峰頻率更高。實(shí)驗(yàn)里讓你“查閱資料回答男聲轉(zhuǎn)女聲的方法”答案里提到的基頻調(diào)整、聲道長(zhǎng)度調(diào)整、諧波調(diào)整本質(zhì)上都是在改這兩個(gè)東西。用 MATLAB 做 FFT 只是第一步真正難的是從頻譜里把基音和共振峰分離出來(lái)——這通常需要倒譜分析或線性預(yù)測(cè)編碼實(shí)驗(yàn)報(bào)告里沒(méi)有展開(kāi)但方向已經(jīng)指出來(lái)了。4. 時(shí)域采樣與重建奈奎斯特間隔不是“經(jīng)驗(yàn)值”而是硬約束4.1 采樣定理的驗(yàn)證為什么 Ts2 一定會(huì)混疊實(shí)驗(yàn)三用升余弦脈沖信號(hào)做采樣和頻譜分析。提示里給了一個(gè)關(guān)鍵推導(dǎo)升余弦脈沖的頻譜大部分集中在[0, 2π/τ]設(shè)截止頻率Ω_m2π/τ代入?yún)?shù)得Ω_m2所以奈奎斯特間隔T_s1/(2f_m)π/2。當(dāng)T_sπ/2時(shí)采樣頻譜剛好不混疊當(dāng)T_s2時(shí)T_s π/2采樣頻譜發(fā)生混疊。這個(gè)推導(dǎo)里有一個(gè)容易忽略的點(diǎn)Ω_m2是角頻率對(duì)應(yīng)的實(shí)際頻率是f_mΩ_m/(2π)1/π。奈奎斯特率是2f_m2/π奈奎斯特間隔是1/(2f_m)π/2。實(shí)驗(yàn)里直接給T_sπ/2是因?yàn)樯嘞颐}沖的頻譜不是嚴(yán)格帶限的Ω_m只是一個(gè)“大部分能量集中”的邊界。如果信號(hào)嚴(yán)格帶限T_sπ/2剛好滿足如果信號(hào)有高頻拖尾即使T_sπ/2也可能有輕微混疊。實(shí)驗(yàn)里判斷“未發(fā)生混疊”是基于頻譜圖的主瓣是否重疊不是嚴(yán)格的數(shù)學(xué)證明。當(dāng)T_s2時(shí)采樣率f_s1/T_s0.5遠(yuǎn)低于奈奎斯特率2/π≈0.637混疊不可避免。頻譜圖上會(huì)看到原本在高頻的副本搬移到了低頻和主瓣重疊。這個(gè)現(xiàn)象在音頻處理里就是“ aliasing ”聽(tīng)起來(lái)像金屬聲或鳥(niǎo)叫。實(shí)驗(yàn)里讓你親眼看到混疊后的頻譜比記公式深刻得多。4.2 信號(hào)重建理想低通濾波器為什么在現(xiàn)實(shí)中做不到P3_2做的是采樣后重建。理想重建是用一個(gè)截止頻率為Ω_cΩ_m的理想低通濾波器對(duì)采樣信號(hào)做濾波恢復(fù)出原始連續(xù)信號(hào)。代碼里用sinc插值實(shí)現(xiàn)。當(dāng)T_sπ/2時(shí)重建誤差很小當(dāng)T_s2且Ω_cΩ_m時(shí)重建誤差大大增加。誤差增大的原因有兩層第一采樣間隔增大導(dǎo)致頻譜混疊高頻分量被“折疊”到低頻這部分信息在重建時(shí)無(wú)法恢復(fù)第二低通濾波器的截止頻率設(shè)為Ω_m但混疊后的信號(hào)里包含了原本高于Ω_m的頻率分量這些分量被濾波器保留了一部分表現(xiàn)為噪聲。實(shí)驗(yàn)報(bào)告里寫“高頻分量噪聲成分增加”指的就是這個(gè)。這里有一個(gè)常見(jiàn)的誤解很多人以為只要采樣率滿足奈奎斯特定理重建就完美。實(shí)際上理想低通濾波器在時(shí)域是sinc函數(shù)它無(wú)限長(zhǎng)、非因果物理不可實(shí)現(xiàn)。實(shí)際重建只能用近似濾波器會(huì)引入相位失真和幅度失真。實(shí)驗(yàn)里用sinc插值是在“模擬”理想重建不是真的實(shí)現(xiàn)了它。理解這個(gè)區(qū)別才能明白為什么實(shí)際音頻系統(tǒng)里過(guò)采樣和抗混疊濾波器那么重要。4.3 采樣間隔與誤差的定量關(guān)系實(shí)驗(yàn)里沒(méi)有給出誤差的定量公式但可以從機(jī)制上推導(dǎo)。對(duì)于升余弦脈沖采樣后的重建誤差主要來(lái)自兩部分混疊誤差和截?cái)嗾`差?;殳B誤差隨T_s增大而增大大致正比于T_s超過(guò)奈奎斯特間隔的部分。截?cái)嗾`差來(lái)自sinc插值的有限長(zhǎng)度在實(shí)驗(yàn)里因?yàn)閠的范圍足夠大可以忽略。如果要在實(shí)驗(yàn)里定量驗(yàn)證可以計(jì)算重建信號(hào)與原始信號(hào)的均方誤差error mean((xa_recon - xa_original).^2);當(dāng)T_sπ/2時(shí)error應(yīng)該很小當(dāng)T_s2時(shí)error會(huì)顯著增大。這個(gè)計(jì)算在實(shí)驗(yàn)報(bào)告里沒(méi)有寫但它是驗(yàn)證采樣定理最直接的方法。很多新手做完實(shí)驗(yàn)只畫圖不量化結(jié)果對(duì)“誤差增大”只有模糊的視覺(jué)印象沒(méi)有數(shù)字概念。加上這一行整個(gè)實(shí)驗(yàn)的嚴(yán)謹(jǐn)性就上了一個(gè)臺(tái)階。5. 避坑與排查信號(hào)處理實(shí)驗(yàn)里最容易翻車的五個(gè)地方5.1 現(xiàn)象conv結(jié)果長(zhǎng)度不對(duì)畫圖時(shí)橫坐標(biāo)對(duì)不上原因conv輸出長(zhǎng)度是LM-1但很多人直接用輸入序列的索引作為輸出索引。如果輸入索引不是從 0 開(kāi)始輸出索引需要平移。解決先確定輸出起始索引n_start n_x_start n_h_start然后生成n_out n_start : n_start length(y) - 1。畫圖時(shí)用stem(n_out, y)不要用stem(y)。5.2 現(xiàn)象fft得到的頻譜和理論譜對(duì)不上幅度差一個(gè)系數(shù)原因MATLAB 的fft不除以N而傅里葉級(jí)數(shù)系數(shù)c_k通常要除以N或T。另外fft的輸出順序是[0, 1, ..., N-1]對(duì)應(yīng)頻率[0, f_s/N, ..., f_s(N-1)/N]不是雙邊譜的順序。解決用fftshift把零頻移到中心用abs(fft(x))/N得到幅度譜。如果要畫雙邊譜頻率軸用(-N/2:N/2-1)*f_s/N。5.3 現(xiàn)象采樣后頻譜圖看起來(lái)“沒(méi)混疊”但重建誤差很大原因頻譜圖的分辨率不夠或者只看了幅度譜沒(méi)看相位譜。混疊可能發(fā)生在高頻段幅度很小肉眼看不出來(lái)但重建時(shí)誤差會(huì)累積。解決把頻譜圖放大到高頻段用對(duì)數(shù)坐標(biāo)20*log10(abs(X))觀察。同時(shí)計(jì)算重建誤差的均方值不要只靠肉眼判斷。5.4 現(xiàn)象吉布斯過(guò)沖的數(shù)值每次運(yùn)行都不一樣原因t的步長(zhǎng)10^(-4)可能沒(méi)有剛好落在不連續(xù)點(diǎn)上導(dǎo)致max(f_N)取到的峰值位置有微小偏移。另外N越大過(guò)沖峰值越靠近跳變點(diǎn)但t的離散化可能錯(cuò)過(guò)真正的峰值。解決用更細(xì)的步長(zhǎng)比如10^(-5)或者在跳變點(diǎn)附近加密采樣。理論上過(guò)沖極限是1.0896實(shí)際計(jì)算值應(yīng)該穩(wěn)定在這個(gè)數(shù)附近。5.5 現(xiàn)象男女聲頻譜對(duì)比時(shí)發(fā)現(xiàn)女聲低頻能量反而更高原因錄音設(shè)備的頻響特性、錄音距離、音量歸一化都會(huì)影響頻譜。如果女聲錄音音量更大低頻能量可能反而更高。解決先對(duì)兩段音頻做歸一化除以最大絕對(duì)值再計(jì)算頻譜。比較時(shí)看的是頻譜形狀基音間隔、共振峰位置不是絕對(duì)幅度。6. 從實(shí)驗(yàn)到工程采樣定理的邊界與重建的取舍做完這四個(gè)實(shí)驗(yàn)最容易留下的印象是“采樣率大于兩倍最高頻率就萬(wàn)事大吉”。但真實(shí)工程里這個(gè)結(jié)論只在理想條件下成立。信號(hào)不是嚴(yán)格帶限的濾波器不是理想低通采樣時(shí)鐘有抖動(dòng)量化有噪聲。奈奎斯特間隔是一個(gè)下界不是舒適區(qū)。實(shí)際系統(tǒng)里通常把采樣率做到奈奎斯特率的 2.5 到 10 倍就是為了給抗混疊濾波器留過(guò)渡帶給時(shí)鐘抖動(dòng)留裕量。另一個(gè)取舍是重建濾波器的設(shè)計(jì)。理想sinc插值在數(shù)學(xué)上完美但時(shí)域無(wú)限長(zhǎng)無(wú)法實(shí)時(shí)實(shí)現(xiàn)。實(shí)際系統(tǒng)用零階保持器加補(bǔ)償濾波器或者用多相濾波器組做插值。這些方法都會(huì)引入幅度失真和相位失真需要在通帶平坦度和阻帶衰減之間做權(quán)衡。實(shí)驗(yàn)里用sinc插值是在“事后”重建沒(méi)有實(shí)時(shí)性要求所以可以用很長(zhǎng)的窗口。如果換成實(shí)時(shí)音頻流就必須用短濾波器誤差會(huì)更大。還有一個(gè)容易被忽略的點(diǎn)采樣定理保證的是“可以從采樣點(diǎn)唯一恢復(fù)原信號(hào)”但前提是原信號(hào)嚴(yán)格帶限。如果信號(hào)有寬帶噪聲即使采樣率很高噪聲也會(huì)混疊。這時(shí)候需要先做抗混疊濾波把噪聲帶寬限制在奈奎斯特頻率以內(nèi)。實(shí)驗(yàn)里沒(méi)有涉及噪聲但實(shí)際做音頻采集或傳感器信號(hào)處理時(shí)抗混疊濾波器是前端電路里最關(guān)鍵的模擬部分。我自己的習(xí)慣是每次做采樣相關(guān)的設(shè)計(jì)先算奈奎斯特率然后問(wèn)三個(gè)問(wèn)題——信號(hào)真的帶限嗎抗混疊濾波器的過(guò)渡帶夠嗎重建濾波器的相位失真在可接受范圍內(nèi)嗎這三個(gè)問(wèn)題答不上來(lái)采樣率再高也是白搭。從那以后我每次看到“采樣率 44.1 kHz 就夠了”這種說(shuō)法都會(huì)先確認(rèn)信號(hào)帶寬和濾波器特性再下結(jié)論。希望這份實(shí)驗(yàn)報(bào)告里的代碼和現(xiàn)象能幫你把采樣定理從公式變成直覺(jué)。本文還有配套的精品資源點(diǎn)擊獲取