2. 电机控制系统建模基础
电机控制系统仿真的第一步，不是搭控制器，而是把电机本体“翻译”成仿真器能计算的数学对象。对三相交流电机而言，这套数学对象由三组方程构成：电压方程描述绕组端电压与电流、磁链的关系，磁链方程描述各绕组磁链由哪些电流和电感产生，转矩方程描述电流与磁链如何生成电磁转矩。三者联立，再配合机械运动方程，才构成完整的电机模型。
一、三相静止坐标系下的原始方程
以定子三相绕组 A、B、C 和转子绕组为例，电压方程可写成矩阵形式：
u = R·i + dψ/dt
其中 u = [uA uB uC]ᵀ，i = [iA iB iC]ᵀ，ψ = [ψA ψB ψC]ᵀ
R = diag(Rs, Rs, Rs)磁链方程则体现绕组间的互感耦合：
ψ = L(θr)·i
L 中包含自感 Ls 与互感 M，且互感随转子位置 θr 周期变化转矩方程由磁场储能对机械角求偏导得到：
Te = (1/2)·iᵀ·(∂L/∂θm)·i这三组方程在物理上完全正确，但直接用于仿真有两个致命问题：一是电感矩阵随 θr 变化，每个仿真步都要重算矩阵；二是三相变量之间存在强耦合，控制器设计无从下手。因此必须做坐标变换。
二、Clarke 变换与 Park 变换
Clarke 变换把三相静止坐标系 (A,B,C) 投影到两相静止坐标系 (α,β)。采用等幅值变换时：
iα = (2/3)·(iA - 0.5·iB - 0.5·iC)
iβ = (2/3)·(√3/2·iB - √3/2·iC)Park 变换再把 (α,β) 旋转到随转子磁场同步旋转的 (d,q) 坐标系：
id = iα·cosθe + iβ·sinθe
iq = -iα·sinθe + iβ·cosθe其中 θe 为转子电角度。变换后，原本时变的电感矩阵变为常数矩阵，磁链方程简化为：
ψd = Ld·id + ψf
ψq = Lq·iq对表贴式永磁同步电机，Ld ≈ Lq = Ls，转矩方程退化为：
Te = (3/2)·p·ψf·iq这意味着：控制 iq 就能线性控制转矩，id 则用于弱磁或维持零。这正是矢量控制的理论根基，也是仿真中必须先把变换关系写对的原因。
参照约定：本节统一采用等幅值 Clarke 变换，d 轴对齐转子永磁磁链方向，q 轴超前 d 轴 90° 电角度，θe 为 d 轴与 α 轴的夹角。全文及图示均按此约定，切换约定需重新推导系数。
三、机电能量转换在仿真中的表达
电机本质是“电—磁—机”能量转换器。仿真中，这一转换通过两个接口变量体现：电磁转矩 Te 从电气侧输出到机械侧，转速 ωm 从机械侧反馈回电气侧。机械运动方程为：
J·dωm/dt = Te - TL - B·ωm
dθm/dt = ωm
θe = p·θm其中 J 为转动惯量，TL 为负载转矩，B 为阻尼系数，p 为极对数。仿真时，电气方程与机械方程通常以固定步长联立求解，步长需小于电气时间常数与机械时间常数中较小者的 1/10。
下面这张图把上述建模链路串起来：从三相电压输入，经电压方程得到电流，经 Clarke/Park 变换进入 dq 轴，再由磁链与转矩方程输出 Te，最后经机械方程反馈转速与角度，形成闭环。
四、仿真实现中的关键细节
把上述方程落到仿真里，有几个容易出错的地方值得提前说明。
第一，变换系数要统一。Clarke 变换有等幅值和等功率两种版本，前者系数为 2/3，后者为 √(2/3)。如果变换用等幅值、转矩公式却用等功率系数，Te 会差 1.5 倍。建议在模型里显式写出系数，不要凭记忆。
第二，角度积分要防溢出。θe 由 ωm 积分得到，长时间仿真会无限增大。实际实现时应对 θe 取模 2π，或使用 sin/cos 查表时自动处理周期性，避免数值精度下降。
第三，磁链微分不能直接对 ψ 求导。电压方程中的 dψ/dt 在仿真中通常改写为：
dψd/dt = ud - Rs·id + ωe·ψq
dψq/dt = uq - Rs·iq - ωe·ψd其中 ωe = p·ωm 为电角速度。这样写把旋转反电势项显式分离出来，数值求解更稳定，也便于观察交叉耦合对电流的影响。
第四，负载转矩要区分符号。TL 为阻碍运动的转矩，在运动方程中取负号；若负载为风机类，TL 与 ωm² 成正比，仿真中需用乘法模块而非常数。
教学示例参数（合理假设，非实测）：表贴式永磁同步电机，Rs = 0.5 Ω，Ld = Lq = 5 mH，ψf = 0.1 Wb，p = 4，J = 0.001 kg·m²，B = 0.0001 N·m·s。给定 ud = 0，uq = 10 V，TL = 0.5 N·m。仿真步长 1e-5 s。可观察到 iq 在约 5 ms 内上升至稳态，Te 随之建立，转速按机械时间常数缓慢上升。若把 Ld 改为 3 mH、Lq 保持 5 mH，则出现磁阻转矩分量，Te 表达式中的 (Ld-Lq)·id·iq 项不再为零。
五、从方程到可运行模型的检查清单
在把模型接入控制器之前，建议按以下顺序自查，避免把建模错误误判为控制问题：
- 电压方程中的电阻压降、磁链微分、旋转反电势三项是否齐全，符号是否与所选 dq 轴方向一致；
- Clarke/Park 变换系数是否与转矩公式中的 3/2 或 √3 系数匹配；
- θe 是否由 θm 乘极对数得到，而非直接使用机械角；
- 机械方程中 TL 的符号是否与“阻碍运动”定义一致；
- 仿真步长是否小于电气时间常数 Ls/Rs 的 1/10，本例中 Ls/Rs = 10 ms，步长 1e-5 s 满足要求；
- 初始条件是否合理：id、iq 初值为零，ωm 初值为零，θe 初值为零。
完成以上检查后，电机本体模型即可作为被控对象接入电流环、速度环。下一阶段的所有控制器设计，都建立在本节这组方程和变换关系之上；方程写错，后面的仿真结果再漂亮也没有意义。