【任何电机可使用的七段式s曲线-----包括matlab代码和simulink仿真】

最新推荐文章于 2024-06-10 08:46:56 发布
最新推荐文章于 2024-06-10 08:46:56 发布
阅读量1.8k 收藏 40
点赞数 19
文章标签： matlab 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44312889/article/details/136301388
版权
本文详细介绍了永磁同步电机的七段式S曲线加减速策略，包括理论分析和MATLAB中的实现过程，以及如何将其应用到Simulink中进行双闭环FOC控制的仿真示例。通过MATLAB代码计算出各阶段的加速度、速度和位置，并展示了在Simulink中模拟的整个过程和关键变量的变化情况。
摘要由CSDN通过智能技术生成

永磁同步电机七段式s曲线

一、前言

S形：加加速(T1)→>匀加速(T2)→减加速(T3)→>匀速(T4 )>加减速(T5)>匀减速(T6)-→>减减速(T7)

二、理论分析

1、这个是7阶段的加速度

在这里插入图片描述

加速度是𝑎max，J是加加速度，𝜏𝑘(𝑘 = 1，2，… ，7)是每个区间时间坐标，𝑡𝑘(k = 0，1，2… ，7)是每区间转折点时刻，𝑇𝑘(𝑘 = 1，2，…. ，7)是每个区间的连续工作时间；

2、这个是7阶段的速度

在这里插入图片描述

3、这个是7阶段的位置

在这里插入图片描述

三、7段式s曲线加减速关系图

在这里插入图片描述

四、matlab代码

% 定义参数

q0 = 0;
q1 = 200;
vmax =15;
amax = 10;
v0 = 0;
v1 = 0;
jmax = 10;
count = 0;

% 计算Ta、Td、Tj1、Tj2、Tv、alima、alimd等参数，他们就是时间点
if (vmax - v0) * jmax < amax^2
    if v0 > vmax
        Tj1 = 0;
        Ta = 0;
        alima = 0;
    else
        Tj1 = ((vmax - v0) / jmax)^0.5;
        Ta = 2 * Tj1;
        alima = Tj1 * jmax;
    end
else
    Tj1 = amax / jmax;
    Ta = Tj1 + (vmax - v0) / amax;
    alima = amax;
end

if (vmax - v1) * jmax < amax^2
    Tj2 = ((vmax - v1) / jmax)^0.5;
    Td = 2 * Tj2;
    alimd = Tj2 * jmax;
else
    Tj2 = amax / jmax;
    Td = Tj2 + (vmax - v1) / amax;
    alimd = amax;
end

Tv = (q1 - q0) / vmax - Ta / 2 * (1 + v0 / vmax) - Td / 2 * (1 + v1 / vmax);

% 动态调整参数，确保轨迹满足条件
if Tv <= 0
    Tv = 0;
    amax_org = amax;
    
    while Ta < 2 * Tj1 || Td < 2 * Tj2
        count = count + 1;
        amax = amax - amax_org * 0.1;
        alima = amax;
        alimd = amax;
        
        delta = (amax^4) / (jmax^2) + 2 * (v0^2 + v1^2) + amax * (4 * (q1 - q0) - 2 * amax / jmax * (v0 + v1));
        Tj1 = amax / jmax;
        Ta = (amax^2 / jmax - 2 * v0 + delta^0.5) / (2 * amax);
        Tj2 = amax / jmax;
        Td = (amax^2 / jmax - 2 * v1 + delta^0.5) / (2 * amax);
    end
end

% 计算轨迹
p = [];
vc = [];
ac = [];
jc = [];

if Tv > 0
    vlim = vmax;
    T = Tv + Ta + Td;
else
    Tv = 0;
    vlim = v0 + (Ta - Tj1) * alima;
    T = Tv + Ta + Td;
end

for t = 0:0.1:T
    if t >= 0 && t < Tj1
        % 段1：加加速度段
        q = q0 + v0 * t + jmax * t^3 / 6;
        p = [p q];
        v = v0 + jmax * t^2 / 2;
        vc = [vc v];
        a = jmax * t;
        ac = [ac a];
        jc = [jc jmax];
    elseif t >= Tj1 && t < (Ta - Tj1)
        % 段2：匀加速度段
        q = q0 + v0 * t + alima / 6 * (3 * t^2 - 3 * Tj1 * t + Tj1^2);
        p = [p q];
        v = v0 + alima * (t - Tj1 / 2);
        vc = [vc v];
        a = alima;
        ac = [ac a];
        jc = [jc 0];
    elseif t >= (Ta - Tj1) && t < Ta
        % 段3：减加速度段
        q = q0 + (vlim + v0) * Ta / 2 - vlim * (Ta - t) + jmax * (Ta - t)^3 / 6;
        p = [p q];
        v = vlim - jmax * (Ta - t)^2 / 2;
        vc = [vc v];
        a = jmax * (Ta - t);
        ac = [ac a];
        jc = [jc -jmax];
    elseif t >= Ta && t < (Ta + Tv)
        % 段4：匀速段
        q = q0 + (vlim + v0) * Ta / 2 + vlim * (t - Ta);
        p = [p q];
        v = vlim;
        vc = [vc v];
        a = 0;
        ac = [ac 0];
        jc = [jc 0];
    elseif t >= (T - Td) && t < (T - Td + Tj2)
        % 段5：加减速度段
        q = q1 - (vlim + v1) * Td / 2 + vlim * (t - T + Td) - jmax * (t - T + Td)^3 / 6;
        p = [p q];
        v = vlim - jmax * (t - T + Td)^2 / 2;
        vc = [vc v];
        a = -jmax * (t - T + Td);
        ac = [ac a];
        jc = [jc -jmax];
    elseif t >= (T - Td + Tj2) && t < (T - Tj2)
        % 段6：匀减速度段
        q = q1 - (vlim + v1) * Td / 2 + vlim * (t - T + Td) - alimd / 6 * (3 * (t - T + Td)^2 - 3 * Tj2 * (t - T + Td) + Tj2^2);
        p = [p q];
        v = vlim - alimd * (t - T + Td - Tj2 / 2);
        vc = [vc v];
        a = -alimd;
        ac = [ac a];
        jc = [jc 0];
    elseif t >= (T - Tj2) && t < T
        % 段7：减减速度段
        q = q1 - v1 * (T - t) - jmax * (T - t)^3 / 6;
        p = [p q];
        v = v1 + jmax * (T - t)^2 / 2;
        vc = [vc v];
        a = -jmax * (T - t);
        ac = [ac a];
        jc = [jc jmax];
    end
end

% 绘制位置曲线
subplot(4, 1, 1);
plot(0:0.1:T, p, 'LineWidth', 2);
xlabel('时间 (秒)');
ylabel('位置');
title('位置随时间变化');

% 绘制速度曲线
subplot(4, 1, 2);
plot(0:0.1:T, vc, 'LineWidth', 2);
xlabel('时间 (秒)');
ylabel('速度');
title('速度随时间变化');

% 绘制加速度曲线
subplot(4, 1, 3);
plot(0:0.1:T, ac, 'LineWidth', 2);
xlabel('时间 (秒)');
ylabel('加速度');
title('加速度随时间变化');

% 绘制加加速度曲线
subplot(4, 1, 4);
plot(0:0.1:T, jc, 'LineWidth', 2);
xlabel('时间 (秒)');
ylabel('加加速度');
title('加加速度随时间变化');

代码结果

在这里插入图片描述

五、simulink七段式s曲线

有小伙伴想在永磁同步电机或者其他电机的simulink仿真上使用，所以我在此修改，得到下面的结果。

整体框图

在这里插入图片描述

simulink结果

从0到400rpm，在从400rpm到0

在这里插入图片描述

s-function内嵌代码

有一个事情非常重要

最大加速度和最大加加速度越大，目标位置越小，可以使采样时间变小

这就说明你想关注 最大速度和采样时间就得注意这一点，增大或者减小最大速度就得变其他两个，而采样时间是使你的simulink仿真时间变短

下面这个函数名字是fun，所以你的simulink中s-function的名字也要设置一致

function [sys,x0,str,ts,simStateCompliance] = fun(t,x,u,flag)



switch flag,
 
  case 0,
    [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes;
 
  case 1,
    sys=mdlDerivatives(t,x,u);
    
  case 2,
    sys=mdlUpdate(t,x,u);
 
  case 3,
    sys=mdlOutputs(t,x,u);
 
  case 4,
    sys=mdlGetTimeOfNextVarHit(t,x,u);
 
 
  case 9,
    sys=mdlTerminate(t,x,u);
 
  otherwise
    DAStudio.error('Simulink:blocks:unhandledFlag', num2str(flag));
 
end
 
 
function [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes
 
sizes = simsizes;
 
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 1;
sizes.NumInputs      = 1;
sizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 1;  
 
sys = simsizes(sizes);
x0  = [];
str = [];
ts  = [0 0];
 
simStateCompliance = 'UnknownSimState';
 
function sys=mdlDerivatives(t,x,u)
 
sys = [];
 
 
function sys=mdlUpdate(t,x,u)
 
sys = [];
 
 
function sys=mdlOutputs(t,x,u)
% 定义参数
q0 = 0;     % 初始位置
q1 = 500;   % 目标位置
vmax = 400;  % 最大速度
amax = 5000;  % 最大加速度
v0 = 0;     % 初始速度
v1 = 0;     % 最终速度
jmax = 5000;  % 最大加加速度
count = 0;  % 计数器

 %% 最大加速度越大，目标位置越小，可以使采样时间变小
 
 %%
 
% 计算Ta、Td、Tj1、Tj2、Tv、alima、alimd等参数
if (vmax - v0) * jmax < amax^2
    if v0 > vmax
        Tj1 = 0;
        Ta = 0;
        alima = 0;
    else
        Tj1 = ((vmax - v0) / jmax)^0.5;
        Ta = 2 * Tj1;
        alima = Tj1 * jmax;
    end
else
    Tj1 = amax / jmax;
    Ta = Tj1 + (vmax - v0) / amax;
    alima = amax;
end

if (vmax - v1) * jmax < amax^2
    Tj2 = ((vmax - v1) / jmax)^0.5;
    Td = 2 * Tj2;
    alimd = Tj2 * jmax;
else
    Tj2 = amax / jmax;
    Td = Tj2 + (vmax - v1) / amax;
    alimd = amax;
end

Tv = (q1 - q0) / vmax - Ta / 2 * (1 + v0 / vmax) - Td / 2 * (1 + v1 / vmax);

% 动态调整参数，确保轨迹满足条件
if Tv <= 0
    Tv = 0;
    amax_org = amax;
    
    while Ta < 2 * Tj1 || Td < 2 * Tj2
        count = count + 1;
        amax = amax - amax_org * 0.1;
        alima = amax;
        alimd = amax;
        
        delta = (amax^4) / (jmax^2) + 2 * (v0^2 + v1^2) + amax * (4 * (q1 - q0) - 2 * amax / jmax * (v0 + v1));
        Tj1 = amax / jmax;
        Ta = (amax^2 / jmax - 2 * v0 + delta^0.5) / (2 * amax);
        Tj2 = amax / jmax;
        Td = (amax^2 / jmax - 2 * v1 + delta^0.5) / (2 * amax);
    end
end

% 计算轨迹
p = []; % 位置数据
vc = []; % 速度数据
ac = []; % 加速度数据
jc = []; % 加加速度数据

if Tv > 0
    vlim = vmax;
    T = Tv + Ta + Td;
else
    Tv = 0;
    vlim = v0 + (Ta - Tj1) * alima;
    T = Tv + Ta + Td;
end

    if u >= 0 && u < Tj1
        % 段1：加加速度段

        sys = v0 + jmax * u^2 / 2;

    elseif u >= Tj1 && u < (Ta - Tj1)
        % 段2：匀加速度段

        sys = v0 + alima * (u - Tj1 / 2);

    elseif u >= (Ta - Tj1) && u < Ta
        % 段3：减加速度段

        sys = vlim - jmax * (Ta - u)^2 / 2;

    elseif u >= Ta && u < (Ta + Tv)
        % 段4：匀速段

        sys = vlim;
    elseif u >= (T - Td) && u < (T - Td + Tj2)
        % 段5：加减速度段

        sys = vlim - jmax * (u - T + Td)^2 / 2;

    elseif u >= (T - Td + Tj2) && u < (T - Tj2)
        % 段6：匀减速度段

        sys = vlim - alimd * (u - T + Td - Tj2 / 2);

    elseif u >= (T - Tj2) && u < T
        % 段7：减减速度段

        sys = v1 + jmax * (T - u)^2 / 2;
else
    sys = 1;
end

 
 
function sys=mdlGetTimeOfNextVarHit(t,x,u)
 
sampleTime = 1;   
sys = t + sampleTime;
 
function sys=mdlTerminate(t,x,u)
 
sys = [];

六、使用双闭环FOC控制-----simulink仿真示意

在这里插入图片描述
输出得到转速图，如下所示：
在这里插入图片描述

小陈IT
关注
  • 19
    点赞
  • 踩
  • 40
    收藏
    觉得还不错? 一键收藏
  • 7
    评论
步进电机S型曲线算法模型
01-04
网上给出的步进电机S型的算法不是在起始和终止的输出值不是很精确，写的太复杂了，我研究了下发现看图像做简单的平移和拉伸就能解决，这里提供两个解决方案
速度曲线规划之S形曲线
Tunny_one
10-12 4736
速度曲线规划之S形曲线
用SCL编写完整程序实现堆垛机S型曲线速度控制，仿真测试可行，基于西门子S7-1500博途V15以上的堆垛机S型曲线速度控制部分完整程序，采用SCL高级编程语言，可仿真测试堆垛机速度曲线S曲线和梯形曲
最新发布
VfdSKTvuXUBx的博客
06-10 440
在博途V15以上版本中，我们可以打开编程环境，进行堆垛机速度控制程序的编译和调试，确保其在实际运行中的稳定性和性能。通过使用西门子博图1500和SCL编写的堆垛机S型曲线速度控制部分完整程序，我们可以实现堆垛机的平稳运动和高效操作。本文将围绕堆垛机速度曲线的S型曲线和梯形曲线展开，以及使用西门子博图1500和SCL编写堆垛机S型曲线速度控制的完整程序。在编写堆垛机的S型曲线速度控制部分程序时，我们可以利用博图1500的编程环境，采用SCL高级编程语言进行编写。堆垛机速度曲线S曲线 梯形曲线。
速度规划(梯形、余弦、多项式、7段S型、7段修正S型、15段S型、31段S型)
木牛的博客
11-04 7283
文章目录 一、7种速度规划对比 1、梯形速度规划 1.1 只有一段 1.1.1 加速段 1.1.2匀速段 1.1.3减速段 1.2 只有两段 1.2.1 加速段+匀速段 1.2.2匀速段 +减速段 1.2.3加速段+减速段 1.3 有三段(加速段+匀速段 +减速段) 2、余弦速度规划 2.1 只有一段 2.1.1 加速段 2.1.2 匀速段 2.1.3
S形曲线加减速--7段
hanmingjunv5的博客
05-22 5996
S加减速–7段 接前文5段S形速度规划，这里推导一下7段S形曲线的计算公式 1. S加减速曲线 定义： 上图中，令t1-t7为各个位置的绝对时间节点，T1-T7为每一段的增量时间节点。 这里的定义和图中有些出入，线下面的计算全部按照此定义为准. 计算公式： 加加速度 j(t)={J0≤t≤t10t1≤t≤t2−Jt2≤t≤t30t3≤t≤t4Jt4≤t≤t50t5≤t≤t6−Jt6≤t≤t7 j(t)=\begin{cases} J & 0\leq t \leq t_1\\ 0 &
分段式S形速度规划算法
02-18 1万+
1 典型七段式S形速度曲线 典型7段式S形曲线位置q、速度v、加速度a、加加速度j的波形如图所示： 位置q和时间t的关系式： 2 分段式S形曲线速度规划 已知初速度V0,末速度Ve,最大速度Vmax，路程距离s，最大加速度amax,最大jerk（加加速度）jmax,要规划位置、速度、加速度、加加速度随时间变化的曲线。具体步骤如下： 1） 假设能够达到最大速度Vmax，...
S形曲线加减速算法流程--7段式
hanmingjunv5的博客
05-24 1万+
7段S形加减速的算法流程. 7段S形加减速相对于5段，增加了匀加速过程，计算流程更加复杂一些，需要进行多次分类讨论。 加减速规划的流程： 曲线加减速规划是否存在匀速段取决于速度达到给定的指令速度口时，加速段和减速段位移之和s'是否小于待插补位移s． 加速及减速过程中是否存在匀加速段或匀减速段则由vs,ve,v,a,j决定。 规划后的S曲线可能出现下图中8中情况： 1.匀加速段存在的条件 若： v−vs>am2/jm v-v_s > a^2_m/j_m v−vs​>am2​/jm​
七段S形速度曲线轨迹规划
xueluowutong的博客
04-07 2万+
S形7段式曲线推导
qq_38601902的博客
08-18 276
c）减加速阶段：对Ta-t积分，速度v、加速度a、位移s通过计算加速段最大值，再减去Ta-t部分得到。d）匀速阶段：对速度曲线计算加速段面积和匀速段面积。7段式轨迹各个阶段公式推导。
步进电机驱动算法——S形加减速算法原理
热门推荐
aaa376070331的博客
06-29 3万+
目录步进电机S形加减速简介七段S形加减速算法原理分析五段S形加减速算法实现算法分析查表法编程实现思路及方法S曲线加减速流程图 步进电机S形加减速简介 一个物体从起点运动到终点要经历加速、匀速、减速的过程，把整个过程中速度随时间变化关系画出来，就是速度曲线。 “S”形加减速就是指速度曲线中加减速阶段的曲线呈现一个英文字母“S”型。 “S”形加减速的曲线并没有限定是什么曲线，所以可以是指数曲线，可以是正弦曲线，这里分析的“S”形加减速是基于加速度匀速变化的曲线。 基于加速度匀速变化的S形加减速可分为： 7段：加
步进电机S型曲线加减速算法.rar
07-07
步进电机的s型曲线的加减速算法，对整个s型曲线有ppt的形式进行详细的分析，希望能帮助大家更好的理解并更好的运用起来。
开环降压-升压转换器的数学建模：降压-升压转换器的数学表示-matlab开发
05-29
转换器的电路主要包括电感L、电容C、开关S、负载电阻Rload和二极管D1。通过基尔霍夫电流定律（KCL）和电压定律（KVL），我们可以得到以下两组方程： 1. 对于电感L： - 在开关S导通时（电感电流增加阶段）：Vin = ...
Hi-School-Differential-Equation:在CIE（剑桥国际考试）中求解高中微积分的一阶非线性ODE（常微分方程）。-matlab开发
05-28
MATLAB中的Simulink是系统级建模和仿真环境，尤其适用于处理动态系统的模型。在这个Hi-School-Differential-Equation项目中，Simulink模型被用来可视化地构建一阶非线性ODE的解决方案。通过拖放组件，连接它们来表示...
matlab仿真实验指导书.docx
09-14
【MATLAB仿真实验指导】 ...通过这两个实验，学生不仅可以掌握S函数的使用，还能深入理解电力系统中三相桥式整流器的工作原理和仿真技术。这样的实践将有助于提升在信号处理、控制系统设计和电力系统分析等方面的能力。
matlab仿真实验指导书.pdf
09-14
总的来说，MATLAB仿真是电气与电子工程领域中一种强大的工具，通过S函数和Simulink库，工程师和研究人员能够准确地模拟和分析复杂的系统行为。这两个实验提供了一个基础，帮助学习者深入理解S函数的使用和实际工程...
matlab仿真实验指导书..pdf
09-14
在本实验中，我们将探讨如何利用MATLAB的Simulink环境来实现S函数，以模拟单摆运动和晶闸管三相桥式整流器的工作。 **S函数详解** S函数，全称System Function，是Simulink中的一个重要组件，它允许用户自定义系统...
S形速度规划算法
01-14 1万+
S形速度规划相对于梯形速度规划其速度曲线会更加平滑，电机运行会更加平稳。常见的S曲线包括7段式加速度曲线，这种曲线计算量大，而且规划起来困难。一种简单的方法是采用Sigmoid对称函数加减速曲线规划法。典型的Sigmoid函数为： 其值域为(0,1)，函数关于横坐标左右对称，关于点（0,0.5）中心对称。其函数图形为： 如要将此曲线应用在步进电机的加、减速过程中，需要将方程在XY坐标系...
【STM32篇】步进电机之S型曲线
we_long的博客
03-14 6588
基于STM32 HR4988驱动控制。实现S曲线控制步进电机完成指定频率和范围的运行，实现”电机振动“。
matlab七段式s曲线表达式
07-29
MATLAB七段式S曲线表达式是一种用于描述连续曲线的函数表达式。它由多段分段线性函数组成，从而形成了平滑的曲线。 具体而言，七段式S曲线表达式由七个线性函数加上平滑因子组成。每个线性函数由两个控制点（起点和终点）决定，这些控制点可以通过用户输入来指定。平滑因子用于调整曲线的平滑度。 这个七段式S曲线表达式可以用方程表示为: y = a1*x^3+b1*x^2+c1*x + d1 (0 ≤ x ≤ x1) y = a2*x^3+b2*x^2+c2*x + d2 (x1 ≤ x ≤ x2) ... y = a7*x^3+b7*x^2+c7*x + d7 (x6 ≤ x ≤ x7) 其中，x是输入变量，y是输出变量，ai、bi、ci和di是每个线性函数的系数。 根据以上表达式，可以通过指定起始点、终止点和平滑度来生成S曲线。首先，需要将整个曲线从起始点到终止点等分为七个部分，每个部分对应一个线性函数。然后，通过求解每个线性函数的系数，可以得到完整的S曲线表达式。 通过这个七段式S曲线表达式，可以得到一条平滑的曲线，适用于各种应用领域，如图像处理、信号处理和控制系统等。在MATLAB中，可以使用曲线拟合和插值函数来实现对S曲线的生成和绘制。
评论 7
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值