lxblijia
2009-12-18, 16:27
m为采样的一段时程，样本个数20480，采样频率312.5.我分别采用psd和ar(取20阶)进行谱分析，为什么ar模型方法得出的谱值会有约只有psd求出的1/200.我初学，请大侠指教一二，不胜感激！！！
[pxx,f]=psd(m,1024,312.5);
semilogy(f,pxx)
for i=1:20 %取ar模型阶数小于等于20
[a,e]=armcov(m,i);
fpe(i)=(20480+i)*e/(20480-i); %用FPE准则判断ar模型的阶数
end
[b,n]=min(fpe);
hold on
[pxx1,f1]=pmcov(m,n,1024,312.5);
plot(f1,pxx1,'r')
[pxx,f]=psd(m,1024,312.5);
semilogy(f,pxx)
for i=1:20 %取ar模型阶数小于等于20
[a,e]=armcov(m,i);
fpe(i)=(20480+i)*e/(20480-i); %用FPE准则判断ar模型的阶数
end
[b,n]=min(fpe);
hold on
[pxx1,f1]=pmcov(m,n,1024,312.5);
plot(f1,pxx1,'r')