【正文】 hd = sin(wc*m) ./ (pi*m)。 n = [0:1:(M1)]。 西安歐亞學 院課程設計 15 附錄 2. 理想低通濾波器計算源程序 function hd = ideal_lp(wc,M)。 grd(2:1:501) 0]。 % grd = [0 grd(1:1:500)。 % grd = diff(pha)。 % pha = unwrap(angle(H))。 db = 20*log10((mag+eps)/max(mag))。 w = (w(1:1:501))39。)。 % db = [0 到 pi 弧度 ]區(qū)間內(nèi)的相對振幅 (db) % mag = [0 到 pi 弧度 ]區(qū)間內(nèi)的絕對振幅 % pha = [0 到 pi 弧度 ]區(qū)間內(nèi)的相位響應 % grd = [0 到 pi 弧度 ]區(qū)間內(nèi)的群遲延 % w = [0 到 pi 弧度 ]區(qū)間內(nèi) 501 個頻率樣本向量 % b = Ha(z)的分子多項式系數(shù) (對 FIR b=h) % a = Ha(z)的分母多項式系數(shù) (對 FIR: a=[1]) % [H,w] = freqz(b,a,1000,39。 西安歐亞學 院課程設計 13 參考文獻 [1] 楊會成， 數(shù)字信號處理 ，國防工業(yè)出版社， 2020 [2] 程佩青．數(shù)字信號處理教程 [M]， 北京：清華大學出版社， 2020 西安歐亞學 院課程設計 14 附錄 附錄 1. 濾波器脈沖響應源程序 function [db,mag,pha,grd,w] = freqz_m(b,a)。 最后，感謝學校能夠開設這樣的課程，如果學校沒有開設這樣的課，我很可能也不會學到這么多理論課以外的知識。三 個臭皮匠能頂一個諸葛亮，團隊合作是一種很好的方式，它可以彌補個人之間的不足， 大眾的智慧才是真正的智慧。如果沒有他們無私的幫助與指導，我想我就不會那么快的得到實驗結果。 在此次的課程設計中除了感謝老師給我們的幫助，給我最大幫助的就是我的同學們了，在此期間，同學們互相的解答彼此的疑惑。 西安歐亞學 院課程設計 12 致謝 在整個設計過程中， 魏老師給了我們很大的幫助，從選題到寫報告她都悉心的給我們講解，在我們遇到困難的時候，總是孜孜不倦的為我們解惑。 課程設計不僅是對前面所學知識的一種檢驗，而且也是對自己能力的一種提高。 平時的理論知識只有通過自己動手做過課題，做過實驗，才能從中發(fā)現(xiàn)問題，解決問題。 經(jīng)過這次的課程設計，我了解到了做任何事都要有耐心、更是要細心做事。有的事不是僅僅靠自己一個人苦思冥想就可以的到結果的，一個人的能力與智慧始終是有限的，要想高效率的完成一項任務，只有通過團隊合作才能盡快得出結果。它加深了老師與同學，同學與同學之間的情感。 在此期間，老師與同學們打成一片，同學之間相互請教、探討、指導。 5 總結與展望 本次 課程 設計 歷時一個星期，雖然只是短短的一周，但卻讓我們學到了許多。 As 和 Rp 是否達標，根本不了解在什么樣的情況下 As 和 Rp 才能符合濾波器的要求； 解決方法：通過請教老師和尋找同學的幫助解決這一問題，慢慢的得到改善，最終符合設計要求。 4 調(diào)試問題與解決方案 FIR 濾波器時得不到理想的濾波器； 0 1 2 3 0 . 200 . 2時間 ( 單位 : s )幅度原始語音信號 x0 1000 2020 3000 4000050010001500頻率 ( 單位 : H z )幅度譜語音信號幅度譜 X0 1 2 3 0 . 200 . 2時間 ( 單位 : s )幅度加干擾后的語音信號 x10 1000 2020 3000 4000050010001500頻率 ( 單位 : H z )幅度譜加干擾語音信號幅度譜 X10 1 2 3 0 . 200 . 2時間 ( 單位 : s )幅度濾波后語音信號 y0 1000 2020 3000 4000050010001500頻率 ( 單位 : H z )幅度譜濾波后語音信號幅度譜 Y圖 34語音信號去噪前后時域圖及頻域幅度圖 西安歐亞學 院課程設計 11 解決方法： 這是因為通帶截止頻率與阻帶起始頻率之間的差值設置的太小或者太大，通過適當?shù)倪x擇參數(shù) fp fs fs2 及 fp2，繪制出來的圖形效果比較明顯，基本符合設計指標。其調(diào)用格式為： sound (x， fs， bits)。在將三個信號的時域波形和頻譜圖比較之后，我們還要通過回放去濾波去噪音樂信號，來跟原信號相比，以檢驗濾波器的效果。 仿真結果與分析 通過觀察濾波前后語音信號波形的變化，即觀察圖 33，我們可以知道，原信號的時域圖與濾波去噪信號的時域圖基本相似；原信號與濾波去噪信號的頻譜圖波形也大致相似。grid on。濾波后語音信號幅度譜 Y39。)。ylabel(39。頻率 (Hz)39。plot(f,Y_fil)。grid on。濾波后語音信號 y39。)。ylabel(39。時間 (s)39。plot(t,y_fil)。grid on。加干擾語音信號幅度譜 X139。)。ylabel(39。頻率 (Hz)39。axis([0 4000 0 70])。 subplot(3,2,4)。)。title(39。幅度 39。)。 xlabel(39。plot(t,y)。grid on。語音信號幅度譜 X39。)。ylabel(39。頻率 (Hz)39。axis([0 4000 0 70])。 subplot(3,2,2)。)。title(39。幅度 39。)。 xlabel(39。plot(t,x)。Y_fil=Y_fil(1:N/2)。具體程序如下： y_fil=fftfilt (h_bs, y)。 所設計的濾波器的幅度響應和脈沖響應如圖 33所示： 0 0 . 0 5 0 . 1 0 . 1 5 0 . 2 6 0 4 0 2 00w / p idB濾波器幅度響應圖0 0 . 0 5 0 . 1 0 . 1 5 0 . 200 . 51w / p i幅度mag濾波器幅度響應圖0 0 . 5 1202w / p i相位pha濾波器相位響應圖0 500 1000 150000 . 51nh(n)濾波器脈沖響應圖圖 33 濾波器的頻率特性圖 西安歐亞學 院課程設計 9 由 33 圖可以看出 Rp 接近于 0，而 As大于開始的設定值（ As=60db)，由此可見滿足了設計要求，達到了設計目的。)。title(39。h(n)39。)。 xlabel(39。stem(n,h_bs)。grid on。濾波器相