It is known that any waveform (signal) can be rewritten as the sum of sinusoidal functions. Fourier Transform is the tool for the realization. By analyzing a signal in the frequency domain, we may know the characteristics of a system.
Fourier Series written as
where T is the fundamental period, and are the coefficients of the Fourier Series. The coefficients are given by
The signal can be written using complex exponential function that
Fourier transform is defined as
Note that Laplace transform has a similar form as Fourier transform that
where is a complex number frequency parameter .
Fourier series coefficients of a time-series signal and its power spectral density can be calculated in MATLAB:
function [f,amp,psdx]=fft_freq(SAMPL,H) % input: SAMPL: signal, H: sampling time % output: f, frequency [1/s], amp: Fourier series coefficients N = length(SAMPL); if(mod(N,2) ~= 0) SAMPL = SAMPL(1:N-1); N = N-1; end f = (0:N/2); f = f/N*1/H; sampl_fft2 = fft(SAMPL,N);% sampl_fft1 = sampl_fft2(1:N/2+1);% nyquist frequency psdx = (1/(1/H*N)) * abs(sampl_fft1).^2; amp = abs(sampl_fft1/N); amp(2:end-1) = 2*amp(2:end-1); psdx(2:end-1)=2*psdx(2:end-1); fonts = 20; figure; h=plot(f,amp,'-'); grid; ylabel('Amplitude','FontSize',fonts); set(gca,'Fontsize',fonts); set(h,'Linewidth',2.5); xlabel('Frequency [Hz]','FontSize',fonts); figure; h=plot(f,10*log10(psdx),'-'); grid; ylabel('PSD [dB/Hz]','FontSize',fonts); set(gca,'Fontsize',fonts); set(h,'Linewidth',2.5); xlabel('Frequency [Hz]','Fontsize',fonts); end
Response Amplitude Operator
In the community of marine science, a RAO is an engineering method used to evaluate the likely behavior of a marine system (such as a ship) when operating at sea. The RAO is only defined when the motion of a marine structure can be assumed to be linear.
where is the linear excitation force complex amplitude per wave height, the added mass, the damping term, and the restoring coefficient.
RAOs are usually obtained as time-series data from model tests in waves, or from CFD. In model tests, as the generated waves would not be perfect, time-series data would include other components. In order to obtain the RAOs, Fourier series coefficients of both wave signal and motion signal are firstly calculated. RAO at the specified frequency can then be calculated.