Unbalanced Optimal Transport and Density Control for Discrete-Time Linear Systems

तकनीकी सारांश: असंतुलित इष्टतम परिवहन और डिस्क्रीट-टाइम रैखिक प्रणालियों के लिए घनत्व नियंत्रण

समस्या विवरण
यह कार्य शास्त्रीय इष्टतम परिवहन (Optimal Transport - OT) की सीमाओं को संबोधित करता है, जिसके लिए वितरणों का कुल द्रव्यमान समान होना आवश्यक है, जो इसे द्रव्यमान निर्माण, क्षय, या लुप्त डेटा वाले परिदृश्यों के लिए अनुपयुक्त बनाता है। लेखक असंतुलित इष्टतम परिवहन (Unbalanced Optimal Transport - UOT) की जांच करते हैं, जो संदर्भ मापों (reference measures) से विचलन को कुलबैक–लीब्लर (KL) डाइवर्जेंस का उपयोग करके दंडित करके द्रव्यमान संरक्षण बाधा को शिथिल करता है। इसके अलावा, यह शोध पत्र डिस्क्रीट-टाइम रैखिक प्रणालियों के लिए असंतुलित घनत्व नियंत्रण (Unbalanced Density Control - UDC) हेतु इस ढांचे का विस्तार करता है। UDC में, लक्ष्य एक प्रारंभिक माप से एक लक्षित माप तक सिस्टम के अवस्था वितरण (state distribution) को निर्देशित करना है, जबकि नियंत्रण प्रयास को कम करना और संदर्भ वितरणों के प्रति निष्ठा के विरुद्ध परिवहन लागत को संतुलित करना है, और यह बिना कुल संभाव्यता द्रव्यमान के संरक्षण की धारणा के किया जाता है।

कार्यप्रणाली
यह शोध पत्र उस मामले पर केंद्रित है जहाँ संदर्भ माप गाऊसी (Gaussian) हैं और सिस्टम डायनेमिक्स रैखिक हैं। मुख्य पद्धतिगत योगदान दोनों UOT और UDC के लिए वैश्विक रूप से इष्टतम, परिमित-आयामी उत्तलता (convex) सूत्रीकरण का व्युत्पन्न है।

1. गाऊसी मापों तक प्रतिबंध:
   लेखक सिद्ध करते हैं कि UOT और UDC दोनों के लिए, एक ऐसा इष्टतम समाधान मौजूद है जहाँ मार्जिनल वितरण (और नियंत्रण के मामले में अवस्था माप) गाऊसी होते हैं।

• UOT के लिए: वे प्रदर्शित करते हैं कि निश्चित द्रव्यमान, माध्य और सहप्रसरण (covariance) के लिए, परिवहन लागत तब न्यूनतम होती है जब मार्जिनल गाऊसी होते हैं (प्रपोजिशन 1), और KL डाइवर्जेंस पद भी समान बाधाओं के तहत गाऊसी मापों द्वारा न्यूनतम किया जाता है (प्रपोजिशन 2)। यह अनंत-आयामी परिवर्तनशील समस्या को मार्जिनल्स के द्रव्यमान ($c$), माध्य ($m$), और सहप्रसरण ($\Sigma$) पर अनुकूलन (optimization) में बदलने की अनुमति देता है।
• UDC के लिए: इसी तरह, वे दिखाते हैं कि किसी भी व्यवहार्य प्रक्रिया के लिए, एक गाऊसी अवस्था माप और एक गाऊसी-एफाइन फीडबैक नीति बनाई जा सकती है जो मूल प्रक्रिया के प्रथम और द्वितीय क्षणों (moments) से मेल खाती है और एक निम्न या समान ऑब्जेक्टिव वैल्यू प्रदान करती है (थ्योरम 5)।

2. अनुकूलन का पृथक्करण (Separation of Optimization):
   परिणामी परिमित-आयामी समस्याएं सभी चरों (विशेष रूप से द्रव्यमान और क्षणों) में संयुक्त रूप से उत्तलीय (jointly convex) नहीं हैं। हालांकि, लेखक इस समस्या की संरचना का लाभ उठाते हुए अनुकूलन को अलग करते हैं:

• UOT सूत्रीकरण: समस्या द्रव्यमान और क्षणों के फलन को न्यूनतम करने तक सीमित हो जाती है। लेखक एक एल्गोरिदम (एल्गोरिदम 1) प्रस्तावित करते हैं जो पहले माध्य और सहप्रसरणों के लिए एक उत्तलीय उप-समस्या को हल करता है (एक इष्टतम मान $p^*$ प्राप्त करता है) और फिर $p^*$ के आधार पर बंद रूप (closed form) में इष्टतम द्रव्यमान $c^*$ की गणना करता है।
• UDC सूत्रीकरण: नियंत्रण समस्या को द्रव्यमान, माध्य प्रक्षेपवक्र (trajectory), सहप्रसरण प्रक्षेपवक्र और फीडबैक मापदंडों पर अनुकूलन के रूप में तैयार किया गया है। बाल्सी और बाकोलास (2022) के मानक कोवेरियेंस स्टीयरिंग पुनर्रूपण (reparametrization) को लागू करके, गैर-उत्तलीय सहप्रसरण बाधाओं को शूर कॉम्प्लीमेंट (Schur complement) के माध्यम से एफाइन बाधाओं और लीनियर मैट्रिक्स इनइक्वेलिटीज (LMIs) में परिवर्तित किया जाता है। यह उप-समस्या को एक सेमिडेफिनेट प्रोग्रामिंग (SDP) में बदल देता है। UOT की तरह, प्रक्षेप पथों के लिए SDP को हल करने के बाद द्रव्यमान को एक क्लोज्ड-फॉर्म स्टेप में अपडेट किया जाता है।

मुख्य योगदान

• गाऊसी UOT: गाऊसी मापों के बीच अनबैलेंस्ड ऑप्टिमल ट्रांसपोर्ट के लिए एक परिमित-आयामी उत्तलीय सूत्रीकरण और एक वैश्विक रूप से इष्टतम एल्गोरिदम का व्युत्पन्न।
• गाऊसी UDC: डिस्क्रीट-टाइम रैखिक प्रणालियों के लिए अनबैलेंस्ड डेंसिटी कंट्रोल का सूत्रीकरण और एक वैश्विक रूप से इष्टतम, SDP-आधारित एल्गोरिदम का प्रस्ताव।
• सैद्धांतिक न्यूनीकरण (Theoretical Reduction): यह प्रमाण कि खोज स्थान को गाऊसी मापों और एफाइन कपलिंग/नीतियों तक सीमित करने से इन विशिष्ट समस्या वर्गों के लिए वैश्विक इष्टतमता से समझौता नहीं होता है।

परिणाम

• एल्गोरिद्मिक दक्षता: प्रस्तावित विधियाँ जटिल परिवर्तनशील समस्याओं को सुलभ उत्तलीय अनुकूलन समस्याओं (विशेष रूप से नियंत्रण मामले के लिए SDP) में परिवर्तित करती हैं जिन्हें वैश्विक रूप से हल किया जा सकता है।
• संख्यात्मक चित्रण: एक एक-आयामी संख्यात्मक प्रयोग UOT समाधान के व्यवहार को प्रदर्शित करता है। परिणाम दिखाते हैं कि जैसे-जैसे दंड पैरामीटर $\gamma$ बढ़ता है, समाधान मानक संतुलित OT योजना की ओर अभिसरित (converge) होता है, और इष्टतम मार्जिनल्स संदर्भ मापों के साथ अधिक निकटता से संरेखित होते हैं। इसके विपरीत, छोटे $\gamma$ के लिए, समाधान परिवहन लागत को कम करने (विकर्ण के पास केंद्रित होने) को प्राथमिकता देता है, जिससे प्रभावी रूप से परिवहन दूरी को कम करने के लिए द्रव्यमान का निर्माण या विनाश संभव हो जाता है।

महत्व और दावे
यह शोध पत्र दावा करता है कि कोवेरियेंस स्टीयरिंग को अनबैलेंस्ड सेटिंग में विस्तारित करके, यह उन प्रणालियों में संभाव्यता घनत्वों को नियंत्रित करने के लिए एक कठोर ढांचा प्रदान करता है जहाँ द्रव्यमान संरक्षित नहीं होता है। इसका महत्व डेटा हानि या क्षय के कारण होने वाले "मिसमैच्ड मास" को एक नियंत्रण-सैद्धांतिक ढांचे के भीतर संभालने की क्षमता में निहित है। लेखक इस कार्य को एक आधारभूत कदम के रूप में प्रस्तुत करते हैं, यह उल्लेख करते हुए कि पूर्ण प्रमाण और आगे के प्रयोगों के साथ एक काफी विस्तारित संस्करण एक अलग प्रकाशन (नाकाशिमा एट अल., 2026) में उपलब्ध है। वे स्पष्ट रूप से कहते हैं कि भविष्य के कार्य गैर-गाऊसी संदर्भ मापों, वैकल्पिक डाइवर्जेंस विकल्पों और वास्तविक दुनिया की गतिशील प्रणालियों के अनुप्रयोगों को संबोधित करेंगे, जो यह संकेत देते हैं कि वर्तमान कार्य गाऊसी रैखिक मामले के लिए एक सैद्धांतिक और एल्गोरिद्मिक प्रूफ-ऑफ-कॉन्सेप्ट है।