




已閱讀5頁(yè),還剩19頁(yè)未讀, 繼續(xù)免費(fèi)閱讀
版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1 數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告 實(shí)驗(yàn)一 頻譜分析及采樣定理實(shí)驗(yàn)一 頻譜分析及采樣定理 班班 級(jí) 級(jí) 姓姓 名 名 學(xué)學(xué) 號(hào) 號(hào) 2 實(shí)驗(yàn)內(nèi)容實(shí)驗(yàn)內(nèi)容 實(shí)驗(yàn)一實(shí)驗(yàn)一 頻譜分析與采樣定理頻譜分析與采樣定理 一 實(shí)驗(yàn)?zāi)康?1 觀察模擬信號(hào)經(jīng)理想采樣后的頻譜變化關(guān)系 2 驗(yàn)證采樣定理 觀察欠采樣時(shí)產(chǎn)生的頻譜混疊現(xiàn)象 3 加深對(duì) DFT 算法原理和基本性質(zhì)的理解 4 熟悉 FFT 算法原理和 FFT 的應(yīng)用 二 實(shí)驗(yàn)原理 采樣定理 要想抽樣后能不失真的還原出原信號(hào) 抽樣頻率 必須大于等于兩倍信號(hào)頻譜的最高頻率 即 fs 2fh 對(duì)給定信號(hào)進(jìn)行采樣 信號(hào)的頻譜出現(xiàn)周期延拓 三 實(shí)驗(yàn)內(nèi)容和步驟 實(shí)驗(yàn)內(nèi)容 在給定信號(hào)為 1 x t cos 100 at 2 x t exp at 3 x t exp at cos 100 at 其中 a 為實(shí)驗(yàn)者的學(xué)號(hào) 記錄上述各信號(hào)的頻譜 表明采 樣條件 分析比較上述信號(hào)頻譜的區(qū)別 實(shí)驗(yàn)步驟 3 1 復(fù)習(xí)采樣理論 DFT 的定義 性質(zhì)和用 DFT 作譜分析的有 關(guān)內(nèi)容 2 復(fù)習(xí) FFT 算法原理和基本思想 3 確定實(shí)驗(yàn)給定信號(hào)的采樣頻率 編制對(duì)采樣后信號(hào)進(jìn)行 頻譜分析的程序 實(shí)驗(yàn)一 頻譜分析與采樣定理 T 0 0001 采樣間隔T 0 0001 F 1 T 采樣頻率F 1 T L 0 10 幾率長(zhǎng)度L 0 10 N L T 幾率點(diǎn)數(shù) t 0 T L a 5 f1 0 F N F f2 F 2 F N F 2 x1 cos 100 pi a t 輸入信號(hào) y1 T abs fft x1 DFT算法 y11 fftshift y1 FFT算法 subplot 3 1 1 plot t x1 title 正弦信號(hào) subplot 3 1 2 stem f1 y1 title 正弦信號(hào)頻譜 subplot 3 1 3 plot f2 y11 title 正弦信號(hào)頻譜 x2 exp a t y2 T abs fft x2 y21 fftshift y2 figure 2 subplot 3 1 1 stem t x2 title 指數(shù)信號(hào) subplot 3 1 2 stem f1 y2 title 指數(shù)信號(hào)頻譜 subplot 3 1 3 plot f2 y21 title 指數(shù)信號(hào)頻譜 x3 x1 x2 y3 T abs fft x3 y31 fftshift y3 figure 3 subplot 3 1 1 stem t x3 title 兩信號(hào)相乘 subplot 3 1 2 stem f1 y3 title 兩信號(hào)相乘頻譜 subplot 3 1 3 plot f2 y31 title 兩信號(hào)相乘頻譜 4 T 1 10000 fs 10000 L 0 10 00 010 020 030 040 050 060 070 080 090 1 1 0 1 兩 兩 兩 兩 010002000300040005000600070008000900010000 0 0 05 0 1 兩 兩 兩 兩 兩 兩 5000 4000 3000 2000 1000010002000300040005000 0 0 05 0 1 兩 兩 兩 兩 兩 兩 00 010 020 030 040 050 060 070 080 090 1 0 0 5 1 兩 兩 兩 兩 010002000300040005000600070008000900010000 0 0 05 0 1 兩 兩 兩 兩 兩 兩 5000 4000 3000 2000 1000010002000300040005000 0 0 05 0 1 兩 兩 兩 兩 兩 兩 5 00 010 020 030 040 050 060 070 080 090 1 1 0 1 兩 兩 兩 兩 兩 010002000300040005000600070008000900010000 0 0 02 0 04 兩 兩 兩 兩 兩 兩 兩 5000 4000 3000 2000 1000010002000300040005000 0 0 02 0 04 兩 兩 兩 兩 兩 兩 兩 Fs 200 N1 N2 1 時(shí) 線性卷積等于圓周卷積 因此可利用 FFT 計(jì)算線性卷積 三 實(shí)驗(yàn)內(nèi)容和步驟 1 給定離散信號(hào) x n 和 h n 用圖解法求出兩者的線性卷積和 圓周卷積 2 編寫(xiě)程序計(jì)算線性卷積和圓周卷積 3 比較不同列長(zhǎng)時(shí)的圓周卷積與線性卷積的結(jié)果 分析原因 下面是實(shí)驗(yàn)具體可運(yùn)行程序 x 3 0 2 1 3 原始序列 y 3 0 2 1 3 直接計(jì)算圓周卷積或線性卷積 z conv x y figure 1 subplot 311 stem x axis 1 9 0 4 title 原始序列x 10 subplot 312 stem y axis 1 9 0 4 title 原始序列y subplot 313 stem z axis 1 9 0 30 title 直接線性卷積z 利用FFT計(jì)算 N 10 N 10 9時(shí) x1 x zeros 1 N length x y1 y zeros 1 N length y X1 fft x1 Y1 fft y1 Z1 X1 Y1 z1 ifft Z1 figure 2 subplot 321 stem x1 title x1 subplot 322 stem real X1 title X1 subplot 323 stem y1 title y1 subplot 324 stem real Y1 title Y1 subplot 325 stem z1 title z1 subplot 326 stem real Z1 title Z1 N 8 N 8 9時(shí) x2 x zeros 1 N length x y2 y zeros 1 N length y X2 fft x2 Y2 fft y2 Z2 X2 Y2 z2 ifft Z2 figure 3 subplot 321 stem x2 title x2 subplot 322 stem real X2 title X2 subplot 323 stem y2 title y2 subplot 324 stem real Y2 title Y2 subplot 325 stem z2 title z2 subplot 326 stem real Z2 title Z2 11 123456789 0 2 4 兩 兩 兩 兩 x 123456789 0 2 4 兩 兩 兩 兩 y 123456789 0 10 20 30 兩 兩 兩 兩 兩 兩 z N 11 N 11 051015 0 2 4 x1 051015 0 5 10 X1 051015 0 2 4 y1 051015 0 5 10 Y1 051015 0 20 40 z1 051015 100 0 100 Z1 N 10 N 10 12 0510 0 2 4 x1 0510 0 5 10 X1 0510 0 2 4 y1 0510 0 5 10 Y1 0510 0 20 40 z1 0510 100 0 100 Z1 N 8 N 8 02468 0 2 4 x2 02468 10 0 10 X2 02468 0 2 4 y2 02468 10 0 10 Y2 02468 0 20 40 z2 02468 100 0 100 Z2 N 6 N 6 13 0246 0 2 4 x2 0246 10 0 10 X2 0246 0 2 4 y2 0246 10 0 10 Y2 0246 0 20 40 z2 0246 100 0 100 Z2 四 實(shí)驗(yàn)設(shè)備 計(jì)算機(jī) Matlab 軟件 五 實(shí)驗(yàn)結(jié)果分析及心得體會(huì) 1 在編寫(xiě)程序的時(shí)候要注意坐標(biāo)軸幅值的選擇 2 LN1 N2 1 時(shí)才可以用線性卷積代替圓周卷積 14 實(shí)驗(yàn)內(nèi)容 實(shí)驗(yàn)三實(shí)驗(yàn)三 IIR 濾波器設(shè)計(jì)實(shí)驗(yàn)濾波器設(shè)計(jì)實(shí)驗(yàn) 一 實(shí)驗(yàn)?zāi)康?1 學(xué)習(xí)模擬 數(shù)字變換濾波器的設(shè)計(jì)方法 2 掌握雙線性變換濾波器的設(shè)計(jì)方法 3 掌握實(shí)現(xiàn)數(shù)字濾波的具體方法 二 實(shí)驗(yàn)原理 1 IIR數(shù)字低通濾波器設(shè)計(jì)時(shí) 可根據(jù)指標(biāo)要求先設(shè)計(jì)一個(gè)模擬低 通濾波器 選取合適算法逼近這一濾波器 再通過(guò)轉(zhuǎn)化方法的選擇 獎(jiǎng)模擬濾波器轉(zhuǎn)換為數(shù)字濾波器 2 通常模擬濾波器與數(shù)字濾波器轉(zhuǎn)換方法和沖擊響應(yīng)不變法 雙線 性變換法利用的是Z域和S域的相互變換 Z e sT 三 實(shí)驗(yàn)內(nèi)容及步驟 1 用雙線性變換法設(shè)計(jì)一個(gè)巴特沃斯低通IIR數(shù)字濾波器 設(shè)計(jì)指 標(biāo)參數(shù)為 在銅帶內(nèi)頻率低于0 2 時(shí) 最大衰減小于1大B 在租代 內(nèi) 0 3 頻率區(qū)間上 最小衰減大于15dB 先設(shè)計(jì)模擬濾波器 再轉(zhuǎn)化數(shù)字濾波器 設(shè)計(jì)參數(shù)指標(biāo) wp 0 2 pi ws 0 3 pi 數(shù)字濾波器的截止頻率 Rp 1 Rs 15 衰減設(shè)置 Ts 0 02 pi Fs 1 Ts 采樣間隔和采樣頻率 15 wp1 2 Ts tan wp 2 ws1 2 Ts tan ws 2 轉(zhuǎn)換為模擬濾波器截止頻率 N Wn buttord wp1 ws1 Rp Rs s 選擇濾波器的最小階數(shù) N 創(chuàng)建butterworth模擬濾波器 Z P K buttap N 歸一化圓形濾波器設(shè)計(jì) N為階數(shù) Z零點(diǎn) P極點(diǎn) K為增 益 B A zp2tf Z P K 零極點(diǎn)增益模型到傳遞函數(shù)模型的轉(zhuǎn)換 B為傳遞函數(shù)分子 A為分母 b a lp2lp B A Wn 截止頻率為Wn的模擬低通濾波器 用雙線性變換法實(shí)現(xiàn)模擬濾波器到數(shù)字濾波器的轉(zhuǎn)換 bz az bilinear b a Fs 數(shù)字低通濾波器的系數(shù) 回執(zhí)頻率響應(yīng)曲線 H W freqz bz az 50 H為幅度 W為相位 L length W 2 1 figure 1 plot W 1 L pi abs H 1 L grid xlabel 角頻率 pi ylabel 頻率響應(yīng)幅度 x 4 8 12 12 10 6 6 6 4 0 0 0 0 0 2 4 0 0 0 2 2 0 0 2 2 2 2 0 4 2 0 4 6 4 2 4 6 6 4 4 6 6 2 6 12 8 0 16 38 60 84 90 66 32 4 2 y filter bz az x 濾波 figure 2 subplot 2 1 1 plot x title 原始信號(hào) subplot 2 1 2 plot y title 濾波后信號(hào) 濾波器的頻率響應(yīng) 16 012345678 0 0 1 0 2 0 3 0 4 0 5 0 6 0 7 0 8 0 9 1 兩 兩 兩 Hz 兩 兩 兩 兩 兩 兩 濾波效果圖 0102030405060 100 50 0 50 兩 兩 兩 兩 0102030405060 100 50 0 50 兩 兩 兩 兩 兩 17 實(shí)驗(yàn)內(nèi)容 實(shí)驗(yàn)四 FIR 濾波器設(shè)計(jì)實(shí)驗(yàn) 一 實(shí)驗(yàn)?zāi)康?1 熟悉濾波器的計(jì)算機(jī)仿真方法 2 掌握用窗函數(shù)法設(shè)計(jì)FIR 數(shù)字濾波器的原理和方法 3 解各種窗函數(shù)對(duì)濾波特性的影響 18 二 實(shí)驗(yàn)原理 1 可通過(guò)加窗把無(wú)限長(zhǎng)序列變成有限長(zhǎng)序列 從而使用 FFT 快速 計(jì)算卷積 提高運(yùn)算效率 但加窗會(huì)對(duì)頻譜造成影響 因此要根據(jù) 需求注意窗函數(shù)的選擇 2 吉布斯現(xiàn)象 改變N只改變窗譜的主瓣寬度 不能改變主瓣和旁 瓣的相對(duì)比例 因而肩峰的相對(duì)值不隨意改變 改變的只是過(guò)渡帶 寬 三 實(shí)驗(yàn)要求 1 設(shè)計(jì)一線性相位FIR 低通濾波器濾波器 給定抽樣頻率為 s 3 10 rad s 通帶截止頻率為 p 3 10 rad s 阻帶起始頻率為 st 6 10 rad s 帶衰減比小于50dB 2 選擇不同的窗函數(shù)設(shè)計(jì)該濾波器 觀察其頻率響應(yīng)函數(shù)有什么變 化 下面是實(shí)驗(yàn)具體可運(yùn)行程序 海明窗 wp 3000 pi 通帶截止角頻率 ws 6000 pi 阻帶截止角頻率 wsam 30000 pi 采樣角頻率 fsam wsam 2 pi 采樣截止頻率 passrad wp ws 2 fsam 截止頻率 A 3 3 wdelta ws wp fsam N ceil 2 pi wdelta A 濾波器的階數(shù) w hamming N 1 用海明窗實(shí)現(xiàn) L N 2 1 19 n 1 1 N 1 hd sin passrad n L pi n L 理想低通濾波器 if N ceil N 2 2 hd L passrad pi end h hd w 加窗 mag rad freqz h omega linspace 0 pi 512 magdb 20 log10 abs mag subplot 111 plot omega pi magdb k 繪制對(duì)數(shù)幅度特性曲線 axis 0 1 100 5 xlabel omega pi ylabel 20lg H e j omega grid on 海明窗 00 511 522 533 5 120 100 80 60 40 20 0 20 布萊克曼窗 wp 3000 pi 通帶截止角頻率 20 ws 6000 pi 阻帶截止角頻率 wsam 30000 pi 采樣角頻率 fsam wsam1 2 pi 采樣截止頻率 passrad wp ws 2 fsam 截止頻率 A 5 5 wdelta ws wp fsam N ceil 2 pi wdelta A 濾波器的階數(shù) w blackman N 1 用布萊克曼窗實(shí)現(xiàn) L N 2 1 n 1 1 N 1 hd sin passrad n L pi n L 理想低通濾波器 if N ceil N 2 2 hd L passrad pi end h hd w 加窗 mag rad freqz h magdb 20 log10 abs mag omega linspace 0 pi 512 subplot 111 plot omega pi magdb k 繪制對(duì)數(shù)幅度特性曲線 axis 0 1 100 5 xlabel omega pi ylabel 20lg H e j omega grid on 布萊克曼窗 21 00 511 522 533 5 140 120 100 80 60 40 20 0 20 凱澤窗 wp1 2 pi 3000 pi 30000 pi 通帶截止角頻率 ws1 2 pi 6000 pi 30000 pi 阻帶截止角頻率 rs 50 wsam1 30000 pi 采樣角頻率 fsam1 wsam1 2 pi 采樣截止頻率 wc ws1 wp1 2 pi 截止頻率 beta 0 1102 rs 8 7 N ceil rs 7 95 2 286 ws1 wp1 1 濾波器的階數(shù) hdn kaiser N beta 用凱澤窗實(shí)現(xiàn) hn fir1 N 1 wc low kaiser N beta 用凱澤窗函數(shù)設(shè)計(jì)低通濾波器 figure 1 subplot 111 stem 0 N 1 hn k 繪制濾波器時(shí)域波形 axis 0 N 1 1 0 5 xlabel n ylabel h n omega linspace 0 pi 512 mag freqz hn 1 omega magdb 20 log10 abs mag figure 2 subplot 111 plot omega pi m
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 三年級(jí)語(yǔ)文下冊(cè)第六單元第23課和時(shí)間賽跑教學(xué)實(shí)錄語(yǔ)文S版
- 樹(shù)林管理合同
- 感應(yīng)移相器及同步器企業(yè)ESG實(shí)踐與創(chuàng)新戰(zhàn)略研究報(bào)告
- 管線管企業(yè)數(shù)字化轉(zhuǎn)型與智慧升級(jí)戰(zhàn)略研究報(bào)告
- 塑料熱成型機(jī)企業(yè)數(shù)字化轉(zhuǎn)型與智慧升級(jí)戰(zhàn)略研究報(bào)告
- 自動(dòng)半自動(dòng)超聲波焊接機(jī)企業(yè)數(shù)字化轉(zhuǎn)型與智慧升級(jí)戰(zhàn)略研究報(bào)告
- 原生鎵(礦產(chǎn)鎵)企業(yè)ESG實(shí)踐與創(chuàng)新戰(zhàn)略研究報(bào)告
- 紅外熱成像儀企業(yè)縣域市場(chǎng)拓展與下沉戰(zhàn)略研究報(bào)告
- 冷藏或保溫掛車企業(yè)數(shù)字化轉(zhuǎn)型與智慧升級(jí)戰(zhàn)略研究報(bào)告
- 郵政代辦服務(wù)企業(yè)數(shù)字化轉(zhuǎn)型與智慧升級(jí)戰(zhàn)略研究報(bào)告
- 氨吹脫塔單元設(shè)計(jì)示例
- 中國(guó)移動(dòng)-安全-L3
- 骨齡評(píng)測(cè)方法課件
- GB/T 42314-2023電化學(xué)儲(chǔ)能電站危險(xiǎn)源辨識(shí)技術(shù)導(dǎo)則
- 人教小學(xué)數(shù)學(xué)五年級(jí)下冊(cè)綜合與實(shí)踐《怎樣通知最快》示范公開(kāi)課教學(xué)課件
- 海陸熱力性質(zhì)差異的說(shuō)課課件
- 科學(xué)院大學(xué)博士入學(xué)復(fù)試匯報(bào)個(gè)人介紹PPT模板
- 四川省稅務(wù)局財(cái)行處土增稅清算復(fù)審指導(dǎo)口徑(2021年6月22日)
- 迎春杯2023年-2023年中高年級(jí)初賽復(fù)賽試題真題整理
- GB/T 6322-1986光滑極限量規(guī)型式和尺寸
- GB/T 31052.5-2015起重機(jī)械檢查與維護(hù)規(guī)程第5部分:橋式和門(mén)式起重機(jī)
評(píng)論
0/150
提交評(píng)論