有没有人能帮我看看自相关函数的功率谱估计的程序啊?程序如下：load'signal1.mat'Fs=1000;xn=s63;subplot(211);plot(xn);window=boxcar(length(xn));nfft=512;cxn=xcorr(xn,'unbiased');CXk=fft(cxn,nfft);Pxx=abs(CXk);index=0:round

Question

有没有人能帮我看看自相关函数的功率谱估计的程序啊?
程序如下：
load 'signal1.mat'
Fs=1000;
xn=s63;
subplot(211);
plot(xn);
window=boxcar(length(xn));
nfft=512;
cxn=xcorr(xn,'unbiased') ;
CXk=fft(cxn,nfft) ;
Pxx=abs(CXk) ;
index=0:round(nfft/2-1) ;
k=index*Fs/nfft;
plot_Pxx =10*log10(Pxx(index+1));
subplot(212);
plot(k,plot_Pxx);
最主要是求绝对值后面啊?到底有什么用啊?.

Answer 1

1,boxcar()是矩形窗.
2,cxn是信号的自相关；自相关的傅里叶变换的幅度分布就是功率谱,因此CXk=fft(cxn,nfft)表示信号的快速傅里叶变换.再取幅度abs,表示功率谱.然后后面的
index=0:round(nfft/2-1) ;
k=index*Fs/nfft;
plot_Pxx =10*log10(Pxx(index+1));
subplot(212);
plot(k,plot_Pxx);
这一段是用来画图的.