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

重头再来

我想用MATLAB/simulink中的S函数做个阵风风速模型。老是提示错误如下：Error, M-File S-function 'gustvelocity' in 'half/S-Function' is expecting 1 parameter(s) while 3 were provided.
0 V4 w  b+ D/ j* `5 _源代码如下：
6 W3 \+ |4 n  T. ifunction [sys,x0,str,ts] = gustvelocity(u,flag,p1,p2,p3)
% _+ L5 B4 l- _' H( \switch flag,
/ b& @6 [- b; Ncase 0,
9 l" {! f& J: N# r' k    [sys,x0,str,ts]=mdlInitializeSizes(p1,p2,p3);
case 3,
+ M7 |% Q$ F+ x* P% M    sys=mdlOutputs(u,p1,p2,p3);
    case{1,2,4,9}
: Y* G; s1 r; K2 I: n        sys=[];
: M; O; w, Y' w: F  otherwise
. F- |2 Y  ^% l4 r8 i# l/ M    error(['Unhandled flag = ',num2str(flag)]);
end
function [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
/ d! A9 R- g3 Ssizes.NumOutputs     = 4;
sizes.NumInputs      = 3;
, N7 b: R& t* @& Jsizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 1;  
sys = simsizes(sizes);
x0  = [];
- [1 }& z" b& g* a( fstr = [];
# L, T6 D/ e" a% N3 ]* k9 wts  = [0 0];
function sys=mdlOutputs(u,p1,p2,p3)
if u<p1
    sys=0;
7 {' S2 p! h/ m* Eelseif u>=p1&&u<=p1+p2
) k. z6 F1 R  L' @9 u  w    sys=p3/2-p3/2*cos(2*pi*((u/p1)-(p2/p1)));
. A9 j, |2 n/ ^9 |3 r9 celse
    sys=0;
end
我不知道哪出了问题，请高手指点一二。

重头再来

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

lhw606

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

重头再来

回复 3# lhw606
/ {, D' J. d, k
! e5 {# h/ L1 v% y" @4 T
/ X" v1 r. _; K    是啊，我掌握的文献中，都是在程序中设定好的风速，这样就有点讨巧。
    我在想是不是可以灵活些，把程序做成模块，我们在模块外部进行设定风速数值。

senyul

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

LJ.G

回复 1# 重头再来