五一数学建模A题不确定性下的最短路径问题CUMT赖增强

2015年暑期数学建模培训第一次模拟承诺书我们仔细阅读了数学建模联赛的竞赛规则。

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与本队以外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其它公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。

如有违反竞赛规则的行为，我们愿意承担由此引起的一切后果。

我们授权数学建模联赛赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。

我们参赛选择的题号为（从A/B/C中选择一项填写）：我们的参赛报名号为：参赛组别（研究生或本科或专科）：本科所属学校（请填写完整的全名）中国矿业大学南湖校区参赛队员 (打印并签名) ：1. 赖增强2. 兰卫旗3. 李康杰日期：2015年8月11日获奖证书邮寄地址：中国矿业大学南湖校区桃4B5032邮政编码：221116收件人姓名：赖增强联系电话：2015年暑期数学建模培训第一次模拟编号专用页竞赛评阅编号（由竞赛评委会评阅前进行编号）：评阅记录裁剪线裁剪线裁剪线竞赛评阅编号（由竞赛评委会评阅前进行编号）：参赛队伍的参赛号码：1（请各参赛队提前填写好）：不确定条件下的最优路径问题摘要本文针对如何在复杂的交通环境下寻找一条可靠、快速、安全的最优路径的问题，考虑到交通堵塞、恶劣天气、路途成本等不确定因素对司机路径选择的影响，建立多个不确定条件下的最优路径模型。

对于问题一，我们在各个路段所用时间服从正态分布N(μ，δ2)的基础上,建立了在不确定条件下求最短路的NP 模型，给每个路段设定一个预留到达的时间t ，为了尽可能准确的到达目的地，选取95%的概率，满足P{T ≤t}?95%,那么最优路径的定义就是预留时间最小的那个路径，将其转换为标准的正态分布，通过标准的正态分布得到了在不确定性条件下车辆从起点到终点预留时间的数学表达式：t=μ+Φ−1δ。

计算得对应的t （绕城）=，t （市区）=，那么最优路径为绕城快速路。

对于问题二，在第一问定义的基础上进一步引入Bool 系数β(a ,k ),在搜集得到的具体的交通网络中，建立了一个从起点到终点路径为∑β(a ,k )n a =1T a link 的正态分布，通过求最小预留时间t (min)=E[T k path ]+Φ−1√Var [T k path ] ,得出最优路径的算法。

其中E [T Kpath ]=∑β(a ,k )n a =1E[T a link ]，Var [T k path ]= ∑β(a ,k )n a =1Var [T a link ]，但Var [T kpath ]的根式不具有线性可加性。

不能用经典的dijkstra 算法求解。

对此采用基于双目标规划的思路，利用第K 短路径算法，分别对E[T k path ]，Var [T kpath ],运用matlab 编程，找出各自前十条最短路径。

之后在其并集中找出最优路径：V1→V3→V4→V8。

由此建立了求最短路的NPK 模型。

最后从时间的渐进性态上分析模型的复杂性和收敛性。

对于问题三，我们只考虑各路段空间上的相关性，并用概率论中的协方差来表示这种耦合关系，建立了NPK 模型。

得出可靠时间的数学表达式t= E [T K path ]+Φ−1(ρ)√∑δ(a ,k )V ar [T a link +∑cov (a −1,a )n a =2n a =1；求解得最优路径：V1->V3->V5->V8。

关键词：NPK 模型，K-短路，预留时间，正态分布，渐进性态。

一.问题重述目前，交通拥挤和事故正越来越严重的困扰着城市交通,如何在复杂的交通环境寻找一条可靠、快速、安全的最优路径，已经成为所有驾驶员的共识。

传统的最优路径就是行驶时间最短的路径，这是基于理想交通情况下分析的，而事实上，在现实生活中，受到很多不确定性因素的影响，例如：交通工具、恶劣天气、突发事件，导致车辆的行驶时间均存在不确定性。

为了减少交通阻塞所耽误的时间，尽可能快的到达目的地，解决不确定性性条件下的最优路径问题，现依次提出以下问题：（一）针对一般的交通网络，假设已知每条路段行驶时间的均值和标准差，请建立相关的数学模型，定量地分析车辆行驶时间的不确定性，然后给出在不确定性条件下车辆从起点到终点的最优路径的定义和数学表达式。

并将此模型应用到图1的例子中会选择哪条道路。

（二）根据第一问的定义，假设已知每条路段行驶时间的均值和标准差，设计算法搜索最优路径，并将该算法应用到具体的交通网络中，用计算结果验证算法的有效性。

从理论上分析算法的收敛性、复杂性等性质。

（三）建立数学模型，描述下游路段发生交通拥堵使车辆减速或者排队，导致上游路段发生拥堵的交通路段之间驾驶时间的相关性，并将此相关性应用到第一问和第二问的最优路径搜索问题中，并设计算法解决考虑相关性的最优路径搜索问题，给出算例验证算法的有效性。

从理论上分析算法的收敛性、复杂性等问题。

二.基本假设符号说明基本假设各个路段所用时间服从正态分布。

假设仅相邻两条路段之间具有相关性。

假设任意两条相邻路段组成的协方差矩阵为一个随机的正半定矩阵。

市区道路均值30分钟，标准中国矿业大学 徐州火符号说明T n link两相邻节点之间的路段T m path从起点到终点的路径β(a,k) Bool系数。

当路径a通过此路段k时为1Var[T]对随机变量T求方差运算Φ−1标准正太分布累计概率函数的逆函数（ρ）cov(t1,t2)随机向量t1,t2协方差四.模型的建立与求解问题一问题分析本问题是对于题设的交通网络，已知每条路段行驶时间的均值μ和标准差δ条件下。

定量分析车辆行驶时间的不确定性，以及给出在不确定性条件下最优路径的定义和数学表达式。

假设各个路段所用时间服从正态分布N(μ，δ2)模型，利用MATLAB 模拟（附录一）两条路径的正态分布图[1]：（图）图给每个路段设定一个预留到达的时间t，为了尽可能准确的到达目的地，选取95%的概率，满足P{T≤t}≥95%。

那么最优路径的定义就是预留时间最小的那个路径。

模型建立与求解已知到达目的地所用时间和预留时间满足：P{T ≤t}≥95%，将其转换为标准的正态分布：P{T δμ-＜ t δμ-}≥95%，得到Φ0（t δμ-）≥95%，t δμ-≥Φ0-1（ρ），（其中ρ为95%）即可解得每条路段的t ≥μ+δ，取t=μ+δ。

我们将上诉模型称之为不确定条件下的NP 模型。

应用已知的数学表达式，将图1所给的数据：μ 1=33， δ1=1；μ 2=30， δ2=15；带入公式计算出：t （绕城）=，t （市区）=，最优路径为绕城快速路。

问题二问题分析 根据第一问中的最优路径定义和相关数学表达式，对于一般交通网络，可以结合K 短路径算法建立NPK 模型，最后从时间的渐进性态上分析其复杂性和收敛性。

模型建立与求解对于一般交通网络，为了方便设计算法找到最优路径，搜集了八个城市之间路段的时间均值和时间标准差，列表如下(表：表我们可以将其转化为图论问题。

将八个城市看做节点构成一个节点集：V={V1,V2,V3,V4,V5,V6,V7,V8}各个城市之间的道路看做边集：E={V1→V2，V1→V3，V1→4，V1→V5，V1→V6，V1→V7，V2→V3，V1→V4，V2→V7，V3→V4，V3→V5，V4→V5，V4→V8，V5→V6，V4→V8，V6→V7，V6→V8，V7→V8}则可将八个城市之间的交通网络图看做一个无向赋权图G （图.）每条路为图中的边。

八个城市之间交通网络数据图图.在第一问定义的基础上，针对每条路段引入Bool 系数β(a ,k ),当该路段被选择时为1，否则为0。

那么从起点到终点的路径可表示为∑β(a ,k )n a =1Τk path ，可知其服从正态分布。

通过求该路径的最小预留时间t (min)=E[T k path ]+Φ−1 （ρ）√Var [T k path ] ,得出最优路径。

[2]对于t （min ）= E [T k path ]+Φ−1 （ρ）√Var [T k path ]，其中：E [T kpath ]= ∑β(a ,k )n a =1 E [T a link ] Var [T kpath ]= ∑β(a ,k )n a =1Var [T a link ]由于Var [T k path ]的根式不具有线性可加性。

不能用经典的dijkstra 算法求解。

对此我们基于双目标规划的思路[3]，分别将E [T k path ]，Var [T kpath]作为边权，运用matlab 编程（附录二），求出前十条最短路径。

表表根据上述两表数据，运用公式:t （min ）=E[T k path ]+Φ−1 （ρ） √Var [T k path ]求出并集中的前16条最优路径的预留时间（表）表为了直观进行对比，将上表用excel 制得如下柱状图：（图）（图）由图可知最优路径为V1→V3→V4→V8。

(图V7V2V3(20,10)(60,10)（图）收敛性分析：两个多项式时间算法之和还为多项式时间算法，其复杂性比列举的低。

当问题的规模趋向无穷大时，时间复杂度的数量级将表现为渐进性态。

即当路径K 趋于无穷时，该模型一定收敛。

问题三 问题分析在问题三中，要求进一步考虑各路段之间行驶时间的相关性。

我们用概率论中的协方差来表示这种耦合关系，并建立推广的NPK 模型。

模型建立与求解在问题二中我们已得出最短预留时间的数学表达式：t （min ）= [T kpath ]+Φ−1（ρ） √Var [T k path]为方便模型的建立与求解。

在此我们假设仅相邻两条路段之间具有相关性。

根据协方差的性质Var （T1+T2）=Var(T1)+Var(T2)+2cov(T1,T2);可以得出t= E [T K path ]+Φ−1(ρ)√∑δ(a ,k )Var [T a link+∑cov (a −1,a )n a =2n a =1]称t 为可靠时间。

[4] 以下图为例：图图为从A 到B 的一条路径。

可靠时间t = E [T k path ]+Φ−1 （ρ） √Var [T k path]，其中E [T kpath ]=E[T1]+E[T2]，√Var [T kpath ]=√Var [T1]+Var [T2]+cov [T1+T2]。

由于问题二中，我们没有给出任意两条路段其时间随机向量（t1，t2）的密度函数。