2. 数学模型（一）：PMSM在三相静止坐标系下的电压方程与磁链方程
各位同学，咱们今天正式进入数学模型的环节。说实话，很多搞电机控制的同行，一看到数学公式就头大。我刚开始做项目那会儿也是这样，总觉得能调通代码就行，数学嘛，差不多得了。
结果呢？有一次调试一个高速主轴电机，转速一上去系统就震荡，我折腾了整整两周，最后才发现是模型里忽略了一个交叉耦合项。嗯，从那以后，我再也不敢轻视数学模型了。
2.1 为什么要从三相静止坐标系开始？
你想想看，电机本体长什么样？定子上三个绕组，空间上相差120度，通入三相交流电。这就是最自然的描述方式——三相静止坐标系（ABC坐标系）。
我个人习惯把这一步叫做"物理模型直译"。就是把电机的物理结构，直接翻译成数学方程。虽然这个模型用起来不方便（后面我们会做坐标变换），但它是最真实、最没有信息损失的描述。
核心观点：三相静止坐标系下的方程是所有后续模型推导的根基。你如果把这个搞透了，后面的Park变换、Clark变换、dq轴模型，本质上都是在做"数学上的等价变形"，物理本质没变。
2.2 电压方程——绕不过去的基尔霍夫
每一相绕组，说白了就是一个电阻串联一个电感线圈。根据基尔霍夫电压定律，我们可以写出A相、B相、C相的电压方程。
以A相为例：
u_A = R_s * i_A + d(ψ_A)/dt其中：
- u_A：A相绕组端电压（V）
- R_s：定子每相电阻（Ω）
- i_A：A相电流（A）
- ψ_A：A相绕组的总磁链（Wb）
B相和C相完全对称，我就不啰嗦了，直接给出三相电压方程的矩阵形式：
[u_A] [R_s 0 0 ] [i_A] d [ψ_A]
[u_B] = [0 R_s 0 ] [i_B] + -- [ψ_B]
[u_C] [0 0 R_s] [i_C] dt [ψ_C]这个方程看起来简单，但有个坑——磁链ψ_A、ψ_B、ψ_C并不是独立的。它们之间通过互感耦合在一起。我在项目中遇到过有人直接把三相当成三个独立的RL回路来算，结果仿真和实测差了30%以上。
注意：三相绕组之间存在互感，而且这个互感是随转子位置变化的！你不能把三相割裂开来看。这也是为什么三相静止坐标系下的方程用起来麻烦——耦合太严重了。
2.3 磁链方程——耦合的根源
磁链方程才是这一节的重头戏。每一相的总磁链由三部分组成：
- 自感磁链：本相电流产生的磁链
- 互感磁链：另外两相电流通过互感耦合过来的磁链
- 永磁体磁链：转子上的永磁体在定子绕组中感应的磁链
写成矩阵形式：
[ψ_A] [L_AA L_AB L_AC] [i_A] [ψ_fA]
[ψ_B] = [L_BA L_BB L_BC] [i_B] + [ψ_fB]
[ψ_C] [L_CA L_CB L_CC] [i_C] [ψ_fC]这里面的电感参数，才是真正让人头疼的地方。
2.4 电感参数——随转子位置变化的"怪物"
为什么说电感参数是怪物？因为对于表贴式PMSM（SPMSM）和内置式PMSM（IPMSM），电感的表现完全不同。
| 参数 | 表达式 | 说明 |
|---|---|---|
| 自感 L_AA | L_s0 + L_s2 * cos(2θ) | 含2次谐波，IPMSM明显，SPMSM近似常数 |
| 自感 L_BB | L_s0 + L_s2 * cos(2θ - 120°) | 相位偏移120° |
| 自感 L_CC | L_s0 + L_s2 * cos(2θ + 120°) | 相位偏移-120° |
| 互感 L_AB | -M_s0 + M_s2 * cos(2θ + 30°) | 互感也有2次谐波 |
| 互感 L_BC | -M_s0 + M_s2 * cos(2θ - 90°) | 依此类推 |
| 互感 L_CA | -M_s0 + M_s2 * cos(2θ + 150°) | 依此类推 |
你看，每个电感都包含一个常数项和一个随转子位置θ变化的2次谐波项。这就是为什么三相静止坐标系下的方程是时变系统——电感矩阵里的每个元素都在随时间变化。
我的经验：如果你在做参数辨识，千万别直接拿三相静止坐标系下的方程去辨识电感。因为L_AA、L_AB这些参数是时变的，你辨识出来的是一个随时间变化的曲线，很难直接用于控制器设计。我建议先变换到dq坐标系再辨识，会简单很多。
2.5 永磁体磁链——转子位置的正弦函数
永磁体在定子绕组中感应的磁链，是转子位置θ的正弦函数：
[ψ_fA] [cos(θ) ]
[ψ_fB] = [cos(θ - 120°) ] * ψ_f
[ψ_fC] [cos(θ + 120°) ]这里的ψ_f是永磁体磁链的幅值，对于给定的电机，它是一个常数（不考虑温度变化时）。
为什么会是余弦函数？因为永磁体产生的磁场在气隙中是正弦分布的，而定子绕组在空间上也是正弦分布的。两个正弦量相乘，自然得到余弦关系。
2.6 完整的数学模型——把方程合起来
把电压方程和磁链方程合在一起，就得到了PMSM在三相静止坐标系下的完整数学模型：
u_ABC = R_s * i_ABC + d/dt (L_ABC(θ) * i_ABC + ψ_fABC(θ))展开后：
u_ABC = R_s * i_ABC + L_ABC(θ) * d(i_ABC)/dt + d(L_ABC(θ))/dt * i_ABC + d(ψ_fABC(θ))/dt这里面有四项：
- 电阻压降：R_s * i_ABC
- 变压器电动势：L_ABC(θ) * d(i_ABC)/dt —— 电流变化引起的
- 运动电动势：d(L_ABC(θ))/dt * i_ABC —— 电感变化引起的（因为转子在转）
- 反电动势：d(ψ_fABC(θ))/dt —— 永磁体旋转引起的
你看，一个看似简单的电机，数学模型写出来其实挺复杂的。这也是为什么我们在实际工程中很少直接使用三相静止坐标系下的模型——太复杂了，而且电感矩阵需要实时更新，计算量很大。
小结一下：三相静止坐标系下的模型，物理意义清晰，但数学形式复杂。它是我们理解电机原理的基础，但不是工程实现的首选。下一节我们会学习Clark变换和Park变换，把这些时变系数变成常数，那才是工程上真正用的东西。
2.7 知识体系结构图
为了帮你理清这一节的知识脉络，我画了一张结构图：
2.8 避坑指南
最后，分享几个我在实际项目中踩过的坑：
- 别忽略互感：我曾经在搭建仿真模型时，为了省事只考虑了自感，结果仿真出来的电流波形和实测完全对不上。三相之间的互感不是小量，不能忽略。
- 注意电感参数的2次谐波：对于IPMSM，如果你用电感表去测量L_AA，你会发现它随着转子位置变化。这不是测量误差，是物理本质。做参数辨识时一定要考虑这个。
- 永磁体磁链的方向：ψ_f的正方向定义要和反电动势的正方向一致。我见过有人把符号搞反了，结果整个模型的反电动势方向都是错的，控制器输出自然也不对。
- 温度影响：R_s和ψ_f都受温度影响。R_s随温度升高而增大，ψ_f随温度升高而减小（永磁体退磁）。如果你在做高精度模型，一定要把温度效应加进去。
好了，这一节的内容就到这里。三相静止坐标系下的模型虽然在实际控制中不常用，但它是理解电机本质的必经之路。把这些方程搞明白，后面学坐标变换和dq轴模型，你会觉得豁然开朗。