Optimal bounds for numerical approximations of infinite horizon problems based on dynamic programming approach

ملخص تقني: الحدود المثلى للتقريبات العددية لمسائل الأفق اللانهائي القائمة على نهج البرم la الديناميكية

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

تشير الأدبيات الحالية، وبالتحديد بالرجوع إلى أعمال فالكوني وفريتي [12] وآخرين [9, 11]، إلى حد خطأ بمقدار $O(k/h)$ لطرق "سي-لاغرانجيان" شبه المعتمدة المتقطعة كلياً، حيث $h$ هو الخطوة الزمنية و $k$ هو مقياس الشبكة المكانية. يشير هذا الحد إلى أنه كلما تناقصت الخطوة الزمنية $h$، تدهور الخطأ بسبب عامل $1/h$. ومع ذلك، يلاحظ المؤلفون أن التجارب العددية في الأدبيات (على سبيل المثال [4]) لا ترصد سلوك $1/h$ هذا، مما يوحي بأن الحد النظري دون الأمثل. تهدف الورقة إلى حل هذا التناقض من خلال مراجعة تحليل الخطأ لإثبات حد أكثر إحكاماً.

المنهجية
يقترح المؤلفون توصيفاً جديداً للتقريب المتقطع كلياً لاستخلاص حدود خطأ محسنة. تسير المنهجية كما يلي:

1. إعادة الصياغة عبر الدوال التكلفية: بدلاً من الاعتماد فقط على مبدأ البرمجة الديناميكية القياسي، يعرّف المؤلفون دالة تكلفة متقطعة كلياً $\hat{J}_{h,k}$. تستخدم هذه الدالة مؤثر استكمال خطي $I_k$ مطبقاً على ديناميكيات النظام وتكلفة الجريان عند كل خطوة زمنية.
2. توصيف الحل: يثبت النظرية 3 أن دالة القيمة المتقطعة كلياً $v_{h,k}$ (المعرفة بمعادلة البرمجة الديناميكية المتقطعة) هي بالضبط الحد الأدنى لهذه الدالة التكلفية الجديدة المتقطعة كلياً عبر فضاء متتاليات التحكم المقبولة. وهذا يماثل التوصيف المستخدم في الحالة شبه المتقطعة في [5].
3. تفكيك الخطأ: يقارن تحليل الخطأ بين دالة التكلفة المستمرة $\hat{J}$ والدالة المتقطعة كلياً $\hat{J}_{h,k}$. ويتم حصر الفرق من خلال تحليل أخطاء المسار بين حل النظام المستمر $y(t)$ والحل المستكمل (interpolated) المتقطع $\tilde{y}(t)$، بالإضافة إلى أخطاء الاستكمال لدوال التكلفة والديناميكيات.
4. افتراضات الانتظام: يأخذ التحليل في الاعتبار افتراضات انتظام مختلفة لمدخلات التحكم $u(t)$:
   • التحكمات المتصلة بانتظام (مع عدد محدود من قفزات الانقطاع).
   • التحكمات ليبتشيتز المستمرة.
   • التحكمات الثابتة جزئياً (تحت فرضية التقعر على مجموعة الديناميكيات/التكاليف).

المساهمات والنتائج الرئيسية

• الحد المحسن للخطأ: المساهمة الأساسية هي إثبات أن الخطأ بين دالة القيمة المستمرة $v$ والتقريب المتقطع كلياً $v_{h,k}$ هو من الرتبة $O(h + k)$، بدلاً من $O(k/h)$ المذكورة سابقاً. وهذا يثبت التقارب من الدرجة الأولى في كل من الزمان والمكان، وهو ما يتوافق مع استخدام تقريب "أويلر" الزمني من الدرجة الأولى والاستكمال المكاني الخطي الجزئي.
• النظرية 4 (الحد العام): تحت افتراضات ليبتشيتز القياسية على ديناميكيات النظام $f$ والتكلفة $g$، وبفرض أن $\lambda > L$ (حيث $L$ ثابت مشتق من ثوابت ليبتشيتز)، يكون الخطأ محصوراً بـ:
  $$|v(y) - v_{h,k}(y)| \leq C_1(h + k) + C_2 M_u$$
  حيث تقيس $M_u$ تذبذب التحكم ضمن الفترات الزمنية.
• النظرية 6 (التحكمات الليبتشيتزية): إذا كانت عناصر التحكم مستمرة وفق شرط ليبتشيتز، فإن $M_u$ تتناسب مع $h$، مما يؤدي إلى حد الخطأ الأمثل:
  $$|v(y) - v_{h,k}(y)| \leq C(h + k)$$
• النظرية 5 (التحكمات المتصلة بانتظام): إذا كانت عناصر التحكم متصلة بانتظام فحسب (تسمح بعدد محدود من القفزات)، فإن الخطأ يتقارب إلى الصفر مع اقتراب $h, k \to 0$، وإن كان المعدل المحدد يعتمد على معامل الاستمرار.
• التحكمات الثابتة جزئياً (النظرية 7): من خلال تكييف حجج [6] (المخصصة أصلاً لمسائل الأفق المحدود) لتناسب حالة الأفق اللانهائي، يثبت المؤلفون التقارب تحت فرضية التقعر (CA) على المجموعة $\{f(y,u), g(y,u)\}$. ومع ذلك، يؤدي هذا التكييف إلى معدل تقارب زمني أضعف لحالة الأفق اللانهائي:
  $$|v(y) - v_{h,k}(y)| \leq C_1(h + k) + C_2 (\log h)^2 h^{\frac{1}{1+\beta}}$$
  حيث $\beta = \frac{\sqrt{n}L_f}{\lambda}$. ويشير المؤلفون إلى أن هذا يمثل خسارة لا تزيد عن نصف رتبة في التقارب الزمني مقارنة بحالة الأفق المحدود، ويُعزى ذلك إلى الأفق الزمني اللانهائي.

الأهمية والادعاءات
تدعي الورقة أن التوصيف الجديد القائم على حجج الأمثلية يفسر بنجاح السلوك العددي الملحوظ في الممارسة العملية، حيث لا يُلاحظ تدهور $1/h$ الذي تنبأ به الحد القديم. ومن خلال إثبات حد $O(h+k)$، يؤكد المؤلفون أن طريقة "سي-لاغرانجيان" المتقطعة كلياً تحقق تقارباً من الدرجة الأولى في كل من الزمان والمكان تحت افتراضات الانتظام القياسية.

يشير المؤلفون بتواضع إلى أن هذا التوصيف الجديد، المتجذر في حجج الأمثلية، يمتلك إمكانات للتطبيق على أنواع أخرى من معادلات هاملتون-جاكوبي-بيل حيث تكون معدلات التقارب الحالية دون الأمثل. ويذكرون صراحة أنه بينما يستعيدون التقارب من الدرجة الأولى للتحكمات الليبتشيتزية، فإن الانتظام الأضعف (مثل الاستمرار من نوع هولدر بمعامل $\alpha$) يؤدي إلى معدلات تقارب من رتبة $O(h^\alpha)$، كما أن نهج التحكم الثابت جزئياً للأفق اللانهائي يتضمن عقوبة لوغاريتمية مقارنة بنتائج الأفق المحدود. لا تقترح هذه الدراسة تجارب جديدة، بل تقدم تصحيحاً نظرياً لتحليل الخطأ الحالي في هذا المجال.