%Samp7_9
clf
f=0:0.1:1;                         %歸一化頻率點數
m=[0 0 1 1 0 0 1 1 1 0 0];          %幅頻特性值
Order=10;                            % 濾波器的階數
b=fir2(Order,f,m,hamming(Order+1));            %設計濾波器
[h,w]=freqz(b,1,128);                  %計算濾波器的頻率響應
subplot(2,1,1)
plot(f,m,w/pi,abs(h),'r:')       %繪制理想幅頻響應和設計的濾波器幅頻響應
legend('理想特性', '實際設計')   %給出圖例
title('Order=10');xlabel('歸一化頻率');ylabel('振幅');
Order=100;
b=fir2(Order,f,m,hamming(Order+1));               %設計階數為 100 的濾波器
[h,w]=freqz(b,1,128);                 %計算濾波器的頻率響應
subplot(2,1,2),plot(f,m,w/pi,abs(h),'r:'); %繪制理想幅頻響應和設計的幅頻響應
ylim([0 1])
legend('理想特性', '實際設計')   %給出圖例
title('Order=100');xlabel('歸一化頻率');ylabel('振幅');
f1=12;f2=36;                % 輸入信號的兩種頻率成分
t=0:1/100:2;                   % 時間序列
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t);     %輸入信號
y=fftfilt(b,x);                  %對輸入信號進行濾波
figure(2)
subplot(2,1,1), plot(t,x),title('輸入信號')           %繪制輸入信號
subplot(2,1,2),plot(t,y)                        %繪制輸出信號
hold on;plot(Order/2/100*ones(1,2),ylim, 'r') %繪制延遲到的時刻
title('輸出信號'),xlabel('時間/s')