版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
1、<p> 課 程 設(shè) 計 報 告</p><p> 課程名稱 數(shù)字信號處理 </p><p> 課題名稱 數(shù)字濾波器設(shè)計及在語音信號分析中的應(yīng)用 </p><p> 專 業(yè) 通信工程 </p><p>
2、班 級 通信工程1101 </p><p> 學(xué) 號 </p><p> 姓 名 </p><p> 指導(dǎo)教師 &l
3、t;/p><p> 2013年 9月 8日</p><p> 課 程 設(shè) 計 任 務(wù) 書</p><p> 課程名稱 數(shù)字信號處理 </p><p> 課 題 數(shù)字信號濾波器設(shè)計及在語音信號中的</p><p> 分析中的應(yīng)用 </p>
4、<p> 專業(yè)班級 通信工程1101 </p><p> 學(xué)生姓名 </p><p> 學(xué) 號 </p><p> 指導(dǎo)老師 </p><p> 審
5、 批 </p><p> 任務(wù)書下達日期 2013 年 9月 1 日</p><p> 任務(wù)完成日期 2013 年 9月 8日</p><p> 《數(shù)字信號處理》課程設(shè)計任務(wù)書</p><p><b> 一 、設(shè)計目的</b></p&g
6、t;<p> 綜合運用數(shù)字信號處理的理論知識進行頻譜分析和濾波器設(shè)計,通過理論推導(dǎo)得出相應(yīng)結(jié)論,再利用 MATLAB 作為編程工具進行計算機實現(xiàn),從而加深對所學(xué)知識的理解,建立概念。</p><p><b> 二、設(shè)計要求</b></p><p> 1、 MATLAB 的使用,掌握 MATLAB 的程序設(shè)計方法。</p><p&
7、gt; 2、 Windows 環(huán)境下語音信號采集的方法。</p><p> 3、數(shù)字信號處理的基本概念、基本理論和基本方法。</p><p> 4、TLAB 設(shè)計 FIR 和 nR 數(shù)字濾波器的方法。</p><p> 5、用 MATLAB 對信號進行分析和處理</p><p> 6、計報告4000以上,含程序設(shè)計說明,用戶使用說明
8、,源程序清單及程序框圖。</p><p><b> 7、機演示。</b></p><p> 8、有詳細的文檔。文檔中包括設(shè)計思路、設(shè)計仿真程序、仿真結(jié)果及相應(yīng)的分析與結(jié)論。</p><p><b> 三、進度安排</b></p><p> 第一周 星期一: 課題講解,查閱資料<
9、/p><p> 星期二: 總體設(shè)計,詳細設(shè)計</p><p> 星期三: 編程,上機調(diào)試、修改程序</p><p> 星期四: 上機調(diào)試、完善程序</p><p><b> 星期五: 答辯</b></p><p> 星期六-星期天:撰寫課程設(shè)計報告</p><
10、p><b> 附:</b></p><p> 課程設(shè)計報告裝訂順序:封面、任務(wù)書、目錄、正文、評分、附件(A4大小的圖紙及程序清單)。 </p><p> 正文的格式:一級標題用3號黑體,二級標題用四號宋體加粗,正文用小四號宋體;行距為22。</p><p> 正文的內(nèi)容:一、課題的主要功能;二、課題的功能模塊的劃分;三、主要功
11、能的實現(xiàn);四、程序調(diào)試;五、總結(jié);六、附件(所有程序的原代碼,要求對程序?qū)懗霰匾淖⑨專?;七、評分表。</p><p><b> 目 錄</b></p><p> 1、 課題的主要功能1</p><p> 1.1 問題描述1</p><p> 1.2 功能要求1</p><p>
12、 2、課題的功能模塊的劃分2</p><p> 3、 濾波器實現(xiàn)的原理3</p><p> 3.1 FIR濾波器3</p><p> 3.2 窗函數(shù)設(shè)計法3</p><p><b> 4、 程序調(diào)試4</b></p><p> 4.1 對原始信號進行時域和頻域的分析4&l
13、t;/p><p> 4.2 含噪語音信號的合成5</p><p> 4.3 數(shù)字濾波器的設(shè)計及濾波8</p><p><b> 5、 總結(jié)16</b></p><p><b> 6、 附件16</b></p><p><b> 7、評分表21<
14、;/b></p><p><b> 課題的主要功能</b></p><p><b> 1.1 問題描述</b></p><p> 語音信號采集(每個同學(xué)必須錄制本人的語音信號)</p><p> 錄制一段課程設(shè)計學(xué)生的語音信號并保存為文件,要求長度不小于10秒,并對錄制的信號進行采樣;
15、錄制時可以使用Windows自帶的錄音機,或者使用其它專業(yè)的錄音軟件,錄制時需要配備錄音硬件(如麥克風(fēng)),為便于比較,需要在安靜、干擾小的環(huán)境下錄音。</p><p><b> 語音信號分析</b></p><p> 使用MATLAB繪出采樣后的語音信號的時域波形和頻譜圖。根據(jù)頻譜圖求出其帶寬,并說明語音信號的采樣頻率不能低于多少赫茲。</p>&l
16、t;p><b> 含噪語音信號合成</b></p><p> 在MATLAB軟件平臺下,給原始的語音信號疊加上噪聲,噪聲類型分為如下幾種:(1)白噪聲;(2)單頻噪色(正弦干擾);(3)多頻噪聲(多正弦干擾);(4)其它干擾,可設(shè)置為低頻、高頻、帶限噪聲,或Chirp干擾、沖激干擾。繪出疊加噪聲后的語音信號時域和頻譜圖,在視覺上與原始語音信號圖形對比,也可通過Windows播放軟件
17、從聽覺上進行對比,分析并體會含噪語音信號頻譜和時域波形的改變。</p><p><b> 1.2 功能要求</b></p><p> 給定濾波器的規(guī)一化性能指標(參考指標,實際中依據(jù)每個同學(xué)所疊加噪聲情況而定)例如:通帶截止頻率wp=0.25*pi, 阻通帶截止頻率ws=0.3*pi; 通帶最大衰減Rp=1 dB; 阻帶最小衰減Rs=15 dB。</p&g
18、t;<p> 采用窗函數(shù)法(至少采用兩種以上不同的窗函數(shù))分別設(shè)計各型FIR濾波器(低通、高通、帶通、帶阻中的至少3種類型)來對疊加噪聲前后的語音信號進行濾波處理,繪出濾波器的頻域響應(yīng),繪出濾波后信號的時域波形和頻譜,并對濾波前后的信號進行對比,分析信號的變化;在相同的性能指標下比較各方法的濾波效果,并從理論上進行分析(或解釋)。</p><p> 2、課題的功能模塊的劃分</p>
19、<p> 根據(jù)設(shè)計的要求,首先自己錄制一段語音信號,對語音信號加入噪聲干擾,再利用矩形窗和哈明窗設(shè)計合理的FIR濾波器。最后用濾波器對干擾后的語音信號進行濾波去噪。具體設(shè)計流程圖如圖所示:</p><p><b> 否</b></p><p><b> 是</b></p><p><b> 圖
20、 1 設(shè)計流程圖</b></p><p> 本次濾波器設(shè)計分為原始語音信號模塊,噪音模塊,濾波器設(shè)計模塊,濾波器濾波模塊。其中噪音模塊加的有高斯白噪聲,單頻余弦信號,多頻余弦信號噪音。濾波器設(shè)計模塊,我選用的是矩形窗的低通和帶通濾波器設(shè)計,哈明窗的低通和帶通濾波器設(shè)計。濾波器濾波模塊經(jīng)過多次調(diào)整,濾波效果已趨完善。</p><p><b> 濾波器實現(xiàn)的原理<
21、;/b></p><p> 3.1 FIR濾波器</p><p> FIR濾波器:有限長單位沖激響應(yīng)濾波器,是數(shù)字信號系統(tǒng)中最基本的元件,它可以在保證任意幅頻特性的同時具有嚴格的線性相頻特性,同時其單位抽樣響應(yīng)是有限長的,因而濾波器是穩(wěn)定的系統(tǒng)。因此,F(xiàn)IR濾波器在通信、圖像處理、模式識別等領(lǐng)域都有著廣泛的應(yīng)用。在設(shè)計一個濾波器之前,必須首先根據(jù)工程實際的需要確定濾波器的技術(shù)指標
22、。在很多實際應(yīng)用中,數(shù)字濾波器常常被用來實現(xiàn)選頻操作。因此,指標的形式一般在頻域中給出幅度和相位響應(yīng)。幅度指標主要以兩種方式給出。第一種是絕對指標。它提供對幅度響應(yīng)函數(shù)的要求,一般應(yīng)用于FIR濾波器的設(shè)計。第二種指標是相對指標。它以分貝值的形式給出要求。</p><p> 確定了技術(shù)指標后,就可以建立一個目標的數(shù)字濾波器模型。通常用理想的數(shù)字濾波器模型。之后,利用數(shù)字濾波器的設(shè)計方法,設(shè)計出一個實際濾波器模型來
23、逼近給定的目標。根據(jù)這個描述就可以分析其頻率特性和相位特性,以驗證設(shè)計結(jié)果是否滿足指標要求;或者利用計算機仿真實現(xiàn)設(shè)計的濾波器,再分析濾波結(jié)果來判斷。FIR濾波器的設(shè)計問題實際上是確定能滿足所要求的轉(zhuǎn)移序列或脈沖響應(yīng)的常數(shù)的問題,設(shè)計方法主要有窗函數(shù)、頻率響應(yīng)法和等波紋最佳逼近法等。</p><p> 3.2 窗函數(shù)設(shè)計法</p><p> 窗函數(shù)設(shè)計法是一種通過截斷和計權(quán)的方法使無限
24、長非因果序列成為有限長脈沖響應(yīng)序列的設(shè)計方法。通常在設(shè)計濾波器之前,應(yīng)該先根據(jù)具體的工程應(yīng)用確定濾波器的技術(shù)指標。在大多數(shù)實際應(yīng)用中,數(shù)字濾波器常常被用來實現(xiàn)選頻操作,所以指標的形式一般為在頻域中以分貝值給出的相對幅度響應(yīng)和相位響應(yīng)。</p><p> 窗函數(shù)設(shè)計法步驟如下:</p><p> ?。?)根據(jù)過渡帶寬及阻帶衰減要求,選擇窗函數(shù)的類型并估計窗口長度N。窗函數(shù)的類型可根據(jù)最小阻
25、帶衰減AS獨立選擇。</p><p> (2)根據(jù)待求濾波器的理想頻率響應(yīng)求出理想單位脈沖響應(yīng)hd(n)。</p><p> ?。?)由性能指標確定窗函數(shù)W(n)和長度N。</p><p> ?。?)求得實際濾波器的單位脈沖響應(yīng)h(n), h(n)即為所設(shè)計FIR濾波器系數(shù)向量b(n)。</p><p><b> 程序調(diào)試<
26、;/b></p><p> 4.1 對原始信號進行時域和頻域的分析</p><p> [x1,fs,bits]=wavread('d:/text');</p><p> %sound(x1,fs,bits);</p><p> y1=fft(x1,32768);</p><p><
27、b> figure(1)</b></p><p> subplot(2,1,1)</p><p><b> plot(x1);</b></p><p> title('原始語音信號時域波形');</p><p> subplot(2,1,2)</p><p&
28、gt; plot(abs(y1));</p><p> title('原始語音信號頻譜')</p><p> 圖 2 原始信號波形</p><p> 4.2 含噪語音信號的合成</p><p> x2=awgn(x1,10);</p><p> %sound(x2,fs,bits)<
29、/p><p> y2=fft(x2,32768);</p><p><b> figure(2)</b></p><p> subplot(2,1,1)</p><p><b> plot(x2);</b></p><p> title('加高斯白噪聲語音信號的
30、時域波形')</p><p> subplot(2,1,2)</p><p> plot(abs(y2));</p><p> title('加高斯白噪聲語音信號的頻譜')</p><p> N=length(x1)-1;</p><p> t=0:1/fs:N/fs;</p&g
31、t;<p> d=[0.8*cos(2*pi*5000*t)]';</p><p> d2=[0.8*cos(2*pi*20*t)]';</p><p><b> x3=x1+d;</b></p><p> x4=x1+d2+d;</p><p> %sound(x3,fs,bit
32、s);</p><p><b> figure(3)</b></p><p> subplot(2,1,1)</p><p><b> plot(x3)</b></p><p> title('加單頻余弦信號的語音信號時域波形')</p><p>
33、y3=fft(x3,32768);</p><p> subplot(2,1,2)</p><p> plot(abs(y3));</p><p> title('加單頻余弦信號的語音信號頻譜')</p><p> figure(12)</p><p> subplot(2,1,1)</
34、p><p><b> plot(x4)</b></p><p> title('加多頻余弦信號的語音信號時域波形')</p><p> y4=fft(x4,32768);</p><p> subplot(2,1,2);</p><p> plot(abs(y4));<
35、/p><p> title('加多頻余弦信號的語音信號頻譜')</p><p> 圖3 高斯白噪聲波形</p><p> 圖 4 加單頻噪聲波形</p><p> 圖 5 加多頻余弦信號波形</p><p> 4.3 數(shù)字濾波器的設(shè)計及濾波</p><p> 矩形窗低
36、通濾波器設(shè)計及濾波結(jié)果</p><p><b> %矩形窗 %低通</b></p><p> wp=0.3*pi;ws=0.4*pi;</p><p><b> Bt=ws-wp;</b></p><p> N=ceil(1.8*pi/Bt);</p><p> w
37、c=(wp+ws)/2/pi;</p><p> hn=fir1(N,wc,boxcar(N+1));</p><p> [H,f]=freqz(hn,1,512,8000); %采用8000Hz的采樣頻率求出頻率響應(yīng)</p><p><b> figure(4)</b></p><p> subplot(2,1
38、,1),plot(f,20*log10(abs(H)))</p><p> xlabel('頻率/Hz');ylabel('振幅/dB');grid on;</p><p> subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))</p><p> xlabel('頻率/Hz'
39、;);ylabel('相位/^o');grid on;</p><p> y4=fftfilt(hn,x3); %采用fftfilt對輸入信號濾波</p><p> %sound(y4,fs,bits);</p><p><b> figure(5)</b></p><p>
40、 subplot(2,1,1), plot(t,x3),title('輸入信號')%繪出輸入信號波形</p><p> subplot(2,1,2),plot(t,y4) %繪出輸出信號波形</p><p> 圖 6 矩形窗低通濾波器</p><p> 圖 7 矩形窗低通濾波結(jié)果</p>&l
41、t;p> 矩形帶通濾波器設(shè)計及濾波結(jié)果</p><p> wls=0.25*pi;wlp=0.45*pi;whs=0.625*pi;whp=0.7*pi;</p><p> delta_w=whp-whs;</p><p> wp=[0.35 0.66];N=48; %通帶邊界頻率(歸一化頻率)和濾波器階數(shù)</p><p>
42、<b> Fs=8000;</b></p><p> b=fir1(N,wp); %設(shè)計FIR帶通濾波器</p><p><b> figure(6)</b></p><p> [H,f]=freqz(b,1,512,Fs); %以50Hz為采樣頻率求出濾波器頻率響應(yīng)</p&
43、gt;<p> subplot(2,1,1),plot(f,20*log10(abs(H)))</p><p> xlabel('頻率/Hz');ylabel('振幅/dB');grid on;</p><p> subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))</p><p
44、> xlabel('頻率/Hz');ylabel('相位/^o');grid on;</p><p> y5=fftfilt(b,x4); %采用fftfilt對輸入信號濾波</p><p> %sound(y5,fs,bits);</p><p><b> figure(7)</b
45、></p><p> subplot(2,1,1), plot(t,x4),title('輸入信號')%繪出輸入信號波形</p><p> subplot(2,1,2),plot(t,y5) %繪出輸出信號波形</p><p> hold on;plot(N/2*0.02*ones(1,2),ylim,
46、'r') %繪制延遲到的時刻</p><p> title('輸出信號'),xlabel('時間/s')</p><p> 圖 8 矩形窗帶通濾波器設(shè)計</p><p> 圖 9 矩形窗帶通濾波器濾波結(jié)果</p><p> 哈明窗低通濾波器設(shè)計及濾波結(jié)果</p>&
47、lt;p><b> %哈明窗 %低通</b></p><p> wp=0.5*pi;ws=0.6*pi;</p><p><b> Bt=ws-wp;</b></p><p> N=ceil(6.6*pi/Bt);</p><p> wc=(wp+ws)/2/pi;</p>
48、;<p> hn=fir1(N,wc,hamming(N+1));</p><p> [H,f]=freqz(hn,1,512,8000); %采用8000Hz的采樣頻率求出頻率響應(yīng)</p><p><b> figure(8)</b></p><p> subplot(2,1,1),plot(f,20*log10(abs
49、(H)))</p><p> xlabel('頻率/Hz');ylabel('振幅/dB');grid on;</p><p> subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))</p><p> xlabel('頻率/Hz');ylabel('相位/^o
50、9;);grid on;</p><p> y6=fftfilt(hn,x3); %采用fftfilt對輸入信號濾波</p><p> %sound(y6,fs,bits);</p><p><b> figure(9)</b></p><p> subplot(2,1,1), plot(t
51、,x3),title('輸入信號')%繪出輸入信號波形</p><p> subplot(2,1,2),plot(t,y6) %繪出輸出信號波形</p><p> title('輸出信號'),xlabel('時間/s')</p><p> 圖 10 哈明窗低通濾波器</p
52、><p> 圖 11 哈明窗低通濾波器濾波結(jié)果</p><p> 哈明窗帶通濾波器設(shè)計及濾波結(jié)果</p><p><b> %哈明窗 %帶通</b></p><p><b> %Samp7_6</b></p><p> wp=[0.35 0.65];N=48;
53、%通帶邊界頻率(歸一化頻率)和濾波器階數(shù)</p><p><b> Fs=8000;</b></p><p> b=fir1(N,wp); %設(shè)計FIR帶通濾波器</p><p> figure(10)</p><p> [H,f]=freqz(b,1,512,Fs); %以50
54、Hz為采樣頻率求出濾波器頻率響應(yīng)</p><p> subplot(2,1,1),plot(f,20*log10(abs(H)))</p><p> xlabel('頻率/Hz');ylabel('振幅/dB');grid on;</p><p> subplot(2,1,2),plot(f,180/pi*unwrap(angl
55、e(H)))</p><p> xlabel('頻率/Hz');ylabel('相位/^o');grid on;</p><p> y7=fftfilt(b,x4); %采用fftfilt對輸入信號濾波</p><p> sound(y7,fs,bits);</p><p> fi
56、gure(11)</p><p> subplot(2,1,1), plot(t,x4),title('輸入信號')%繪出輸入信號波形</p><p> subplot(2,1,2),plot(t,y7) %繪出輸出信號波形</p><p> hold on;plot(N/2*0.02*ones(1,2),yl
57、im, 'r') %繪制延遲到的時刻</p><p> title('輸出信號'),xlabel('時間/s')</p><p> 圖 12 哈明窗帶通濾波器</p><p> 圖13 哈明窗帶通濾波器濾波結(jié)果</p><p><b> 總結(jié)</b><
58、;/p><p> 這次課設(shè)是我自己辛辛苦苦做出的,雖然限于時間短暫濾波器還有一些問題,低通濾波器濾波結(jié)果很不理想,但是那種努力換得成功的感覺還是讓我著實讓我開心了許久??偟母杏X:首先從以前對濾波器設(shè)計,語音信號降噪幾乎完全不懂慢慢有了一定了解,感覺它不再那么神秘和難以理解,自信心增長了不少;其次,對于matlab的使用更加得心應(yīng)手,增強了專業(yè)與工具的聯(lián)系。最后,真心感覺很多課設(shè)是我自己可以做出來的,只是你在成功之前
59、就已經(jīng)放棄,也就是說你過不了自己這關(guān),而不是課設(shè)。</p><p><b> 附件</b></p><p><b> 源程序代碼:</b></p><p> [x1,fs,bits]=wavread('d:/text');</p><p> %sound(x1,fs,bits
60、);</p><p> y1=fft(x1,32768);</p><p><b> figure(1)</b></p><p> subplot(2,1,1)</p><p><b> plot(x1);</b></p><p> title('原始語音信
61、號時域波形');</p><p> subplot(2,1,2)</p><p> plot(abs(y1));</p><p> title('原始語音信號頻譜')</p><p> x2=awgn(x1,10);</p><p> %sound(x2,fs,bits)</p&
62、gt;<p> y2=fft(x2,32768);</p><p><b> figure(2)</b></p><p> subplot(2,1,1)</p><p><b> plot(x2);</b></p><p> title('加高斯白噪聲語音信號的時域波
63、形')</p><p> subplot(2,1,2)</p><p> plot(abs(y2));</p><p> title('加高斯白噪聲語音信號的頻譜')</p><p> N=length(x1)-1;</p><p> t=0:1/fs:N/fs;</p>
64、<p> d=[0.8*cos(2*pi*5000*t)]';</p><p> d2=[0.8*cos(2*pi*20*t)]';</p><p><b> x3=x1+d;</b></p><p> x4=x1+d2+d;</p><p> %sound(x3,fs,bits);
65、</p><p><b> figure(3)</b></p><p> subplot(2,1,1)</p><p><b> plot(x3)</b></p><p> title('加單頻余弦信號的語音信號時域波形')</p><p> y3=
66、fft(x3,32768);</p><p> subplot(2,1,2)</p><p> plot(abs(y3));</p><p> title('加單頻余弦信號的語音信號頻譜')</p><p> figure(12)</p><p> subplot(2,1,1)</p&g
67、t;<p><b> plot(x4)</b></p><p> title('加多頻余弦信號的語音信號時域波形')</p><p> y4=fft(x4,32768);</p><p> subplot(2,1,2);</p><p> plot(abs(y4));</p&
68、gt;<p> title('加多頻余弦信號的語音信號頻譜')</p><p><b> %矩形窗 %低通</b></p><p> wp=0.3*pi;ws=0.4*pi;</p><p><b> Bt=ws-wp;</b></p><p> N=ceil
69、(1.8*pi/Bt);</p><p> wc=(wp+ws)/2/pi;</p><p> hn=fir1(N,wc,boxcar(N+1));</p><p> [H,f]=freqz(hn,1,512,8000); %采用8000Hz的采樣頻率求出頻率響應(yīng)</p><p><b> figure(4)</b&g
70、t;</p><p> subplot(2,1,1),plot(f,20*log10(abs(H)))</p><p> xlabel('頻率/Hz');ylabel('振幅/dB');grid on;</p><p> subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))</p&g
71、t;<p> xlabel('頻率/Hz');ylabel('相位/^o');grid on;</p><p> y4=fftfilt(hn,x3); %采用fftfilt對輸入信號濾波</p><p> %sound(y4,fs,bits);</p><p><b> figur
72、e(5)</b></p><p> subplot(2,1,1), plot(t,x3),title('輸入信號')%繪出輸入信號波形</p><p> subplot(2,1,2),plot(t,y4) %繪出輸出信號波形</p><p> title('輸出信號'),xlabel
73、('時間/s')</p><p><b> %矩形窗 %帶通</b></p><p> wls=0.25*pi;wlp=0.45*pi;whs=0.625*pi;whp=0.7*pi;</p><p> delta_w=whp-whs;</p><p> wp=[0.35 0.66];N=48;
74、 %通帶邊界頻率(歸一化頻率)和濾波器階數(shù)</p><p><b> Fs=8000;</b></p><p> b=fir1(N,wp); %設(shè)計FIR帶通濾波器</p><p><b> figure(6)</b></p><p> [H,f]=freqz(b,1,
75、512,Fs); %以50Hz為采樣頻率求出濾波器頻率響應(yīng)</p><p> subplot(2,1,1),plot(f,20*log10(abs(H)))</p><p> xlabel('頻率/Hz');ylabel('振幅/dB');grid on;</p><p> subplot(2,1,2),plot(
76、f,180/pi*unwrap(angle(H)))</p><p> xlabel('頻率/Hz');ylabel('相位/^o');grid on;</p><p> y5=fftfilt(b,x4); %采用fftfilt對輸入信號濾波</p><p> %sound(y5,fs,bits);<
77、/p><p><b> figure(7)</b></p><p> subplot(2,1,1), plot(t,x4),title('輸入信號')%繪出輸入信號波形</p><p> subplot(2,1,2),plot(t,y5) %繪出輸出信號波形</p><p&
78、gt; hold on;plot(N/2*0.02*ones(1,2),ylim, 'r') %繪制延遲到的時刻</p><p> title('輸出信號'),xlabel('時間/s')</p><p><b> %哈明窗 %低通</b></p><p> wp=0.5*pi;ws
79、=0.6*pi;</p><p><b> Bt=ws-wp;</b></p><p> N=ceil(6.6*pi/Bt);</p><p> wc=(wp+ws)/2/pi;</p><p> hn=fir1(N,wc,hamming(N+1));</p><p> [H,f]=fr
80、eqz(hn,1,512,8000); %采用8000Hz的采樣頻率求出頻率響應(yīng)</p><p><b> figure(8)</b></p><p> subplot(2,1,1),plot(f,20*log10(abs(H)))</p><p> xlabel('頻率/Hz');ylabel('振幅/dB
81、9;);grid on;</p><p> subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))</p><p> xlabel('頻率/Hz');ylabel('相位/^o');grid on;</p><p> y6=fftfilt(hn,x3); %采用fftf
82、ilt對輸入信號濾波</p><p> %sound(y6,fs,bits);</p><p><b> figure(9)</b></p><p> subplot(2,1,1), plot(t,x3),title('輸入信號')%繪出輸入信號波形</p><p> subplot(2,1,2)
83、,plot(t,y6) %繪出輸出信號波形</p><p> title('輸出信號'),xlabel('時間/s')</p><p><b> %哈明窗 %帶通</b></p><p><b> %Samp7_6</b></p><
84、;p> wp=[0.35 0.65];N=48; %通帶邊界頻率(歸一化頻率)和濾波器階數(shù)</p><p><b> Fs=8000;</b></p><p> b=fir1(N,wp); %設(shè)計FIR帶通濾波器</p><p> figure(10)</p><p> [H,f
85、]=freqz(b,1,512,Fs); %以50Hz為采樣頻率求出濾波器頻率響應(yīng)</p><p> subplot(2,1,1),plot(f,20*log10(abs(H)))</p><p> xlabel('頻率/Hz');ylabel('振幅/dB');grid on;</p><p> subplot(
86、2,1,2),plot(f,180/pi*unwrap(angle(H)))</p><p> xlabel('頻率/Hz');ylabel('相位/^o');grid on;</p><p> y7=fftfilt(b,x4); %采用fftfilt對輸入信號濾波</p><p> sound(y7,fs
87、,bits);</p><p> figure(11)</p><p> subplot(2,1,1), plot(t,x4),title('輸入信號')%繪出輸入信號波形</p><p> subplot(2,1,2),plot(t,y7) %繪出輸出信號波形</p><p> ho
88、ld on;plot(N/2*0.02*ones(1,2),ylim, 'r') %繪制延遲到的時刻</p><p> title('輸出信號'),xlabel('時間/s')</p><p><b> 7、評分表</b></p><p> 計算機與通信學(xué)院課程設(shè)計評分表</p&
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 眾賞文庫僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 數(shù)字信號處理課程設(shè)計--數(shù)字濾波器
- 數(shù)字信號處理課程設(shè)計--數(shù)字濾波器
- 數(shù)字信號處理課程設(shè)計---數(shù)字濾波器設(shè)計及濾波
- 數(shù)字信號課程設(shè)計---fir數(shù)字濾波器
- 數(shù)字信號處理課程設(shè)計---數(shù)字濾波器的設(shè)計及實現(xiàn)
- 數(shù)字信號處理課程設(shè)計--fir數(shù)字濾波器的設(shè)計
- 數(shù)字信號處理課程設(shè)計-- iir數(shù)字濾波器的設(shè)計
- 數(shù)字信號處理課程設(shè)計報告--- 數(shù)字濾波器的設(shè)計
- 數(shù)字信號課程設(shè)計報告--數(shù)字濾波器設(shè)計
- 數(shù)字信號處理課程設(shè)計--有限沖擊響應(yīng)數(shù)字濾波器設(shè)計
- 語音信號去噪 ——數(shù)字濾波器的設(shè)計
- 數(shù)字信號處理課程設(shè)計-橢圓高通iir數(shù)字濾波器設(shè)計
- 數(shù)字信號處理課程設(shè)計報告--fir數(shù)字濾波器設(shè)計及matlab實現(xiàn)
- 數(shù)字信號處理課程設(shè)計--基于dsp 的 fir 數(shù)字濾波器的設(shè)計
- 基于matlab數(shù)字濾波器設(shè)計——語音信號去噪
- 數(shù)字信號處理課程設(shè)計-- 用窗函數(shù)法設(shè)計fir數(shù)字濾波器
- 用窗函數(shù)法設(shè)計fir數(shù)字濾波器-數(shù)字信號處理課程設(shè)計
- 數(shù)字信號處理課程設(shè)計--語音信號的處理與濾波
- 數(shù)字信號處理課程設(shè)計--語音信號濾波去噪
- 數(shù)字信號課程設(shè)計---語音信號的處理與濾波
評論
0/150
提交評論