冲击道路建模

以正弦函数为例建立减速带模型，仅供参考

%冲击道路建模，241111
clc
clear
close all
warning off
%% 原始数据
dt=0.01;%仿真间隔时间
t_end=10;%仿真总时长
t=0:dt:t_end;
u=36;%车速,km/h
u=u/3.6;%车速转化为m/s
%% Hump道路输入建模
kuandu=0.4;%减速带宽度，米
gaodu=0.10;%减速带高度，米%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%振幅

q=zeros(1,length(t));
q(102)=sin(2*pi/kuandu/2*0.1)*gaodu;
q(103)=sin(2*pi/kuandu/2*0.2)*gaodu;
q(104)=sin(2*pi/kuandu/2*0.3)*gaodu;

dq=zeros(1,length(t));
dq(102)=cos(2*pi/kuandu/2*0.1)*gaodu*2*pi/kuandu/2;
dq(103)=cos(2*pi/kuandu/2*0.2)*gaodu*2*pi/kuandu/2;
dq(104)=cos(2*pi/kuandu/2*0.3)*gaodu*2*pi/kuandu/2;

figure(1)
length_road=(0:dt:10)*u;
plot(length_road,q(1:length(length_road)))
xlabel('行驶位移/m')
ylabel('路面位移/m')

Sim_q=[t',q'];
Sim_dq=[t',dq'];