MATLAB轨道剖面PSD绘图与文献不符问题排查请求

我正在MATLAB/SIMULINK中进行列车建模，参考文献中轨道剖面采用功率谱密度（PSD）表示，给出了对应方程，核心参数为：

• Av = 0.0339（轨道状态标识量，单位：cm²·rad/m）
• ωc = 0.8245（轨道不平顺特征角频率，单位：rad/m）

文献附有对应的PSD参考绘图，但我编写的MATLAB代码生成的绘图与文献完全不符，请求协助排查错误：

clc
clear all
%% PSD of track profile
f = 0:0.5:100;            %Hz
% Av=0.0339;              %Track condition identifier (cm2rad/m)
Av=0.0339*(10e-4)/(2*pi); %m2Hz/m
%wc=0.8245;               %Constant pulsation for rail irregularity (rad/m)
fc=0.8245/(2*pi);         %(Hz/m)   

S=(0.25*Av.*fc.^2)./((f.^2+fc.^2).*f.^2);
semilogx(f, S)
xlabel('frequency [Hz]')
ylabel('power spectral density of track profile [m2/Hz]')
grid on

当前生成的绘图在半对数坐标系下，低频段出现异常数值尖峰，整体曲线形态与文献给出的PSD图完全不一致。

内容的提问来源于stack exchange，提问作者Hakan