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

重头再来

我想用MATLAB/simulink中的S函数做个阵风风速模型。老是提示错误如下：Error, M-File S-function 'gustvelocity' in 'half/S-Function' is expecting 1 parameter(s) while 3 were provided.
源代码如下：
1 o9 r; H  n, X: qfunction [sys,x0,str,ts] = gustvelocity(u,flag,p1,p2,p3)
0 E7 A' q, @) M6 qswitch flag,
case 0,
, ?* a* B* L  D5 M/ w% y. n0 n    [sys,x0,str,ts]=mdlInitializeSizes(p1,p2,p3);
case 3,
    sys=mdlOutputs(u,p1,p2,p3);
    case{1,2,4,9}
  J! e' _) G4 w: U& Q. I        sys=[];
  otherwise
4 a2 ^" @4 ]. j" O  P" t    error(['Unhandled flag = ',num2str(flag)]);
; T: @; l- A- B/ Q+ nend
function [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
4 z. A( e3 x+ S5 @* \( \& |/ d  i1 ~sizes.NumContStates  = 0;
/ M+ V6 A1 R- m# S* G6 e- hsizes.NumDiscStates  = 0;
# d/ K# X% x, q" @, ?sizes.NumOutputs     = 4;
0 v7 X; J  ]- @sizes.NumInputs      = 3;
2 P4 O! _( D! A) Y- Wsizes.DirFeedthrough = 1;
3 W( i) H6 Z4 S. x0 rsizes.NumSampleTimes = 1;  
sys = simsizes(sizes);
" @- \1 t* e* U5 X) N* E4 Cx0  = [];
& B7 C& f/ X9 h0 j+ Astr = [];
& b; O0 y' X5 t2 m* f; Ots  = [0 0];
7 E) m% M' M% Z* ]6 m, Lfunction sys=mdlOutputs(u,p1,p2,p3)
( o2 v# E* k7 M' }; I. X: W* Xif u<p1
5 ?" @- Q; P; O, G- X( Y    sys=0;
6 c& d1 w- J( R5 _- t4 E7 t/ nelseif u>=p1&&u<=p1+p2
9 m' I3 ^1 {; w  m& i  Y+ j9 b  y6 f, ^    sys=p3/2-p3/2*cos(2*pi*((u/p1)-(p2/p1)));
else
5 t8 `' I1 N  d$ j    sys=0;
( N6 T! E% V) Q' Pend
我不知道哪出了问题，请高手指点一二。

重头再来

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

lhw606

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

重头再来

回复 3# lhw606

% G$ b" v; L$ b( u0 p: A
/ K7 s$ M( w( n0 K6 z    是啊，我掌握的文献中，都是在程序中设定好的风速，这样就有点讨巧。
    我在想是不是可以灵活些，把程序做成模块，我们在模块外部进行设定风速数值。

senyul

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

LJ.G

回复 1# 重头再来