第一部分：基础物理概念与坐标系
在推导公式前，必须建立坐标系的概念，这是现代电机控制的基石。
1. 坐标系定义
- 三相静止坐标系 ($abc$)：电机定子侧真实的物理坐标系，三相互差120度。
- 两相静止坐标系 ($\alpha\beta$)：通过Clark变换得到的静止坐标系，通常$\alpha$轴与$a$轴重合。
- 两相旋转坐标系 ($dq$)：这是最核心的坐标系。$d$轴（直轴）方向与转子磁链方向重合，$q$轴（交轴）超前$d$轴90度（电角度）。所有的解耦控制、MTPA、弱磁都在这个坐标系下进行。
2. 关键物理量符号
- $u_d, u_q$：定子电压在$d, q$轴上的分量。
- $i_d, i_q$：定子电流在$d, q$轴上的分量。
- $\psi_f$：永磁体磁链。
$L_d, L_q$：$d$轴电感和$q$轴电感。
- SPMSM (表贴式)：$L_d = L_q$，无凸极效应。
- IPMSM (内置式)：$L_d < L_q$，存在磁阻转矩，这是MTPA的基础。
- $R_s$：定子电阻。
- $\omega_e$：电角速度。
- $\omega_r$：机械角速度。
- $p_n$：极对数。
第二部分：核心数学模型
这是所有高级控制算法的“公理”。
1. 坐标变换公式
Clark变换 ($abc \to \alpha\beta$)（等幅值变换）：
$$ \begin{bmatrix} i_\alpha \\ i_\beta \end{bmatrix} = \frac{2}{3} \begin{bmatrix} 1 & -\frac{1}{2} & -\frac{1}{2} \\ 0 & \frac{\sqrt{3}}{2} & -\frac{\sqrt{3}}{2} \end{bmatrix} \begin{bmatrix} i_a \\ i_b \\ i_c \end{bmatrix} $$
Park变换 ($\alpha\beta \to dq$)：
$$ \begin{bmatrix} i_d \\ i_q \end{bmatrix} = \begin{bmatrix} \cos\theta & \sin\theta \\ -\sin\theta & \cos\theta \end{bmatrix} \begin{bmatrix} i_\alpha \\ i_\beta \end{bmatrix} $$
(注：$\theta$ 为转子电角度位置)
2. 电压方程
这是描述电机电气特性的核心方程：
$$ \begin{cases} u_d = R_s i_d + \frac{d\psi_d}{dt} - \omega_e \psi_q \\ u_q = R_s i_q + \frac{d\psi_q}{dt} + \omega_e \psi_d \end{cases} $$
在稳态或忽略电阻压降时，简化为：
$$ \begin{cases} u_d = -\omega_e \psi_q \\ u_q = \omega_e \psi_d \end{cases} $$
3. 磁链方程
描述电流如何产生磁场：
$$ \begin{cases} \psi_d = L_d i_d + \psi_f \\ \psi_q = L_q i_q \end{cases} $$
4. 转矩方程
这是MTPA控制的目标函数来源：
$$ T_e = \frac{3}{2} p_n (\psi_d i_q - \psi_q i_d) $$
将磁链方程代入，得到含磁阻转矩的完整转矩方程：
$$ T_e = \frac{3}{2} p_n [\psi_f i_q + (L_d - L_q) i_d i_q] $$
解析：
- 第一项 $\psi_f i_q$ 是永磁转矩。
- 第二项 $(L_d - L_q) i_d i_q$ 是磁阻转矩（只有IPMSM才有，因为 $L_d \neq L_q$）。
5. 运动方程
连接电磁转矩与机械转速：
$$ J \frac{d\omega_r}{dt} = T_e - T_L - B \omega_r $$
- $J$: 转动惯量
- $T_L$: 负载转矩
- $B$: 摩擦系数
第三部分：MTPA（最大转矩电流比）控制模型
MTPA的目标是：在给定转矩需求下，使定子电流幅值最小（即效率最优）。
1. 数学推导依据
构建拉格朗日函数或直接求极值。电流幅值平方：
$$ i_s^2 = i_d^2 + i_q^2 $$
目标：在 $T_e$ 恒定条件下，求 $i_s$ 的最小值。
2. MTPA 电流轨迹方程
对转矩方程求导并令导数为零，推导出 $i_d$ 和 $i_q$ 的关系：
$$ i_d = \frac{\psi_f - \sqrt{\psi_f^2 + 4(L_d - L_q)^2 i_q^2}}{2(L_d - L_q)} $$
或者用电流角 $\beta$（电流矢量与 $q$ 轴夹角）表示：
$$ \beta = \arcsin\left( \frac{-\psi_f + \sqrt{\psi_f^2 + 8(L_d - L_q)^2 i_s^2}}{4(L_d - L_q) i_s} \right) $$
- 控制策略：在低速区，根据目标转矩，通过查表或计算上述公式，给出 $i_d^{ref}$ 和 $i_q^{ref}$。
第四部分：弱磁控制模型
弱磁控制的目标是：突破基速，扩展电机运行范围。
1. 电压极限圆
逆变器输出电压有上限 $u_{max}$（通常为 $U_{dc}/\sqrt{3}$）。忽略电阻，稳态电压方程近似为：
$$ u_s = \sqrt{u_d^2 + u_q^2} \le u_{max} $$
代入磁链方程，得到电压极限椭圆方程：
$$ (L_d i_d + \psi_f)^2 + (L_q i_q)^2 \le \left(\frac{u_{max}}{\omega_e}\right)^2 $$
- 概念：随着转速 $\omega_e$ 增加，椭圆半径缩小，电流工作点被“挤压”向 $-i_d$ 方向。
2. 电流极限圆
$$ i_d^2 + i_q^2 \le i_{max}^2 $$
3. 弱磁控制策略
当电机转速升高，反电动势接近母线电压，电压饱和。此时需要：
- 控制 $u_s = u_{max}$（电压闭环）。
调节 $i_d$ 向负方向增加（去磁电流），以抵消永磁体磁链，降低反电动势。
$$ i_d = -\frac{\psi_f}{L_d} + \frac{1}{L_d} \sqrt{\left(\frac{u_{max}}{\omega_e}\right)^2 - (L_q i_q)^2} $$
第五部分：滑模控制 (SMC) 相关模型
SMC用于提高控制器的鲁棒性，常用于速度环或位置环，替代传统的PI控制。
1. 状态空间模型
定义状态变量为转速误差：$x = \omega_{ref} - \omega_r$。
对时间求导，代入运动方程：
$$ \dot{x} = \dot{\omega}_{ref} - \dot{\omega}_r = 0 - \frac{1}{J}(T_e - T_L - B\omega_r) $$
代入转矩方程，得到误差状态方程（这是设计滑模面的基础）。
2. 滑模面设计
通常采用积分滑模面，以减小抖振：
$$ s = e + c \int_0^t e \, dt $$
其中 $e$ 为误差（如转速误差），$c > 0$。
3. 滑模控制律
设计控制律使其满足可达性条件（Lyapunov稳定性）：
$$ \dot{s} \cdot s < 0 $$
典型的指数趋近律：
$$ \dot{s} = -\epsilon \cdot \text{sgn}(s) - k s $$
- $\epsilon$: 切换增益，需大于系统扰动（负载变化）的上界。
- $\text{sgn}(s)$: 符号函数。
- $k$: 指数趋近速率。
总结：你需要掌握的知识图谱
为了掌握这些高级算法，你需要按以下顺序攻克数学模型：
- 基础层：Park变换、电压方程、转矩方程。这是FOC（矢量控制）的基础。
- 进阶层 (MTPA)：理解 $L_d \neq L_q$ 带来的磁阻转矩，掌握 $i_d, i_q$ 的最优配比推导。
- 高速层：深刻理解电压极限圆和电流极限圆的几何意义，理解为什么需要负的 $i_d$。
- 算法层 (SMC)：掌握 Lyapunov 稳定性理论，理解符号函数 $\text{sgn}(s)$ 的作用及抖振问题。
建议你将这些公式手写推导一遍，特别是转矩方程的分解和电压极限圆的推导，这是理解电机控制从“理论”到“实践”跨越的关键。