s函数报错

勇往直前11

s函数报错“during flag=i call must be a real vector of length n” ，怎么解决？
代码如下：
; b0 G8 G+ ~- o. _# @6 p) {function [sys,x0,str,ts] = PMSM_chaos_controller(t,x,u,flag)
switch flag,
+ r2 }( K: e( u( u! c, _$ Y    case 0,
        [sys,x0,str,ts] = mdlIniticializeSizes;
/ Z( i$ D7 C; b* }: ]' R" D    case 3,
        sys = mdlOutputs(t,x,u);
    case {1,2,4,9},
& e- u( x! R: ?9 e        sys = [];
    otherwise,
        DAStudio.error('Simulink:blocks:unhandledFlag',num2str(flag));
end
8 N2 S; b) O. C! b: I8 Y  v
function [sys,x0,str,ts] = mdlIniticializeSizes
/ G: M* ^! N  u# R7 U' g% H, Z; o$ U3 Asizes = simsizes;
4 k8 W8 ~. @  o0 W( esizes.NumContStates  = 3;    %连续状态个数
sizes.NumDiscStates  = 0;    %离散状态个数
sizes.NumOutputs     = 3;    %输出变量个数
" c6 |4 r8 y: m, k5 A. `! D- Esizes.NumInputs      = 9;    %输入信号个数
2 W4 n! j. d' q5 q9 ^6 Gsizes.DiRFeedthrough = 1;    %输入直接传入输出信号否
sizes.NumSampleTimes = 0;    % at least one sample time is neededWO  一般来说为1个
  l" ?# Q8 b) Y3 Y( D( Vsys = simsizes(sizes);
9 I& u8 U& {9 Z% X# Y+ xx0 = [0.001 0.001 0.001];
str = [];
ts = [];

function sys=mdlOutputs(t,x,u)
e1=u(1);e2=u(2);e3=u(3);
s21=u(4);s22=u(5);s1=u(6);
es_w=u(9);
  r6 d& w5 U& w) j6 B3 \) d% Xalpha = 0.8;
7 s, G! b0 s, B: L+ H5 g& _# dalpha1 = 5/7; beta1 = 7/5;
alpha2 = 7/9; beta2 = 9/7;
k1 = 1;
k2 = 1.2;
! G+ L9 T# W4 g6 I" @5 ~' k; H$ W$ kmiu = 10;
" B+ _' T; ~: O/ k( Jdelta = 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;
    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;
else
4 O% Q) a5 I, n    u1 = 0;
' s3 V% E9 ?$ Y/ x( K% w8 ^  f& {    u2 = 0;
6 ^& B, J) z. R! O8 A    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;
8 l) L3 K4 e% d2 }! Qend
' d% l$ \6 c! S# I3 y, psys(1) = u1;
sys(2) = u2;
- _" {2 K6 E7 s, Qsys(3) = u3;
复制代码

zzz.dan

zaiyiaaaa

你上传一下simulink仿真图？

小小鲁班

你这个没有用到状态变量吧

Uifhjvv

帮你顶一下