DiggleGatesPointProcess[μ,ρ,d]
表示一个 Diggle–Gates 点过程,其固定强度为 μ 且交互半径为
中的 ρ.
DiggleGatesPointProcess
DiggleGatesPointProcess[μ,ρ,d]
表示一个 Diggle–Gates 点过程,其固定强度为 μ 且交互半径为
中的 ρ.
更多信息
- DiggleGatesPointProcess 模拟点布局,其中的点在相互距离半径 ρ 的范围内有一个软核递减互斥配对作用,否则均匀分布.
-

- Diggle–Gates 点过程可定义为关于其强度 μ 和配对势能
或配对作用
的 GibbsPointProcess,皆可通过 ρ 按照下列方式参数化: -
![phi(r)= 2 log(sin(pi r/(2 rho))) r<rho; 1 TemplateBox[{True, paclet:ref/True}, RefLink, BaseStyle -> {2ColumnTableMod}]; phi(r)= 2 log(sin(pi r/(2 rho))) r<rho; 1 TemplateBox[{True, paclet:ref/True}, RefLink, BaseStyle -> {2ColumnTableMod}];](Files/DiggleGatesPointProcess.zh/1.png)
配对势能 ![h(r)= sin^2(pi r/(2 rho)) r<rho; 1 TemplateBox[{True, paclet:ref/True}, RefLink, BaseStyle -> {2ColumnTableMod}]; h(r)= sin^2(pi r/(2 rho)) r<rho; 1 TemplateBox[{True, paclet:ref/True}, RefLink, BaseStyle -> {2ColumnTableMod}];](Files/DiggleGatesPointProcess.zh/2.png)
配对作用 - 在区域 reg 中,Diggle–Gates 点过程 DiggleGatesPointProcess[μ,ρ,d] 中的点布局
的关于 PoissonPointProcess[1,d] 的密度函数
与
成比例. - 将点
添加到点布局
的 Papangelou 条件密度
为
. - DiggleGatesPointProcess 允许 μ 和 ρ 为任意正数,且 d 为任意正整数.
- DiggleGatesPointProcess 是 GibbsPointProcess 的特殊情况.
- 当
时,DiggleGatesPointProcess 简化为 PoissonPointProcess. - RandomPointConfiguration 中,StraussPointProcess 的可能 Method 设置为:
-
"MCMC" 马尔科夫链蒙特卡洛方法初始和消亡 "Exact" 从过去耦合 - EstimatedPointProcess 中,DiggleGatesPointProcess 可能的 PointProcessEstimator 设置为:
-
Automatic 自动选择参数估计量 "MaximumPseudoLikelihood" 最大化伪似然度 - DiggleGatesPointProcess 可与诸如 RipleyK 和 RandomPointConfiguration 这样的函数一起使用.
范例
打开所有单元 关闭所有单元基本范例 (2)
domain = Rectangle[{0, 0}, {10, 10}];pts = RandomPointConfiguration[DiggleGatesPointProcess[10, 2, 2], domain]Show[RegionPlot[domain], ListPlot[pts]]reg = Entity["Country", "Denmark"]["Polygon"]pts = RandomPointConfiguration[DiggleGatesPointProcess[Quantity[3/1000, 1/"Kilometers"^2], Quantity[3, "Kilometers"], 2], reg, 1]GeoListPlot[pts]范围 (2)
proc = DiggleGatesPointProcess[50, .1, 2];
reg = Disk[];SeedRandom[1];
pts = RandomPointConfiguration[proc, reg, 2]ListPlot[pts]Clear[μ, ρ, d];
EstimatedPointProcess[pts, DiggleGatesPointProcess[μ, ρ, d]]生成两个在地球表面的 Diggle–Gates 点过程的实现过程:
μ = Quantity[.3, "Kilometers" ^ -2];
ρ = Quantity[10., "Kilometers"];proc = DiggleGatesPointProcess[μ, ρ, 2];
reg = GeoDisk[Entity["City", {"BuenosAires", "BuenosAires", "Argentina"}], Quantity[30, "Kilometers"]];pts = RandomPointConfiguration[proc, reg, 2]GeoListPlot[pts]Clear[μ, ρ, d];
EstimatedPointProcess[pts, DiggleGatesPointProcess[μ, ρ, d]]选项 (3)
Method (3)
proc = DiggleGatesPointProcess[10, .1, 2];
reg = Disk[];RandomPointConfiguration[proc, reg, Method -> "MCMC"]RandomPointConfiguration[proc, reg, Method -> {"MCMC", MaxRecursion -> 6}]RandomPointConfiguration[proc, reg, Method -> {"MCMC", "LengthOfRun" -> 5000}]proc = DiggleGatesPointProcess[10, .1, 2];
reg = Disk[];RandomPointConfiguration[proc, reg, Method -> {"MCMC", "InitialState" -> RandomPoint[reg, 100]}]proc = DiggleGatesPointProcess[10, .1, 2];reg = Rectangle[{0, 0}, {10, 10}];pts = RandomPointConfiguration[proc, reg, Method -> {"Exact", MaxRecursion -> 20}]Show[RegionPlot[reg], ListPlot[pts]]可能存在的问题 (1)
默认情况下,模拟会在点的数量收敛到一个稳定状态之前或达到默认迭代次数之前持续运行:
proc = DiggleGatesPointProcess[10000, .04, 2];
reg = Disk[];RandomPointConfiguration[proc, reg]RandomPointConfiguration[proc, reg, Method -> {"MCMC", MaxRecursion -> 6}]RandomPointConfiguration[proc, reg, Method -> {"MCMC", "LengthOfRun" -> 5 * 10 ^ 4}]相关指南
-
▪
- 空间点过程
文本
Wolfram Research (2020),DiggleGatesPointProcess,Wolfram 语言函数,https://reference.wolfram.com/language/ref/DiggleGatesPointProcess.html.
CMS
Wolfram 语言. 2020. "DiggleGatesPointProcess." Wolfram 语言与系统参考资料中心. Wolfram Research. https://reference.wolfram.com/language/ref/DiggleGatesPointProcess.html.
APA
Wolfram 语言. (2020). DiggleGatesPointProcess. Wolfram 语言与系统参考资料中心. 追溯自 https://reference.wolfram.com/language/ref/DiggleGatesPointProcess.html 年
BibTeX
@misc{reference.wolfram_2026_digglegatespointprocess, author="Wolfram Research", title="{DiggleGatesPointProcess}", year="2020", howpublished="\url{https://reference.wolfram.com/language/ref/DiggleGatesPointProcess.html}", note=[Accessed: 15-September-2026]}
BibLaTeX
@online{reference.wolfram_2026_digglegatespointprocess, organization={Wolfram Research}, title={DiggleGatesPointProcess}, year={2020}, url={https://reference.wolfram.com/language/ref/DiggleGatesPointProcess.html}, note=[Accessed: 15-September-2026]}