cadzow.m

基于累计量的奇异值－总体最小二乘法求AR参数 用奇异值－总体最小二乘法求AR参数 一般最小二乘法求AR参数 根据AR参数和自相关函数以及AR阶数用Cadzow谱估计子求出频谱密度

function Cadzow(a,Rx,p,N)
%N=128;
Ro=Rx;
Ro(N)=0.5*Ro(N);
[p,m]=size(a);
p=p-1;
n=zeros(1,p+1);
for k=1:p+1
    n(k)=a(1:p+1)'*Ro(k+N:-1:k-p+N)';
end
wfreq=[0.3:0.002:0.5]*pi;
%wfreq=[1:N]/N*pi;
z=exp(sqrt(-1)*wfreq);
Pw=polyval(n,z)./polyval(a,z)+polyval(n,z.^-1)./polyval(a,z.^-1);
Pw=abs(Pw);
Hpsd = dspdata.psd(Pw, wfreq);
plot(Hpsd);