samp7_7.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 18 行

%Samp7_7
clf;N=1000;Fs=1000;    %数据总数和采样频率
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阶的低通滤波器，归一化截止频率据6-19式
yfft=fftfilt(b,x,256);           %对数据进行滤波
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')