用MATLAB中S函数仿真模型参考自适应遇到的问题

maggie19811117

我在做感应电机参数辨识,采取模型参考自适应的方法,采用参数误差逼近的方法.用S函数编写控制模型时,遇到一些问题,我的程序如下:
function [sys,x0,str,ts] = ac_model(t,x,u,flag)

switch flag,
   
  case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
   
  case 1,
    sys=mdlDerivatives(t,x,u);
   
case 3,
    sys=mdlOutputs(t,x,u);
   
case{2,4,9}
    sys=[];

  otherwise
    error(['Unhandled flag = ',num2str(flag)]);

end

function [sys,x0,str,ts]=mdlInitializeSizes(x0);
sizes = simsizes;

sizes.NumContStates  = 7;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 7;
sizes.NumInputs      = 10;
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1;  

sys = simsizes(sizes);
x0  = [0;0;0;0;0;0;20];
str = [];
ts  = [0 0];

function sys=mdlDerivatives(t,x,u)
g1=100;g2=500;g3=500;g4=10;g5=10;k=100;
dx(1)=x(3)*u(1)-u(5)*u(2)+x(4)*u(7)+x(5)*u(9)+x(6)*u(5)*u(8)+x(7)*(u(5)*(u(10)+u(3)))-k*u(6);
dx(2)=x(3)*u(2)+u(5)*u(1)+x(4)*u(7)+x(5)*u(10)-x(6)*u(5)*u(7)+x(7)*(u(4)-u(5)*u(9));
dx(3)=-g1*u(6)*u(1);
dx(4)=-g2*u(6)*u(7);
dx(5)=-g3*u(6)*u(9);
dx(6)=-g4*u(6)*u(5)*u(8);
dx(7)=-g5*u(6)*(u(5)*u(10)+u(3));

sys=dx;

function sys=mdlOutput(t,x,u)

tspan=[t0;tf];
options=odeset('RelTol',1e-3,'AbsTol',1e-3);
[t,x]=ode45('ac_model',tspan,x0,options);

if abs(u(6))>1e-8
    sys=[];
else
   sys = [x(1);x(2);x(3);x(4);x(5);x(6);x(7)];
end

老是出现 flag = 3 (output), at time 0. MATLAB error message:
Undefined function or variable 'mdlOutputs'为什么呢?
望不舍赐教,谢谢

cdwxg

你定义的是mdlOutput
但你上面用的是mdlOutputs
如红笔所示
function [sys,x0,str,ts] = ac_model(t,x,u,flag)

switch flag,
   
  case 0,
    [sys,x0,str,ts]=mdlInitializeSizes;
   
  case 1,
    sys=mdlDerivatives(t,x,u);
   
case 3,
    sys=mdlOutputs(t,x,u);
   
case{2,4,9}
    sys=[];

  otherwise
    error(['Unhandled flag = ',num2str(flag)]);

end

function [sys,x0,str,ts]=mdlInitializeSizes(x0);
sizes = simsizes;

sizes.NumContStates  = 7;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 7;
sizes.NumInputs      = 10;
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1;  

sys = simsizes(sizes);
x0  = [0;0;0;0;0;0;20];
str = [];
ts  = [0 0];

function sys=mdlDerivatives(t,x,u)
g1=100;g2=500;g3=500;g4=10;g5=10;k=100;
dx(1)=x(3)*u(1)-u(5)*u(2)+x(4)*u(7)+x(5)*u(9)+x(6)*u(5)*u(8)+x(7)*(u(5)*(u(10)+u(3)))-k*u(6);
dx(2)=x(3)*u(2)+u(5)*u(1)+x(4)*u(7)+x(5)*u(10)-x(6)*u(5)*u(7)+x(7)*(u(4)-u(5)*u(9));
dx(3)=-g1*u(6)*u(1);
dx(4)=-g2*u(6)*u(7);
dx(5)=-g3*u(6)*u(9);
dx(6)=-g4*u(6)*u(5)*u(8);
dx(7)=-g5*u(6)*(u(5)*u(10)+u(3));

sys=dx;

function sys=mdlOutput(t,x,u)

tspan=[t0;tf];
options=odeset('RelTol',1e-3,'AbsTol',1e-3);
[t,x]=ode45('ac_model',tspan,x0,options);

if abs(u(6))>1e-8
    sys=[];
else
   sys = [x(1);x(2);x(3);x(4);x(5);x(6);x(7)];
end

xiAozhI407

细心呀～～