简介面向通信工程研究者与蜂窝网络仿真学习者这份资源围绕泊松点过程随机几何与蒙特卡洛方法提供了一套蜂窝通信系统性能分析仿真代码。包内共33个文件以28个m脚本为主体覆盖PPP分布生成、用户基站空间建模、干扰与覆盖概率计算、蒙特卡洛多次随机试验及结果可视化另含3个fig图表与2张jpg示意图便于对照分析仿真输出。资源压缩包仅135KB轻量便携适合快速部署运行。已有347人学习下载可作为理解随机几何在蜂窝网络中应用的入门与进阶工具。通过运行这些脚本读者可量化分析不同基站密度下的覆盖概率、吞吐量与干扰水平探索PPP与蒙特卡洛结合的网络规划优化方法也可在此基础上扩展多径、阴影衰落等复杂无线环境研究。1. 蜂窝仿真里的PPP与蒙特卡洛为什么随机几何成了覆盖率分析的标准答案做蜂窝网络仿真的人十有八九都掉进过一个坑用规则六边形网格假设基站位置仿真出来的覆盖率高得不敢信一落地就翻车。PPP泊松点过程把基站位置当成随机几何里的随机点集蒙特卡洛方法再把这些随机实现统计成覆盖率曲线这套组合最近几年几乎成了蜂窝通信仿真里做覆盖率评估的默认做法。这篇笔记直接拆一个典型的“PPP仿真”工程该怎么组织泊松点过程怎么生成、SINR怎么算、蒙特卡洛要跑多少次、参数怎么调才不出玄学结果。适合正在做无线网络仿真、或者准备用随机几何做蜂窝覆盖研究的同学照着跑一遍就能上手。2. 泊松点过程建模从PPP定义到蜂窝仿真里的两种采样思路2.1 为什么蜂窝仿真不用六边形网格PPP的随机性来自哪里传统教材里蜂窝网络画成整齐的六边形基站间距相等、小区大小一致算出来的信干噪比SINR只跟用户落在小区哪个位置有关。问题是现实的站址从来不是这样楼顶租赁条件、地价、物业协调、既有铁塔资源每一条都让基站的最终落点偏离规划网格。你把模型建得越整齐仿真结果就越像“理想国”跟实测覆盖一对比就露馅。随机几何里的泊松点过程恰恰反过来它假设基站在二维平面上完全随机地出现每个位置出现的概率相等任意两个基站的位置互不影响。PPP的数学定义很便宜就两个要素。一个是强度 λ代表单位面积里基站的期望个数单位是“个/平方米”或者“个/平方公里”另一个是点的分布在有限区域内点数服从泊松分布每个点在该区域内均匀独立分布。生成起来也简单先按泊松分布抽一个点数 N再抽 N 个均匀分布坐标完事。密度 λ 是蜂窝仿真里最核心的一个旋钮城区宏站稀疏、微站密集反映到模型里就只是把 λ 从 1e-6 调到 1e-5 的区别。PPP在蜂窝仿真里还有一个特别好用的性质叫 Sliwnyak 定理的通俗版在 PPP 条件下观察“典型用户”给定用户位置后其余基站仍然构成一个同强度的 PPP。这意味着你把手机放在仿真区域正中央、还是放在角落看干扰的统计特性基本是一样的。蒙特卡洛实现里这意味着用户位置不必反复挪只需要保证仿真区域足够大、边界效应被控制住就能把接收机固定在中心位置跑统计。当然工程上还是建议在每次迭代里随机化用户位置原因放在第4章的避坑记录里讲。2.2 蒙特卡洛在PPP仿真里的角色从一次随机实现到覆盖率曲线有了 PPP 点过程下一步要回答的问题是“用户信干噪比超过某个门限的概率有多大”。这个概率理论上是个积分对所有可能的基站布局、所有可能的信道衰落求平均。问题在于被积函数没有解析形式——路径损耗指数 α 不是整数、加了阴影衰落、再加多天线波束赋形闭式解就彻底写不出来了。这时候蒙特卡洛登场它不碰积分直接把期望写成样本均值随机生成 N 次基站布局和信道实现每次判断用户 SINR 是否超过阈值 T覆盖率就是“超过次数 / 总次数”。通信仿真里蒙特卡洛几乎是默认做法原因就一条模型复杂度随便加。你想加对数正态阴影衰落在每次实现里乘一个 lognormal 随机数就行想加两天线接收分集把信号功率的衰落系数换成两个独立指数分布变量的均值就行想加基站负载和关断在每次迭代里按负载概率丢弃部分基站即可。这些操作对解析方法来说每一步都让公式复杂度翻倍对蒙特卡洛来说只是循环体里多一行代码。代价是收敛速度估计值的标准差按 1/√N 下降想把误差缩到原来的三分之一迭代次数要加到九倍。所以蒙特卡洛 PPP 仿真的核心矛盾不是“怎么算”而是“要算多少次才能信”这一条贯穿整篇笔记。2.3 泊松簇过程和硬核过程的取舍什么时候别用纯PPP纯 PPP 有个内在假设基站位置完全独立最小站间距可以无限小。真实网络里两个基站离个三五米的情况极少于是有人引入硬核过程Matérn 硬核或者泊松簇过程Thomas 簇过程当进阶模型。硬核过程先把 PPP 撒出来再按最小站间距删除冲突点得到的结果比纯 PPP 更均匀泊松簇过程适合模拟家庭基站或微站成簇部署的场景比如一个商场里密集放十几个小站。取舍标准其实很直接如果仿真目标是看宏站覆盖、评估干扰态势纯 PPP 足够因为它在“随机”和“可解析”之间取了最佳平衡论文里也有大量理论基准可对照如果仿真目标是微站组网、基站分簇部署带来的容量差异再考虑升级成簇过程或硬核过程。我一般建议先跑通纯 PPP 的蒙特卡洛框架把干扰统计、覆盖率计算、参数标定的流程全部打通再谈换点过程模型——框架本身不换只换点过程生成函数。这个“先框架后模型”的顺序能省掉大量返工。3. 跑通PPP蒙特卡洛仿真Python实现、SINR统计与覆盖率曲线3.1 工程结构与依赖一个可复现的PPP仿真工程怎么组织拿到一个名为“PPP仿真”的工程包最常见的组织方式是四个脚本加一个配置文件点过程生成、信道与SINR计算、蒙特卡洛主循环、结果绘图配置独立成文件。这样拆的核心目的是把随机性来源隔离——PPP点集生成是一类随机快衰落是另一类随机两种随机混在一个文件里后面想固定种子复现时会很痛苦。工程里各文件的职责大致如下表。文件职责关键接口ppp_gen.py泊松点过程生成generate_ppp(lambda_b, L, rng)channel.py路径损耗、瑞利衰落、SINR计算compute_sinr(user, bs, alpha, noise_power, rng)monte_carlo.py蒙特卡洛主循环与收敛判断monte_carlo_coverage(...)config.py所有仿真参数集中管理常量定义禁止在其它文件里改参数plot_results.py覆盖率曲线与收敛曲线绘图输出 SVG/PNG依赖只有三个numpy 负责随机数与向量运算scipy 可选用于特殊统计检验matplotlib 出图。Python 3.8 以上就行不需要 GPU也不需要任何通信仿真专用库——这个方案的全部优势就在于它足够裸每一行都在你眼皮底下。配置集中管理这条我后面还会强调因为参数散落在各文件里是蒙特卡洛仿真翻车的第一大源头。3.2 核心代码拆解PPP生成、SINR计算与蒙特卡洛主循环第一步是生成泊松点过程。实现就三行但里面藏着 PPP 的关键假设点数服从泊松分布位置服从均匀分布。import numpy as np def generate_ppp(lambda_b, L, rng): 在 [-L/2, L/2] x [-L/2, L/2] 区域内生成强度为 lambda_b 的泊松点过程。 lambda_b: 基站强度单位 每平方米 L: 仿真区域边长单位 米 rng: numpy 随机数生成器用于复现 area L * L n_points rng.poisson(lambda_b * area) # 期望点数 强度 x 面积 points rng.uniform(-L / 2, L / 2, size(n_points, 2)) return points逻辑说明n_points 是每次实现里基站的数量它本身是个随机变量均值是 λL²方差也是 λL²这是泊松分布的特性不要手滑写成固定 int(lambda_b * area)那会把 PPP 变成确定性网格。points 的每一行是一个基站的 (x, y) 坐标。rng 必须由外部传入而不是在函数内部新建否则无法固定随机种子。参数方面L 的选择直接决定边界效应强弱经验值是最小取 5 倍“最近基站期望距离”λ1e-5 时这个距离约 158 米L 至少取 800 米λ 越密L 可以适当缩小但要按公式走。第二步是核心的 SINR 计算。信号功率取最近基站的路径增益干扰功率累加其余所有基站快衰落按瑞利模型独立生成。def compute_sinr(user, bs, alpha, noise_power, rng, rayleighTrue): user: 用户坐标形如 (x, y) bs: 基站坐标矩阵形状 (N, 2) alpha: 路径损耗指数典型值 3~4 noise_power: 噪声功率线性值非 dB rng: 随机数生成器 dist np.linalg.norm(bs - user, axis1) dist np.maximum(dist, 1e-6) # 防止基站与用户重合导致除零 path_gain dist ** (-alpha) # 路径损耗转成功率增益 serve_idx np.argmax(path_gain) # 最近基站 最大路径增益 signal_power path_gain[serve_idx] if rayleigh: # 瑞利衰落指数分布均值固定为1 # 信号与干扰必须使用独立衰落系数否则相关性会虚增覆盖率 signal_power * rng.exponential(1.0) interf_power 0.0 for i in range(len(path_gain)): if i ! serve_idx: interf_power path_gain[i] * rng.exponential(1.0) else: interf_power path_gain.sum() - signal_power # 无衰落纯路径损耗 sinr signal_power / (interf_power noise_power) return sinr逻辑说明argmax 选最近基站是一个关键决策。实测网络里用户不一定总是接入最近基站——还可能接入信号最强、负载最轻或偏置后的邻站——但对于覆盖率基准仿真“最近基站”是最干净也最保守的假设改接入规则不会影响框架结构只会改变 serve_idx 的选择逻辑。干扰累加里特别注意信号衰落和干扰衰落必须是独立抽取的指数随机变量互不相关。rayleighFalse 的纯路径损耗模式适合先做调试跑通后再开衰落。第三步是蒙特卡洛主循环。每次迭代生成一份新的 PPP 布局和一组新的快衰落实现记录阈值判断结果。def monte_carlo_coverage(lambda_b, L, alpha, noise_power, threshold_db, n_iter, seed42, rayleighTrue): rng np.random.default_rng(seed) threshold 10 ** (threshold_db / 10) # 阈值从 dB 转线性值 coverage_count 0 for _ in range(n_iter): bs generate_ppp(lambda_b, L, rng) user np.zeros(2) # 典型用户固定在中心 sinr compute_sinr(user, bs, alpha, noise_power, rng, rayleigh) if sinr threshold: coverage_count 1 return coverage_count / n_iter逻辑说明每次迭代里基站布局和快衰落全部重新抽取保证样本独立。coverage_count / n_iter 就是覆盖率的蒙特卡洛估计值大数定律保证它随 n_iter 增大收敛到真实覆盖概率。seed 固定后同一份代码跑出来的曲线完全一致这在调参和排错时是底线保障。threshold_db 默认给一个单阈值标量实际画曲线时要在一系列阈值上重复整个蒙特卡洛过程也就是把 for threshold in [-10, -5, 0, 5, 10, 15, 20] 包在外层。3.3 结果怎么画才有效SINR互补CDF曲线是标准输出单点覆盖率数值没有说服力行业里默认画的是 SINR 互补 CDF 曲线横轴是 SINR 阈值dB纵轴是“SINR 超过该阈值的概率”。这条曲线从左到右单调下降横轴取 -10 dB 到 20 dB纵轴从接近 1 掉到接近 0。画曲线时每条配置一条线比如密集城区 λ1e-5 一条、郊区 λ1e-6 一条线之间的间距就是密度带来的覆盖差异。绘图代码不复杂但要输出矢量图而不是截图放进论文。plt.plot(thresholds_db, coverage_list, markero)plt.xlabel(SINR threshold (dB))plt.ylabel(P(SINR threshold))plt.grid(True, whichboth, alpha0.3)。需要注意纵轴单位不是百分比而是概率横轴必须是 dB 值而非线性值。另外每条曲线至少要跑三组随机种子取平均画误差带否则曲线上的“小锯齿”会被审稿人追问。4. PPP蒙特卡洛仿真的避坑记录边界、衰落与种子的4个教训4.1 覆盖率总比论文高一截边界效应把干扰算少了现象仿真跑出来的覆盖率明显高于同参数下的论文曲线或实测参考值尤其是在低 SINR 阈值区间差值能到 5~10 个百分点。原因仿真区域边长 L 取小了。用户固定在中心时理想情况是干扰来自无穷平面上的所有基站但有限区域内离用户很远的那些基站根本不在点集里干扰功率被系统性截断。干扰少了SINR 虚高覆盖率整体上移。这个问题在 λ 小、L 相对更小时尤其严重。解决把 L 设为最近基站期望距离的 5 倍以上。PPP 下典型用户到最近基站的期望距离是 1/(2√λ)L 取这个值的 5 倍时边界截断造成的干扰误差已经低于蒙特卡洛本身的统计波动。另一种做法是周期边界wrap-around把有限区域折成环面点集在边界处左右翻转对齐用户收到来自“镜像”基站的干扰这样做的好处是效率高缺点是区域太小时会出现同一个基站以多个镜像身份重复加干扰的伪影。先加 L再考虑 wrap-around。4.2 加了瑞利快衰落覆盖率掉得离谱信号与干扰的衰落系数没有独立现象从纯路径损耗模式切换到瑞利衰落模式后覆盖率曲线掉得特别猛中高阈值区间甚至腰斩怎么调都回不去。原因两种典型的手误。第一种是把衰落系数的均值写错比如 rng.exponential(0.5)指数分布的均值正好是 scale 参数写 0.5 等于平均把功率砍了一半第二种是把同一个衰落系数同时乘到信号和干扰上这样信号与干扰完全相关SINR 的分子分母被同一个随机数缩放方差被严重低估高 SINR 区间的概率被抬高或压低都不自知。解决衰落系数统一写成 rng.exponential(1.0)均值保真为 1。更重要的是信号衰落和每个干扰基站的衰落都必须独立抽取不能复用同一个随机变量。可以用一个快速自检生成 10 万次指数随机数打印均值确认接近 1.0再分别在 rayleighTrue 和 rayleighFalse 下跑同一组参数确认衰落只会让曲线的中段变缓、低阈值段略升不会让整条曲线的形状面目全非。4.3 换了机器曲线对不上伪随机种子没固定现象同一份代码在自己笔记本上跑出的覆盖率曲线到服务器上重跑变成另一条完全不同的曲线。更尴尬的是改了一个参数想对比差异结果曲线整体漂移根本看不出是参数变化引起的还是随机波动引起的。原因蒙特卡洛结果本质上是随机数序列的函数。numpy 默认每次运行时从系统熵池取种子两次运行的随机序列完全不同点过程布局和衰落系数全部换了曲线自然对不上。解决在 monte_carlo.py 入口处固定 rng np.random.default_rng(42)并把 seed 作为参数暴露在 config.py 里。想对比两组参数时务必用同一个 seed想换一组随机性验证稳健性时用 seed43、seed44 各跑一遍。更稳妥的做法是把每次迭代的 PPP 点集用 np.save 存盘这样即使随机数版本升级导致序列变化numpy 2.x 和 1.x 的 default_rng 算法保持兼容但也保不齐后续变化点集还在历史结果永远可复现。4.4 迭代次数加了一倍曲线还在抖用户位置没参与随机化现象把 n_iter 从 1000 加到 5000覆盖率曲线还是明显抖动滑动平均的波动幅度没有随迭代次数下降好像钱花在了刀刃以外的地方。原因每次迭代只重新抽取基站布局和衰落用户位置固定在中心零点。虽然 PPP 的平移不变性让中心用户和任意位置用户的干扰统计理论上相同但单次实现的有限区域边界对干扰截断的影响与用户位置有关固定中心等于把“用户位置”这个随机维度完全丢掉了方差被低估收敛速度反而变慢。解决每次迭代时在仿真区域内随机化用户位置让用户位置和基站布局互相独立。更标准的做法是借助 Sliwnyak 定理在每次迭代里把用户放在区域中心、然后把整个点集做一个随机的整体平移等效于随机化用户位置而不改变 PPP 的分布属性。实现上就是在 generate_ppp 之后给所有基站坐标加一个随机偏移量。这个改动的成本是零收益是覆盖率曲线更早收敛第 5 章的收敛判断也更可信。5. 把仿真调到可信λ、路径损耗指数、噪声功率与迭代次数的联动5.1 四个必调参数含义、默认值与曲线走向蒙特卡洛 PPP 仿真里参数不少但真正决定曲线形状的只有四个基站强度 λ、路径损耗指数 α、噪声功率 σ²、迭代次数 n_iter。仿真区域 L 和阈值 T 属于“必须合理但不必细调”的参数。每个参数调大以后覆盖率曲线的走向是调参时最该背下来的东西。参数典型范围调大后的影响备注λ每平方米1e-6 ~ 1e-5干扰受限场景下覆盖率几乎不变噪声受限场景下低阈值段明显上升λ 翻倍信号与干扰同步增强α无量纲3.0 ~ 4.0覆盖率整体上升尤其在干扰受限区α 越大远端干扰衰减越快σ²线性功率1e-13 ~ 1e-12低阈值段下降高阈值段几乎不动对应 -130 dBm 量级的接收机底噪n_iter次数500 ~ 5000曲线波动变小估值更稳收敛后继续加次数收益递减λ 的影响最反直觉在无噪声的纯干扰场景λ 变大以后最近基站距离变小、信号变强但其余基站也变近、干扰同步变大SINR 的分布整体不变。这个结论在随机几何里有严格推导工程上意味着如果你仿真的场景是“密集城区整体加密基站”提升的只是覆盖概率中的“被服务”部分而不是“SINR 变好”的部分真正让 SINR 变好的是α、是加装基站时的功率偏置而不是无脑加密。α 的影响相对直白路径损耗指数变大远处干扰衰减更快信号从最近的基站来、距离损失小SINR 自然受益。城市峡谷里 α3.5、开阔区域 α4.0这是仿真里常见的取值。噪声功率 σ² 的标定要结合实际接收机灵敏度。灵敏度 -100 dBm、系统带宽 20 MHz 时底噪约 -120 dBm 量级但仿真里用的是线性功率要把 dBm 转成瓦特再参与运算。很多人在这一步翻车直接用 dBm 数值当线性功率相加结果噪声功率比信号功率大十几个数量级覆盖率直接归零。参数配置里写明“σ² 是线性值非 dB”并在代码里加一行 assert noise_power 0。5.2 用收敛曲线决定迭代次数别等跑到天亮蒙特卡洛仿真的迭代次数不是拍脑袋定的也不是越大越好。常用的做法是画“覆盖率随迭代次数的滑动平均”曲线观察它是否进入平稳区。这里给一个带提前终止的收敛判断实现。def run_until_converge(lambda_b, L, alpha, noise_power, threshold_db, seed42, window200, tol0.005, max_iter5000): rng np.random.default_rng(seed) threshold 10 ** (threshold_db / 10) hits [] for k in range(1, max_iter 1): bs generate_ppp(lambda_b, L, rng) user np.zeros(2) sinr compute_sinr(user, bs, alpha, noise_power, rng, rayleighTrue) hits.append(1 if sinr threshold else 0) if k window: recent np.mean(hits[-window:]) std np.std(hits[-window:]) if std tol: print(fconverged at {k} iters, coverage{recent:.4f}) return recent, k return np.mean(hits), max_iter逻辑说明hits 列表里每个元素是单次迭代的 0/1 判断结果window 是滑动窗口长度tol 是窗口内标准差的阈值。标准差低于 tol 时认为估计值已经平稳。窗口长度取 200、tol 取 0.005 是一个保守但可靠的组合对应覆盖率约 50% 时标准误约 0.035窗口 200 会把这个量级压到 0.005 附近说明曲线确实进入平台期。参数说明tol 不是越小越好取 0.001 可能让程序跑到 max_iter 还不触发提前终止浪费时间取 0.01 又太早停曲线可能还在缓慢漂移。调试期建议把 window 调小到 100、tol 调大到 0.01快速看到趋势最终出图时再换保守参数。另外这个收敛判断针对的是单阈值点的收敛画整条曲线时每个阈值点都要独立跑收敛。5.3 参数组合的常见翻车场景干扰受限与噪声受限判断错了方向调参时最怕的场景是λ 翻了两倍覆盖率曲线纹丝不动于是怀疑代码有 bug。第 5.1 节已经说了这不是 bug而是场景进入了干扰受限区。反过来在城市边缘、基站稀疏的郊区场景λ 变大时低阈值段覆盖率明显上升这是噪声受限区增加基站有效。怎么快速判断当前仿真在哪个区做法是做一个“无噪声极限”实验把 noise_power 设为 0重跑覆盖率。如果结果跟有噪声时的曲线几乎重合说明噪声对 SINR 的贡献可以忽略场景是干扰受限如果低阈值段明显上移说明噪声原来压低了弱信号用户的 SINR场景是噪声受限。这个实验只改一个参数、跑一遍蒙特卡洛十秒出结果却能让后续所有调参都踩在正确的方向上。还有一个常被忽略的参数联动α 和 λ。α 变大时干扰受限区会更早出现因为远端干扰的衰减抵消了基站密度上升带来的干扰增强。也就是说在城市峡谷 α4 的场景里加密基站覆盖率收益比在开阔地 α3 的场景更小。想验证这种联动是否合理可以把两组参数跑出来叠在同一张图里曲线的形状差异一目了然也方便向同事或导师解释为什么“基站加密了但覆盖率没涨”。6. 进阶用法用对偶变量压方差让覆盖率曲线平滑一个数量级蒙特卡洛覆盖率估计的方差主要来自两块PPP 点过程的空间随机性和快衰落的乘性随机性。想在不增加迭代次数的前提下把曲线跑平滑一个非常实用的技巧是“对偶变量”。思路特别简单同一次迭代里用同一份 PPP 点集生成两套独立的快衰落实现分别计算 SINR 并做阈值判断取两次判断的平均作为这一个样本点的值。两套随机数天然负相关——一次衰落让 SINR 偏高时另一次往往偏低——平均之后方差被显著压低。def coverage_antithetic(bs, user, alpha, noise_power, threshold_lin, rng, num_antithetic2): # 同一份基站布局生成多套独立快衰落取平均降低方差 decisions [] for _ in range(num_antithetic): sinr compute_sinr(user, bs, alpha, noise_power, rng, rayleighTrue) decisions.append(1 if sinr threshold_lin else 0) return np.mean(decisions)逻辑说明compute_sinr 内部会独立抽取衰落系数所以循环里连续调用两次得到的是两份条件独立样本。把两次判断取平均单次迭代的输出不再是 0/1而是 0/0.5/1 三档方差比原来的伯努利变量小。实际跑下来同样的 2000 次迭代覆盖率曲线的高频抖动能明显减小等效于迭代次数翻倍却只花了一半的运行时间。更激进的做法是用逆变换采样构造严格负相关的衰落系数对比如对指数分布用 U 和 1-U 分别做逆变换但工程上两次独立采样的收益已经足够没必要把简单代码复杂化。另一个值得做的验证是和理论基准比对。在无噪声、纯干扰受限、瑞利衰落的条件下PPP 蜂窝网络的覆盖率存在近似解析解随机几何教材里都有这个闭式。把蒙特卡洛曲线和理论曲线叠在一起画如果中段偏离超过 2~3 个百分点几乎可以断定是代码有系统性偏差——常见偏差源就是第 4 章的边界效应或衰落均值没标定对。这个交叉验证比任何“跑更多迭代”都更能说明仿真可信。最后说个我自己踩过的教训早先跑仿真一看到曲线抖就先加迭代次数从 1000 加到 5000再不行加 10000跑一次要等十几分钟结果曲线还是带毛刺。后来把种子固定、用户位置随机化、加上对偶变量三件事一起做了同样 2000 次迭代的曲线平滑程度远超从前。从那以后我养成了习惯先固定种子和随机维度再谈加大次数先做无噪声极限实验再调具体参数。这个顺序能帮每一个刚上手 PPP 蒙特卡洛仿真的人少走大半个晚上的弯路。希望帮到你。本文还有配套的精品资源点击获取