9.利用频率采样法设计FIR滤波器1%1)频率采样技术:低通,最优法,T1=0.5%wp=0.2pi,ws=
0.3pi,Rp=0.25dB,As=50dBM=40;alpha=(M-1)/2;l=0:M-1;wl
=(2pi/M)l;Hrs=[ones(1,5),0.5,zeros(1,29),0.5,ones(1,4)];Hd
r=[1,1,0,0];wdl=[0,0.25,0.25,1];k1=0:floor((M-1)/2);k2=
floor((M-1)/2)+1:M-1;angH=[-alpha(2pi)/Mk1,alpha(2pi)/M
(M-k2)];H=Hrs.exp(jangH);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'',wdl,Hdr);axis([0,
1,-0.1,1.1]);title(''频率样本:M=40,T1=0.5'')xlabel('''');ylabel(''Hr(k
)'')set(gca,''XTickMode'',''manual'',''XTick'',[0,0.2,0.3,1])set(gca,''
YTickMode'',''manual'',''YTick'',[0,0.5,1]);gridsubplot(2,2,2);stem
(l,h);axis([-1,M,-0.1,0.3])title(''脉冲响应'');ylabel(''h(n)'');text(M
+1,-0.1,''n'')subplot(2,2,3);plot(ww/pi,Hr,wl(1:21)/pi,Hrs(1:21),
''o'');axis([0,1,-0.1,1.1]);title(''振幅响应'')xlabel(''频率（单位：pi）'');yl
abel(''Hr(w)'')set(gca,''XTickMode'',''manual'',''XTick'',[0,0.2,0.3,1])
set(gca,''YTickMode'',''manual'',''YTick'',[0,0.5,1]);gridsubplot(2,
2,4);plot(w/pi,db);axis([0,1,-100,10]);gridtitle(''幅度响应'');xlab
el(''频率（单位：pi）'');ylabel(''分贝数'');set(gca,''XTickMode'',''Manual'',''XTi
ck'',[0;0.2;0.3;1]);set(gca,''YTickMode'',''Manual'',''YTick'',[-30;0])
;set(gca,''YTickLabelMode'',''manual'',''YTickLabels'',[''30'';''0''])%
2)频率采样技术:低通,最优法,T1&T2%wp=0.2pi,ws=0.3pi,Rp=0
.25dB,As=50dB%T1=0.5925,T2=0.1099M=60;alpha
=(M-1)/2;l=0:M-1;wl=(2pi/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.2,0.3,1];k1=0:floor((M-1)/2);k2=floor((M-1)/2)+1:M-1;ang
H=[-alpha(2pi)/Mk1,alpha(2pi)/M(M-k2)];H=Hrs.exp(ja
ngH);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:3
1)/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,''XTic
kMode'',''manual'',''XTick'',[0,0.2,0.3,1])set(gca,''YTickMode'',''manua
l'',''YTick'',[0,0.109,0.59,1]);gridsubplot(2,2,2);stem(l,h);axi
s([-1,M,-0.1,0.3])title(''脉冲响应'');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(''振幅响应'')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(''幅度响应'');xlabel(
''频率（单位：pi）'');ylabel(''分贝数'');set(gca,''XTickMode'',''Manual'',''XTick''
,[0;0.2;0.3;1]);set(gca,''YTickMode'',''Manual'',''YTick'',[-63;0]);s
et(gca,''YTickLabelMode'',''manual'',''YTickLabels'',[''63'';''0''])%3)频
率采样技术:带通,最优法T1&T2%ws1=0.2pi,wp1=0.35pi,wp2=0.
65pi,%ws2=0.8pi,Rp=1dB,As=60dB%T2=
0.59417456,T1=0.109021M=40;alpha=(M-1)/2;l=0:M-1;wl=
(2pi/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.2,0.35,0.65,0.8,1];k1=0:floor(
(M-1)/2);k2=floor((M-1)/2)+1:M-1;angH=[-alpha(2pi)/Mk1,
alpha(2pi)/M(M-k2)];H=Hrs.exp(jangH);h=real(ifft(H,M))
;[db,mag,pha,grd,w]=freqz_m(h,1);[Hr,ww,a,L]=Hr_Type2(h);s
ubplot(1,1,1)subplot(2,2,1);plot(wl(1:21)/pi,Hrs(1:21),''o'',wdl,H
dr);axis([0,1,-0.1,1.1]);title(''带通:M=40,T1=0.5941,T2=0.109'')
xlabel('''');ylabel(''Hr(k)'')set(gca,''XTickMode'',''manual'',''XTick''
,[0,0.2,0.35,0.65,0.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(''脉冲响应'');ylabel(''h(n)'');text(M+1,-0.4,''n'')subplot(2,2
,3);plot(ww/pi,Hr,wl(1:21)/pi,Hrs(1:21),''o'');axis([0,1,-0.1,1.1
]);title(''振幅响应'')xlabel(''频率（单位：pi）'');ylabel(''Hr(w)'')set(gca,''X
TickMode'',''manual'',''XTick'',[0,0.2,0.35,0.65,0.8,1])set(gca,''YTic
kMode'',''manual'',''YTick'',[0,0.109,0.59,1]);gridsubplot(2,2,4);pl
ot(w/pi,db);axis([0,1,-100,10]);gridtitle(''幅度响应'');xlabel(''频率（
单位：pi）'');ylabel(''Decibels'');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''])%4)频率采样技术:高通,最优法T1%ws=0.6pi,wp=0.8pi,Rp=
1dB,As=50dB%M=33,T1=0.1095;T2=0.598;M=33;
alpha=(M-1)/2;l=0:M-1;wl=(2pi/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.6,0.8,1];k1=0:floor((M-1)/2);k2=floo
r((M-1)/2)+1:M-1;angH=[-alpha(2pi)/Mk1,alpha(2pi)/M(M-k
2)];H=Hrs.exp(jangH);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)subpl
ot(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;.6;.8;1])set(gc
a,''XTickLabelMode'',''manual'',''XTickLabels'',[''0'';''.6'';''.8'';''1''])
set(gca,''YTickMode'',''manual'',''YTick'',[0,0.109,0.59,1]);gridsubp
lot(2,2,2);stem(l,h);axis([-1,M,-0.4,0.4])title(''脉冲响应'');ylabe
l(''h(n)'');text(M+1,-0.4,''n'')subplot(2,2,3);plot(ww/pi,Hr,wl(1:1
7)/pi,Hrs(1:17),''o'');axis([0,1,-0.1,1.1]);title(''振幅响应'')xlabel(
''频率（单位：pi）'');ylabel(''Hr(w)'')set(gca,''XTickMode'',''manual'',''XTick
'',[0;.6;.8;1])set(gca,''XTickLabelMode'',''manual'',''XTickLabels'',[''
0'';''.6'';''.8'';''1''])set(gca,''YTickMode'',''manual'',''YTick'',[0,0.10
9,0.59,1]);gridsubplot(2,2,4);plot(w/pi,db);axis([0,1,-100,10]
);gridtitle(''幅度响应'');xlabel(''频率（单位：pi）'');ylabel(''分贝数'');set(gc
a,''XTickMode'',''manual'',''XTick'',[0;.6;.8;1])set(gca,''XTickLabelMo
de'',''manual'',''XTickLabels'',[''0'';''.6'';''.8'';''1''])set(gca,''YTickM
ode'',''Manual'',''YTick'',[-50;0]);set(gca,''YTickLabelMode'',''manual''
,''YTickLabels'',[''50'';''0''])10.利用频率采样法设计FIR滤波器2频率采样技术应用(最优法):参照
实验内容9的程序或自编程序来实现1)设计低通滤波器:wp=0.3pi,ws=0.4pi,Rp=0.25
dB,As=50dBT1=0.5925,T2=0.10992)设计高通滤波器:wp=
0.4pi,ws=0.3pi,Rp=0.25dB,As=50dBT1=0.5925,T2=0.1099
3)设计带通滤波器:ws1=0.3pi,wp1=0.35pi,wp2=0.40pi,ws2
=0.45pi,Rp=1dB,As=60dBT2=0.59417456,T1=0.1090214)设计
带阻滤波器:ws1=0.35pi,wp1=0.30pi,wp2=0.45pi,ws2=0.4
pi,Rp=1dB,As=60dBT2=0.59417456,T1=0.109021四、实验报告要求
1.简述实验目的及实验原理；2.记录实验步骤3、6、8、10的（程序和要求）的图形(图形的标题及横纵轴要标示明确)。3.总结
FIR滤波器的设计步骤。实验七基于MATLAB的FIRDF设计一、实验目的1.练习并掌握窗函数法设计FIR数字滤
波器的方法；2.练习并掌握频率采样法设计FIR数字滤波器的方法；二、实验原理与方法1.线性相位FIR数字滤波器的条件和
特点满足第一类线性相位的条件是：h(n)是实序列且对(N-1)/2偶对称，即h(n)=h(N-n-1)
满足第二类线性相位的条件是：h(n)是实序列且对(N-1)/2奇对称，即h(n)=-h(N-n-1)窗函
数：1、rectwin；2、hamming；3、hann；4、blackman；5、kaiser；6、bartlet
t；7、triang数学函数：1、ceil；2、round；3、floor；4、real显示函数：1、
window；2、set滤波器函数：1、fir1；2、fir2其他函数：wvtool相关函数幅度特性:1
)h(n)=h(N-n-1),N=奇数2)h(n)=h(N-n-1),N=偶数3)h(n)=-h(N-n-1),
N=奇数4)h(n)=-h(N-n-1),N=偶数2.常用的窗函数1)三角形窗(BartlettWindow)2)
汉宁(Hanning)窗，又称升余弦窗3)汉明(Hamming)窗，又称改进的升余弦窗4)布莱克曼(Blankman)窗，
又称二阶升余弦5)凯泽(Kaiser)窗6)矩形(Rectangle)窗六种窗函数的基本参数-8010pi/N或10pi
/N-57凯塞窗a=7.865-7411pi/N或12pi/N-57布莱克曼窗-536.6pi/N或8pi/N-
41哈明窗-446.2pi/N或8pi/N-31汉宁窗-254.2pi/N或8pi/N-25三角形窗-211
.8pi/N或4pi/N-13矩形窗阻带最小衰减/dB过渡带宽或主瓣宽度旁瓣峰值幅度/dB窗函数3.用窗函数设计FI
R滤波器的步骤由允许的过渡带宽度(近似等于窗函数主瓣宽度)△ω及阻带衰减查表选择合适的窗函数并估计△ω=A/N确定滤波器的对称
中心延时值1)2)3)4)5)6)4.利用频率采样法设计FIR滤波器4.相关函数-系统函数:1)窗函数n=
50;w1=bartlett(n);w2=blackman(n);w3=hamming(n);w4=hann(n);w5
=rectwin(n);w6=triang(n);wvtool(w1,w2,w3);wvtool(w4,w5,w6);w
7=kaiser(n,beta);2)window-窗函数选择w=window(fhandle,n);%fhandle指定函
数句柄的n点的window窗口%函数句柄:@bartlett,@hamming,@hann,@blackman,%函数句
柄:@rectwin,@triang,@hann,@blackman,3)fir1-设计加窗有限冲激响应滤波器b=fir
1(n,wn);%0high’或ftype=‘stop’b=fir1(n,wn,window);%window指定具体的窗函数b=fir1(n,
wn,’ftype’,window);4)fir2-设计频率取样的有限冲激响应滤波器b=fir2(n,f,m);
%如f=[00.60.61],m=[1100]b=fir2(n,f,m,npt);%npt指定内插的点数
b=fir2(n,f,m,npt,window);%window指定具体的窗函数function[Hr,w,a,L]=
Hr_Type1(h);%计算第一种低通滤波器设计的振幅响应Hr(w)%-----------------------
------------------------------------%[Hr,w,a,L]=Hr_Type1(h)%
Hr=振幅响应%w=在[0pi]区间内计算Hr的500个频率点%a=第一种低通滤波器的系数%
L=Hr的阶次%h=第一种低通滤波器的频率响应%N=length(h);L=(N-1)/2;
a=[h(L+1)2h(L:-1:1)];%1乘(L+1)行向量rowvectorn=[0:1:L];
%(L+1)乘1列向量w=[0:1:500]''pi/500;Hr=cos(wn)a
'';1)计算第一种低通滤波器设计的振幅响应Hr(w)4.相关函数-扩展函数:h(n)=h(N-n-1),N=奇数2)计算第
二种低通滤波器设计的振幅响应Hr(w)h(n)=h(N-n-1),N=偶数function[Hr,w,b,L]=Hr_T
ype2(h);%计算第二种低通滤波器设计的振幅响应Hr(w)%----------------------------
-------------------------------%[Hr,w,b,L]=Hr_Type2(h)%Hr
=振幅响应%w=在[0pi]区间内计算Hr的频率点%b=第二种低通滤波器的系数%L=Hr的
阶次%h=第二种低通滤波器的频率响应%N=length(h);L=N/2;b=2[h(L:
-1:1)];n=[1:1:L];n=n-0.5;w=[0:1:500]''pi/500;Hr=cos
(wn)b'';3)计算第三种低通滤波器设计的振幅响应Hr(w)h(n)=-h(N-n-1),N=奇数function[
Hr,w,c,L]=Hr_Type3(h);%计算第三种低通滤波器设计的振幅响应Hr(w)%------------
-----------------------------------------------%[Hr,w,c,L]=Hr
_Type3(h)%Hr=振幅响应%w=在[0pi]区间内计算Hr的频率点%c=第三种低通滤波
器的系数%L=Hr的阶次%h=第三种低通滤波器的频率响应%N=length(h);L=(N-
1)/2;c=[2h(L+1:-1:1)];n=[0:1:L];w=[0:1:500]''pi/500;
Hr=sin(wn)c'';4)计算第四种低通滤波器设计的振幅响应Hr(w)h(n)=-h(N-n-1),N=偶数f
unction[Hr,w,d,L]=Hr_Type4(h);%计算第四种低通滤波器设计的振幅响应Hr(w)%---
------------------------------------------------%[Hr,w,d,L]=H
r_Type4(h)%Hr=振幅响应%w=在[0pi]区间内计算Hr的500个频率点%d=第四
种低通滤波器的系数%L=Hr的阶次%h=第四种低通滤波器的频率响应%N=length(h);L
=N/2;d=2[h(L:-1:1)];n=[1:1:L];n=n-0.5;w=[0:1:50
0]''pi/500;Hr=sin(wn)d'';三、实验内容及步骤1.fir1函数应用练习1:wn=[0.350
.65];N=34;%汉明窗b=fir1(N,wn,hamming(N+1));figure(1);freqz(b,
1,512);title(''使用汉明窗'')%海宁窗b=fir1(N,wn,hanning(N+1));figure(2)
;freqz(b,1,512);title(''使用海宁窗'')%矩形窗b=fir1(N,wn,boxcar(N+1));
figure(3);freqz(b,1,512);title(''使用矩形窗'')%三角窗b=fir1(N,wn,tr
iang(N+1));figure(4);freqz(b,1,512);title(''使用三角窗'')%布莱克窗b=f
ir1(N,wn,blackman(N+1));figure(5);freqz(b,1,512);title(''使用布莱
克窗'')%凯泽窗b=fir1(N,wn,kaiser(N+1,5.6));figure(6);freqz(b,1,512
);title(''使用凯泽窗,beta=6.5'')b=fir1(n,wn,window);wn=[0.480.68]
;N=34;%汉明窗b=fir1(N,wn,''stop'',hamming(N+1));figure(1)freqz(b,
1,512);title(''使用汉明窗'')%海宁窗b=fir1(N,wn,''stop'',hanning(N+1));fi
gure(2)freqz(b,1,512);title(''使用海宁窗'')%矩形窗b=fir1(N,wn,''stop'',
boxcar(N+1));figure(3)freqz(b,1,512),title(''使用矩形窗'')wn=[0.480.
68];N=34;%三角窗b=fir1(N,wn,''stop'',triang(N+1));figure(4)freq
z(b,1,512),title(''使用三角窗'')%布莱克窗b=fir1(N,wn,''stop'',blackman(N+
1));figure(5)freqz(b,1,512),title(''使用布莱克窗'')%凯泽窗b=fir1(N,wn,
''stop'',kaiser(N+1,5.6));figure(6)freqz(b,1,512);title(''使用凯泽窗,
beta=6.5'')2.fir1函数应用练习2:b=fir1(n,wn,’ftype’,window);3.fir1函数应用
练习1)利用b=fir1(n,wn,’ftype’,window);设计高通滤波器.2)利用b=fir1(n,wn,win
dow);设计带通滤波器.自己指定n和wn按下面要求设计滤波器,并总结规律.4.窗函数应用练习---注意不同窗函数的时频特性
n=50;w1=bartlett(n);w2=blackman(n);w3=hamming(n);w4=hann(n);
w5=rectwin(n);w6=triang(n);wvtool(w1,w2,w3);wvtool(w4,w5,w6);
w7=kaiser(n,beta);5.扩展函数应用练习1已知:h(n)={-4,1,-1,-2,5,6,5,-2,-
1,1,-4},确定第一种线性相位滤波器的振幅响应Hr(ω).h=[-4,1,-1,-2,5,6,5,-2,-1,1,-4]
;N=length(h);n=0:N-1;[Hr,w,a,L]=Hr_Type1(h);amax=max(
a)+1;amin=min(a)-1;subplot(2,2,1);stem(n,h);axis([-12L+1
aminamax])xlabel(''n'');ylabel(''h(n)'');title(''ImpulseRespons
e'')subplot(2,2,3);stem(0:L,a);axis([-12L+1aminamax])xlab
el(''n'');ylabel(''a(n)'');title(''a(n)coefficients'')subplot(2,2,
2);plot(w/pi,Hr);gridxlabel(''frequencyinpiunits'');ylabel(''Hr
'')title(''Type-1AmplitudeResponse'')subplot(2,2,4);zplane(h,1)
6.扩展函数应用练习21)已知:h(n)={-4,1,-1,-2,5,5,-2,-1,1,-4},参照实验内容5的程序或自
编程序确定并绘出第二种线性相位滤波器的单位冲激响应h(n)和振幅响应Hr(ω).2)已知:h(n)={-4,1,-1,-
2,5,0,-5,2,1,-1,4},参照实验内容5的程序或自编程序确定第三种线性相位滤波器的单位冲激响应h(n)和振幅响
应Hr(ω).3)已知:h(n)={-4,1,-1,-2,5,-5,2,1,-1,4},参照实验内容5的程序或自编程序
确定第四种线性相位滤波器的单位冲激响应h(n)和振幅响应Hr(ω).7.利用窗函数设计FIR滤波器11)设计具有指标ωp
=0.2π,Rp=0.25dB,ωs=0.3π,As=50dB的低通数字FIR滤波器。选择合适的窗函数，确定冲激响应，并画出
滤波器的频率响应。%Lowpassfilterdesign-Hammingwindowwp=0.2pi;w
s=0.3pi;tr_width=ws-wpM=ceil(6.6pi/tr_width)+1n=[0
:1:M-1];wc=(ws+wp)/2hd=ideal_lp(wc,M);w_ham=(hamming(M))
'';h=hd.w_ham;[db,mag,pha,grd,w]=freqz_m(h,[1]);delta_w
=2pi/1000;Rp=-(min(db(1:1:wp/delta_w+1)))%Passband
RippleAs=-round(max(db(ws/delta_w+1:1:501)))%MinStopbandat
tenuation%plotssubplot(1,1,1)subplot(2,2,1);stem(n,hd);titl
e(''IdealImpulseResponse'')axis([0M-1-0.10.3]);xlabel(''n'');
ylabel(''hd(n)'')subplot(2,2,2);stem(n,w_ham);title(''HammingWind
ow'')axis([0M-101.1]);xlabel(''n'');ylabel(''w(n)'')subplot(2,2
,3);stem(n,h);title(''ActualImpulseResponse'')axis([0M-1-0.1
0.3]);xlabel(''n'');ylabel(''h(n)'')subplot(2,2,4);plot(w/pi,db);
title(''MagnitudeResponseindB'');gridaxis([01-10010]);xlabe
l(''frequencyinpiunits'');ylabel(''Decibels'')set(gca,''XTickMode
'',''manual'',''XTick'',[0,0.2,0.3,1])set(gca,''YTickMode'',''manual'',''Y
Tick'',[-50,0])set(gca,''YTickLabelMode'',''manual'',''YTickLabels'',[''
50'';''0''])2)设计具有指标ωp1=0.35π,Rp=1dB,ωp2=0.65π,Rp=1dB;ωs1=0
.2π,As=60dB,ωs2=0.8π,As=60dB的数字带通FIR滤波器。选择合适的窗函数，确定并画出滤波器的冲激响应
和频率响应。%Bandpassfilterdesign-Blackmanwindowws1=0.2pi;
wp1=0.35pi;wp2=0.65pi;ws2=0.8pi;As=60;tr_w
idth=min((wp1-ws1),(ws2-wp2))M=ceil(11pi/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=2pi/1000;Rp
=-min(db(wp1/delta_w+1:1:wp2/delta_w))%Actua;PassbandRipple
As=-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''])1)设计具有指标ωp=0.4π,Rp=0.25dB,ωs=0.2π,As=50dB的高通数字FIR滤波器。选择合适的窗函数，确定并画出滤波器的冲激响应和频率响应。8.利用窗函数设计FIR滤波器2:参照实验内容7的程序或自编程序来实现2)设计具有指标ωs1=0.35π,Rp=1dB,ωs2=0.65π,Rp=1dB;ωp1=0.2π,As=60dB,ωp2=0.8π,As=60dB的数字带阻FIR滤波器。选择合适的窗函数，确定并画出滤波器的冲激响应和频率响应。 |
|
