功率谱估计周期图法
clear all;
Fs=1000;
%第一种情况:N=256
N=256;Nfft=256;%数据长度和FFT所用的数据长度
n=0:N-1;t=n/Fs;%采用的时间序列
xn=load('C:\Users\Administrator\Desktop\shuju\0.5xielou\0.5xielou1-1.txt');
Pxx=10*log10(abs(fft(xn,Nfft).^2)/N);
%Fourier振幅谱平方的平均值,并装换为dB
f=(0:length(Pxx)-1)*Fs/length(Pxx);%给出频率序列
subplot(211);plot(f,Pxx);%绘制功率谱曲线
xlabel('频率/Hz');ylabel('功率谱/dB');
title('周期图N=256');
grid on;
%第二种情况:N=1024
Fs=1000;
N=1024;Nfft=1024;%数据长度和FFT所用的数据长度
n=0:N-1;t=n/Fs;%采用的时间序列
xn=load('C:\Users\Administrator\Desktop\shuju\0.5xielou\0.5xielou1-1.txt');
Pxx=10*log10(abs(fft(xn,Nfft).^2)/N);
%Fourier振幅谱平方的平均值,并装换为dB
f=(0:length(Pxx)-1)*Fs/length(Pxx);%给出频率序列
subplot(211);plot(f,Pxx);%绘制功率谱曲线
xlabel('频率/Hz');ylabel('功率谱/dB');
title('周期图N=1024');
grid on;