【正文】 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))。title(39。IIR低通濾波器39。)。legend(39。用butter設(shè)計39。)。grid。附錄4Ft=8000。Fp=4000。Fs=3500。wp1=tan(pi*Fp/Ft)。 %高通到低通濾波器參數(shù)轉(zhuǎn)換ws1=tan(pi*Fs/Ft)。wp=1。ws=wp1*wp/ws1。[n13,wn13]=cheb1ord(wp,ws,1,50,39。s39。)。 %求模擬的低通濾波器階數(shù)和截止頻[b13,a13]=cheby1(n13,1,wn13,39。s39。)。 %求S域的頻率響應(yīng)的參數(shù)[num,den]=lp2hp(b13,a13,wn13)。 %將S域低通參數(shù)轉(zhuǎn)為高通的[num13,den13]=bilinear(num,den,)。 %利用雙線性變換實現(xiàn)S域到Z域轉(zhuǎn)[h,w]=freqz(num13,den13)。plot(w*21000*,abs(h))。title(39。IIR高通濾波器39。)。legend(用cheby1設(shè)計39。)。附錄5[y,fs,bits]=wavread(39。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。幅值 A39。)subplot(2,1,2)。plot(f,abs(y_p(1:n/2)))。title(39。原始語音信號采樣后的頻譜圖39。)。xlabel(39。頻率Hz39。)。ylabel(39。頻率幅值39。)。noise=1*sin(2*pi*3000*n)。y_z=y+noise。sound(y_z,fs)L=length(y_z)。y_zp=fft(y_z,L)。f=fs*(0:L/21)/L。figure(2)subplot(2,1,1)。plot(y_z)。title(39。加噪語音信號時域波形39。)。xlabel(39。時間軸39。)ylabel(39。幅值 A39。)subplot(2,1,2)。plot(f,abs(y_zp(1:L/2)))。title(39。加噪語音信號頻譜圖39。)。xlabel(39。頻率Hz39。)。ylabel(39。頻率幅值39。)。附錄6[y,fs,bits]=wavread(39。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。)。附錄7[y,fs,bits]=wavread(39。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)附錄8Ft=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。用butter設(shè)計39。)。grid。[y,fs,nbits]=wavread (39。39。)。n = length (y) 。 %求出語音信號的長度noise=*randn(n,2)。 %隨機(jī)函數(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。)。