永磁同步电机自抗扰无位置传感器控制仿真，同时实现自抗扰和基于eso扩张状态观测器的无位置控制仿真。

咱们搞电机控制的都清楚，无位置传感器控制就像在玩"盲人摸象"，特别是当系统存在未知扰动时更是头疼。最近折腾了半个月的ADRC+ESO组合方案，实测效果确实有点东西，今天就把仿真过程和关键代码掰开了说说。

先说核心思路：用ESO实时观测系统总扰动（包括外部干扰和模型误差），同时在速度环用自抗扰控制器替代传统PI。重点在于让观测器能同时估计转速和位置，这就得搞点数学魔法了。

先看滑模观测器的关键实现代码片段：

function [theta_est, omega_est] = SMO_Observer(v_alpha, v_beta, i_alpha, i_beta)
    persistent z_alpha z_beta;
    if isempty(z_alpha)
        z_alpha = 0; z_beta = 0;
    end
    
    L = 0.05;  % 观测器增益
    k = 100;    % 滑模增益
    Ts = 1e-4;  % 采样时间
    
    e_alpha = i_alpha - z_alpha;
    e_beta = i_beta - z_beta;
    
    s_alpha = sign(e_alpha);
    s_beta = sign(e_beta);
    
    dz_alpha = (v_alpha - R*i_alpha)/L + k*s_alpha;
    dz_beta = (v_beta - R*i_beta)/L + k*s_beta;
    
    z_alpha = z_alpha + dz_alpha*Ts;
    z_beta = z_beta + dz_beta*Ts;
    
    theta_est = atan2(z_beta, z_alpha);
    omega_est = (diff(theta_est))/Ts;
end

这段代码实现了滑模观测器的核心算法，重点注意sign函数带来的抖振问题。实际调试时发现把sign换成饱和函数sat(x)=x/(|x|+δ)效果更好，δ取0.01左右能明显平滑波形。

接下来是ADRC控制器的灵魂——ESO部分：

function [omega_hat, disturbance] = ESO_Update(i_q, u_q)
    persistent z1 z2;
    if isempty(z1)
        z1 = 0; z2 = 0;
    end
    
    beta1 = 100;  % 带宽参数
    beta2 = 3000;
    Ts = 1e-4;
    
    e = z1 - i_q;
    
    dz1 = z2 - beta1*e + b*u_q;
    dz2 = -beta2*e;
    
    z1 = z1 + dz1*Ts;
    z2 = z2 + dz2*Ts;
    
    omega_hat = z1;
    disturbance = z2/b;  % 扰动估计
end

这个三阶ESO的设计精髓在于带宽参数的选取。beta1控制跟踪速度，beta2决定抗扰能力，但过大会引入噪声。经验值是让beta2≈(3~5)*beta1^2，调试时用阶跃响应观察估计曲线是否紧贴真实值。

仿真结果对比传统PI控制，在突加负载时转速波动从±15rpm降到±3rpm，位置估计误差在1000rpm时小于0.5机械角度。不过要注意电机参数敏感性，特别是定子电阻变化超过20%时需要重新整定ESO参数。

最后给个调参小技巧：先关闭ADRC只用ESO做开环观测，把估计波形调准后再接入闭环。遇到高频振荡别急着降带宽，先检查PWM频率与采样周期是否匹配，很多时候问题出在离散化方法选错（推荐用Tustin变换）。

这套方案在医疗离心机控制中实测效果拔群，下次可以聊聊如何移植到DSP时避免浮点运算爆内存的骚操作。搞控制就像炒菜，火候到了自然香，参数整定急不得。