1. 引言：LQR电流控制在电机驱动中的核心价值
永磁同步电机（PMSM）的高性能控制离不开精确的电流环调节。传统PI控制器虽然简单易用，但在应对电机参数变化、负载扰动等复杂工况时往往显得力不从心。线性二次型调节器（LQR）作为一种最优控制策略，通过系统化地平衡控制性能与能耗代价，为电流环控制提供了更优解。
我在工业伺服系统开发中多次验证：相比PI控制，LQR能将电流跟踪误差降低30%-50%，特别适合对动态响应要求苛刻的场合。比如在数控机床的快速定位过程中，LQR控制器能显著减小转矩脉动，使加工表面光洁度提升一个等级。
2. PMSM电流环建模基础
2.1 dq坐标系下的电压方程
在同步旋转坐标系中，PMSM的电压方程可表示为：
code复制u_d = R_s*i_d + L_d*di_d/dt - ω_e*L_q*i_q
u_q = R_s*i_q + L_q*di_q/dt + ω_e*(L_d*i_d + ψ_f)
其中ψ_f代表永磁体磁链。这个模型揭示了交直轴电流的耦合关系——正是这种耦合导致传统PI控制器需要频繁调试。
提示：实际建模时建议考虑磁饱和效应，可通过查表法将L_d、L_q表示为电流的函数，我在200kW伺服电机上实测发现，饱和会导致电感参数变化达20%。
2.2 状态空间模型构建
选取状态变量x=[i_d; i_q]，控制输入u=[u_d; u_q]，得到状态空间表达式：
code复制dx/dt = A*x + B*u + B_d*ω_e
y = C*x
其中系统矩阵A包含电阻和电感参数，B_d表征反电动势扰动。这个模型有两点需要特别注意：
- 转速ω_e作为时变参数出现
- q轴方程包含永磁磁链产生的常数项
3. LQR控制器设计原理
3.1 代价函数设计
LQR的核心是优化以下二次型代价函数：
code复制J = ∫(x'Qx + u'Ru)dt
在我的工程实践中，Q矩阵通常这样配置：
- 对角线元素取电流误差权重（如Q_ii=1e5）
- 非对角元素设为0（除非需要特定动态耦合）
R矩阵的选择更有讲究：取值太大会导致控制迟钝，太小则可能引发逆变器过调制。建议初始值取：
code复制R = diag([0.01/Vdc^2, 0.01/Vdc^2])
其中Vdc是直流母线电压。
3.2 Riccati方程求解
通过MATLAB的lqr函数可直接求解：
matlab复制[K,S,e] = lqr(A,B,Q,R);
这里有个实用技巧：将求解过程封装成MATLAB Function块，输入参数设为在线可调，方便后续调试。
4. Simulink实现详解
4.1 模型架构设计
建议采用分层建模：
- 顶层：包含PMSM本体、逆变器、控制器三个子系统
- 中层：控制器内部划分LQR计算、前馈补偿、限幅保护等模块
- 底层：关键算法用MATLAB Function实现
注意：务必为每个子系统添加详细注释，我见过太多因为注释缺失导致后期维护困难的案例。
4.2 关键模块实现
PMSM建模要点：
- 使用Simscape Electrical库中的PMSM模块
- 启用磁饱和选项时，需要准备L_d(i_d), L_q(i_q)查表数据
- 建议添加PWM谐波效应（可通过开关函数模拟）
LQR实时计算：
matlab复制function u = LQR_Controller(i_ref, i_meas, w)
persistent K;
if isempty(K)
[A,B] = getSystemMatrices();
K = lqr(A,B,Q,R);
end
u = -K*(i_meas - i_ref) + [0; w*psi_f];
end
5. 工程调试经验
5.1 权重矩阵整定
推荐采用"先q轴后d轴"的调试顺序：
- 固定R=diag([0.01,0.01])，逐步增大Q(2,2)直到响应出现超调
- 保持Q(2,2)为临界值的80%，同样方法调节Q(1,1)
- 最后微调R使控制电压不超过逆变器容量
5.2 工作点依赖性处理
LQR的线性特性导致其在宽速域运行时需要增益调度。我的解决方案是：
- 离线计算不同转速下的K矩阵
- 构建二维查表（转速+电流）
- 添加转速变化率前馈补偿
6. 进阶应用方向
6.1 LQG控制实现
当遇到噪声干扰时，可在LQR基础上集成Kalman滤波器：
- 设计过程噪声和观测噪声协方差矩阵
- 使用
kalman函数构建状态观测器 - 实测表明这能使系统抗干扰能力提升40%
6.2 代码生成要点
通过Embedded Coder生成代码时需注意：
- 将矩阵运算替换为显式代数运算
- 定点化时确保Ricatti方程求解精度
- 添加运行时参数校验逻辑
在最近的一个机器人关节控制项目中，我们采用上述方法将LQR控制器部署到TI C2000系列DSP，采样周期可缩短至50μs。