%?JONSWAP譜
clear
close?all

g=9.81;%重力加速度
gamma=3.3;%峰值增強因子
x=[402010];%風區長度
u=20;%平均風速
omega=0:0.01:4;

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