电液伺服系统ADRC控制：Simulink实现与参数整定

1. 电液伺服控制中的ADRC实战：从原理到Simulink实现

电液伺服系统在工业自动化领域扮演着重要角色，但其强非线性、参数时变和未知干扰等特性，使得传统PID控制器往往难以取得理想的控制效果。自抗扰控制(ADRC)技术因其独特的扰动处理能力，成为解决这类问题的利器。本文将带您深入理解ADRC的核心原理，并手把手教您在Simulink中搭建可自由切换版本的三阶ADRC控制器。

提示：本文所有代码和模型均基于MATLAB R2021a开发，但设计思路适用于各版本Simulink环境

1.1 为什么ADRC适合电液伺服系统？

电液伺服系统的位移跟踪控制面临三大挑战：

1. 非线性特性：液压油的压缩性、阀口的流量-压力特性等导致系统呈现显著非线性
2. 参数时变：油温变化、密封件磨损等使系统参数随时间漂移
3. 未知扰动：负载突变、外部力干扰等不可预测因素

传统PID控制在这些挑战面前表现不佳，主要是因为：

• 线性组合的PID难以应对非线性系统
• 固定参数无法适应时变特性
• 扰动出现后才开始反应，存在滞后

ADRC的核心优势在于其"总扰动"概念和实时估计补偿机制，通过扩张状态观测器(ESO)将所有不确定性和扰动视为一个总扰动项进行实时估计和补偿，相当于为系统安装了一个"智能减震器"。

2. ADRC核心模块深度解析

2.1 三阶ADRC架构设计

标准三阶ADRC包含三个关键子系统：

2.1.1 跟踪微分器(TD)

matlab复制function [v1, v2] = TD(v, h, r)
    persistent v1_prev v2_prev
    if isempty(v1_prev)
        v1_prev = 0;
        v2_prev = 0;
    end
    fh = fhan(v1_prev-v, v2_prev, r, h);
    v1 = v1_prev + h*v2_prev;
    v2 = v2_prev + h*fh;
    v1_prev = v1;
    v2_prev = v2;
end

TD的作用是对给定信号进行平滑处理并提取微分信号，其核心参数r决定跟踪速度。对于电液伺服系统，建议初始值设为系统响应时间的倒数。

2.1.2 扩张状态观测器(ESO)

ESO是ADRC的核心创新，其Matlab实现如下：

matlab复制function [z1, z2, z3] = ESO(u, y, h, beta1, beta2, beta3, b0)
    persistent z1_prev z2_prev z3_prev
    if isempty(z1_prev)
        z1_prev = y;
        z2_prev = 0;
        z3_prev = 0;
    end
    e = z1_prev - y;
    dz1 = z2_prev - beta1*e;
    dz2 = z3_prev - beta2*fal(e,0.5,h) + b0*u;
    dz3 = -beta3*fal(e,0.25,h);
    % 状态更新
    z1 = z1_prev + h*dz1;
    z2 = z2_prev + h*dz2;
    z3 = z3_prev + h*dz3;
    % 保持状态
    z1_prev = z1;
    z2_prev = z2;
    z3_prev = z3;
end

关键点说明：

• z3是扩张状态，专门用于估计总扰动
• fal()函数实现非线性跟踪，比线性函数具有更好的鲁棒性
• b0是系统先验知识补偿项，即使设为0也能工作，但适当设置可提高收敛速度

2.1.3 非线性状态误差反馈(NLSEF)

matlab复制function u = NLSEF(e1, e2, alpha, delta)
    u = -fal(e1,alpha(1),delta) - fal(e2,alpha(2),delta);
end

function f = fal(e,alpha,delta)
    if abs(e) > delta
        f = abs(e)^alpha * sign(e);
    else
        f = e / (delta^(1-alpha));
    end
end

NLSEF采用非线性组合取代PID的线性组合，参数α决定非线性强度，δ为线性区间阈值。

2.2 参数整定方法论

ADRC参数整定可遵循"带宽参数化"原则：

现场调试技巧：

1. 先调ESO参数确保扰动估计准确
2. 再调NLSEF参数优化跟踪性能
3. 最后微调TD使指令过渡平滑

注意：调试时应先施加小幅度阶跃信号，观察各状态变量是否正常收敛后再逐步增大输入幅度

3. Simulink实现与版本切换

3.1 模块化设计架构

在Simulink中，我们采用分层封装的设计方法：

code复制ADRC_Controller (顶层)
├── TD_Subsystem
├── ESO_Subsystem
├── NLSEF_Subsystem
└── Version_Switch

3.1.1 版本切换实现方案

通过Model Variants实现多版本共存：

1. 创建Controller_Lib.slx库文件
2. 为每个版本创建子系统变体
3. 使用Switch模块和回调函数实现运行时切换

切换逻辑代码示例：

matlab复制function switch_ADRC_version(version)
    blk = gcb;
    switch version
        case 'Standard'
            set_param([blk '/Controller'], 'Variant', 'Standard_ADRC');
            load('ADRC_std_params.mat');
        case 'Reduced_Order'
            set_param([blk '/Controller'], 'Variant', 'Reduced_ADRC');
            load('ADRC_red_params.mat');
        case 'Cascade'
            set_param([blk '/Controller'], 'Variant', 'Cascade_ADRC');
            load('ADRC_cas_params.mat');
    end
end

3.2 各版本特点对比

4. 实战调试与性能分析

4.1 典型问题排查指南

4.2 实测性能对比

在电液伺服试验台上进行阶跃响应测试：

调试心得：

1. ESO的z3输出是系统健康的"晴雨表"，调试时应重点监控
2. 采样周期h的选择很关键，一般取系统响应时间的1/10~1/20
3. 实际系统中，ADRC的b0参数可以离线辨识获得，能显著提升初始性能

5. 高级技巧与扩展应用

5.1 调试模式开发

在模型中添加Debug开关：

matlab复制if strcmp(get_param('ADRC_Controller/Debug','sw'), 'on')
    disturbance_comp = z3/b0;
    assignin('base','disturbance_comp',disturbance_comp);
end

这样可以在不修改主代码的情况下观察扰动补偿量。

5.2 离散化实现要点

对于实际数字控制器，需注意：

1. 使用Tustin变换保持非线性特性
2. 添加抗积分饱和逻辑
3. 实现输出限幅保护

离散ESO示例：

matlab复制function [z1, z2, z3] = discrete_ESO(u, y, Ts, beta1, beta2, beta3, b0)
    persistent z1 z2 z3
    if isempty(z1)
        z1 = y;
        z2 = 0;
        z3 = 0;
    end
    e = z1 - y;
    z1 = z1 + Ts*(z2 - beta1*e);
    z2 = z2 + Ts*(z3 - beta2*fal(e,0.5,Ts) + b0*u);
    z3 = z3 - Ts*beta3*fal(e,0.25,Ts);
end

5.3 自适应ADRC进阶

对于参数变化剧烈的系统，可以实现参数在线调整：

matlab复制function adjust_params()
    persistent error_integral
    if isempty(error_integral)
        error_integral = 0;
    end
    error = ref - actual;
    error_integral = error_integral + error*Ts;
    if abs(error_integral) > threshold
        beta1 = beta1 * (1 + sign(error_integral)*0.05);
        % 其他参数调整类似
        error_integral = 0;
    end
end

在实际液压伺服系统调试中，我发现ADRC对采样噪声比较敏感，建议在ESO输入端添加适当的低通滤波。一个实用的技巧是将滤波时间常数设为采样周期的2-3倍，这样既能有效抑制噪声，又不会引入明显相位滞后。