传输理论

传输理论（tansport theory、transportation theory），又称为运输理论，是数学、经济学等学科中研究最优运输和资源配置的理论。该问题最早由法国数学家加斯帕尔·蒙日于1781年提出。
1920年代，A·N·托尔斯泰是最早运用数学方法研究传输问题的学者之一。1930年，他在苏联国家交通部编纂的《运输规划》第一卷中发表了题为《寻找太空货物运输的最小千公里方法》的论文。
第二次世界大战期间，苏联数学家、经济学家列昂尼德·坎托罗维奇在该领域取得了重要进展。因此，这一问题有时也被称为蒙日-坎托罗维奇运输问题（Monge–Kantorovich transportation problem）。 该问题的线性规划形式也被称为-库普曼斯运输问题。
背景.
矿山与工厂.
假设有formula_1个开采铁矿石的矿山以及formula_2个工厂使用这些矿山生产的铁矿石，且这些矿山和工厂构成欧几里得平面formula_3中两个不相交的子集formula_4和formula_5。同时假设存在一个成本函数formula_6， 其中formula_7表示将一批矿石从formula_8运送到formula_9的成本。为简化起见，此处忽略运输所需的时间。我们还假定每个矿山只能供应一家工厂（不能拆分运输），并且每家工厂需要恰好一批货物才能运营（工厂不能半负荷或双倍负荷运转）。基于上述条件，一个传输计划可以看作是一个双射formula_10 。换句话说，每个矿井formula_11仅供应一个目标工厂formula_12，而每个工厂也只由一个矿山供货。我们希望找到最优传输计划formula_13，使得总成本
在所有formula_4到formula_5的传输计划中是最小的。该问题是传输问题的一个特例，可看成一个任务分配问题。更具体地说，它等价于在二分图中寻找最小权重匹配。
书籍移动：成本函数的重要性.
下面这个简单的例子说明了成本函数在确定最优传输计划中的重要性。假设我们有formula_2本宽度相等的书摆放在书架上（可以看成具像化的实数线），形成连续一排书。我们希望将它们重新排列，在保持其连续性的同时将整体向右移动一本书的宽度。针对这个问题，有两个显而易见的最优传输候选方案：
如果成本函数与欧几里得距离成正比（即 formula_20，其中formula_21 ），那么这两种候选方案都是最优的。而如果我们选择与欧几里得距离的平方成比例的严格凸成本函数（即 formula_22，其中formula_21 ），则“许多小步”的方案则是唯一的最优解。
需要注意的是，上述成本函数仅考虑书籍本身移动的水平距离，而没有考虑拿起每本书并将其移动到位的设备所行进的水平距离。如果考虑后者，那么在两种传输计划中，第二种方案对于欧几里得距离始终是最优的，而第一种方案则对于平方欧几里得距离是最优的（至少有三本书的情况下）。
希区柯克问题.
以下传输问题的表述由弗兰克·劳伦·希区柯克提出：
佳林·库普曼斯也为运输经济学与资源分配问题的表述作出了贡献。
问题的抽象表述.
蒙日和坎托罗维奇形式.
由于黎曼几何和测度论的发展，在现代或更加技术性的文献中传输问题的表述有所不同。不过，上述矿山与工厂的简单例子还是可以作为考虑抽象形式时一个有用的参照。此时我们可以考虑并非所有矿山和工厂都营业的情况，并允许一个矿山向多家工厂供货，而一家工厂也可以从多个矿山接受矿石。
假设formula_35和formula_36为两个可分度量空间，使得formula_35（或者formula_36 ) 上的每个概率测度都是拉东测度，并假设formula_39是一个博雷尔可测函数。给定formula_35上的概率测度formula_41和formula_36上的概率测度formula_43，最优传输问题的蒙日形式是指寻找一个传输映射formula_44，使得下式中的下确界成立：
其中formula_46由formula_13推进formula_41的向前推进算子。如果一个映射formula_13达到这个下确界，则该映射被称为“最优传输映射”。
最优传输问题的蒙日形式有其局限性，因为有时可能不存在满足formula_50的映射formula_13。例如，当formula_41是狄拉克测度而formula_43不是时，就会出现这种情况。
此时可以通过采用最优传输问题的坎托罗维奇形式来克服这一局限，即寻找formula_54上的一个概率测度formula_55，使得下式中的下确界成立：
其中formula_57表示formula_54上所有概率测度的集合并满足边缘分布formula_41和formula_43。可以证明，当成本函数formula_61是下半连续并且formula_62是紧测度集合（拉东空间formula_35和formula_36蕴含该条件）时，该问题总是存在最小值（另见沃瑟斯坦度量）。、史蒂文·哈克尔（Steven Haker）与提出了蒙日–坎托罗维奇问题解的梯度下降表述。
对偶形式.
坎托罗维奇问题的最小值等于
其中上确界遍历所有有界且连续、满足
的函数对。
经济学解释.
将以上表述中的符号翻转更利于从经济学角度解释这一问题。假设formula_67表示工人特征的向量，formula_68表示企业特征的向量，formula_69则表示工人formula_8与企业formula_9配对所创造的经济产出。令formula_72、formula_73，蒙日–坎托罗维奇问题可以重新表述为：
它的对偶形式为：
其中下确界遍历所有有界且连续的函数formula_76和formula_77。如果对偶问题有解，则有：
可以将formula_79解释为formula_8类型工人的均衡工资 ，并将formula_81解释为formula_9类型企业的均衡利润。
问题求解.
一维连续情形.
对于formula_83， 假设formula_84表示formula_85上所有formula_86阶矩有限的概率测度的集合。设formula_87且formula_88，其中formula_89是一个凸函数。
拉切夫（Rachev）与吕申多夫（Rüschendorf）于1998年给出了对此的证明。
离散情形与线性规划.
在边缘分布formula_41和formula_43是离散的情形下，令formula_98和formula_99分别是分配给formula_100和formula_101的概率质量，而formula_102则是formula_103的分配概率。原始坎托罗维奇问题中的目标函数为
并满足约束条件
为了将这一问题作为线性规划问题处理，我们需要将矩阵formula_102向量化，可以通过堆叠其行或列来完成，我们用formula_108来表示这一操作。在的情况下，上述约束条件可改写为
其中formula_111是克罗内克积，formula_112是一个大小为formula_113、所有元素为1的矩阵，而formula_114是大小为formula_2的单位矩阵。设formula_116，该问题的线性规划形式为
这一问题可以很容易地通过大规模线性规划求解器计算。
半离散情形.
在半离散情况下，令formula_118，且formula_41是formula_120上的连续分布，formula_121则是分配概率质量formula_122到formula_123的离散分布。此时，坎托罗维奇问题的原始形式为：
其中formula_125满足formula_126和formula_127 。
而其对偶形式则为
还可以写为：
这是一个有限维凸优化问题，可以通过梯度下降等方法求解。
当formula_130时，可以证明分配给特定formula_131的formula_100集合是一个凸多面体，而得到的配置称为。
二次正态情形.
假设一个特殊情形formula_133、formula_134 且formula_135，其中formula_136是可逆矩阵。此时有
加利雄（Galichon）于2016年证明了该情况下的解。
可分希尔伯特空间.
令formula_35是一个可分希尔伯特空间，定义formula_141为formula_35上所有formula_86阶矩有限的概率测度的集合，formula_144则表示其中高斯正则的测度集合，即如果formula_145是formula_35上任何严格正的且满足formula_147 ，则formula_148也成立。
假设formula_149、formula_150，并且formula_151，其中formula_152 。则坎托罗维奇问题存在一个唯一解formula_153 ，并且该解对应一个最优传输映射：即存在一个博雷尔映射formula_154，使得
此外，如果formula_43具有有界支撑，那么对于formula_41-几乎所有的formula_158，存在局部利普希茨、 formula_61-凹和最大坎托罗维奇势formula_160，使得 
其中formula_162表示formula_160的加托导数。
熵正则化.
考虑上述离散问题的一个变体：在原始问题的目标函数中添加一个熵正则化项
相应的对偶问题为
相较于不含正则化项的问题，原先对偶问题中的硬约束（ formula_166 ）被替换为了软约束，即惩罚项formula_167。对偶问题的最优条件可以表示为
令formula_136为formula_171的矩阵，其中元素formula_172。此时对偶问题的求解等价于寻找两个对角正矩阵formula_173和formula_174，它们的大小分别为formula_175和formula_176 ，使得formula_177、formula_178。矩阵formula_173和formula_174的存在性是辛克宏定理的推广，可以使用辛克宏-诺普算法进行求解。该算法通过迭代求解中的formula_181和中的formula_182实现。因此，辛克宏-诺普算法相当于对偶正则问题的坐标下降法。
应用.
蒙日-坎托罗维奇运输问题已广泛运用于许多领域，例如：