sinrls.cpp

Go to the documentation of this file.

#include <fstream>
#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
#include "RLSCanceller.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 int nsamples = 2000;
        const double period = 35.2;
        const double PI = acos(-1.0);
        const double noiselevel = 0.5;
        const int shift = 3;
        RLSCanceller Filter(filterlength); 
        UniformRNG Random;
        Filter.SetLambda(lambda);
        Filter.SetDelta(delta);
        
        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);
        Canceller.FilterData();
        
        
        
        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;
        
        
        
        /*
        weightfile.open((noisefilename+".weights").c_str());
        for (i = 0; i < Canceller.WeightHistory.size(); ++i)
        {
                for (j = 0; j < Configuration.filterlength*ntimeseries; ++j)
                {
                        weightfile << j << " " << i << " " << Canceller.WeightHistory.at(i).at(j) << " " << endl;
                        //weightfile << Canceller.WeightHistory.at(i).at(j) << " ";
                }
                //weightfile << endl;
        }
        weightfile.close();*/
}

---