自抗扰控制算法

这个一个能用的自抗扰程序，初学者好用，里面有微分器，观测器，非线性反馈模块！

function [sysx0strts]=eso（txuflag）

switch flag
  case 0                                                
    [sysx0strts] = mdlInitializeSizes;
  case 2                                                
    sys = mdlUpdates（xu）;
   case 3                                                
    sys = mdlOutputs（x）;  
  case {149}                                         
    sys = []; % do nothing    
  otherwise 
     error（[‘Unhandled flag=‘num2str（flag）]）; 
end
function [sysx0strts]=mdlInitializeSizes
    sizes=simsizes;
    sizes.NumContStates=0;
    sizes.NumDiscStates=3;
    sizes.NumOutputs=3;
    sizes.NumInputs=2;
    sizes.DirFeedthrough=0;
    sizes.NumSampleTimes=1;
    sys=simsizes（sizes）;
    x0=[0;0;0];
    str=[];
    ts=[-1 0];
function sys=mdlUpdates（xu）
bet=[19  2 1.5];b=0.1;h=0.1;d=9*h;% d取大一点的值减小了波动幅度 bet的第二个参数调小，波动减小
    e=x（1）-u（2）;
    sys（11）=x（1）+h*（x（2）-bet（1）*e）;
    sys（21）=x（2）+h*（x（3）-bet（2）*fal（e0.5d）+b*u（1））;
    sys（31）=x（3）-h*bet（3）*fal（e0.25d）;
function sys=mdlOutputs（x）
    sys=x;
function f=fal（ead）
    if abs（e）<=d
        f=e*d^（a-1）;
    else f=（abs（e））^a*sign（e）;
    end
    
        
    
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1240  2010-11-21 01:25  自抗扰控制（改为h0）\d_eso50.m

     文件       1207  2010-11-21 01:25  自抗扰控制（改为h0）\d_nlsef50.m

     文件       1371  2015-08-21 07:38  自抗扰控制（改为h0）\d_td50.m

     文件        143  2010-11-21 01:25  自抗扰控制（改为h0）\fal.m

     文件       3411  2010-11-21 01:25  自抗扰控制（改为h0）\model_zhengqifashengqi.m

     文件      53642  2010-11-21 01:25  自抗扰控制（改为h0）\zikangrao.mdl

     文件      32178  2010-11-21 01:25  自抗扰控制（改为h0）\蒸汽发生器模型参数.jpg

     目录          0  2018-02-02 20:42  自抗扰控制（改为h0）

----------- ---------  ---------- -----  ----

                93192                    8