sigadd.m

来自「西交大刘树棠写的《数字信号处理》那本书的所有源代码」· M 代码 · 共 12 行

function [y,n]=sigadd(x1,n1,x2,n2)
%y(n)=x1(n)+x2(n)
%------------------------
%y(n)=sigadd(x1,n1,x2,n2)
%
n=min(min(n1),min(n2)):max(max(n1),max(n2));
y1=zeros(1,length(n));y2=y1;
y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;
y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;
y=y1+y2;