windSpeed_kts = 25;
;
windSpeed_mps = windSpeed_kts * 1852.0 / 3600;

g = 9.8;
U195 = windSpeed_mps;
T1 = 10;

f =  0.02:.001:0.5;
w = 2 * pi * f;
alpha = 8.1e-3;
beta = 0.74;
w0 = g / U195;
S_pm = (2*pi)*((alpha * g^2) ./ (w.^5)) .* exp((-1*beta)*((w0./w).^4));

Hchar = 2.5;
Tavg = 10;
A = 172.75 * (Hchar^2 / Tavg^4);
B = 691 / Tavg^4;
S_ittc = (2*pi)*(A ./ w.^5) .* exp(-B./w.^4);

plot(f,S_pm,f,S_ittc)
legend('PM Spectrum', 'ITTC Spectrum')

