40 return ((std::sqrt(2.0) + 1.0) / std::sqrt(8.0 * M_PI * std::sqrt(2.0))) *
41 (1.0 / std::sqrt(x)) * std::exp(-(2.0 - std::sqrt(2.0)) * x);
48 double EbN0 = sinr * 22000000.0 / 1000000.0;
49 double ber = 0.5 * std::exp(-EbN0);
50 return std::pow((1.0 - ber),
static_cast<double>(nbits));
57 double EbN0 = sinr * 22000000.0 / 1000000.0 / 2.0;
59 return std::pow((1.0 - ber),
static_cast<double>(nbits));
68 double EbN0 = sinr * 22000000.0 / 1375000.0 / 4.0;
69 double sep = SymbolErrorProb16Cck(4.0 * EbN0 / 2.0);
70 return std::min(1.0, std::pow(1.0 - sep, nbits / 4.0));
72 NS_LOG_WARN(
"Running a 802.11b CCK Matlab model less accurate than GSL model");
86 double a1 = 5.3681634344056195e-001;
87 double a2 = 3.3092430025608586e-003;
88 double a3 = 4.1654372361004000e-001;
89 double a4 = 1.0288981434358866e+000;
90 ber = a1 * std::exp(-std::pow((sinr - a2) / a3, a4));
92 return std::min(1.0, std::pow((1.0 - ber),
static_cast<double>(nbits)));
103 double EbN0 = sinr * 22000000.0 / 1375000.0 / 8.0;
104 double sep = SymbolErrorProb256Cck(8.0 * EbN0 / 2.0);
105 return std::min(1.0, std::pow(1.0 - sep, nbits / 8.0));
107 NS_LOG_WARN(
"Running a 802.11b CCK Matlab model less accurate than GSL model");
121 double a1 = 7.9056742265333456e-003;
122 double a2 = -1.8397449399176360e-001;
123 double a3 = 1.0740689468707241e+000;
124 double a4 = 1.0523316904502553e+000;
125 double a5 = 3.0552298746496687e-001;
126 double a6 = 2.2032715128698435e+000;
127 ber = (a1 * sinr * sinr + a2 * sinr + a3) /
128 (sinr * sinr * sinr + a4 * sinr * sinr + a5 * sinr + a6);
130 return std::min(1.0, std::pow((1.0 - ber),
static_cast<double>(nbits)));
136IntegralFunction(
double x,
void* params)
138 double beta = ((FunctionParameters*)params)->beta;
139 double n = ((FunctionParameters*)params)->n;
140 double IntegralFunction = std::pow(2 * gsl_cdf_ugaussian_P(x + beta) - 1, n - 1) *
141 std::exp(-x * x / 2.0) / std::sqrt(2.0 * M_PI);
142 return IntegralFunction;
146DsssErrorRateModel::SymbolErrorProb16Cck(
double e2)
151 FunctionParameters params;
152 params.beta = std::sqrt(2.0 * e2);
155 gsl_integration_workspace* w = gsl_integration_workspace_alloc(1000);
158 F.function = &IntegralFunction;
161 gsl_integration_qagiu(&F, -params.beta, 0, 1e-7, 1000, w, &sep, &error);
162 gsl_integration_workspace_free(w);
172DsssErrorRateModel::SymbolErrorProb256Cck(
double e1)
174 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.