%Samp7_12
clf;n=30;                      %濾波器的階數
f=[0 0.4 0.5 1];
a=[1 1 0 0];
weit=[1 10];
b=firls(n,f,a,weit);          %采用 firls 設計濾波器
[h,w1]=freqz(b);                %計算濾波器的頻率響應
bb=remez(n,f,a,weit);                  %采用 remez 設計濾波器
[hh,w2]=freqz(bb);                   % 計算濾波器的頻率響應
figure(1)
plot(w1/pi,abs(h),'b.',w2/pi, abs(hh),'r:',f,a,'m');%繪濾波器幅頻響應
xlabel('歸一化頻率');ylabel('振幅');
legend('firls','remez','理想特性');               %給出圖例