jindianpuguji.m

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

% Fs=200;
% NFFT=280;%指定FFT运算所采用的点数
%基于经典谱估计的系统辨识
Fs=2000;
NFFT=256;
n=0:1/Fs:1;
% load 100a.txt
% x=X100a;
% 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)');
x=randn(size(n));
b=ones(1,5)/5;
y=filter(b,1,x);
[h,f]=tfe(x,y,NFFT,Fs,256,128,'none');
h0=freqz(b,1,f,Fs);
subplot(2,1,1)
plot(f,abs(h0));
title('理想传递函数的幅度');
xlabel('频率(Hz）');
subplot(2,1,2)
plot(f,abs(h));
title('估计的传递函数幅度');
xlabel('频率(Hz）');