离散优化模型乐及算法设计.pptVIP

第九章 离散优化模型及算法设计 浙江大学数学建模基地  第九章结束 The end （设计方法1） 注意到Pk每行（或列）元素之和均为1，故不管如何指派开关的使用时间（即不论如何取λk），矩阵 均具有某些特殊的性质，例如其行和（及列和）均为同一常数。这样的矩阵构成一个线性空间（参见逻辑模型第一节 Dürer魔方），为减少开关模式的种类，可取此空间的一组基底作为开关模式。在使用这种开关模式时，无论T的元素tij怎么取，通讯卫星对每一发（收）点的开通时间总和是恒定的。在这种开关模式下，可按如下方式指派各开关模式的使用时间： 步1 先将T改变为 ， 满足： （1） ≥T （2）记 ， 步2 用Pk表示 ，即将 分解为 （r为空间的维数） 将T化为 的方法一般有无穷多种，如可如下化法： 令 事实上， ，（即通讯卫星传送总时间的下界）。 令 其中 用这种方法化例9.39中的T，得到 的任一行（或列）中元素之和均为7。 定义9.13 称行和、列和均相等的矩阵为双随机矩阵（Doubly stochastic matrix） 定理9.24 （Birkhoff定理，1944）任一n阶双随机矩阵均可写成至多 （n－1）2+1个置换矩阵的非线性组合。 的分解可如下进行： 步1 选取由Pij0可推出 0的置换矩阵P 步2 确定 步3 取 ，用 － 代替 步4 若 =0，停；否则，返回步1。 例9.33 为方便起见，我们来分解一个元素均为非负整数的3阶双随机矩阵， （由Birkhoff定理，r≤5） 解：取 ，λ=min {1, 3, 3 } =1 分解成 ，再取 因min{ 5, 5, 3} = 3，又有 ，取 于是又有 易得分解结果为： 尚需解决的问题是如何求P，使得Pij0必有 。读者不难发现，此问题可以通过求解一个两分图上的最大流（或最大匹配）来实现，计算量为O(n4)，是多项式时间可解的。具体方法为：作一两分图，若 ，则作边（i, j），令边容量为1，这样，可作出P的充要条件是该最大流问题的最大流量为n。对例9.33，n=3。由于所有 ，先取 ， －P1为 相应的两分图为： 于是又可求得 ，相应的两分图为： 又可得 ，…，如此下去，直到作不出P为至， 由于 的特殊性质及Birkhoff定理，上述分解必能在不超过r= (n－1)2 + 1步内终止。 上述开关设计方法要求在通讯卫星上设置(n－1)2 + 1种不同的开关模式（即Pk），当n稍大时，(n－1)2 + 1仍显得太大而使得使用时不便。例如，当n=41时，(n－1)2 + 1=1601 。为实用方便，人们研究了限止开关模式个数的相应问题。 若要求r≤n，即要求通讯卫星上至多设置n种开关模式，则问题化为令r≤n，求不超过n个置换矩阵Pk及λk，使之满足： min S.t 为了使任意一对发射站与接收站之间的传送均为可能实现的，自然应要求 Pk满足 （9.5） （9.6） （右面的矩阵有n2个值为1的分量，每一Pk 恰有n个1分量）故r=n。 容易看出，（9.5）隐含着T的每一元素只能被唯一的P复盖，即T的元素在分解中是不可分割的，这当然是一个好性质，使实际操作时较为方便，但可惜的是对一般的双随机矩阵，分解很可能无解。 例9.34 若取 ， ， （注意：T已是双随机矩阵，行和列和均为10） 则min S.t 的解为λ1=3，λ2=4，λ3=5。 （大于10）而 但等号经常并不成立。1985年，F·Rendel证明，在给定满足（9.6）的置换矩阵P1,…,Pn后，求解问题（9.5）是NP难的，从而不可能存在多项式时间算法，除非P=NP。 现要求r≤2n 一种自然而方便的开关设置为引入两组各有n个开关模式的置换矩阵P1,…,Pn，Q1,…,Qn，满足 例如，当n=3时，可令： （注：这种设置方法保持了其内在的对称性，不失为一种明智的做法。） 现在，我们来分解例9.33中的双随机矩阵