example_ceo_model_fusion.cpp

 
#include <Pds/Ra>
#include <Pds/It>
 
int main(void)
{
    unsigned int L=100000;
    
    // Dados do canal
    double ps=0.5;
    Pds::Vector f0={1,2,3,4,5};
    Pds::Vector p;
    
    unsigned int N=24;
    Pds::Vector Rho=Pds::GeomSpace(0.01,0.1,N);
    Pds::Vector BER(N);
    
    for(unsigned int n=0;n<N;n++)
    {
        if(n==0)    L=200/1.0e-5;
        else        L=200/BER[n-1];
        
        p=f0*Rho[n];
        
        // Source Us
        Pds::Vector Us(L);
        Us.FillRandC(ps);
        
        // Input channel U
        Pds::Matrix U=Pds::ChannelModel::Bsc(Us,p);
        
        // Second fusion method: NaiveBayes 
        Pds::Vector Unb=Pds::Ceo::Binary::Fusion::NaiveBayes(ps,p,U);
        U.MakeEmpty();
        
        // BER
        BER.SetRaw(n,Us.Xor(Unb).Mean());
        Unb.MakeEmpty();
        
        std::cout<<"p: "<<Rho[n]<<"\tBERm: "<<BER[n]<<"\n";
    }
    
    Pds::Octave::XLabel="Rho";
    Pds::Octave::YLabel="BER";
    Pds::Octave::Legend={"Rho","BER"};
    Pds::Octave::Plot::LogLogXYXY(Rho,Rho,Rho,BER,"temporal.m","example_ceo_model_fusion.png");
    
    return 0;
}