glpmd.rar_petyw5_power spectral_psd_功率谱 曲线_功率谱曲线

Fs=200; NFFT=1024; n=0:1/Fs:1; x=load('优化前.txt'); 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');