ex0812.m

数字信号处理（机械工业出版社）的源码

% Chapter 08: Example 8.12 Design a filter
f=0:0.002:1;
m(1:201)=1;
m(202:301)=0;
m(302:351)=0.5;
m(352:401)=0;
m(402:501)=1;
b = fir2(30,f,m);
freqz(b,1,128); %frequency response.
figure
[h,w] = freqz(b,1,128);
plot(f,m,w/pi,abs(h))
legend('Ideal','fir2 Designed')