sinamrls.cpp

Go to the documentation of this file.

#include <fstream>
#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
#include "AMRLSCanceller.h"
#include "ApplyFilter.h"
#include "TimeSeriesComponent.h"
#include "UniformRNG.h"
using namespace std;

int main()
{
        ofstream weightfile;

        
        const int filterlength = 5;
        const double lambda = 0.9999;
        const double delta = 1.1;
        const double alpha = 0.001;
        const int nsamples = 2000;
        const double period = 35.2;
        const double PI = acos(-1.0);
        const double noiselevel = 0.5;
        const int shift = 3;
        AMRLSCanceller Filter(filterlength,delta,lambda,alpha); 
        UniformRNG Random;
        
        ApplyFilter Canceller(Filter);
        
        TimeSeriesComponent  reference, input;
        
        for (int i = 0; i < nsamples; ++i)
        {
                reference.GetData().push_back(2.0*cos( 2.0 * PI * i/period)+Random.GetNumber(-noiselevel,noiselevel));
                input.GetData().push_back(sin( 2.0 * PI * i/period));
        }

        Canceller.AddReferenceChannel(reference);
        Canceller.AddInputChannel(input);
        Canceller.SetShift(shift);
        cout << "Initial Lambda: " << lambda << endl;
        Canceller.FilterData();
        cout << "Final Lambda: " << Filter.GetLambda() << endl;
        
        
        ofstream cleanfile("test.clean");
        for (int i = 0; i < nsamples - filterlength-shift; ++i)
                cleanfile << i << " " << Canceller.GetOutChannels().front()->GetData().at(i) << " " << reference.GetData().at(i) << " " 
                                         << 2.0*cos( 2.0 * PI * i/period) << endl;
        //copy(Canceller.GetOutChannels().front()->GetData().begin(),Canceller.GetOutChannels().front()->GetData().end(),ostream_iterator<double>(cleanfile,"\n"));
        ofstream difffile("test.diff");
        for (int i = 0; i < nsamples - filterlength-shift; ++i)
                difffile << i << " " << Canceller.GetOutChannels().front()->GetData().at(i) -2.0*cos( 2.0 * PI * i/period)  
                          << " " << reference.GetData().at(i) - 2.0*cos( 2.0 * PI * i/period) << endl;
        ofstream epsfile("test.eps");
        copy(Canceller.GetEpsValues().front().begin(),Canceller.GetEpsValues().front().end(),ostream_iterator<double>(epsfile,"\n"));
        
        cout << "Last Weights: " << endl;
        Filter.PrintWeights(cout);
        cout << endl;
        
}

---