ittc1.asv

ittc标准海浪谱的matlab实现

%S(w)=A/w^5*exp(-B/w^4)  P-M谱

   g=9.81; %重力加速度 kg.m/s^2    
   
   H31; %有义波高  m
   %V; %航速  m/s
   
   A=0.0081*g^2; %单参数谱
   B=3.11/(H31^2);
   w=0.25;
   deltaw=0.08;
   for i=1:50
       omiga(i)=w+(i-1)*deltaw;
       s(i)=((A/omiga(i)^5)*exp(-B/omiga(i)^4))*(omiga(i)^4/g^2);
       kesai(i)=sqrt(2*s(i)*deltaw);
       rand('seed');
       emslong(i)=2*pi*rand(1); 
   end
   %plot(omiga,s)
   for t=1:6000
       wave(t)=0;
       for i=1:50
           wave(t)=wave(t)+kesai(i)*cos(omiga(i)*(t-1)/10+emslong(i));
       end
   end
%plot(wave)
for i=1:6000
    alpha(i)=wave(i)*180/pi;
end
plot([0:0.1:599.9],alpha)
t=0:0.1:599.9;
alphabqj=[t;alpha];
save wave_slope alphabqj;