直流电机H∞控制实战：从状态建模到鲁棒控制器设计

简介：本资源是一份面向控制工程领域研究生、科研人员及工程师的H∞鲁棒控制实战资料，聚焦直流电机在参数不确定性与外部扰动下的高性能闭环控制问题，系统覆盖状态空间建模、广义被控对象构建、权函数设计、H∞控制器综合与MATLAB仿真验证全流程。压缩包共590个文件，含482个核心m脚本（实现控制器设计、仿真主程序与分析函数）、26个mat数据文件（存储系统参数与仿真结果）、20个fig图形文件（直观呈现阶跃响应、正弦跟踪及抗扰对比曲线），辅以C/C++底层接口（.c/.cpp/.dll/.mex*）支持混合仿真，整体7.02MB，结构完整、模块解耦清晰。已有127人学习下载，提供从理论推导到代码落地的全链路支撑，包括可复现的完整仿真案例、权函数选型依据说明、控制器参数敏感性分析及工程调试建议，助读者深入掌握H∞控制在机电系统中的工程实现路径。

1. H∞控制不是“万能滤波器”，而是给直流电机装上可量化的抗扰铠甲

很多工程师第一次接触H∞控制时，会下意识把它当成一种“高级PID”——参数调得越小，系统越稳；加个权函数，就能压住所有干扰。但实际跑通直流电机仿真后才发现：阶跃响应超调压不下去，正弦跟踪相位滞后严重，甚至控制器增益突然爆表。问题出在哪？根本原因在于H∞不是靠经验试凑，而是通过状态空间模型的广义被控对象重构 + 权函数对性能/鲁棒性边界的显式刻画，把“系统能承受多大参数摄动”和“对外部扰动能抑制到什么程度”变成可计算、可验证的数学不等式。本文提供的MATLAB实现，正是围绕直流电机电枢回路与机械运动耦合特性，构建包含电阻变化±20%、负载转矩突变、测量噪声三类不确定性的真实闭环场景，用hinfsyn求解器直接输出满足γ<1.8的最小范数控制器。适合已掌握状态方程建模、能手推传递函数到状态空间转换的控制工程师，尤其适用于伺服驱动器、精密定位平台等对参数漂移敏感的工业现场。

2. 从直流电机物理方程到广义被控对象：状态空间建模的三个关键跃迁

H∞控制器设计成败，70%取决于广义被控对象（Generalized Plant）的构建质量。它不是简单把电机传递函数塞进sys = ss(A,B,C,D)，而需完成三次本质性建模跃迁：第一，将非线性摩擦、电感饱和等忽略项转化为结构化不确定性Δ；第二，把控制目标（如跟踪精度、抗扰带宽）和硬件约束（如PWM限幅、传感器带宽）映射为权重矩阵W1、W2、W3；第三，确保所有不确定性通道在LFT（Linear Fractional Transformation）框架下可解耦。下面以永磁直流电机为例，逐步展开。

2.1 基础状态空间模型：从物理定律到线性化方程

直流电机经典二阶模型由电枢电压方程和运动方程导出：

di_a/dt = (-R/L)·i_a - (K_e/L)·ω + (1/L)·u_a dω/dt = (K_t/J)·i_a - (B/J)·ω

其中i_a为电枢电流，ω为角速度，R=2.5Ω、L=0.015H、K_e=K_t=0.15N·m/A、J=0.002kg·m²、B=0.001N·m·s/rad。将上述方程写成标准状态空间形式：

A0 = [-R/L, -K_e/L; K_t/J, -B/J]; B0 = [1/L; 0]; C0 = [0, 1]; % 仅测量角速度 D0 = 0; sys_nominal = ss(A0, B0, C0, D0);

注意：此处C矩阵设为[0,1]表示只反馈角速度，若需电流反馈则改为[1,0]或[1,1]。D0=0是因理想电机无直通路径，但实际中需考虑ADC采样延迟，此时D应设为e^(-sT_s)的Padé近似。

2.2 不确定性建模：用Δ块封装参数摄动与外部扰动

H∞框架要求所有不确定性必须显式表达为Δ，且满足||Δ||_∞ ≤ 1。针对直流电机，我们定义三类Δ：

• 参数不确定性Δ1：电枢电阻R在[2.0,3.0]Ω间变化 → 构造R = R0 + δ_R·R0，δ_R∈[-0.2,0.2] → Δ1 = diag(δ_R)
• 外部扰动Δ2：负载转矩突变 → 在运动方程右侧添加-1/J·d(t)→ 将d(t)作为外部输入w1，对应Bw1=[0;-1/J]
• 测量噪声Δ3：编码器量化误差 → 在输出端叠加v(t) → Cz=[0,1;0,0]，Dzw=[0;1]，使z=[ω;v]

最终广义被控对象P(s)维度为[z; y] = P(s)·[w; u]，其中z为加权后的性能输出（含跟踪误差、控制量），w为扰动输入（含d(t), v(t)），y为实际测量输出（ω）。MATLAB中构造如下：

% 定义权重矩阵（核心！） W1 = makeweight(1, 0.1, 10); % 跟踪误差权重：低频增益1，高频衰减 W2 = makeweight(0.01, 0.001, 100); % 控制量权重：抑制u幅值 W3 = tf(1,[0.01,1]); % 测量噪声权重：强调高频抑制 % 构建广义被控对象P P = augw(sys_nominal, W1, W2, W3); % 此时P为4输入2输出系统：w=[d;v], u=controller output, z=[e;u], y=omega

提示：makeweight生成的权重函数形状直接影响控制器带宽。若W1低频增益过小（如0.1），会导致阶跃响应稳态误差增大；若W2高频衰减太慢，则控制量会出现剧烈抖动。实际调试中需反复调整W1、W2的转折频率。

2.3 广义被控对象验证：检查LFT结构与维度匹配

构建完P后，必须验证其LFT分解是否正确。H∞控制器K满足z = P11·w + P12·u,y = P21·w + P22·u，其中P11、P12等子块需符合物理意义。使用lftdata函数提取：

[P11,P12,P21,P22] = lftdata(P); size(P11) % 应为2×2：z有2维（误差+控制量），w有2维（扰动+噪声） size(P12) % 应为2×1：u为单输入 size(P21) % 应为1×2：y为单输出 size(P22) % 应为1×1：标量

若尺寸不匹配，常见原因是augw输入的sys_nominal未正确指定输入输出通道。例如，若原始sys_nominal的输入端口未标记为'u'，输出未标记为'y'，augw会默认按顺序连接，导致P22维度错误。解决方案是显式定义：

sys_nominal.InputName = 'u'; sys_nominal.OutputName = 'y'; P = augw(sys_nominal, W1, W2, W3);

3. H∞控制器求解与参数整定：hinfsyn的四个不可跳过的配置项

hinfsyn是MATLAB Robust Control Toolbox中求解标准H∞问题的核心函数，但它绝非“一键生成”。其输出控制器K的质量，高度依赖四个关键配置项的设置。本节结合直流电机案例，逐项解析参数含义与调试逻辑。

3.1 指定性能指标γ与迭代容差：平衡鲁棒性与可行性

H∞问题本质是求解最小γ使得||F_l(P,K)||_∞ < γ成立，其中F_l为闭环传递函数。hinfsyn默认尝试γ=1开始迭代，但直流电机因存在右半平面零点（由电感L引入），常需γ>1.5才能收敛。必须显式设置初始γ和容差：

options = hinfsynOptions('MaxIter', 100, 'RelTol', 1e-5); [K, CL, gamma, info] = hinfsyn(P, 1, 1, options); % 第二、三参数分别为控制输入数、测量输出数

若gamma > 2.5，说明权重设计过于激进（如W1低频增益太大），需降低W1增益或拓宽W2带宽；若info.Stabilizing == 0，表明控制器无法稳定P，此时应检查P的极点是否含不稳定模态（如未补偿的机械谐振）。

3.2 处理控制器降阶：避免高阶控制器引发实时性危机

hinfsyn默认输出控制器阶数与P相同。本例中P为4阶，但实际嵌入式控制器（如TI C2000系列）通常要求≤3阶。必须启用降阶选项：

options.ReduceOrder = 'on'; options.Order = 2; % 强制降为2阶 [K_red, CL_red, gamma_red, info_red] = hinfsyn(P, 1, 1, options);

降阶后需对比bode(K)与bode(K_red)：若降阶控制器在100Hz以上增益骤降，则可能削弱高频抗扰能力。此时应检查info_red.Rank——若小于设定阶数，说明原控制器本身已接近低阶，强行降阶会损失性能。

3.3 权重灵敏度分析：用sigma plot定位瓶颈环节

权重选择是否合理，不能仅看仿真曲线，而需分析闭环系统的奇异值（σ）图。执行：

sigma(CL); grid on; title('Closed-loop singular values'); legend('sigma_{max}(CL)','sigma_{min}(CL)');

关键观察点：

• 低频段（<1Hz）σ_max应低于0dB，确保阶跃无稳态误差；
• 中频段（1~50Hz）σ_max需快速下降，反映跟踪带宽；
• 高频段（>100Hz）σ_max应趋近于γ值，验证鲁棒裕度。

若σ_max在10Hz处出现尖峰，说明W2抑制不足，控制量在此频段共振；若整体σ_max高于γ，表明hinfsyn未收敛，需重新运行。

3.4 控制器离散化：ZOH与Tustin方法对相位的影响

仿真用连续控制器，但实际部署需离散化。c2d函数提供多种方法，对直流电机影响最大的是相位特性：

Ts = 1e-4; % 10kHz采样率 Kd_zoh = c2d(K_red, Ts, 'zoh'); % 零阶保持：相位滞后小，但幅频响应有混叠 Kd_tus = c2d(K_red, Ts, 'tustin'); % 双线性变换：幅频保真度高，但相位超前易引发振荡

实测发现：当电机机械谐振频率为250Hz时，tustin法在125Hz处引入额外相位超前，与谐振峰叠加导致不稳定；而zoh法虽在200Hz以上增益衰减，但相位单调滞后，更易稳定。因此本例选用zoh，并验证离散闭环脉冲响应：

CLd = feedback(Pd*Kd_zoh, 1); % Pd为离散化后的广义对象 impulse(CLd, 0.1); % 观察前0.1秒冲击响应是否收敛

4. 仿真验证三要素：阶跃响应、正弦跟踪、抗扰测试的量化判据

控制器设计完成，必须通过三项硬性测试验证其鲁棒性。每项测试均需定义明确的量化判据，而非仅观察曲线形态。本节提供可直接复用的MATLAB脚本框架与阈值标准。

4.1 阶跃响应测试：聚焦超调量、调节时间与稳态误差

在simulink中搭建闭环系统，施加单位阶跃参考信号，记录角速度响应。关键判据：

• 超调量σ% ≤ 8%：超过此值说明相位裕度不足，需增强W2在中频段的抑制；
• 调节时间t_s（2%准则）≤ 0.15s：反映带宽，若超限则扩大W1的高频衰减带宽；
• 稳态误差e_ss ≤ 0.002 rad/s：对应W1低频增益≥100。

% 生成阶跃响应数据 t = 0:Ts:1; r = ones(size(t)); % 单位阶跃 [y,tout,x] = lsim(CL, r, t); % CL为闭环传递函数 % 计算指标 ys = y(end); % 稳态值 ess = 1 - ys; peak = max(y); sigma_pct = (peak - 1)/1 * 100; ts_idx = find(abs(y - 1) < 0.02*1, 1, 'first'); ts = t(ts_idx); fprintf('超调:%.2f%%, 调节时间:%.3fs, 稳态误差:%.4f\n', sigma_pct, ts, ess);

4.2 正弦跟踪测试：评估幅值衰减与相位滞后

输入频率为10Hz、20Hz、50Hz的正弦信号，测量输出幅值比与相位差。合格标准：

• 10Hz时幅值比 ≥ 0.95，相位滞后 ≤ 15°：保证低频跟踪精度；
• 50Hz时幅值比 ≥ 0.7，相位滞后 ≤ 45°：验证带宽余量。

freq_list = [10, 20, 50]; for i = 1:length(freq_list) f = freq_list(i); t_sin = 0:Ts:0.5; r_sin = sin(2*pi*f*t_sin); [y_sin,~,~] = lsim(CL, r_sin, t_sin); % 提取稳态段（最后0.1秒） idx_steady = find(t_sin > 0.4, 1, 'first'):end; Y_fft = fft(y_sin(idx_steady)); R_fft = fft(r_sin(idx_steady)); amp_ratio(i) = abs(Y_fft(2))/abs(R_fft(2)); % 基频幅值比 phase_diff(i) = angle(Y_fft(2)) - angle(R_fft(2)); % 相位差（rad） end

4.3 抗扰性能测试：量化扰动抑制比（DSR）

在t=0.3s时刻注入幅值为额定转矩20%的阶跃负载扰动d(t)，测量角速度波动峰值Δω_max与恢复时间t_recovery。定义扰动抑制比DSR = |Δω_max / d_step|，要求DSR ≤ 0.05（即5%转速波动）：

% 构建含扰动的仿真模型 P_dist = P; % P已含扰动输入w1 P_dist.InputName{1} = 'd'; % 显式命名扰动通道 CL_dist = feedback(P_dist*Kd_zoh, 1, 1, 1); % 1,1指d为第一输入，u为第一输出 % 施加扰动 t_dist = 0:Ts:0.6; d_step = [zeros(1,3000), 0.2*ones(1,1000)]; % 0.3s后注入20%扰动 [y_dist,~,~] = lsim(CL_dist, d_step, t_dist); delta_omega = y_dist - y_dist(1); % 扰动引起的转速变化 dsr = max(abs(delta_omega)) / 0.2; fprintf('扰动抑制比DSR=%.4f\n', dsr);

提示：若DSR超标，优先检查W3（测量噪声权重）是否过大——过大的W3会弱化控制器对扰动的响应速度。可尝试将W3的DC增益降低至0.1，重新运行hinfsyn。

5. 工程落地关键技巧：从仿真到DSP的三步参数迁移法

仿真验证通过后，控制器参数需迁移到实际硬件。直接复制K的分子分母系数常导致性能劣化，根本原因是仿真中忽略了ADC采样延迟、PWM死区、电流环带宽等非理想因素。本节提供经某伺服驱动器项目验证的三步迁移法。

5.1 步骤一：用d2c反向离散化，获取连续域控制器结构

即使部署离散控制器，也需先还原其连续等效形式，以便分析零极点分布：

K_cont = d2c(Kd_zoh, 'tustin'); % 必须用tustin，zoh反变换不唯一 [z,p,k] = zpk(K_cont); fprintf('连续控制器零点:%.3f±j%.3f, 极点:%.3f±j%.3f\n', ... real(z(1)), imag(z(1)), real(p(1)), imag(p(1)));

若发现控制器零点实部>0（右半平面零点），说明存在非最小相位特性，需在硬件中增加前置滤波器补偿。

5.2 步骤二：注入真实延迟模型，重优化离散控制器

在P中加入1.5个采样周期的纯延迟（对应ADC+DSP处理延时）：

delay_T = 1.5 * Ts; P_delay = P * ss(1,'InputName','w','OutputName','w_delay'); % 占位 P_delay = series(P_delay, exp(-s*delay_T)); % 添加延迟 [K_final,~,~,~] = hinfsyn(P_delay, 1, 1, options);

此步骤使控制器主动补偿延迟，避免现场调试时反复修改。

5.3 步骤三：量化系数适配定点DSP，避免溢出

TI C2000系列常用Q15格式（15位小数），需将控制器系数缩放：

% 获取离散控制器系数 [num,den] = tf(K_final); % 计算最大系数绝对值 max_coef = max([abs(num), abs(den)]); scale_factor = 2^15 / max_coef; num_q15 = round(num * scale_factor); den_q15 = round(den * scale_factor); % 生成C代码头文件 fprintf(fid, '#define K_NUM_COEF {%d,%d}\n', num_q15); fprintf(fid, '#define K_DEN_COEF {%d,%d,%d}\n', den_q15);

注意：缩放后必须验证闭环稳定性，方法是将Q15系数代入tf(num_q15/scale_factor, den_q15/scale_factor)，再执行margin查看相位裕度是否>45°。若裕度下降超10°，需增大scale_factor或改用Q31格式。

最终，在某数控机床主轴驱动器上实测：H∞控制器相较传统PI，在电机电阻升高30%时，转速波动从±12rpm降至±1.8rpm；面对50ms阶跃负载，恢复时间缩短42%。这印证了——鲁棒性不是玄学指标，而是可通过状态空间建模、权重函数设计、闭环奇异值分析精确调控的工程参数。

本文还有配套的精品资源，点击获取