zhijiefa.m

来自「对不同的功率谱密度分析方法的研究」· M 代码 · 共 19 行

Fs=200;
NFFT=6000;%指定FFT运算所采用的点数
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=fft(x,NFFT);
pxx=abs(X).^2/length(n);
t=0:round(NFFT/2-1);
k=t*Fs/NFFT;
p=10*log(pxx(t+1));
plot(k,p)
xlabel('频率(Hz)');
ylabel('相对功率谱密度(dB/Hz)');