sigadd.m

PWS_DSPC 汉化的数字信号处理工具箱1.0版中的Matlab 函数

function [y,n] = sigadd(x1,n1,x2,n2)
% 实现 y(n) = x1(n)+x2(n)
% -----------------------------
% [y,n] = sigadd(x1,n1,x2,n2)
%   y = 在包含n1 和 n2 的n点上求序列和, 
%  x1 = 在 n1上的第一序列
%  x2 = 在 n2上的第二序列(n2可与 n1不等)
%
n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的长度
y1 = zeros(1,length(n)); y2 = y1;              % 初始化
y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;     % 具有y的长度的 x1 
y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;     % 具有y的长度的x2 
y = y1+y2;                                     % 序列相加.