example_channel_model_bsc.cpp

 
#include <Pds/Ra>
#include <Pds/It>
    
int main(void)
{
    unsigned int L=1000000;
    unsigned int M=5;
    unsigned int N=24;
    unsigned int n;
    
    Pds::Matrix U;
    Pds::Vector Up;
    Pds::Vector p(M);
    
    Pds::Vector Rho=Pds::GeomSpace(0.01,0.5,N);
    Pds::Vector BER(N);
    
    for(n=0;n<N;n++)
    {
        if(n==0)    L=200/1.0e-5;
        else        L=200/BER[n-1];
        
        Pds::Vector Us(L);
        p.Fill(Rho[n]);
        Us.FillRandC(0.5);
        
        U=Pds::ChannelModel::Bsc(Us,p);
        Up=U.MeanInRows();
        U.MakeEmpty();
        
        BER.SetRaw(n,Us.Xor(Up).Mean());
        Up.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_channel_model_bsc.png");
    
    Pds::Vector BERteo=Pds::Ceo::Binary::Ber::SymetricModel(Rho,M);
    Pds::Octave::XLabel="BERteo";
    Pds::Octave::YLabel="BER";
    Pds::Octave::Legend={"BERteo","BER"};
    Pds::Octave::Plot::LogLogXYXY(Rho,BERteo,Rho,BER,"temporal.m","example_channel_model_bsc_theoric.png");
    return 0;
}