Commit f567ef73 authored by thomas.forbriger's avatar thomas.forbriger
Browse files

ts/croposp [WP]: provide two sets of test series

parent 09c56e80
......@@ -42,6 +42,11 @@ clean:
TS_DT=0.01
# number of samples
TS_N=1000000
# number of samples to shift series for test of coherent time series
TS_SHIFT=1000
#
#
# total number of samples to be produce
#
# set random number generator type
# --------------------------------
......@@ -52,33 +57,44 @@ RNGSEED=0
# create to raw time series
rawseries.bin:
randomseries -v -o -t bin -rngtype $(RNGTYPE) -nseries 2 \
-nsamples $(TS_N) -dt $(TS_DT) -seed $(RNGSEED) \
-nsamples $(shell dc -e "$(TS_N) $(TS_SHIFT) + p") \
-dt $(TS_DT) -seed $(RNGSEED) \
-std 1. -mean 0. $@
# place the two time series in separate files
series1.bin series2.bin: \
%.bin: rawseries.bin
any2any --verb --over --itype bin --otype bin \
$@ $< t:$(patsubst series%.bin,%,$@)
tijerasxx --verbose --overwrite --iformat bin --oformat bin \
--samplesl $(TS_SHIFT) $@ $< t:$(patsubst series%.bin,%,$@)
# ----------------------------------------------------------------------
# create shifted versions of series1
# these are time series which are perfectly coherent in our understanding of
# the corresponding continues time signal, but which may appear incoherent
# due to windowing of segments in the discrete analysis
# create shifted version
# the shifted and the not shifted version are time series which are perfectly
# coherent in our understanding of the corresponding continues time signal,
# but which may appear incoherent due to windowing of segments in the
# discrete analysis
TS_SHIFT=1000
seriesA.bin: series1.bin
series_shifted.bin: rawseries.bin
tijerasxx --verbose --overwrite --iformat bin --oformat bin \
--samplesf $(TS_SHIFT) xxx$@ $<
--samplesf $(TS_SHIFT) xxx$@ $< t:1
sehefixx xxx$@ $@ -verbose -overwrite -itype bin -otype bin \
-st $(shell sigval -format %HT -type bin $<)
seriesB.bin: series1.bin
tijerasxx --verbose --overwrite --iformat bin --oformat bin \
--samplesl $(TS_SHIFT) $@ $<
-st $(shell sigval -format %HT -type bin $< t:1)
# ----------------------------------------------------------------------
# pairs of time series
# --------------------
# incoherent pair (as far as rng produces incoherent random numbers)
seriesI1.bin: series1.bin; /bin/cp -fv $< $@
seriesI2.bin: series2.bin; /bin/cp -fv $< $@
# coherent series
seriesC1.bin: series1.bin; /bin/cp -fv $< $@
seriesC2.bin: series_shifted.bin; /bin/cp -fv $< $@
# ----------------------------------------------------------------------
# provide a time series plot to confirm incoherence
rawsig.ps: series1.bin series2.bin seriesA.bin seriesB.bin
rawsig.ps: seriesI1.bin seriesI2.bin seriesC1.bin seriesC2.bin
stuplox -d $@/cps -s x -i -ty bin $^
# ======================================================================
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment