s函数报错

勇往直前11

s函数报错“during flag=i call must be a real vector of length n” ，怎么解决？
代码如下：
+ f' e8 d, |$ X$ u; Z- yfunction [sys,x0,str,ts] = PMSM_chaos_controller(t,x,u,flag)
' [: q% F4 N: ]( K/ p+ d8 `+ |switch flag,
    case 0,
8 ~3 W& F+ q3 _% C  x        [sys,x0,str,ts] = mdlIniticializeSizes;
    case 3,
' ]' {) c, S* Y* N* Y  a        sys = mdlOutputs(t,x,u);
! P; T( \& w! C0 \3 q+ k4 u    case {1,2,4,9},
        sys = [];
7 e) e! e: l7 j2 M; x& \    otherwise,
        DAStudio.error('Simulink:blocks:unhandledFlag',num2str(flag));
' n( G3 _' r# g2 N8 o  z8 n2 nend
5 v( I* e/ a; S
9 e8 u) S) N& v3 t, rfunction [sys,x0,str,ts] = mdlIniticializeSizes
0 d) h4 I( D! W/ w$ t8 Isizes = simsizes;
sizes.NumContStates  = 3;    %连续状态个数
sizes.NumDiscStates  = 0;    %离散状态个数
3 y+ o3 \" L* q9 s( d+ V& Hsizes.NumOutputs     = 3;    %输出变量个数
sizes.NumInputs      = 9;    %输入信号个数
sizes.DiRFeedthrough = 1;    %输入直接传入输出信号否
9 ?$ h$ Y& |: d3 y2 k$ Esizes.NumSampleTimes = 0;    % at least one sample time is neededWO  一般来说为1个
; Y5 e# W$ @0 y! D" K4 E( {sys = simsizes(sizes);
x0 = [0.001 0.001 0.001];
str = [];
6 o- E3 n+ D% T+ ^; y& J* h9 Uts = [];
. I+ F% N8 t  a% x
function sys=mdlOutputs(t,x,u)
e1=u(1);e2=u(2);e3=u(3);
s21=u(4);s22=u(5);s1=u(6);
% Z+ U: S5 i5 X& t( Ves_w=u(9);
4 n8 A1 J! `3 ?alpha = 0.8;
& r, e; r  S8 h  valpha1 = 5/7; beta1 = 7/5;
alpha2 = 7/9; beta2 = 9/7;
k1 = 1;
k2 = 1.2;
miu = 10;
6 J  W. Z" j; t" s7 v% Ndelta = 5.46;
' Z% b: n2 U" q! R4 M/ [; R! Q%assert(all(imag(u)==0), 'u is imaginary or nan');不起作用
if e3==0
9 Z- K$ _* Y( E" p8 \$ I    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;
& w$ D) G$ l* {0 @& \    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;
) s- l9 N2 ?- }2 }else
5 X9 p6 {# j- {% B0 b# Y    u1 = 0;
4 H+ L5 b0 y" O5 ^1 ~    u2 = 0;
    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;
end
sys(1) = u1;
& F2 _: ^% ^- B7 m  J; Psys(2) = u2;
  s4 a' w, s! rsys(3) = u3;
4 r4 Q) s) K3 R% c8 x2 Q复制代码

Uifhjvv

帮你顶一下

小小鲁班

  \' T" R3 j) t) _7 }- o你这个没有用到状态变量吧

zaiyiaaaa

你上传一下simulink仿真图？

zzz.dan