MATLAB中的fft问题

Question

我的导师告诉我功率谱密度函数应该按照下面这个图中的公式算

所以我根据公式把代码改成如下所示：figure;N =length(a);Fs = 500; U = fft(a);P2 = abs(U.^2/2piNFs);P = P2(1:N/2+1);P(2:end-1) = 2P(2:end-1);f = Fs*(0:(N/2))/N;semilogx(f, P);
但是画出来的图纵坐标PSD非常大，如下图已经圈出来了（达到了10的17次方，也太大了）：

但是实际上纵坐标PSD应该是如下图这样的几千的范围内：

那么问题是，在导师给的公式没错的前提下，我应该怎么改代码或者设置图像能使我的图的纵坐标PSD恢复到像最后一个图这样几千几千的范围？请根据导师给的公式帮我改一下代码解决我的问题。

Answer 1

abs(U.^2/2piNFs)你这里是不是写成了

abs((U.^2/2)*pi*N*Fs)

的意思了呀？
应该是

abs(U.^2/(2*pi*N*Fs))