2023年全國碩士研究生考試考研英語一試題真題(含答案詳解+作文范文)_第1頁
已閱讀1頁,還剩24頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權(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)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論