correl.m

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

function [a,S,R]=correl(x,P)

L=length(x);
x=[x(:);zeros(P,1)];
X=x;
for p=1:P, X=[X cshift(x, -p)]; end;
X=X((P+1):L,:);
R=X'*X;
a=inv(R)*[1;zeros(P,1)];
S=1/(a(1));
a=a/a(1);