伺服电机正反转Simulink仿真建模与PI参数整定详解

简介：基于MATLAB/Simulink的伺服电机正反转控制仿真脚本，面向自动化设备、机器人系统相关专业的学生与工程师，旨在解决缺少实物时无法直观理解伺服电机建模、闭环控制及换向逻辑的问题。压缩包内仅含1个m脚本，整体约1KB，脚本将电机数学模型、PID控制器、输入信号生成和动画显示代码整合在一起，仿真动画能够直观呈现电机在正反转指令下的速度、位置与转向变化。已有300人学习下载，读者可借此脚本快速搭建仿真学习环境，观察不同控制参数对响应曲线的影响，并通过修改信号或参数验证正反转控制策略；它也能作为课程设计或项目预研的起点，在轻量代码基础上扩展电机负载、扰动等场景，或进一步研究PID参数整定与换向时序。

1. 一个 zip 包背后的伺服电机仿真链路

打开一个 t6rey.zip 这类伺服电机仿真包，里面多半是 .slx 模型、参数脚本和一堆 .mat 结果文件。跑通的人很多，但真正把“正反转”这个动作做干净的人很少：仿真正反转不是在速度给定上换个负号就完事，逆变器换相、电流环反向、角度反馈符号，任何一环没对齐，转速就会在零点附近爬台阶。下面直接把最关键的环节捋一遍——伺服电机数学模型怎么进 MATLAB/Simulink、电流环和速度环怎么搭、正反转信号时序怎么给，以及最后用哪些波形确认换向正确。适合正在写课程设计、毕业论文，以及要在产线验证阶段快速拿出可交付结果的人。

2. 伺服电机仿真的数学基础：电角度、转矩方向与求解器选择

2.1 三个方程决定正反转能不能转起来

伺服电机的电气动态和机械动态，在 MATLAB 仿真里最终都收敛到三个方程。无论你用的是 Simulink 里的 PMSM 模块，还是自己写 S 函数，下面这三组式子都绕不开。

定子电压方程在 d-q 旋转坐标系下写成：

ud = Rs*id + Ld*d(id)/dt - ωe*Lq*iq uq = Rs*iq + Lq*d(iq)/dt + ωe*(Ld*id + ψf)

ωe 是电角速度，与机械角速度 ωm 的关系是 ωe = p*ωm，p 为极对数。这个换算直接决定角度反馈进坐标变换前要不要乘 p——很多正反转异常都出在这个系数上。

电磁转矩方程为：

Te = 1.5*p*(ψf*iq + (Ld - Lq)*id*iq)

对表贴式永磁同步电机，Ld 与 Lq 近似相等，转矩就基本正比于 iq。给负的 iq，转矩反向，电机才能反转。所以正反转的本质不是“倒换电源线”，而是电流环要把 iq 从正控到负。

机械方程：

J*d(ωm)/dt = Te - B*ωm - TL

J 是转子惯量，B 是粘滞摩擦系数，TL 是负载转矩。仿真里最容易造假的就是 B：设成 0 虽然转速曲线更“好看”，但过零时的动态响应会明显失真，正反转落点会和实际台架对不上。

2.1.1 为什么不能只改三相相序

有变频器/逆变器背景的人会问：工频电机反转换两根相线就行，仿真里直接把三相电压相序倒过来不行吗？

问题在于伺服电机是电子换向。PMMSM 模块内部根据转子位置实时计算换向逻辑，你从外部直接反相序，只会让电流环的角度解算与真实电角度错位，转矩方向与期望方向互相拉扯。正确做法是改变转速参考的符号，让速度环输出负 iq，再由电流环和逆变器完成内部换向。这部分对应标题里“正反转”的真实语义，也是后续所有配置的出发点。

2.2 建模前先定一套能用的电机参数

没有实测数据时，我一般直接用一套典型 PMSM 小功率参数起模型。这套参数接近教材和课程设计中常见的额定 4000 rpm、额定转矩约 1.2 N·m 的小型伺服：

这套参数可以直接输入 Simulink 的 PMSM 模块，也可以方便地换算成传递函数做线性分析。摩擦系数 B 这个值用来让仿真稳定即可，真实值需要通过空载减速曲线辨识，后面第 5 章会展开。

2.2.1 阶跃响应先过一遍

参数定完不要急着搭闭环，先用阶跃响应验证模型本身。把 q 轴电压突加一个恒定值，观察转速上升曲线，机械时间常数大约等于 J 除以等效阻尼。这是判断“模型参数是否合理”最快的手段，比直接跑完整闭环更容易定位问题。在 MATLAB 里可以对简化模型做线性化：[num, den] = ss2tf(A, B, C, D)，然后step(tf(num, den))直接看到机械环节的阶跃响应时长。如果你连闭环电流环一起阶跃，还能顺便识别出电流环带宽是否符合预期。

2.3 求解器与步长：逆变器模型决定取舍

Simulink 中伺服模型选连续求解器还是离散求解器，取决于你是否把逆变器功率电路建出来。

只做控制策略验证，PMSM 模块用平均值模型或直接把逆变器理想化，ode45 就够了，计算快、波形平滑。一旦加入 PWM 载波或 IGBT 开关管模型，开关频率 10 kHz 意味着仿真里需要分辨 1e-4 秒级的开关事件，此时 ode15s 或定步长离散求解器更稳。我一般设载波 10 kHz、采样时间 1e-5 秒，一个 0.6 秒的正反转工况仿真大约几十秒能跑完。

提示：变步长求解器跑带开关模型时，如果看到“仿真速率过慢”或波形在开关沿附近抖动，先把容差从 1e-3 调到 1e-4，再不行就换定步长。定步长值取开关周期的 1/10 是常见起点，即 10 kHz 对应 1e-5 秒。

求解器选型没有绝对正确，判断标准是你关心什么：关心转矩脉动和换向电流尖峰，必须保留开关细节；只关心速度环正反转曲线，平均值模型加 ode45 完全够用。

3. 用 Simulink 搭建可复现的伺服电机仿真模型

3.1 最小完整链路：从电源到测量

一个能跑正反转的控制模型，最少需要这几个环节串联：直流电源、三相逆变器、PMSM、负载转矩、测量与反馈。模块可以从 Simscape Electrical 的 Specialized Power Systems 库中拖出，近几年无论是 R2023b 还是新装用户关注的 2026b，路径基本稳定在“Electric Drives > Machines > Permanent Magnet Synchronous Machine”。

推荐的模块连接顺序：

DC Voltage Source → Universal Bridge → PMSM ↓ Torque Input (TL) ↓ Bus Creator / Voltage-Current Measurement → Scope/To Workspace

这里的 Universal Bridge 用三相桥臂，开关器件选 IGBT/Diodes，载波频率设 10 kHz。PMSM 模块的 Mechanical input 选 Torque TL，用一个常量或阶跃信号接在输入端作为负载。千万不要把机械负载直接挂在转速端口上——那是错误的物理含义。

3.2 PMSM 模块对话框里必须逐项确认的 8 个参数

在模块双击对话框中，参数要和第 2 章的表格一一对应。容易出错的是两个地方：

第一，Initial rotor angle 初始转子角度。仿真开始时如果为 0，对齐坐标系没问题；但如果你想模拟任意启动位置，就要在外部用一个常量连接到模块的 angle 输入端口。正反转测试时建议从 0 开始，减少坐标变换初始跳变对判断的干扰。

第二，Torque 端口的方向。模块默认正方向是“电磁转矩驱动转子正向旋转”，负载转矩输入正值代表阻力矩。反转工况下，电磁转矩要克服的是负方向阻力矩，所以负载转矩信号的符号不要随便改，保持物理方向一致。

电磁转矩和机械角速度的测量端口，我习惯用 Demux 分开，或者用 Bus Selector 挑出需要的两路再接 To Workspace。这样第 4 章的运行脚本才能干净地拿到数据。

3.2.1 反电动势波形选正弦还是梯形

PMSM 模块里通常有 Permanent magnet flux 和 Back EMF waveform 两个选项。永磁同步伺服电机基本都选正弦波反电动势；无刷直流电机才选梯形波。标题既然是伺服电机，这里不要选错，否则电流波形会多出明显的高频分量，正反转特点会被这些谐波掩盖。

3.3 电流环 PI 参数：用零极点抵消一次算到位

电流环结构常见做法是：把三相电流 Clarke/Park 变换到 d-q 轴，分别用两个 PI 控制器控 id 和 iq，输出再接逆变换到三相电压。这个环节在 Simulink 里可以用现成的 Current Controller 模块，也可以手搭 PI 加限幅。

PI 参数我用零极点抵消法，先定 q 轴电流环带宽 αi：

Kp = αi * Lq Ki = Kp * Rs / Lq

取 αi = 1500 rad/s，对应第 2 章的参数：

Kp = 1500 * 0.00525 = 7.875 Ki = 7.875 * 0.958 / 0.00525 ≈ 1437

这里 Kp 的量纲是 V/A，Ki 是 V/(A·s)。αi 的选择有约束：不能超过 PWM 开关频率的十分之一，10 kHz 载波对应 αi 上限约 6000 rad/s，1500 到 3000 是稳妥区间。αi 太大，数字控制一拍延迟带来的相位滞后会让电流环振荡；太小，电流响应慢，正反转时 iq 跟不上给定，转速过零会拖出明显圆弧。

d 轴电流环参数与 q 轴相同。对表贴式电机，id 给定直接设 0，这个环的主要作用是抑制外部扰动。

3.4 用初始化脚本把参数一次性写进工作区

手动在模块对话框一个个填参数既慢又容易错。常见做法是写一个参数脚本，在启动 Simulink 模型前先运行。模型的模块参数里直接填变量名即可。

% servo_param.m % PMSM 参数，先运行本脚本再打开模型 Rs = 0.958; Ld = 0.00525; Lq = 0.00525; flux = 0.1827; % Wb J = 0.0006329; % kg.m^2 B = 0.0001; % N.m.s p = 4; % 极对数 Vdc = 300; % 直流母线电压，逆变器输入端 % 电流环 PI（零极点抵消） alpha_i = 1500; % 电流环带宽 rad/s Kp_i = alpha_i * Lq; Ki_i = Kp_i * Rs / Lq; % 速度环初始值，仅作为起步，后续按阶跃响应微调 alpha_s = 80; % 速度环带宽 rad/s Kp_s = alpha_s * J * alpha_i / (1.5 * p * flux); Ki_s = alpha_s * Kp_s / 5; % 限幅 iq_max = 10; % q 轴电流限幅 A u_max = Vdc / 2; % 电压限幅 V

脚本的逻辑说明：变量名与模型对话框里填的名字必须完全一致，Simulink 在仿真开始时从 MATLAB 工作区解析这些变量。速度环 Kp 的公式来自“电流环闭环后近似一阶惯性环节 αi”的前提，把速度环被控对象简化成：转矩常数 1.5·p·ψf，再除以 Js。Kp_s 表达式里除以 5 是预留的相位裕度折减，不同模型可能需要调整。脚本里同时定义了电流限幅，Simulink 的饱和模块直接引用 iq_max 和 -iq_max，这样后面测正反转时，反向电流能一眼看出有没有碰到限幅边界。

4. 正反转的信号时序：给定、限幅与换向验证

4.1 先分清转速环里“谁在反转”

伺服系统里的反转，在仿真中可以发生在三个层面。

第一层是参考输入反转，即你在外部给速度环一个负转速，这是我们要验证的主要目标。第二层是速度环输出的 iq 参考值从正变负，它说明控制器已经决定反向出力。第三层是电机实际的电角度换向，这一层由逆变器和坐标变换共同完成，仿真中通常不会出错，但如果你的编码器模型漏掉了方向位，反馈角速度符号不对，控制器永远追不上目标。

排错时先看 iq_ref 有没有变号，再去看实际转速。如果 iq_ref 变号了但转速不反向，问题在电机/负载环节；如果 iq_ref 都不变号，问题一定在反馈或控制器环节。

4.2 用阶跃和梯形信号构造双向测试工况

正反转测试不能只给一个负阶跃，常见的做法是给一段含“正转—停止—反转—停止”的时序。Simulink 里可以简单用 Signal Builder 或 Signal Editor 画梯形速度曲线，也可以用 Step 模块叠加。下面是一组我常测试的工况：

注意 0.3 秒的负阶跃，速度环要给负的 iq 参考才能把转速从正向拉到反向。如果这里速度环 PI 或饱和限幅配置错误，最常见的现象是转速先停住，然后停在零点附近振荡很久才反向，这就是“换向爬台阶”。

阶跃测试在 MATLAB 里还可以直接对简化模型求闭环阶跃响应，用stepinfo看超调量和调节时间，把这两个指标当正反转工况的验收基线，比人眼判断波形可靠得多。

4.3 跑一次仿真并导出正反转波形的完整脚本

模型里把给定转速接到一个 Constant 模块，Constant 的值设为变量 n_ref，然后把输出的转速和转矩接到 To Workspace 模块，Save format 设为 Timeseries。运行脚本如下：

% run_reversal.m % 用 SimulationInput 覆盖工作区变量，跑 0.6s 正反转仿真 in = Simulink.SimulationInput('servo_pmsm_model'); in = in.setModelParameter('StopTime', '0.6'); in = in.setVariable('n_ref', 1000); % 覆盖 Constant 里的 n_ref out = sim(in); t = out.tout; wm = out.yout{1}.Values.Data; % 机械角速度 rad/s te = out.yout{2}.Values.Data; % 电磁转矩 N.m rpm = wm * 60 / (2 * pi); % rad/s 转 rpm figure subplot(2,1,1) plot(t, rpm); ylabel('转速 rpm'); grid on; title('正反转转速响应') subplot(2,1,2) plot(t, te); ylabel('转矩 N.m'); xlabel('时间 s'); grid on

逐行说明：setModelParameter 覆盖停止时间；setVariable 把 n_ref 注入工作区并覆盖同名变量，这是 Simulink 推荐的参数化仿真方式，比在脚本里直接n_ref=1000; sim('模型名')更干净，因为不会污染全局工作区。out.yout 的序号顺序取决于模型中 To Workspace 模块的排列顺序，不要靠猜，跑完后先out.yout.getElementNames看一眼再取数。转速换算因子 60/(2π) 是 rad/s 到 rpm 的标准换算。

这个脚本跑完后，重点看 subplot(2,1,1) 在 0.3 秒附近的斜率是否连续，以及 subplot(2,1,2) 的转矩有没有在过零时产生来回跳变。阻尼不足的 PI 参数会让转速在零点附近来回穿越，产生“反反复复正转反转”的假象。

4.4 三种经典反转失败模式

4.4.1 指令反转但 iq 不反向

现象：给定转速已经变负，但 iq_ref 仍然是正饱和值。检查速度环输出限幅——如果限幅模块的上下界用了变量但脚本没定义，Simulink 会在仿真开始时报错，但如果你写成常量限幅且只允许正向输出，就会出现这种“静默失败”。正反转测试前，把限幅改成对称的 ±iq_max。

4.4.2 角度反馈少了极对数

现象：空载正转正常，一给负给定就振荡甚至发散。原因是坐标变换使用的电角度写成了机械角度，或者反向时角度符号没有翻转。用 Scope 同时观察 rotor angle 输入和实际转速，如果角度斜率方向与给定方向不一致，问题就在这一环。

4.4.3 逆变器死区时间设得过大

现象：过零时转矩波形出现很短的凹坑，转速响应有一小段“滞后”。逆变器死区一般设 2~5 μs，对 10 kHz 载波来说不能超过开关周期的 5%。正反转过零处电流接近 0，死区引起的电压误差在这个区间最明显。这不是控制参数问题，改小死区或改用平均值逆变器模型即可验证。

5. 让正反转更稳的三个技巧：摩擦辨识、S 函数与换向自检

5.1 摩擦模型选对，过零点才真实

第 2 章里把 B 设为小常数，只是起步。正反转过零最敏感的是静摩擦和库仑摩擦。Simulink 里可以用 MATLAB Function 模块写一个改进摩擦转矩输入到负载端口，代替常值负载：

function Tf = friction(wm) Tc = 0.05; % 库仑摩擦 N.m Ts = 0.12; % 静摩擦 N.m ws = 1e-3; % 零速切换阈值 rad/s if abs(wm) < ws Tf = Ts * sign(wm + eps); % 零速附近用静摩擦 else Tf = Tc * sign(wm); end end

注意这个模型在零速附近的符号切换会造成仿真抖动，实际测试时可以在 sign 函数输出后串联一个低通滤波，或者用双曲正切近似。摩擦模型的意义在于：速度环的输出在过零瞬间要先克服静摩擦，正反转的“爬台阶”现象有一部分确实是真实物理特性，不要全部当成控制器问题。

5.2 用 S 函数扩展非线性负载和扰动观测器

如果要在仿真里加入带传动间隙、弹性联轴器这类非线性环节，Simulink 基础模块会变得很乱，我一般写一个 C MEX S 函数或 Level-2 MATLAB S 函数。在 S 函数的 Derivative 部分把机械方程扩展成二阶：把负载惯量分成电机侧和负载侧，中间用扭转弹簧连接。

S 函数的好处是可以直接加扰动观测器逻辑，实时估计负载转矩并前馈补偿。正反转时扰动的方向不变或缓慢变化，但速度给定在快速反向，这时前馈补偿的价值最明显。这一步做扎实以后，把参数交给产线同事做台架验证时，仿真和实测的转速曲线能对到相当高的置信度。

5.3 方向一致性自检：让模型自己判断换向是否正确

人眼看波形容易漏掉短暂的符号不一致。我习惯在模型里嵌入一个 MATLAB Function 模块，实时检查转速实际方向与指令方向：

function flag = dir_check(n_ref, wm) tol = 10; % 容差 rpm if n_ref > tol && wm < 0 flag = 1; % 实际反转，期望正转 elseif n_ref < -tol && wm > 0 flag = 2; % 实际正转，期望反转 else flag = 0; end end

这个模块的输出接一个 Scope，仿真结束后统计 flag 非零区间的时长。期望值与实际值符号不一致的总时长远小于整个换向过程时，可以判断为基本正常；一旦过长，说明速度环响应滞后明显，或者电机模块的机械端口配置错误。后续如果要更精细地整定速度环 PI，可以试试用优化工具箱里现成的fmincon对stepinfo的超调量和调节时间做离线寻优，把人工调参的来回次数降到最低。这样整套仿真从正反转验证到参数交付，链路就完整了。

本文还有配套的精品资源，点击获取