On maximal positive invariant set computation for rank-deficient linear systems

إليك ملخص تقني مفصل للورقة البحثية بعنوان "حساب مجموعة الثبات الموجب القصوى للأنظمة الخطية ذات الرتبة الناقصة" (On maximal positive invariant set computation for rank-deficient linear systems).

1. بيان المشكلة

تتناول الورقة البحثية التحدي الحسابي المتمثل في تحديد المجموعة الموجبة الثابتة القصوى (MPI) للأنظمة الخطية ذات الزمن المنفصل حيث تكون مصفوفة الديناميكيات في الحلقة المغلقة ناقصة الرتبة (rank-deficient).

• السياق: في التحكم التنبئي بالنموذج (MPC)، تعد مجموعة الـ MPI ضرورية كقيد نهائي لضمان الجدوى التكرارية والاستقرار. ويتم حسابها عادةً عبر إجراء الوصول الخلفي (backward reachability) الذي يتضمن التقاطع التكراري لمجموعات السوابق: $\Omega_{k+1} = A_\circ^{-1}\Omega_k \cap \mathcal{X}$.
• المشكلة: غالبًا ما تؤدي استراتيجيات تخليق التحكم الأمثل (مثل حل معادلة ريكاتي الجبرية المنفصلة - DARE) أو استراتيجيات التحكم "الضربي الميت" (dead-beat control) إلى مصفوفات حلقة مغلقة ($A_\circ = A + BK$) ذات قيم ذاتية عند الصفر.
• النتيجة: عندما تكون $A_\circ$ منفردة (ناقصة الرتبة)، تصبح علاقة التكرار القياسية غير مستقرة عددياً أو غير محددة هندسياً. فالديناميكيات تعمل فعلياً على إسقاط فضاء الحالة على فضاءات فرعية ذات أبعاد أقل خلال الخطوات المتوسطة، مما يؤدي إلى فشل الخوارزميات القياسية (التي تستخدم تمثيلات متعددة السطوح أو الزونوتوب) أو جعلها مكلفة حسابياً بسبب انفجار القيود والتكرارات المطلوبة للتقارب.

2. المنهجية

يقترح المؤلفون خوارزمية قوية تستغل تحليل شور الحقيقي (Real Schur Decomposition) لمصفوفة الحلقة المغلقة للتعامل مع حالات التفرد بشكل صريح. الفكرة الجوهرية هي فصل ديناميكيات النظام إلى فضاء فرعي مستقر وغير منفرد، وفضاء فرعي "نيلي" (nilpotent - ذو قيم ذاتية صفرية).

أ. الأساس الرياضي

1. تحليل شور: يتم تفكيك مصفوفة الحلقة المغلقة إلى $A_\circ = USU^\top$، حيث $U$ مصفوفة متعامدة و $S$ مصفوفة مثلثية علوية.
2. التقسيم الكتلي: يتم تقسيم المصفوفة $S$ بحيث تُعزل القيم الذاتية الصفرية (والكتل جوردان المرتبطة بها) في الكتلة السفلية اليمنى ($S_{22}$)، بينما توجد القيم الذاتية غير الصفرية في الكتلة العلوية اليسرى ($S_{11}$):
   $$S = \begin{bmatrix} S_{11} & S_{12} \\ 0 & S_{22} \end{bmatrix}$$
   هنا، $S_{22}$ هي مصفوفة نيلية من الدرجة $p+1$ (أي أن $S_{22}^{p+1} = 0$).
3. تحويل الإحداثيات: يتم تحويل متجه الحالة إلى $y = U^\top x$، مما يفصل الديناميكيات إلى مكونين ($y_1$ و $y_2$). المكون $y_2$ يتلاشى بعد $p$ من الخطوات بسبب الطبيعة النيلية.

ب. الخوارزمية المقترحة

تحسب الخوارزمية مجموعة الـ MPI على مرحلتين وتتقاطع النتائج بينهما:

1. المرحلة 1 (قيود المرحلة الانتقالية): حساب تقاطع القيود لأول $p$ من التكرارات ($k < p$). خلال هذه المرحلة، لا تزال الديناميكيات النيلية نشطة، ويجب فرض قيود الحالة الكاملة.
2. المرحلة 2 (الديناميكيات المختزلة): لـ $k \geq p$، تتقلص ديناميكيات النظام فعلياً إلى $z_{k+1} = S_{11}z_k$ في فضاء فرعي أقل أبعاداً ($\mathbb{R}^{d_1}$).
   • يحسب المؤلفون مجموعة الـ MPI ($\Phi_z$) لهذا النظام المختزل كامل الرتبة.
   • ثم يتم "رفع" هذه المجموعة المختزلة إلى الفضاء الأصلي باستخدام مصفوفة مساعدة $T$ (مشتقة من $S_{11}, S_{12}, S_{22}$) لفرض القيود على الإحداثيات الأصلية.
3. التقاطع النهائي: مجموعة الـ MPI النهائية هي تقاطع قيود المرحلة الانتقالية (المرحلة 1) والمجموعة المختزلة المرفوعة (المرحلة 2).

ج. التمثيلات

تمت صياغة الطريقة لتمثيلين للمجموعات:

• المجموعات متعددة السطوح (Polyhedral Sets): باستخدام متباينات نصف الفضاء ($Fx \leq \theta$).
• الزونوتوبات المقيدة (Constrained Zonotopes - CZ): باستخدام مصفوفات المولدات وقيود المساواة. يستفيد المؤلفون من خاصية أن الـ CZ مغلقة تحت عملية التقاطع والتحويل الخطي، مما يسمح بإجراء التكرار بدقة دون تقريب.

3. المساهمات الرئيسية

• التعامل مع الديناميكيات المنفردة: تقدم الورقة أول إطار نظري وخوارزمية لحساب مجموعات الـ MPI للأنظمة ذات القيم الذاتية الصفرية في الحلقة المغلقة، وهي حالة تفشل فيها طرق الوصول الخلفي القياسية.
• تقليل الأبعاد: من خلال استخدام شكل شور، يتم تنفيذ الجزء الأكثر استهلاكاً للوقت في التكرار (التكرار ذو الأفق اللانهائي) في أبعاد مختزلة ($d_1 < n$)، مما يقلل التعقيد بشكل كبير.
• دعم التمثيل المزدوج: الخوارزمية قابلة للتطبيق على كل من التمثيلات متعددة السطوح والزونوتوبات المقيدة، مما يوفر مرونة لمختلف تطبيقات التحكم.
• المتانة: تتجنب الطريقة عدم الاستقرار العددي المرتبط بعكس المصفوفات المنفردة من خلال العمل مباشرة مع البنية النيلية.

4. النتائج والتقييم

تحقق المؤلفون من صحة نهجهم باستخدام مثالين:

1. مثال عددي ($\mathbb{R}^6$):

   • نظام ذو مصفوفة حلقة مغلقة تمتلك القيم الذاتية $\{0, 0, 0, 0.2, 0.5, 0.7\}$.
   • نجحت الخوارزمية في حساب مجموعة الـ MPI في 0.4 ثانية، حيث تعاملت بشكل صحيح مع نقص الرتبة ($d_2=3$) عن طريق تقليص المشكلة إلى فضاء فرعي ثلاثي الأبعاد.

2. تجربة النابض المزدوج (Coupled Spring Experiment - CSE):

   • نظام قابل للتوسع مع عدد متغير من الكتل ($\ell$)، حيث يزدير أبعاد الحالة حتى $\mathbb{R}^{20}$.
   • المقارنة: تمت مقارنة طريقة "نقص الرتبة" المقترحة مع خوارزمية قياسية (التي تفشل أو تصبح غير فعالة في الحالات المنفردة).
   • مكاسب الأداء:
     • التكرارات: تطلبت الطريقة القياسية زيادة هائلة في التكرارات مع نمو الأبعاد، بينما أظهرت الطريقة المقترحة نمواً معتدلاً.
     • وقت الحساب: بالنسبة لـ $\ell \geq 8$، أصبح الوقت المطلوب للطريقة القياسية غير مجدٍ حسابياً، بينما ظلت الطريقة المقترحة فعالة.
     • التعقيد: كان عدد المتباينات المطلوبة لوصف المجموعة متعددة السطوح الناتجة أقل بكثير للطريقة المقترحة (على سبيل المثال، 704 مقابل 1496 متباينة لـ $\ell=10$).

5. الأهمية

يحل هذا العمل عقبة حرجة في تصميم التحكم التنبئي بالنموذج (MPC) للأنظمة التي تستخدم قوانين تحكم مثلى تسبب بطبيعتها نقصاً في الرتبة (مثل التحكم القائم على معادلة ريكاتي أو التحكم الضربي الميت).

• تمكين الـ MPC القوي: يسمح للمهندسين باستخدام متحكمات مثالية نظرياً تضع الأقطاب عند الصفر دون التضحية بالقدرة على حساب المجموعات الثابتة النهائية اللازمة لضمان الاستقرار.
• الكفاءة الحسابية: من خلال تقليل أبعاد مشكلة الوصول، تجعل هذه الطريقة حساب المجموعات الثابتة أمراً ممكناً للأنظمة ذات الأبعاد العالية حيث قد تفشل الطرق القياسية.
• الصرامة النظرية: تجسر هذه الورقة الفجوة بين الهندسة الجبرية (أشكال شور/جوردان) والتحسين المحدب (Set Invariance)، مما يوفر حلاً سليماً رياضياً للأنظمة المنفردة.

في الختام، تقدم هذه الورقة تقدماً ضرورياً لنظرية التحكم الحديثة، مما يضمن إمكانية حساب خاصية "القصوى" للمجموعات الثابتة بشكل موثوق حتى عندما تكون ديناميكيات النظام الأساسية منفردة.