system: [SYSTEM]

user: 3. 不确定随机网络带资源约束的最短路问题新型互熵优化模型及实验 3.1. 问题描述 带资源约束的最短路径问题本质上是一个最短路问题，它的目标是寻找一条从源节点1到目的节点n的满足资源约束的路径中最短的路径。 为了方便研究，在介绍模型之前，我们首先做一些假设： 1) 在本文中我们研究单源节点带资源约束的路径问题，节点1被设置为唯一的源节点，节点n被设置为唯一的终点。 2) 假定不确定随机网络 ( N ， U ， R ， W ， C ) 有向且无环，其中 N 为网络中节点的集合，一般设置为n个节点，即 N = { 1 ， 2 ， ⋯ ， n } ， U 为网络中不确定边的集合，即 U = { ( i ， j ) | ( i ， j ) 是 不 确 定 边 } ， R 网络中随机边的集合，即 R = { ( i ， j ) | ( i ， j ) 是 随 机 边 } ， W 为网络中所有的边的权重的集合，即 W = { w i j | ( i ， j ) ∈ U ∪ R } ， C 为资源成本矩阵(可设定为消耗时间，建设费用等，本文统一设置为成本)，并给定一个成本约束值d，被选择的路径还需要满足 ∑ ( i ， j ) ∈ U ∪ R x i j c i j ≤ d 这一条件的约束，记约束条件为 c ( x ) ≤ d 。 3) 假设研究的网络 ( N ， U ， R ， W ， C ) 中所有的边 ( i ， j ) 的权重为正则不确定随机变量，也就是说边的权重 w i j ，要么是正则不确定变量要么是正的随机变量，不存在其他情况。 4) 网络中存在的所有不确定随机变量相互独立，即所有的不确定变量相互独立，且所有的随机变量也相互独立。 5) 为了方便建模，在本文中使用 L = { x i j | ( i ， j ) ∈ U ∪ R ， c ( x ) ≤ d } (13) 来代表不确定随机网络 ( N ， U ， R ， W ) 中的一条从节点1到节点 n 带资源约束的路径，其中 x i j ∈ { 0 ， 1 } ， x i j = 1 代表边 ( i ， j ) 出现在此条路径当中，反之，若 x i j = 0 ，则说明 ( i ， j ) 没有出现在此路径中。若 ξ i j 代表 ( N ， U ， R ， W ) 中边的权重(通常为距离、消费、时间等的一种表现形式)，则 ( N ， U ， R ， W ， C ) 中一条路径 L 的权值为 ∑ ( i ， j ) ∈ U ∪ R x i j ξ i j (14) 另外可证明，若 L 为网络中的一条带资源约束的路径，当且仅当： { ∑ ( i ， j ) ∈ U ∪ R x i j c i j ≤ d ∑ ( i ， j ) ∈ U ∪ R x i j − ∑ ( i ， j ) ∈ U ∪ R x j i { 1 ， i = 1 0 ， 2 ≤ i ≤ n − 1 − 1 ， i = n x i j ∈ { 0 ， 1 } ， ∀ ( i ， j ) ∈ U ∪ R (15) 3.2. 新型互熵优化模型 定理3：(不确定随机带资源约束的最短路径问题的理想路径机会分布函数)假定不确定随机网络为 ( N ， U ， R ， W ， C ) ，其中的不确定边 ( i ， j ) ∈ U 的权重 w i j 对应的不确定分布为 ϒ i j ，其中的随机边 ( i ， j ) ∈ R 的权重 w i j 对应的概率分布为 Y i j ，网络连接的成本矩阵 C = ( c i j ) n × n ，限定的资源约束最大值为d，记路径 ( x | x i j ， ( i ， j ) ∈ U ∪ R ) 的连接成本 ∑ ( i ， j ) ∈ L x i j c i j 为 c ( x ) ，则从源节点到目的节点的资源约束最短路径 L = ( x | x i j ， ( i ， j ) ∈ U ∪ R ， c ( x ) ≤ d ) 的机会分布函数为 Φ ( x ) = ∫ 0 + ∞ ⋯ ∫ 0 + ∞ G ( x ; y i j ， ( i ， j ) ∈ R ， c ( x ) ≤ d ) ∏ ( i ， j ) ∈ R d Y i j ( y i j ) (16) 其中 G ( x ; y i j ， ( i ， j ) ∈ R ， c ( x ) ≤ d ) 由其逆不确定分布决定 G − 1 ( α ; y i j ， ( i ， j ) ∈ R ， c ( x ) ≤ d ) = g ( s i j ， ( i ， j ) ∈ U ∪ R ， c ( x ) ≤ d ) (17) s i j = { ϒ − 1 ( α ) ， y i j ， ( i ， j ) ∈ U ( i ， j ) ∈ R (18) 根据给定的不同的 α ， g 可以根据Flyod算法计算。 定理4：若 L = { x i j | ( i ， j ) ∈ U ∪ R ， c ( x ) ≤ d } 为不确定随机网络 ( N ， U ， R ， W ， C ) 中的一条约束路径，对应的机会分布函数为 Φ L ( x ) ，此网络的理想资源约束路径的机会分布函数为 Φ ( x ) ，则对于网络中的任意路径 L ，有 Φ L ( x ) ≤ Φ ( x ) 恒成立。 定义6：(带资源约束的新型互熵最小的最短路径)假设一个不确定随机网络 ( N ， U ， R ， W ， C ) 的理想资源约束路径对应的函数为 机会分布函数为 Φ ( x ) ，网络中存在的其他路径 L ，对应的机会分布函数为 Φ L ( x ) 。若存在某条路径 L * ，对应的函数为 ξ * ，使得其与理想路径的新型互熵 最小，即对 ∀ x ∈ R ，和所有的资源约束路径 L ，有下式成立： ∫ − ∞ + ∞ | Φ ( x ) − Φ L * ( x ) | d x ≤ ∫ − ∞ + ∞ | Φ ( x ) − Φ L ( x ) | d x (19) 则称路径 L * 为此不确定随机网络的带资源约束的新型互熵最小的最短路径。 根据定义6，如果一条路径 L * = { x i j | ( i ， j ) ∈ U ∪ R } 是一条新型互熵最小资源约束最短路径，则它应该是满足下述模型的最优解： { min { x i j ， ( i ， j ) ∈ U ∪ R } ∫ − ∞ + ∞ | Φ ( x ) − Φ L ( x | x i j ， ( i ， j ) ∈ U ∪ R ， c ( x ) ≤ d ) | d x subjectto: ∑ ( i ， j ) ∈ U ∪ R x i j c i j ≤ d ∑ j : ( i ， j ) ∈ U ∪ R x i j − ∑ j : ( i ， j ) ∈ U ∪ R x j i = { 1 ， i = 1 0 ， 2 ≤ i ≤ n − 1 − 1 ， i = n x i j = { 0 ， 1 } ， ∀ ( i ， j ) ∈ U ∪ R (20) 其中 Φ ( x ) 是理想资源约束路径的机会分布函数， Φ L ( x | ( i ， j ) ∈ U ∪ R ， c ( x ) ≤ d ) = C h { ∑ ( i ， j ) ∈ U x i j τ i j + ∑ ( i ， j ) ∈ R x i j ξ i j ≤ x } (21) 为网络中其他资源约束路径的机会分布函数。 定理5：给定一个不确定随机网络 ( N ， U ， R ， ξ ， C ) ， Φ ( x ) 是此网络的理想资源约束路径的机会分布函数， ( x | x i j ， ( i ， j ) ∈ U ∪ R ， c ( x ) ≤ d ) 为此网络中从源节点1到目的节点n的一条路径，网络中不确定边，即 ( i ， j ) ∈ U 的权重为不确定变量 τ i j ，其对应的不确定分布为 γ i j ，网络中随机边，即 ( i ， j ) ∈ R 的权重为随机变量 ξ i j ，其对应的概率分布为 Y i j 。则模型(20)可以转化为下面的模型： { min { x i j ， ( i ， j ) ∈ U ∪ R } ∫ − ∞ + ∞ | Φ ( x ) − ∫ − ∞ + ∞ γ ( x − ∑ ( i ， j ) ∈ R x i j y i j ) d Y ( y ) | d x subjectto: ∑ ( i ， j ) ∈ U ∪ R x i j c i j ≤ d ∑ j : ( i ， j ) ∈ U ∪ R x i j − ∑ j : ( i ， j ) ∈ U ∪ ​ R x j i = { 1 ， i = 1 0 ， 2 ≤ i ≤ n − 1 − 1 ， i = n x i j = { 0 ， 1 } ， ∀ ( i ， j ) ∈ U ∪ R (22) 其中 Φ ( x ) = ∫ 0 + ∞ ⋯ ∫ 0 + ∞ G ( x ; y i j ， ( i ， j ) ∈ R ， c ( x ) ≤ d ) ∏ ( i ， j ) ∈ R d Y i j ( y i j ) (23) Y ( y ) = ∫ y 1 + y 2 + ⋯ + y m ≤ y d Y 1 ( y 1 ) d Y 2 ( y 2 ) ⋯ d Y m ( y m ) (24) ϒ ( r ) = sup ∑ ( i ， j ) ∈ U x i j r i j = r ( min ( i ， j ) ∈ U ϒ i j ( x i j r i j ) ) (25) 证明：根据定理3可得，网络的理想路径机会分布函数为 Φ ( x ) = ∫ 0 + ∞ ⋯ ∫ 0 + ∞ G ( x ; y i j ， ( i ， j ) ∈ R ， c ( x ) ≤ d ) ∏ ( i ， j ) ∈ R d Y i j ( y i j ) (26) 假设 L * = ( x | x i j ， ( i ， j ) ∈ U ∪ R ， c ( x ) ≤ d ) 为此网络中从源节点1到目的节点n的一条资源约束路径，则 Φ L ( x ) = C h { ∑ ( i ， j ) ∈ U x i j τ i j + ∑ ( i ， j ) ∈ R x i j ξ i j ≤ x } = ∫ R M { ∑ ( i ， j ) ∈ U x i j τ i j + y ≤ x } Y ( y ) (27) 其中 Y ( y ) 是 ∑ ( i ， j ) ∈ R x i j ξ i j 的概率分布函数。 由于所有的不确定随机变量 ξ i j 彼此相互独立，故 Y ( y ) = ∫ y 1 + y 2 + ⋯ + y m ≤ y d Y 1 ( y 1 ) d Y 2 ( y 2 ) ⋯ d Y m ( y m ) (28) 又由于所有的不确定随机变量相互独立，故根据例1有 M { ∑ ( i ， j ) ∈ U x i j τ i j ≤ x − y } = sup ∑ ( i ， j ) ∈ U x i j r i j = x − y ( min ( i ， j ) ∈ U ϒ i j ( x i j r i j ) ) = ϒ ( x − ∑ ( i ， j ) ∈ R x i j y i j ) . (29) 结合式子(27) (28) (29)带入目标函数方程，得证。 3.3. 新型互熵优化模型算法及数值实验 由于在模型求解过程中都涉及到理想路径机会分布，但对于此类问题而言，直接给出机会分布函数的显式表达式是很困难的，于是在本文中我们设计了一个线性插值的办法用来拟合机会分布函数，步骤如下： 算法1：(用于求解理想带资源约束路径的机会分布)。 步骤1：对于不确定随机网络中的每个随机变量，我们根据随机变量的概率分布函数，使用蒙特卡洛模拟的方法生成1000个数，随机变量的取值便为蒙特卡洛模拟生成的依数值大小排序后的随机数值。 步骤2：设置矩阵 T = [ 0 ， 0.01 ， ⋯ ， 0.99 ， 1 ] ， α 依次从T中取值，根据上述过程生成的数值计算 G − 1 ( α ) ，得到不确定分布 G ( x ) 的离散形式。 步骤3：对步骤2生成的离散形式的 G ( x ) ，进行线性插值，并带入(4)式计算理想路径的机会分布 Φ ( x ) 。 基于3.2节提出模型——不确定随机网络下的带资源约束的最短路径新型互熵优化模型(22)，我们设计了下面算法用来求解。 算法2：(适用于不确定随机网络的带资源约束的最短路径新型互熵优化模型(22)的求解)。 步骤1：结合定理1并利用上述提出的算法1计算理想路径的机会分布函数 Φ ( x ) 。 步骤2：利用智能算法( A * 算法、广度优先算法、蚁群算法、遗传算法等均可)搜寻网络中的所有路径并储存。 步骤3：根据成本矩阵计算每条路径的成本，去掉不符合条件约束的路径。 步骤4：计算步骤3中储存所有路径与理想路径之间的新型互熵值。 步骤5：比较上一步得到的新型互熵值，选择新型互熵值最小的路径，即为所求。 下面给出一个数值实验，来验证我们提出的算法与模型的有效性。 例3给定一个具有8个节点12条边的不确定随机网络图 ( N ， U ， R ， W ， C ) 如图2所示，假设不确定权重为 τ 23 ， τ 25 ， τ 38 ， τ 45 ， τ 46 ， τ 58 ， τ 67 ， τ 78 它们为相互独立的正则不确定变量，且对应的正则不确定分布为 ϒ 23 ， ϒ 25 ， ϒ 38 ， ϒ 45 ， ϒ 46 ， ϒ 58 ， ϒ 67 ， ϒ 78 ，随机边的权重为 ξ 12 ， ξ 14 ， ξ 16 ，它们为相互独立的不确定变量，对应的概率分布函数分别为 Κ 12 ， Κ 14 ， Κ 16 ，假设此不确定随机网络图 ( N ， U ， R ， W ， C ) 的权重与成本如表1所示，其中 U ( a ， b ) 代表线形不确定分布， Z ( a ， b ， c ) 代表Z字形不确定分布。给定成本约束值d = 200。 图2. 8节点12条边不确定随机网络 Table 1 The edge weights and connection costs of the uncertain random network in Figure ( i ， j ) ( i ， j ) 权值 边 ( i ， j ) 成本 ( i ， j ) ( i ， j ) 权值 边 ( i ， j ) 成本 ( i ， j ) ( i ， j ) 权值 边 ( i ， j ) 成本 ( 2 ， 3 ) L ( 50 ， 74 ) 60 ( 4 ， 5 ) Z ( 29 ， 37 ， 49 ) 62 ( 1 ， 2 ) U ( 45 ， 59 ) 70 ( 2 ， 5 ) L ( 27 ， 43 ) 77 ( 4 ， 7 ) Z ( 33 ， 42 ， 45 ) 49 ( 1 ， 4 ) U ( 34 ， 52 ) 90 ( 3 ， 8 ) L ( 31 ， 45 ) 39 ( 5 ， 8 ) Z ( 6 ， 17 ， 25 ) 16 ( 1 ， 6 ) U ( 60 ， 90 ) 75 ( 4 ， 6 ) L ( 52 ， 79 ) 48 ( 7 ， 8 ) Z ( 15 ， 37 ， 42 ) 30 ( 6 ， 7 ) L ( 48 ， 63 ) 70 表1. 图2的不确定随机网络的边权值与连接成本表 根据定理3可得 Φ ( x ) = ∫ 0 + ∞ ∫ 0 + ∞ ∫ 0 + ∞ G ( x ; y 12 ， y 14 ， y 16 ， c ( x ) ≤ 200 ) d Y 12 ( y 12 ) d Y 14 ( y 14 ) d Y 16 ( y 16 ) (30) 其中 G ( x ; y 12 ， y 14 ， y 16 ， c ( x ) ≤ 200 ) 可以根据其逆不确定分布计算。利用深度优先算法搜索网络发现一共有6条路径，但通过计算成本发现其中有一条路径 L ( 1 → 4 → 6 → 7 → 8 ) 的成本值为238 > 200，故在第三步中完成后，提前删除此路径，将剩余的五条路径进行算法的第四步运算计算新型互熵值，计算结果列在表2中。 通过比较新型互熵的数值结果，我们发现路径3 ( 1 → 4 → 5 → 8 ) 新型互熵值最小，而且通过图3的机会分布图也可以验证我们的计算结果的有效性——路径3是最接近理想约束路径的，即 ( 1 → 4 → 5 → 8 ) 为此不确定随机网络下的带资源约束的新型互熵最小的最短路径。 Table 2 The calculation results about the newly cross entropy mode 路径 新型互熵计算的数值结果 1 → 2 → 3 → 8 35.9699 1 → 2 → 5 → 8 44.4151 1 → 4 → 5 → 8 4.95090 1 → 4 → 7 → 8 83.3799 1 → 6 → 7 → 8 47.6151 表2. 新型互熵计算结果 图3. 理想资源约束路径与不确定随机网络中的其他约束路径的机会分布