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)
, [) Q8 j8 @* vswitch flag,
% q. P5 Y" Z( D- W; b    case 0,
        [sys,x0,str,ts] = mdlIniticializeSizes;
    case 3,
- W1 X" G; k+ H# M( a/ O        sys = mdlOutputs(t,x,u);
    case {1,2,4,9},
9 E: T3 j/ L+ |" F* l4 {        sys = [];
6 x5 F9 u- L' g( T    otherwise,
% y! X  b4 b# d' Q        DAStudio.error('Simulink:blocks:unhandledFlag',num2str(flag));
: {5 O' E4 o4 c& T4 J* J, c% x( {& Rend

function [sys,x0,str,ts] = mdlIniticializeSizes
1 [2 K* H/ {. c; Z& O4 csizes = simsizes;
9 b- K5 Y* d# Z+ E% `% Q, msizes.NumContStates  = 3;    %连续状态个数
- z* K, A: i6 Msizes.NumDiscStates  = 0;    %离散状态个数
sizes.NumOutputs     = 3;    %输出变量个数
sizes.NumInputs      = 9;    %输入信号个数
sizes.DiRFeedthrough = 1;    %输入直接传入输出信号否
/ v7 X9 `/ a" _! ?8 k5 c! @sizes.NumSampleTimes = 0;    % at least one sample time is neededWO  一般来说为1个
. E$ R- X. T9 \6 D/ R1 m* tsys = simsizes(sizes);
x0 = [0.001 0.001 0.001];
str = [];
ts = [];

function sys=mdlOutputs(t,x,u)
e1=u(1);e2=u(2);e3=u(3);
  e  Q6 C! {7 {s21=u(4);s22=u(5);s1=u(6);
es_w=u(9);
) F- a  s) I  Q# kalpha = 0.8;
alpha1 = 5/7; beta1 = 7/5;
5 v2 d' C. S8 C5 ]1 h+ g0 ialpha2 = 7/9; beta2 = 9/7;
k1 = 1;
2 H+ i2 F. @- ck2 = 1.2;
miu = 10;
delta = 5.46;
%assert(all(imag(u)==0), 'u is imaginary or nan');不起作用
) b, f$ P& M$ \9 Dif 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;
    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;
  O1 z7 E- y4 e9 {    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
& C6 @0 s& Q" `9 M( x! n* c" L    u1 = 0;
    u2 = 0;
" L- R( p; s  A1 i) _9 `    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;
: ?' |1 b0 ~6 z* f: V$ jend
sys(1) = u1;
sys(2) = u2;
  M$ P) B0 Q% r0 D: f4 J  psys(3) = u3;
* U7 C. N9 X" G0 P, J9 E5 {复制代码

Uifhjvv

帮你顶一下

小小鲁班

" K) X' m7 P4 D你这个没有用到状态变量吧

zaiyiaaaa

你上传一下simulink仿真图？

zzz.dan