1. SVPWM在永磁同步电机控制系统中的定位与工程意义

空间矢量脉宽调制（Space Vector Pulse Width Modulation，SVPWM）并非一种孤立的波形生成技术，而是PMSM（永磁同步电机）矢量控制系统中连接上层控制算法与底层功率执行机构的关键枢纽。其核心作用在于：将由FOC（磁场定向控制）环路计算出的、以αβ坐标系表示的期望电压矢量 $ \mathbf{V} {ref} = [U \alpha, U_\beta]^T $，精确地映射为三相全桥逆变器所能理解的八种开关状态组合，并通过时间维度上的动态加权，最终在电机定子绕组上合成出幅值与相位连续可调的旋转电压矢量。

这一映射过程的必要性源于电力电子器件的物理本质。一个标准的三相两电平电压源型逆变器（VSI）由三个桥臂、六个功率开关管（通常为IGBT或MOSFET）构成。由于同一桥臂上下两个开关管必须遵循“互补导通、死区保护”的安全约束，以避免直通短路，因此其所有合法的开关组合仅有 $ 2^3 = 8 $ 种。这8种组合在电机定子三相坐标系中，对应着8个离散的电压空间矢量：6个非零矢量 $ \mathbf{U}_1 $ 至 $ \mathbf{U}_6 $，以及2个零矢量 $ \mathbf{U}_0 $（全上桥臂关断）和 $ \mathbf{U}_7 $（全上桥臂导通）。

若直接将这8个固定矢量按固定时长顺序施加于电机，所合成的定子磁链轨迹将是一个正六边形。从电磁理论出发，定子绕组上的电压 $ \mathbf{u}(t) $ 与由此产生的磁链 $ \boldsymbol{\psi}(t) $ 满足基本关系：
$$
\frac{d\boldsymbol{\psi}(t)}{dt} = \mathbf{u}(t) - R_s \mathbf{i}(t)
$$
在忽略电阻压降 $ R_s \mathbf{i}(t) $ 的高速运行近似下，该式可简化为 $ d\boldsymbol{\psi} \approx \mathbf{u} \cdot dt $。这意味着，一个恒定电压矢量 $ \mathbf{U}_k $ 作用时间 $ \Delta t $，将在磁链平面上产生一个大小为 $ |\mathbf{U}_k| \cdot \Delta t $、方向与 $ \mathbf{U}_k $ 一致的增量 $ \Delta \boldsymbol{\psi} $。当依次轮换施加 $ \mathbf{U}_1 $ 至 $ \mathbf{U}_6 $，每个作用 $ \Delta t $，其首尾相接的磁链增量便勾勒出一个正六边形轨迹。

这个六边形轨迹是系统性能的瓶颈所在。理想的磁场定向控制要求定子磁链以恒定角速度、恒定幅值作圆周运动，从而在气隙中建立一个平滑、均匀的旋转磁场，驱动转子实现无脉动的平稳旋转。而六边形轨迹意味着磁链幅值在 $ \frac{2}{3}U_{dc} $ 与 $ \frac{2}{3}U_{dc} \cdot \cos(30^\circ) $ 之间周期性波动，这种幅值不恒定性直接导致了电磁转矩的脉动（Torque Ripple），表现为电机振动、噪声增大及效率下降。SVPWM的工程价值，正在于它提供了一套严谨、可计算的数学框架，用以解决“如何用有限的8个离散电压矢量，去逼近无限多个连续旋转电压矢量”这一根本性问题，从而将六边形的磁链轨迹“柔化”为高保真度的圆形轨迹。

2. 电压空间矢量的数学建模与几何分析

2.1 三相桥臂开关状态与线电压的映射

为建立清晰的物理图像，需从逆变器拓扑出发，严格推导开关状态与输出电压的关系。定义三相桥臂的开关函数 $ S_a, S_b, S_c $，其中 $ S_x = 1 $ 表示上桥臂导通、下桥臂关断；$ S_x = 0 $ 表示上桥臂关断、下桥臂导通。则三相相对于直流母线中点 $ O $ 的相电压为：
$$
U_{aO} = S_a \cdot U_{dc}, \quad U_{bO} = S_b \cdot U_{dc}, \quad U_{cO} = S_c \cdot U_{dc}
$$
三相线电压 $ U_{ab}, U_{bc}, U_{ca} $ 可由相电压差得到：
$$
U_{ab} = U_{aO} - U_{bO}, \quad U_{bc} = U_{bO} - U_{cO}, \quad U_{ca} = U_{cO} - U_{aO}
$$

以开关状态 $ [S_a, S_b, S_c] = [1, 0, 0] $ 为例，此时 $ U_{aO} = U_{dc}, U_{bO} = U_{cO} = 0 $，故 $ U_{ab} = U_{dc}, U_{bc} = 0, U_{ca} = -U_{dc} $。根据Clarke变换（αβ变换）的定义，电压空间矢量 $ \mathbf{U} $ 在αβ坐标系中的分量为：
$$
U_\alpha = \frac{2}{3}(U_{aO} - \frac{1}{2}U_{bO} - \frac{1}{2}U_{cO}), \quad U_\beta = \frac{2}{3}(\frac{\sqrt{3}}{2}U_{bO} - \frac{\sqrt{3}}{2}U_{cO})
$$
将 $ [1,0,0] $ 代入，得 $ U_\alpha = \frac{2}{3}U_{dc}, U_\beta = 0 $，即 $ \mathbf{U} 1 = \frac{2}{3}U {dc} \angle 0^\circ $。

对全部8种开关状态进行相同计算，可得完整的电压矢量表：

开关状态 $[S_a,S_b,S_c]$	矢量名称	幅值	相位角	物理含义
$[0,0,0]$	$ \mathbf{U}_0 $	$ 0 $	—	全上桥臂关断，三相悬空
$[1,0,0]$	$ \mathbf{U}_1 $	$ \frac{2}{3}U_{dc} $	$ 0^\circ $	A相接正，B/C相接负
$[1,1,0]$	$ \mathbf{U}_2 $	$ \frac{2}{3}U_{dc} $	$ 60^\circ $	A/B相接正，C相接负
$[0,1,0]$	$ \mathbf{U}_3 $	$ \frac{2}{3}U_{dc} $	$ 120^\circ $	B相接正，A/C相接负
$[0,1,1]$	$ \mathbf{U}_4 $	$ \frac{2}{3}U_{dc} $	$ 180^\circ $	B/C相接正，A相接负
$[0,0,1]$	$ \mathbf{U}_5 $	$ \frac{2}{3}U_{dc} $	$ 240^\circ $	C相接正，A/B相接负
$[1,0,1]$	$ \mathbf{U}_6 $	$ \frac{2}{3}U_{dc} $	$ 300^\circ $	A/C相接正，B相接负
$[1,1,1]$	$ \mathbf{U}_7 $	$ 0 $	—	全上桥臂导通，三相悬空

2.2 六扇区划分与磁链轨迹的几何直观

将上述6个非零矢量 $ \mathbf{U}_1 $ 至 $ \mathbf{U}_6 $ 绘制在复平面上，它们首尾相连，恰好将整个平面划分为6个互不重叠的扇形区域（Sector），每个扇区夹角为 $ 60^\circ $。扇区编号 $ N $ 通常按逆时针方向从1开始排列，其中扇区1覆盖 $ 0^\circ \leq \theta < 60^\circ $，扇区2覆盖 $ 60^\circ \leq \theta < 120^\circ $，依此类推。

此划分具有深刻的物理意义。对于任意一个位于扇区 $ N $ 内的期望电压矢量 $ \mathbf{V} {ref} $，根据向量空间的线性组合原理，它必然可以且仅可以由该扇区相邻的两个非零矢量 $ \mathbf{U}_x $ 和 $ \mathbf{U}_y $（例如扇区1中为 $ \mathbf{U}_1 $ 和 $ \mathbf{U}_2 $）以及零矢量 $ \mathbf{U}_0 $ 或 $ \mathbf{U}_7 $ 线性合成。这种合成并非简单的算术加法，而是在时间维度上的加权平均。设 $ \mathbf{V} {ref} $ 的幅值为 $ V_{ref} $，相位为 $ \theta $，PWM开关周期为 $ T_s $，则有：
$$
\mathbf{V} {ref} \cdot T_s = \mathbf{U}_x \cdot T_x + \mathbf{U}_y \cdot T_y + \mathbf{U}_0 \cdot T_0
$$
其中 $ T_x + T_y + T_0 = T_s $。该方程在二维平面上提供了两个独立的标量方程（α轴与β轴分量），恰好可以求解出 $ T_x $、$ T_y $ 和 $ T_0 $ 这三个未知数。零矢量 $ \mathbf{U}_0 $ 或 $ \mathbf{U}_7 $ 的引入，不仅提供了求解自由度，更重要的是，它允许我们通过调节 $ T_0 $ 来控制合成矢量的有效作用时间，从而实现对 $ \mathbf{V} {ref} $ 幅值的连续调节，这是SVPWM实现调压功能的核心机制。

3. 扇区判断与作用时间计算的工程实现

3.1 基于 $ U_\alpha $、$ U_\beta $ 的原始判断法

扇区判断是SVPWM算法的第一步，其目标是确定 $ \mathbf{V} {ref} $ 所在的60°区间。最直接的方法是利用三角函数关系。以扇区1为例，其边界条件为：
- $ U \beta > 0 $ （位于上半平面）
- $ U_\alpha > 0 $ （位于右半平面）
- $ \tan\theta = \frac{U_\beta}{U_\alpha} < \tan(60^\circ) = \sqrt{3} $，即 $ U_\beta < \sqrt{3} U_\alpha $

同理，可推导出其余扇区的判据。然而，这种基于浮点除法和三角函数的判断在资源受限的MCU上计算开销较大，且易受数值精度影响。因此，工程实践中普遍采用一种巧妙的代数变换，定义三个中间变量 $ X, Y, Z $：
$$
X = U_\beta, \quad Y = \frac{\sqrt{3}}{2}U_\alpha - \frac{1}{2}U_\beta, \quad Z = -\frac{\sqrt{3}}{2}U_\alpha - \frac{1}{2}U_\beta
$$
这三个变量的物理意义是：它们分别代表了 $ \mathbf{V}_{ref} $ 在三个互成120°的坐标轴（即对应于 $ \mathbf{U}_2, \mathbf{U}_4, \mathbf{U}_6 $ 方向）上的投影。扇区判断逻辑可简化为对 $ X, Y, Z $ 符号的组合查询：

扇区 $ N $	$ X $	$ Y $	$ Z $	判据说明
1	>0	>0	<0	$ \mathbf{V}_{ref} $ 在 $ \mathbf{U}_1 $ 与 $ \mathbf{U}_2 $ 之间
2	<0	>0	<0	$ \mathbf{V}_{ref} $ 在 $ \mathbf{U}_2 $ 与 $ \mathbf{U}_3 $ 之间
3	<0	>0	>0	$ \mathbf{V}_{ref} $ 在 $ \mathbf{U}_3 $ 与 $ \mathbf{U}_4 $ 之间
4	<0	<0	>0	$ \mathbf{V}_{ref} $ 在 $ \mathbf{U}_4 $ 与 $ \mathbf{U}_5 $ 之间
5	>0	<0	>0	$ \mathbf{V}_{ref} $ 在 $ \mathbf{U}_5 $ 与 $ \mathbf{U}_6 $ 之间
6	>0	<0	<0	$ \mathbf{V}_{ref} $ 在 $ \mathbf{U}_6 $ 与 $ \mathbf{U}_1 $ 之间

此方法仅需三次比较操作，无需乘除与三角函数，极大提升了实时性与鲁棒性。

3.2 作用时间的解析推导与通用公式

在确定扇区后，即可计算基础矢量的作用时间。以扇区1为例，$ \mathbf{V} {ref} $ 由 $ \mathbf{U}_1 $ 和 $ \mathbf{U}_2 $ 合成。将 $ \mathbf{V} {ref} $ 分解到 $ \mathbf{U} 1 $ 和 $ \mathbf{U}_2 $ 方向上，利用几何关系（见图1）：
- $ \mathbf{U}_1 $ 与 $ \mathbf{U}_2 $ 夹角为 $ 60^\circ $。
- $ \mathbf{V} {ref} $ 在垂直于 $ \mathbf{U} 1 $ 方向的分量为 $ V {ref} \sin\theta $，该分量完全由 $ \mathbf{U} 2 $ 提供，故 $ V {ref} \sin\theta \cdot T_s = |\mathbf{U} 2| \cdot T_2 \cdot \sin(60^\circ) $。
- 代入 $ |\mathbf{U}_2| = \frac{2}{3}U {dc} $，得 $ T_2 = \frac{3}{2} \cdot \frac{U_\beta}{U_{dc}} \cdot T_s $。

同理，在 $ \mathbf{U} 1 $ 方向的分量满足：
$$
V {ref} \cos\theta \cdot T_s = |\mathbf{U} 1| \cdot T_1 + |\mathbf{U}_2| \cdot T_2 \cdot \cos(60^\circ)
$$
解得 $ T_1 = \frac{3}{2} \cdot \frac{U \alpha - \frac{1}{\sqrt{3}}U_\beta}{U_{dc}} \cdot T_s $。

将 $ T_1 $、$ T_2 $ 用中间变量 $ X, Y, Z $ 表示，可得到一套高度统一的公式。观察发现，对于任意扇区，其 $ T_x $ 和 $ T_y $ 总是等于 $ X, Y, Z $ 中某两个变量的绝对值与 $ T_s / U_{dc} $ 的乘积。例如，在扇区1中，$ T_1 = \frac{T_s}{U_{dc}} \cdot |Y| $，$ T_2 = \frac{T_s}{U_{dc}} \cdot |X| $。这种模式在整个六扇区中保持一致，使得代码实现可以高度模块化，只需一个查表（或switch-case）即可获取当前扇区对应的变量映射关系。

3.3 过调制处理与最大输出能力

SVPWM的理论最大输出电压幅值并非 $ U_{dc} $，而是内切于六边形的圆的半径，即 $ V_{ref}^{max} = \frac{2}{3}U_{dc} $。当控制器计算出的 $ V_{ref} $ 超过此值时，即进入过调制（Over-modulation）区域。此时，若仍强行按线性关系计算 $ T_1 $、$ T_2 $，会导致 $ T_1 + T_2 > T_s $，违反时间守恒定律，算法失效。

工程上，过调制处理通常分为两个阶段：
1. 线性调制区 ($ V_{ref} \leq \frac{2}{3}U_{dc} $)：按前述公式正常计算。
2. 六拍波区 ($ \frac{2}{3}U_{dc} < V_{ref} \leq U_{dc} $)：将 $ \mathbf{V}_{ref} $ “钳位”至六边形边界，此时输出波形退化为方波，但能提供更高的基波电压。具体实现是将 $ T_1 $、$ T_2 $ 按比例缩放，使 $ T_1 + T_2 = T_s $，并用零矢量填充剩余时间。更高级的方案会采用特定的过调制算法（如七段式过调制），在保证输出电压连续性的前提下，最大化利用直流母线电压。

4. 七段式SVPWM波形生成与定时器配置

4.1 七段式发波策略的原理与优势

SVPWM的波形生成，本质上是将一个PWM周期 $ T_s $ 划分为若干时间段，并在每个时间段内设置逆变器的开关状态。七段式（Seven-Segment）发波是一种被广泛采用的优化策略，其核心思想是：在一个 $ T_s $ 周期内，安排7个时间片段，依次使用 $ \mathbf{U}_0 $、$ \mathbf{U}_x $、$ \mathbf{U}_y $、$ \mathbf{U}_7 $、$ \mathbf{U}_y $、$ \mathbf{U}_x $、$ \mathbf{U}_0 $ 这7个矢量。其中，零矢量 $ \mathbf{U}_0 $ 和 $ \mathbf{U}_7 $ 的作用时间各为 $ T_0/2 $，基础矢量 $ \mathbf{U}_x $ 和 $ \mathbf{U}_y $ 的作用时间各为 $ T_x/2 $ 和 $ T_y/2 $。

这种对称结构带来了两大显著优势：
- 开关损耗最小化 ：由于每次切换只改变一个桥臂的状态（即只有一个比特翻转），整个周期内开关管的总动作次数仅为6次，远低于五段式或非对称式的12次。这直接降低了开关损耗与EMI噪声。
- 谐波特性最优 ：对称的波形结构使得输出电压的偶次谐波被自然抵消，主要谐波集中在 $ 6k \pm 1 $（k为整数）次，即5、7、11、13次等，这些谐波易于通过电机绕组电感滤除，从而获得更纯净的基波电流。

4.2 STM32高级定时器TIM1的中央对齐模式配置

在STM32平台上，利用高级定时器（如TIM1、TIM8）的中央对齐（Center-Aligned）模式是实现七段式SVPWM的理想选择。该模式下，计数器 $ CNT $ 并非单向递增，而是先从0计数至自动重装载寄存器 $ ARR $，再反向递减回0，构成一个完整的PWM周期 $ T_s = 2 \times ARR \times T_{clk} $。此模式天然支持对称的七段波形。

关键寄存器配置如下：
- ARR ：设置为 $ T_s/2 $ 对应的计数值。例如，若系统主频为168MHz，PWM频率为10kHz（$ T_s = 100\mu s $），则 $ ARR = (168 \times 10^6 \times 100 \times 10^{-6}) / 2 = 8400 $。
- CCRx ：捕获/比较寄存器，用于设定PWM通道的占空比。其值决定了 $ CNT $ 在上升沿和下降沿的切换点。
- CCMRx ：捕获/比较模式寄存器，需配置为PWM模式1或模式2。
- CCER ：捕获/比较使能寄存器，使能相应通道。

4.3 PWM模式1与模式2的切换点计算

SVPWM的最终输出是三相PWM信号，每相的高电平持续时间（即有效作用时间）由 $ T_x $、$ T_y $ 和 $ T_0 $ 共同决定。以A相为例，其在一个 $ T_s $ 周期内的导通时刻（Switching Point）$ T_{p1} $、$ T_{p2} $、$ T_{p3} $ 可由矢量作用顺序推导得出。对于扇区1的七段式序列 $ [\mathbf{U}_0, \mathbf{U}_1, \mathbf{U}_2, \mathbf{U}_7, \mathbf{U}_2, \mathbf{U}_1, \mathbf{U}_0] $，A相的开关状态变化如下：
- $ \mathbf{U}_0 $：$ S_a=0 $，A相下管导通，输出低电平。
- $ \mathbf{U}_1 $：$ S_a=1 $，A相上管导通，输出高电平。
- $ \mathbf{U}_2 $：$ S_a=1 $，A相上管导通，输出高电平。
- $ \mathbf{U}_7 $：$ S_a=1 $，A相上管导通，输出高电平。
- 后续同理。

因此，A相的第一次上升沿发生在 $ T_{p1} = T_0/2 $ 时刻，第一次下降沿发生在 $ T_{p2} = T_0/2 + T_1/2 + T_2/2 + T_7/2 = T_0/2 + T_1/2 + T_2/2 + T_0/2 = T_0 + T_1/2 + T_2/2 $ 时刻。最终，A相的高电平持续时间为 $ T_{onA} = T_{p2} - T_{p1} = T_0/2 + T_1/2 + T_2/2 $。

在PWM模式2（Active High on Compare Match）下，当 $ CNT $ 从0开始递增，到达 $ CCRA $ 时，通道输出高电平；当 $ CNT $ 从 $ ARR $ 开始递减，再次到达 $ CCRA $ 时，通道输出低电平。因此，要使A相在 $ T_{p1} $ 上升、$ T_{p2} $ 下降，需设置 $ CCRA = T_{p1} $。

在PWM模式1（Active High on Compare Match, Inverted）下，当 $ CNT $ 小于 $ CCRA $ 时，通道输出高电平；否则输出低电平。此时，$ CCRA $ 应设置为A相高电平的“中心点”，即 $ CCRA = T_{onA}/2 $。

两种模式在硬件上均能正确复现所需的波形，选择取决于开发者习惯与已有代码库的兼容性。模式2的逻辑更直观，与“切换点”的概念直接对应；模式1则与某些SDK（如ST MC SDK）的实现方式一致。

5. SVPWM算法的工程化封装与实战经验

5.1 算法流程的模块化封装

一个工业级的SVPWM函数，其接口应简洁、鲁棒且可移植。典型的C语言函数签名如下：

void SVPWM_Generate(uint16_t *pDutyA, uint16_t *pDutyB, uint16_t *pDutyC,
                     float U_alpha, float U_beta, float U_dc, uint16_t T_s);

其内部实现可划分为清晰的模块：
1. 输入预处理 ：对 $ U_\alpha $、$ U_\beta $ 进行限幅，防止过调制导致计算溢出。
2. 扇区判断 ：计算 $ X, Y, Z $，通过查表或逻辑运算确定扇区号 $ N $。
3. 时间计算 ：根据扇区号 $ N $，从预定义的映射表中取出对应的变量索引，计算 $ T_x $、$ T_y $、$ T_0 $。
4. 切换点计算 ：依据七段式规则，计算出三相的六个切换点 $ T_{p1} $ 至 $ T_{p6} $。
5. 输出映射 ：将切换点转换为定时器比较寄存器 $ CCRx $ 的值，并写入指针参数。

这种封装方式将复杂的数学计算与底层硬件细节完全解耦，上层FOC算法只需提供 $ U_\alpha $、$ U_\beta $，即可获得可直接载入定时器的占空比数据，极大地提升了代码的可维护性与复用性。

5.2 实战中的关键调试技巧与常见陷阱

在实际项目调试中，SVPWM是故障率较高的模块之一，以下几点经验至关重要：
- 死区时间（Dead Time）的插入 ：SVPWM计算出的切换点是理想值。在将其写入定时器之前，必须通过 HAL_TIMEx_ConfigDeadTime() 函数为每个通道配置死区。典型值为1~2μs。 未配置死区是导致IGBT直通炸机的最主要原因 。死区会轻微畸变波形，因此在计算 $ T_x $、$ T_y $ 时，应预留足够的裕量（例如，将 $ T_s $ 设为略小于理论值）。
- 坐标系一致性检查 ：务必确认Clarke变换的符号约定与SVPWM模块一致。常见的错误是 $ U_\beta $ 的符号取反，这会导致扇区判断完全错误，电机剧烈抖动。一个快速验证方法是：给FOC环路注入一个微小的 $ I_q^ $ 给定，观察电机是否平稳旋转；若反转或抖动，则首先检查 $ U_\alpha $、$ U_\beta $ 的符号。
- 过调制的平滑过渡 ：在线性调制区与六拍波区的交界处，若处理不当会产生明显的扭矩冲击。推荐的做法是：当 $ V_{ref} $ 接近 $ \frac{2}{3}U_{dc} $ 时，启动一个缓慢的软切换过程，逐步将算法从线性插值过渡到六拍波逻辑，而非硬切换。
- ADC采样与PWM更新的同步 *：为了获得最佳的电流环性能，ADC对相电流的采样时刻应严格对齐在PWM周期的中点（即 $ CNT = ARR $ 时刻），此时上下桥臂的占空比均为50%，电流纹波最小。这需要配置ADC的触发源为定时器的更新事件（Update Event）。

我在一个3kW伺服驱动器项目中，曾因忽略了死区时间的温度漂移特性，在高温环境下连续烧毁了三块驱动板。后来改为在固件中加入温度传感器读数，并根据温度查表动态调整死区时间，问题才彻底解决。这提醒我们，SVPWM不仅是数学游戏，更是与物理世界紧密耦合的工程实践。