从Bianchi模型到实践：深入解析IEEE 802.11 DCF饱和吞吐量计算

1次阅读

共计 3013 个字符，预计需要花费 8 分钟才能阅读完成。

在无线局域网（WLAN）中，IEEE 802.11 标准定义了分布式协调功能（DCF）作为介质访问控制（MAC）层的基础协议。DCF 采用载波侦听多路访问 / 冲突避免（CSMA/CA）机制，通过随机退避来协调多个节点共享无线信道。理解 DCF 的性能特征，尤其是饱和吞吐量（即网络在重负载条件下的最大吞吐量），对于网络协议开发者至关重要。

Bianchi 在 2000 年提出的马尔可夫模型是分析 DCF 饱和吞吐量的经典方法。该模型通过数学建模，量化了竞争节点数、冲突概率和吞吐量之间的关系。本文将详细解析 Bianchi 模型的理论基础，并提供完整的 Python 实现，帮助开发者掌握这一核心分析方法。

Bianchi 模型基于三个核心假设：

饱和条件 ：每个节点始终有数据帧等待发送，不考虑空队列情况。
理想信道 ：忽略信道错误，仅考虑因冲突导致的传输失败。
固定节点数 ：竞争节点数在分析过程中保持不变。

模型将节点的退避过程建模为二维马尔可夫链，其中状态由退避阶段 $i$ 和退避计数器 $k$ 决定。稳态概率 $b_{i,k}$ 表示节点处于状态 $(i,k)$ 的概率。通过求解平衡方程，可以得到稳态概率的闭式表达式：

$$
\tau = \frac{2(1-2p)}{(1-2p)(W+1) + pW(1-(2p)^m)}
$$

其中，$\tau$ 为节点在随机时隙发送帧的概率，$p$ 为冲突概率，$W$ 为最小竞争窗口大小，$m$ 为最大退避阶段。

饱和吞吐量 $S$ 定义为成功传输的负载与总时间的比值：

$$
S = \frac{P_s P_{tr} E[P]}{(1-P_{tr})\sigma + P_{tr} P_s T_s + P_{tr} (1-P_s) T_c}
$$

其中：
– $P_{tr} = 1 – (1-\tau)^n$：至少一个节点发送的概率
– $P_s = \frac{n\tau(1-\tau)^{n-1}}{P_{tr}}$：发送成功的条件概率
– $E[P]$：帧平均长度
– $T_s$, $T_c$：成功传输和冲突的持续时间
– $\sigma$：时隙长度

以下是基于 Python 的 Bianchi 模型实现，完整代码包含冲突概率计算、稳态概率求解和吞吐量计算。

import numpy as np
import matplotlib.pyplot as plt

def bianchi_throughput(n, W=32, m=5, payload=1500, ack=14, 
                      slot_time=9e-6, SIFS=16e-6, DIFS=34e-6):
    """
    计算 IEEE 802.11 DCF 饱和吞吐量

    参数:
        n: 竞争节点数
        W: 最小竞争窗口大小 (默认 32)
        m: 最大退避阶段 (默认 5)
        payload: 数据帧长度 (字节)
        ack: ACK 帧长度 (字节)
        slot_time: 时隙时间 (秒)
        SIFS: SIFS 时间 (秒)
        DIFS: DIFS 时间 (秒)

    返回:
        饱和吞吐量 (Mbps)
    """
    # 物理层参数 (802.11a)
    PHY_header = 16  # 物理层头 (字节)
    MAC_header = 34  # MAC 层头 (字节)
    R = 54e6  # 数据速率 (bps)
    R_control = 6e6  # 控制帧速率 (bps)

    # 帧传输时间计算
    T_data = (PHY_header + MAC_header + payload) * 8 / R
    T_ack = (PHY_header + ack) * 8 / R_control

    # 成功和冲突时间
    T_s = DIFS + T_data + SIFS + T_ack
    T_c = DIFS + T_data + SIFS

    # 迭代求解 tau 和 p
    tau = 1/(W+1)  # 初始猜测
    p = 1 - (1-tau)**(n-1)

    for _ in range(10):  # 固定迭代次数
        tau = 2*(1-2*p) / ((1-2*p)*(W+1) + p*W*(1-(2*p)**m))
        p = 1 - (1-tau)**(n-1)

    # 计算吞吐量
    P_tr = 1 - (1-tau)**n
    P_s = n*tau*(1-tau)**(n-1) / P_tr if P_tr > 0 else 0

    E_P = payload * 8  # 比特数
    S = (P_tr * P_s * E_P) / ((1-P_tr)*slot_time + P_tr*P_s*T_s + P_tr*(1-P_s)*T_c)

    return S * 1e-6  # 转换为 Mbps

# 绘制吞吐量曲线
node_counts = range(1, 50)
throughputs = [bianchi_throughput(n) for n in node_counts]

plt.figure(figsize=(10, 6))
plt.plot(node_counts, throughputs, 'b-', linewidth=2)
plt.xlabel('Number of Competing Nodes')
plt.ylabel('Saturation Throughput (Mbps)')
plt.title('IEEE 802.11 DCF Saturation Throughput (Bianchi Model)')
plt.grid(True)
plt.show()