理论如下：
式（13-7）中等式左边的路面位移随时间的导数可用差分法近似，即
整理可得
代码如下：
C级路面

clc
clear
close all

%% 参数定义
dt=0.01;%仿真间隔时间
t_end=10;%仿真总时长
t=0:dt:t_end;
n00=0.011;%下截止频率
u=60;%车速,km/h
u=u/3.6;%车速转化为m/s
f0=2* pi *