%Samp7_10
clf;n=20;                  %濾波器階數
f=[0 0.4 0.5 1];            %頻率向量
a=[1 1 0 0];                 %振幅向量
b=firls(n,f,a);              %采用 firls 設計濾波器
[h,w1]=freqz(b);                %計算其頻率響應
bb=remez(n,f,a);             %采用 remez 設計濾波器
[hh,w2]=freqz(bb);              %計算濾波器的頻率響應
figure(1)
plot(w1/pi,abs(h),w2/pi,abs(hh),'r:');         %繪制濾波器的幅頻響應
xlabel('歸一化頻率');ylabel('振幅');
legend('firls','remez');                       %給出圖例
grid on;
figure(2)
fs=50;t=0:1/fs:2;  %采樣頻率和時間序列
f1=5;f2=15;  %輸入信號頻率
x1=sin(2*pi*f1*t)+8.*cos(2*pi*f2*t);
subplot(2,1,1),plot(t,x1)
title('原始信號')
y1=filter(b,1,x1);
y2=filter(bb,1,x1);
subplot(2,1,2),plot(t,y1,t,y2,'r:')
legend('firls','remez');                       %給出圖例
title('輸出信號')
xlabel('時間/s')