当前位置：首页 > news >正文

基于Simulink的鲁棒H∞整流器控制器设计

news 2026/5/12 11:50:51

手把手教你学Simulink

——基于Simulink的鲁棒H∞整流器控制器设计

一、问题背景

二、系统建模：dq 坐标系下的状态空间

三、加权函数设计（Weighting Functions）

1. 跟踪误差权重 (W_1(s))

2. 控制量权重 (W_2(s))

3. 扰动权重 (W_d(s))（可选）

四、广义被控对象构建（P(s)）

五、H∞ 控制器综合

六、Simulink 闭环系统搭建

第一步：导出控制器到 Simulink

第二步：构建完整整流系统

第三步：注入扰动测试鲁棒性

七、仿真结果与分析

性能指标对比（H∞ vs PI）

频域特性（Bode 图）

八、工程挑战与应对

九、总结

十、动手建议

手把手教你学Simulink

——基于Simulink的鲁棒H∞整流器控制器设计

一、问题背景

在高性能三相 PWM 整流器中，传统 PI 控制虽结构简单，但面临以下挑战：

模型不确定性：电感、电阻随温度/频率变化
外部扰动：电网电压谐波、不平衡、跌落
多目标冲突：既要快速动态响应，又要强抗扰性

H∞ 控制（H-infinity Control）是一种频域鲁棒控制方法，通过优化系统对扰动的最坏情况增益（即 ( |T_{zd}|_\infty )），实现：

对模型摄动的鲁棒稳定性
对外部扰动的鲁棒性能
显式处理多目标权衡（如带宽 vs 噪声抑制）

本教程将在 Simulink + MATLAB 中完成：

整流器小信号建模
加权函数设计
H∞ 控制器综合（使用hinfsyn）
闭环系统仿真与鲁棒性验证

先决条件：需安装Robust Control Toolbox和Simulink Control Design

二、系统建模：dq 坐标系下的状态空间

考虑 d 轴电流子系统（q 轴同理），忽略交叉耦合项（或已前馈补偿），其线性化模型为：

[
L \frac{di_d}{dt} = -R i_d + v_d - v_{gd}
]

整理为标准状态空间形式：

[
\begin{aligned}
\dot{x} &= A x + B_1 w + B_2 u \
z &= C_1 x + D_{11} w + D_{12} u \
y &= C_2 x + D_{21} w + D_{22} u
\end{aligned}
]

其中：

状态 (x = i_d)
控制输入 (u = v_d)
外部扰动 (w = [v_{gd},, n]^\top)（电网电压扰动 + 电流测量噪声）
被控输出 (z = [W_1(s)(i_d^* - i_d),, W_2(s) v_d]^\top)（跟踪误差 + 控制量）
测量输出 (y = i_d)

取参数（标称值）：

(L = 2,\text{mH}), (R = 0.1,\Omega)

则：
[
A = -\frac{R}{L} = -50,\quad
B_2 = \frac{1}{L} = 500,\quad
B_1 = \begin{bmatrix} -\frac{1}{L} & 0 \end{bmatrix} = \begin{bmatrix} -500 & 0 \end{bmatrix}
]

三、加权函数设计（Weighting Functions）

H∞ 设计的核心是合理构造加权函数，将时域指标转化为频域约束。

1. 跟踪误差权重 (W_1(s))

要求低频高增益（稳态无静差），高频衰减（抑制噪声）：

[
W_1(s) = \frac{s/M + \omega_b}{s + \omega_b \epsilon}
]

(\omega_b)：期望带宽（如 1000 rad/s）
(M)：最大峰值（如 1.5，限制超调）
(\epsilon)：低频增益（如 0.01，保证稳态精度）

wb = 1000; % 带宽 (rad/s) M = 1.5; % 最大灵敏度 eps = 0.01; % 低频增益 W1 = tf([1/M, wb], [1, wb*eps]);

2. 控制量权重 (W_2(s))

限制高频控制能量（避免放大噪声、保护器件）：

[
W_2(s) = \frac{s}{\omega_u} + 1
]

(\omega_u)：控制量转折频率（如 5000 rad/s）

wu = 5000; W2 = tf([1/wu, 1], 1);

3. 扰动权重 (W_d(s))（可选）

若已知扰动频谱（如电网 5th 谐波），可加入 (W_d) 提升针对性。

四、广义被控对象构建（P(s)）

将原系统与加权函数互联，形成广义 plant (P(s))：

┌───────────────┐ w ──►│ │──► z │ P(s) │ u ──►│ │──► y └───────────────┘

在 MATLAB 中使用augw或手动构建：

% 标称 plant G(s) = 1/(Ls + R) s = tf('s'); G = 1/(0.002*s + 0.1); % 构建 P = [W1*G, W1; W2, 0] （简化版，忽略扰动通道细节） P = augw(G, W1, W2, []);

更严谨的方式是显式定义所有通道：

A = -50; B1 = [-500; 0]; % [vgd; noise] B2 = 500; C1 = [1; 0]; % [tracking error; control effort] C2 = 1; % measurement D11 = [0; 0]; D12 = [0; 1]; D21 = [0, 1]; % assume noise directly adds to measurement D22 = 0; P = ss(A, [B1 B2], [C1; C2], [D11 D12; D21 D22]); P = append(W1, W2, 1) * P; % apply weights