1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
/***********************************************/
/**
* @file parametrizationTemporalFourier.cpp
*
* @brief Fourier expansion.
* @see ParametrizationTemporal
*
* @author Torsten Mayer-Guerr
* @date 2011-12-28
*
*/
/***********************************************/
#include "base/import.h"
#include "config/config.h"
#include "classes/timeSeries/timeSeries.h"
#include "classes/parametrizationTemporal/parametrizationTemporal.h"
#include "parametrizationTemporalFourier.h"
/***********************************************/
ParametrizationTemporalFourier::ParametrizationTemporalFourier(Config &config)
{
try
{
TimeSeriesPtr timeSeries;
readConfig(config, "fourierDegree", order, Config::MUSTSET, "", "");
readConfig(config, "interval", timeSeries, Config::DEFAULT, "", "intervals of fourier series");
readConfig(config, "includeLastTime", includeLastTime, Config::DEFAULT, "0", "");
if(isCreateSchema(config)) return;
times = timeSeries->times();
isInterval = (times.size() != 0);
if(!isInterval)
times = {Time(), date2time(2500,1,1)};
idxStart = 0;
idxEnd = times.size()-1;
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
Bool ParametrizationTemporalFourier::setInterval(const Time &timeStart, const Time &timeEnd, Bool estimatePerArc)
{
try
{
const UInt idxStartOld = idxStart;
const UInt idxEndOld = idxEnd;
if(estimatePerArc && !isInterval)
times = {timeStart, timeEnd};
idxStart = 0;
while((idxStart+1<times.size()) && (timeStart>=times.at(idxStart+1)))
idxStart++;
idxEnd = idxStart;
while((idxEnd<times.size()-1) && (timeEnd>times.at(idxEnd)))
idxEnd++;
return (idxStartOld != idxStart) || (idxEndOld != idxEnd);
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
void ParametrizationTemporalFourier::parameterName(std::vector<ParameterName> &name) const
{
try
{
for(UInt i=idxStart; i<idxEnd; i++)
for(UInt m=1; m<=order; m++)
{
name.push_back(ParameterName("", "", "fourier.cos("+m%"%i"s+"*x)", times.at(i), times.at(i+1)));
name.push_back(ParameterName("", "", "fourier.sin("+m%"%i"s+"*x)", times.at(i), times.at(i+1)));
}
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
void ParametrizationTemporalFourier::factors(const Time &time, UInt startIndex, std::vector<UInt> &index, std::vector<Double> &factor) const
{
try
{
if((time < times.at(idxStart)) || (time > times.at(idxEnd)) || (!includeLastTime && (time == times.at(idxEnd))))
return;
// find index (interval)
UInt idx = idxStart;
while((idx+1 < idxEnd) && (time >= times.at(idx+1)))
idx++;
const Double angle = 2*PI*(time-times.at(idx)).mjd()/(times.at(idx+1)-times.at(idx)).mjd();
for(UInt m=1; m<=order; m++)
{
index.push_back(2*order*(idx-idxStart)+2*m-2+startIndex);
index.push_back(2*order*(idx-idxStart)+2*m-1+startIndex);
factor.push_back(cos(m*angle));
factor.push_back(sin(m*angle));
}
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
|
