s函数报错

勇往直前11

s函数报错“during flag=i call must be a real vector of length n” ，怎么解决？
8 @( c' ^% j' U& Z, D2 C. T代码如下：
function [sys,x0,str,ts] = PMSM_chaos_controller(t,x,u,flag)
$ [8 C! ~" {7 d% v# tswitch flag,
4 M8 r8 D# a4 l7 h. A: y6 B% c    case 0,
0 S0 Z! a8 ~$ Z, Q8 m  J) ~        [sys,x0,str,ts] = mdlIniticializeSizes;
    case 3,
, f7 D3 L7 w9 R( H- e' p" Q        sys = mdlOutputs(t,x,u);
1 B5 q4 \" G9 W; v    case {1,2,4,9},
1 Z+ L! b$ d; ?( J2 i8 y        sys = [];
    otherwise,
        DAStudio.error('Simulink:blocks:unhandledFlag',num2str(flag));
- h# l0 j0 C: K6 |* {" u- g- bend

4 O* D  |( O' s! ~/ u6 `function [sys,x0,str,ts] = mdlIniticializeSizes
& Y& t: A4 n1 t; j2 K! G+ N. Msizes = simsizes;
* D# D2 F0 q/ Osizes.NumContStates  = 3;    %连续状态个数
sizes.NumDiscStates  = 0;    %离散状态个数
sizes.NumOutputs     = 3;    %输出变量个数
sizes.NumInputs      = 9;    %输入信号个数
" Z4 x" [! Z  h9 V! _sizes.DiRFeedthrough = 1;    %输入直接传入输出信号否
" H8 V( P% B7 [0 }sizes.NumSampleTimes = 0;    % at least one sample time is neededWO  一般来说为1个
# w6 K% p- y# h$ qsys = simsizes(sizes);
4 P6 G" P/ q+ d, S+ [' Tx0 = [0.001 0.001 0.001];
  \1 J/ c2 R: e) W& n& I" ]str = [];
( q6 d" ?% m9 _ts = [];
5 }) Q* ?, }3 G7 j4 R. Z/ K
7 |$ p+ ~' o" y5 }& R) q. |function sys=mdlOutputs(t,x,u)
+ P2 J8 V7 L, ~5 p$ T: Ne1=u(1);e2=u(2);e3=u(3);
! X( V) p9 o5 b9 w$ [. W5 Ps21=u(4);s22=u(5);s1=u(6);
es_w=u(9);
alpha = 0.8;
alpha1 = 5/7; beta1 = 7/5;
8 k3 Z" V. p. }- Y7 ealpha2 = 7/9; beta2 = 9/7;
1 {+ `6 x# v: g6 P7 Pk1 = 1;
% Q0 A0 s* c; s/ ?k2 = 1.2;
miu = 10;
' r3 {+ m( ?. k& z7 q* I5 ~) }2 Idelta = 5.46;
7 e4 a; l6 |8 {% L) l6 P%assert(all(imag(u)==0), 'u is imaginary or nan');不起作用
if e3==0
& Q, D8 u% g. |# L  M    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;
; ~# m) {( Y! y, Q1 q; u    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
( N/ l: Y9 g9 D- o2 {3 l! t7 c: w( z    u1 = 0;
    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;
sys(2) = u2;
5 v# x) z# G- F' v7 ?) bsys(3) = u3;
复制代码

Uifhjvv

帮你顶一下

小小鲁班

9 F2 Y7 [+ ~: b' f  j你这个没有用到状态变量吧

zaiyiaaaa

你上传一下simulink仿真图？

zzz.dan