%Samp7_7
clf;N=1000;Fs=1000;    %數(shù)據(jù)總數(shù)和采樣頻率
fc=200;
n=[0:N-1];t=n/Fs;            %時間序列
f1=50;f2=250;
x=sin(2*pi*f1*t)+sin(2*pi*f2*t);    %輸入信號
b=fir1(40,fc*2/Fs);     %設計40階的低通濾波器，歸一化截止頻率據(jù)6-19式
yfft=fftfilt(b,x,256);           %對數(shù)據(jù)進行濾波
n1=81:241;t1=t(n1);             %選擇采樣點間隔
x1=x(n1);           %與采樣點對應的輸入信號
subplot(2,1,1);plot(t1,x1); grid on;       %繪制輸入信號
title('輸入信號');
n2=n1-40/2;t2=t(n2);   %輸出信號，扣除了相位延遲N/2
y2=yfft(n2);
subplot(2,1,2);plot(t2,y2);          %繪制輸出信號
title('輸出信號');
grid on; xlabel('時間/s')