基于matlab的海浪模拟程序

function[key]=wave（）
beta1=[input（‘请输入遭遇角：‘）];
V1=[input（‘请输入航速（节）：‘）];
beta=beta1*pi/180;
V=V1*1852/3600;
g=9.8;
A=0.0081;
B=10;
N=25;

for t=1:1:200
    a1=0;
    for w=0.4:0.05:1.6
        kt=1-0.0833*w+0.0032*w^2-0.0000044*w^6;
        kb=1+0.0049*