16#include <gsl/gsl_cdf.h>
17#include <gsl/gsl_integration.h>
18#include <gsl/gsl_math.h>
19#include <gsl/gsl_sf_bessel.h>
36 return ((M_SQRT2 + 1.0) / std::sqrt(8.0 * M_PI * M_SQRT2)) * (1.0 / std::sqrt(x)) *
37 std::exp(-(2.0 - M_SQRT2) * x);
44 double EbN0 = sinr * 22000000.0 / 1000000.0;
45 double ber = 0.5 * std::exp(-EbN0);
46 return std::pow((1.0 - ber),
static_cast<double>(nbits));
53 double EbN0 = sinr * 22000000.0 / 1000000.0 / 2.0;
55 return std::pow((1.0 - ber),
static_cast<double>(nbits));
64 double EbN0 = sinr * 22000000.0 / 1375000.0 / 4.0;
65 double sep = SymbolErrorProb16Cck(4.0 * EbN0 / 2.0);
66 return std::min(1.0, std::pow(1.0 - sep, nbits / 4.0));
68 NS_LOG_WARN(
"Running a 802.11b CCK Matlab model less accurate than GSL model");
82 double a1 = 5.3681634344056195e-001;
83 double a2 = 3.3092430025608586e-003;
84 double a3 = 4.1654372361004000e-001;
85 double a4 = 1.0288981434358866e+000;
86 ber = a1 * std::exp(-std::pow((sinr - a2) / a3, a4));
88 return std::min(1.0, std::pow((1.0 - ber),
static_cast<double>(nbits)));
99 double EbN0 = sinr * 22000000.0 / 1375000.0 / 8.0;
100 double sep = SymbolErrorProb256Cck(8.0 * EbN0 / 2.0);
101 return std::min(1.0, std::pow(1.0 - sep, nbits / 8.0));
103 NS_LOG_WARN(
"Running a 802.11b CCK Matlab model less accurate than GSL model");
117 double a1 = 7.9056742265333456e-003;
118 double a2 = -1.8397449399176360e-001;
119 double a3 = 1.0740689468707241e+000;
120 double a4 = 1.0523316904502553e+000;
121 double a5 = 3.0552298746496687e-001;
122 double a6 = 2.2032715128698435e+000;
123 ber = (a1 * sinr * sinr + a2 * sinr + a3) /
124 (sinr * sinr * sinr + a4 * sinr * sinr + a5 * sinr + a6);
126 return std::min(1.0, std::pow((1.0 - ber),
static_cast<double>(nbits)));
132IntegralFunction(
double x,
void* params)
134 double beta = ((FunctionParameters*)params)->beta;
135 double n = ((FunctionParameters*)params)->n;
136 double IntegralFunction = std::pow(2 * gsl_cdf_ugaussian_P(x + beta) - 1, n - 1) *
137 std::exp(-x * x / 2.0) / std::sqrt(2.0 * M_PI);
138 return IntegralFunction;
142DsssErrorRateModel::SymbolErrorProb16Cck(
double e2)
147 FunctionParameters params;
148 params.beta = std::sqrt(2.0 * e2);
151 gsl_integration_workspace* w = gsl_integration_workspace_alloc(1000);
154 F.function = &IntegralFunction;
157 gsl_integration_qagiu(&F, -params.beta, 0, 1e-7, 1000, w, &sep, &error);
158 gsl_integration_workspace_free(w);
168DsssErrorRateModel::SymbolErrorProb256Cck(
double e1)
170 return 1.0 - std::pow(1.0 - SymbolErrorProb16Cck(e1 / 2.0), 2.0);
static double GetDsssDqpskSuccessRate(double sinr, uint64_t nbits)
Return the chunk success rate of the differential encoded QPSK.
static const double WLAN_SIR_PERFECT
WLAN perfect.
static double GetDsssDbpskSuccessRate(double sinr, uint64_t nbits)
Return the chunk success rate of the differential BPSK.
static double GetDsssDqpskCck5_5SuccessRate(double sinr, uint64_t nbits)
Return the chunk success rate of the differential encoded QPSK for 5.5Mbps data rate.
static double GetDsssDqpskCck11SuccessRate(double sinr, uint64_t nbits)
Return the chunk success rate of the differential encoded QPSK for 11Mbps data rate.
static const double WLAN_SIR_IMPOSSIBLE
WLAN impossible.
static double DqpskFunction(double x)
A function DQPSK.
#define NS_LOG_COMPONENT_DEFINE(name)
Define a Log component with a specific name.
#define NS_LOG_DEBUG(msg)
Use NS_LOG to output a message of level LOG_DEBUG.
#define NS_LOG_FUNCTION_NOARGS()
Output the name of the function.
#define NS_LOG_WARN(msg)
Use NS_LOG to output a message of level LOG_WARN.
Every class exported by the ns3 library is enclosed in the ns3 namespace.