附頁第頁陜西科技大學實驗報告第頁班級信工082學號16姓名劉剛實驗組別實驗日期室溫報告日期成績報告內容：(目的和要求，原理，步驟，數(shù)據(jù)，計算，小結等)實驗七語音處理系統(tǒng)仿真（綜合濾波器設計）一代碼：（1）主程序：x=wavread('login.wav');n=0:1/22050:(size(x)-1)/22050;x_dft=fft(x);%DFT[x_dtftw]=freqz(x);%DTFTsubplot(2,2,1);plot(n,x);title('原語音信號x1');subplot(2,2,2);plot(n,abs(x_dft));title('x的DFT變換結果x2');gridon;subplot(2,2,3);plot(w/pi,abs(x_dtft));title('x2的幅度相位響應');gridon;subplot(2,2,4);plot(w/pi,unwrap(angle(x_dtft)));title('x2的相位相位響應');gridon;%用FIR濾波器進行濾波[y_lowb_low]=lowpass_fir_filter(x);%調用FIR低通濾波器函數(shù)[y_bandb_band]=bandpass_fir_filter(x);%調用FIR帶通濾波器函數(shù)[y_highb_high]=highpass_fir_filter(x);%調用FIR高通濾波器函數(shù)[hh_lowww_low]=freqz(b_low,1,512);[hh_bandww_band]=freqz(b_band,1,512);[hh_highww_high]=freqz(b_high,1,512);[h_loww_low]=freqz(y_low);[h_bandw_band]=freqz(y_band);[h_highw_high]=freqz(y_high);figure;subplot(3,3,1);plot(ww_low/pi,20*log10(abs(hh_low)));title('FIR低通濾波器');gridon;subplot(3,3,2);plot(n,y_low);title('經(jīng)FIR低通濾波后的信號x1');gridon;subplot(3,3,3);plot(w_low/pi,abs(h_low));title('x1的幅度相位響應');gridon;subplot(3,3,4);plot(ww_band/pi,20*log10(abs(hh_band)));title('FIR帶通濾波器');gridon;subplot(3,3,5);plot(n,y_band);title('經(jīng)FIR帶通濾波后的信號x2');gridon;subplot(3,3,6);plot(w_band/pi,abs(h_band));title('x2的幅度相位響應');gridon;subplot(3,3,7);plot(ww_high/pi,20*log10(abs(hh_high)));title('FIR高通濾波器');gridon;subplot(3,3,8);plot(n,y_high);title('經(jīng)FIR高通濾波后的信號x3');gridon;subplot(3,3,9);plot(w_high/pi,abs(h_high));title('x3的幅度相位響應');gridon;%用IIR濾波器進行濾波[y_lowbz_lowaz_low]=lowpass_iir_filter(x);%調用IIR低通濾波器函數(shù)[y_bandbz_bandaz_band]=bandpass_iir_filter(x);%調用IIR帶通濾波器函數(shù)[y_highbz_highaz_high]=highpass_iir_filter(x);%調用IIR高通濾波器函數(shù)[hh_lowww_low]=freqz(bz_low,az_low,256,1);[hh_bandww_band]=freqz(bz_band,az_band,256,1);[hh_highww_high]=freqz(bz_high,az_high,256,1);[h_loww_low]=freqz(y_low);[h_bandw_band]=freqz(y_band);[h_highw_high]=freqz(y_high);figure(3);subplot(3,3,1);plot(ww_low,abs(hh_low));title('IIR低通濾波器');gridon;subplot(3,3,2);plot(n,y_low);title('經(jīng)IIR低通濾波后的信號x1');gridon;subplot(3,3,3);plot(w_low/pi,abs(h_low));title('x1的幅度相位響應');gridon;subplot(3,3,4);plot(ww_band,abs(hh_band));title('IIR帶通濾波器');gridon;subplot(3,3,5);plot(n,y_band);title('經(jīng)IIR帶通濾波后的信號x2');gridon;subplot(3,3,6);plot(w_band/pi,abs(h_band));title('x2的幅度相位響應');gridon;subplot(3,3,7);plot(ww_high,abs(hh_high));title('IIR高通濾波器');gridon;subplot(3,3,8);plot(n,y_high);title('經(jīng)IIR高通濾波后的信號x3');gridon;subplot(3,3,9);plot(w_high/pi,abs(h_high));title('x3的幅度相位響應');gridon;（1）主程序：1)FIR低通濾波器function[yb]=bandpass_fir_filter(x)w1=0.4*pi;w2=0.7*pi;wn=[w1w2];wdelta=w2-w1;M=ceil(3.32*pi/wdelta);N=2*M+1;win=hamming(N);b=fir1(N-1,wn/pi,'bandpass',win);y=filter(b,1,x);2)FIR高通濾波器function[yb]=highpass_fir_filter(x)wp=0.6*pi;ws=0.9*pi;wdelta=ws-wp;M=ceil(3.32*pi/wdelta);N=2*M+1;wc=(ws+wp)/2;win=hamming(N);b=fir1(N-1,wc/pi,'high',win);y=filter(b,1,x);3)FIR帶通濾波器function[yb]=lowpass_fir_filter(x)wp=0.22*pi;ws=0.46*pi;wdelta=ws-wp;M=ceil(3.32*pi/wdelta);N=2*M+1;%′°wc=(ws+wp)/2;win=hamming(N);b=fir1(N-1,wc/pi,win);y=filter(b,1,x);4)IIR低通濾波器function[ybzaz]=lowpass_iir_filter(x)rp=1;rs=15;wp=0.25*pi;ws=0.4*pi;wap=tan(wp/2);was=tan(ws/2);[n,wn]=buttord(wap,was,rp,rs,'s');[z,p,k]=buttap(n);[bp,ap]=zp2tf(z,p,k);[bs,as]=lp2lp(bp,ap,wap);[bz,az]=bilinear(bs,as,1/2);y=filter(bz,az,x);5)IIR高通濾波器function[ybzaz]=highpass_iir_filter(x)rp=1;rs=15;wp=0.7*pi;ws=0.85*pi;wap=tan(wp/2);was=tan(ws/2);[n,wn]=buttord(wap,was,rp,rs,'s');[z,p,k]=buttap(n);[bp,ap]=zp2tf(z,p,k);[bs,as]=lp2hp(bp,ap,wap);[bz,az]=bilinear(bs,as,1/2);y=filter(bz,az,x);6)IIR帶通濾波器function[ybzaz]=bandpass_iir_filter(x)rp=1;rs=15;wp=0.35*