MT Tools/Time Series Noise Removal/testlmscancel.cpp

Go to the documentation of this file.

#include <fstream>
#include <iostream>
#include <vector>
#include <string>
#include <algorithm>
#include "LMSCanceller.h"
#include "RLSCanceller.h"
#include "ApplyFilter.h"
#include "UniformRNG.h"
#include "TimeSeriesComponent.h"
#include "VecMat.h"
namespace ublas = boost::numeric::ublas;
using namespace std;


int main()
{
        
        
        const int filterlength = 2;
        const int inputsize = 2;
        const int outputsize = 1;
        const int maxiter = 1;
        const double mu = 0.2;
        LMSCanceller Filter(filterlength); 
        Filter.SetMu(mu);
        UniformRNG Random;      
        ApplyFilter Canceller(Filter);
        
        gplib::rvec currentinput(inputsize), currdesired(outputsize), currout(outputsize);
        TimeSeriesComponent Input1, Input2, Desired;
        Input1.GetData().assign(3*inputsize,0.0);
        Input2.GetData().assign(3*inputsize,0.0);
        Desired.GetData().assign(3*inputsize,0.0);
        Canceller.AddInputChannel(Input1);
        Canceller.AddInputChannel(Input2);
        Canceller.AddReferenceChannel(Desired);
        for (int iter = 0; iter < maxiter; ++iter)
        {
                for (int i = 0; i < 3*inputsize; ++i)
                {
                        Input1.GetData().at(i) = Random.GetNumber();
                        Input2.GetData().at(i) = Random.GetNumber();
                        Desired.GetData().at(i) = Random.GetNumber();
                }
                Canceller.FilterData();
                cout << "Input1: ";
                copy(Input1.GetData().begin(),Input1.GetData().end(),
                        ostream_iterator<double>(cout," "));
                cout << endl;
                cout << "Input2: ";
                copy(Input2.GetData().begin(),Input2.GetData().end(),
                        ostream_iterator<double>(cout," "));
                cout << endl;
                cout << "Output: ";
                copy(Canceller.GetOutChannels().front()->GetData().begin(),Canceller.GetOutChannels().front()->GetData().end(),
                        ostream_iterator<double>(cout," "));
                cout << endl;
                cout << "Desired: ";
                copy(Desired.GetData().begin(),Desired.GetData().end(),
                        ostream_iterator<double>(cout," "));
                cout << endl;
                cout << "Eps: ";
                copy(Canceller.GetEpsValues().front().begin(),Canceller.GetEpsValues().front().end(),
                        ostream_iterator<double>(cout," "));
                cout << endl << endl;
        }
}

---