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)
switch flag,
    case 0,
        [sys,x0,str,ts] = mdlIniticializeSizes;
    case 3,
, n# c1 a* A+ ?7 c% k% p2 O        sys = mdlOutputs(t,x,u);
    case {1,2,4,9},
+ K0 [" n3 f! H        sys = [];
    otherwise,
: `$ V9 \4 b4 C+ q8 x        DAStudio.error('Simulink:blocks:unhandledFlag',num2str(flag));
- {# h  _/ `$ z9 `  t: ^end
2 \8 @6 U' Z! B2 z
6 \/ v; J4 n5 Ufunction [sys,x0,str,ts] = mdlIniticializeSizes
. z9 H+ `/ s' G% a: ]sizes = simsizes;
4 M( j7 ?! n( r+ o' Jsizes.NumContStates  = 3;    %连续状态个数
sizes.NumDiscStates  = 0;    %离散状态个数
sizes.NumOutputs     = 3;    %输出变量个数
* I4 A) s$ _+ |( A1 esizes.NumInputs      = 9;    %输入信号个数
  V$ _: t( D3 r, s1 gsizes.DiRFeedthrough = 1;    %输入直接传入输出信号否
* y! i) t8 `) \! xsizes.NumSampleTimes = 0;    % at least one sample time is neededWO  一般来说为1个
6 a0 |; A, ~  I9 e; M3 p: k7 K7 ?sys = simsizes(sizes);
x0 = [0.001 0.001 0.001];
str = [];
ts = [];
5 r# T" D! M% o
function sys=mdlOutputs(t,x,u)
& Y* {( |) d8 y6 }5 i; Le1=u(1);e2=u(2);e3=u(3);
s21=u(4);s22=u(5);s1=u(6);
es_w=u(9);
alpha = 0.8;
! a) V6 b$ T6 S8 F# B4 balpha1 = 5/7; beta1 = 7/5;
5 u! h) b% B5 x6 {9 X1 _0 ealpha2 = 7/9; beta2 = 9/7;
k1 = 1;
k2 = 1.2;
' B; k9 ]  C% G3 c% j  L- q3 p) A" `miu = 10;
4 H8 ^( e* Q+ O" s9 ]# cdelta = 5.46;
%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;
: i# P% W: \7 R! H1 v    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;
. S; u) G! ]4 p! X0 i+ y& d4 q    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;
$ ^2 y* J, Z' relse
9 g  v$ H* G( R% R4 d+ t2 T! X    u1 = 0;
% R! D+ o7 a1 D6 S0 e2 @; X    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;
) a6 w1 d# o3 F+ |0 |! o* b4 esys(2) = u2;
. `9 n% N) ]* }* m% Jsys(3) = u3;
复制代码

Uifhjvv

帮你顶一下

小小鲁班

0 _) H( l; H4 q9 [你这个没有用到状态变量吧

zaiyiaaaa

你上传一下simulink仿真图？

zzz.dan