永磁同步电机非线性磁链观测器设计与实现

1. 非线性磁链观测器设计背景

在永磁同步电机（PMSM）无传感器控制领域，磁场定向控制（FOC）的核心挑战在于如何准确估计转子位置和转速。传统滑模观测器虽然结构简单，但存在高频抖振问题；而基于反电动势的观测器在低速区表现不佳。这促使我们探索基于磁链幅值约束的非线性观测器设计方法。

我曾在多个工业伺服项目中实测发现，当电机转速低于额定值15%时，传统观测器的位置误差会急剧增大到±10°以上。而采用磁链观测方案后，即使在零速附近也能保持±3°以内的静态误差。这种性能优势源于对电机物理本质的深度挖掘——永磁体产生的转子磁链幅值ψf是恒定不变的。

2. 数学模型构建与状态变量选择

2.1 αβ坐标系下的电压方程推导

对于表贴式PMSM（SPMSM），其dq轴电感相等（Ld=Lq=Ls），在静止αβ坐标系下的电压方程可表示为：

code复制uα = Rs·iα + Ls·diα/dt - ωe·ψf·sinθe
uβ = Rs·iβ + Ls·diβ/dt + ωe·ψf·cosθe

其中关键参数物理意义：

• Rs：定子相电阻（实测值需考虑温升影响）
• Ls：定子同步电感（饱和效应会导致非线性变化）
• ψf：永磁体磁链（随温度每升高100°C约减小4%）
• ωe：电角速度（与机械转速关系ωe=P·ωm）

实际工程中需注意：Ls参数辨识时应考虑电流饱和效应，建议在额定电流的20%-120%范围内多点测量。

2.2 定子磁链的物理构成

通过坐标变换得到的定子磁链方程揭示了其双重来源：

code复制ψα = Ls·iα + ψf·cosθe  
ψβ = Ls·iβ + ψf·sinθe

• 电流产生分量（Ls·iαβ）：与负载电流成正比，体现电枢反应
• 永磁耦合分量（ψf·cosθe）：转子位置的正余弦函数，包含待估计信息

2.3 状态变量的优选策略

选择定子磁链作为状态变量的优势在于：

1. 直接包含转子位置信息θe
2. 导数项可通过电压方程消去未知量
3. 满足李雅普诺夫稳定性分析要求

具体状态定义：

code复制x = [ψα, ψβ]^T = [Ls·iα + ψf·cosθe, Ls·iβ + ψf·sinθe]^T

3. 状态动态方程的建立与简化

3.1 状态导数推导过程

对状态变量求导得到：

code复制dψα/dt = Ls·diα/dt - ωe·ψf·sinθe
dψβ/dt = Ls·diβ/dt + ωe·ψf·cosθe

将电压方程中的电流导数项代入：

code复制Ls·diα/dt = uα - Rs·iα + ωe·ψf·sinθe
Ls·diβ/dt = uβ - Rs·iβ - ωe·ψf·cosθe

3.2 关键简化步骤

通过代数运算神奇地消去了包含ωe和θe的项：

code复制dψα/dt = (uα - Rs·iα) + [ωe·ψf·sinθe - ωe·ψf·sinθe] = uα - Rs·iα
dψβ/dt = (uβ - Rs·iβ) - [ωe·ψf·cosθe - ωe·ψf·cosθe] = uβ - Rs·iβ

最终得到简洁形式：

code复制dx/dt = y = [uα - Rs·iα, uβ - Rs·iβ]^T

注意：Rs的准确性直接影响方程有效性。建议采用递推最小二乘法在线辨识，更新周期≤100ms。

4. 转子磁链约束与观测器设计

4.1 转子磁链的提取方法

定义向量函数η(x) = x - Ls·i，代入状态变量可得：

code复制η(x) = ψf·[cosθe, sinθe]^T

其几何意义为：

• 幅值恒等于ψf（关键约束条件）
• 方向与转子d轴重合（包含位置信息）

4.2 幅值误差函数的构造

设计误差函数衡量估计偏差：

code复制ε(x̂) = 0.5*(ψf² - ||η(x̂)||²)

当||η(x̂)||=ψf时，ε=0，此时估计准确。该函数具有以下特性：

• 连续可微，适合梯度优化
• 全局最小值唯一
• 对幅值偏差敏感

4.3 梯度流校正项设计

采用梯度下降法最小化ε²，得到校正项：

code复制dx̂/dt|correction = γ·ε·η(x̂)

其中增益γ的选择原则：

• 过大导致超调振荡
• 过小收敛缓慢
• 建议取值范围：0.5ψf⁻² ~ 2ψf⁻²

5. 完整观测器实现方案

5.1 观测器动态方程

结合开环预测和闭环校正：

code复制dx̂/dt = y + (λ/2)·η(x̂)·[ψf² - ||η(x̂)||²]

参数λ的整定方法：

1. 初始取λ=1/(2τ)，τ为期望收敛时间
2. 通过阶跃响应测试调整
3. 最终值通常为50~200

5.2 位置信息提取算法

通过反正切计算估计位置：

code复制θ̂e = atan2(ψ̂β - Ls·iβ, ψ̂α - Ls·iα)

实际实现时需要处理的问题：

• 象限判断（atan2函数自动解决）
• 相位补偿（考虑数字控制延迟）
• 滤波处理（截止频率≥10倍带宽）

5.3 锁相环(PLL)速度估算

采用二阶PLL结构：

code复制εθ = ψ̂rβ·cosθ̂e - ψ̂rα·sinθ̂e
ω̂e = Kp·εθ + Ki·∫εθ dt
θ̂e = ∫ω̂e dt

参数整定经验：

• 带宽设为系统带宽的1/5~1/3
• Kp=2ξωn, Ki=ωn²（ξ=0.7~1.0）
• 加入抗饱和措施

6. 工程实现关键问题

6.1 参数敏感性分析

主要影响因素及应对措施：

6.2 离散化实现要点

采用双线性变换离散化：

code复制x̂[k+1] = x̂[k] + T·(y[k] + λ/2·η(x̂[k])·[ψf² - ||η(x̂[k])||²])

采样周期T的选择原则：

• 满足T ≤ 1/(10·ωb)，ωb为系统带宽
• 典型值：50-100μs（对应10-20kHz PWM）

6.3 启动策略设计

特殊处理方案：

1. 初始位置检测（高频注入或DC激励）
2. 开环启动阶段（固定加速度斜坡）
3. 观测器软启动（λ从0渐增至额定值）
4. 切换逻辑（速度误差<5%时切换）

7. 实测性能优化记录

在某400W伺服电机上的调试数据对比：

关键优化手段：

1. 引入自适应λ：λ=λ0·(1+0.2|ω̂e|/ωrated)
2. 增加谐波补偿项
3. 采用滑动平均滤波（窗口5-10个周期）

经过三个月现场运行验证，该方案在0-3000rpm范围内位置估计误差始终保持在±3°以内，完全满足工业伺服定位要求。特别是在低速区段的稳定性，使得机床主轴定向停止精度提升40%以上。