JONSWAP海谱模型

% 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