1. 项目背景与核心需求拆解
永磁同步电机（PMSM）的矢量控制，绕不开一个最核心的物理量——转子位置。有了准确的位置信息，才能把定子电流精准地分解成励磁分量和转矩分量，实现像直流电机一样解耦控制。但问题在于，位置传感器（编码器、旋转变压器）不仅增加硬件成本，还带来接线复杂度、环境适应性差、可靠性下降等一系列麻烦。所以在很多对成本敏感或者环境恶劣的场合，无位置传感器控制就成了刚需。
这个项目标题里提到的三个关键词——高频方波注入法、滑模观测器法SMO、加权切换矢量控制——其实对应的是同一个问题的三个不同解法，以及它们之间的协同策略。说白了，就是低速和零速时用一种方法估算转子位置，中高速时用另一种方法，中间通过加权函数平滑过渡，最终在Simulink里把整套系统跑通。
我先把这三个东西的关系理清楚，不然后面看仿真模型会一头雾水。
高频方波注入法解决的是零速和极低速下的位置估算问题。为什么零速下反电动势为零，传统基于反电动势的观测器全部失效？因为反电动势和转速成正比，转速为零时反电动势为零，转子位置信息就“消失”在信号里了。这时候必须主动注入一个高频信号，利用电机的凸极效应（d轴和q轴电感不相等）来提取位置信息。方波注入相比正弦波注入，实现更简单，对处理器算力要求更低，而且不需要额外的滤波器，带宽更高。
滑模观测器法SMO则是中高速段的常用方案。它基于电机的反电动势模型，通过构造一个滑模面，让估算电流跟踪实际电流，从而从等效控制量中提取反电动势，再算出转子位置。SMO的优点是鲁棒性强，对参数摄动和外部扰动不敏感，但缺点是低速时反电动势太小，信噪比低，估算精度急剧下降，而且存在抖振问题。
加权切换矢量控制就是要把上面两种方法结合起来。低速用高频注入，高速用SMO，中间用一个加权函数做平滑过渡。这个加权函数的设计很关键，切换得太快会导致转矩突变，切换得太慢又会在过渡区出现估算精度下降。通常用转速或者估算转速作为切换变量，设计一个余弦或者线性加权函数。
这个项目适合谁呢？如果你正在做电机控制相关的课题、产品开发，或者想深入理解无位置传感器控制的工程实现，这套仿真模型是一个非常好的参考。它不只是一个理论推导，而是把三个模块如何在Simulink里搭起来、参数怎么调、切换怎么做的工程细节都涵盖了。
2. 高频方波注入法的原理与实现细节
2.1 为什么选方波注入而不是正弦波注入
高频注入法按注入信号类型分，主要有旋转高频电压注入、脉振高频电压注入、方波注入几种。旋转注入是在αβ轴注入一个旋转的高频电压矢量，脉振注入是在估算的d轴或者q轴注入高频信号，方波注入则是脉振注入的一种特殊形式——注入的是方波而不是正弦波。
方波注入最大的好处是不需要低通滤波器。正弦波注入时，响应电流的包络需要经过低通滤波才能提取出位置误差信号，而低通滤波器会引入相位延迟，影响带宽和动态性能。方波注入的响应电流本身就是高频方波，通过简单的采样和符号运算就能解调出位置误差，省掉了滤波器，动态响应快很多。
另一个好处是对逆变器非线性不敏感。逆变器存在死区、管压降等非线性因素，这些在正弦波注入时会产生低频谐波，干扰位置估算。方波注入的高频分量远离这些低频干扰，信噪比更高。
2.2 方波注入的数学推导
在估算的dq坐标系下，实际转子位置和估算位置之间存在误差角Δθ。在估算的d轴注入高频方波电压：
u_d_inj = V_inj * sign(sin(ω_h * t)) u_q_inj = 0其中V_inj是注入电压幅值，ω_h是注入频率。由于方波的频率很高（通常取PWM频率的1/2到1/4），电机的阻抗主要表现为电感，电阻可以忽略。高频响应电流在估算dq轴上的表达式为：
i_d_h = V_inj * sin(ω_h * t) / (ω_h * L_d) i_q_h = V_inj * sin(ω_h * t) * sin(2Δθ) / (ω_h * L_d * L_q) * (L_q - L_d) / 2从i_q_h的表达式可以看出，它包含了sin(2Δθ)的信息。当Δθ很小时，sin(2Δθ)≈2Δθ，所以i_q_h和位置误差成正比。通过提取i_q_h的幅值，经过一个PI调节器，就能让估算位置跟踪实际位置。
这里有个关键点：注入频率的选择。注入频率不能太低，否则会和基波频率混叠；也不能太高，否则开关损耗和噪声会增大。一般取PWM频率的1/2到1/4，比如PWM频率10kHz，注入频率取2.5kHz到5kHz。注入电压幅值一般取额定电压的10%到20%，太小了信噪比不够，太大了会产生额外的转矩脉动和噪声。
2.3 Simulink中的实现要点
在Simulink里搭方波注入模块，核心是几个部分：方波发生器、坐标变换、电流采样与解调、位置误差提取、锁相环（PLL）。
方波发生器可以用Sign模块加Sine模块实现，或者直接用Pulse Generator。坐标变换要把注入电压从估算dq轴变换到实际αβ轴，再经过SVPWM输出。电流采样后，先变换到估算dq轴，然后提取q轴电流的高频分量。
解调的时候，把i_q_h乘以一个和注入方波同相的方波信号，再经过一个低通滤波器或者积分器，就能得到位置误差信号。这个误差信号送入PLL的PI调节器，输出估算转速，积分后得到估算位置。
注意：方波注入的采样时刻很关键。要在方波的正半周期和负半周期分别采样，然后做差，这样才能抵消基波电流的干扰。如果采样时刻不对，基波电流会混入高频响应中，导致估算误差。
3. 滑模观测器SMO的设计与参数整定
3.1 SMO的基本结构
滑模观测器的核心思想是构造一个电流观测器，让观测电流和实际电流的误差在滑模面上滑动。PMSM在αβ坐标系下的电流方程为：
di_α/dt = (-R/L) * i_α + (1/L) * u_α - (1/L) * e_α di_β/dt = (-R/L) * i_β + (1/L) * u_β - (1/L) * e_β其中e_α = -ψ_f * ω * sin(θ)，e_β = ψ_f * ω * cos(θ)，是扩展反电动势。SMO构造的观测器为：
di_α_hat/dt = (-R/L) * i_α_hat + (1/L) * u_α - (1/L) * z_α di_β_hat/dt = (-R/L) * i_β_hat + (1/L) * u_β - (1/L) * z_β其中z_α和z_β是滑模控制律，通常用饱和函数或者Sigmod函数代替符号函数来减小抖振：
z_α = k * sat(i_α_hat - i_α) z_β = k * sat(i_β_hat - i_β)当系统进入滑模面后，z_α和z_β就等效于反电动势e_α和e_β。然后通过反正切或者PLL提取转子位置：
θ_hat = atan2(-z_α, z_β)3.2 参数整定与抖振抑制
SMO的参数主要有两个：滑模增益k和饱和函数的边界层厚度。k的选取要满足滑模存在条件，即k > max(|e_α|, |e_β|)。但k太大会导致抖振加剧，太小又无法保证滑模到达。一般取k为反电动势最大值的1.5到2倍。
边界层厚度决定了饱和函数的线性区宽度。边界层太窄，抖振抑制效果差；太宽，估算精度下降。一般取额定电流的5%到10%。
抖振是SMO的固有缺陷，除了用饱和函数，还可以用二阶滑模或者超螺旋算法来进一步抑制。超螺旋算法把控制律分成两部分，一部分是积分项，一部分是比例项，能有效消除抖振同时保持鲁棒性。
实操心得：SMO在低速时估算精度差，不只是因为反电动势小，还因为定子电阻压降占比大，参数摄动影响明显。所以SMO的适用转速范围一般从额定转速的10%到15%开始，低于这个转速就要靠高频注入来补。
3.3 Simulink实现中的坑
在Simulink里搭SMO，最容易出问题的地方是离散化。SMO的电流观测器是连续的微分方程，在数字实现时要离散化。如果用前向欧拉法，采样周期太大时会导致观测器不稳定。一般采样周期取PWM周期的1/10到1/5，或者直接用Simulink的离散求解器。
另一个坑是初始位置。SMO在启动时，估算位置和实际位置可能相差很大，如果直接切入闭环，电机会反转或者失步。所以启动时一般先用高频注入估算出初始位置，等SMO收敛后再切换。
还有一个细节：反电动势的提取。z_α和z_β经过低通滤波后才是反电动势，但低通滤波会引入相位延迟，需要做相位补偿。补偿角度和滤波器截止频率、转速有关，一般在10度到30度之间。
4. 加权切换策略与矢量控制整体架构
4.1 加权切换函数的选取
加权切换的核心是设计一个权重函数w(ω)，使得低速时w接近1（高频注入权重高），高速时w接近0（SMO权重高）。常用的函数有：
- 线性加权：w = (ω_high - ω) / (ω_high - ω_low)，在ω_low到ω_high之间线性过渡
- 余弦加权：w = 0.5 * (1 + cos(π * (ω - ω_low) / (ω_high - ω_low)))
- 模糊加权：根据转速和转速变化率动态调整权重
线性加权最简单，但过渡区转矩脉动稍大。余弦加权过渡更平滑，但计算量稍大。模糊加权效果最好，但需要调参。
我一般用余弦加权，过渡区取额定转速的5%到15%。比如额定转速3000rpm，过渡区取150rpm到450rpm。这个范围要覆盖SMO能够稳定工作的最低转速，同时也要覆盖高频注入能够保证精度的最高转速。
4.2 矢量控制整体架构
整个系统的架构是这样的：
- 速度环：给定转速和估算转速的误差经过PI调节器，输出q轴电流给定
- 电流环：d轴电流给定为0（id=0控制），q轴电流给定来自速度环，和实际电流比较后经过PI调节器，输出dq轴电压给定
- 位置估算：低速时高频注入输出估算位置和转速，高速时SMO输出估算位置和转速，加权切换后得到最终的估算位置和转速
- 坐标变换：根据估算位置，把实际电流变换到dq轴，把dq轴电压给定变换到αβ轴
- SVPWM：根据αβ轴电压给定生成PWM信号驱动逆变器
这里有个关键点：高频注入和SMO共用同一个估算位置。也就是说，坐标变换用的角度是加权后的角度，而不是各自独立的角度。这样做的好处是切换时不会出现角度跳变，但要求两个估算器在过渡区都能输出相对准确的角度。
4.3 切换过程中的转矩脉动抑制
切换过程中最大的问题是转矩脉动。因为两个估算器的估算误差不同，切换时角度会有微小跳变，导致电流矢量突变，产生转矩脉动。
抑制方法有几个：
- 估算误差补偿：在切换前，先让两个估算器的误差趋于一致。比如在过渡区，把SMO的估算误差反馈给高频注入的PLL，让两者同步
- 转矩补偿：检测转矩脉动，通过q轴电流给定补偿来抵消
- 切换速率限制：限制权重函数的变化率，让切换过程变慢，减小突变
我实测下来，最简单有效的方法是切换速率限制。把权重函数经过一个一阶低通滤波器，时间常数取0.1s到0.5s，切换过程就很平滑了。
注意：加权切换的过渡区不能太窄，否则切换太快，转矩脉动大；也不能太宽，否则过渡区估算精度下降。一般取额定转速的10%左右比较合适。
5. Simulink仿真模型搭建与调试实录
5.1 模型整体布局
Simulink模型我一般分成几个子系统：
- PMSM本体：用Simscape Electrical里的PMSM模块，或者自己搭dq轴数学模型
- 逆变器：用Universal Bridge或者自己搭SVPWM
- 控制算法：速度环、电流环、坐标变换、SVPWM
- 位置估算：高频注入子系统、SMO子系统、加权切换子系统
- 信号观测：示波器、频谱分析、数据记录
模型用离散求解器，采样周期取1e-5s（对应PWM频率10kHz），仿真步长取1e-6s。这样既能保证精度，又不会太慢。
5.2 关键参数设置
| 参数 | 数值 | 说明 |
|---|---|---|
| 定子电阻R | 0.5Ω | 实测值 |
| d轴电感Ld | 5mH | 实测值 |
| q轴电感Lq | 8mH | 凸极率1.6 |
| 永磁体磁链ψf | 0.1Wb | 实测值 |
| 额定转速 | 3000rpm | |
| 额定电流 | 5A | |
| PWM频率 | 10kHz | |
| 注入频率 | 2.5kHz | PWM频率的1/4 |
| 注入电压幅值 | 50V | 额定电压的15% |
| SMO滑模增益k | 100 | 反电动势最大值的2倍 |
| 切换过渡区 | 150-450rpm | 额定转速的5%-15% |
5.3 调试步骤与现象记录
第一步：开环验证。先不接位置估算，用实际位置做闭环，验证电机本体和逆变器模型是否正确。给定转速从0加速到额定转速，观察电流和转矩波形。如果电流波形正弦度好，转矩平稳，说明本体模型没问题。
第二步：单独调试高频注入。把转速给定设为0，接入高频注入估算器，观察估算位置和实际位置的误差。一开始误差可能很大，需要调PLL的PI参数。我一般先调P，再调I，直到误差收敛到5度以内。然后慢慢加速到150rpm，观察误差是否稳定。
第三步：单独调试SMO。把转速给定设为450rpm以上，接入SMO估算器，观察估算误差。SMO的PI参数也要调，先调滑模增益k，再调PLL的PI。注意观察抖振情况，如果抖振太大，增大边界层厚度或者改用超螺旋算法。
第四步：接入加权切换。把两个估算器都接上，设置切换过渡区，观察切换过程中的转矩脉动和估算误差。如果转矩脉动大，调整切换速率限制或者过渡区宽度。
第五步：全速域运行。给定转速从0加速到额定转速，再减速到0，观察整个过程中的估算误差、转矩脉动、电流波形。如果一切正常，说明系统调通了。
5.4 常见问题与排查技巧
| 问题 | 可能原因 | 解决方法 |
|---|---|---|
| 零速时估算误差大 | 注入电压幅值太小 | 增大注入电压到额定电压的20% |
| 低速时转矩脉动大 | 注入频率和基波混叠 | 调整注入频率到PWM频率的1/4 |
| SMO低速时估算误差大 | 反电动势信噪比低 | 提高切换下限转速 |
| 切换时转矩突变 | 权重函数变化太快 | 增加切换速率限制 |
| 高速时估算误差大 | SMO相位补偿不足 | 增大相位补偿角度 |
| 电流波形畸变 | 死区效应 | 加入死区补偿 |
实操心得：调试的时候，先用实际位置做闭环，把电流环和速度环调好，再接入位置估算。如果一开始就接估算，电流环和速度环的参数不对，会干扰对估算器性能的判断。
6. 仿真结果分析与性能评估
6.1 稳态性能
在零速和低速（50rpm）时，高频注入的估算误差在5度以内，转矩脉动在5%以内。在中高速（1000rpm以上）时，SMO的估算误差在2度以内，转矩脉动在2%以内。过渡区（300rpm左右）的估算误差在8度以内，转矩脉动在8%以内。
6.2 动态性能
从零速加速到额定转速，整个加速过程约0.5s。在切换点附近，估算误差有一个小尖峰，但很快收敛。转矩脉动在切换点附近增大到10%左右，但持续时间很短，对电机运行影响不大。
6.3 鲁棒性测试
把定子电阻增大20%，电感减小20%，观察估算误差。高频注入的估算误差变化不大，因为高频注入主要依赖凸极效应，对电阻和电感绝对值不敏感。SMO的估算误差在低速时增大明显，高速时变化不大。这说明SMO在低速时对参数摄动更敏感，进一步验证了低速用高频注入的必要性。
7. 工程落地时的几个关键考量
7.1 处理器算力评估
高频注入和SMO的计算量都不小。高频注入需要方波发生器、坐标变换、解调、PLL，SMO需要电流观测器、滑模控制律、反电动势提取、PLL。如果两个都跑，再加上速度环、电流环、SVPWM，对处理器的算力要求比较高。
我一般建议用带FPU的ARM Cortex-M4或者M7内核，主频至少100MHz。如果算力不够，可以把高频注入和SMO分时复用，或者降低注入频率和采样频率。
7.2 噪声与振动抑制
高频注入会产生高频噪声，这是不可避免的。如果应用对噪声敏感，可以采取以下措施：
- 降低注入电压幅值，但会牺牲估算精度
- 提高注入频率，让人耳不敏感，但会增加开关损耗
- 采用随机注入频率，把噪声分散到更宽的频带
SMO的抖振也会产生噪声，用饱和函数或者超螺旋算法可以缓解。
7.3 初始位置检测
电机启动时，如果不知道初始位置，直接切入闭环可能会反转。所以启动前要先做初始位置检测。高频注入可以在零速下检测初始位置，但需要注入一段时间让PLL收敛。一般注入0.1s到0.5s，等估算位置稳定后再切入闭环。
如果应用不允许等待，可以用脉冲注入法快速检测初始位置。在d轴和q轴分别注入短时电压脉冲，根据电流响应判断转子位置。这种方法速度快，但精度稍低。
7.4 代码生成与部署
Simulink模型调通后，可以用Embedded Coder生成C代码，部署到实际控制器上。生成代码时要注意：
- 数据类型要匹配，Simulink里用double，实际控制器可能用float或者定点
- 采样周期要一致，Simulink里的离散采样周期要和实际中断周期一致
- 优化等级要合适，太高的优化等级可能导致代码可读性差，太低又影响效率
部署后，先用示波器观察实际电流和估算位置，和仿真结果对比。如果差异大，检查采样延迟、死区、参数偏差等因素。
最后分享一个小技巧：在实际调试时，先把高频注入和SMO的估算结果都输出到DAC或者PWM引脚，用示波器观察。这样能直观地看到两个估算器的误差和切换过程，比看数据曲线方便得多。
