R语言最大流最小割定理和最短路径算法分析交通网络流量拥堵问题

今天早上,我们使用一些论文中提到的示例,使用最大流最小割定理将流量拥塞降至最低, 应用了最短路径分析了交通瓶颈

由Kaizong Ye,Liao Bao撰写

我们可以在下面看到

 
map=openp(map)
plot(map)
points(t(m[3:2,]),col="black", pch=19, cex=3
×

割(CUT)是网络中顶点的划分,它把网络中的所有顶点划分成两个顶点的集合源点S和汇点T。记为CUT(S,T)。

如下图:源点:s=1;汇点:t=5。框外是容量,框内是流量

 

如下图是一个图的割。顶点集合S={1,2,3}和T={4,5}构成一个割。

           

如果一条弧的两个顶点分别属于顶点集S和T那么这条弧称为割CUT(S,T)的一条割边。  从S指向T的割边是正向割边,从T指向S的割边是逆向割边。如上图正向割边:

1->2;3->5      逆向割边:2->3。

割CUT(S,T)中所有正向割边的容量和称为割的容量。不同的个的容量不同。如上图割的容量为4+4=8;割的正向流量:4+2=6    逆向割的流量:1。

定理一:如果f是网络中的一个流,CUT(S,T)是任意一个割,那么流量f的值等于正向割边的流量与负向割边的流量之差。

推理提示:如果V既不是源点也不是汇点,那么:f({V},S∪T)-f(S∪T,{V})=0; 任何一个点,流入的与流出的量相等。

        结论:f= f(S,T)- f(T,S)<=f(S,T)<=割CUT(S,T)的容量 。 

推论1:如果f是网络中的一个流,CUT(S,T)是一个割,那么f的值不超过割CUT(S,T)的容量。

推论2:网络中的最大流不超过任何割的容量 


定理2: 在任何网络中,如果f是一个流,CUT(S,T)是一个割,且f的值等于割CUT(S,T)的容量,那么f是一个最大流,CUT(S,T)是一个最小割(容量最小的割)。



定理3:最大流最小割定量: 在任何的网络中,最大流的值等于最小割的容量。

结论1:最大流时,最小割cut(S,T)中,正向割边的流量=容量,逆向割边的流量为0。否则还可以增广。

结论2:在最小割中cut(S,T)中:
 
                        ① 源点s∈S。
 
                        ② 如果i∈S,结点j满足:有弧<i,j>,并且c[I,j]>f[I,j]  或者有弧<j,i>并且f[j,i]>0,那么j∈S。//否则不是最小割
 
                       即从s出发能找到的含有残留的点组成集合S。其余的点组成集合T。



提取有关边缘容量

要提取有关边缘容量的信息,在该网络上使用以下代码,该代码将从论文中提取三个表

 extract_tab(location)

在Windows中,要先下载另一个软件包


课程

R语言数据分析挖掘必知必会

从数据获取和清理开始,有目的的进行探索性分析与可视化。让数据从生涩的资料,摇身成为有温度的故事。

立即参加

library(devtools)
 
extract_tab(locatio

现在我们可以得出具有容量的数据框

B1=as.data.frame(out[[2]])
B2=as.data.frame(out[[3
 
capacity=as.character(B2$V3[-1])
capacity[6]="843"
ic(capacity)

我们可以在地图上添加这些边

plot(map)
points(t(m[3:2,]),col="black", pch=1
 
for(i in 1:nrow(E)){
i1=which(B$i==as.character(E$from
]))
segments(B[i1,"x"],B[i1,"y"],B[i2,
 
text(t(m[3:2,]),c("s",1:10,"t"),col="white")

要获得具有容量的图形,可以使用另一种方法

有容量的网络

为了更好地了解道路通行能力,使用

plot(g, layout=as.matrix(B[,c("x","y")]),
edge.width=E$capacity/200)

通过具有容量的网络,目标是确定该网络上从源到宿的最大流量。可以使用R

$value
[1] 2571
 
$flow
[1] 10 142 130 23 0 2

我们的最大流量为2571,这与两篇论文中的最大流量最小割定理以及 最短路径的应用中都实际要求的不同   ,因为表格和图表上的值不同。

E$flux1=m$flow
plot(g, layout=as.matrix(B[,c("x","y")]),

考虑采用更简单的流程,但是相同的全局值

 
 
E(g)$label=E$flux2
plot(g, layout=as.matrix(B[,c("x","y")]),
edge.width=E$flux2/200)

实际上,有可能在同一城市的另一篇论文中做同样的事情,这是道路网络的交通拥堵问题

 
 
dim(out[[3]])
B1=a
ame(from=B1[2:61,"V2"],
to=B1[2:6
as.numeric(
as.characte
data_frame(E)
m=max_flow(graph=g,
source="S",
 
E$flux1=m$flow
E(g)$label=E
 
edge.width=E$flux1/200,
edge.arrow.size=0.15)

此处的最大流量值为4017,就像原始论文中发现的那样


可下载资源

关于作者

Kaizong Ye拓端研究室(TRL)的研究员。

本文借鉴了作者最近为《R语言数据分析挖掘必知必会 》课堂做的准备。

​非常感谢您阅读本文,如需帮助请联系我们!

 
QQ在线咨询
售前咨询热线
15121130882
售后咨询热线
0571-63341498

关注有关新文章的微信公众号


永远不要错过任何见解。当新文章发表时,我们会通过微信公众号向您推送。

技术干货

最新洞察

This will close in 0 seconds