%Fir2函數：設計具有任意頻率響應的FIR濾波器
f=[0 0.2 0.2 0.4 0.4 0.6 0.6 1];
m=[1 1 0 0 1 1 0 0];
b=fir2(40,f,m);
[H,W]=freqz(b,1,256);
plot(f,m,'r:',W/pi,abs(H),'b-');
legend('期望幅頻響應','實際幅頻響應');