seqadd.m

来自「matlab数字信号处理代码」· M 代码 · 共 18 行

function [y,n] = seqadd(x1,n1,x2,n2)
%《数字信号处理教程——MATLAB释义与实现》
% 序列相加子程序
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% [y,n] = seqadd(x1,n1,x2,n2)
% -----------------------------
% 实现y(n) = x1(n)+x2(n)
%  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))))=x1;     % 具有y的长度的x1 
y2(find((n>=min(n2))&(n<=max(n2))))=x2;     % 具有y的长度的x2 
y = y1+y2;                                     % 序列相加.