Nonlinear MPC for Feedback-Interconnected Systems: a Suboptimal and Reduced-Order Model Approach

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

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

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

يركز المؤلفون على فئة محددة من الأنظمة المترابطة بالتغذية الراجعة في الزمن المنفصل والتي تتميز بوجود زمنين مختلفين:

1. ديناميكيات الهدف ($x_t$): الحالة الأساسية المراد التحكم فيها (الديناميكيات البطيئة).
2. الديناميكيات الإضافية ($\xi_t$): الديناميكيات الأسرع، والتي قد تمثل تأثيرات غير مبرمجة، أو أنظمة فرعية عالية التردد، أو معوضات منخفضة المستوى.

التحدي الجوهري هو تصميم بنية تحكم تحقق ما يلي:

• استخدام نموذج منخفض الرتبة (بتجاهل الديناميكيات السريعة $\xi_t$) لتبسيط عملية الأمثلة.
• العمل دون المثالية (عبر إجراء عدد محدود فقط من تكرارات الأمثلة بدلاً من الحل للوصول إلى المثالية التامة) لتقليل الحمل الحسابي.
• ضمان استقرار النظام الكامل ذي الحلقة المغلقة رغم اختزال النموذج وإجراءات التحكم دون المثالية.

2. المنهجية

أ. نمذجة النظام

يتم نمذجة النظام كترابط تغذية راجعة في الزمن المنفصل:
$$\begin{aligned} x_{t+1} &= f(x_t, \xi_t, u_t, \delta) \\ \xi_{t+1} &= g(\xi_t, x_t, u_t, \delta) \end{aligned}$$
حيث $\delta > 0$ هو معامل يمثل زمن أخذ العينات. يستفيد المؤلفون من فصل المقاييس الزمنية، حيث يعمل $\delta$ كمعامل صغير يعمل على تسريع الديناميكيات السريعة بالنسبة للديناميكيات البطيئة.

• الافتراض 1: الديناميكيات السريعة هي مستمرة ليبشيتز (Lipschitz continuous) بالنسبة للحالات البطيئة والمدخلات، ومدرجة بمقياس $\delta$.
• الافتراض 2: يمتلك النظام الفرعي السريع سطح اتزان مستقر أسياً عالمياً $\xi_{eq}(x, u)$، بشكل موحد بالنسبة لـ $(x, u)$.

ب. استراتيجية التحكم

تجمع بنية التحكم المقترحة بين استراتيجيتين للكفاءة:

1. نموذج منخفض الرتبة: يستخدم مُحسن الـ MPC نموذجاً مبسطاً $f_R$ يفترض أن الحالة السريعة $\xi$ موجودة فورياً عند اتزانها $\xi_{eq}(x, u)$.
   $$f_R(x, u, \delta) := f(x, \xi_{eq}(x, u), u, \delta)$$
2. الأمثلة دون المثالية: بدلاً من حل مسألة التحكم الأمثل ذات الأفق المحدود للوصول إلى المثالية العالمية، يقوم الحل (solver) بإجراء عدد ثابت ومحدود من التكرارات (على سبيل المثال، خطوة واحدة من طريقة التدرج المسقط). تتطور حالة المُحسن $z_t$ وفقاً لديناميكيات $A(z_t, x_t)$، وصولاً إلى الحل الأمثل $z^*(x_t)$.

ج. إطار تحليل الاستقرار

يستخدم المؤلفون نظرية الاضطراب المنفرد (فصل المقاييس الزمنية) المعدلة للأنظمة في الزمن المنفصل لإثبات الاستقرار. يقوم التحليل بتفكيك نظام الحلقة المغلقة إلى:

1. نظام الطبقة الحدية (النظام الفرعي السريع): يحلل تقارب الحالة السريعة $\xi_t$ وحالة المُحسن $z_t$ إلى أسطح الاتزان الخاصة بهما، بافتراض ثبات الحالة البطيئة $x_t$.
2. النظام المختزل (النظام الفرعي البطيء): يحلل استقرار الحالة البطيئة $x_t$ بافتراض أن الحالات السريعة في حالة اتزان تام وأن المُحسن مثالي.
3. النظام المترابط: يجمع بين النتائج أعلاه لإظهار أنه إذا كان زمن أخذ العينات $\delta$ صغيراً بما يكفي، فإن "الانحراف" الناتج عن عدم التطابق بين الديناميكيات السريعة الفعلية والاتزان المثالي سيكون صغيراً بما يكفي للحفاظ على الاستقرار الأسي العالمي.

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

• بنية موحدة: تقترح الورقة أول إطار عمل يدمج في آن واحد بين التحكم دون المثالي واختزال رتبة النموذج للأنظمة المترابطة بالتغذية الراجعة غير الخطية.
• إثبات صارم للاستقرار: تثبت الورقة الاستقرار الأسي العالمي لنظام الحلقة المغلقة الناتج عن الترابط بين نبات (plant) كامل الرتبة ومُحسن الـ MPC منخفض الرتبة ودون المثالي.
• زمن أخذ العينات كمعامل استقرار: يحدد المؤلفون رسمياً زمن أخذ العينات $\delta$ ليس فقط كمعامل تقسيم (discretization parameter)، بل كـ معامل حرج لفصل المقاييس الزمنية. ويثبتون أنه يوجد حد $\bar{\delta}$ بحيث لأي $\delta < \bar{\delta}$، يظل النظام مستقراً.
• التعامل مع عدم تطابق النموذج: يأخذ النظرية في الاعتبار صراحةً عدم التطابق الهيكلي بين النموذج منخفض الرتبة المستخدم في المُحسن والواقع كامل الرتبة للنبات.

4. النتائج

تم التحقق من النتائج النظرية من خلال عمليات محاكاة عددية على نظام ميكاتروني: بندول يتم تشغيله بواسطة محرك تيار مستمر (DC motor).

• إعداد النظام: زاوية البندول والسرعة ($x$) هما الحالات البطيئة، بينما يمثل تيار العضو الدوار ($\xi$) الديناميكيات السريعة.
• سيناريوهات المحاكاة: قارن المؤلفون بين نهجهم "دون المثالي ومنخفض الرتبة" وبين:
  1. MPC الأمثل مع نبات مثالي (بدون عدم تطابق).
  2. MPC دون المثالي مع نبات مثالي.
  3. النهج المقترح (دون المثالي + منخفض الرتبة) بالتفاعل مع النبات كامل الرتبة.
• النتائج الرئيسية:
  • $\delta$ صغير (0.01 ثانية): حقق النهج المقترح تقارباً أسياً مشابهاً للمعايير المثالية، مما يؤكد صحة نظرية فصل المقاييس الزمنية.
  • $\delta$ متوسط (0.12 ثانية): تدهورت أداء النهج المقترح عند تفاعله مع النبات كامل الرتبة، بينما ظلت المعايير المثالية مستقرة. سلط هذا على ضرورة شرط فصل المقاييس الزمنية.
  • $\delta$ كبير (0.2 ثانية): فشل النهج المقترح في تحقيق استقرار النظام، في حين ظل الـ MPC الأمثل على النبات المثالي يعمل. أظهر هذا أنه بدون فصل كافٍ للمقاييس الزمنية، تؤدي الأخطاء المجتمعة من اختزال النموذج والأمثلة دون المثالية إلى عدم الاستقرار.

5. الأهمية

هذا العمل مهم للتطبيق العملي للتحكم التنبئي بالنماذج (MPC) في التطبيقات ذات الموارد المحدودة أو السرعات العالية (مثل الروبوتات، الفضاء، إلكترونيات القدرة) حيث:

1. القيود الحسابية تمنع حل مسائل الأمثلة كاملة الرتبة والدقيقة في الوقت الفعلي.
2. تعقيد النظام يتضمن ديناميكيات سريعة غير مبرمجة أو متحكمات منخفضة المستوى يصعب نمذجتها بدقة.

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