用matlab生成一个七段的S型加减速曲线

以下是一个简单的七段的S型加减速曲线的matlab代码：

% 定义插补参数
t_total = 10;  % 插补总时间
v_max = 10;  % 最大速度
a_max = 5;  % 最大加速度
j_max = 2;  % 最大加加速度
dt = 0.01;  % 时间间隔

% 初始化变量
t = 0:dt:t_total;  % 时间序列
s = zeros(size(t));  % 位置序列
v = zeros(size(t));  % 速度序列
a = zeros(size(t));  % 加速度序列
j = zeros(size(t));  % 加加速度序列

% 计算各段时间
t1 = v_max / a_max;
s1 = 0.5 * a_max * t1^2;
t5 = t_total - 2 * t1 - 2 * s1 / v_max;
s5 = v_max * t5;
t7 = t1;
s7 = s1;

% 计算各段加速度、速度、位置、加加速度
for i = 1:length(t)
    if t(i) <= t1
        a(i) = a_max;
        j(i) = j_max;
        v(i) = a_max * t(i);
        s(i) = 0.5 * a_max * t(i)^2;
    elseif t(i) <= t1 + t5/2
        a(i) = 0;
        j(i) = 0;
        v(i) = v_max;
        s(i) = s1 + v_max * (t(i)-t1);
    elseif t(i) <= t1 + t5
        a(i) = -a_max;
        j(i) = -j_max;
        v(i) = v_max - a_max*(t(i)-t1-t5/2);
        s(i) = s1 + v_max*t5/2 + v_max*(t(i)-t1-t5/2) - 0.5*a_max*(t(i)-t1-t5/2)^2;
    elseif t(i) <= t_total - t7
        a(i) = 0;
        j(i) = 0;
        v(i) = v_max;
        s(i) = s1 + s5 + v_max*(t(i)-t1-t5);
    elseif t(i) <= t_total - t1
        a(i) = -a_max;
        j(i) = -j_max;
        v(i) = v_max - a_max*(t_total-t(i)-t1-t5);
        s(i) = s1 + s5 + v_max*(t_total-t1-t5-t(i)) - 0.5*a_max*(t_total-t(i)-t1-t5)^2;
    else
        a(i) = -a_max;
        j(i) = -j_max;
        v(i) = v_max - a_max*t7;
        s(i) = s1 + s5 + s7 + v_max*t7 - 0.5*a_max*t7^2 + v_max*(t(i)-t_total+t1+t5+t7);
    end
end

% 绘制速度-时间图和位置-时间图
subplot(2,1,1);
plot(t, v);
xlabel('时间 (s)');
ylabel('速度 (m/s)');
title('速度-时间图');

subplot(2,1,2);
plot(t, s);
xlabel('时间 (s)');
ylabel('位置 (m)');
title('位置-时间图');

你可以根据需要修改插补参数和绘图方式。