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

重头再来

我想用MATLAB/simulink中的S函数做个阵风风速模型。老是提示错误如下：Error, M-File S-function 'gustvelocity' in 'half/S-Function' is expecting 1 parameter(s) while 3 were provided.
源代码如下：
function [sys,x0,str,ts] = gustvelocity(u,flag,p1,p2,p3)
/ X% T% i. H1 Y0 L( V+ Uswitch flag,
8 C/ k; o8 n. Scase 0,
    [sys,x0,str,ts]=mdlInitializeSizes(p1,p2,p3);
case 3,
3 M: E: ^' n" E! P$ r0 {    sys=mdlOutputs(u,p1,p2,p3);
    case{1,2,4,9}
( I! i/ d- U6 A$ k$ a        sys=[];
  otherwise
    error(['Unhandled flag = ',num2str(flag)]);
end
function [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
) h% b; `9 p0 y. k% H+ ?sizes.NumContStates  = 0;
6 A5 Y, _7 L5 L* ssizes.NumDiscStates  = 0;
2 @" c/ [3 K8 l2 c6 c( {8 z7 Ssizes.NumOutputs     = 4;
sizes.NumInputs      = 3;
* N0 ]/ C: P( Z5 u3 Rsizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 1;  
sys = simsizes(sizes);
$ n: z6 i. L0 S1 T0 K) K# S' s+ bx0  = [];
- Z# T! n- T3 q2 r+ Zstr = [];
ts  = [0 0];
( Y7 r! |1 T! B0 U/ a! o7 ]0 _function sys=mdlOutputs(u,p1,p2,p3)
if u<p1
    sys=0;
7 x  q  Y, _. \elseif u>=p1&&u<=p1+p2
    sys=p3/2-p3/2*cos(2*pi*((u/p1)-(p2/p1)));
( c2 M) d. t3 ]% S$ s: d4 kelse
    sys=0;
end
7 e' m4 }( \5 J; N我不知道哪出了问题，请高手指点一二。

重头再来

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

lhw606

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

重头再来

回复 3# lhw606
. x% w! i% X) M6 d) F. b$ X
9 b* o0 ?5 @' Y3 o: I- L: ^% k* \
: E6 |: L" w2 f2 y5 g6 a    是啊，我掌握的文献中，都是在程序中设定好的风速，这样就有点讨巧。
8 P( [5 v" b' c3 k  U- M2 }    我在想是不是可以灵活些，把程序做成模块，我们在模块外部进行设定风速数值。

senyul

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

LJ.G

回复 1# 重头再来