功率谱估计周期图法
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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;