A Semi-smooth Newton Method for the Constrained Optimal Control of Continuous-Time Linear Systems

ملخص تقني: طريقة نيوتن شبه سلسة للتحكم الأمثل المقيد للأنظمة الخطية في الزمن المستمر

بيان المشكلة
يتناول البحث تحدي حل مسائل التحكم الأمثل (OCPs) المقيدة للأنظمة الخطية في الزمن المستمر مباشرة ضمن فضاء الدوال ذي البعد اللانهائي. يسعى صياغة المشكلة المحددة إلى تقليل دالة التكلفة التي تتكون من تكلفة نهائية $J(x(T))$ وتكامل لتكلفة متزايدة $\ell(x(t), u(t))$، مع الخضوع لديناميكيات خطية $\dot{x} = Ax + Bu$، وشروط أولية، وقيود نقطية غير متساوية $c(x(t), u(t)) \leq 0$. يركز المؤلفون على المسائل التربيعية الخطية ذات الأفق المحدود، حيث تنتمي مسارات الحالة والمدخلات إلى فضاءات سوبوليف وليبيج ($W^{1,\infty}$ و $L^\infty$)، مما يضمن سلامة المسألة تحت افتراضات التحدب القياسية وتأهيل القيود.

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

1. إعادة صياغة KKT عبر دالات NCP: يتم دمج شروط كاروش-كون-تاكر (KKT) الضرورية للتحسين داخل دالة تكامل غير سلسة. وتحديداً، يتم استبدال شروط التكامل التي تتضمن المتغيرات المزدوجة $\mu$ والقيود $c$ بدالة مسألة تكامل غير خطية (NCP) $\phi(a, b) = 0$، والتي تكافئ $a, b \geq 0, ab=0$. يستخدم المؤلفون إما دالة "الحد الأدنى" (min) أو دالة فيشر-بورميستر. يؤدي هذا إلى تحويل شروط التحسين إلى معادلة غير سلسة $F(z) = 0$، حيث تمثل $z = (x, u, \lambda, \mu)$ الحالة، والمدخل، والمتغير المصاحب، والمتغيرات المزدوجة.

2. تكرار نيوتن شبه السلس (Semi-Smooth Newton Iteration): يتم حل مسألة إيجاد الجذر $F(z)=0$ باستخدام طريقة نيوتن شبه السلسة. تولد هذه العملية التكرارية تحديثاً $\tilde{z}$ عن طريق حل معادلة خطية تتضمن جاكوبي كلايك (Clarke generalized Jacobian) للدالة $F$. وقد ثبت أن الطريقة تتقارب بشكل فائق (superlinearly) نحو الحل تحت شروط انتظام طفيفة.

3. حل خطوات نيوتن القائم على ريكاتي: المساهمة التقنية المركزية هي اشتقاق طريقة فعالة حاسوبياً لحل معادلة خطوة نيوتن. يوضح المؤلفون أن مسألة القيمة الحدية ذات النقطتين والخطية الناشئة عن التخطي يمكن إعادة صياغتها كمسألة منظم لترتيب خطي-تربيعي (LQR) متغير زمنياً وغير مقيد.

   • يتم دمج قيود عدم التساوي بفعالية في هذه المسألة الفرعية كحدود تكلفة "حاجز" (barrier)، حيث يتم إعادة وزن مصفوفات التكلفة ($Q, R, S$) والحدود الخطية ($q, r$) في كل تكرار بناءً على مضاعفات القيود الحالية والمشتقات المعممة لدالة NCP.
   • تسمح هذه إعادة الصياغة بحساب تحديث نيوتن عبر حل معادلة ريكاتي التفاضلية إلى الخلف في الزمن، يتبعها تكامل أمامي لمسار الحالة. وهذا يستفيد من حلول المعادلات التفاضلية العادية (ODE) القياسية وعالية الكفاءة بدلاً من اشتراط مخططات تقطيع معقدة.

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

• طريقة غير مباشرة في فضاء الدوال: يفصل البحث خوارزمية غير مباشرة مبتكرة تعمل مباشرة في فضاء الدوال في الزمن المستمر، متجنبة نهج "التقطيع ثم التفاضل" الشائع في الطرق المباشرة.
• تكامل ريكاتي للمسائل المقيدة: المساهمة الأساسية هي إثبات أن مسألة القيمة الحدية المعقدة ذات النقطتين المتأصلة في طرق نيوتن شبه السلسة لمسائل التحكم الأمثل المقيدة يمكن حلها تحليلياً عبر معادلات ريكاتي التفاضلية. هذا يربط بين خوارزميات نيوتن شبه السلسة النظرية (مثل Gerdts, 2008) واستراتيجيات التنفيذ عالية الكفاءة المستخدمة في الطرق غير المقيدة مثل PRONTO.
• إعادة وزن التكلفة التكيفية: تفسر الطريقة خطوة نيوتن كبحث تكيفي حيث يتم تحديث الأوزان الخطية والتربيعية لدالة التكلفة تكرارياً بناءً على مضاعفات القيود، مما يعالج القيود بفعالية دون الحاجة إلى تحديد المجموعة النشطة (active-set) أو حواجز النقطة الداخلية بالمعنى التقليدي.

النتائج
يتحقق المؤلفون من صحة الطريقة من خلال محاكاة عددية لمناورة تغيير المسار في مسار الطوارئ لنموذج مركبة خطي (يمثل سيارة BMW 740i موديل 2017). يخضع النظام لقيود على زوايا الانزلاق الأمامية/الخلفية وزاوية التوجيه.

• التقارب: تظهر الخوارزمية تقارباً فائقاً. حيث تنخفض قيمة معيار $L^2$ لمتجه البواقي بأكثر من أربعة مراتب عشرية خلال أول 10 تكرارات.
• الاستقرار: تتقارب الطريقة حتى تصل إلى حد التسامح الخاص بحلال المعادلات التفاضلية العادية (دقة الآلة)، مما يؤكد أن معيار التوقف محدد جيداً.
• المقارنة: أظهرت المحاكاة باستخدام كل من دالتي "الحد الأدنى" ودالة فيشر-بورميستر سلوكيات تقارب متشابهة، وإن كانت دالة فيشر-بورميستر قد استُخدمت في رسومات المسار الرئيسية.

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