【正文】 N = length (y) 。 %求出語音信號的長度noise=*randn(N,1)。 %隨機函數(shù)產(chǎn)生噪聲Si=y+noise。 %語音信號加入噪聲 sound(Si)。subplot(2,1,1)。plot(Si)。title(39。加噪語音信號的時域波形39。)。xlabel(39。時間39。)。ylabel(39。幅度39。)。S=fft(Si)。 %傅里葉變換subplot(2,1,2)。plot(abs(S))。 title(39。加噪語音信號的頻域波形39。)。 xlabel(39。時間39。)。ylabel(39。幅度39。)。程序結(jié)果如下圖：4 數(shù)字濾波器設(shè)計 用窗函數(shù)法設(shè)計IIR帶通濾波器程序如下：wp=*pi。 ws=*pi。 %性能指標wdelta=wswp。 %過渡帶寬度M=ceil(*pi/wdelta)。 %濾波器長度，朝正無窮方向舍入N=2*M+1。 %窗口長度wc=(ws+wp)/2。 %截止頻率win=hamming(N)。 %因為衰減不小于50dB，所以選擇海明窗，這里得到海明窗的時域響應(yīng)b=fir1(N1,wc/pi,win)。n=0:1:N。[hi t]=impz(b,1,n)。%得到脈沖響應(yīng)[hf w]=freqz(b,1,512)。 %得到頻率響應(yīng)subplot(3,1,1)。 stem(n,hi)。xlabel(39。n39。)。 ylabel(39。幅度39。)。 title(39。脈沖響應(yīng)39。)。subplot(3,1,2)。 plot(w/pi,20*log10(abs(hf)))。xlabel(39。頻率(Hz)39。)。 ylabel(39。幅度(dB)39。)。title(39。頻率響應(yīng)39。)。subplot(3,1,3)。 plot(w/pi,180/pi*unwrap(angle(hf)))。xlabel(39。頻率(Hz)39。)。 ylabel(39。相位(degrees)39。)。title(39。頻率響應(yīng)39。)。程序結(jié)果如下圖： 用窗函數(shù)法設(shè)計FIR低通濾波器程序如下：Ft=8000。 Fp=1000。 Fs=1200。wp=2*Fp/Ft。ws=2*Fs/Ft。rp=1。rs=50。p=110.^(rp/20)。s=10.^(rs/20)。fpts=[wp,ws]。mag=[1,0]。dev=[p,s]。[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev)。b21=fir1(n21,wn21,kaiser(n21+1,beta))。[h,w]=freqz(b21,1)。plot(w/pi,abs(h))。title(39。FIR低通濾波器39。)。grid。其程序結(jié)果如下圖：5信號處理用設(shè)計好的數(shù)字濾波器對含噪聲的語音信號進行濾波，在MATLAB中FIR濾波器利用函數(shù)fftfilt對信號進行濾波，IIR濾波器利用函數(shù)filter對信號進行濾波，在一個窗口同時畫出濾波前后的波形及頻譜如下圖。 IIR帶通濾波程序如下：[y,fs,nbits]=wavread(39。39。)。N=length(y)。Noise=*randn(N,1)。Si=y+Noise。 %語音信號加入噪聲 nn=0:1:50。y=filter(b,1,Si)。subplot(2,1,1)。 stem(Si)。xlabel(39。n39。)。 ylabel(39。幅度39。)。title(39。濾波前信號頻譜39。)。 axis([0 50 2 2])。subplot(2,1,2)。 stem(out)。xlabel(39。n39。)。 ylabel(39。幅度39。)。title(39。濾波后信號頻譜39。)。 axis([0 50 1 1])。sound(y)。程序結(jié)果如下圖： FIR低通濾波程序如下：[y,fs,nbits]=wavread(39。39。)。N = length (y) 。 noise=*randn(N,1)。 Si=y+noise。 S=fft(Si)。z21=fftfilt(b21,Si)。sound(z21)。m21=fft(z21)。subplot(2,2,1)。plot(abs(S),39。g39。)。title(39。濾波前信號的頻譜39。)。xlabel(39。時間39。)。ylabel(39。幅度39。)。grid。subplot(2,2,2)。plot(abs(m21),39。r39。)。title(39。濾波后信號頻譜39。)。xlabel(39。時間39。)。ylabel(39。幅度39。)。grid。subplot(2,2,3)。plot(S)。title(39。濾波前信號波形39。)。xlabel(39。時間39。)。ylabel(39。幅度39。)。grid。subplot(2,2,4)。plot(z21)。title(39。濾波后信號波形39。)。xlabel(39。時間39。)。ylabel(39。幅度39。)。grid。程序結(jié)果如下圖：6心得體會通過這次課程設(shè)計，我學到了很多東西。進一步加深了對MATLAB的了解。在設(shè)計的過程中我遇到了不少的問題，比如剛開始要處理的是WAV格式的語音信號，可是我發(fā)現(xiàn)不是所有的WAV文件都可以用WAVREAD函數(shù)讀取，通過軟件轉(zhuǎn)化了才可以使用。在編寫程序的時候遇到了很多錯誤，最后經(jīng)過查找一些資料和自己的理解都一一糾正了。總的來說，通過這次的課程設(shè)計我對語音信號有的更多的認識，對數(shù)字信號處理的知識有了更多的理解，鍛煉了我的理論與實際相結(jié)合的能力，并且熟練的應(yīng)用MATLAB也加深了我對課程的理解。這次設(shè)計使我更深刻的了解MATLAB的使用方法，同時我相信，進一步加強對MATLA