model.zip

%习题7-7 wn=[0.4,0.6,0.7,0.8]; %设置多通带 N=38; b=fir1(N,wn,'dc-1'); %令第一个频带为通带 freqz(b) figure(2) stem(b,'.'); xlabel('n'),ylabel('h(n)'); f=0:0.01:1 m(1:21)=1;m(22:41)=0;m(42:61)=1;m(62:81)=0;m(82:101)=1; %设置多通带 plot(f,m,'k:');hold on N=30; b=fir2(N,f,m); %令第一个频带为通带 [h,f1]=freqz(b); f1=f1./pi; plot(f1,abs(h)); legend('理想滤波器','设计滤波器'); xlabel('归一化频率'),ylabel('幅值');