數(shù)字信號處理課程設(shè)計報告--基于matlab的語音去噪處理-資料下載頁
2025-01-21 17:21本頁面
【正文】 (39。幅值A(chǔ)39。)subplot(2,1,2)。plot(f,abs(y_zp(1:n/2)))。title(39。加噪語音信號頻譜圖39。)。xlabel(39。頻率Hz39。)。ylabel(39。頻率幅值39。)。附錄3[y,fs,bits]=wavread(C:\Users\Administrator\Desktop\數(shù)字信號處理課程設(shè)計\39。)。sound(y,fs)n=length(y)y_p=fft(y,n)。f=fs*(0:n/21)/n。figure(1)subplot(2,1,1)。plot(y)。title(39。原始語音信號采樣后的時域波形39。)。xlabel(39。時間軸39。)ylabel(39。幅值A(chǔ)39。)subplot(2,1,2)。plot(f,abs(y_p(1:n/2)))。title(39。原始語音信號采樣后的頻譜圖39。)。xlabel(39。頻率Hz39。)。ylabel(39。頻率幅值39。)。L=length(y)noise=*randn(L,2)。y_z=y+noise。sound(y_z,fs)n=length(y)。y_zp=fft(y_z,n)。f=fs*(0:n/21)/n。figure(2)subplot(2,1,1)。plot(y_z)。title(39。加噪語音信號時域波形39。)。xlabel(39。時間軸39。)ylabel(39。幅值A(chǔ)39。)subplot(2,1,2)。plot(f,abs(y_zp(1:n/2)))。title(39。加噪語音信號頻譜圖39。)。xlabel(39。頻率Hz39。)。ylabel(39。頻率幅值39。)。Ft=5000。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。)。x=fftfilt(b21,y_z)。X=fft(x,n)。figure(4)。subplot(2,2,1)。plot(f,abs(y_zp(1:n/2)))。title(39。濾波前信號的頻譜39。)。subplot(2,2,2)。plot(f,abs(X(1:n/2)))。title(39。濾波后信號的頻譜39。)。subplot(2,2,3)。plot(y_z)。title(39。濾波前信號的時域波形39。)subplot(2,2,4)。plot(x)。title(39。濾波后信號的時域波形39。)sound(x,fs,bits)附錄4Ft=8000。Fp=1000。Fs=1200。wp=2*pi*Fp/Ft。ws=2*pi*Fs/Ft。fp=2*Ft*tan(wp/2)。fs=2*Fs*tan(wp/2)。[n11,wn11]=buttord(wp,ws,1,50,39。s39。)。 %求低通濾波器的階數(shù)和截止頻率[b11,a11]=butter(n11,wn11,39。s39。)。 %求S域的頻率響應(yīng)的參數(shù)[num11,den11]=bilinear(b11,a11,)。 %雙線性變換實現(xiàn)S域到Z域的變換[h,w]=freqz(num11,den11)。 %根據(jù)參數(shù)求出頻率響應(yīng)plot(w*8000*,abs(h))。legend(39。用切比雪夫lchoby1設(shè)計39。)。grid。[y,fs,nbits]=wavread (39。C:\Users\Administrator\Desktop\數(shù)字信號處理課程設(shè)計\39。)。n = length (y) 。 %求出語音信號的長度noise=*randn(n,2)。 %隨機函數(shù)產(chǎn)生噪聲s=y+noise。 %語音信號加入噪聲S=fft(s)。 %傅里葉變換z11=filter(num11,den11,s)。sound(z11)。m11=fft(z11)。 %求濾波后的信號subplot(2,2,1)。plot(abs(S),39。g39。)。title(39。濾波前信號的頻譜39。)。grid。subplot(2,2,2)。plot(abs(m11),39。r39。)。title(39。濾波后信號的頻譜39。)。grid。subplot(2,2,3)。plot(s)。title(39。濾波前信號的波形39。)。grid。subplot(2,2,4)。plot(z11)。title(39。濾波后的信號波形39。)。第23頁
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1