import to data.dat
#
# import
#
# file-format string should be composed from the following letters:
# t ... time
# o ... observed
# a ... adjusted
# c ... calculated
# w ... point weight
# u ... point error
# r ... residuals(observed)
# R ... residuals(adjusted)
# 1 ... attribute #1
# 2 ... attribute #2
# 3 ... attribute #3
# 4 ... attribute #4
# 5 ... weight(attribute #1)
# 6 ... weight(attribute #2)
# 7 ... weight(attribute #3)
# 8 ... weight(attribute #4)
# i ... ignore
#
# note:
# if you import point errors/weights ('u', 'w') the weights
# will be automatically activated for use in calculations!
#
loadproject myproject.p04
#
# loadproject
#
# load project from file
#
fourier 0 nyq o y
#
# fourier
#
# ... lower frequency limit
# ... higher frequency limit or 'nyq' (=computes spectrum up to Nyquist freq)
# ... choose one of these letters:
# o ... observed
# a ... adjusted
# r ... residuals(observed)
# R ... residuals(adjusted)
# 1 ... spectral window 1
# 2 ... spectral window 2
# ... choose y or n
#
fit o 1
fit o 2
#
# fit
#
# ... choose one of these letters:
# o ... observed
# a ... adjusted
# ... calculation mode
# 1 ... improve amplitudes & phases
# 2 ... improve frequencies, amplitudes & phases
#
addharmonics 15
#
# addharmonics
#
# Adds harmonics of the frequency that was found
# in the previous step. This command only defines the additional frequencies,
# no least squares fit is performed at this stage.
#
addsubharmonic 2
#
# addsubharmonic
#
# Add a specific subharmonic (1/) of the frequency that was found
# in the previous step. This command only registers the subharmonic,
# no least squares fit is performed at this stage.
#
savefreqs data.per
#
# savefreqs
#
# save the current frequencies, amplitudes & phases in
#
importfreqs data.per
#
# importfreqs
#
# import frequencies, amplitudes & phases from
#
savefourier 0 nyq r n data.fou
#
# savefourier
#
# compute fourier spectrum and save it in
# ... any number
# ... any number or 'nyq' (=computes spectrum up to Nyquist freq)
# ... choose one of these letters:
# o ... observed
# a ... adjusted
# r ... residuals(observed)
# R ... residuals(adjusted)
# 1 ... spectral window 1
# 2 ... spectral window 2
# ... choose y or n
#
#
saveproject myproject.p04
#
# saveproject
#
# save project in
#
createartificial d 1000.0 1500.0 1.0 0.0 artificial.dat a
#
# createartificial d
#
# generate artificial data based on the current frequency solution
# ... desired start time
# ... desired end time
# ... time step
# ... additional leading/trailing time
# ... output file containing the generated data
# ... choose one of these letters:
# a ... append to outputfile if it already exists
# n ... create new file
#
# createartificial f
#
# generate artificial data based on the current frequency solution
# ... input file with one column containing the times used for data generation
# ... output file containing the generated data
# ... choose one of these letters:
# a ... append to outputfile if it already exists
# n ... create new file
#
#
calcsnr 2.0 0 r snr.txt
#
# calcsnr
#
# compute signal to noise ratio for active frequencies
#
# ... size of the frequency range to be used for the noise calculation:
# [frequency-box/2, frequency+box/2]
# ... defines the step rate
# 0 ... high
# 1 ... medium
# 2 ... low
# ... choose one of these letters:
# o ... observed
# a ... adjusted
# r ... residuals(observed)
# R ... residuals(adjusted)
# 1 ... spectral window 1
# 2 ... spectral window 2#
# ... output file containing the results
#
calcerrors a uncertainties.txt
calcerrors ls uncertainties.txt o 2 t0
calcerrors mc uncertainties.txt 1000 t0 r
#
# calcerrors a
#
# calculate analytically computed uncertainties
# ... output file containing the results
#
# calcerrors ls
# calcerrors ls t0
#
# calculate uncertaintes derived from least-squares fit
# ... output file containing the results
# ... choose one of these letters:
# o ... observed
# a ... adjusted
# ... calculation mode
# 1 ... improve amplitudes & phases
# 2 ... improve frequencies, amplitudes & phases
# t0 ... use original zeropoint in time
#
# calcerrors mc
# calcerrors mc t0
# calcerrors mc r
#
# calculate uncertaintes derived by monte carlo simulation
# ... output file containing the results
# ... number of monte carlo experiments
# t0 ... use original zeropoint in time
# r ... use system time to initialize random generator
#
exit
#
# exit Period04
#