想做个风速变化模型遇到困难，请高手指点。

重头再来

我想用MATLAB/simulink中的S函数做个阵风风速模型。老是提示错误如下：Error, M-File S-function 'gustvelocity' in 'half/S-Function' is expecting 1 parameter(s) while 3 were provided.
/ y, l4 m; k3 p1 k) l  r源代码如下：
& F- G2 j( i& m1 f  ]1 afunction [sys,x0,str,ts] = gustvelocity(u,flag,p1,p2,p3)
switch flag,
4 I9 [9 g$ W. r, u* ^& F/ j, mcase 0,
    [sys,x0,str,ts]=mdlInitializeSizes(p1,p2,p3);
3 d) Q( s. m) E8 [case 3,
    sys=mdlOutputs(u,p1,p2,p3);
    case{1,2,4,9}
+ p( k  [3 e6 S* s* C$ t        sys=[];
  otherwise
& q: e, U8 w. O! f' k. M- H    error(['Unhandled flag = ',num2str(flag)]);
: q* h3 S5 C+ E9 {: A) \end
3 A1 ~) ~4 P  n/ ~" V7 o5 Yfunction [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates  = 0;
( b0 i6 e/ H9 f5 gsizes.NumDiscStates  = 0;
1 A; [6 B; `5 Z2 nsizes.NumOutputs     = 4;
sizes.NumInputs      = 3;
6 m/ X8 n% d) L# {0 h* @9 F, Qsizes.DirFeedthrough = 1;
3 S, M+ b  x+ Y- z4 osizes.NumSampleTimes = 1;  
sys = simsizes(sizes);
x0  = [];
' O4 k6 W( {. w7 astr = [];
ts  = [0 0];
$ f4 Z% l( f7 W' Q+ ^0 f5 |/ nfunction sys=mdlOutputs(u,p1,p2,p3)
  r0 }0 L8 E% w3 G% Z7 yif u<p1
    sys=0;
elseif u>=p1&&u<=p1+p2
7 ^5 g, S0 |& K, Y8 w    sys=p3/2-p3/2*cos(2*pi*((u/p1)-(p2/p1)));
+ }& y) u+ T$ X7 e5 s$ j& z2 Qelse
    sys=0;
end
2 K. a- {' N: O' `我不知道哪出了问题，请高手指点一二。

重头再来

解释下：p1表示启动时间；p2表示阵风周期；p3最大阵风风速；

lhw606

是用S 函数编写的牙~~~~~~~~~~~~~~

重头再来

回复 3# lhw606
% `4 C8 ?- P4 |/ Y
8 o" Q6 |+ c3 J
    是啊，我掌握的文献中，都是在程序中设定好的风速，这样就有点讨巧。
    我在想是不是可以灵活些，把程序做成模块，我们在模块外部进行设定风速数值。

senyul

遇到同样的问题啊，我是用S函数写，但也没成功，还在思考中。。。。

LJ.G

回复 1# 重头再来