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

تقليل رتبة النموذج الموضعي للأنظمة الهاميلتونية ذات المنافذ

المعلومات الأساسية

  • معرّف الورقة: 2203.07751
  • العنوان: تقليل رتبة النموذج الموضعي للأنظمة الهاميلتونية ذات المنافذ
  • المؤلفون: Silke Glas, Mir Mamunuzzaman, Hongliang Mu, Hans Zwart
  • التصنيف: math.OC (التحسين والتحكم)، math.CA (التحليل الكلاسيكي)
  • وقت النشر: مارس 2022 (نسخة أولية من arXiv)
  • رابط الورقة: https://arxiv.org/abs/2203.07751

الملخص

تقدم هذه الورقة طريقة جديدة لتقليل رتبة النموذج مع الحفاظ على البنية (MOR) للأنظمة الهاميلتونية ذات المنافذ (pH) الخطية غير المتغيرة مع الزمن. يهدف البحث إلى بناء نظام pH مختزل يحافظ على قابلية التفسير في المجال الفيزيائي للنموذج الكامل الأصلي. بشكل محدد، إذا كان النظام عالي الأبعاد الأولي عبارة عن دائرة كهربائية، فنحن نرغب في أن يظل النموذج المختزل قابلاً للتفسير كدائرة إلكترونية. بالنسبة لأنظمة الكتلة-الزنبرك-التخميد (MSD)، توجد طرق MOR تضمن الحفاظ على هذه البنية المحددة. تُظهر هذه الورقة أن طريقة MOR الجديدة المحافظة على البنية بناءً على طرق MOR الموضعية تسترجع طريقة Arnoldi من الدرجة الثانية المعروفة في حالة أنظمة MSD. ومع ذلك، بالنسبة لنماذج pH الدوائر الكهربائية وغيرها ذات البنية الكتلية المماثلة، تنتج الطريقة تقنيات تقليل نموذج جديدة.

الخلفية البحثية والدافع

وصف المشكلة

الأنظمة الهاميلتونية ذات المنافذ هي تعميم للإطار الهاميلتوني الكلاسيكي وتكون مفيدة للغاية في المشاكل الهندسية، لأنها تجمع بين الطرق الهاميلتونية الكلاسيكية والإطار القوي لنمذجة الأنظمة الفيزيائية المفتوحة. ومع ذلك، في تصميم المتحكمات، عادة ما يكون للمتحكم نفس رتبة النظام الأصلي، مما يدفع الحاجة إلى نماذج مختزلة يمكنها التقاط السلوك الأساسي للنظام الأصلي.

حدود الطرق الموجودة

  1. مشكلة فقدان البنية: تطبيق طرق MOR التقليدية مباشرة يؤدي إلى فقدان بنية pH
  2. غياب القابلية للتفسير الفيزيائي: الطرق الموجودة لا تضمن قابلية التفسير الفيزيائي للنموذج المختزل
  3. تدمير البنية الكتلية: بالنسبة لأنظمة pH ذات البنية الكتلية الخاصة، لا تستطيع الطرق الموجودة الحفاظ على هذه البنية

الدافع البحثي

عند إعادة صياغة النماذج الفيزيائية في شكل pH، تمتلك المصفوفات J و R و H عادة بنية كتلية خاصة. الهدف من هذه الورقة هو الحفاظ على هذه البنية الكتلية أثناء عملية التقليل، مما يضمن أن النموذج المختزل لا يزال قابلاً للتفسير كنظام فيزيائي من نفس نوع النموذج الكامل.

المساهمات الأساسية

  1. إدخال طريقة Symplectic-pH: اقتراح طريقة جديدة لـ MOR لأنظمة pH بناءً على MOR الموضعية، قابلة للتطبيق على البنية الكتلية (2)، بحيث يحافظ النموذج المختزل على نفس شكل نظام pH
  2. الربط النظري: إثبات أنه في الحالة الخاصة لأنظمة MSD، تسترجع الطريقة المقترحة طريقة Arnoldi من الدرجة الثانية (SO-Arnoldi)
  3. التحقق العددي: عرض النتائج العددية على الدوائر الإلكترونية (دائرة LRCR) مع المقارنة بطرق أخرى معروفة من الأدبيات
  4. مزايا الاستقرار: اكتشاف أن طريقة 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 يحقق:

  1. عدم التماثل: Ω(v₁, v₂) = -Ω(v₂, v₁)
  2. عدم التحلل: لجميع v₁ ∈ V، Ω(v₁, v₂) = 0 إذا وفقط إذا كان v₂ = 0

التعريف 2 (الخريطة الموضعية): تُسمى المصفوفة Q بـ (J₂ₙ, J₂ₖ)-موضعية إذا:

Q*J₂ₙQ = J₂ₖ

تقليل رتبة النموذج الموضعي للحالة بدون تبديد

النظرية 4: بالنسبة لنظام pH بدون تبديد (R = 0)، إذا كانت Q ∈ C²ⁿˣ²ᵏ مصفوفة (J₂ₙ⁻¹, J₂ₖ⁻¹)-موضعية، وتحقق:

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

فإن نظام pH المختزل يساوي النظام الأصلي عند التردد المعطى.

تقليل رتبة النموذج الموضعي للحالة مع التبديد

النظرية 5: بالنسبة لنظام pH مع التبديد، نقوم بإنشاء مصفوفة موسعة:

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

حيث Q₂ هي (J₂ₙ, J₂ₖ)-موضعية، و Q₁* هي معكوس يساري لـ Q₂. تحت الشروط المناسبة، يحقق النظام المختزل شروط الاستيفاء.

خوارزمية الحفاظ على البنية الكتلية

الخوارزمية 1 (Symplectic-pH MOR):

  1. بناء V₁ بناءً على بيانات الاستيفاء
  2. حساب V₂ := HV₁
  3. تقسيم V₁ و V₂ حسب البنية الكتلية
  4. بناء Ψ₁, Ψ₂, Φ₁, Φ₂ من خلال الصيغة (27)
  5. بناء Q₁ و Q₂
  6. تطبيق إسقاط Petrov-Galerkin

صيغ البناء:

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

نقاط الابتكار التقني

  1. الحفاظ على البنية الموضعية: الحفاظ على البنية الهندسية الأساسية لنظام هاميلتوني من خلال الخرائط الموضعية
  2. صيانة البنية القطرية الكتلية: بناء خاص لـ Q₁ و Q₂ يضمن الحفاظ على البنية الكتلية ذات المعنى الفيزيائي
  3. شروط الاستيفاء: ضمان أن النموذج المختزل يطابق النموذج الأصلي بدقة عند نقاط تردد محددة مسبقاً
  4. إطار موحد: تضمين طريقة SO-Arnoldi لأنظمة MSD كحالة خاصة

إعداد التجارب

مثال عددي: دائرة LRCR

النظر في دائرة LRCR الموضحة في الشكل 1، تحتوي على 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ω)‖₂

طرق المقارنة

  1. pH-Arnoldi: طريقة الحفاظ على البنية بناءً على النظرية 2
  2. SO-Arnoldi: طريقة Arnoldi من الدرجة الثانية
  3. Symplectic-pH: الطريقة المقترحة في هذه الورقة

اختيار نقاط الاستيفاء

اختيار نقاط الاستيفاء بشكل موحد على مقياس لوغاريتمي في الفترة 10⁻², 10²، على سبيل المثال 6 نقاط استيفاء: 10⁻²i, -10⁻²i, 10⁰i, -10⁰i, 10²i, -10²i

نتائج التجارب

النتائج الرئيسية

من الشكل 3 يمكن ملاحظة:

  1. Symplectic-pH مقابل pH-Arnoldi: توفر Symplectic-pH خطأ H₂ و H∞ أقل في معظم الحالات
  2. Symplectic-pH مقابل SO-Arnoldi:
    • الأداء متشابهة عندما تكون رتبة التقليل <24
    • عند رتبة التقليل =32 تنخفض أداء SO-Arnoldi بشكل حاد
    • عند رتب التقليل 28 و 36 تنتج SO-Arnoldi أنظمة غير مستقرة

تحليل الاستقرار

يوضح الشكل 4 الأقطاب غير المستقرة التي تنتجها طريقة SO-Arnoldi عند r=28 و r=36، بينما تحافظ طريقة Symplectic-pH على الاستقرار في جميع أبعاد الاختبار.

التحقق النظري: تكافؤ أنظمة MSD

النظرية 7 تثبت أنه بالنسبة لأنظمة MSD، طريقة Symplectic-pH مكافئة لطريقة SO-Arnoldi من خلال تحويل وحدوي: يوجد مصفوفة وحدوية Θ بحيث zred(t) = Θẑred(t)

الأعمال ذات الصلة

تطور MOR لأنظمة pH

  1. الأعمال المبكرة: طرق الأنظمة السلبية بناءً على دوال النقل الموجبة الحقيقية
  2. مطابقة اللحظات: أعمال Polyuga و van der Schaft
  3. تقريب بنية Dirac: مسار مختلف للحفاظ على البنية
  4. طرق التحلل الطيفي: طرق بناءً على دالة Popov
  5. طرق التحسين: طرق بناءً على تحسين المعاملات والتحسين الريماني

تطور MOR الموضعي

  1. الأنظمة الهاميلتونية المستقلة: العمل الرائد لـ Peng و Mohseni
  2. الطرق غير الخطية: تقليل الرتبة الموضعي بناءً على المتشعبات
  3. الطرق المدفوعة بالبيانات: استدلال المشغلين وتحديد النظام

الخلاصة والنقاش

الاستنتاجات الرئيسية

  1. تطوير ناجح لطريقة MOR تحافظ على بنية pH والبنية الكتلية
  2. توحيد نظري لطرق MOR المختلفة لأنظمة MSD
  3. التحقق العددي من فعالية الطريقة ومزايا الاستقرار

القيود

  1. طريقة الاستيفاء: حالياً تطابق فقط قيم دالة النقل، لم تأخذ في الاعتبار مطابقة المشتقات
  2. قيود البنية الكتلية: الطريقة مصممة خصيصاً لبنية كتلية محددة
  3. التعقيد الحسابي: تحتاج تكلفة حسابية لبناء المصفوفات الموضعية إلى تحليل إضافي

الاتجاهات المستقبلية

  1. مطابقة المشتقات: التوسع لمطابقة مشتقات دالة النقل لتحسين جودة الملاءمة
  2. الاختيار التكيفي: تطوير استراتيجيات اختيار تكيفية لنقاط الاستيفاء
  3. التوسع غير الخطي: التوسع إلى أنظمة pH غير الخطية
  4. التطبيقات واسعة النطاق: تحسين الخوارزمية للتعامل مع أنظمة أكبر حجماً

التقييم المتعمق

المزايا

  1. الصرامة النظرية: أساس رياضي متين بناءً على الهندسة الموضعية
  2. الحفاظ على البنية: الحفاظ المتزامن على بنية pH والبنية الكتلية ذات المعنى الفيزيائي
  3. التوحيد: توحيد الطرق الموجودة كحالات خاصة في إطار واحد
  4. الجدوى العملية: الطريقة لها تطبيق خوارزمي واضح والتحقق العددي

أوجه القصور

  1. نطاق التطبيق: الطريقة مخصصة لبنية كتلية محددة، الشمولية محدودة
  2. الكفاءة الحسابية: تحليل التعقيد الحسابي لبناء المصفوفات الموضعية غير كافٍ
  3. التحقق واسع النطاق: نقص التحقق العددي للأنظمة الأكبر حجماً

التأثير

  1. المساهمة النظرية: توفير منظور نظري جديد لـ MOR لأنظمة pH
  2. القيمة العملية: لها قيمة مباشرة في التطبيقات الهندسية مثل تحليل الدوائر
  3. الأهمية المنهجية: توضيح الدور المهم للطرق الهندسية في MOR

السيناريوهات المناسبة

  1. تحليل الدوائر: النمذجة المختزلة للدوائر الإلكترونية واسعة النطاق
  2. الأنظمة الميكانيكية: تصميم التحكم لأنظمة الكتلة-الزنبرك-التخميد
  3. الأنظمة متعددة الفيزياء: الأنظمة الفيزيائية المعقدة ذات خصائص الحفاظ على الطاقة

المراجع

تتضمن الورقة 31 مرجعاً، تغطي أعمالاً مهمة في عدة مجالات بما فيها نظرية MOR، نمذجة أنظمة pH، الهندسة الموضعية ونظرية التحكم، مما يوفر أساساً نظرياً متيناً للبحث.