Characterization of stability radii for robustly asymptotically stable dissipative Hamiltonian differential-algebraic systems

तकनीकी सारांश: सुदृढ़ रूप से स्पर्शोन्मुखी रूप से स्थिर विपाश्विक हैमिल्टोनियन विभेदक-बीजगणितीय प्रणालियों के लिए स्थिरता त्रिज्याओं का अभिलक्षणन

समस्या विवरण
यह शोध पत्र रैखिक समय-अपरिवर्तनीय विपाश्विक हैमिल्टोनियन विभेदक-बीजगणितीय समीकरण (dHDAE) प्रणालियों के रूप में $E \dot{x} = (J - R)Qx$ के सुदृढ़ रूप से स्पर्शोन्मुखी स्थिरता (robust asymptotic stability) के विषय को संबोधित करता है। जबकि यह स्थापित है कि ऐसी प्रणालियों के परिमित आइगेन मान (eigenvalues) बंद बाएं अर्ध-तल (closed left half-plane) में स्थित होते हैं और विशुद्ध काल्पनिक आइगेन मान (शून्य को छोड़कर) अर्ध-सरल (semisimple) होते हैं, ये प्रणालियाँ अनिवार्य रूप से सुदृढ़ रूप से स्पर्शोन्मुखी स्थिर नहीं होती हैं। गुणांक मैट्रिसेस (coefficient matrices) में छोटे विचलन नियमितता (regularity) को नष्ट कर सकते हैं, विभेदक-बीजगणितीय सूचकांक (index) को बढ़ा सकते हैं, या आइगेन मानों को काल्पनिक अक्ष की ओर ले जा सकते हैं, जिससे अस्थिरता उत्पन्न हो सकती है।

मुख्य समस्या यह मात्रा निर्धारित करना है कि एक दी गई सुदृढ़ रूप से स्पर्शोन्मुखी स्थिर dHDAE प्रणाली और ऐसी प्रणालियों के समुच्चय की सीमा के बीच की "दूरी" कितनी है। विशेष रूप से, लेखक उन न्यूनतम संरचित विचलनों (perturbations) का लक्षण वर्णन करने का प्रयास करते हैं जो निम्नलिखित कारणों से उत्पन्न होते हैं:

1. विलक्षणता (singularity) (नियमितता की हानि)।
2. विभेदक-बीजगणितीय सूचकांक का एक से अधिक होने तक बढ़ना।
3. विशुद्ध काल्पनिक आइगेन मानों का प्रकट होना (स्पर्शोन्मुखी स्थिरता की हानि)।

विश्लेषण दो प्रकार के संरचित विचलनों पर विचार करता है:

• $S_d$ (विपाश्विक/Dissipative): विचलन जो $E$ और $R$ की चिह्न-निश्चितता (sign-definiteness) को बनाए रखते हैं (अर्थात $\Delta E \le 0, \Delta R \le 0$) और $J$ की विषम-हर्मिटियन (skew-Hermitian) संरचना को बनाए रखते हैं।
• $S_i$ (अनिश्चित/Indefinite): विचलन जो हर्मिटियन/विषम-हर्मिटियन संरचनाओं को बनाए रखते हैं लेकिन विचलनों की चिह्न-निश्चितता बाधाओं को अनिवार्य रूप से सुरक्षित नहीं रखते हैं, केवल यह सुनिश्चित करते हैं कि परिणामी मैट्रिसेस धनात्मक अर्ध-निश्चित (positive semidefinite) बने रहें।

कार्यप्रणाली
लेखक संरचनात्मक मैट्रिक्स सिद्धांत, स्पेक्ट्रल विश्लेषण और अनुकूलन तकनीकों के संयोजन का उपयोग करते हैं:

1. संरचनात्मक अभिलक्षणन: शोध पत्र क्रोनेकर कैनोनिकल रूप (Kronecker canonical form) और एक विशिष्ट "सीढ़ीनुमा रूप" (staircase form) (लेम्मा 8) का उपयोग करता है, जिसे विभेदक-बीजगणितीय पेंसिल के विश्लेषण के लिए यूनिटरी रूपांतरणों के माध्यम से प्राप्त किया गया है। सुदृढ़ स्पर्शोन्मुखी स्थिरता को तीन शर्तों द्वारा अभिलक्षित किया जाता है: नियमितता, सूचकांक $\le 1$, और एक ऋणात्मक स्पेक्ट्रल एब्सिसा (spectral abscissa) (प्रमेय 11)।
2. मानक रूपों में न्यूनीकरण: विश्लेषण अक्सर यह मान लेता है कि $Q$ व्युत्क्रमणीय (invertible) है, जिससे प्रणाली को $E \dot{x} = (J - R)x$ के रूप में रूपांतरित किया जा सकता है जहाँ $E=E^* \ge 0$, $R=R^* \ge 0$, और $J=-J^*$ है।
3. मैपिंग समस्याएँ: स्थिरता त्रिज्या की गणना करने के लिए, लेखक एक न्यूनतम मानदंड विचलन $(\Delta E, \Delta J, \Delta R)$ खोजने के रूप में समस्या को तैयार करते हैं जो विशिष्ट बीजगणितीय बाधाओं (जैसे, किसी मैट्रिक्स को विलक्षण बनाना या काल्पनिक अक्ष पर आइगेन मान बनाना) को संतुष्ट करता है। वे हर्मिटियन और विषम-हर्मिटियन मैपिंग समस्याओं के न्यूनतम मानदंड समाधानों (प्रमेय 15 और 16) का लाभ उठाते हैं।
4. अनुकूलन सूत्र (Optimization Formulations): स्थिरता दूरियों को सामान्यीकृत रेले कोशेंट (Rayleigh quotients) के माध्यम से अनुकूलन समस्याओं के रूप में व्यक्त किया गया है। उन मामलों के लिए जहाँ $E$ स्थिर है (आंशिक विचलन), सटीक अभिलक्षण प्राप्त किए गए हैं। पूर्ण विचलनों (जहाँ $E$ परिवर्तित होता है) के लिए, लेखक सीमाएँ (bounds) प्राप्त करते हैं।
   • विलक्षणता और उच्च सूचकांक की दूरी में रैंक स्थितियों के अधीन वर्ग मानदंडों के योग को न्यूनतम करना शामिल है।
   • काल्पनिक आइगेन मानों की दूरी को आवृत्ति पैरामीटर $\omega$ और एक वेक्टर $x$ पर न्यूनतम करने के रूप में दिया गया है।
5. संख्यात्मक कार्यान्वयन: प्राप्त अनुकूलन समस्याओं को दो-स्तरीय रणनीति का उपयोग करके हल किया जाता है: एक आंतरिक लूप जो लेवल-शिफ्टेड सेल्फ-कंसिस्टेंट फील्ड (SCF) पुनरावृत्तियों के माध्यम से एक नॉनलीनियर आइगेन मान समस्या (NEPv) को हल करता है, और एक बाहरी लूप जो मानक अनुकूलन विधियों (जैसे fminsearch, fmincon) का उपयोग करके आवृत्ति पैरामीटर $\omega$ पर न्यूनतमकरण करता है।

प्रमुख योगदान और परिणाम

• आंशिक विचलनों के लिए सटीक अभिलक्षणन: जब मैट्रिक्स $E$ स्थिर रहता है और केवल $J$ और $R$ को विचलित किया जाता है, तो यह शोध पत्र विलक्षणता, उच्च सूचकांक और काल्पनिक आइगेन मानों के लिए संरचित दूरियों के सटीक सूत्र प्रदान करता है (प्रमेय 24–26)।
  • $S_i$ विचलनों के लिए, जिसमें $R > 0$ है, विलक्षणता की दूरी ठीक $\sqrt{\lambda_{\min}(N^*(R^2 - J^2)N)}$ है, जहाँ $N$, $E$ के नुल स्पेस (null space) को दर्शाता है।
  • $S_i$ के लिए काल्पनिक आइगेन मानों की दूरी $\inf_{\omega} \sigma_{\min}([R, i\omega E - J])$ द्वारा दी गई है।
• पूर्ण विचलनों के लिए सीमाएँ: जब सभी मैट्रिसेस ($E, J, R$) में विचलन की अनुमति होती है, तो लेखक स्पष्ट निचली और ऊपरी सीमाएँ प्राप्त करते हैं:
  • विलक्षणता की दूरी: $S_d$ के लिए सटीक सूत्र (प्रमेय 20) और $S_i$ के लिए सीमाएँ (प्रमेय 21) प्रदान की गई हैं।
  • उच्च सूचकांक की दूरी: एक दो-चरणीय विचलन दृष्टिकोण (पहले $E$ को नुल स्पेस के आयाम को बढ़ाने के लिए विचलित करना, फिर $J$ और $R$ को विचलित करना) का उपयोग करके ऊपरी सीमाएँ प्राप्त की गई हैं (प्रमेय 22)। लेखक नोट करते हैं कि द्वितीय-क्रम के विचलन पदों के कारण, ये वर्तमान में केवल ऊपरी सीमाएँ हैं।
  • काल्पनिक आइगेन मानों की दूरी: $S_d$ और $S_i$ दोनों के लिए निचली सीमाएँ प्राप्त की गई हैं (प्रमेय 18)। इनमें विशिष्ट मैट्रिक्स पेंसिल के सामान्यीकृत रेले कोशेंट और आइगेन मानों पर जटिल अनुकूलन शामिल है।
• एल्गोरिदम ढांचा: मानक रैखिक बीजगणिक अपघटन (SVD, स्पेक्ट्रल अपघटन) का उपयोग करके सुदृढ़ स्पर्शोन्मुखी स्थिरता की जाँच करने और संबंधित दूरी सीमाओं की गणना करने के लिए एक चरण-दर-चरण एल्गोरिदम (अनुभाग 2.7) प्रस्तुत किया गया है।

महत्व और दावे
यह शोध पत्र दावा करता है कि यह संरचना-संरक्षण विचलनों के तहत dHDAE प्रणालियों के लिए स्थिरता त्रिज्याओं का पहला व्यापक अभिलक्षणन प्रदान करता है। जबकि सामान्य असंरचित DAEs के लिए स्थिरता त्रिज्या की गणना करना कठिन ज्ञात है, लेखक प्रदर्शित करते हैं कि विपाश्विक हैमिल्टोनियन प्रणालियों के विशिष्ट ज्यामितीय और बीजगणितीय गुण स्पष्ट सूत्रों और गणना योग्य सीमाओं के व्युत्पन्न को सक्षम करते हैं।

इसका महत्व भौतिक मॉडलों (जैसे संरचनात्मक यांत्रिकी, द्रव गति विज्ञान और सरंध्र माध्यम प्रवाह) की मॉडलिंग और विविक्तकरण (discretization) त्रुटियों के प्रति सहनशीलता को मापने की क्षमता में निहित है। निकटतम विलक्षण, उच्च-सूचकांक, या अस्थिर प्रणाली की दूरी निर्धारित करके, इंजीनियर और वैज्ञानिक अपने dHDAE मॉडलों की विश्वसनीयता का आकलन कर सकते हैं।

लेखक पूर्ण विचलनों (विशेष रूप से उच्च सूचकांक की दूरी के लिए) के लिए सीमाओं की सटीकता के संबंध में विनम्र हैं, और स्वीकार करते हैं कि जटिल द्वितीय-क्रम के विचलन पदों के कारण वर्तमान परिणाम केवल ऊपरी सीमाएँ हैं। वे भविष्य के अनुसंधान के लिए आवश्यक दिशाओं के रूप में सटीक सूत्रों के निर्धारण और अधिक कुशल संख्यात्मक अनुकूलन विधियों के विकास की पहचान करते हैं।