33 for(
int i=0;
i<
N;
i++) x[
i] =
rnd.Gaus();
34 for(
int i=0;
i<
N;
i++) w[
i] =
rnd.Gaus();
35 for(
int i=0;
i<
N;
i++) y[
i] =
rnd.Gaus();
49 r.
add(tfmap,const_cast<char*>(
"hchannel"));
52 r.
add(A,const_cast<char*>(
"witness"));
61 cout <<
"Ratio rms: (" << clean.
rms() <<
"/" << H.
rms() <<
")= " << clean.
rms()/H.
rms() << endl;
63 cout <<
"x : " << x.
mean() <<
" " << x.
rms() << endl;
64 cout <<
"clean : " << clean.
mean() <<
" " << clean.
rms() << endl;
66 cout <<
"clean-x : " << clean.
mean() <<
" " << clean.
rms() << endl;
75 plot.
plot(eigen,const_cast<char*>(
"alp"),1);
76 plot.
graph[0]->SetTitle(
"Eigen-values of all layers");
84 g->SetTitle(
"Ranking for all layers");
86 c1->SaveAs(
"Test2.png");
static double g(double e)
void setMatrix(double edge=0., double f=1.)
virtual void rate(double r)
std::vector< TGraph * > graph
size_t add(WSeries< double > &target, char *name, double fL=0., double fH=0.)
virtual void start(double s)
void plot(wavearray< double > &, char *=NULL, int=1, double=0., double=0., bool=false, float=0., float=0., bool=false, float=0., bool=false)
virtual size_t size() const
void apply(double threshold=0., char c='a')
wavearray< double > getRank(int n)
void solve(double th, int nE=0, char c='s')
wavearray< double > vfreq
wavearray< double > getVEIGEN(int n=-1)
wavearray< double > getClean()
WDM< double > wdtf(lev, 2 *lev, 6, 10)
virtual double mean() const
virtual void stop(double s)
void Forward(int n=-1)
param: wavelet - n is number of steps (-1 means full decomposition)
virtual void resize(unsigned int)