corrcomp.m

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

function [ry,rz,ryz,rzy] = corrcomp(n,y,z,ip,iq) 
      
for k=0:ip 
	r = corr(n,k,y,y);
	ry(k+1)=r;
end
for k=0:iq 
	r = corr(n,k,z,z);
	rz(k+1)=r;
end
for k=0:ip-1 
	r = corr(n,k,y,z);
	ryz(k+1)=r;
end
for k=1:iq-1 
	r = corr(n,k,z,y);
	rzy(k+1)=r;
end