5 #define nLAYERS 32 // layers used in TF WDM transform
7 #define DURATION 2 // sec
9 #define TIME_PIXEL_INDEX 100
10 #define FREQ_PIXEL_INDEX 20
25 cout <<
"ts size = " <<
ts.
size() <<
" ts rate = " <<
ts.
rate() << endl;
29 cout <<
"tf size = " << tf.
size() << endl;
32 double dT = 1./(2*dF);
34 cout <<
"rate(hz) : " <<
RATE <<
"\t layers : " <<
nLAYERS
35 <<
"\t dF(hz) : " << dF <<
"\t dT(ms) : " << dT*1000. << endl;
40 int index = (levels+1)*itime+ifreq;
44 double freq = (ifreq>0) ? ifreq*dF : dF/4;
46 cout <<
"PIXEL TIME = " << time <<
" sec " << endl;
47 cout <<
"PIXEL FREQ = " << freq <<
" Hz " << endl;
76 cout <<
"mean time : " << tavr << endl;
81 trms+=x[
i]*x[
i]*pow(t-tavr,2);
85 cout <<
"rms time : " << trms*1000 <<
" (ms) " << endl;
101 double e = x[2*
i]*x[2*
i]+x[2*
i+1]*x[2*
i+1];
106 cout <<
"mean frequency : " << favr << endl;
109 for(
int i=0;
i<x.
size()/2;
i++) {
111 double e = x[2*
i]*x[2*
i]+x[2*
i+1]*x[2*
i+1];
112 frms+=e*pow(f-favr,2);
116 cout <<
"rms frequency : " << frms <<
" (Hz) " << endl;
wavearray< double > t(hp.size())
virtual size_t size() const
int getBaseWave(int m, int n, SymmArray< double > &w)
virtual void rate(double r)
cout<< endl;cout<< "ts size = "<< ts.size()<< " ts rate = "<< ts.rate()<< endl;tf.Forward(ts, wdm);int levels=tf.getLevel();cout<< "tf size = "<< tf.size()<< endl;double dF=tf.resolution();double dT=1./(2 *dF);cout<< "rate(hz) : "<< RATE<< "\t layers : "<< nLAYERS<< "\t dF(hz) : "<< dF<< "\t dT(ms) : "<< dT *1000.<< endl;int itime=TIME_PIXEL_INDEX;int ifreq=FREQ_PIXEL_INDEX;int index=(levels+1)*itime+ifreq;double time=itime *dT;double freq=(ifreq >0)?ifreq *dF:dF/4;cout<< endl;cout<< "PIXEL TIME = "<< time<< " sec "<< endl;cout<< "PIXEL FREQ = "<< freq<< " Hz "<< endl;cout<< endl;wavearray< double > x
WDM< double > wdm(nLAYERS, nLAYERS, 6, 10)
gwavearray< double > * gx
wavearray< double > ts(N)
void Forward(int n=-1)
param: wavelet - n is number of steps (-1 means full decomposition)
virtual void resize(unsigned int)