scientimate.wavepropfrompsd¶
Hm0, fp, Tp, Tm01, Tm02, m0, m1, m2 = scientimate.wavepropfrompsd(Syy, f, fcL=0, fcH=None, dispout='no')
Description¶
Calculate wave properties from a power spectral density
Inputs¶
- Syy
- Power spectral density (m^2/Hz)
- f
- Frequency (Hz)
- fcL=0
- Low cut-off frequency, between 0*fs to 0.5*fs (Hz)
- fcH=max(f)
- High cut-off frequency, between 0*fs to 0.5*fs (Hz)
- dispout=’no’
- Define to display outputs or not (‘yes’: display, ‘no’: not display)
Outputs¶
- Hm0
- Zero-Moment Wave Height (m)
- fp
- Peak wave frequency (Hz)
- Tp
- Peak wave period (second)
- Tm01
- Wave Period from m01 (second), Mean Wave Period
- Tm02
- Wave Period from m02 (second), Mean Zero Crossing Period
- m0
- Zero-Moment of the power spectral density (m^2)
- m1
- First-Moment of the power spectral density (m^2)
- m2
- Second-Moment of the power spectral density (m^2)
Examples¶
import scientimate as sm
import numpy as np
N=2**11
fs=8
df=fs/N #Frequency difference
f=np.arange(0,fs/2+df,df) #Frequency vector
f[0]=f[1]/2 #Assign temporarily non-zero value to fisrt element of f to prevent division by zero
Syy=0.016*9.81**2/((2*np.pi)**4*(f**5))*np.exp(-1.25*(0.33/f)**4) #Calculating Spectrum
f[0]=0
Syy[0]=0
Hm0,fp,Tp,Tm01,Tm02,m0,m1,m2=sm.wavepropfrompsd(Syy,f,0,np.max(f),'yes')