testcomplexica.cpp

Go to the documentation of this file.

#include "ComplexICA.h"
#include "UniformRNG.h"
#include "miscfunc.h"
#include <fstream>
int main()
{
        const int nchannels = 2;
        const int nsamples = 2000;
        const float minrand = -0.5;
        const float maxrand = 0.5;
        UniformRNG Random;
        
        //Convolve(Master,Corr,Output);
        
        cmat orig_source(nchannels,nsamples);
        cmat input(nchannels,nsamples);
        for (int j = 0; j < nsamples; ++j)
        {
                        orig_source(0,j) = std::complex<double>(sin(j/125.23),0.0);
                        orig_source(1,j) = std::complex<double>(Random.GetNumber(minrand,maxrand),Random.GetNumber(minrand,maxrand)); 
        }
        cmat orig_mix(nchannels,nchannels);
        for (int i = 0; i < nchannels; ++i)
                for (int j = 0; j < nchannels; ++j)
                        orig_mix(i,j) = std::complex<double>(Random.GetNumber(minrand,maxrand),Random.GetNumber(minrand,maxrand)); 
        for (int i = 0; i < nsamples; ++i)
                ublas::column(input,i) = ublas::prod(orig_mix,ublas::column(orig_source,i));
        cmat source(nchannels,nsamples);
        cmat mix(nchannels,nchannels);
        std::cout << "input: " << input << std::endl << std::endl;
        ComplexICA(input,source,mix);
        std::cout << std::endl << std::endl;
        //std::cout << "orig_source: " << orig_source << std::endl << std::endl;
        std::cout << "orig_mix: " << orig_mix << std::endl<< std::endl;
        //std::cout << "est_source: " << source << std::endl<< std::endl;
        std::cout << "est_mix: " << mix << std::endl<< std::endl; 
        std::ofstream outfile("compica.out");
        for (int i = 0; i < nsamples; ++i)
                outfile << i << " " <<  source(0,i).real() << " " << source(1,i).real()<< " " <<  source(0,i).imag() << " " << source(1,i).imag() << std::endl;
}

---