seqmult.m

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

function [y,n] = seqmult(x1,n1,x2,n2)
%《数字信号处理教程——MATLAB释义与实现》子程序
% 序列相乘子程序
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% [y,ny] = seqmult(x1,n1,x2,n2)
% -----------------------------
% 实现 y(n) = x1(n)*x2(n)
%  y = 在n区间上的乘积序列,ny 包含 n1 和 n2
%  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,y2初始置零
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;                                   % 序列相乘