Script de la fonction JONSWAP sur Matlab
function [S,Hs,Ts] = jonswapF(fr,U,F)
%JONSWAPFU Summary of this function goes here
%
fonction qui traite le spectre de jonswap sur les houles
% une approche par le fetch et la vitesse du vent au large
% OUNIS Faiza Sabrina 5éme année ENSMAL
% sabrinafaizaounis@gmail.com
%input
%fetch en m et U en m/s
g=9.81; %accéleration de la gravité m/s
f_pic=3.5*(g/U)*(g*F/U^2)^(-0.33); %peak frequency
alpha= 0.076*(g*F/U^2)^(-0.22); %facteur alpha
gamma=7*(g*F/U^2)^(-0.143); %facteur gamma
if fr sigma=0.07;
else
sigma=0.09;
end
a=exp(-((fr-f_pic).^2)/((2*(sigma^2)*f_pic^2))); %puissance a en fonction de
la frequence
S=(alpha*g^2/(2*pi)^4)*fr.^(-5).*exp(
-1.25*(f_p./fr).^(4)).*gamma.^a;
%spectre d'énérgie
plot(fr,S);
figure
plot(fr,S)
xlabel('frequence (Hz)');
ylabel('Spectrum (m^2/Hz)');
grid;
m_0=trapz(fr,S); %moment m0
Hm=4*sqrt(m_0);
disp(Hm);
Y_m1=fr.*S;
m_1=trapz(fr,Y_m1); %moment m1
disp(m_1);
T_01=m_0/m_1;
Y_m2=(fr.^2).*S;