example2_3.m

经典《信号与系统》教程的matlab例程,对深入理解信号与系统相关概念有很大帮助

% 程序名：example2_3.m
clear;close all
n0=-25;
n1=25;
t=n0:n1;
x=0.5.^(t).*u(t);
h=u(t);
y=conv(x,h);%卷积计算
%将卷积的结果取得和两个输入的信号一样长，以便在相同的坐标中绘制其图形并方便比较。
l=length(y);
m=2*n0:n0-1;
n=n1+1:2*n1;
m1=zeros(1,abs(n0));
p1=zeros(1,abs(n1));
lm=length(m1);
ln=length(p1);
y=y([lm+1:l-ln]);
c=u(-t);
%计算移位长度
hl=u(-t-15);
c=u(-t-15);
l=length(hl);
%以下程序绘制激励信号和系统冲激响应的波形
clf;
subplot(4,2,1);stem(t,x,'.');hold on;axis([-10,10,0,1.5])
subplot(4,2,2);
stem(t,h,'r.');hold on;axis([-10,10,0,1.5])
%下段程序实现移位动画
b=31;
for i=1:b
    for j=l-1:-1:1
        hl(j+1)=hl(j);
        c(j+1)=c(j);
    end
    c(1)=1;
    hl(1)=1;    
    subplot(4,1,3);
    stem(t,hl,'r.'); axis([-15,15,0,1.5])
    subplot(4,1,2);stem(t,x,'.');axis([-15,15,0,1.5])
    hold off;
    subplot(4,1,4);
    stem(t,c.*y,'r.');axis([-15,15,0,2.5]);    
    pause;
end