%Samp7_13
clf;n=21;              %濾波器的階數(shù)為奇數(shù)
f=[0 0.1 0.2 1];a=[0 0 1 1];          %理想濾波器的幅頻特性
b=firls(n,f,a,'h');              %采用 firls 設(shè)計奇對稱系數(shù)濾波器
[h,w]=freqz(b,1,512);            %計算濾波器的脈沖響應(yīng)
bb=remez(n, f,a,'h');           %采用 remez 設(shè)計奇對稱系數(shù)濾波器
[hh,w]=freqz(b,1,512);           %計算濾波器的脈沖響應(yīng)
figure(1)
plot (w/pi,abs(h),'b.',w/pi, abs(hh),'r--',f,a,'m');     %繪濾波器幅頻響應(yīng)
xlabel('歸一化頻率');ylabel('振幅');
legend('firls','remez','理想特性');                %給出圖例
t=0:1/50:3;x=sin(2*pi*5*t);    %時間序列和輸入信號  
figure(2)
subplot(2,1,1),plot(t(1:100),x(1:100))   %繪制輸入信號的前100個樣本
title('輸入信號')
y1=filter(b,1,x);   %運用firls設(shè)計的濾波器進行濾波
y2=filter(bb,1,x);  %運用remez設(shè)計的濾波器進行濾波
subplot(2,1,2),
plot(t([1:100]+20/2),y1([1:100]+20/2),t([1:100]+20/2),y2([1:100]+20/2),'r:')
%繪制與輸入信號對應(yīng)的輸出信號，考慮了延遲效應(yīng)
legend('firls', 'remez')
xlabel('時間/s')
title('輸出信號')