随机波浪P-M谱

随机海浪往往具有统计特征，组成频率会呈现出某一频率集中的特征。由此而衍生出的海浪谱多种多样。P-M谱由于其结构简单，便于计算，运用于许多领域当中。
用圆频率表示的P-M谱有如下形式：
$S(\omega )=\frac{0.78}{{\omega }^5}\exp [-\frac{3.11}{{\omega }^4{H}_s^2}]$
由$\frac{\partial S(\omega )}{\partial \omega }=0$可求得谱峰频率为：
${\omega }_m=1.253/\sqrt{H_s}$.

% Pierson-Moscowitz (PM) Spectral
% Designed by: JN-Cui 
% Modified on 12/09/2019
%% DEFINITIONS
% H_s - significant wave height;
% g - gravity acceleration;
%% IMPUT PARAMETERS
% H_s - significant wave height; dm - the calculation interval of frequency 
%% FUNCTION

function [S,Omega,omega_m,T_m]=PM_Spectral(H_s,T_s,dm)
i=1;
omega_m=1.253/sqrt(H_s);
f_m=omega_m/2/pi;
T_m=1/f_m;
Omega=0:dm:omega_m*4;
S=length(0:dm:omega_m*4);
for omega=0:dm:omega_m*4
    S(i)=0.78/(omega^5)*exp(-3.11/(omega^4*H_s^2));
    i=i+1;
end
Omega=Omega(2:end);
S=S(2:end);
end