clear all;
close all;
dt=0.01;
t=0:dt:10;
xt=0.1*cos(0.15*pi*t)+0.5*cos(4*pi*t);
[f,xf]=T2F(t,xt);
%抽样信号,抽样频率为4Hz
fs=4;
sdt=1/fs;% 抽样间隔
t1=0:sdt:10;
st=0.1*cos(0.15*pi*t1)+0.5*cos(4*pi*t1);
[f1,sf]=T2F(t1,st);
%恢复原始信号
t2=-50:dt:50;
gt=sinc(fs*t2);%以dt为间隔,理想低通滤波器的冲激响应
stt=sigexpand(st,sdt/dt);%将以sdt为间隔的st扩展为以dt为间隔,与gt一致,注意到10时,也扩展这么多,为10+ sdt-dt
xt_t=conv(stt,gt);%抽样信号通过理想低通滤波器,恢复出原始信号
figure(1);
subplot(311);
plot(t,xt);
title('原始信号');
subplot(312);
plot(t1,st);
title('抽样信号');
subplot(313);
t3=-50:dt:60+sdt-dt;%卷积的长度(左坐标之和到右坐标之和),见stt的解释
plot(t3,xt_t);
title('抽样信号恢复');
axis([0 10 -1 1]);