pmpy.m

来自「这是一个用于语音信号处理的工具箱」· M 代码 · 共 16 行

function [Z,IZ]=pmpy(X,IX,Y,IY)
%
if IX*IY <= 0
	IZ=0;
	Z=0;
	return;
end
IZ=IX+IY-1;
Z=zeros(1,IZ);
for I=1:IX
	for J=1:IY
		K=I+J-1;
		Z(K)=X(I)*Y(J)+Z(K);
	end
end