sigmult.m

数字信号处理的Matlab实现的随书光盘

function [y,n] = sigmult(x1,n1,x2,n2)
% 执行信号相乘：y(n) = x1(n)*x2(n)
% 调用方式：[y,n] = sigmult(x1,n1,x2,n2)
%  y：输出的信号的乘积序列，序列号为n, 包括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(find((n>=min(n1))&(n<=max(n1))==1))=x1;    %将x1赋给y对应位置的元素
y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;    %将x2赋给y对应位置的元素
y = y1 .* y2;                               %信号相乘