gusucode.com > 基于matlab的阵风风速模型源码程序 > 基于matlab的阵风风速模型源码程序/code/zhen_feng.m
function vws=zhen_feng(in) t1g=in(1); t2g=in(2); vgmax=in(3); t=in(4); if t<t1g vws=0; elseif t>=t1g&&t<=(t1g+t2g) vws=0.5*vgmax*(1-cos(2*pi*(t/t2g-t1g/t2g))); else vws=0; end