sigadd.m

数字信号处理—基于计算机的方法教材上的所有例题源程序

function [y,n]=sigadd(x1,n1,x2,n2)
n=min(n1(1), n2(1)): max(n1(end), n2(end));
y1=zeros(1,length(n)); y2=y1;

find((n>=n1(1)) & (n<=n1(end)))
find((n>=n2(1)) & (n<=n2(end)))

y1(find((n>=n1(1)) & (n<=n1(end))==1))=x1;
y2(find((n>=n2(1)) & (n<=n2(end))==1))=x2;
y=y1 + y2;