jianjiefa.m

对不同的功率谱密度分析方法的研究

Fs=200;
NFFT=6000;%指定FFT运算所采用的点数
% Fs=2000;
% NFFT=1024;
n=0:1/Fs:1;
% load 001.txt
% x=X0001;
load 001II.txt
x=X001II;
% x=sin(2*pi*100*n)+4*sin(2*pi*500*n)+randn(size(n));
Cx=xcorr(x,'unbiased');%计算序列的自相关函数
Cxk=fft(Cx,NFFT);%计算FFT
pxx=abs(Cxk);%求解PSD
t=0:round(NFFT/2-1);
k=t*Fs/NFFT;
p=10*log(pxx(t+1));
plot(k,p)
xlabel('频率(Hz)');
ylabel('相对功率谱密度(dB/Hz)');