Observatoire des Sciences de l’Univers Terre Homme Environnement Temps Astronomie

Accueil > Recherche > Logiciels libres

SigmaTheta Software

par François Vernotte - publié le , mis à jour le

SigmaTheta is a collection of numerical programs for time and frequency metrology.

SigmaTheta is free software, you can redistribute it and/or modify it under the terms of the CeCILL License.

The CeCILL License does not permit this software to be redistributed in proprietary programs.


Download

 Unix/Max OS X
 Windows


Contents

 1col2col

Usage : 1col2col SOURCE TARGET [x]

Transforms a 1 column file into a 2 column file.

The input file SOURCE contains a N line / 1 column table.

The output file TARGET contains a N line / 2 column table with time values (dates) in the first column and the data of SOURCE in the second column.

The optional argument x is the sampling step in seconds. If x is omitted, the sampling rate is assumed to be 1 s.

 X2Y

Usage : X2Y SOURCE TARGET

Transforms a time error sequence x(t) into a normalized frequency deviation sequence Yk.

The input file SOURCE contains a N line / 2 column table with time values (dates) in the first column and time error samples in the second column.

The output file TARGET contains a N-1 line / 2 column table with time values (dates) in the first column and normalized frequency samples in the second column.

 DriRem

Usage : DriRem SOURCE TARGET

Removes the linear drift of a sequence of normalized frequency deviation measurements.

The input file SOURCE contains a 2-column table with time values (dates) in the first column and normalized frequency deviation measurements in the second column.

The output file TARGET contains a 2-column table with the same time values in the first column and the normalized frequency deviation measurements from which the linear drift was removed in the second column.

The coefficients a and b of the drift (equation y=a.x+b) are given in the header of the TARGET file.

 SigmaTheta

Usage : SigmaTheta SOURCE TARGET

Computes the (modified) Allan Deviations of a sequence of normalized frequency deviation measurements, the 95 % confidence intervals (2 sigma), the 68 % confidence intervals (1 sigma), the asymptotes and plots a graph as a postscript file.

The input file SOURCE contains a 2-column table with time values (dates) in the first column and normalized frequency deviation measurements in the second column.

A 7-column table is sent to the standard output with :

  • 1st column : tau values
  • 2nd column : (modified) Allan deviation estimate
  • 3rd column : unbiased estimate
  • 4th column : 2.5 % bound
  • 5th column : 16 % bound
  • 6th column : 84 % bound
  • 7th column : 97.5 % bound.

The file TARGET.gnu is generated for invoking gnuplot. The configuration file ".SigmaTheta.conf" is taken into account (e.g. selection of the modified Allan varance).
The file TARGET.pdf is the pdf file of the gnuplot graph (if the PDF option has been chosen in the configuration file).

If the option ’-m’ is selected, the modified Allan variance is invoked.

If the option ’-h’ is selected, the Hadamard variance is invoked.

If the option ’-p’ is selected, the Parabolic variance is invoked.

Otherwise, the classical Allan variance is invoked.

 ADev

Usage : ADev SOURCE

Computes the Allan Deviations of a sequence of normalized frequency deviation measurements.

The input file SOURCE contains a 2-column table with time values (dates) in the first column and normalized frequency deviation measurements in the second column.

A 2-column table containing tau values (integration time) in the first column and Allan deviation measurement in the second column is sent to the standard output.

A redirection should be used for saving the results in a TARGET file : ADeV SOURCE > TARGET

 MDev

Usage : MDev SOURCE

Computes the modified Allan Deviations of a sequence of normalized frequency deviation measurements.

The input file SOURCE contains a 2-column table with time values (dates) in the first column and normalized frequency deviation measurements in the second column.

A 2-column table containing tau values (integration time) in the first column and modified Allan deviation measurement in the second column is sent to the standard output.

A redirection should be used for saving the results in a TARGET file : MDeV SOURCE > TARGET

  HDev

Usage : HDev SOURCE

Computes the modified Hadamard Deviations of a sequence of normalized frequency deviation measurements.

The input file SOURCE contains a 2-column table with time values (dates) in the first column and normalized frequency deviation measurements in the second column.

A 2-column table containing tau values (integration time) in the first column and Hadamard deviation measurement in the second column is sent to the standard output.

A redirection should be used for saving the results in a TARGET file : HDeV SOURCE > TARGET

 PDev

Usage : PDev SOURCE

Computes the Parabolic Deviations of a sequence of normalized frequency deviation measurements.

The input file SOURCE contains a 2-column table with time values (dates) in the first column and normalized frequency deviation measurements in the second column.

A 2-column table containing tau values (integration time) in the first column and Parabolic deviation measurement in the second column is sent to the standard output.

A redirection should be used for saving the results in a TARGET file : PDeV SOURCE > TARGET

 uncertainties

Usage : uncertainties SOURCE TARGET

Computes the 95 % confidence intervals (2 sigma), the 68 % confidence intervals (1 sigma) of a sequence of Allan Deviations, the asymptotes and plots a graph as a postscript file.

The input file SOURCE contains a 2-column table with tau values (integration time) in the first column and Allan deviation measurement in the second column.

The first tau value is assumed to be equal to the sampling step.
The last tau value is assumed to be equal to the half of the whole time sequence duration.

The output file TARGET contains a 4-column table with :

  • 1st column : tau values
  • 2nd column : Allan deviation measurement
  • 3rd column : inferior bound of the 95 % confidence interval
  • 4th column : superior bound of the 95 % condidence interval.

The file TARGET.gnu is generated for invoking gnuplot.
The file TARGET.ps is the postscript file of the gnuplot graph.

 STshell

STshell is an example of bourne shell chaining ADev and uncertainties for obtaining the same result as SigmaTheta.

Same input and output format as SigmaTheta.

 Asymptote

Usage : Asymptote SOURCE [edfFILE]

Computes the tau^-1, tau^-1/2, tau^0, tau^1/2 and tau asymptotes of a sequence of Allan Deviations.

The input file SOURCE contains a 2-column table with tau values (integration time) in the first column and Allan deviation measurement in the second column.

The optional input file edfFILE contains a 2-column table with tau values in the first column and the equivalent degrees of freedom (edf) of the Allan deviation measurement in the second column.

The Allan deviation estimates are weigted by the tau-values in the default case and by the inverse of the "edf" if the optional file edfFILE has been used.

If the configuration file ".SigmaTheta.conf" contains the line "Unbiased estimates : ON", the fit is performed over the unbiased estimates.

The 5 asymptote coefficients are sent to the standard output separated by a tabulation.

A redirection should be used for saving the results in a TARGET file : Asymptote SOURCE > TARGET

 RaylConfInt

Usage : RaylConfInt ’value’

Computes the mean and the 95 % confidence interval of a chi distribution with ’value’ degrees of freedom, normalized by the square root of ’value’.

The input ’value’ is a floating point number.

The inferior bound, the logarithmic mean and the superior bound are sent to the standard output separated by a tabulation.

 Asym2Alpha

Usage : Asym2Alpha ADevFILE FitFILE

Finds the dominating power law noise (alpha) versus tau.

The input file ADevFILE contains a 2-column table with tau values (integration time) in the first column and Allan deviation measurement in the second column.

The input file FitFILE contains the 5 asymptote coefficients (from tau^-1 to tau^+1) in a 1-line 5-column table.

A a 2-column table with tau values in the first column and the dominating power law noise (from alpha=-2 to +2) in the second column is sent to the standard output.

A redirection should be used for saving the results in a TARGET file : Asym2Alpha ADevFILE FitFILE > TARGET

 AVarDOF

Usage : AVarDOF SOURCE

Computes the degrees of freedom of a sequence of tau values.

The input file SOURCE contains a 2-column table with tau values (integration time) in the first column and the exponent of the power law of the dominating noise in the second column.

The first tau value is assumed to be equal to the sampling step.
The last tau value is assumed to be equal to the half of the whole time sequence duration.

A 2-column table containing tau values (integration time) in the first column and the equivalent degrees of freedom in the second column is sent to the standard output.

A redirection should be used for saving the results in a TARGET file : AVarDOF SOURCE > TARGET

 ADUncert

Usage : ADUncert AdevFILE EdfFILE

Computes the 95 % (2 sigma) and the 68% (1 sigma) confidence intervals of a sequence of Allan Deviations.

The input file ADevFILE contains a 2-column table with tau values (integration time) in the first column and Allan deviation measurement in the second column.

The input file EdfFILE contains a 2-column table with tau values in the first column and the equivalent degrees of freedom (edf) of the Allan deviation measurement in the second column.

A 7-column table is sent to the standard output with :

  • 1st column : tau values
  • 2nd column : Allan deviation estimate
  • 3rd column : unbiased estimate
  • 4th column : 2.5 % bound
  • 5th column : 16 % bound
  • 6th column : 84 % bound
  • 7th column : 97.5 % bound.

A redirection should be used for saving the results in a TARGET file : ADUncert AdevFILE EdfFILE > TARGET

 ADGraph

Usage : ADGraph ADevFILE FitFILE

Plots the graph of the Allan Deviation estimates versus tau.

The input file ADevFILE contains a 2-column table with tau values (integration time) in the first column and Allan deviation measurement in the second column.

The input file FitFILE contains the 5 asymptote coefficients (from tau^-1 to tau^+1) in a 1-line 5-column table.

The ADevFILE is given to store :
- the file ADevFILE.gnu for invoking gnuplot,
- the file ADevFILE.pdf which is the pdf file of the gnuplot graph (if the PDF option has been chosen in the configuration file).

 PSDGraph

Usage : PSDGraph yktFILE

Plots the graph of the Power Spectrum Density (PSD) of normalized frequency deviation versus the frequency.

The input file yktFILE contains a 2-column table with dates (in s) in the first column and frequency deviation samples (yk) in the second column.

- The data of the PSD versus the frequency are stored in the output file yktFILE.psd,
- the file yktFILE.psd.gnu is generated for invoking gnuplot,
- the file yktFILE.psd.pdf is the pdf file of the gnuplot graph (if the PDF option has been chosen in the configuration file).

 YkGraph

Usage : YkGraph yktFILE

Plots the graph of the normalized frequency deviation samples versus time.

The input file yktFILE contains a 2-column table with dates (in s) in the first column and frequency deviation samples (yk) in the second column.

The yktFILE is given to store :
- the file yktFILE.gnu for invoking gnuplot,
- the file yktFILE.pdf which is the pdf file of the gnuplot graph (if the PDF option has been chosen in the configuration file).

  XtGraph

Usage : XtGraph xttFILE

Plots the graph of the time error data versus time.

The input file xttFILE contains a 2-column table with dates (in s) in the first column and time error data x(t) (in s) in the second column.

The xttFILE is given to store :
- the file xttFILE.gnu for invoking gnuplot,
- the file xttFILE.pdf which is the pdf file of the gnuplot graph (if the PDF option has been chosen in the configuration file).

 bruiteur

Usage : bruiteur TARGET

Computes a sequence of simulated time error x(t) and/or normalized frequency deviation yk samples.

A 2-column table containing dates in the first column and time errors in the second column is saved in the file TARGET.xtt and/or a 2-column table containing dates in the first column and frequency deviations in the second column is saved in the file TARGET.ykt.

A redirection should be used for loading the input parameters from a SOURCE file : bruiteur TARGET < SOURCE

 GCoDev

Usage : GCoDev SOURCE1 SOURCE2

Computes the Groslambert CoDeviations (see arXiv:1904.05849) of 2 sequences of normalized frequency deviation measurements.

The input files SOURCE1 and SOURCE2 contain a 2-column table with time values (dates) in the first column and normalized frequency deviation measurements in the second column.

A 2-column table containing tau values (integration time) in the first column and Groslambert codeviation measurement in the second column is sent to the standard output.

A redirection should be used for saving the results in a TARGET file : GCoDeV SOURCE1 SOURCE2 > TARGET.

 Aver

Usage : AVer FILE1 FILE2 FILE3

Computes the average of the 2nd columns of FILE1, FILE2 and FILE3 (suitable to compute the measurement noise of a three-cornered hat system thanks to the closure relationship).

The input files FILE1, FILE2 and FILE3 contain a N line / 2 column table, the 1st columns are assumed identical.

A 2-column table containing the 1st column of FILE1 as the first column of the output table and the average of the 2nd columns of FILE1, FILE2, FILE3 as the the second column of the output table is sent to the standard output.

A redirection should be used for saving the results in a TARGET file : AVer FILE1 FILE2 > TARGET.

 DevGraph

Usage : DevGraph OutFILE DevFILE1 [Dev FILE2] ... [DevFILE8]

Plots the graph of 1 to 8 2-columns Dev files (suitable to plot the 3 ADEV computed by GCoDev on a three-cornered hat system and the ADEV of the measurement noise computed by Aver).

If a file name is preceded by a dash (-DevFILEn), the contain of the second column of this file is multiplied by -1.

The OutFILE is given to store :
- the file OutFILE.gnu for invoking gnuplot,
- the file OutFILE.pdf which is the pdf file of the gnuplot graph (if the PDF option has been chosen in the configuration file).

 STshell2

STshell2 is an example of bourne shell chaining ADev, Asymptote, Asym2Alpha, AVarDOF, ADUncert and ADGraph for obtaining the same result as SigmaTheta.

Same input and output format as SigmaTheta.

 STshell3

STshell3 is an example of bourne shell chaining MDev, Asymptote, Asym2Alpha, AVarDOF, ADUncert and ADGraph for obtaining the same result as SigmaTheta when the modified Allan variance is selected from the configuration file.

Same input and output format as SigmaTheta.


Unix/Mac OS X : The latest Unix/Mac Os X version is available on
https://sourcesup.renater.fr/sigmatheta/SigmaTheta_latest.tar.gz
or directly here

SigmaTheta-4.1 (UNIX/Mac OS X version, tar.gz file)

Windows :

  • exe and dll binaries
    SigmaTheta 2.1.3 (exe and dll windows, zip file)
    This zip file contains .exe and .dll binaries directly usable under Microsoft Windows.
  • source files

Contributors : SigmaTheta is maintained as a git repository on sourcesup. Just type :

git clone https://git.renater.fr/sigmatheta.git

Contriburors are very welcome !


F. Vernotte / A. Kinali / F. Meyer / C. Plantard - August 2019

francois.vernotte_at_obs-besancon.fr