写一个可以嵌入MATLAB Function的梯形冲击路面激励模型

function [F] = trapezoidal_impact_road_model(t, duration, amplitude, rise_time, fall_time)

% t: 时间向量
% duration: 冲击持续时间
% amplitude: 冲击幅值
% rise_time: 上升时间
% fall_time: 下降时间

% 计算梯形冲击路面激励的力信号

F = zeros(size(t));
for i = 1:length(t)
if t(i) < duration
if t(i) < rise_time
F(i) = (amplitude / rise_time) * t(i);
elseif t(i) < duration - fall_time
F(i) = amplitude;
else
F(i) = (amplitude / fall_time) * (duration - t(i));
end
end
end

end