s函数报错

勇往直前11

s函数报错“during flag=i call must be a real vector of length n” ，怎么解决？
代码如下：
( K  Z+ C4 T/ x5 @4 X4 Sfunction [sys,x0,str,ts] = PMSM_chaos_controller(t,x,u,flag)
; J( j) k# v% ?0 Yswitch flag,
    case 0,
        [sys,x0,str,ts] = mdlIniticializeSizes;
    case 3,
        sys = mdlOutputs(t,x,u);
7 ]5 N8 L# l( i% ]    case {1,2,4,9},
        sys = [];
, b# I2 B0 l" a) t* d% w: [+ W$ x7 D    otherwise,
: E0 h, l! I6 W6 P1 d+ w  [6 y7 w        DAStudio.error('Simulink:blocks:unhandledFlag',num2str(flag));
end

function [sys,x0,str,ts] = mdlIniticializeSizes
sizes = simsizes;
8 K- w, A0 u- |  x: Osizes.NumContStates  = 3;    %连续状态个数
% o& `. J, [. Gsizes.NumDiscStates  = 0;    %离散状态个数
' Q5 V5 X0 \! o9 i3 t0 f) xsizes.NumOutputs     = 3;    %输出变量个数
sizes.NumInputs      = 9;    %输入信号个数
sizes.DiRFeedthrough = 1;    %输入直接传入输出信号否
7 }* S( @( t  D9 ~+ Dsizes.NumSampleTimes = 0;    % at least one sample time is neededWO  一般来说为1个
sys = simsizes(sizes);
) U- ?) w' o1 E% m) p$ l; Zx0 = [0.001 0.001 0.001];
str = [];
+ F7 c  }4 Y" {4 R: z# m! cts = [];

/ x, S3 q$ t: Gfunction sys=mdlOutputs(t,x,u)
0 Y$ G; Y/ B2 u1 W& p& Ze1=u(1);e2=u(2);e3=u(3);
. t8 U0 N1 D! `; ss21=u(4);s22=u(5);s1=u(6);
es_w=u(9);
, k) A( N' Q) @- m; @: g) r- ?: ~$ Palpha = 0.8;
alpha1 = 5/7; beta1 = 7/5;
alpha2 = 7/9; beta2 = 9/7;
6 B2 V- ~- R$ J- ~- b- m6 R8 \k1 = 1;
k2 = 1.2;
miu = 10;
/ s; Y, F+ W4 o( G+ vdelta = 5.46;
( t# ?$ t8 _$ s7 ]  i: D! J%assert(all(imag(u)==0), 'u is imaginary or nan');不起作用
if e3==0
    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;
6 @% ]2 e) V  d) J    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;
& T$ E& V% R* j. Q: u1 E/ {: h! L    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;
else
  A9 T- }3 I  A( X' g    u1 = 0;
    u2 = 0;
4 M+ ]' \7 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;
end
9 D; Q9 z8 a2 Q: gsys(1) = u1;
$ S$ P3 x: m$ Q7 K, H" wsys(2) = u2;
sys(3) = u3;
复制代码

Uifhjvv

帮你顶一下

小小鲁班

. s% h7 d4 w) o; H你这个没有用到状态变量吧

zaiyiaaaa

你上传一下simulink仿真图？

zzz.dan