Lagrange Interpolation with time interval of
day(s) First date (optional)
Produce file of the selected parameters
Draw data
Spectral analysis (FFT, complex for 2D signal)
Amplitude
PSD
min. frequency or period
max. frequency or period
x linear scale x log scale
y linear scale y log scale
frequency spectrum (in cycle / unit of time)
period spectrum (in unit of time)
First derivative
Produce data
Plot
Weighted least square fit of periodic components (periods in unit of time)
Positive periods
Polynomial of degree
Negative periods
Weighted least square
(take negative periods only for 2 dimensional signal)
Draw residuals
Draw fit and input data
Print residuals
In-phase and out-of-phase terms (a, b) are estimated,
as well as amplitude A and phase φ :
1-D :$$\small X = A \cos[2\pi/T (t-t_0) + \phi] = a \cos[2\pi/T (t-t_0)] + b \cos[2\pi/T (t-t_0)]$$
2-D : $$\small X +i Y = A e^{i [2\pi/T (t-t_0) + \phi]$$
with the reference epoch $$\small t_0$$ = 1/1/2000 0hUT that is :
$$\small X = a \cos[2\pi/T (t-t_0)] - b \sin[2\pi /T (t-t_0)$$ $$\small Y = b cos[2\pi/T (t-t_0)] + a sin[2\pi/T (t-t_0) ]$$
with $$\small a = A \cos\phi$$ $$\small b = A\sin\phi$$
Vondrak low/high pass filter
Panteleev band pass filter (for 2D signal)
Produce data file
Draw filtered data and envelope / phase referred to 2πfct
cycle/time unit Band width f0=
cycle/time unit
This band pass filter was designed by Russian astronomer and gravimetrist V. L. Panteleev. Its frequency transfer function is given by
$$\small T(f) = \frac{f_0^4}{ ( f - f_c)^4 + f_0^4 }$$. At the edges of the window $$\small |f - f_c| = f_0$$ and $$\small T = 0.5$$.
Singular Spectral Analysis (SSA)
- Zoom between
and
The extracted components are decorellated over time windows of
(in the time unit) with the interpolation lag
(in the time unit). Firt step consists in the determination of the
eigenvalues and eigenvectors printed by decreasing weight. Then
5 singular components are reconstructed according to the following
combinations of eigenvectors, to be stated from the analysis of the
eigenvalues .
RC1
Reconstructed Component (RC) based upon eigenvectors N1 and N2
RC2
RC3
RC4
RC5
produce time series (date, signal, RC1,RC2,RC3,RC4,RC5,residuals)
draw
Graphic dimension x
output graphics
png
pdf
ps
Partial Interface with the C-Fortran Libraries SLAVA (C. Bizouard) & MIMOSA (S. Lambert).
Thank you for bringing to our knowledge any possible mistake, mail to : christian.bizouard at obspm.fr
|