MATLAB作業(yè)第一題已知:h(n)={-4,1,-1,-2,5,5,-2,-1,1,-4},參照實(shí)驗(yàn)內(nèi)容5的程序或自編程序確定并繪出第二種線性相位濾波器的單位沖激響應(yīng)h(n)和振幅響應(yīng)Hr(ω).h=[-4,1,-1,-2,5,5,-2,-1,1,-4];N=length(h)n=0:N-1;[Hr,w,b,L]=Hr_Type2(h);amax=max(b)+1;;amin=min(b)-1;subplot(221);stem(n,h);axis([-12*L+1aminamax])xlabel('n');ylabel('h(n)');title('ImpulseResponse')subplot(223);stem(1:L,b);axis([-12*L+1aminamax])xlabel('n');ylabel('a(n)');title('a(n)coefficients')subplot(222);plot(w/pi,Hr);gridxlabel('Frequencyinpiunits');ylabel('Hr')title('Type-1AmplitudeResponse')subplot(224);zplane(h,1);已知:h(n)={-4,1,-1,-2,5,0,-5,2,1,-1,4},參照實(shí)驗(yàn)內(nèi)容5的程序或自編程序確定第三種線性相位濾波器的單位沖激響應(yīng)h(n)和振幅響應(yīng)Hr(ω).h=[-4,1,-1,-2,5,0,-5,2,1,-1,4];N=length(h);n=0:N-1;[Hr,w,c,L]=Hr_Type3(h);amax=max(c)+1;;amin=min(c)-1;subplot(2,2,1);stem(n,h);axis([-12*L+1aminamax])xlabel('n');ylabel('h(n)');title('ImpulseResponse')subplot(223);stem(0:L,c);axis([-12*L+1aminamax])xlabel('n');ylabel('c(n)');title('a(n)coefficients')subplot(222);plot(w/pi,Hr);gridxlabel('Frequencyinpiunits');ylabel('Hr')title('Type-1AmplitudeResponse')subplot(224);zplane(h,1);已知:h(n)={-4,1,-1,-2,5,-5,2,1,-1,4},參照實(shí)驗(yàn)內(nèi)容5的程序或自編程序確定第四種線性相位濾波器的單位沖激響應(yīng)h(n)和振幅響應(yīng)Hr(ω).h=[-4,1,-1,-2,5,-5,2,1,-1,4];N=length(h);n=0:N-1;[Hr,w,d,L]=Hr_Type4(h);amax=max(d)+1;;amin=min(d)-1;subplot(2,2,1);stem(n,h);axis([-12*L+1aminamax])xlabel('n');ylabel('h(n)');title('ImpulseResponse')subplot(223);stem(1:L,d);axis([-12*L+1aminamax])xlabel('n');ylabel('d(n)');title('a(n)coefficients')subplot(222);plot(w/pi,Hr);gridxlabel('Frequencyinpiunits');ylabel('Hr')title('Type-1AmplitudeResponse')subplot(224);zplane(h,1);第二題利用窗函數(shù)設(shè)計(jì)FIR濾波器2:參照實(shí)驗(yàn)內(nèi)容7的程序或自編程序來實(shí)現(xiàn)設(shè)計(jì)具有指標(biāo)ωp=0.4π,Rp=0.25dB,ωs=0.2π,As=50dB的高通數(shù)字FIR濾波器。選擇合適的窗函數(shù)，確定并畫出濾波器的沖激響應(yīng)和頻率響應(yīng)。%Highpassfilterdesign-Hammingwindowwp=0.4*pi;ws=0.2*pi;tr_width=wp-ws;M=ceil(6.6*pi/tr_width)+1;n=[0:1:M-1];wc=(ws+wp)/2;hd=ideal_lp(wc,M);w_ham=(hamming(M))';h=hd.*w_ham;[db,mag,pha,grd,w]=freqz_m(h,[1]);delta_w=2*pi/1000;Rp=-(min(db(1:1:wp/delta_w+1)))%PassbandRippleAs=-round(max(db(ws/delta_w+1:1:501)))%MinStopbandattenuation%plotssubplot(1,1,1)subplot(2,2,1);stem(n,hd);title('IdealImpulseResponse')axis([0M-1-0.10.3]);xlabel('n');ylabel('hd(n)')subplot(2,2,2);stem(n,w_ham);title('HammingWindow')axis([0M-101.1]);xlabel('n');ylabel('w(n)')subplot(2,2,3);stem(n,h);title('ActualImpulseResponse')axis([0M-1-0.10.3]);xlabel('n');ylabel('h(n)')subplot(2,2,4);plot(w/pi,db);title('MagnitudeResponseindB');gridaxis([01-10010]);xlabel('frequencyinpiunits');ylabel('Decibels')set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1])set(gca,'YTickMode','manual','YTick',[-50,0])set(gca,'YTickLabelMode','manual','YTickLabels',['50';'0'])設(shè)計(jì)具有指標(biāo)ωs1=0.35π,Rp=1dB,ωs2=0.65π,Rp=1dB;ωp1=0.2π,As=60dB,ωp2=0.8π,As=60dB的數(shù)字帶阻FIR濾波器。選擇合適的窗函數(shù)，確定并畫出濾波器的沖激響應(yīng)和頻率響應(yīng)。%Bandstopfilterdesign-Blackmanwindowws1=0.35*pi;wp1=0.2*pi;wp2=0.8*pi;ws2=0.65*pi;As=60;tr_width=min((ws1-wp1),(wp2-ws2))M=ceil(11*pi/tr_width)+1%;M=68n=[0:1:M-1];wc1=(ws1+wp1)/2;wc2=(wp2+ws2)/2;hd=ideal_lp(wc2,M)-ideal_lp(wc1,M);w_bla=(blackman(M))';h=hd.*w_bla;[db,mag,pha,grd,w]=freqz_m(h,[1]);delta_w=2*pi/1000;Rp=-min(db(wp1/delta_w+1:1:wp2/delta_w))%Actua;PassbandRippleAs=-round(max(db(ws2/delta_w+1:1:501)))%MinStopbandAttenuation%plotssubplot(1,1,1);subplot(2,2,1);stem(n,hd);title('IdealImpulseResponse')axis([0M-1-0.40.5]);xlabel('n');ylabel('hd(n)')subplot(2,2,2);stem(n,w_bla);title('BlackmanWindow')axis([0M-101.1]);xlabel('n');ylabel('w(n)')subplot(2,2,3);stem(n,h);title('ActualImpulseResponse')axis([0M-1-0.40.5]);xlabel('n');ylabel('h(n)')subplot(2,2,4);plot(w/pi,db);%set(gca,'FontName','cmr12');title('MagnitudeResponseindB');grid;xlabel('frequencyinpiunits');ylabel('Decibels')axis([01-15010]);set(gca,'XTickMode','manual','XTick',[0,0.2,0.35,0.65,0.8,1])set(gca,'YTickMode','manual','YTick',[-60,0])set(gca,'YTickLabelMode','manual','YTickLabels',['60';'0']；第三題.利用頻率采樣法設(shè)計(jì)FIR濾波器2頻率采樣技術(shù)應(yīng)用(最優(yōu)法):參照實(shí)驗(yàn)內(nèi)容9的程序或自編程序來實(shí)現(xiàn)1)設(shè)計(jì)低通濾波器:wp=0.3pi,ws=0.4pi,Rp=0.25dB,As=50dBT1=0.5925,T2=0.1099M=60;alpha=(M-1)/2;l=0:M-1;wl=(2*pi/M)*l;Hrs=[ones(1,7),0.5925,0.11,zeros(1,43),0.11,0.5925,ones(1,6)];Hdr=[1,1,0,0];wdl=[0,0.3,0.4,1];k1=0:floor((M-1)/2);k2=floor((M-1)/2)+1:M-1;angH=[-alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2)];H=Hrs.*exp(j*angH);h=real(ifft(H,M));[db,mag,pha,grd,w]=freqz_m(h,1);[Hr,ww,a,L]=Hr_Type2(h);subplot(1,1,1)subplot(2,2,1);plot(wl(1:31)/pi,Hrs(1:31),'o',wdl,Hdr);axis([0,1,-0.1,1.1]);title('低通:M=60,T1=0.59,T2=0.109')xlabel('');ylabel('Hr(k)')set(gca,'XTickMode','manual','XTick',[0,0.3,0.4,1])set(gca,'YTickMode','manual','YTick',[0,0.109,0.59,1]);gridsubplot(2,2,2);stem(l,h);axis([-1,M,-0.1,0.3])title('脈沖響應(yīng)');ylabel('h(n)');text(M+1,-0.1,'n')subplot(2,2,3);plot(ww/pi,Hr,wl(1:31)/pi,Hrs(1:31),'o');axis([0,1,-0.1,1.1]);title('振幅響應(yīng)')xlabel('頻率（單位：pi）');ylabel('Hr(w)')set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1])set(gca,'YTickMode','manual','YTick',[0,0.109,0.59,1]);gridsubplot(2,2,4);plot(w/pi,db);axis([0,1,-100,10]);gridtitle('幅度響應(yīng)');xlabel('頻率（單位：pi）');ylabel('分貝數(shù)');set(gca,'XTickMode','Manual','XTick',[0;0.2;0.3;1]);set(gca,'YTickMode','Manual','YTick',[-63;0]);set(gca,'YTickLabelMode','manual','YTickLabels',['63';'0'])2)設(shè)計(jì)高通濾波器:wp=0.4pi,ws=0.3pi,Rp=0.25dB,As=50dBT1=0.5925,T2=0.1099M=33;alpha=(M-1)/2;l=0:M-1;wl=(2*pi/M)*l;T1=0.1095;T2=0.598;Hrs=[zeros(1,11),T1,T2,ones(1,8),T2,T1,zeros(1,10)];Hdr=[0,0,1,1];wdl=[0,0.3,0.4,1];k1=0:floor((M-1)/2);k2=floor((M-1)/2)+1:M-1;angH=[-alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2)];H=Hrs.*exp(j*angH);h=real(ifft(H,M));[db,mag,pha,grd,w]=freqz_m(h,1);[Hr,ww,a,L]=Hr_Type1(h);subplot(1,1,1)subplot(2,2,1);plot(wl(1:17)/pi,Hrs(1:17),'o',wdl,Hdr);axis([0,1,-0.1,1.1]);title('高通:M=33,T1=0.1095,T2=0.598')xlabel('');ylabel('Hr(k)')set(gca,'XTickMode','manual','XTick',[0;.3;.4;1])set(gca,'XTickLabelMode','manual','XTickLabels',['0';'.6';'.8';'1'])set(gca,'YTickMode','manual','YTick',[0,0.109,0.59,1]);gridsubplot(2,2,2);stem(l,h);axis([-1,M,-0.4,0.4])title('脈沖響應(yīng)');ylabel('h(n)');text(M+1,-0.4,'n')subplot(2,2,3);plot(ww/pi,Hr,wl(1:17)/pi,Hrs(1:17),'o');axis([0,1,-0.1,1.1]);title('振幅響應(yīng)')xlabel('頻率（單位：pi）');ylabel('Hr(w)')set(gca,'XTickMode','manual','XTick',[0;.3;.4;1])set(gca,'XTickLabelMode','manual','XTickLabels',['0';'.6';'.8';'1'])set(gca,'YTickMode','manual','YTick',[0,0.109,0.59,1]);gridsubplot(2,2,4);plot(w/pi,db);axis([0,1,-100,10]);gridtitle('幅度響應(yīng)');xlabel('頻率（單位：pi）');ylabel('分貝數(shù)');set(gca,'XTickMode','manual','XTick',[0;.3;.4;1])set(gca,'XTickLabelMode','manual','XTickLabels',['0';'.3';'.4';'1'])set(gca,'YTickMode','Manual','YTick',[-50;0]);set(gca,'YTickLabelMode','manual','YTickLabels',['50';'0'])3)設(shè)計(jì)帶通濾波器:ws1=0.3pi,wp1=0.35pi,wp2=0.40pi,ws2=0.45pi,Rp=1dB,As=60dBT2=0.59417456,T1=0.109021M=40;alpha=(M-1)/2;l=0:M-1;wl=(2*pi/M)*l;T1=0.109021;T2=0.59417456;Hrs=[zeros(1,5),T1,T2,ones(1,7),T2,T1,zeros(1,9),T1,T2,ones(1,7),T2,T1,zeros(1,4)];Hdr=[0,0,1,1,0,0];wdl=[0,0.3,0.35,0.4,0.45,1];k1=0:floor((M-1)/2);k2=floor((M-1)/2)+1:M-1;angH=[-alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2)];H=Hrs.*exp(j*angH);h=real(ifft(H,M));[db,mag,pha,grd,w]=freqz_m(h,1);[Hr,ww,a,L]=Hr_Type2(h);subplot(1,1,1)subplot(2,2,1);plot(wl(1:21)/pi,Hrs(1:21),'o',wd