covariance.m

来自「wav声音文件的通过多脉冲激励编码器」· M 代码 · 共 19 行

function lina=covariance(p,N2,y)
M=N2/2-1;
n=1+M;
for i=1:2*M+1
    for j=1:p
        if(n-M+i-1-j<=0)
            R(i,j)=0;
        else
        R(i,j)=y(n-M+i-1-j);
        end
    end
end
for i=1:2*M+1
    r(i)=y(n-M+i-1);
end
lina=zeros(p,1);%建立一个p元列向量
     lina=pinv(R)*r';