65 double w,
double t,
size_t skip)
73 cout<<
"wavecor::init() invalid arguments"<<endl;
94 short**
P = (
short **)malloc(m*
sizeof(
short*));
95 for(i=0; i<
m; i++) P[i] = (
short*)malloc(m*
sizeof(
short));
106 for(l=-
int(n); l<=
int(n); l++){
116 sign = x>pa[
j] ? 1 : -1;
117 sign *= y>pb[
j] ? 1 : -1;
138 if(i+m+abs(l)>=N)
continue;
145 sign = x>pa[i+j+1] ? 1 : -1;
146 sign *= y>pb[i+j+1] ? 1 : -1;
156 xcor *= 1./double(m*(m-1));
157 for(i=0; i<
m; i++)
delete P[i];
167 double w,
double t,
size_t skip)
175 cout<<
"wavecor::init() invalid arguments"<<endl;
214 for(j=N-n+i; j<
N; j++) x.
data[j]=0.;
222 var = float(x.
data[j*skip]);
231 for(j=N-n+i; j<
N; j++) x.
data[j]=0.;
238 var = float(x.
data[j*skip]);
245 xcor *= sqrt(
double(m))/2.;
256 if(T <= 0.)
return 1.;
275 if(w<=0. || pw==NULL)
return 0.;
283 size_t nL = N-
int(nM+1);
291 for(i=0; i<=
n; i++) {
293 if(x.
data[i] > 0.) count++;
304 if(x.
data[last] > 0.) count--;
306 if(x.
data[last++] > 0.) count++;
wavearray< double > t(hp.size())
virtual size_t size() const
virtual double select(double)
virtual void rate(double r)
wavearray< double > a(hp.size())
std::list< vector_int > cList
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
virtual void init(wavearray< double > &, wavearray< double > &, double, double, size_t=0)
param: two wavearrays, integration window, interval for time lag analysis and skip parameter for runn...
virtual double median(size_t=0, size_t=0) const
virtual void start(double s)
virtual double mean() const
virtual void kendall(wavearray< double > &, wavearray< double > &, double, double, size_t=0)
param: two wavearrays, integration window, interval for time lag analysis and skip parameter for runn...
virtual double coincidence(double, wavecor *)
param: coincidence window, pointer to wavecor object
double fabs(const Complex &x)
virtual void resize(unsigned int)
void cpf(const wavearray< DataType_t > &, int=0, int=0, int=0)
wavecor & operator=(const wavecor &)