基于matlab的fir濾波器語音信號去噪畢業(yè)設計-資料下載頁
2025-06-27 18:17本頁面
【正文】 ot(x)。 %做原始語音信號的時域圖形title(39。原始語音信號的時域圖39。)subplot(212)。 plot(x1)。 %繪制加噪語音信號的時域圖形title(39。加噪語音信號的時域圖39。)。 xlabel(39。t39。)。ylabel(39。x139。)。figure(2)。y1=fft(x1,1024)。 %對加噪信號做1024點FFT變換f=fs*(0:511)/512。subplot(211)。plot(f,abs(y(1:512)))。title(39。原始語音信號頻譜圖39。)。xlabel(39。Hz39。)。ylabel(39。fudu39。)。grid on。Subplot(212)。plot(f,abs(y1(1:512)))。 %繪制原始語音信號頻率響應圖title(39。加噪語音信號頻譜圖39。)。xlabel(39。Hz39。)。ylabel(39。fudu39。)。grid on。%語音信號濾波去噪程序設計fs=22050。 %語音信號采集頻率為22050[x,fs,bits]=wavread(39。D:\39。)。 %加隨機噪聲L=length(x)。noise=*randn(L,2)。x1=x+noise。 %利用kaiser濾波器對語音信號濾波168。fp=1000。fm=1200。rs=100。Fs=8000。 %濾波器設計wp=2*pi*fp/Fs。ws=2*pi*fm/Fs。Bt=wswp。 %計算過渡帶寬度alph=*()。%計算kaiser窗的控制參數alphM=ceil((rs8)/)。%計算kaiser窗所需階數Mwc=(wp+ws)/2/pi。hn=fir1(M,wc,kaiser(M+1,alph))。%調用kaiser計算低通FIDF的h(n)figure(1)。freqz(hn)。y=fft(x)。y1=fft(x1)。x2=fftfilt(hn,x1)。 %利用kaiser濾波器對語音信號濾波繪圖y=fft(x,1024)。 %對原始信號做1024點FFT變換y1=fft(x1,1024)。 %對加噪信號做1024點FFT變換y2=fft(x2,1024)。 %對加噪濾波后的信號做1024點FFT變換figure(1)。freqz(hn)。grid on figure(2) %畫出時域波形subplot(311)。 axis([0 7 1 1])； plot(x)。title(39。原始語音信號的時域波形39。)。subplot(312)。plot(x1)title(39。加噪后語音信號的時域波形39。)。subplot(313)。 plot(x2)。title(39。加噪語音濾波后信號的時域波形39。)。f=fs*(0:511)/512。figure(3) %畫出頻譜圖y1=fft(x1,1024)。 subplot(311)。 plot(f,abs(y(1:512)))。%畫出原始語音信號頻譜圖title(39。原始語音信號頻譜圖39。)xlabel(39。Hz39。)。ylabel(39。fuzhi39。)。subplot(312)。plot(f,abs(y1(1:512)))。 %畫出加噪后語音信號頻譜圖title(39。加噪后語音信號頻譜圖39。)xlabel(39。Hz39。)。ylabel(39。fuzhi39。)。subplot(313)。 plot(f,abs(y2(1:512)))。 %畫出濾波后的信號頻譜圖title(39。加噪語音濾波后的信號頻譜圖39。)xlabel(39。Hz39。)。ylabel(39。fuzhi39。)。致謝本次設計是在陳老師的悉心指導下完成的，在整個課程設計的完成過程中，陳老師給予了精心的指導，每一個環(huán)節(jié)都注入了導師大量的心血。導師嚴謹的治學態(tài)度、大公無私的品德和為科學獻身的精神都是我學習的榜樣，尤其是她的認真塌實、實事求是、精益求精、刻苦耐勞的工作作風使我獲益良多。學生在老師身上學到的東西不僅僅是科學知識，更多的是學到了研究的方法及做人的道理，使我終生難忘，在此向導師致以崇高的敬意和衷心的感謝!
語文相關推薦
鄂ICP備17016276號-1