%?JONSWAP譜
clear
close?all

g=9.81;%重力加速度
gamma=3.3;%峰值增強(qiáng)因子
x=[402010];%風(fēng)區(qū)長(zhǎng)度
u=20;%平均風(fēng)速
omega=0:0.01:4;

for?i=1:3;
????xx=g.*x（i）.*1000/power（u2）;
????a=0.076.*power（xx-0.22）;%無(wú)因次常數(shù)
????omega0=22.*g.*power（xx-0.33）./u;%譜峰頻率
????y=-power（omega-omega02）;
????delta=0.07.*（omega<=omeg