ns3 802.11b PHY model

I use the ubuntu and do not install the chinse input.

The Code: c file requires gnu gsl library, it can be installed easily because many tutorial.

Although the code style is poor, it can be clear for you to read and can be copy to your edit tool.

gcc -Wall -I/usr/local/include -c 80211b.c

gcc -Wall -L/usr/local/lib 80211b.o lgsl 0lgslcblas -lm -o 80211b
./80211b >80211b.txt

 // Copyright 2009, The Boeing Company

 #include "math.h"
 #include "stdlib.h"
 #include "stdio.h"

 #include <gsl/gsl_math.h>
 #include <gsl/gsl_integration.h>
 #include <gsl/gsl_cdf.h>
 #include <gsl/gsl_sf_bessel.h>

 #define min(a,b) ((a)<(b) ? (a) : (b))
 #define max(a,b) ((a)>(b) ? (a) : (b))
 //page 67 802.15 that pdf kanghl
 #define WLAN_SIR_perfect 10.0 // if SIR > 10dB, perfect reception
 #define WLAN_SIR_impossible 0.1 // if SIR < -10dB, impossible to receive

 typedef struct fn_parameter_t
 {
   double beta;
   double n;
 }fn_parameters;

 double QFunction (double x)
 {
   return 0.5 * erfc(x/sqrt(2.0));
 }
 //not understand
 double f(double x, void *params)
 {
   double beta = ((fn_parameters *) params)->beta;
   double n = ((fn_parameters *) params)->n;
   double f = pow( *gsl_cdf_ugaussian_P (x+ beta) - , n-)
     * exp (-x*x/2.0) / sqrt (2.0 * M_PI);

   return f;
 }

 double p_e2(double e2)
 {
     double sep;
     double error;

     fn_parameters params;
     params.beta = sqrt (2.0*e2);
     params.n = 8.0;

     gsl_integration_workspace * w = gsl_integration_workspace_alloc ();

     gsl_function F;
     F.function = &f;
     F.params = &params;

     gsl_integration_qagiu(&F,-params.beta,
     , 1e-, , w, &sep, &error);
     gsl_integration_workspace_free (w);
     if (error == 0.0) sep = 1.0;

 return 1.0 - sep;
 }

 double p_e1(double e1)
 {
     return 1.0 - pow( 1.0 - p_e2 (e1/2.0), 2.0);
 }

 double DbToNoneDb (double x)
 {
     return pow(10.0, x/10.0);
 }

 double NoneDbToDb (double x)
 {
     return 10.0 * log10 (x) ;
 }

 double DQPSKFunction (double x)
 {
     double pi = acos (-1.0);
     return ( (sqrt(2.0) + 1.0) / sqrt(8.0*pi*sqrt(2.0)))
     *(1.0/sqrt(x))
     *exp( - (2.0 - sqrt(2.0)) * x) ;
 }
 //P1MBPS-SYMBOL BER=SER
 double Get80211bDsssDbpskBerIeee(double EcNc)
 {
     double ber;
     if(EcNc > WLAN_SIR_perfect) ber = ;
     else if(EcNc < WLAN_SIR_impossible) ber = 0.5;
     else
       ber = min(QFunction(sqrt(11.0*EcNc)),0.5);
     return ber;
 }

 double Get80211bDsssDbpskBer(double sinr)
 {
   double EbN0 = sinr * 22000000.0 / 1000000.0;
   double ber = 0.5 * exp(-EbN0);
   return ber;
 }

 double Get80211bDsssDqpskBerIeee(double EcNc)
 {
     double ber;
     if (EcNc > WLAN_SIR_perfect) ber = ;
     else if(EcNc < WLAN_SIR_impossible) ber = 0.5;
     else
       ber = min(QFunction(sqrt(5.5*EcNc)),0.5);
     return ber;
 }

 double Get80211bDsssDqpskBer(double sinr)
 {
   // 2 bits per symbol, 1 MSPS
   double EbN0 = sinr * 22000000.0 / 1000000.0 / 2.0;
   double ber = DQPSKFunction(EbN0);
   return ber;
 }

 double Get80211bDsssDqpskCCK5_5BerIeee(double EcNc)
 {
   double ber;
   if(EcNc > WLAN_SIR_perfect) ber = 0.0 ;
   else if(EcNc < WLAN_SIR_impossible) ber = 0.5;
   else
   {
       double pew = 14.0*QFunction(sqrt(EcNc*8.0)) + QFunction(sqrt(EcNc*16.0));
       pew = min(pew, 0.99999);
       ber = 8.0/15.0 * pew;
   }
 return ber;
 }

 double Get80211bDsssDqpskCCK11BerIeee(double EcNc)
 {
   double ber;
   if(EcNc > WLAN_SIR_perfect) ber = 0.0 ;
   else if(EcNc < WLAN_SIR_impossible) ber = 0.5;
   else
   {
     double pew = 24.0*QFunction(sqrt(EcNc*4.0)) +
       16.0*QFunction(sqrt(EcNc*6.0)) +
       174.0*QFunction(sqrt(EcNc*8.0)) +
       16.0*QFunction(sqrt(EcNc*10.0)) +
       24.0*QFunction(sqrt(EcNc*12.0)) +
       QFunction(sqrt(EcNc*16.0));
       pew = min(pew, 0.99999);
       ber = 128.0/255.0 * pew;
   }
   return ber;
 }

 int main (int argc, char *argv[])
 {
   double rss, sinr;
   double totalPkt = 200.0;
   //double noise = 1.552058;
   double noise = ;
   double EcNc, EbN01, EbN02, EbN05, EbN011;
   double ieee1,ieee2,ieee5,ieee11;
   double numBits = (. + . + .) * .;
   double dbpsk,dqpsk,cck16,cck256,sepcck16,sepcck256;

   noise = DbToNoneDb(noise) * 1.3803e-23 * 290.0 * ;
   for (rss=-102.0; rss <= -80.0; rss += 0.1)
   {
     sinr = DbToNoneDb(rss)/1000.0/noise;
     EcNc = sinr * 22000000.0 / 11000000.0; // IEEE sir
     EbN01 = sinr * 22000000.0 / 1000000.0;
     // 2 bits per symbol, 1 MSPS
     EbN02 = sinr * 22000000.0 / 1000000.0 / 2.0;
     EbN05 = sinr * 22000000.0 / 1375000.0 / 4.0;
     EbN011 = sinr * 22000000.0 / 1375000.0 / 8.0;
     // 1=rss, 2=EcNc, 3=EbN01, 4=EbN02, 5=EBN05, 6=EbN011
     printf("%g %g %g %g %g %g ", rss, NoneDbToDb(EcNc),
     NoneDbToDb(EbN01),NoneDbToDb(EbN02),
     NoneDbToDb(EbN05),NoneDbToDb(EbN011));

     ieee1 = Get80211bDsssDbpskBerIeee (EcNc);
     ieee2 = Get80211bDsssDqpskBerIeee (EcNc);
     ieee5 = Get80211bDsssDqpskCCK5_5BerIeee (EcNc);
     ieee11 = Get80211bDsssDqpskCCK11BerIeee (EcNc);
     // 7=ber_ieee1, 8=ber_ieee2, 9=ber_ieee5, 10=ber_ieee11
     printf(" %g %g %g %g ", ieee1, ieee2,ieee5,ieee11);

     ieee1 = totalPkt*pow(-ieee1, numBits);
     ieee2 = totalPkt*pow(-ieee2, numBits);
     ieee5 = totalPkt*pow(-ieee5, numBits);
     ieee11 = totalPkt*pow(-ieee11, numBits);
     // 11=pkt_ieee1, 12=pkt_ieee2, 13=pkt_ieee5, 14=pkt_ieee11
     printf(" %g %g %g %g ", ieee1, ieee2,ieee5,ieee11);

     dbpsk = Get80211bDsssDbpskBer (sinr);
     dqpsk = Get80211bDsssDqpskBer (sinr);
     cck16 = max(, 8.0/15.0*p_e2(4.0*EbN05/2.0));
     cck256 = max(, 128.0/255.0*p_e1(8.0*EbN011/2.0));
     // 15=ber_dbpsk, 16=ber_dqpsk, 17=ber_cck16, 18=ber_cck256
     printf(" %g %g %g %g ", dbpsk, dqpsk,cck16,cck256);

     dbpsk = totalPkt*pow(-dbpsk,numBits);
     dqpsk = totalPkt*pow(-dqpsk,numBits);
     sepcck16 = p_e2(4.0*EbN05/2.0);
     sepcck256 = p_e1(8.0*EbN011/2.0);
     cck16 = totalPkt*pow(1.0-sepcck16,numBits/4.0);
     cck256 = totalPkt*pow(1.0-sepcck256,numBits/8.0);
     // 19=pkt_dbpsk, 20=pkt_dqpsk, 21=pkt_cck16, 22=pkt_cck256
     printf(" %g %g %g %g ", dbpsk, dqpsk,cck16,cck256);
     // 23=sinr
     printf(" %g \n",NoneDbToDb(sinr));
   }
   return ;
 }

80211b.c

gnuplot code

 set term postscript eps color enh "Times-BoldItalic"
 set output '80211b.ieee.pkt.eps'
 set xlabel "RSS (dBm)"
 set ylabel "Packet Received"
 set yrange [:]
 set xrange [-:-]
 plot "80211b.txt" using : title '1M IEEE', \
 "80211b.txt" using : title '2M IEEE', \
 "80211b.txt" using : title '5.5M IEEE', \
 "80211b.txt" using : title '11M IEEE'
 set term postscript eps color enh "Times-BoldItalic"
 set output '80211b.ns3.pkt.eps'
 set xlabel "RSS (dBm)"
 set ylabel "Packet Received"
 set yrange [:]
 set xrange [-:-]
 plot "80211b.txt" using : title '1M DBPSK', \
 "80211b.txt" using : title '2M DQPSK', \
 "80211b.txt" using : title '5.5M CCK16', \
 "80211b.txt" using : title '11M CCK256'
 set term postscript eps color enh "Times-BoldItalic"
 set output '80211b.ieee.sir.eps'
 set xlabel "SIR"
 set ylabel "BER"
 set yrange [10e-:]
 set xrange [-:]
 set logscale y
 plot "80211b.txt" using : title '1M IEEE', \
 "80211b.txt" using : title '2M IEEE', \
 "80211b.txt" using : title '5.5M IEEE', \
 "80211b.txt" using : title '11M IEEE'

plot80211b

reference: http://www.nsnam.org/~pei/80211b.pdf