Аннотация:
Рассматривается оптимальная транспортировка мер на метрических и топологических пространствах в случае, когда функция стоимости и маргинальные распределения зависят от параметра со значениями в метрическом пространстве. Расстояние Хаусдорфа между множествами вероятностных мер с заданными проекциями оценивается через расстояния между самими проекциями. Эта оценка используется для доказательства непрерывности стоимости оптимальной транспортировки относительно параметра в случае непрерывной зависимости функции стоимости и маргинальных распределений от этого параметра. Установлено существование приближенных оптимальных планов, непрерывных относительно параметра. Показано, что оптимальный план непрерывен по параметру в случае единственности. Однако построены примеры, когда не существует непрерывного выбора оптимальных планов. Другое применение оценки для расстояния Хаусдорфа связано с дискретными приближениями транспортных задач. Наконец, доказан общий результат о сходимости оптимальных отображений Монжа.
Библиография: 46 названий.
Ключевые слова:
задача Канторовича, задача Монжа, расстояние Хаусдорфа, каплинг, слабая сходимость, непрерывность по параметру.
Напомним, что задача Канторовича оптимальной транспортировки имеет дело с тройкой (μ,ν,h), где μ и ν – борелевские вероятностные меры на топологических пространствах X и Y соответственно и h⩾0 – борелевская функция на X×Y. Задача заключается в минимизации интеграла
∫hdσ
по всем мерам σ из множества Π(μ,ν), состоящего из борелевских вероятностных мер на X×Y с проекциями μ и ν на сомножители, т.е. σ(A×Y)=μ(A) и σ(X×B)=ν(B) для всех борелевских множеств A⊂X и B⊂Y. Меры μ и ν называются маргинальными распределениями или маргиналами, а функция h называется функцией стоимости. Меры из Π(μ,ν) называют каплингами мер μ и ν или транспортными планами Канторовича. В общем случае имеется лишь инфимум Kh(μ,ν) (который может быть бесконечным), но если функция h непрерывна (или полунепрерывна снизу) и ограниченна, а меры μ и ν радоновы, то достигается минимум, а меры, на которых он достигается, называются оптимальными мерами или оптимальными планами Канторовича. Более того, ограниченность h можно заменить предположением, что существует мера из Π(μ,ν), относительно которой функция h интегрируема. Задача имеет смысл также и в чисто теоретико-множественной постановке, но здесь мы рассматриваем топологический случай, более того, рассматриваемые пространства будут вполне регулярными, в некоторых результатах метрическими. Аналогично вводится мультимаргинальная задача Канторовича с маргиналами μ1,…,μn на пространствах X1,…,Xn. Общую информацию о задачах Монжа и Канторовича можно найти в работах [3], [4], [12], [17], [25], [36], [39], [42], [43].
Мы изучаем оптимальную транспортировку мер на метрических и топологических пространствах в случае, когда функция стоимости ht и маргинальные распределения μt и νt зависят от параметра t со значениями в метрическом пространстве. Здесь естественно возникают вопросы о непрерывности по t оптимальной стоимости Kht(μt,νt), а также о возможности выбора оптимального плана в Π(μt,νt), непрерывного по параметру. Кроме того, множество всех транспортных планов Π(μt,νt) также зависит от параметра, так что можно задать вопрос о его непрерывности, когда пространство множеств мер наделяется метрикой Хаусдорфа, порожденной некоторой метрикой на пространстве мер. Наш первый основной результат дает простую оценку для расстояния Хаусдорфа между множествами каплингов через расстояния Хаусдорфа между их маргиналами. Качественный результат такого рода был установлен в [9], позже другой подход был развит в [40]; в недавней работе [26] более короткое обоснование непрерывности соответствия (μ,ν)↦Π(μ,ν) было предложено с помощью метрики Прохорова (см. [26; теорема 1]). Мы показываем, что расстояние Канторовича позволяет получить более простую и информативную оценку.
Задачи Канторовича, зависящие от параметра, исследовались в ряде работ, см. [13]–[15], [18], [23], [30], [32], [43], [46], где рассматривались вопросы измеримости. Свойства непрерывности также представляют большой интерес и имеют фундаментальное значение; в частности, они могут быть полезны в изучении дифференциальных уравнений и включений в пространствах мер (см. [21]), в регуляризации оптимальных транспортировок (см. [22] и [31]), в построении приближений дискретными транспортными задачами (см. [6], результат такого рода дан ниже), а также в других приложениях (см., например, [1], [10], [41], [45]). Как уже было упомянуто, непрерывная зависимость от маргиналов рассматривалась в [9], [40] и [26]. Близкие задачи возникают для нелинейных функционалов стоимости, так называемых слабых транспортных стоимостей, рассмотренных в [27], [2], [7], [8].
Основные результаты настоящей работы состоят в следующем.
1. Расстояние Хаусдорфа между множествами вероятностных мер Π(μ1,ν1) и Π(μ2,ν2) оценено через расстояния между μ1 и μ2 и между ν1 и ν2 (см. теорему 1). Аналогичный результат верен для p-метрики Канторовича Wp. Эта оценка применяется для построения дискретных приближений задач транспортировки. Фактически установлена более общая оценка для псевдометрик Канторовича и Хаусдорфа в случае общих вполне регулярных пространств.
2. Стоимость оптимальной транспортировки непрерывна по параметру в случае непрерывной зависимости функции стоимости и маргинальных распределений от этого параметра (см. теорему 2 и теорему 5 в случае нелинейного функционала стоимости).
3. Оптимальный план непрерывен по параметру в случае единственности (см. следствие 2).
4. Существуют приближенные оптимальные планы, непрерывные по параметру (см. теорему 3).
5. Построены примеры, когда нет непрерывного выбора оптимальных планов на [0,1]2 с функцией стоимости, непрерывной по параметру, и маргинальными распределениями, равными мере Лебега.
Наконец, доказан общий результат о сходимости оптимальных отображений Монжа в духе известных результатов из [5] и [24]: согласно следствию 4 оптимальные отображения Монжа, переводящие μn в νn, сходятся по мере μ0, если μn→μ0 по вариации, νn→ν0 слабо, функции стоимости hn сходятся к функции стоимости h0 равномерно на компактах, а соответствующие оптимальные отображения Монжа и оптимальные планы Канторовича единственны. Часть этих результатов была анонсирована в заметке [19].
§ 2. Обозначения и терминология
Напомним, что неотрицательная радоновская мера на топологическом пространстве X есть ограниченная борелевская мера μ⩾0 такая, что для каждого борелевского множества B и каждого ε>0 найдется компакт K⊂B, для которого μ(B∖K)<ε (см. [11]). Если X – полное сепарабельное метрическое пространство, то все борелевские меры радоновы. Ограниченная знакопеременная борелевская мера μ называется радоновской, если ее полная вариация |μ| радонова.
Пространство Mr(X) знакопеременных ограниченных радоновских мер на X может быть наделено слабой топологией, порожденной полунормами
μ↦|∫fdμ|,
где f – ограниченная непрерывная функция. Подмножество Pr(X) вероятностных радоновских мер наделяется индуцированной топологией.
Для двух заданных мер μ∈Pr(X) и ν∈Pr(Y) множество Π(μ,ν) состоит из всех мер σ∈Pr(X×Y) с проекциями μ и ν на сомножители.
Если h⩾0 – борелевская функция на X×Y (называемая функцией стоимости), то для мер μ∈Pr(X) и ν∈Pr(Y)стоимость по Канторовичу определяется равенством
Kh(μ,ν)=infσ∈Π(μ,ν)∫X×Yhdσ,
где возможно также значение +∞.
Множество M неотрицательных радоновских мер на пространстве X называется равномерно плотным, если для всякого ε>0 существует такой компакт K⊂X, что μ(X∖K)<ε для всех μ∈M.
Пусть (X,dX) и (Y,dY) – метрические пространства. Пространство X×Y наделяется метрикой
d((x1,y1),(x2,y2))=dX(x1,x2)+dY(y1,y2).
Слабая топология на пространствах радоновских вероятностных мер Pr(X), Pr(Y), Pr(X×Y) метризуема соответствующими метриками Канторовича–Рубинштейна dKR (называемыми также метриками Форте–Мурье, см. [12]), определенными формулой
dKR(μ,ν)=sup{∫fd(μ−ν):f∈Lip1,|f|⩽1},
где Lip1 – множество 1-липшицевых функций. Если X полно, то (Pr(X),dKR) также полно, а если X – польское пространство, то Pr(X) также оказывается польским.
Подмножества P1r(X), P1r(Y), P1r(X×Y), состоящие из мер, относительно которых все функции вида x↦d(x,x0) интегрируемы, наделяются метрикой Канторовича
dK(μ,ν)=sup{∫fd(μ−ν):f∈Lip1}.
Для этой метрики используется также термин “Wasserstein distance”, но мы не будем использовать это исторически некорректное название.
Поскольку мы рассматриваем вероятностные меры, в формуле для dK супремум может вычисляться по функциям f с дополнительным условием f(x0)=0 для фиксированной точки x0. Значит, для пространства, содержащегося в шаре радиуса 1, выполнено равенство dK=dKR.
Заметим, что неограниченное метрическое пространство (X,d) может быть наделено ограниченной метрикой ˜d=min(d,1), порождающей исходную топологию. Для новой метрики имеем ˜dK=˜dKR. Более того,
2−1˜dK⩽dKR⩽2˜dK.
В самом деле, если |f|⩽1 и f – 1-липшицева функция по исходной метрике, то относительно новой метрики f является 2-липшицевой. С другой стороны, всякая функция f, которая 1-липшицева по новой метрике и равна нулю в точке x0, удовлетворяет оценке |f|⩽1 и является 2-липшицевой по исходной метрике, поскольку ˜d(x,y)=d(x,y) при d(x,y)⩽1, а при d(x,y)>1 нужное неравенство вытекает из оценки |f|⩽1.
Стоит отметить, что если X=Y и в качестве функции стоимости взято расстояние, то выполнено равенство Kd(μ,ν)=dK(μ,ν) на P1r(X), называемое формулой двойственности Канторовича.
Аналогично для всякого p∈[1,+∞) подпространство Ppr(X) в Pr(X), состоящее из мер, относительно которых функция x↦d(x,x0)p интегрируема для некоторой (тогда и для всякой) фиксированной точки x0, может быть наделено p-метрикой Канторовича
Wp(μ,ν)=Kdp(μ,ν)1/p.
Напомним, что расстояние Хаусдорфа между ограниченными замкнутыми подмножествами A и B метрического пространства (M,d) определяется формулой
H(A,B)=max{supx∈Ad(x,B),supy∈Bd(y,A)}.
Это расстояние будет рассматриваться для подмножеств пространства вероятностных мер Pr(X×Y) с метрикой Канторовича–Рубинштейна dKR (порожденной метрикой на X×Y, введенной выше) или его подпространства P1r(X×Y) с метрикой Канторовича dK, что дает соответствующие расстояния Хаусдорфа HKR и HK. Благодаря оценке (2.1) можно иметь дело с последним случаем.
Когда Ppr(X×Y) наделено метрикой Wp, получаем соответствующее расстояние Хаусдорфа Hp на пространстве замкнутых подмножеств Ppr(X×Y).
Аналогичные конструкции вводятся в более общем случае вполне регулярных пространств X и Y. Топологии таких пространств можно задать семействами псевдометрик ΨX и ΨY (напомним, что псевдометрика d отличается от метрики тем, что разрешается равенство d(x,y)=0 при x≠y). Тогда топология произведения X×Y порождается псевдометриками
Для псевдометрики d∈ΨX указанным выше способом определяются псевдометрики Канторовича и Канторовича–Рубинштейна dK,d и dKR,d на пространствах PΨr(X) и Pr(X), где первое состоит из мер, относительно которых функции x↦d(x,x0) интегрируемы для всех псевдометрик из ΨX. Легко проверить, что слабая топология на Pr(X) порождается семейством таких псевдометрик dKR,d (см. работу [1], где изучаются топологии Канторовича и Канторовича–Рубинштейна на пространствах мер на вполне регулярных пространствах).
Двойственность Канторовича Kd(μ,ν)=dK(μ,ν) имеет место также для псевдометрик (см. [37] или [17; теорема 1.3.1]).
Задача Монжа для той же самой тройки (μ,ν,h) состоит в нахождении борелевского отображения T:X→Y, переводящего μ в ν, т.е. удовлетворяющего равенству ν=μ∘T−1, где (μ∘T−1)(B)=μ(T−1(B)) для всех борелевских множеств B⊂Y, для которого минимален интеграл
∫h(x,T(x))μ(dx).
Опять в общем случае есть только инфимум h(μ,ν) этого интеграла (возможно, бесконечный), но во многих интересных случаях существуют оптимальные отображения Монжа. В любом случае Kh(μ,ν)⩽Mh(μ,ν), но для безатомических сепарабельных радоновских мер и непрерывной функции стоимости верно равенство Kh(μ,ν)=Mh(μ,ν); см. [16], [35]. Из этого равенства вытекает, что если есть единственное решение T задачи Монжа, то образ μ при отображении x↦(x,T(x)) является оптимальным планом Канторовича.
Нам понадобится ниже так называемая лемма о склейке, см., например, лемму 1.1.6 в [17] или лемму 3.3.1 в [12] (в последней формулировка дана для метрических пространств, но доказательство фактически приведено для радоновских мер на вполне регулярных пространствах). Пусть X1, X2, X3 – вполне регулярные пространства и μ1,2, μ2,3 – радоновские вероятностные меры на X1×X2 и X2×X3 соответственно, причем их проекции на X2 совпадают. Тогда существует такая радоновская вероятностная мера μ на X1×X2×X3, что ее проекция на X1×X2 есть μ1,2 и проекция на X2×X3 есть μ2,3.
§ 3. Основные результаты
Мы начнем со следующей общей оценки для мер на вполне регулярных пространствах. Для функций α1,…,αn на пространствах X21,…,X2n (где X2i=Xi×Xi) положим
Теорема 1. Пусть μ1,μ2∈Pr(X), ν∈Pr(Y), и пусть α и β – непрерывные неотрицательные функции на X2 и Y2 соответственно, где β(y,y)=0. Тогда для всякой меры σ1∈Π(μ1,ν) найдется такая мера σ2∈Π(μ2,ν), что
Kα⊕β(σ1,σ2)⩽Kα(μ1,μ2).
Если ν1,ν2∈Pr(Y) и α,β – псевдометрики, то для всякой меры σ1∈Π(μ1,ν1) найдется такая мера σ2∈Π(μ2,ν2), что
dK,α⊕β(σ1,σ2)⩽dK,α(μ1,μ2)+dK,β(ν1,ν2).
Значит, для соответствующих псевдометрик Канторовича и Хаусдорфа имеем
Аналогичное утверждение верно для n маргиналов: если μi,νi∈Pr(Xi), i=1,…,n, αi – непрерывные псевдометрики на Xi, σ∈Π(μ1,…,μn), то существует такая мера π∈Π(ν1,…,νn), что
dK,α1⊕⋯⊕αn(π,σ)⩽dK,α1(μ1,ν1)+⋯+dK,αn(μn,νn).
Доказательство. Возьмем такую меру η∈Π(μ1,μ2), что
∫X×Xα(x1,x2)η(dx1dx2)=Kα(μ1,μ2).
Согласно упомянутой выше лемме о склейке, применяемой к X1=X2=X, существует такая мера λ∈Pr(X1×X2×Y), что ее проекция на X1×Y есть σ1 и ее проекция на X1×X2 есть η. В качестве σ2 возьмем проекцию меры λ на произведение X2×Y. Проекции меры σ2 на X и Y равны μ2 и ν соответственно. В самом деле, проекция на X получается проектированием меры λ сначала на X2×Y, затем на X2, что совпадает с композицией операторов проектирования на X1×X2 и X2, но эта композиция переводит λ в μ2, поскольку λ сначала отображается в η и затем в μ2. Проекция меры σ2 на Y получается проектированием на X1×Y и затем на Y, т.е. равна проекции меры σ1 на Y.
Для доказательства желаемой оценки рассмотрим меру ζ, равную образу λ при отображении
Замечание 1. Если вместо непрерывности α и β потребовать только их борелевскую измеримость (или универсальную измеримость), то рассуждение, использованное выше, покажет, что для всякого ε>0 найдется мера σε2, для которой
Kα⊕β(σ1,σε2)⩽Kα(μ1,μ2)+ε.
Для этого возьмем такую меру ηε∈Π(μ1,μ2), что
∫X×Xα(x1,x2)ηε(dx1dx2)⩽Kα(μ1,μ2)+ε.
Тогда последнее равенство в доказательстве с величиной Kα(μ1,μ2) заменится на неравенство с величиной Kα(μ1,μ2)+ε в правой части.
Непосредственная модификация доказательства дает следующее утверждение для p-метрики Канторовича Wp.
Следствие 1. Предположим, что X – метрическое пространство, p∈[1,+∞) и μ1,μ2,ν1,ν2∈Ppr(X). Тогда для всякой меры σ1∈Π(μ1,ν1) существует такая мера σ2∈Π(μ2,ν2), что
Wp(σ1,σ2)⩽Wp(μ1,μ2)+Wp(ν1,ν2).
Значит, для расстояния Хаусдорфа Hp имеем
Hp(Π(μ1,ν1),Π(μ2,ν2))⩽Wp(μ1,μ2)+Wp(ν1,ν2).
Доказательство. Повторяя приведенные выше рассуждения с h=dp, сначала найдем меру σ0∈Π(μ2,ν1) с Wp(σ1,σ0)⩽Wp(μ1,μ2), а затем возьмем меру σ2∈Π(μ2,ν2) с Wp(σ0,σ2)⩽Wp(ν1,ν2).
Напомним, что вполне регулярное пространство называется секвенциально прохоровским, если всякая последовательность вероятностных радоновских мер на этом пространстве, слабо сходящаяся к радоновской мере, равномерно плотна. Например, полные метрические пространства секвенциально прохоровские (см. [11] или [12] об этом свойстве).
Теорема 2. Пусть X и Y – вполне регулярные пространства. Предположим, что последовательность мер μn∈Pr(X) слабо сходится к мере μ∈Pr(X), последовательность мер νn∈Pr(Y) слабо сходится к мере ν∈Pr(Y), обе последовательности равномерно плотны (что выполнено автоматически, если оба пространства секвенциально прохоровские), а непрерывные функции hn:X×Y→[0,+∞) сходятся к функции h:X×Y→[0,+∞) равномерно на компактах. Предположим также, что даны такие неотрицательные борелевские функции an∈L1(μn) и bn∈L1(νn), что
В частности, это верно, если supn(‖ для некоторого p>1.
Доказательство. Сначала рассмотрим случай, когда h_n\equiv h\leqslant 1. Поскольку функция h непрерывна на компактах (на всем пространстве она не обязательно является непрерывной), то существуют оптимальные меры \sigma_n\in \Pi(\mu_n,\nu_n) для функции h (см. [17; комментарии после теоремы 1.2.1]). По нашему предположению последовательности мер \mu_n и \nu_n равномерно плотны, что влечет равномерную плотность последовательности мер \sigma_n. Перейдем к подпоследовательности индексов n_i, для которой числа K_h(\mu_{n_i},\nu_{n_i}) стремятся к \liminf_{n\to\infty} K_h(\mu_n,\nu_n). Эта последовательность имеет некоторую предельную точку \sigma в слабой топологии, более того, \sigma\in \Pi(\mu,\nu). Интеграл от h относительно меры \sigma является предельной точкой интегралов по мерам \sigma_{n_i}, т.е. чисел K_h(\mu_{n_i},\nu_{n_i}). Значит, выполнено неравенство K_h(\mu,\nu)\leqslant \liminf_{n\to\infty} K_h(\mu_n,\nu_n).
Покажем, что K_h(\mu,\nu)\geqslant \limsup_{n\to\infty} K_h(\mu_n,\nu_n). Пусть \varepsilon>0. Перейдя к подпоследовательности, можно считать, что
Есть несколько известных способов это сделать. Один таков. Вложим наши пространства X и Y гомеоморфно в локально выпуклые пространства X_1 и Y_1 (или просто наделим их равномерностями, порождающими исходные топологии). Тогда функция h равномерно непрерывна на S. Она имеет равномерно непрерывное продолжение g\colon X_1\times Y_1\to [0,1] (это верно для всякого подмножества, см. [29]; но для компактного подмножества достаточно взять непрерывное продолжение на компактификацию Стоуна–Чеха пространства X_1\times Y_1). Функция
поэтому |g(x,y)-g(x_1,y_1)|\leqslant \alpha(x,x_1)+\beta(y,y_1).
Рассмотрим псевдометрики Канторовича d_{\mathrm K,\alpha} и d_{\mathrm K,\beta} на X и Y, порожденные этими псевдометриками \alpha и \beta. Для всех достаточно больших n имеем
Теперь рассмотрим случай различных h_n, но с общей оценкой h_n\leqslant 1. Пусть \varepsilon>0. Используя то же самое компактное множество S, что и выше, найдем такое N, что |h_n-h|\leqslant \varepsilon для всех (x,y)\in S и n\geqslant N. Тогда
поэтому разность между интегралами от h_n и h относительно меры \eta не превосходит 2\varepsilon. Это дает наше утверждение в случае равномерно ограниченных h_n. Общий случай легко сводится к этому случаю, так как для функций \min (h_n,R) имеем
для всех мер \eta\in\Pi(\mu_n,\nu_n) и аналогично для тройки (h,\mu,\nu). Благодаря этой оценке и (3.4) можно выбрать R столь большим, что правая часть предыдущего неравенства будет меньше заданного числа \varepsilon. Следовательно,
Теорему 2 можно сравнить с теоремой 5.20 из [43], где пространства польские, маргиналы \mu_n и \nu_n слабо сходятся к \mu и \nu соответственно, функции стоимости h_n равномерно сходятся к h, а заключение состоит в том, что последовательность оптимальных планов \pi_n для (\mu_n,\nu_n,h_n) содержит подпоследовательность, слабо сходящуюся к оптимальному плану для (\mu,\nu,h).
Следствие 2. Если в ситуации теоремы 2 оптимальные планы \sigma_n и \sigma для троек (\mu_n,\nu_n,h_n) и (\mu,\nu,h) единственны, то меры \sigma_n слабо сходятся к \sigma.
Доказательство. Последовательность мер \sigma_{n} равномерно плотна в силу равномерной плотности маргиналов этих мер. Значит, она имеет слабо сходящуюся поднаправленность (но не обязательно подпоследовательность, ибо пространства не предполагаются метризуемыми). Заметим, что если ее поднаправленность \{\sigma_\alpha\} слабо сходится к мере \sigma_0, то \sigma_0 оптимальна для h. В самом деле, покажем, что
так как меры \sigma_n равномерно плотны и функции \min(h_n,R) сходятся к \min(h,R) равномерно на компактах.
Правая часть (3.6) равна K_{h}(\mu,\nu) по теореме 2. Значит, мера \sigma_0 оптимальна, а тогда \sigma_0=\sigma в силу единственности. Таким образом, последовательность \{\sigma_n\} не имеет отличных от \sigma предельных точек. Следовательно, меры \sigma_{n} слабо сходятся к \sigma.
Следствие доказано.
Для несчетных семейств мер \mu_t, \nu_t и функций a_t\in L^1(\mu_t), b_t\in L^1(\nu_t) аналог (3.4) записывается как
для некоторого p>1. Более общее достаточное условие таково: найдется неограниченная возрастающая функция V>0 на [0,\infty), для которой интегралы от a_t V(a_t) относительно \mu_t и интегралы от b_t V(b_t) относительно \nu_t равномерно ограничены.
Следствие 3. Пусть X и Y – секвенциально прохоровские вполне регулярные пространства и T – топологическое пространство. Предположим, что отображения
секвенциально непрерывны и (t,x,y)\mapsto h_t(x,y), T\times X\times Y\to [0,+\infty) является непрерывной функцией. Предположим также, что существуют такие неотрицательные борелевские функции a_t\in L^1(\mu_t) и b_t\in L^1(\nu_t), что выполнено (3.7). Тогда функция t\mapsto K_{h_t}(\mu_t,\nu_t) секвенциально непрерывна.
Для доказательства достаточно воспользоваться тем, что h_{t_n}(x,y)\to h_t(x,y) равномерно на компактах, если t_n\to t.
Ясно, что прохоровское свойство можно заменить равномерной плотностью семейств \{\mu_t\} и \{\nu_t\}.
Предположим теперь, что X, Y, T – метрические пространства, X и Y полны и для каждого t\in T даны меры \mu_t\in \mathcal{P}_r(X) и \nu_t\in \mathcal{P}_r(Y), причем отображения t\mapsto \mu_t и t\mapsto \nu_t непрерывны в слабой топологии (что равносильно непрерывности по метрике Канторовича–Рубинштейна). Предположим также, что имеется непрерывная неотрицательная функция (t,x,y)\mapsto h_t(x,y). Возникает вопрос: можно ли выбрать оптимальный план, непрерывно зависящий от параметра t? Оказывается, что такой выбор не всегда возможен, как показывают примеры ниже. Однако положение улучшается для приближенных оптимальных планов или в случае единственных оптимальных планов. Для заданного \varepsilon>0 меру \sigma\in \Pi(\mu,\nu) будем называть \varepsilon-оптимальной для функции стоимости h, если
Теорема 3. Предположим, что для каждого t существуют такие неотрицательные борелевские функции a_t\in L^1(\mu_t) и b_t\in L^1(\nu_t), что выполнено (3.7). Тогда можно выбрать \varepsilon-оптимальные меры \sigma_t^\varepsilon\in \Pi(\mu_t,\nu_t) для функций стоимости h_t так, что они будут непрерывными по t в слабой топологии для всякого фиксированного \varepsilon>0.
Если для каждого t имеется единственный оптимальный план \sigma_t, то он непрерывен по t.
Доказательство. Наше предположение влечет включение h_t\in L^1(\sigma) для всех \sigma\in \Pi(\mu_t,\nu_t). Множество
выпукло и компактно в \mathcal{P}_r(X\times Y). В самом деле, выпуклость очевидна, а компактность вытекает из того факта, что M_t замкнуто в компактном множестве \Pi(\mu_t,\nu_t), что проверяется следующим образом. Предположим, что меры \pi_n\in M_t слабо сходятся к мере \pi. Для сходимости интегралов от непрерывных функций h_t по мерам \pi_n к интегралам по мере \pi достаточно проверить (см. [12; теорема 2.7.1]), что
Чтобы найти непрерывную выборку, мы проверим предположения классической теоремы Майкла о селекции, применяемой к многозначному отображению t\mapsto \Psi_t с выпуклыми компактными значениями в полном метризуемом подмножестве M=\mathcal{P}_r(X\times Y) локально выпуклого пространства \mathcal{M}_r(X\times Y). В статье [33] соответствующий результат установлен для множества M в метризуемом локально выпуклом пространстве (см. также [38; с. 41, теорема 1.5] про такой случай), его можно применить, если пространство \mathcal{M}_r(X\times Y) наделить нормой Канторовича–Рубинштейна. Однако теорема 1.2 из [34] покрывает случай, когда M лежит в полном локально выпуклом пространстве, что в нашей ситуации можно применить к пополнению пространства \mathcal{M}_r(X\times Y) со слабой топологией. Для проверки выполнения условий теоремы Майкла нам надо еще показать, что для каждого открытого множества U в \mathcal{P}_r(X\times Y) множество
Можно считать, что это неравенство строгое, поскольку в противном случае можно взять выпуклую линейную комбинацию с оптимальной мерой \sigma_0\in M_w, что для \alpha\in (0,1) даст меру (1-\alpha) \sigma + \alpha \sigma_0\in U с меньшим интегралом от функции h_w. Итак,
Далее, из слабой непрерывности \mu_t и \nu_t по t следует, что последовательности \{\mu_{w_n}\} и \{\nu_{w_n}\} равномерно плотны (напомним, что здесь пространства метризуемы). Это влечет равномерную плотность объединения множеств \Pi(\mu_{w_n}, \nu_{w_n}). Значит, имеется такой компакт S\subset X\times Y, что
В силу непрерывности h на компактном множестве (\{w_n\}\cup \{w\})\times S существует липшицева функция (t,x,y)\mapsto L_t(x,y) на T\times X\times Y со значениями в [0,R], для которой
Пусть L>0 – постоянная Липшица этой функции. Поскольку компактное множество \Pi(\mu_t, \nu_t) зависит от t непрерывно в метрике Хаусдорфа, то для всех достаточно больших n имеем
Поскольку K_{h_{w_n}}(\mu_{w_n},\nu_{w_n})\to K_{h_{w}}(\mu_{w},\nu_{w}), приходим к противоречию с (3.8). Итак, множество \{t\colon \Psi(t) \cap U \ne \varnothing\} открыто.
Наконец, в случае единственности оптимальных планов наш результат вытекает из следствия 2.
Можно было бы сказать, что непрерывность оптимальной стоимости по параметру есть нечто ожидаемое, но как насчет непрерывности оптимальных планов? В общем случае нет непрерывной выборки оптимальных планов.
Пример. Пусть \mu и \nu – такие борелевские вероятностные меры на полном сепарабельном метрическом пространстве X, что существуют по крайней мере два непрерывных отображения f, g\colon X\to X, переводящих \mu в \nu, для которых образы \mu при отображениях x\mapsto (x,f(x)) и x\mapsto (x,g(x)) различны. Например, в качестве \mu и \nu можно взять меру Лебега на [0,1]. В качестве пространства T возьмем последовательность точек 1/n и 0. Положим
Пусть \mu_t=\mu, \nu_t=\nu. При t=0 все меры в \Pi(\mu,\nu) оптимальны для функции стоимости h_0. При t=(2n-1)^{-1} единственная оптимальная мера для функции стоимости h_t равна образу меры \mu при отображении x\mapsto (x,f(x)), а для t=(2n)^{-1} единственная оптимальная мера для функции стоимости h_t равна образу \mu при отображении x\mapsto (x,g(x)). Тогда указанные оптимальные меры не имеют предела при t\to 0. Например, если \mu=\nu есть мера Лебега на [0,1], то можно взять f(x)=x, g(x)=1-x. Вот другой явный пример с теми же самыми мерами на [0, 1]:
Тогда при t>0 оптимальный план сосредоточен на диагонали x=y, при t<0 оптимальный план сосредоточен на диагонали x+y=1, во всех случаях он является нормированной линейной мерой Лебега. При t=0 нет единственности: имеется много оптимальных планов, сосредоточенных на объединении диагоналей, некоторые из них не порождаются отображениями Монжа. Однако левые и правые пределы оптимальных планов порождаются отображениями Монжа.
Стоит отметить, что близкая, но не эквивалентная задача была рассмотрена в [5] и [24] (см. также [42]), где была доказана сходимость по мере оптимальных отображений Монжа T_\varepsilon\to T_0 при \varepsilon\to 0 для некоторой специальной функции стоимости h_\varepsilon. В нашей ситуации соответствующий вопрос ставился бы о некотором виде непрерывности оптимальных отображений Монжа T_t в случае, когда мера \mu_t не зависит от t или все меры \mu_t абсолютно непрерывны относительно некоторой базовой меры \lambda. Конечно, в случае неединственности T_t вопрос возникает о непрерывных выборках оптимальных отображений. Приведем общий результат в духе цитированных работ в случае, когда оптимальные планы порождены отображениями. Частная ситуация, рассмотренная в цитированных работах (см. также [42; теорема 2.53]), имеет дело с постоянными маргиналами и функцией стоимости h_t(x,y)=|x-y|+t |x-y|^2 или более общим семейством функций с тем свойством, что h_t(x,y)=|x-y|+t \varphi(|x-y|)+o(t) со строго выпуклой функцией \varphi. Однако эта специальная структура нужна, чтобы гарантировать существование и единственность решения Монжа. В следующем абстрактном результате существование “отображений Монжа” является частью предположений.
Напомним, что сходимость по мере для отображений во вполне регулярные пространства определяется следующим образом. Пусть \{d_\alpha\} – семейство псевдометрик, задающих топологию вполне регулярного пространств X и \mu\in\mathcal{P}_r(X). Борелевские отображения T_n\colon X\to X сходятся по мере \mu к борелевскому отображению T\colon X\to X, если для каждого \alpha и каждого \delta>0 имеем \mu(x\colon d_\alpha(T_n(x),T(x))\geqslant \delta)\to 0 при n\to \infty.
Предложение. Пусть X – вполне регулярное пространство, \mu_n\in\mathcal{P}_r(X) для n\in \mathbb{Z}^+ и \mu_n\to\mu_0 по вариации. Предположим, что имеются такие борелевские отображения T_n\colon X\to X, n\geqslant 0, что меры \sigma_n, равные образам мер \mu_n при отображениях x\mapsto (x,T_n(x)), слабо сходятся к \sigma_0. Кроме того, предположим, что меры \mu_0 и \nu_0:=\mu_0\circ T_0^{-1} сосредоточены на счетных объединениях метризуемых компактов (что выполнено автоматически, если X – суслинское пространство). Тогда отображения T_n сходятся к T_0 по мере относительно \mu_0.
Доказательство. Пусть \psi – непрерывная псевдометрика на X. Можно считать, что 0\leqslant \psi\leqslant 1. Предположим сначала, что T_{0} непрерывно. Тогда функция \psi(T_{0}(x),y) непрерывна на произведении X^2, значит, в силу слабой сходимости \sigma_n\to \sigma_{0} имеем
Значит, по определению T_n\to T_{0} по мере \mu_0.
В общем случае можно вложить X гомеоморфно в подходящую степень прямой и считать, что X=\mathbb{R}^\tau. Для заданного \varepsilon>0 можно найти компакт K с \mu_0(K)>1-\varepsilon, на котором T_{0} непрерывно. Для этого найдем метризуемый компакт Q с \nu_0(Q)>1-\varepsilon и в борелевском множестве T_0^{-1}(Q) выберем метризуемое компактное множество Q_1 с \mu_0(Q_1)>1-\varepsilon, что возможно, так как \mu_0(T_0^{-1}(Q))=\nu_0(Q)>1\,{-}\,\varepsilon. Остается применить теорему Лузина к борелевскому отображению T_0 между метризуемыми компактами Q_1 и Q. На самом деле для применения обобщения теоремы Лузина достаточно, чтобы только мера \nu_0 была сконцентрирована на метризуемых компактах (см. [11; теорема 7.1.13]). Далее, существует непрерывное отображение S, которое совпадает с T_{0} на K (здесь используется, что X=\mathbb{R}^\tau, поэтому достаточно продолжать компоненты T_0). Повторяя то же самое рассуждение с функцией \psi(S(x),y), получаем, что интегралы от \psi(S(x),T_n(x)) относительно \mu_0 стремятся к интегралу от \psi(S(x), T_0(x)), который оценивается через \varepsilon, поскольку \psi(S(x),T_{0}(x))=0 на K. Это дает сходимость по мере в общем случае.
Предложение доказано.
Заметим, что сходимость отображений T_n к T_0 по мере \mu_0 также достаточна для слабой сходимости мер \sigma_n.
Следствие 4. Предположим, что X – полное сепарабельное метрическое пространство, меры \mu_n\in \mathcal{P}_r(X) сходятся по вариации к мере \mu_0, меры \nu_n\in \mathcal{P}_r(X) слабо сходятся к мере \nu_0, непрерывные функции h_n\geqslant 0 на X\times X сходятся к функции h равномерно на компактах и выполнено (3.4). Предположим также, что для троек (\mu_n,\nu_n,h_n), n\geqslant 0, существуют единственные оптимальные планы Канторовича \sigma_n, которые порождаются единственными оптимальными отображениями Монжа T_n. Тогда отображения T_n сходятся к T_0 по мере \mu_0.
Заметим, что следствие 5.20 в [43] содержит аналогичный результат в случае, когда все меры \mu_n совпадают с \mu_0 и X локально компактно.
Замечание 2. Как показано выше, в ситуации, когда мы имеем дело с задачей оптимальной транспортировки для троек (\mu_n,\nu_n,h_n), предположение, что оптимальные планы \sigma_n слабо сходятся к \sigma_0, выполнено, если оптимальные меры единственны, меры \nu_n сходятся слабо к \nu_0 и равномерно плотны, а функции h_n непрерывны и сходятся к h_0 равномерно на компактах.
Если не предполагается слабая сходимость оптимальных планов, но последовательность \{\nu_n\} слабо сходится и равномерно плотна, то заключение остается в силе для подпоследовательности в \{T_n\}, выбранной так, что соответствующие планы сходятся, если известно, что все оптимальные планы для h_0 (или хотя бы лежащие в замыкании рассматриваемой подпоследовательности планов) порождаются отображениями Монжа.
Частный случай постоянных маргиналов \mu и \nu не намного проще, так как в любом случае надо обеспечить сходимость планов и нужно существование отображений Монжа. Было бы интересно изучить приближенные решения Монжа, непрерывно зависящие от параметра. Для этого можно было бы проанализировать конструкции из работ [16] и [35].
где H – борелевская функция на X\times Y\times \mathcal{P}_r(X\times Y) или на X\times Y\times \Pi(\mu,\nu). Функционалы такого типа недавно изучались рядом авторов; см., например, [27], [2], [7], где можно найти дополнительные ссылки. Для упрощения можно предполагать, что функция H ограниченна, более того, поскольку мы обсуждаем свойства непрерывности, можно также считать, что функция H непрерывна. При таких предположениях для всех мер \mu\in\mathcal{P}_r(X) и \nu\in\mathcal{P}_r(Y) существует решение нелинейной задачи Канторовича для функционала J_H, т.е. план \sigma\in \Pi(\mu,\nu), для которой величина J_H(\sigma) достигает минимума K_H(\mu,\nu) на \Pi(\mu,\nu); см. [20], [13].
Однако следует отметить, что непрерывность по \sigma не выполняется в типичных примерах из цитированных работ, где H(x,y,\sigma)=H(x,\sigma^x) и (\sigma^x)_{x\in X} есть дезинтегрирование меры \sigma относительно ее проекции \sigma_X на X (равной \mu, если \sigma\in\Pi(\mu,\nu)). Напомним, что дезинтегрирование есть такое борелевское отображение x\mapsto \sigma^x из X в \mathcal{P}_r(Y), что \sigma (dx\, dy)=\sigma^x(dy)\, \sigma_X(dx) в смысле равенства
Общие нелинейные функционалы стоимости с параметром будут изучаться в другой работе, а здесь мы представим лишь следующий результат о непрерывности по параметру.
Пусть T – топологическое пространство, X, Y – секвенциально прохоровские пространства (скажем, метрические пространства). Предположим, что
является ограниченной непрерывной функцией. Пусть H_t(x,y,\sigma)\,{:=}\,H(x,y,\sigma,t). Для заданных мер \mu\in \mathcal{P}_r(X) и \nu\in \mathcal{P}_r(Y) положим
Сначала рассмотрим случай, когда функция H не зависит от t.
Теорема 4. Предположим, что меры \mu_n\in\mathcal{P}_r(X) слабо сходятся к мере \mu\in\mathcal{P}_r(X), меры \nu_n\in\mathcal{P}_r(Y) слабо сходятся к мере \nu\in\mathcal{P}_r(Y). Тогда K_H(\mu,\nu)=\lim_{n\to\infty} K_H(\mu_n,\nu_n).
Кроме того, всякая мера \sigma_0\in \mathcal{P}_r(X\times Y), являющаяся предельной точкой последовательности оптимальных мер \sigma_n для маргиналов \mu_n, \nu_n, служит оптимальной для маргиналов \mu,\nu.
Доказательство. Можно считать, что 0\leqslant H\leqslant 1. Так как по предположению имеет место слабая сходимость \mu_{n}\to\mu и \nu_{n}\to\nu, а пространства X и Y секвенциально прохоровские, то обе последовательности равномерно плотны. Поэтому равномерно плотно объединение всех множеств \Pi(\mu_{n}, \nu_{n}) и \Pi(\mu, \nu). Значит, оно содержится в компактном равномерно плотном множестве S в \mathcal{P}_r(X\times Y). Возьмем оптимальные меры \sigma_n для пар (\mu_n,\nu_n) и какую-либо предельную меру \sigma_0 последовательности \{\sigma_n\} в компакте S. Ясно, что \sigma\in \Pi(\mu,\nu). Покажем, что J_H(\sigma_0)=K_H(\mu,\nu)=\lim_{n\to\infty}K_H(\mu_n,\nu_n). Сначала проверим последнее равенство.
Пусть \varepsilon>0. Возьмем такой компакт K в X\times Y, что
\begin{equation*}
\sigma(K)>1-\varepsilon \quad \forall\, \sigma\in S.
\end{equation*}
\notag
Найдутся такие ограниченные непрерывные псевдометрики d_1 на X и d_2 на Y, что существует функция G\colon X\times Y\times \mathcal{P}_r(X\times Y)\to [0,1], липшицева с некоторой постоянной L относительно псевдометрики {d_1\oplus d_2\oplus d_{\mathrm K}}, где d_{\mathrm K} – псевдометрика Канторовича, порожденная d_1\oplus d_2, что
Это можно сделать с помощью теоремы Вейерштрасса, ибо класс функций такого вида есть алгебра, разделяющая точки компакта K\times S. В силу слабой сходимости имеем d_{\mathrm K,d_1}(\mu_n,\mu)\to 0, d_{\mathrm K,d_2}(\nu_n,\nu)\to 0. Поэтому при некотором N для всех n\geqslant N справедливы неравенства
Поэтому K_G(\mu_n,\nu_n) \leqslant K_G(\mu,\nu)+2\varepsilon. Аналогично получаем оценку K_G(\mu,\nu) \leqslant K_G(\mu_n,\nu_n)+2\varepsilon, взяв оптимальный план \pi_n' для тройки (G,\mu_n,\nu_n) и подобрав \pi'\in \Pi(\mu,\nu) с d_{\mathrm K, d_1\oplus d_2}(\pi',\pi_n')\leqslant \varepsilon/(L+1).
Для исходной функции H выполнены оценки |J_H(\sigma)-J_G(\sigma)|\leqslant 3\varepsilon для всех \sigma\in S, поэтому
при n\geqslant N. Итак, |K_H(\mu,\nu)-K_H(\mu_n,\nu_n)|\leqslant 8\varepsilon, что в силу произвольности \varepsilon доказывает равенство K_H(\mu,\nu)=\lim_{n\to\infty} K_H(\mu_n,\nu_n).
Теперь покажем, что \sigma_0 – оптимальный план для маргиналов \mu, \nu. Пусть \varepsilon>0. Возьмем тот же компакт K\subset X\times Y, что и выше, а также те же псевдометрику d_1\oplus d_2 и L-липшицеву по этой псевдометрике функцию стоимости G. Мера \sigma_0 является предельной для последовательности \{\sigma_n\} в слабой топологии, следовательно, есть подпоследовательность \sigma_{n_k}, сходящаяся к \sigma_0 по указанной псевдометрике. Без потери общности можно считать, что такова вся исходная последовательность. Значит, найдется такое N, что d_{\mathrm K,d_1\oplus d_2}(\sigma_{n},\sigma_0)\leqslant \varepsilon/(L+1) при n\geqslant N. Тогда при таких n имеем
т.е. |K_H(\mu_n,\nu_n)\,{-}\,J_H(\sigma_0)|\leqslant 8\varepsilon. Поэтому |K_H(\mu,\nu)\,{-}\,J_H(\sigma_0)|\leqslant 8\varepsilon, что завершает доказательство.
Вернемся к рассмотрению функции стоимости H, зависящей от t.
Теорема 5. Предположим, что t\mapsto \mu_t и t\mapsto \nu_t есть непрерывные отображения со значениями в \mathcal{P}_r(X) и \mathcal{P}_r(Y) соответственно. Тогда функция t\mapsto K(t,\mu_t,\nu_t) секвенциально непрерывна.
Доказательство. Сведем наше утверждение к теореме 4. При этом воспользуемся тем же начальным построением. Можно считать, что 0\leqslant H\leqslant 1. Пусть t_n\to t_0. Так как по предположению имеет место слабая сходимость \mu_{t_n}\to\mu_{t_0} и \nu_{t_n}\to\nu_{t_0}, а пространства X и Y секвенциально прохоровские, то обе последовательности равномерно плотны. Поэтому равномерно плотно объединение множеств мер \Pi(\mu_{t_n}, \nu_{t_n}), n\geqslant 0. Значит, оно содержится в компактном равномерно плотном множестве S в \mathcal{P}_r(X\times Y). Пусть \varepsilon>0. Опять возьмем такой компакт K в X\times Y, что
\begin{equation*}
\sigma(K)>1-\varepsilon \quad \forall\, \sigma\in S.
\end{equation*}
\notag
В силу непрерывности H и компактности K\times S найдется такое N, что
Как и выше, этот результат легко переносится на случай неограниченных неотрицательных функций стоимости с подходящей равномерной интегрируемостью. Например, достаточно иметь оценку
Заметим, что в случае функционала стоимости, порожденного функцией вида H(x,\sigma_x) (упомянутого выше), в работе [7] доказано, что функционал J_H(\sigma) полунепрерывен снизу по \sigma при условии, что функция H совместно полунепрерывна снизу и выпукла по второму аргументу. Как показано в [2] (см. также [13]), это может быть неверно без выпуклости.
§ 5. Приложение к дискретным приближениям
Пусть X и Y – компактные метрические пространства с борелевскими вероятностными мерами \mu и \nu и h – липшицева функция на X\times Y. Пусть L – ее постоянная Липшица. Полученная выше оценка может быть использована для изучения дискретных приближений задачи Канторовича для тройки (\mu,\nu,h). Рассмотрим дискретные приближения маргиналов \mu и \nu мерами вида
соответственно, где A_1,\dots,A_n и B_1,\dots,B_n – разбиения X и Y на дизъюнктные борелевские множества диаметра не более заданного числа \varepsilon, a_i\in A_i, b_i\in B_i. Поскольку для всяких функций f\in \mathrm{Lip}_1(X), g\in \mathrm{Lip}_1(Y) справедливы неравенства |f(x)-f(a_i)|\leqslant \varepsilon для всех x\in A_i и |g(y)-g(b_i)|\leqslant \varepsilon для всех y\in B_i, то получаем оценки
Взяв в качестве \sigma_2 оптимальную меру для тройки (h,\mu_n,\nu_n) и найдя с помощью нашей оценки каплинг \sigma_1\in \Pi(\mu,\nu) с \|\sigma_1-\sigma_2\|_{\mathrm K}\leqslant 2\varepsilon, получаем
Аналогично доказываем, что K_h(\mu_n,\nu_n)-K_h(\mu,\nu)\leqslant 2L \varepsilon, взяв сначала план \sigma_1 оптимальным для тройки (h,\mu,\nu) и найдя подходящий каплинг \sigma_2\in \Pi(\mu_n,\nu_n). Следовательно,
Задача Канторовича для (\mu_n,\nu_n,h) конечномерна. Соответствующая матрица есть (h(a_i,b_j))_{i,j\leqslant n}. Если \sigma_n – оптимальный план для этой задачи, то наша оценка показывает, что имеется каплинг \pi_n\in \Pi(\mu,\nu) для исходной задачи с \|\sigma_n-\pi_n\|_{\mathrm K}\leqslant 2\varepsilon. Непосредственно проверяется, что каплинг \pi_n является 4L\varepsilon-приближенным для исходной задачи. Конечно, важно контролировать n: наименьшее возможное n зависит от метрической энтропии X и Y. Например, если X и Y лежат в \mathbb{R}^d, то n имеет порядок \varepsilon^{-d}. Скорость приближения мер \mu и \nu дискретными мерами есть, конечно, отдельный вопрос (см. [28]). Если h не липшицева, то аналогичная оценка может быть получена в терминах модуля непрерывности h. Наш результат может быть использован в аппроксимационных схемах, обсуждаемых в [6].
Наконец, случай некомпактных пространств в принципе сводится к рассмотренному, если имеется некоторая информация о компактных множествах, на которых \varepsilon-сосредоточены маргиналы. Компактность была нужна лишь для того, чтобы получить унифицированные дискретные приближения маргиналов и использовать расстояние d_{\mathrm K} вместо d_{\mathrm{KR}}. Дискретные приближения мультимаргинальной задачи Канторовича строятся аналогично, наше неравенство дает довольно универсальные оценки соответствующих стоимостей и оптимальных планов. Недавние результаты о дискретных мультимаргинальных задачах Канторовича см. в работе [44].
Такой же подход работает для нелинейного функционала стоимости, рассмотренного нами в § 4. Предположим, что h – неотрицательная функция на X\times Y\times \mathcal{P}_r(X\times Y), липшицева с постоянной L. Тогда нелинейный функционал
Следовательно, используя те же самые дискретные \varepsilon-приближения маргиналов \mu и \nu, что и выше, получаем оценку |K_h(\mu,\nu)-K_h(\mu_n,\nu_n)|\leqslant 4\varepsilon L.
Этот же метод работает для некомпактных метрических пространств X и Y при условии, что имеются дискретные \varepsilon-приближения маргиналов в метрике Канторовича–Рубинштейна d_{\mathrm{KR}} и функция h является L-липшицевой и ограничена числом L. Тогда заключительная оценка остается той же самой. Если \sigma_n – оптимальный план для (\mu_n,\nu_n), то наша оценка обеспечивает план \pi_n\in \Pi(\mu,\nu) с d_{\mathrm{KR}}(\sigma_n, \pi_n)\leqslant 2\varepsilon. Этот план является 8L\varepsilon-приближенным для пары (\mu,\nu).
Благодарность
Благодарим рецензента за полезные замечания.
Список литературы
1.
K. A. Afonin, V. I. Bogachev, “Kantorovich type topologies on spaces of measures and convergence of barycenters”, Commun. Pure Appl. Anal., 22:2 (2023), 597–612
2.
J.-J. Alibert, G. Bouchitté, T. Champion, “A new class of costs for optimal transport planning”, European J. Appl. Math., 30:6 (2019), 1229–1263
3.
L. Ambrosio, E. Brué, D. Semola, Lectures on optimal transport, Unitext, 130, La Mat. per il 3+2, Springer, Cham, 2021, ix+250 pp.
4.
L. Ambrosio, N. Gigli, “A user's guide to optimal transport”, Modelling and optimisation of flows on networks, Lecture Notes in Math., 2062, Fond. CIME/CIME Found. Subser., Springer, Heidelberg, 2013, 1–155
5.
L. Ambrosio, A. Pratelli, “Existence and stability results in the L^1 theory of optimal transportation”, Optimal transportation and applications (Martina Franca, 2001), Lecture Notes in Math., 1813, Springer, Berlin, 2003, 123–160
6.
M. L. Avendaño-Garrido, J. R. Gabriel-Argüelles, L.-T. Quintana, J. González-Hernández, “An approximation scheme for the Kantorovich–Rubinstein problem on compact spaces”, J. Numer. Math., 26:2 (2018), 63–75
7.
J. Backhoff-Veraguas, M. Beiglböck, G. Pammer, “Existence, duality, and cyclical monotonicity for weak transport costs”, Calc. Var. Partial Differential Equations, 58:6 (2019), 203, 28 pp.
8.
J. Backhoff-Veraguas, G. Pammer, “Applications of weak transport theory”, Bernoulli, 28:1 (2022), 370–394
9.
J. Bergin, “On the continuity of correspondences on sets of measures with restricted marginals”, Econom. Theory, 13:2 (1999), 471–481
10.
S. Bobkov, M. Ledoux, One-dimensional empirical measures, order statistics, and Kantorovich transport distances, Mem. Amer. Math. Soc., 261, no. 1259, Amer. Math. Soc., Providence, RI, 2019, v+126 pp.
11.
В. И. Богачев, Основы теории меры, т. 1, 2, изд. 3-е, испр. и доп., НИЦ “Регулярная и хаотическая динамика”, Ин-т компьютерных исследований, М.–Ижевск, 2021, 584 с., 688 с.; англ. пер.: V. I. Bogachev, Measure theory, т. I, II, Springer-Verlag, Berlin, 2007, xviii+500 pp., xiv+575 с.
12.
V. I. Bogachev, Weak convergence of measures, Math. Surveys Monogr., 234, Amer. Math. Soc., Providence, RI, 2018, xii+286 pp.
13.
В. И. Богачев, “Задача Канторовича оптимальной транспортировки мер: новые направления исследований”, УМН, 77:5(467) (2022), 3–52; англ. пер.: V. I. Bogachev, “Kantorovich problem of optimal transportation of measures: new directions of research”, Russian Math. Surveys, 77:5 (2022), 769–817
14.
В. И. Богачев, “Задачи Канторовича с параметром и ограничениями на плотности”, Сиб. матем. журн., 63:1 (2022), 42–57; англ. пер.: V. I. Bogachev, “Kantorovich problems with a parameter and density constraints”, Siberian Math. J., 63:1 (2022), 34–47
15.
В. И. Богачев, А. Н. Доледенок, И. И. Малофеев, “Задача Канторовича с параметром и ограничениями на плотность”, Матем. заметки, 110:6 (2021), 922–926; англ. пер.: V. I. Bogachev, A. N. Doledenok, I. I. Malofeev, “The Kantorovich problem with a parameter and density constraints”, Math. Notes, 110:6 (2021), 952–955
16.
В. И. Богачев, А. Н. Калинин, С. Н. Попова, “О равенстве значений в задачах Монжа и Канторовича”, Вероятность и статистика. 25, Посвящается памяти Владимира Николаевича Судакова, Зап. науч. сем. ПОМИ, 457, ПОМИ, СПб., 2017, 53–73; англ. пер.: V. I. Bogachev, A. N. Kalinin, S. N. Popova, “On the equality of values in the Monge and Kantorovich problems”, J. Math. Sci. (N.Y.), 238:4 (2019), 377–389
17.
В. И. Богачев, А. В. Колесников, “Задача Монжа–Канторовича: достижения, связи и перспективы”, УМН, 67:5(407) (2012), 3–110; англ. пер.: V. I. Bogachev, A. V. Kolesnikov, “The Monge–Kantorovich problem: achievements, connections, and perspectives”, Russian Math. Surveys, 67:5 (2012), 785–890
18.
V. I. Bogachev, I. I. Malofeev, “Kantorovich problems and conditional measures depending on a parameter”, J. Math. Anal. Appl., 486:1 (2020), 123883, 30 pp.
19.
В. И. Богачев, С. Н. Попова, “О задаче Канторовича с параметром”, Докл. РАН. Матем., информ., проц. упр., 507 (2022), 26–28; англ. пер.: V. I. Bogachev, S. N. Popova, “On Kantorovich problems with a parameter”, Dokl. Math., 106:3 (2022), 426–428
20.
В. И. Богачев, А. В. Резбаев, “Существование решений нелинейной задачи Канторовича оптимальной транспортировки”, Матем. заметки, 112:3 (2022), 360–370; англ. пер.: V. I. Bogachev, A. V. Rezbaev, “Existence of solutions to the nonlinear Kantorovich transportation problem”, Math. Notes, 112:3 (2022), 369–377
21.
B. Bonnet, H. Frankowska, “Differential inclusions in Wasserstein spaces: the Cauchy–Lipschitz framework”, J. Differential Equations, 271 (2021), 594–637
22.
C. Clason, D. A. Lorenz, H. Mahler, B. Wirth, “Entropic regularization of continuous optimal transport problems”, J. Math. Anal. Appl., 494:1 (2021), 124432, 22 pp.
23.
J. Dedecker, C. Prieur, P. Raynaud De Fitte, “Parametrized Kantorovich–Rubinštein theorem and application to the coupling of random variables”, Dependence in probability and statistics, Lect. Notes Stat., 187, Springer, New York, 2006, 105–121
24.
L. De Pascale, J. Louet, F. Santambrogio, “The Monge problem with vanishing gradient penalization: vortices and asymptotic profile”, J. Math. Pures Appl. (9), 106:2 (2016), 237–279
25.
A. Figalli, F. Glaudo, An invitation to optimal transport, Wasserstein distances, and gradient flows, EMS Textbk. Math., EMS Press, Berlin, 2021, vi+136 pp.
26.
M. Ghossoub, D. Saunders, “On the continuity of the feasible set mapping in optimal transport”, Econ. Theory Bull., 9:1 (2021), 113–117
27.
N. Gozlan, C. Roberto, P.-M. Samson, P. Tetali, “Kantorovich duality for general transport costs and applications”, J. Funct. Anal., 273:11 (2017), 3327–3405
28.
S. Graf, H. Luschgy, Foundations of quantization for probability distributions, Lecture Notes in Math., 1730, Springer-Verlag, Berlin, 2000, x+230 pp.
29.
M. Katětov, “On real-valued functions in topological spaces”, Fund. Math., 38 (1951), 85–91; “Correction”, 40 (1953), 203–205
30.
S. Kuksin, V. Nersesyan, A. Shirikyan, “Exponential mixing for a class of dissipative PDEs with bounded degenerate noise”, Geom. Funct. Anal., 30:1 (2020), 126–187
31.
D. A. Lorenz, P. Manns, C. Meyer, “Quadratically regularized optimal transport”, Appl. Math. Optim., 83:3 (2021), 1919–1949
32.
И. И. Малофеев, “Измеримая зависимость условных мер от параметра”, Докл. РАН, 470:1 (2016), 13–17; англ. пер.: I. I. Malofeev, “Measurable dependence of conditional measures on a parameter”, Dokl. Math., 94:2 (2016), 493–497
33.
E. Michael, “Continuous selections. I”, Ann. of Math. (2), 63:2 (1956), 361–382
34.
E. Michael, “A selection theorem”, Proc. Amer. Math. Soc., 17 (1966), 1404–1406
35.
A. Pratelli, “On the equality between Monge's infimum and Kantorovich's minimum in optimal mass transportation”, Ann. Inst. Henri Poincaré Probab. Stat., 43:1 (2007), 1–13
36.
S. T. Rachev, L. Rüschendorf, Mass transportation problems, v. I, Probab. Appl. (N.Y.), Theory, Springer-Verlag, New York, 1998, xxvi+508 pp. ; v. II, Applications, xxvi+430 pp.
37.
D. Ramachandran, L. Rüschendorf, “A general duality theorem for marginal problems”, Probab. Theory Related Fields, 101:3 (1995), 311–319
38.
D. Repovš, P. V. Semenov, Continuous selections of multivalued mappings, Math. Appl., 455, Kluwer Acad. Publ., Dordrecht, 1998, viii+356 pp.
39.
F. Santambrogio, Optimal transport for applied mathematicians. Calculus of variations, PDEs, and modeling, Progr. Nonlinear Differential Equations Appl., 87, Birkhäuser/Springer, Cham, 2015, xxvii+353 pp.
40.
A. Savchenko, M. Zarichnyi, “Correspondences of probability measures with restricted marginals”, Proc. Intern. Geom. Center, 7:4 (2014), 34–39
41.
A. M. Vershik, P. B. Zatitskiy, F. V. Petrov, “Geometry and dynamics of admissible metrics in measure spaces”, Cent. Eur. J. Math., 11:3 (2013), 379–400
42.
C. Villani, Topics in optimal transportation, Grad. Stud. Math., 58, Amer. Math. Soc., Providence, RI, 2003, xvi+370 pp.
43.
C. Villani, Optimal transport. Old and new, Grundlehren Math. Wiss., 338, Springer, New York, 2009, xxii+973 pp.
44.
D. Vögler, “Geometry of Kantorovich polytopes and support of optimizers for repulsive multi-marginal optimal transport on finite state spaces”, J. Math. Anal. Appl., 502:1 (2021), 125147, 31 pp.
45.
Feng-Yu Wang, Jie-Xiang Zhu, “Limit theorems in Wasserstein distance for empirical measures of diffusion processes on Riemannian manifolds”, Ann. Inst. Henri Poincaré Probab. Stat., 59:1 (2023), 437–475
46.
Xicheng Zhang, “Stochastic Monge–Kantorovich problem and its duality”, Stochastics, 85:1 (2013), 71–84
Образец цитирования:
В. И. Богачев, С. Н. Попова, “Расстояния Хаусдорфа между каплингами и оптимальная транспортировка с параметром”, Матем. сб., 215:1 (2024), 33–58; V. I. Bogachev, S. N. Popova, “Hausdorff distances between couplings and optimal transportation”, Sb. Math., 215:1 (2024), 28–51
\RBibitem{BogPop24}
\by В.~И.~Богачев, С.~Н.~Попова
\paper Расстояния Хаусдорфа между каплингами и оптимальная транспортировка с параметром
\jour Матем. сб.
\yr 2024
\vol 215
\issue 1
\pages 33--58
\mathnet{http://mi.mathnet.ru/sm9920}
\crossref{https://doi.org/10.4213/sm9920}
\mathscinet{http://mathscinet.ams.org/mathscinet-getitem?mr=4741221}
\zmath{https://zbmath.org/?q=an:1543.49037}
\adsnasa{https://adsabs.harvard.edu/cgi-bin/bib_query?2024SbMat.215...28B}
\transl
\by V.~I.~Bogachev, S.~N.~Popova
\paper Hausdorff distances between couplings and optimal transportation
\jour Sb. Math.
\yr 2024
\vol 215
\issue 1
\pages 28--51
\crossref{https://doi.org/10.4213/sm9920e}
\isi{https://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcApp=Publons&SrcAuth=Publons_CEL&DestLinkType=FullRecord&DestApp=WOS_CPL&KeyUT=001224793300002}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85193420317}
Образцы ссылок на эту страницу:
https://www.mathnet.ru/rus/sm9920
https://doi.org/10.4213/sm9920
https://www.mathnet.ru/rus/sm/v215/i1/p33
Эта публикация цитируется в следующих 3 статьяx:
Светлана Попова, “Непрерывная выборка приближенных решений Монжа в задаче Канторовича с параметром”, Функц. анализ и его прил., 58:2 (2024), 137–156; Svetlana Popova, “Continuous selection of approximate Monge solutions in the Kantorovich problem with a parameter”, Funct. Anal. Appl., 58:2 (2024), 212–227
С. Н. Попова, “О нелинейных задачах Канторовича для функций стоимости специального вида”, Алгебра и анализ, 36:4 (2024), 165–194
А. Г. Ченцов, Д. А. Серков, “Непрерывная зависимость множеств в пространстве мер и задача на программный минимакс”, Тр. ИММ УрО РАН, 30, № 2, 2024, 277–299; A. G. Chentsov, D. A. Serkov, “Continuous dependence of sets in a space of measures and a program minimax problem”, Proc. Steklov Inst. Math. (Suppl.), 325, suppl. 1 (2024), S76–S98