corr.m

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

function [r] = corr(n,lag,x,y) 
      
sum=0.;
for j=1:n-lag 
	sum=sum+x(j)*y(j+lag);
end
r=sum/n;