xiefangchaguji.asv

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

% Fs=200;
% NFFT=280;%指定FFT运算所采用的点数
Fs=200;
% NFFT=1024;
% n=0:1/Fs:1;

% load 0001.txt
% x=X0001;
% load 100a.txt
% x=X100a;
% load 008.txt
% x=X008;
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)');
% h=fir1(20,0.3);
% r=randn(1024,1);
% x=filter(h,1,r);
[P1,f]=pcov(x,20,[],Fs);
[P2,f]=pmcov(x,20,[],Fs);
Pxx1=10*log10(P1);
Pxx2=10*log10(P2);
plot(f,Pxx1,':',f,Pxx2,'-.')
title('(a)MIT数据');
title('(a)MIT数据');
ylabel('相对幅度(dB)');
xlabel('功率谱估计)');
legend('协方差法','改进的协方差法');