cadzow.m

来自「基于累计量的奇异值－总体最小二乘法求AR参数 用奇异值－总体最小二乘法求AR参数」· M 代码 · 共 21 行

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);