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