ms5_5.m

《数字信号处理实验(MATLAB版)》这

clf;
nf1=0:20;                    %f1的时间向量
f1=0.8.^nf1;
lf1=length(f1);
nf2=0:10;                    %f2的时间向量
lf2=length(nf2);             %取f2时间向量的长度
f2=ones(1,lf2);
lmax=max(lf2,lf1);
if lf2>lf1 nf2=0;nf1=lf2-lf1; 
    elseif lf2<lf1 nf1=0;nf2=lf1-lf2;
    else nf2=0;lf1=0;
end
lt=lmax;
u=[zeros(1,lt),f2,zeros(1,nf2),zeros(1,lt)];
t1=(-lt+1:2*lt);
f1=[zeros(1,2*lt),f1,zeros(1,nf1)];
hf1=fliplr(f1);
N=length(hf1); 
y=zeros(1,3*lt);
for k=0:2*lt
    p=[zeros(1,k),hf1(1:N-k)];
    y1=u.*p
    yk=sum(y1);
    y(k+lt+1)=yk;
    subplot(4,1,1);stem(t1,u);
    subplot(4,1,2);stem(t1,p);  
    subplot(4,1,3);stem(t1,y1);  
    subplot(4,1,4);stem(k,yk);
    axis([-20,50,0,5]);hold on
    pause(1);
end