武漢理工大學(xué)基于矩形窗、三角窗、海明窗、漢寧窗、布拉克曼窗的FIR數(shù)字濾波器設(shè)計(jì)

基于矩形窗的FIR數(shù)字低通濾波器設(shè)計(jì)參數(shù)如下：主程序：%FIR數(shù)字低通濾波器設(shè)計(jì)wp=2*pi*1.75/15;ws=2*pi*3.25/15w=(wp+ws)/2;N=ceil(4*pi/(ws-wp));n=0:1:(N-1);alpha=(N-1)/2;m=n-alpha+eps;hd=sin(w*m)./(pi*m);%得到抱負(fù)低通濾波器B=boxcar(N);h=hd.*(B)”;%FIR數(shù)字濾波器[H,w]=freqz(h,[1],1024);%求其頻率響應(yīng)db=20*log10(abs(H)+eps);%得到幅值pha=angle(H);%得到相位figure(1);subplot(2,2,1);axis([0,N-1,-0.1,0.35]);實(shí)際低通濾波器的text((0.3*N),0.275,string);subplot(2,2,2);axis([0,1,-100,5]);grid;subplot(2,2,3);plot(w,pha);holdon;”””axis([0,3.2,-4,4]);subplot(2,2,4);plot(w,abs(H));””頻率”axis([0,3.15,0,1.5]);text(0.9,1.3,string);基于矩形窗的帶通濾波器的設(shè)計(jì)參數(shù)如下：子程序：functionhd=ideal_bs(Wcl,Wch,N)alpha=(N-1)/2;n=0:1:N-1;m=n-alpha+eps;hd=[sin(Wch*m)-sin(Wcl*m)]./(pi*m);functionhd=ideal_bs(Wcl,Wch,N)alpha=(N-1)/2;n=0:1:N-1;m=n-alpha+eps;hd=[sin(Wch*m)-sin(Wcl*m)]./(pi*m);主程序：clearall;Wph=2*pi*4.75/15;Wpl=2*pi*2.5/15;Wsl=2*pi/15;Wsh=2*pi*6.25/15;tr_width=min((Wpl-Wsl),(Wsh-Wph)); %過(guò)渡帶寬度N=ceil(4*pi/tr_width); %濾波器長(zhǎng)度n=0:1:N-1;Wcl=(Wsl+Wpl)/2; %抱負(fù)濾波器的截止頻率Wch=(Wsh+Wph)/2;hd=ideal_bs(Wcl,Wch,N); %抱負(fù)濾波器的單位沖擊響應(yīng)w_ham=(boxcar(N))”;string=[”矩形窗”,”N=”,num2str(N)];h=hd.*w_ham; %截取取得實(shí)際的單位脈沖響應(yīng)[db,mag,pha,w]=freqz_m2(h,[1]); %計(jì)算實(shí)際濾波器的幅度響應(yīng)subplot(3,2,1);stem(n,hd);title(”抱負(fù)脈沖響應(yīng)hd(n)”)axis([-1,N,-0.5,0.8]);xlabel(”n”);ylabel(”hd(n)”);subplot(3,2,2);stem(n,w_ham);axis([-1,N,0,1.1]);xlabel(”n”);ylabel(”w(n)”);text(1.5,1.3,string);subplot(3,2,3);stem(n,h);title(”實(shí)際脈沖響應(yīng)h(n)”);axis([0,N,-1.4,1.4]);xlabel(”n”);ylabel(”h(n)”);subplot(3,2,4);plot(w,pha);title(”相頻特性”);axis([0,3.15,-4,4]);xlabel(”頻率〔rad〕”);ylabel(”相位〔?〕”);subplot(3,2,5);plot(w/pi,db);title(”幅度特性〔dB〕”);axis([0,1,-80,10]);xlabel(”頻率〔pi〕”);ylabel(”分貝數(shù)”);subplot(3,2,6);plot(w,mag);title(”頻率特性”)axis([0,3.15,0,1.5]);xlabel(”頻率〔rad〕”);ylabel(”幅值”);fs=15000t=(0:100)/fsx=sin(2*pi*t*750)+sin(2*pi*t*5000)+sin(2*pi*t*6100)q=filter(h,1,x)[a,f1]=freqz(x)f1=f1/pi*fs/2[b,f2]=freqz(q)f2=f2/pi*fs/2figure(2)subplot(2,1,1)plot(f1,abs(a))title(”輸入波形頻譜圖”)xlabel(”頻率”)ylabel(”幅度”)subplot(2,1,2)plot(f2,abs(b))title(”輸出波形頻譜圖”)xlabel(”頻率”)ylabel(”幅度”)基于矩形窗的FIR數(shù)字帶阻濾波器設(shè)計(jì)參數(shù)：主程序：clearall;Wph=2*pi*6.25/15;Wpl=2*pi/15;Wsl=2*pi*2.5/15;Wsh=2*pi*4.75/15;tr_width=min((Wsl-Wpl),(Wph-Wsh));%過(guò)渡帶寬度N=ceil(4*pi/tr_width);%濾波器長(zhǎng)度n=0:1:N-1;Wcl=(Wsl+Wpl)/2;%抱負(fù)濾波器的截止頻率Wch=(Wsh+Wph)/2;hd=ideal_bs(Wcl,Wch,N);%抱負(fù)濾波器的單位沖擊響應(yīng)w_ham=(boxcar(N))”;string=[”矩形窗”,”N=”,num2str(N)]h=hd.*w_ham;%截取取得實(shí)際的單位脈沖響應(yīng)[db,mag,pha,w]=freqz_m2(h,[1]);%計(jì)算實(shí)際濾波器的幅度響應(yīng)delta_w=2*pi/1000;subplot(3,2,1);stem(n,hd);title(hd(n)”)axis([-1,N,-0.5,0.8]);xlabel(”n”);ylabel(”hd(n)”);subplot(3,2,2);stem(n,w_ham);axis([-1,N,0,1.1]);xlabel(”n”);ylabel(”w(n)”);text(1.5,1.3,string);subplot(3,2,3);stem(n,h);title(”實(shí)際脈沖響應(yīng)h(n)”)axis([0,N,-1.4,1.4]);xlabel(”n”);ylabel(”h(n)”);subplot(3,2,4);plot(w,pha);title(”相頻特性”)axis([0,3.15,-4,4]);xlabel(”頻率〔rad〕”)ylabel(”相位〔Φ〕”)subplot(3,2,5);plot(w/pi,db);title(”幅度特性〔dB〕”)axis([0,1,-80,10]);xlabel(”頻率〔pi〕”)ylabel(”分貝數(shù)”)subplot(3,2,6);plot(w,mag);title(”頻率特性”)axis([0,3,0,2]);xlabel(”頻率〔rad〕”)ylabel(”幅值”)fs=15000;t=(0:100)/fs;x=sin(2*pi*t*750)+sin(2*pi*t*3000)+sin(2*pi*t*6100);q=filter(h,1,x);[a,f1]=freqz(x);f1=f1/pi*fs/2;[b,f2]=freqz(q);f2=f2/pi*fs/2;figure(2);subplot(2,1,1);plot(f1,abs(a));title(”輸入波形頻譜圖”)xlabel(”頻率”)ylabel(”幅度”)subplot(2,1,2);plot(f2,abs(b));title(”輸出波形頻譜圖”)xlabel(”頻率”)ylabel(”幅度”)；1：functionhd=ideal_bs(Wcl,Wch,m);alpha=(m-1)/2;n=[0:1:(m-1)];m=n-alpha+eps;hd=[sin(m*pi)+sin(Wcl*m)-sin(Wch*m)]./(pi*m)調(diào)用程序2：function[db,mag,pha,w]=freqz_m2(b,a)[H,w]=freqz(b,a,1000,”whole”);H=(H(1:1:501))”;w=(w(1:1:501))”;mag=abs(H);db=20*log10((mag+eps)/max(mag));pha=angle(H);基于三角窗的FIR數(shù)字高通濾波器設(shè)計(jì)參數(shù)：子程序：functionhd=ideal_hp1(Wc,N)alpha=(N-1)/2;n=0:1:N-1;m=n-alpha+eps;%hd=[sin(pi*m)-sin(Wc*m)]/(pi*m);hd=[sin(pi*m)-sin(Wc*m)]./(pi*m);function[db,mag,pha,w]=freqz_m2(b,a)[H,w]=freqz(b,a,1000,”whole”);H=(H(1:1:501))”;w=(w(1:1:501))”;mag=abs(H);db=20*log10((mag+eps)/max(mag));pha=angle(H);主程序：clearall;Ws=6.5/15*piWp=3.5/15*pitr_width=Ws-Wp;%過(guò)渡帶寬度n=0:1:N-1;Wc=(Ws+Wp)/2;%抱負(fù)濾波器的截止頻率hd=ideal_hp1(Wc,N);%抱負(fù)濾波器的單位沖擊響應(yīng)w_ham=(triang(N))”;string=[”三角窗”,”N=”,num2str(N)];h=hd.*w_ham;%截取得到的實(shí)際的單位脈沖響應(yīng)delta_w=2*pi/1000;Ap=-(min(db(Ws/delta_w+1:1:501)));%實(shí)際通帶波浪hd(n)”);axis([0,N-1,-0.5,0.5]);xlabel(”n”);ylabel(”hd(n)”);subplot(3,2,2);stem(n,w_ham);axis([0,N-1,0,1.1]);xlabel(”n”);ylabel(”w(n)”);text(1.5,1.3,string);subplot(3,2,3);stem(n,h);title(”實(shí)際脈沖響應(yīng)hd(n)”);axis([0,N,-1.4,1.4]);xlabel(”n”);ylabel(”hd(n)”);subplot(3,2,4);plot(w,pha);title(”想頻特性”);axis([0,3.15,-4,4]);xlabel(”頻率(pi)”);ylabel(”相位(o)”);subplot(3,2,5);plot(w/pi,db);title(”幅度特性(dB)”);(pi)”);ylabel(”(分貝數(shù))”);subplot(3,2,6);plot(w,mag);title(”頻率特性”);(rad)”);ylabel(”幅值”);fs=12023;t=(0:100)/fs;x=sin(2*pi*t*750)+sin(2*pi*t*1500)+sin(2*pi*t*3000);q=filter(h,1,x);[a,f1]=freqz(x);f1=f1/pi*fs/2;[b,f2]=freqz(q);f2=f2/pi*fs/2;figure(2);subplot(2,1,1);plot(f1,abs(a));title(”輸入波形頻譜圖”);”);ylabel(”幅度”)subplot(2,1,2);plot(f2,abs(b));