samp10_1.m

数字信号处理的Matlab实现的随书光盘

%Samp10_1
t=0:0.01:1.49;   %时间信号，采样间隔为0.01s
sig=sin(2*pi*45*t);   %原始信号
echo=0.5*[zeros(1,20),sig(1:length(t)-20)];
 %在信号前面补加20个零，并使其振幅衰减一半作为回声信号
sigecho=sig+echo;
c=cceps(sigecho);    %求原始信号与回声信号叠加的复时谱
subplot(2,1,1),plot(t,sigecho);   %绘制原始信号与回声的叠加信号
xlabel('时间/s');axis([0 1.5 -1 1]);grid on
title('时间域信号')
subplot(2,1,2),plot(t,c);   %绘制其复时谱
xlabel('时间/s');axis([0 1.5 -1 1]);grid on
title('时谱域信号')