s函数报错

勇往直前11

s函数报错“during flag=i call must be a real vector of length n” ，怎么解决？
代码如下：
function [sys,x0,str,ts] = PMSM_chaos_controller(t,x,u,flag)
1 F- m% {+ C$ U. m$ V! E2 {' Hswitch flag,
' x3 k- U! M- ?' x. m. \( v) r    case 0,
        [sys,x0,str,ts] = mdlIniticializeSizes;
    case 3,
        sys = mdlOutputs(t,x,u);
    case {1,2,4,9},
        sys = [];
6 B. j; l$ w* [9 N' s    otherwise,
        DAStudio.error('Simulink:blocks:unhandledFlag',num2str(flag));
; o; _; N+ v" n- Y& c/ m  Wend
) N$ L6 v/ R6 X9 l: G$ J0 l4 u
/ ~. S. h, @1 h) w- \function [sys,x0,str,ts] = mdlIniticializeSizes
sizes = simsizes;
' j3 a$ }& W/ i4 B- Q4 ssizes.NumContStates  = 3;    %连续状态个数
sizes.NumDiscStates  = 0;    %离散状态个数
sizes.NumOutputs     = 3;    %输出变量个数
sizes.NumInputs      = 9;    %输入信号个数
sizes.DiRFeedthrough = 1;    %输入直接传入输出信号否
sizes.NumSampleTimes = 0;    % at least one sample time is neededWO  一般来说为1个
sys = simsizes(sizes);
( E5 F3 S# [( ]' gx0 = [0.001 0.001 0.001];
str = [];
1 o) x1 U; E" ets = [];

# U2 F- ?8 Q2 S# l: l4 P+ Gfunction sys=mdlOutputs(t,x,u)
0 X1 @% S/ l6 ~/ G" Ue1=u(1);e2=u(2);e3=u(3);
" c( X8 P; p; N9 j% ks21=u(4);s22=u(5);s1=u(6);
8 ^# }$ z: F0 J4 d+ N; Mes_w=u(9);
alpha = 0.8;
alpha1 = 5/7; beta1 = 7/5;
6 P# y$ C) s9 o+ @7 Xalpha2 = 7/9; beta2 = 9/7;
+ D' V" T% c7 X9 W, E2 n6 q: \k1 = 1;
k2 = 1.2;
. T# P; a0 b( P/ c" Z0 ~* Zmiu = 10;
delta = 5.46;
- B' Q& l2 w. j3 `%assert(all(imag(u)==0), 'u is imaginary or nan');不起作用
if e3==0
2 B8 V- j! y/ Q" P    u1 = e1 - es_w*e1 - 10*((abs(e1))^alpha)*sign(e1) - 10*e1 - k2*((abs(s21))^alpha2)*sign(s21) -k2*((abs(s21))^beta2)*sign(s21) - miu*s21;
    u2 = e2 - es_w*e2 - 10*((abs(e2))^alpha)*sign(e2) - 10*e2 - k2*((abs(s22))^alpha2)*sign(s22) -k2*((abs(s22))^beta2)*sign(s22) - miu*s22;
    u3 = -delta*(e2-e3) - 10*((abs(e3))^alpha)*sign(e3) - 10*e3 - k1*((abs(s1))^alpha1)*sign(s1) -k1*((abs(s1))^beta1)*sign(s1) - miu*s1;
7 v2 F8 d  X  F8 q: xelse
  N" k  r+ F) N  e+ p4 `1 r    u1 = 0;
: d/ W5 ]& H3 b* n; E: W& F6 m( H    u2 = 0;
! \3 X; b$ y. y" d9 \' h    u3 = -delta*(e2-e3) - 10*((abs(e3))^alpha)*sign(e3) - 10*e3 - k1*((abs(s1))^alpha1)*sign(s1) -k1*((abs(s1))^beta1)*sign(s1) - miu*s1;
& Z5 }* j  S; h1 x; Fend
. l) S+ I2 L3 |1 Z( Q6 t+ V; v, Bsys(1) = u1;
sys(2) = u2;
sys(3) = u3;
) k2 ]+ h2 I$ Z$ @复制代码

Uifhjvv

帮你顶一下

小小鲁班

# C9 q4 o9 c# j" S你这个没有用到状态变量吧

zaiyiaaaa

你上传一下simulink仿真图？

zzz.dan