lemiicatest.cpp

Go to the documentation of this file.

#include <iostream>
#include <itpp/base/fastica.h>
#include <itpp/base/mat.h>
#include <string>
#include <fstream>
#include <vector>
#include <iterator>
using namespace std;

void WriteToFile(const itpp::Vec<double> &in,const  itpp::Vec<double> &out,const string &base)
{
        
        ofstream ints((base+".in").c_str());    
        ofstream outts((base+".out").c_str());
        
        const int size = in.length();
        for (int i = 0; i < size; ++i)
        {
                ints << in(i) <<  endl;
                outts << out(i) <<  endl;
        }
}


int main()
{
        const int nchannels = 2;
        
        
        string name1, name2;
        cout << "Infile 1: ";
        cin >> name1;
        cout << "Infile 2: ";
        cin >> name2;
        ifstream file1(name1.c_str());
        ifstream file2(name2.c_str());
        vector<double> in1((istream_iterator<double>(file1)), istream_iterator<double>());
        vector<double> in2((istream_iterator<double>(file2)), istream_iterator<double>());
        const size_t tslength = min(in1.size(),in2.size());
        cout << tslength <<endl;
    itpp::mat Input(nchannels,tslength);
    for (size_t i = 0; i < tslength; ++i)
    {
        Input(0,i) = in1.at(i);
        Input(1,i) = in2.at(i);
    }
        
        itpp::Fast_ICA fastica(Input);
        fastica.set_nrof_independent_components(nchannels);
        //fastica.set_non_linearity(FICA_NONLIN_TANH);
        //fastica.set_approach(FICA_APPROACH_DEFL);
        fastica.separate();
        itpp::mat Output = fastica.get_independent_components();
        WriteToFile(Input.get_row(0),Output.get_row(0),"1");
        WriteToFile(Input.get_row(1),Output.get_row(1),"2");
        cout << "Mix Mat: " << endl << fastica.get_mixing_matrix() << endl;
        cout << "Whit Mat: " << endl << fastica.get_whitening_matrix() << endl;
        
}

---