2025-11-12T23:07:10.120463

Symplectic model order reduction of port-Hamiltonian systems

Glas, Mamunuzzaman, Mu et al.

This work proposes a novel structure-preserving model order reduction (MOR) method for linear, time-invariant port-Hamiltonian (pH) systems. Our goal is to construct a reduced order pH system, which can still be interpreted in the physical domain of the full order model. By this we mean, that if an electrical circuit is the initial high-dimensional pH system, we want the reduced order model to be still interpretable as an electronic circuit. In the case of the well-known mass spring damper (MSD) system, there are MOR methods available, which already guarantee the preservation of this particular structure. Moreover, we show that our new structure-preserving MOR method, which is based on symplectic MOR methods, will recover the known second-order Arnoldi method in the case of MSD systems. However, for the example of an electrical circuit pH model (and more models of similar block structure), our method yields a novel model reduction method. We present numerical results on the aforementioned electronic circuit model, highlighting the advantages of the proposed method.

academic

पोर्ट-हैमिल्टनियन सिस्टम्स का सिम्प्लेक्टिक मॉडल ऑर्डर रिडक्शन

मूल जानकारी

पेपर ID: 2203.07751
शीर्षक: Symplectic model order reduction of port-Hamiltonian systems
लेखक: Silke Glas, Mir Mamunuzzaman, Hongliang Mu, Hans Zwart
वर्गीकरण: math.OC (अनुकूलन और नियंत्रण), math.CA (शास्त्रीय विश्लेषण)
प्रकाशन समय: मार्च 2022 (arXiv प्रीप्रिंट)
पेपर लिंक: https://arxiv.org/abs/2203.07751

सारांश

यह पेपर रैखिक समय-अपरिवर्तनीय पोर्ट-हैमिल्टनियन (pH) सिस्टम्स के लिए एक नवीन संरचना-संरक्षक मॉडल क्रम में कमी (MOR) विधि प्रस्तुत करता है। अनुसंधान का लक्ष्य एक कम क्रम वाली pH प्रणाली का निर्माण करना है जो मूल पूर्ण-क्रम मॉडल के भौतिक डोमेन में व्याख्या योग्य रहे। विशेष रूप से, यदि प्रारंभिक उच्च-आयामी pH प्रणाली एक विद्युत परिपथ है, तो हम चाहते हैं कि कम क्रम वाला मॉडल अभी भी एक इलेक्ट्रॉनिक परिपथ के रूप में व्याख्या योग्य हो। द्रव्यमान-वसंत-अवमंदन (MSD) सिस्टम्स के लिए, MOR विधियाँ पहले से ही इस विशिष्ट संरचना के संरक्षण की गारंटी देती हैं। यह पेपर दर्शाता है कि सिम्प्लेक्टिक MOR विधि पर आधारित नई संरचना-संरक्षक MOR विधि MSD प्रणाली के मामले में ज्ञात द्वितीय-क्रम Arnoldi विधि को पुनः प्राप्त कर सकती है। हालांकि, विद्युत परिपथ pH मॉडल्स जैसे समान ब्लॉक संरचना वाले मॉडल्स के लिए, यह विधि नई मॉडल क्रम में कमी तकनीकें उत्पन्न करती है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या विवरण

पोर्ट-हैमिल्टनियन सिस्टम्स शास्त्रीय हैमिल्टनियन ढांचे का सामान्यीकरण हैं, जो इंजीनियरिंग समस्याओं में अत्यंत उपयोगी हैं क्योंकि वे शास्त्रीय हैमिल्टनियन विधि को खुली भौतिक प्रणालियों के मॉडलिंग के शक्तिशाली ढांचे के साथ जोड़ते हैं। हालांकि, नियंत्रण डिजाइन में, नियंत्रक आमतौर पर मूल प्रणाली के समान क्रम का होता है, जो मूल प्रणाली के मौलिक व्यवहार को पकड़ने वाले कम क्रम वाले मॉडल की आवश्यकता को प्रेरित करता है।

मौजूदा विधियों की सीमाएं

संरचना हानि समस्या: पारंपरिक MOR विधियों को सीधे लागू करने से pH संरचना खो जाती है
भौतिक व्याख्या की कमी: मौजूदा विधियाँ भौतिक डोमेन में कम क्रम वाले मॉडल की व्याख्या योग्यता की गारंटी नहीं दे सकती हैं
ब्लॉक संरचना विनाश: विशेष ब्लॉक संरचना वाली pH प्रणालियों के लिए, मौजूदा विधियाँ इस संरचना को बनाए नहीं रख सकती हैं

अनुसंधान प्रेरणा

जब भौतिक मॉडल्स को pH रूप में पुनः तैयार किया जाता है, तो मैट्रिक्स J, R और H आमतौर पर विशेष ब्लॉक संरचना रखते हैं। इस पेपर का लक्ष्य क्रम में कमी की प्रक्रिया के दौरान इस ब्लॉक संरचना को बनाए रखना है, यह सुनिश्चित करते हुए कि कम क्रम वाला मॉडल अभी भी पूर्ण-क्रम मॉडल के समान प्रकार की भौतिक प्रणाली के रूप में व्याख्या योग्य है।

मुख्य योगदान

सिम्प्लेक्टिक-pH विधि का परिचय: सिम्प्लेक्टिक MOR पर आधारित pH सिस्टम्स के लिए एक नवीन MOR विधि प्रस्तावित की गई है, जो ब्लॉक संरचना रूप (2) के लिए उपयुक्त है, जिससे कम क्रम वाला मॉडल समान रूप की pH प्रणाली बना रहता है
सैद्धांतिक संबंध: MSD सिस्टम्स के विशेष मामले में, यह साबित किया गया है कि प्रस्तावित विधि द्वितीय-क्रम Arnoldi विधि (SO-Arnoldi) को पुनः प्राप्त करती है
संख्यात्मक सत्यापन: इलेक्ट्रॉनिक परिपथ (LRCR परिपथ) पर संख्यात्मक परिणाम प्रदर्शित किए गए हैं, साहित्य में अन्य ज्ञात विधियों के साथ तुलना की गई है
स्थिरता लाभ: पाया गया है कि SO-Arnoldi विधि कुछ मामलों में अस्थिर कम क्रम वाले मॉडल्स उत्पन्न करती है, जबकि यह विधि स्थिरता बनाए रखती है

विधि विवरण

कार्य परिभाषा

दी गई रैखिक समय-अपरिवर्तनीय pH प्रणाली:

ẋ(t) = (J - R)Hx(t) + Bu(t)
y(t) = B*Hx(t)

जहाँ x(t) ∈ C²ⁿ स्थिति है, u(t) ∈ Cᵖ इनपुट है, y(t) ∈ Cᵖ आउटपुट है। लक्ष्य समान रूप की लेकिन छोटे स्थिति-स्पेस आयाम वाली प्रणाली खोजना है, जिससे दोनों ट्रांसफर फंक्शन पूर्वनिर्धारित आवृत्ति बिंदुओं पर समान हों।

मुख्य सैद्धांतिक ढांचा

सिम्प्लेक्टिक ज्यामिति आधार

परिभाषा 1 (सिम्प्लेक्टिक रूप): मान लीजिए V, C पर एक सम-आयामी सदिश स्पेस है, सिम्प्लेक्टिक रूप Ω: V × V → C निम्नलिखित को संतुष्ट करता है:

विषम-सममितता: Ω(v₁, v₂) = -Ω(v₂, v₁)
अ-विकृति: सभी v₁ ∈ V के लिए Ω(v₁, v₂) = 0 यदि और केवल यदि v₂ = 0

परिभाषा 2 (सिम्प्लेक्टिक मानचित्र): मैट्रिक्स Q को (J₂ₙ, J₂ₖ)-सिम्प्लेक्टिक कहा जाता है, यदि:

Q*J₂ₙQ = J₂ₖ

अ-अपव्यय स्थिति में सिम्प्लेक्टिक MOR

प्रमेय 4: अ-अपव्यय pH प्रणाली (R = 0) के लिए, यदि Q ∈ C²ⁿˣ²ᵏ एक (J₂ₙ⁻¹, J₂ₖ⁻¹)-सिम्प्लेक्टिक मैट्रिक्स है, और निम्नलिखित को संतुष्ट करता है:

ran(B) ⊆ ran(Q)
(sₘI₂ₙ - J₂ₙH)⁻¹Bμₘ ∈ ran(Q)

तो कम क्रम वाली pH प्रणाली दिए गए आवृत्ति पर मूल प्रणाली के बराबर है।

अपव्यय स्थिति में सिम्प्लेक्टिक MOR

प्रमेय 5: अपव्यय वाली pH प्रणाली के लिए, विस्तारित मैट्रिक्स का निर्माण करें:

QE := [Q₁  0  ]
      [0   Q₂]

जहाँ Q₂ (J₂ₙ, J₂ₖ)-सिम्प्लेक्टिक है, Q₁* Q₂ का बायां प्रतिलोम है। उपयुक्त शर्तों के तहत, कम क्रम वाली प्रणाली अंतर्वेशन शर्तों को संतुष्ट करती है।

ब्लॉक संरचना-संरक्षक एल्गोरिदम

एल्गोरिदम 1 (सिम्प्लेक्टिक-pH MOR):

अंतर्वेशन डेटा के अनुसार V₁ का निर्माण करें
V₂ := HV₁ की गणना करें
V₁ और V₂ को ब्लॉक संरचना के अनुसार विभाजित करें
सूत्र (27) के माध्यम से Ψ₁, Ψ₂, Φ₁, Φ₂ का निर्माण करें
Q₁ और Q₂ का निर्माण करें
Petrov-Galerkin प्रक्षेपण लागू करें

निर्माण सूत्र:

Ψ₂ = orth(V₂₂)
Ψ₁ = V₁₂(Ψ₂*V₁₂)⁻¹
Φ₂ = V₂₁(Ψ₂*Ĵₙ*V₂₁)⁻¹Ĵₖ*
Φ₁ = V₁₁(Φ₂*V₁₁)⁻¹

तकनीकी नवाचार बिंदु

सिम्प्लेक्टिक संरचना संरक्षण: सिम्प्लेक्टिक मानचित्र के माध्यम से हैमिल्टनियन प्रणाली की मौलिक ज्यामितीय संरचना को बनाए रखना
ब्लॉक विकर्ण संरचना रखरखाव: Q₁ और Q₂ का विशेष निर्माण यह सुनिश्चित करता है कि भौतिक अर्थ की ब्लॉक संरचना बनी रहे
अंतर्वेशन शर्तें: कम क्रम वाली प्रणाली को पूर्वनिर्धारित आवृत्ति बिंदुओं पर मूल प्रणाली से सटीक रूप से मेल खाने की गारंटी देना
एकीकृत ढांचा: MSD सिस्टम्स की SO-Arnoldi विधि को एक विशेष मामले के रूप में एक ढांचे में शामिल करना

प्रायोगिक सेटअप

संख्यात्मक उदाहरण: LRCR परिपथ

चित्र 1 में दिखाए गए LRCR परिपथ पर विचार करें, जिसमें n=50 इकाइयाँ हैं, पैरामीटर सेटिंग:

संधारित्र: Cᵢ = 1 pF
प्रेरकत्व: Lᵢ = 1 nH
प्रतिरोध: Rᵢ,₁ = 1 mΩ, Rᵢ,₂ = 1 kΩ

मूल्यांकन मेट्रिक्स

H₂ त्रुटि:

H₂error = √(1/2π ∫ Trace((G(iω) - Gred(iω))*(G(iω) - Gred(iω)))dω)

H∞ त्रुटि:

H∞error = max‖G(iω) - Gred(iω)‖₂

तुलना विधियाँ

pH-Arnoldi: प्रमेय 2 पर आधारित संरचना-संरक्षक विधि
SO-Arnoldi: द्वितीय-क्रम Arnoldi विधि
सिम्प्लेक्टिक-pH: इस पेपर द्वारा प्रस्तावित विधि

अंतर्वेशन बिंदु चयन

अंतराल 10⁻², 10² में लॉगरिदमिक पैमाने पर समान रूप से अंतर्वेशन बिंदु चुनें, उदाहरण के लिए 6 अंतर्वेशन बिंदु: 10⁻²i, -10⁻²i, 10⁰i, -10⁰i, 10²i, -10²i

प्रायोगिक परिणाम

मुख्य परिणाम

चित्र 3 से निम्नलिखित देखा जा सकता है:

सिम्प्लेक्टिक-pH बनाम pH-Arnoldi: सिम्प्लेक्टिक-pH अधिकांश मामलों में कम H₂ और H∞ त्रुटि प्रदान करता है
सिम्प्लेक्टिक-pH बनाम SO-Arnoldi:
- कम क्रम आयाम <24 होने पर प्रदर्शन समान है
- कम क्रम आयाम =32 होने पर SO-Arnoldi का प्रदर्शन तेजी से गिरता है
- कम क्रम आयाम 28 और 36 होने पर SO-Arnoldi अस्थिर प्रणाली उत्पन्न करता है

स्थिरता विश्लेषण

चित्र 4 SO-Arnoldi विधि द्वारा r=28 और r=36 पर उत्पन्न अस्थिर ध्रुवों को दर्शाता है, जबकि सिम्प्लेक्टिक-pH विधि सभी परीक्षित आयामों पर स्थिरता बनाए रखती है।

सैद्धांतिक सत्यापन: MSD प्रणाली समतुल्यता

प्रमेय 7 साबित करता है कि MSD सिस्टम्स के लिए, सिम्प्लेक्टिक-pH विधि एकात्मक रूपांतर के माध्यम से SO-Arnoldi विधि के समतुल्य है: एकात्मक मैट्रिक्स Θ मौजूद है जिससे zred(t) = Θẑred(t)

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

pH संरचना और ब्लॉक संरचना को संरक्षित करने वाली MOR विधि का सफल विकास
MSD सिस्टम्स की विभिन्न MOR विधियों का सैद्धांतिक एकीकरण
विधि की प्रभावशीलता और स्थिरता लाभ का संख्यात्मक सत्यापन

सीमाएं

अंतर्वेशन विधि: वर्तमान में केवल ट्रांसफर फंक्शन मान मिलाता है, व्युत्पन्न मिलान पर विचार नहीं करता है
ब्लॉक संरचना प्रतिबंध: विधि विशेष ब्लॉक संरचना के लिए डिज़ाइन की गई है
कम्प्यूटेशनल जटिलता: सिम्प्लेक्टिक मैट्रिक्स निर्माण की कम्प्यूटेशनल लागत का आगे विश्लेषण आवश्यक है

भविष्य की दिशाएं

व्युत्पन्न मिलान: ट्रांसफर फंक्शन के व्युत्पन्न मिलान तक विस्तार करना ताकि फिटिंग गुणवत्ता में सुधार हो
अनुकूली चयन: अंतर्वेशन बिंदुओं के अनुकूली चयन रणनीति विकसित करना
अरैखिक विस्तार: अरैखिक pH सिस्टम्स तक विस्तार
बड़े पैमाने पर अनुप्रयोग: बड़े पैमाने की प्रणालियों को संभालने के लिए एल्गोरिदम को अनुकूलित करना

गहन मूल्यांकन

लाभ

सैद्धांतिक कठोरता: सिम्प्लेक्टिक ज्यामिति पर आधारित ठोस गणितीय आधार
संरचना संरक्षण: pH संरचना और भौतिक अर्थ की ब्लॉक संरचना दोनों को संरक्षित करना
एकीकृतता: मौजूदा विधियों को एक ढांचे में विशेष मामलों के रूप में एकीकृत करना
व्यावहारिकता: विधि में स्पष्ट एल्गोरिदम कार्यान्वयन और संख्यात्मक सत्यापन है

कमियाँ

प्रयोज्यता सीमा: विधि विशेष ब्लॉक संरचना के लिए है, सामान्यता सीमित है
कम्प्यूटेशनल दक्षता: सिम्प्लेक्टिक मैट्रिक्स निर्माण की कम्प्यूटेशनल जटिलता विश्लेषण अपर्याप्त है
बड़े पैमाने पर सत्यापन: बड़ी प्रणालियों के संख्यात्मक सत्यापन की कमी

प्रभाव

सैद्धांतिक योगदान: pH प्रणाली MOR के लिए नया सैद्धांतिक दृष्टिकोण प्रदान करता है
व्यावहारिक मूल्य: विद्युत परिपथ विश्लेषण जैसे इंजीनियरिंग अनुप्रयोगों में सीधा मूल्य है
पद्धति महत्व: MOR में ज्यामितीय विधियों की महत्वपूर्ण भूमिका को प्रदर्शित करता है

प्रयोज्य परिदृश्य

विद्युत परिपथ विश्लेषण: बड़े पैमाने के इलेक्ट्रॉनिक परिपथों का क्रम में कमी मॉडलिंग
यांत्रिक प्रणाली: द्रव्यमान-वसंत-अवमंदन प्रणालियों का नियंत्रण डिजाइन
बहु-भौतिकी प्रणाली: ऊर्जा संरक्षण विशेषताओं वाली जटिल भौतिक प्रणालियाँ

संदर्भ

पेपर में 31 संदर्भ हैं, जो MOR सिद्धांत, pH प्रणाली मॉडलिंग, सिम्प्लेक्टिक ज्यामिति और नियंत्रण सिद्धांत सहित कई क्षेत्रों के महत्वपूर्ण कार्यों को शामिल करते हैं, जो अनुसंधान के लिए एक ठोस सैद्धांतिक आधार प्रदान करते हैं।