超螺旋高阶滑模+双观测器：Buck变换器负载突变抑制方案

简介：面向电力电子工程师、自动化控制研究人员及高校相关专业师生，提供Buck型变换器高阶滑模控制方法的完整研究资料。内容涵盖Buck变换器工作原理与数学模型建立、传统一阶滑模/积分滑模/终端滑模控制性能对比、二阶滑模抖振抑制设计，以及负载观测器和非线性干扰观测器的具体应用。全文结合详细可运行的MATLAB仿真代码和逐段中文解释，帮助读者理解滑模面构造、控制律平滑处理（如以tanh替换sign）、观测器极点配置与干扰估计补偿等关键环节。仿真结果表明，在±40%输入电压波动工况下输出电压偏差仍小于2%，较传统方法动态响应更快、抗扰能力更强。包内为1个PDF文档，大小695KB，浓缩了理论推导、代码复现、参数整定规则及硬件实现注意事项。已有76人学习，适合需要高精度、强鲁棒性Buck变换器控制方案的实践者参考。 Buck变换器控制，见过太多人卡在同一个问题上：负载一突变，输出电压先砸下去一截，然后慢慢悠悠爬回来。普通PI最怕这种工况，带宽拉高了怕噪声，拉低了又嫌恢复慢；换成滑模控制倒是恢复快，可高频抖振一上来，听着电感“吱吱”叫，心里更没底。这篇文章，我直接讲清楚我自己在12V/5A样机上用的一套方案——用超螺旋高阶滑模把抖振压下去，再叠加上负载观测器和非线性干扰观测器把扰动“吃干净”，整个过程带完整Matlab代码和参数整定经验，适合做数字电源、毕设选题里碰到的Buck控制场景，也适合已经学过经典滑模、想往高阶滑模方向落地的朋友。

1. 本文的对标问题：负载突变下的电压跌落与控制困局

先还原一个典型测试画面：Buck输出12V，静态带上4Ω负载，电压纹波挺好。你拿起电子负载，把负载从4Ω瞬间切到2Ω，示波器上输出电压先是掉下去一二百毫伏，然后经过几个毫秒慢慢回到12V。如果这个跌落幅度超过后端芯片的允许范围，下一级设备可能直接复位。

为什么线性控制器这么“肉”？因为电压环本质上是靠输出电容放电来扛住第一波负载冲击的，环路只能事后补偿。把PI参数调激进一些，跌落确实小一些，但相位裕度下降，动态过程中容易振荡；而且PWM占空比变化太猛，电磁干扰问题也随之而来。

一阶滑模控制（SMC）看起来是解药——它对参数摄动和外部扰动天然不敏感。但用过的人都知道，经典滑模里的符号函数项会直接进占空比，控制量高频切换，稳态时占空比波形像锯齿，实测效率下降，电磁噪声也大。根本原因在于：一阶滑模为了保证鲁棒性，切换增益必须大于扰动上界，而扰动上界又往往被夸大，结果就是抖振。

所以现实中最合理的路线，不是单纯提高切换增益去“硬抗”扰动，而是先用观测器把扰动估计出来，通过前馈抵消，让滑模切换增益维持在很低水平；再用高阶滑模算法把不连续项放到控制量的导数里去，从机制上消除抖振。这是这篇项目里一以贯之的设计逻辑。

2. Buck变换器的平均模型与“集总扰动”的数学定义

2.1 连续导通模式下的状态平均方程

先建立Buck在CCM（连续导通模式）下的平均模型。设输入电压为V_in，开关占空比为u，电感L，电容C，负载电阻R。电感电流i_L和电容电压v_C为状态变量：

L·(di_L/dt) = V_in·u - v_C

C·(dv_C/dt) = i_L - v_C/R

这个模型忽略了开关纹波，把每个开关周期内的平均值当作对象。对控制设计来说这是合理的，因为电压环穿越频率通常只有开关频率的十分之一到二十分之一，开关纹波是高频分量，环路本来就不该去响应它。

重点在于，我们要控制的是输出电压v_C，自然希望得到一个以v_C为输出的二阶方程。把上面两式合并，并假设负载R实际是时变的，可以得到误差动态的完整形式。

2.2 把不确定性全部塞进一个“集总扰动”

定义输出电压误差e = V_ref - v_C，其中V_ref为参考电压。假设参考电压为常量，则有：

ë = v_C/(LC) - ė/(R·C) - (V_in·u)/(LC) + d_total

这里的d_total就是“集总扰动”，它包含了真实设计时最头疼的所有东西：

• 负载电阻R的变化及其变化率；
• 输入电压V_in的波动，因为方程里的V_in如果偏离标称值，差值部分就进入d_total；
• 电感L、电容C的参数漂移，比如电感饱和、电容老化；
• PWM死区、开关延迟、采样噪声带来的非理想效应。

在控制设计时，我先把d_total当作一个未知但有界的量，后面用观测器去估计，而不是靠高增益滑模去压制。划个重点：这里把“负载突变”这种强扰动从模型里显式分离出来，剩下的交给观测器，是整个方案能压住电压跌落的关键。

3. 超螺旋滑模控制：如何把不连续项藏进导数里

3.1 普通滑模哪里疼

经典一阶滑模，通常选滑模面s = ė + c·e，其中c>0决定误差收敛速率。控制律由等效控制u_eq和切换控制u_sw组成，切换项一般是这类形式：

u_sw = -k·sign(s)

问题一眼就能看出来：sign(s)直接作用在占空比u上，而占空比本身是物理量，没法无限频率切换。数字控制器里虽然有开关频率上限，但这个符号函数会让占空比在高频段来回跳，实际表现就是电感啸叫、输出纹波恶化、EMI超标。

3.2 超螺旋算法：二阶滑模的最实用形态

真正工程上用得最多的高阶滑模，是超螺旋算法（Super-Twisting Algorithm, STA）。它属于二阶滑模控制：不连续项不再直接出现在控制量u中，而是先进入一个中间变量w的导数里，再由w积分后形成控制量。结构如下：

s = ė + c·e

w = -λ·|s|^0.5·sign(s) + ξ

dξ/dt = -W·sign(s)

其中λ和W是待整定的正增益。从公式能直观看到：sign(s)出现在ξ的导数里，而控制量u用的是经过积分后的w。积分本身是低通特性，所以u是连续的，抖振从根源上被削掉了。

很多新手会问：这和“把符号函数后面加个低通滤波”有什么区别？区别在于超螺旋算法有严格的可证明的有限时间收敛性，它保证在滑模面上误差的高阶量也收敛到零，增益条件有理论支撑。而我实际调试中感受到的差异是：超螺旋不需要滑模变量的导数信息，只用到s本身，离散实现非常简单，一个采样周期就能跑完。

3.3 稳定性条件简述

超螺旋算法有限时间稳定的充分条件，大致可以写成：

λ > 0，W > 某个关于扰动上界的正值

完整的Lyapunov证明在这里不展开，但整定逻辑要说透：λ决定收敛速度，W决定对扰动的抑制强度。如果W太小，系统在滑模面上会出现“慢性子”的偏差；如果W太大，占空比里又能看到低频振荡。我实际调试中通常是先固定W，再一点一点加λ，直到占空比波形干净为止。

4. 双观测器协作：负载观测器和非线性干扰观测器的具体设计

4.1 负载观测器：用物理关系直接辨识R

负载观测器的思路，是利用电容电流方程反推负载电流。电容电流关系是：

i_C = C·(dv_C/dt) = i_L - v_C/R

所以负载电流i_R = v_C/R = i_L - C·(dv_C/dt)。只要知道电感电流i_L、输出电压v_C及其导数，就能算出当前负载电流，进而得到负载电阻估计值：

R_hat = v_C / max(i_R_est, ε)

这里的max()是保护，防止启动阶段或者空载时分母接近零导致R_hat爆掉。

实际工程里，dv_C/dt不能直接差分，噪声会被放大到没法用。我建议的做法是先对v_C做一阶低通滤波，再差分；或者更讲究一点，用一个跟踪微分器。在本项目的Matlab仿真中，我直接用仿真内部的真实电压变化率来模拟理想传感器信号，代码注释里标明了实际替换方案。

负载观测器对负载突变这种“大尺度慢变量”非常有效。负载从4Ω变到2Ω，R_hat会在几个周期内迅速跟踪到接近2Ω，从而修正控制器里的模型参数。

4.2 非线性干扰观测器：吃掉残余的动态扰动

负载观测器只能针对性地估计R，而对输入电压突变、电感参数漂移、负载变化率冲击这类“杂七杂八”的扰动，需要更通用的工具。非线性干扰观测器（NDOB）就是干这个的。

思路很简单：建立误差系统的状态方程

ë = f0 + b0·u + d_total

其中 f0 = v_C/(LC) - ė/(R_hat·C)，b0 = V_in/(LC) 为标称输入通道增益。NDOB用一个辅助变量ζ来避开对d_total导数的直接测量：

ζ̇ = -l·(ζ + p) - l·(f0 + b0·u)

p = l·ė

d_hat = ζ + p

l是观测器增益，本质上决定观测带宽。这个结构的收敛性很好理解：如果l足够大，d_hat的动态远快于d_total的真实变化，它就能实时逼近d_total。代价是l太大会放大测量噪声，这个矛盾第6章会专门讲。

4.3 为什么两个观测器要一起用

有人会问：既然NDOB什么都能估计，还要负载观测器做什么？我实际对比过，只靠NDOB也能工作，但有两个缺点：一是NDOB把负载突变识别成“扰动”，本质上是靠高带宽去追一个阶跃信号，速度再快也有滞后；二是负载突变瞬间d_total里包含R的变化率脉冲，NDOB很难完全逼近，表现在波形上就是那个“尖刺”没法完全削平。

负载观测器提前把R_hat的慢变部分修正好，NDOB只需要盯着残余快变量，压力小很多。这是典型的“慢变参数观测器+快变扰动观测器”分工协作，比单靠一个观测器效果稳。

5. 从公式到仿真：完整Matlab代码与逐段解释

5.1 代码框架和参数说明

下面这套代码我拆成了三块：参数初始化、主循环（含控制器和双观测器）、结果绘图。仿真用的是CCM平均模型，控制频率等于开关频率20kHz，这个方法验证算法完全够用。若想看开关纹波，可以在Simulink里搭功率级，控制器代码可以直接移植。

%% Buck_HOSMC_NDOB_demo.m % 超螺旋高阶滑模 + 负载观测器 + 非线性干扰观测器 % 说明: 仿真采用CCM平均模型, 控制频率=开关频率=20kHz clear; clc; close all; %% 系统参数 V_in = 24; % 输入电压 V V_ref = 12; % 输出电压参考 V L = 330e-6; % 电感 H C = 470e-6; % 电容 F R_nom = 4; % 标称负载 Ohm f_sw = 20e3; % 开关/控制频率 Hz Ts = 1/f_sw; % 控制周期 s %% 控制器与观测器参数 c = 4000; % 滑模面斜率 lambda_w = 0.5; % 超螺旋比例增益 W_h = 1200; % 超螺旋积分增益 l_ndob = 2000; % 非线性干扰观测器带宽增益 alpha_R = 0.05; % 负载观测器滤波系数 %% 仿真时间 T_end = 0.06; N = round(T_end / Ts); t = (0:N-1) * Ts; %% 初始状态 iL = V_ref / R_nom; % 稳态初始: 电感电流 vC = V_ref; % 稳态初始: 输出电压 R_hat = R_nom; % 负载观测器初值 d_hat = 0; % 非线性干扰观测器初值 z_d = 0; % NDOB辅助状态 xi = 0; % 超螺旋积分项 u_old = V_ref / V_in; % 上一周期占空比(初始为标称值) % 存储 u_log = zeros(1,N); vC_log = zeros(1,N); e_log = zeros(1,N); dhat_log = zeros(1,N); Rhat_log = zeros(1,N); for k = 1:N time = (k-1) * Ts; % 工况: 0.02s时负载从4Ohm突变到2Ohm % 0.035s时输入电压从24V突变到20V if time < 0.02 R_load = R_nom; else R_load = 2.0; end if time > 0.035 V_in_k = V_in - 4; else V_in_k = V_in; end % 电压变化率(实际工程用滤波微分, 这里取理想导数演示算法) dvC_dt = (iL - vC / R_load) / C; dedt = -dvC_dt; % 误差导数 e = V_ref - vC; % 输出电压误差

5.2 双观测器更新

先更新负载观测器，把R_hat送到控制器模型里；再更新NDOB，估计集总扰动的残余部分。注意这里的f0_hat要和第2章的误差方程严格对应，符号错了，整个观测器会失效。

%% 负载观测器: 用物理关系估计负载电流和电阻 i_R_est = iL - C * dvC_dt; % 估计负载电流 R_raw = vC / max(i_R_est, 1e-3); % 估计负载电阻, 加分母保护 R_raw = min(max(R_raw, 0.5), 50); % 限幅避免极端值 R_hat = (1 - alpha_R) * R_hat + alpha_R * R_raw; % 低通平滑 %% 非线性干扰观测器 NDOB % 误差方程: dedt_dot = f0 + b0*u + d_total f0_hat = vC/(L*C) - dedt/(R_hat*C); b0 = V_in/(L*C); p = l_ndob * dedt; z_d = z_d - l_ndob * (z_d + p) * Ts - l_ndob * (f0_hat + b0 * u_old) * Ts; d_hat = z_d + p;

5.3 超螺旋控制律计算

先计算滑模面s，然后离散化超螺旋积分项xi，最后合成控制占空比u。占空比必须限幅，这是所有数字电源的底线。

%% 超螺旋高阶滑模控制律 s = dedt + c * e; % 滑模面 xi = xi - W_h * sign(s) * Ts; % 积分项(不连续项先到导数里) w_htc = -lambda_w * sqrt(abs(s)) * sign(s) + xi; % 超螺旋控制量 % 控制律: 模型前馈 + 双观测器补偿 u = vC/V_in + (L/V_in) * (c - 1/(R_hat*C)) * dedt ... + (L*C/V_in) * (d_hat - w_htc); % 占空比限幅 u = min(max(u, 0.05), 0.95); u_old = u; % 存给下一个周期的NDOB使用

5.4 状态更新与结果保存

平均模型用中点法更新，比向前欧拉的数值精度好不少，跑60ms仿真几乎不担心稳定性问题。

%% 状态更新(中点法) diL_dt = (V_in_k * u - vC) / L; dvc_dt = (iL - vC / R_load) / C; % 中点预测 iL_mid = iL + 0.5 * Ts * diL_dt; vC_mid = vC + 0.5 * Ts * dvc_dt; diL_mid = (V_in_k * u - vC_mid) / L; dvc_mid = (iL_mid - vC_mid / R_load) / C; % 更新 iL = iL + Ts * diL_mid; vC = vC + Ts * dvc_mid; %% 保存数据 u_log(k) = u; vC_log(k) = vC; e_log(k) = e; dhat_log(k) = d_hat; Rhat_log(k) = R_hat; end %% 绘图 figure; subplot(3,1,1); plot(t, vC_log, 'LineWidth', 1.2); hold on; yline(V_ref, 'r--', 'LineWidth', 1); xlabel('Time (s)'); ylabel('v_C (V)'); title('Buck输出电压'); legend('v_C','V_{ref}','Location','best'); grid on; subplot(3,1,2); plot(t, u_log, 'LineWidth', 1.2); xlabel('Time (s)'); ylabel('占空比'); title('PWM占空比'); grid on; subplot(3,1,3); plot(t, Rhat_log, 'LineWidth', 1.2); hold on; plot([0 0.02], [4 4], 'k--'); plot([0.02 0.06], [2 2], 'k--'); xlabel('Time (s)'); ylabel('R_{hat} (Ohm)'); title('负载观测器估计值'); legend('R_{hat}', '真实R', 'Location','best'); grid on;

想要看NDOB效果，可以再画一张dhat_log曲线；想要做对比实验，就把控制律里的d_hat强制置零、R_hat强制回到标称值，电压跌落幅度立刻变大，这套方案的差异一目了然。

6. 参数整定顺序与一组可复现的仿真结果

6.1 核心增益的整定顺序

参数整定最重要的一条原则：先让系统稳定，再追求抗扰性能。我推荐这个顺序：

第一步，整定滑模面斜率c。c决定了误差收敛速率，量级上取开关角频率的十分之一到五分之一。本文20kHz开关频率，c在2000~6000都合理。c太小，恢复慢；c太大，数字延迟容易激发出振荡。

第二步，固定W_h，从小到大调lambda_w。lambda_w是超螺旋比例项，调大能让误差快速归零。但过大会导致占空比在瞬态时打得太满，限幅频繁触发。

第三步，再微调W_h。W_h的作用是保证扰动下s不会“滑出”滑模面。经验值从lambda_w的十分之一起步，逐步加到占空比波形出现轻微低频波动为止，然后回退20%。

第四步，调NDOB增益l_ndob。我习惯让l_ndob比控制环带宽高3~5倍。本文控制环带宽约1kHz，l_ndob取2000左右就够。再大会让观测器把采样噪声也当成扰动估计出来，占空比波形会变毛糙。

第五步，调负载观测器滤波系数alpha_R。alpha_R越小越平滑，但对负载突变的跟踪越慢。0.05在大多数Bluck负载突变场景下兼顾了两头。

以本文参数为例，我跑出来的典型结果是：0.02s负载从4Ω切到2Ω，输出电压最大跌落约60mV，恢复时间约1.5ms；0.035s输入电压从24V突变到20V，输出电压跌落更小，约30mV以内。如果把两个观测器都关掉，同样参数下负载突变跌落会达到200mV以上，恢复时间翻倍还不止。这说明观测器补偿不是“优化锦上添花”，而是这套方案能不能压住动态的胜负手。

6.2 启动过程怎么处理

上面的仿真初值我直接设在了稳态点，是为了验证抗扰性能。如果你要跑启动过程，有两个办法：一是给参考电压加斜率限制，让V_ref在5ms内从0斜坡上升到12V；二是限制占空比变化率。直接让零状态启动，输出电压误差很大，滑模控制很快就打到占空比上限，和开环硬充电没有本质区别，观测器前几毫秒还会被带偏，没必要一上来就给自己上难度。

7. 工程落地的几个教训：采样、限幅与噪声处理

7.1 误差导数不要直接差分

仿真里我可以取理想导数，但样机上直接对v_C差分，结果就是一场灾难。输出电压纹波通常是几十毫伏，控制周期50微秒，差分出来的导数噪声能达到真实值的几倍。我实际用的是“低通滤波+差分”的组合：先对v_C做截止频率3~5kHz的一阶低通，再差分会好很多；进阶方案是用一个滑模微分器或者线性扩张状态观测器同时估计导数，效果更干净。代价是微分器本身也有带宽限制，相位滞后会影响c的上限。

7.2 数字控制器的执行顺序必须固定

采样、观测器、控制器、PWM更新，这个顺序不能乱。我吃过一次亏：先算控制律再更新观测器，结果观测器用的还是“未来”的控制量，一个采样周期的混乱导致突发振荡。正确顺序是：ADC采样iL和vC，更新负载观测器和NDOB，计算滑模控制律，更新占空比寄存器。NDOB里用的u应该是“上一周期实际施加的占空比”，而不是本周期的，代码里我已经用u_old做了这个错拍。

7.3 限幅之外，还要给观测器加“安全带”

占空比限幅是底线，但观测器也有自己的安全隐患。负载观测器的R_hat如果不受限幅，短路工况下会估计出0.1Ω甚至负值，后面NDOB和控制律全乱。代码里min(max(R_raw,0.5),50)这个限幅不是随便写的，0.5Ω对应满负载电流极限，50Ω对应接近空载。NDOB的d_hat同样要限幅，我一般限制在额定占空比修正量的正负30%以内。一旦超过，优先怀疑工况异常而不是真扰动，先把输出稳住再说。

7.4 死区和延迟会让高增益失效，别盲目加l_ndob

数字电源里PWM死区、开关延迟、ADC采样延迟，本质上都是高频不匹配项。NDOB增益l_ndob调高，看似观测误差收敛快了，但延迟会让观测器在扰动到来时“冲过头”，反而在占空比里引入振荡。我的实测经验是：l_ndob超过4000以后，抗扰提升微乎其微，而电感啸叫的几率明显上升。如果你的系统有严重的延迟，用Smith预估器先补偿延迟，再提高观测带宽，这才是正路。

最后再说一条我踩出来的经验：高阶滑模和观测器这套组合，最大的价值不是某个单一性能指标的压倒性优势，而是把“鲁棒性”和“低抖振”这两个原本对立的要求同时满足了。调试的时候，不要一上来就追求电压跌落小于30mV，先把占空比波形调到“稳态安静、瞬态有力”，再去调观测器带宽，你会发现一切都水到渠成。

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