%Samp7_14
clf;n=30;                 %濾波器階數
f=[0 0.9];a=[0 0.9];          %理想濾波器響應
b=firls(n,f,a,'d');          %采用 firls 設計濾波器
[h,w]=freqz(b);           %計算濾波器頻率響應
bb=remez(n,f,a,'d');       % 采用 remez 設計濾波器
[hh,w]=freqz(bb);           %計算濾波器脈沖響應
figure(1)
plot(w/pi,abs(h),'b',w/pi, abs(hh),'r--',f,a,'m:'); %繪幅頻響應
xlabel('歸一化頻率');ylabel('振幅');
legend('firls', 'remez','理想特性',4);  %繪制圖例，4表示圖例位置在右下角
grid on ;
f1=5;                     %輸入信號頻率
t=0:1/1000:1;               %時間序列
x=sin(2* pi*f1* t);         %輸入正弦信號
y=fftfilt(bb,x);          %濾波
figure(2)
subplot(2,1,1), plot(t,x),title ('輸入信號')% 繪輸入信號
subplot(2,1,2), plot(t,y),            %繪輸出信號
ylim([-0.01 0.01]),xlabel('時間/s'),title('輸出信號')