語音信號處理與濾波課程設(shè)計(jì)-畢業(yè)設(shè)計(jì)-在線瀏覽
2025-03-24 02:12本頁面
【正文】 ead(‘e:\’)。 y=fft(x1,1024)。 plot(x1)。 原始語音采樣后的的時(shí)域信號波形圖： （ 2）對語音信號進(jìn)行頻譜分析： 在 MATLAB 中，利用函數(shù) fft 對信號進(jìn)行快速傅里葉變換，得到信號的頻譜特性。 程序代碼： [x1,fs,bits]=wavread(‘D:\’)。 f=fs*(0:511)/1024。 plot(f,abs(1:512)))。 xlabel(39。)。幅值 39。 subplot(2,1,2)。 title(‘原始語音信號 FFT 頻譜’ )。點(diǎn)數(shù) N39。 ylabel(39。)。 首先用窗函數(shù)法（ 矩形窗 (Rectangular window)、三角窗 (Triangular window)、漢寧窗 (Hanning window)、海明窗 (Hamming window)、布拉克曼窗 (Blackman window)、切比雪夫窗 (Chebyshev window)、巴特里特窗 (Bartlett window)及凱塞窗 (Kaiser window)。在 MATLAB 中，利用函數(shù) fir1設(shè)計(jì) FIR 濾波器；然后再用雙線性變換法設(shè)計(jì)上面要求的三種濾波器，利用函數(shù)butte、 cheby1 和 ellip 設(shè)計(jì) IIR 濾波器；最后，利用 MATLAB 中的函數(shù) freqz 畫出各種濾波器的頻率響應(yīng)。 fc=1200。 Ap=1。 wc=2*fc/fs。 N=ceil（ ()/(*(wcwp)/2))+1。 Win=Kaiser(N+1,beta)。 freqz(b,1,512,fs)。 fc=1200。 Ap=1。 wc=2*fc/fs。 [n,wn]=ellipord(wp,wc,Ap,As)。 freqz(b,a,512,fs)。 fc=4800。 Ap=1。 wc=2*fc/fs。 [n,wn]=ellipord(wp,wc,Ap,As)。 freqz(b,a,512,fs)。fp2=3000。 fc2=3200。 Ap=1。 wc=[2*fc1/fs,2*fc2/fs]。 [n,wn]=ellipord(wp,wc,Ap,As)。 freqz(b,a,512,fs)。在MATLAB 中， FIR 濾波器利用函數(shù) fftfilt 對信號進(jìn)行濾波， IIR 濾波器利用函數(shù)filter 對信號進(jìn)行濾波。其調(diào)用格式為 sound(x,fs,bits)，可以感覺濾波前后的聲音有變化。 界面如下圖所示： 四、調(diào)試分析 A．用窗函數(shù)法設(shè)計(jì)低通濾波器中，語句 N=ceil()/(*(wcwp)/2))+1。即可。s39。改成wp=[2*fp1/fs,2*fp2/fs]。 濾波前后的頻譜分析：由于所取的采樣點(diǎn)數(shù)比較大，濾波前后的頻譜比較相近，但仔細(xì)對比仍然可以看出下圖的結(jié)果濾掉了高頻的成分。 濾波前后的頻譜分析 :由圖明顯可以看出，此濾波器將低頻成分濾掉了。 濾波前后的頻譜分析：該濾波器將低頻成分和高頻成分都濾掉了，只剩下通帶里的頻率。在雙線性變換法的設(shè)計(jì)中選用橢圓函數(shù)設(shè)計(jì)，濾波器的幅度和相位均滿足設(shè)計(jì)的要求。我深深的體會(huì)到了知識的重要性，不管是為了應(yīng)付考試還是怎么樣也好，對于自身的幫助還是很大的，整整兩天坐在電腦面前的功夫沒有白費(fèi)，我還是收獲了很多，對于 Matlab 軟件的使用也更加熟悉了。學(xué)會(huì)團(tuán)結(jié)，才能把事情做得更好。西安：西安電子科技大學(xué)出版社， 2021 [2]懷琛 數(shù)字信號處理教程 —— MATLAB釋疑與實(shí)現(xiàn)。北京：清華大學(xué)出版社， 2021 [4]劉順蘭，吳杰 數(shù)字信號處理。請選擇選項(xiàng) 39。原始信號采樣后時(shí)域圖和頻譜圖 39。FIR 濾波器 39。IIR 濾波器39。退出 39。 if b==4 b==0。請選擇選項(xiàng) 39。播放原始語音 39。原始語音時(shí)域圖和頻率響應(yīng)圖 39。原始語音頻譜圖及做 FFT變換后頻譜圖 39。返回 39。 if temp==1 voice。 else if temp==2 pingpu。 else main。 main。 main。 main。 end else if b==3 temp=menu(39。,39。,39。,39。,39。)。 main。 main。 main。 end end fs=22050。D:\39。 wp1=2*pi*1200/fs。 ws1=2*pi*1000/fs。 Rp=1。 wn=(wp1+ws1)/2。 wdelta=wp1ws1。 wn=[wp,ws]。ba
試題試卷相關(guān)推薦
鄂ICP備17016276號-1