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;                               %信號相乘