Synthesis of Limit Cycles and Reference Tracking via Switching Affine Systems

ملخص تقني: تخليق الدورات الحدية وتتبع المرجع عبر الأنظمة التآلفية التبديلية

بيان المشكلة
تعد السلوكيات الدورية ظاهرة منتشرة في مجالات تتراوح من البيولوجيا إلى الهندسة، ومع ذلك يظل التحليل الرياضي للمذبذبات غير الخطية (مثل كوراماتو أو فان دير بول) مقيداً. إن توصيف الدورات الحدية، لا سيال فيما يتعلق بالوحدانية والاستقرار، غالباً ما يكون ممكناً فقط في حالات خاصة. كما أن التقنيات الحالية القائمة على البيانات، مثل تعلم الآلة، يمكنها تقريب الديناميكيات الدورية ولكنها تفتقر عادةً إلى الخصائص الهيكلية المطلوبة للتحقق الرسمي. علاوة على ذلك، فإن النهج السابقة لتقريب الدورات الحدية باستخدام الأنظمة التبديلية أو الأنظمة التآلفية المجزأة (PAS) كانت مقيدة إلى حد كبير في الإعدادات المستوية ($\mathbb{R}^2$) مع تقسيمات بسيطة مكونة من منطقتين، أو أنها تفشل في تقديم ضمانات للاستقرار العالمي في الأبعاد الأعلى.

الهدف الرئيسي لهذا العمل مزدوج:

1. تطوير طريقة تخليق تقرب الدورات الحدية للمعادلات التفاضلية العادية (ODEs) غير الخطية باستخدام الأنظمة التآلفية التبديلية في أبعاد تعسفية ($\mathbb{R}^{n_x}$)، مع ضمان أن تكون الدورة الحدية الناتجة مستقرة عالمياً.
2. تصميم متحكم لتتبع المرجع لهذه الأنظمة يحقق التقارب التقاربي نحو مسار دوري، حتى عندما يفتقر النظام المرجعي نفسه إلى ضمانات الاستقرار العالمي.

المنهجية

تقسيم فضاء الحالة وتخليق الديناميكيات
تبدأ الطريقة المقترحة بمجموعة من نقاط الحالة المأخوذة $\mathcal{F}$ التي تم جمعها على طول دورة حدية مستقرة لنظام غير خطي. يتم تقسيم فضاء الحالة إلى عدد محدود من المتعددات السطحية (polytopes) $P_i$. وخلافاً للطرق المستوية السابقة، يتعامل هذا النهج مع أبعاد تعسفية ($n_x = 2, 3, >3$) عن طريق تحديد نقطة مركزية مشتركة $x_s$ ومجموعة من المستويات الحدودية (أو الأشعة في 2D) التي تتقاطع عند $x_s$.

• بالنسبة لـ $n_x=2$، يتم رسم خطوط من $x_s$ عبر نقاط عينات مختارة لتحديد المناطق.
• بالنسبة لـ $n_x \geq 3$، يتم أولاً تحديد مستوٍ فائق (hyperplane) لتقليل المسافة إلى نقاط العينات، مما يحدد خط تقاطع $\Gamma$ يمر عبر $x_s$. ثم يتم تشكيل المناطق بواسطة مستويات تحتوي على $\Gamma$ ونقاط عينات محددة.
  تتم نمذجة الديناميكيات داخل كل منطقة $P_i$ كـ $\dot{x} = A_i x + b_i + Bu(t)$. يتم صياغة التخليق كمسألة تحسين عددي مقيدة. الهدف هو تقليل الخطأ بين مسار النموذج ونقاط البيانات المأخوذة. ومن الأهمية بمكان أن عملية التحسين تفرض الانكماش (contractivity) لضمان الاستقرار العالمي:

1. الاستمرارية: يجب أن تكون الحقول المتجهة مستمرة عبر حدود التبديل.
2. دالة ليابونوف مشتركة: يجب أن توجد مصفوفة موجبة محددة $Q$ بحيث يكون $A_i^T Q + Q A_i \prec 0$ لجميع المناطق. يضمن هذا أن النظام انكماشي، مما يضمن أن أي مسار سيتقارب نحو دورة حدية فريدة لها نفس دورة الإشارة المدخلة $u(t)$.

للتعامل مع عدم خطية الأس لـ (matrix exponential) في دالة التكلفة، تستخدم الطريقة تقريبات متسلسلة تايلور وإجراءات مضادة (مثل ضمان أوقات انتقال صغيرة بين المناطق وتحديد القيم المفردة لـ $A_i$) للحفاظ على الدقة دون تعقيد حسابي مفرط.

تتبع المرجع
بالنسبة لمشكلة التتبع، يعالج البحث السيناريوهات التي قد لا يكون فيها المسار المرجعي الدوري المحدد مستقراً عالمياً (على سبيل المثال، إذا تم إغفال قيود الاستقرار أثناء عملية التحديد لتحسين دقة التقريب). يتم تصميم قانون التحكم كـ $u_c(t) = u_r(t) + K_i x_c(t) + w_i - (K_j x_r(t) + w_j)$.
لتحقيق التتبع التقاربي ($\lim_{t\to\infty} ||x_c(t) - x_r(t)|| = 0$)، يقترح المؤلفون مخطط تخليق باستخدام دوال ليابونوف متعددة بدلاً من دالة واحدة مشتركة.

• يتم تعريف دالة ليابونوف تربيعية مجزأة $V_i$ لكل منطقة تحتوي على المسار المرجعي.
• يتم اشتقاق شروط الاستقرار لضمان تناقص دالة ليابونوف بشكل صارم داخل المناطق وبقائها محدودة عبر حدود التبديل، بشرما أن وقت المكوث في كل منطقة يتجاوز حداً أدنى $T_{min}$.
• يقلل هذا النهج من متطلبات وجود مصفوفة $Q$ مشتركة، مما يسمح بنتائج تقارب أقل تحفظاً مقارنة بالطرق التي تعتمد على دالة ليابونوف واحدة.

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

1. التخليق عالي الأبعاد: يوسع البحث عملية تقريب الدورة الحدية عبر الأنظمة التآلفية التبديلية من التقسيمات المستوية ذات المنطقتين إلى فضاءات حالات عالية الأبعاد ($\mathbb{R}^{n_x}$) مع تقسيمات عامة.
2. ضمان الاستقرار العالمي: من خلال الاستفادة من مبدأ الانكماش وفرض دالة ليابونوف مشتركة أثناء تحسين التخليق، تضمن الطريقة أن الدورة الحدية المولدة مستقرة عالمياً وفريدة.
3. تتبع المرجع باستخدام دوال ليابونوف متعددة: يقدم العمل متحكماً لتتبع المرجع للأنظمة التآلفية التبديلية الدورية يستخدم دوال ليابونوف متعددة. ويُعد هذا أول تطبيق لهذا التقنية لتتبع المرجع الدوري في هذا السياق، مما يوفر بديلاً أقل تحفظاً لنهج الدوال المشتركة.
4. قواعد التقسيم البنائية: يوفر البحث قواعد هندسية محددة لتقسيم فضاءات الحالة في 2D و3D والأبعاد الأعلى بناءً على البيانات المأخوذة، مما يضمن توافق المتعددات السطحية الناتجة مع قيود التخليق.

النتائج
تم إثبات فعالية الطرق المقترحة من خلال أمثلة عددية:

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

الأهمية والادعاءات
يزعم المؤلفون أن هذا العمل يعالج فجوة حرجة في الاستقصاء المنهجي للظواهر التذبذبية، لا سيما في تطبيقات مثل الإيقاعات البيولوجية (circadian rhythms) حيث يعد فهم الاستقرار والتزامن أمراً لا غنى عنه. ومن خلال توفير نموذج قابل للتطبيق (الأنظمة التآلفية التبديلية) يحافظ على الخصائص الديناميكية الأساسية (الاستقرار، الوحدانية) مع السماح بالتحقق الصارم، تسهل هذه الطريقة تحليل المذبذبات غير الخطية المعقدة.

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