如何写好数学建模竞赛答案及初等模型精要.pptVIP

尚需解决的问题是如何求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=| 60 |。为实用方便，人们研究了限止开关模式个数的相应问题。 若要求r≤n，即要求通讯卫星上至多设置n种开关模式，则问题化为令r≤n，求不超过n个置换矩阵Pk及λk，使之满足： min S.t 为了使任意一对发射法与接收站之间的传送均为可能实现的，自然应要求 Pk满足 （5.1） （5.2） （右面的矩阵有n2个值为1的分量，每一Pk 恰有n个1分量）故r=n。 容易看出，（5.1）隐含着T的每一元素只能被唯一的P复盖，即T的元素在分解中是不可分割的，这当然是一个好性质，使实际操作时较为方便，但可惜的是对一般的双随机矩阵，分解很可能无解。 例3 若取 ， ， （注意：T已是双随机矩阵，行和列和均为10） 则min S.t 的解为λ1=3，λ2=4，λ3=5。 （大于10）而 但等号经常并不成立。1985年，F·Rendel证明，在给定满足（5.2）的置换矩阵P1,…,Pn后，求解问题（5.1）是NP难的，从而不可能存在多项式时间算法，除非P=NP。 现要求r≤2n 一种自然而方便的开关设置为引入两组各有n个开关模式的置换矩阵P1,…,Pn，Q1,…,Qn，满足下面的（5.3)式： 例如，当n=3时，可令： （注：这种设置方法保持了其内在的对称性，不失为一种明智的做法。） 现在，我们来分解例9.33中的双随机矩阵 ，令 = ，得方程组 求出各对角线与反对角线上的三个元素之和，并作一些简单的消去运算； 将矩阵的所有元素相加，可得下面的方程组： 注意到（5.3），易证空间 的维数为 5， 故 之一可任取，（稍加注意即可保持非负性）， 例如，令μ3=0，求得 ，故有 读者不难验证，上述方法可推广到n是奇数的一般情况。事实，由各对角线元素之和可导出n－1个方程，由各反对角线元素之和又可导出n－1个方程，加上矩阵所有元素之和导出的等式，共计可导出2 n－1个方程，并易知它们是独立的。另一方面空间 的维数恰为2 n－1，故 之一可任取，而通过方程组解得所有的 ，（只须注意保持其非负性即可） 但当n为偶数时，情况就不大相同了。让我们先来观察一下n=4的情况。 当n=4时， 易见， 具有非常特殊的结构，一般的偶数阶双随机矩 阵，即使其元素是非负整数，也无法用Pk、Qk来分解。 当 具有上述结构时，能否用Pk和Qk来分解呢？易见，由各对角线元素之和可导出： 另外，由反对角线元素之和又可导出 上述方程中只有6个是独立的，且已不可能再得出新的独立方程，（读者可自行分析之）故可选取其中2个的值，进而可解出其余。例如，若令 4= 3=0，可得 2=1， 1=0，进而可求得 1=2， 4=3， 3=3及 2=4， 已达到下界。 易见，P1 + P3 = Q2 + Q4，P2 + P4 = Q1 + Q3，故空间 的维数为6，与上面的分析是一致的。 读者可将上述讨论推广到n为一般偶数的情况，分析方法是完全类似的。 当n是偶数时，我们虽无法将一般的双随机矩阵分解为Pk 、Qk的非负组合，但上述讨论仍然是十分