eeva/jamulus
1
0
Fork 0
Code Issues Pull requests Releases Wiki Activity

implement new value bound for time variance measurement class

Browse source
This commit is contained in:
Volker Fischer 2009-08-15 14:16:51 +00:00
parent d6ec1b8978
commit e354d985bc
1 changed files with 639 additions and 619 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

36
src/util.h
View file

@ -562,16 +562,22 @@ public:
class CCycleTimeVariance class CCycleTimeVariance
{ {
public: public:
CCycleTimeVariance() : iBlockLengthAtSystemSampleRate ( 0 ) {} CCycleTimeVariance() : iBlockLengthAtSystemSampleRate ( 0 ),
dIntervalTime ( 0.0 ), iNewValueBoundFactor ( 0 ) {}
virtual ~CCycleTimeVariance() {} virtual ~CCycleTimeVariance() {}
void Init ( const int iNewBlockLengthAtSystemSampleRate, void Init ( const int iNewBlockLengthAtSystemSampleRate,
const int iNewSystemSampleRate, const int iNewSystemSampleRate,
const int iHistoryLengthTime ) const int iHistoryLengthTime,
const int iNewNewValueBoundFactor = 4 )
{ {
// store block size and sample rate // store block size and new value bound factor
iBlockLengthAtSystemSampleRate = iNewBlockLengthAtSystemSampleRate; iBlockLengthAtSystemSampleRate = iNewBlockLengthAtSystemSampleRate;
iSystemSampleRate = iNewSystemSampleRate; iNewValueBoundFactor = iNewNewValueBoundFactor;
// calculate interval time
dIntervalTime = static_cast<double> (
iBlockLengthAtSystemSampleRate ) * 1000 / iNewSystemSampleRate;
// calculate actual moving average length and initialize buffer // calculate actual moving average length and initialize buffer
RespTimeMoAvBuf.Init ( iHistoryLengthTime * RespTimeMoAvBuf.Init ( iHistoryLengthTime *
@ -593,10 +599,23 @@ public:
// we want to calculate the standard deviation (we assume that the mean // we want to calculate the standard deviation (we assume that the mean
// is correct at the block period time) // is correct at the block period time)
const double dCurAddVal = ( (double) ( CurTime - TimeLastBlock ) - const double dCurAddVal =
( iBlockLengthAtSystemSampleRate * 1000 / iSystemSampleRate ) ); static_cast<double> ( CurTime - TimeLastBlock ) - dIntervalTime;
RespTimeMoAvBuf.Add ( dCurAddVal * dCurAddVal ); // add squared value if ( iNewValueBoundFactor > 0 )
{
// check if new value is in range
if ( fabs ( dCurAddVal ) < ( iNewValueBoundFactor * dIntervalTime ) )
{
// add squared value
RespTimeMoAvBuf.Add ( dCurAddVal * dCurAddVal );
}
}
else
{
// new value bound is not used, add new value (add squared value)
RespTimeMoAvBuf.Add ( dCurAddVal * dCurAddVal );
}
// store old time value // store old time value
TimeLastBlock = CurTime; TimeLastBlock = CurTime;
@ -613,7 +632,8 @@ protected:
CMovingAv<double> RespTimeMoAvBuf; CMovingAv<double> RespTimeMoAvBuf;
int TimeLastBlock; int TimeLastBlock;
int iBlockLengthAtSystemSampleRate; int iBlockLengthAtSystemSampleRate;
int iSystemSampleRate; double dIntervalTime;
int iNewValueBoundFactor;
}; };
#endif /* !defined ( UTIL_HOIH934256GEKJH98_3_43445KJIUHF1912__INCLUDED_ ) */ #endif /* !defined ( UTIL_HOIH934256GEKJH98_3_43445KJIUHF1912__INCLUDED_ ) */