2025-11-22T16:07:16.564646

Unstable optimal transport maps

Letrouit
The stability of optimal transport maps with respect to perturbations of the marginals is a question of interest for several reasons, ranging from the justification of the linearized optimal transport framework to numerical analysis and statistics. Under various assumptions on the source measure, it is known that optimal transport maps are stable with respect to variations of the target measure. In this note, we focus on the mechanisms that can, on the contrary, lead to instability. We identify two of them, which we illustrate through examples of absolutely continuous source measures $ρ$ in $\mathbb{R}^d$ for which optimal transport maps are less stable, or even very unstable. We first show that instability may arise from the unboundedness of the density: we exhibit a source density on the unit ball of $\mathbb{R}^d$ which blows up superpolynomially at two points of the boundary and for which optimal transport maps are highly unstable. Then we prove that even for uniform densities on bounded open sets, optimal transport maps can be rather unstable close enough to configurations where uniqueness of optimal plans is lost.
academic

অস্থিতিশীল সর্বোত্তম পরিবহন মানচিত্র

মৌলিক তথ্য

  • কাগজ ID: 2510.13265
  • শিরোনাম: অস্থিতিশীল সর্বোত্তম পরিবহন মানচিত্র
  • লেখক: সিরিল লেট্রুইট (ইউনিভার্সিটে প্যারিস-স্যাকলে, CNRS)
  • শ্রেণীবিভাগ: math.OC (গাণিতিক অপ্টিমাইজেশন এবং নিয়ন্ত্রণ)
  • প্রকাশনার সময়: ২০২৫ সালের ১৫ অক্টোবর
  • কাগজের লিঙ্ক: https://arxiv.org/abs/2510.13265

সারসংক্ষেপ

সর্বোত্তম পরিবহন মানচিত্রের প্রান্তিক বিতরণ বিঘ্নের প্রতি স্থিতিশীলতা একটি গুরুত্বপূর্ণ সমস্যা, যা রৈখিক সর্বোত্তম পরিবহন কাঠামোর যথার্থতা, সংখ্যাসূচক বিশ্লেষণ এবং পরিসংখ্যান সহ একাধিক ক্ষেত্র জড়িত। উৎস পরিমাপের বিভিন্ন অনুমানের অধীনে, সর্বোত্তম পরিবহন মানচিত্র লক্ষ্য পরিমাপের পরিবর্তনের প্রতি স্থিতিশীল বলে পরিচিত। এই কাগজটি বিপরীত প্রক্রিয়া—অস্থিতিশীলতার কারণগুলিতে ফোকাস করে। লেখক দুটি প্রক্রিয়া চিহ্নিত করেছেন এবং Rd\mathbb{R}^d এ পরম ক্রমাগত উৎস পরিমাপ ρ\rho এর উদাহরণের মাধ্যমে সর্বোত্তম পরিবহন মানচিত্রের অস্থিতিশীলতা এমনকি উচ্চ অস্থিতিশীলতা প্রদর্শন করেছেন। প্রথমে প্রমাণ করা হয়েছে যে ঘনত্বের সীমাহীনতা অস্থিতিশীলতার দিকে পরিচালিত করতে পারে: Rd\mathbb{R}^d এর একক গোলকে উৎস ঘনত্ব তৈরি করা হয়েছে, যা সীমানার দুটি বিন্দুতে অতি-বহুপদী বিস্ফোরণ ঘটায়, সংশ্লিষ্ট সর্বোত্তম পরিবহন মানচিত্র অত্যন্ত অস্থিতিশীল। তারপর প্রমাণ করা হয়েছে যে এমনকি সীমাবদ্ধ খোলা সেটে সমান ঘনত্বের জন্যও, যখন সর্বোত্তম পরিকল্পনার অনন্যতা হারানোর কনফিগারেশনের কাছাকাছি যায়, সর্বোত্তম পরিবহন মানচিত্র বেশ অস্থিতিশীল হতে পারে।

গবেষণা পটভূমি এবং প্রেরণা

সমস্যার বর্ণনা

সর্বোত্তম পরিবহন তত্ত্বের একটি মূল সমস্যা হল সর্বোত্তম পরিবহন মানচিত্র TμT_\mu (স্থির উৎস পরিমাপ ρ\rho থেকে লক্ষ্য পরিমাপ μ\mu পর্যন্ত) এর লক্ষ্য পরিমাপ বিঘ্নের প্রতি স্থিতিশীলতা অধ্যয়ন করা। নির্দিষ্টভাবে, ধ্রুবক C,α>0C, \alpha > 0 এবং p1p \geq 1 বিদ্যমান কিনা তা অধ্যয়ন করা, যেমন সংক্ষিপ্ত সেট YRdY \subset \mathbb{R}^d এ যেকোনো সম্ভাব্যতা পরিমাপ μ,ν\mu, \nu এর জন্য: TμTνL2(ρ)CWp(μ,ν)α\|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha

গুরুত্ব

এই সমস্যার গুরুত্ব একাধিক দিক থেকে প্রতিফলিত হয়:

  1. রৈখিক সর্বোত্তম পরিবহন কাঠামোর যথার্থতা: মানচিত্র μTμ\mu \mapsto T_\mu Wasserstein স্থান (P(Y),W2)(P(Y), W_2) থেকে Hilbert স্থান L2(ρ,Rd)L^2(\rho, \mathbb{R}^d) এ একটি এমবেডিং প্রদান করে
  2. সংখ্যাসূচক বিশ্লেষণ: যখন লক্ষ্য পরিমাপ μ\mu শুধুমাত্র আনুমানিক μ^\hat{\mu} এর মাধ্যমে পাওয়া যায়, তখন নিশ্চিত করতে হবে যে গণনা করা Tμ^T_{\hat{\mu}} প্রকৃত TμT_\mu এর কাছাকাছি
  3. পরিসংখ্যানগত অনুমান: পরিসংখ্যানগত সর্বোত্তম পরিবহনে, নমুনা থেকে সর্বোত্তম পরিবহন মানচিত্র অনুমান করতে হয়

বিদ্যমান গবেষণার সীমাবদ্ধতা

যদিও উৎস পরিমাপ ρ\rho এর বিভিন্ন অনুমানের অধীনে পরিমাণগত স্থিতিশীলতা অসমতা প্রতিষ্ঠিত হয়েছে, তবে এর আগে কখনো প্রমাণ করা হয়নি যে সর্বোত্তম পরিবহন মানচিত্র অস্থিতিশীল হতে পারে। বিদ্যমান ফলাফল প্রধানত স্থিতিশীলতার ইতিবাচক ফলাফলের উপর দৃষ্টি নিবদ্ধ করে, অস্থিতিশীলতা প্রক্রিয়ার গভীর বোঝার অভাব রয়েছে।

মূল অবদান

  1. প্রথমবার অত্যন্ত অস্থিতিশীল সর্বোত্তম পরিবহন মানচিত্র তৈরি করা: প্রমাণ করা হয়েছে যে পরম ক্রমাগত উৎস পরিমাপ বিদ্যমান, যেমন কোনো ফর্মের পরিমাণগত স্থিতিশীলতা অসমতা বৈধ নয়
  2. দুটি অস্থিতিশীলতা প্রক্রিয়া চিহ্নিত করা:
    • ঘনত্বের অতি-বহুপদী বিস্ফোরণ
    • সর্বোত্তম পরিকল্পনার অ-অনন্যতার কাছাকাছি কনফিগারেশন
  3. নির্দিষ্ট প্রতিউদাহরণ নির্মাণ প্রদান করা: নির্ভুল গাণিতিক নির্মাণের মাধ্যমে অস্থিতিশীলতা ঘটনা প্রদর্শন করা
  4. মানচিত্র μTμ\mu \mapsto T_\mu এর নিয়মিততা সীমাবদ্ধতা প্রকাশ করা: প্রমাণ করা হয়েছে যে কিছু ক্ষেত্রে এই মানচিত্র 12\frac{1}{2}-Hölder ক্রমাগত নয়

পদ্ধতির বিস্তারিত ব্যাখ্যা

কাজের সংজ্ঞা

ρ,μP2(Rd)\rho, \mu \in P_2(\mathbb{R}^d) দেওয়া, যেখানে ρ\rho পরম ক্রমাগত, Brenier উপপাদ্য নিশ্চিত করে যে অনন্য সর্বোত্তম পরিবহন মানচিত্র TμL2(ρ)T_\mu \in L^2(\rho) ρ\rho থেকে μ\mu পর্যন্ত বিদ্যমান। TμT_\mu এর μ\mu বিঘ্নের প্রতি স্থিতিশীলতা অধ্যয়ন করা।

উপপাদ্য 1.1 এর নির্মাণ পদ্ধতি

উৎস পরিমাপ নির্মাণ

A=(1,0,,0)A = (1, 0, \ldots, 0), A=(1,0,,0)A' = (-1, 0, \ldots, 0), E={A,A}E = \{A, A'\} সেট করুন। ফাংশন সংজ্ঞায়িত করুন: f(r)=rdmin(1,(logr)2)f(r) = r^{-d} \min(1, (\log r)^{-2})

ঘনত্ব তৈরি করুন: ρ(x)=c0f(dist(x,E))\rho(x) = c_0 f(\text{dist}(x, E)) যেখানে c0>0c_0 > 0 যাতে ρ\rho একক গোলক BRd(0,1)B_{\mathbb{R}^d}(0,1) এ সম্ভাব্যতা ঘনত্ব হয়।

লক্ষ্য পরিমাপ ক্রম

Bθ=(Rsinθ,Rcosθ,0,,0)B_\theta = (R\sin\theta, R\cos\theta, 0, \ldots, 0) এবং Bθ=(Rsinθ,Rcosθ,0,,0)B'_\theta = (-R\sin\theta, -R\cos\theta, 0, \ldots, 0) বিবেচনা করুন, সংজ্ঞায়িত করুন: μθ=12(δBθ+δBθ)\mu_\theta = \frac{1}{2}(\delta_{B_\theta} + \delta_{B'_\theta})

মূল অনুমান

জ্যামিতিক বিশ্লেষণের মাধ্যমে প্রমাণ করুন: TμθTμ0L2(ρ)2c0c1σd1R2log(θ/4)\|\|T_{\mu_\theta} - T_{\mu_0}\|\|^2_{L^2(\rho)} \geq \frac{c_0 c_1 \sigma_{d-1} R^2}{|\log(\theta/4)|}

যখন Wp(μ0,μθ)Rθ/2W_p(\mu_0, \mu_\theta) \sim R\theta/2, তাই স্থিতিশীলতা অসমতা যেকোনো α>0\alpha > 0 এর জন্য ব্যর্থ হয়।

উপপাদ্য 1.2 এর নির্মাণ পদ্ধতি

সমর্থন সেট নির্মাণ

অসীম সংখ্যক আয়তাকার সমান্তরাল পাইপের জোড়ার ইউনিয়ন তৈরি করুন: X=i=1+SiX = \bigcup_{i=1}^{+\infty} S_i যেখানে Si=T+(Ai+,i,ri)T(Ai,i,ri)S_i = T^+(A_i^+, \ell_i, r_i) \cup T^-(A_i^-, \ell_i, r_i)

প্যারামিটার নির্বাচন

ক্রম নির্বাচন করুন সন্তুষ্ট করতে:

  • ri=c02ir_i = c_0 2^{-i} (অতি-বহুপদী হ্রাস)
  • i=wi=c0c1i2\ell_i = w_i = c_0 c_1 i^{-2}
  • কোষের মধ্যে দূরত্ব কোষের মধ্যে দূরত্বের চেয়ে অনেক বেশি

অস্থিতিশীলতা প্রমাণ

নির্দিষ্ট বিঘ্নের অধীনে সর্বোত্তম পরিবহন মানচিত্র পরিবর্তন বিশ্লেষণের মাধ্যমে, প্রমাণ করুন: TμTνiL2(ρ)2Wp(μ,νi)2αCwi32αpri12αp2α\frac{\|\|T_\mu - T_{\nu_i}\|\|^2_{L^2(\rho)}}{W_p(\mu, \nu_i)^{2\alpha}} \geq C w_i^{3-\frac{2\alpha}{p}} r_i^{1-\frac{2\alpha}{p}-2\alpha}

যখন α>p2(p+1)\alpha > \frac{p}{2(p+1)}, ডান দিক অসীমে প্রবণ হয়।

পরীক্ষামূলক সেটআপ

তাত্ত্বিক যাচাইকরণ

এই কাগজটি বিশুদ্ধ তাত্ত্বিক কাজ, কঠোর গাণিতিক প্রমাণের মাধ্যমে ফলাফল যাচাই করে:

  1. জ্যামিতিক বিশ্লেষণ: সর্বোত্তম পরিবহন মানচিত্রের জ্যামিতিক কাঠামো বিস্তারিত বিশ্লেষণ
  2. পরিমাপ অনুমান: বিভিন্ন অঞ্চলের পরিমাপ নির্ভুলভাবে গণনা করা
  3. দূরত্ব গণনা: Wasserstein দূরত্ব এবং L2L^2 দূরত্ব সঠিকভাবে অনুমান করা

নির্মাণ যাচাইকরণ

নিম্নলিখিত পদক্ষেপের মাধ্যমে নির্মাণের সঠিকতা যাচাই করুন:

  1. প্রমাণ করুন যে নির্মিত ঘনত্ব সত্যিই সম্ভাব্যতা ঘনত্ব
  2. সর্বোত্তম পরিবহন মানচিত্রের স্পষ্ট অভিব্যক্তি যাচাই করুন
  3. জ্যামিতিক কনফিগারেশন প্রয়োজনীয় বিচ্ছিন্নতা শর্ত পূরণ করে তা নিশ্চিত করুন

পরীক্ষামূলক ফলাফল

প্রধান ফলাফল

উপপাদ্য 1.1 ফলাফল

নির্মিত উৎস ঘনত্ব ρ\rho এর জন্য (একক গোলক সীমানার দুটি বিন্দুতে অতি-বহুপদী বিস্ফোরণ), প্রমাণ করা হয়েছে যে যেকোনো গোলক Y=BRd(0,R)Y = B_{\mathbb{R}^d}(0,R), যেকোনো C,α>0C, \alpha > 0 এবং p1p \geq 1 এর জন্য, অসমতা: μ,νP(Y),TμTνL2(ρ)CWp(μ,ν)α\forall \mu, \nu \in P(Y), \|\|T_\mu - T_\nu\|\|_{L^2(\rho)} \leq C W_p(\mu, \nu)^\alpha বৈধ নয়।

উপপাদ্য 1.2 ফলাফল

সীমাবদ্ধ খোলা সেটে সমান সম্ভাব্যতা পরিমাপ ρ\rho এর জন্য, প্রমাণ করা হয়েছে যে যেকোনো C>0C > 0, p1p \geq 1 এবং α>p2(p+1)\alpha > \frac{p}{2(p+1)} (বিশেষত, α=12\alpha = \frac{1}{2}), স্থিতিশীলতা অসমতা ব্যর্থ হয়।

প্রযুক্তিগত আবিষ্কার

  1. ঘনত্ব বিস্ফোরণের সমালোচনামূলকতা: বহুপদী বিস্ফোরণ rd+δr^{-d+\delta} এখনও স্থিতিশীলতা বজায় রাখে, কিন্তু অতি-বহুপদী বিস্ফোরণ সম্পূর্ণ অস্থিতিশীলতার দিকে পরিচালিত করে
  2. Hölder নিয়মিততা সীমাবদ্ধতা: মানচিত্র μTμ\mu \mapsto T_\mu কিছু ক্ষেত্রে 12\frac{1}{2}-Hölder ক্রমাগত নয়, এমনকি 13\frac{1}{3}-Hölder ক্রমাগত নয়
  3. অ-অনন্যতার পূর্বাভাস: সর্বোত্তম পরিকল্পনার অনন্যতা হারানোর আগে স্থিতিশীলতা অবনতি ঘটে

সম্পর্কিত কাজ

স্থিতিশীলতা ইতিবাচক ফলাফল

  • Gigli (2011): মানচিত্র μTμ\mu \mapsto T_\mu এর 12\frac{1}{2}-Hölder ক্রমাগত তত্ত্বের উপরের সীমা
  • Letrouit & Mérigot (2024): John ডোমেনে সীমাবদ্ধ ঘনত্বের স্থিতিশীলতা অসমতা
  • অন্যান্য কাজ: লগ-অবতল ঘনত্ব, বহুপদী বিস্ফোরণ ঘনত্ব ইত্যাদির স্থিতিশীলতা ফলাফল

অস্থিতিশীলতা গবেষণা

এর আগে শুধুমাত্র Kantorovich সম্ভাব্যতা ফাংশনের অস্থিতিশীলতা উদাহরণ ছিল, এই কাগজটি প্রথমবার অস্থিতিশীল পরিবহন মানচিত্র তৈরি করেছে।

উপসংহার এবং আলোচনা

প্রধান উপসংহার

  1. সর্বোত্তম পরিবহন মানচিত্র অত্যন্ত অস্থিতিশীল হতে পারে, যেকোনো পরিমাণগত স্থিতিশীলতা অসমতা ব্যর্থ হতে পারে
  2. ঘনত্বের অতি-বহুপদী বিস্ফোরণ অস্থিতিশীলতার জন্য একটি যথেষ্ট শর্ত
  3. এমনকি সমান ঘনত্বের জন্যও, অ-অনন্য কনফিগারেশনের কাছাকাছি অস্থিতিশীলতা ঘটে

সীমাবদ্ধতা

  1. উপপাদ্য 1.1 সীমাহীন ঘনত্ব প্রয়োজন, বাস্তব প্রয়োগে অবাস্তব হতে পারে
  2. উপপাদ্য 1.2 এর নির্মাণ অসীম সংখ্যক সংযুক্ত উপাদান প্রয়োজন, জ্যামিতিক কাঠামো জটিল
  3. সীমাবদ্ধ ঘনত্ব উপপাদ্য 1.1 এর শক্তিশালী অস্থিতিশীলতা অর্জন করতে পারে কিনা তা নির্ধারণ করা হয়নি

ভবিষ্যত দিকনির্দেশনা

  1. সীমাবদ্ধ ঘনত্ব ক্ষেত্রে সর্বোচ্চ অস্থিতিশীলতা অধ্যয়ন করা
  2. অস্থিতিশীল উৎস পরিমাপের জন্য পরিসংখ্যানগত অনুমান পদ্ধতি বিকাশ করা
  3. অস্থিতিশীলতা এবং জ্যামিতিক কাঠামোর গভীর সংযোগ অন্বেষণ করা

গভীর মূল্যায়ন

সুবিধা

  1. তাত্ত্বিক যুগান্তকারী: প্রথমবার সর্বোত্তম পরিবহন মানচিত্রের অস্থিতিশীলতা প্রমাণ করা, তাত্ত্বিক শূন্যতা পূরণ করা
  2. নির্মাণ সূক্ষ্ম: জ্যামিতিক নির্মাণ এবং প্যারামিটার নির্বাচনের মাধ্যমে তাত্ত্বিক লক্ষ্য অর্জন করা
  3. বিশ্লেষণ কঠোর: গাণিতিক প্রমাণ সম্পূর্ণ, প্রযুক্তিগত বিবরণ যথাযথভাবে পরিচালিত
  4. অন্তর্দৃষ্টি গভীর: স্থিতিশীলতা এবং অনন্যতা, ঘনত্ব নিয়মিততার অভ্যন্তরীণ সংযোগ প্রকাশ করা

অপূর্ণতা

  1. ব্যবহারিক সীমাবদ্ধতা: নির্মিত উদাহরণ বাস্তব প্রয়োগে বিরল হতে পারে
  2. খোলা প্রশ্ন: কিছু প্রাকৃতিক প্রশ্ন (যেমন সীমাবদ্ধ ঘনত্বের শক্তিশালী অস্থিতিশীলতা) এখনও অমীমাংসিত
  3. সংখ্যাসূচক যাচাইকরণ অনুপস্থিত: তাত্ত্বিক কাজ হিসাবে, সংখ্যাসূচক পরীক্ষা যাচাইকরণ অনুপস্থিত

প্রভাব

  1. তাত্ত্বিক মূল্য: সর্বোত্তম পরিবহন তত্ত্বের জন্য গুরুত্বপূর্ণ প্রতিউদাহরণ এবং অন্তর্দৃষ্টি প্রদান করা
  2. প্রয়োগ সতর্কতা: অনুশীলনকারীদের কিছু কনফিগারেশনে অস্থিতিশীলতা ঝুঁকি সম্পর্কে সতর্ক করা
  3. গবেষণা অনুপ্রেরণা: পরবর্তী গবেষণার জন্য নতুন দিকনির্দেশনা এবং সরঞ্জাম প্রদান করা

প্রযোজ্য দৃশ্যকল্প

  1. সর্বোত্তম পরিবহন তত্ত্বের গভীর গবেষণা
  2. রৈখিক সর্বোত্তম পরিবহন কাঠামোর তাত্ত্বিক বিশ্লেষণ
  3. পরিসংখ্যানগত সর্বোত্তম পরিবহনে স্থিতিশীলতা বিশ্লেষণ
  4. সংখ্যাসূচক সর্বোত্তম পরিবহন অ্যালগরিদমের স্থিতিশীলতা মূল্যায়ন

তথ্যসূত্র

কাগজটি সর্বোত্তম পরিবহন ক্ষেত্রের গুরুত্বপূর্ণ সাহিত্য উদ্ধৃত করেছে, যার মধ্যে Brenier এর ভিত্তিমূলক কাজ, Gigli এর নিয়মিততা ফলাফল এবং সাম্প্রতিক পরিমাণগত স্থিতিশীলতা গবেষণা রয়েছে, যা গবেষণার জন্য একটি দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করে।