Hill-Type Stability Analysis of Periodic Solutions of Fractional-Order Differential Equations

ملخص تقني: تحليل استقرار نوع هيل للحلول الدورية للمعادلات التفاضلية ذات الرتب الكسرية

بيان المشكلة
تتناول الورقة البحثية تحليل استقرار الحلول الدورية في المعادلات التفاضلية ذات الرتب الكسرية (FODEs) غير الخطية. ثمة قيد جوهري في معادلات Caputo-type التقليدية، وهو أنها لا تسمح بوجود حلول دورية تماماً؛ إذ تكون الدوال الدورية تحت مشتقات Caputo دورية تقاربية فقط. وبناءً على ذلك، لا يمكن تطبيق أدوات تحليل الاستقرار القياسية، مثل نظرية فلوكيت (Floquet theory) وطريقة مصفوفة هيل (Hill matrix method)، بشكل مباشر على ديناميكيات الاضطراب الخطي للحلول الدورية في أنظمة Caputo. وبينما تسمح معادلات Liouville-Weyl-type (المعرفة بحد سفلي لانهائي، $-\infty$) بوجود حلول دورية تماماً، إلا أنه لم يتم وضع نظرية فلوكيت صارمة لهذه الأنظمة بعد. وقد اعتمدت المحاولات السابقة لتوسيع نظرية فلوكيت لتشمل الأنظمة ذات الرتب الكسرية على تقريبات (مثل افتراض الهويات الضربية لدوال ميتغ-ليفلر) أو فشلت في توصيف حدود الاستقرار وأنواع التفرع (bifurcation) بشكل شامل.

المنهجية
يقترح المؤلفون إطار عمل يعتمد على معادلات Liouville-Weyl ذات الرتب الكسرية لدراسة الحلول الدورية. وتتضمن المنهجية الجوهرية ما يلي:

1. صياغة النظام: تعريف فئة من معادلات Liouville-Weyl التفاضلية غير الخطية ذات الحد السفلي اللانهائي ($-\infty$). يتم صياغة مسألة القيمة الأولية باستخدام دالة تهيئة $x_0(t)$ معرفة على الفترة $(-\infty, t_0]$، تنتمي إلى فضاء دالي محدد $X(-\infty, t_0]$ يضمن وجود حد الإثارة المستمد من تاريخ النظام.
2. الخطيّة (Linearization): يتم تحليل استقرار الحل الدوري $x_p(t)$ عن طريق خطيّة النظام حوله، مما ينتج عنه معادلة تفاضلية كسرية ذات زمن دوري خطي (LTP-FODE) مع مصفوفة جاكوبي (Jacobian matrix) دورية $J(t)$.
3. خصائص مشتقة Liouville-Weyl: باستخدام خاصية أن مشتقة Liouville-Weyl للدالة الأسية $e^{st}$ هي $s^\alpha e^{st}$ (لـ $\text{Re}(s) \ge 0$)، يحاول المؤلفون تمديد صيغة فلوكيت (Floquet ansatz) $y(t) = e^{\lambda t} p(t)$ للأنظمة ذات الرتب الكسرية.
4. مصفوفة هيل الكسرية: من خلال تعويض صيغة فلوكيت في LTP-FODE واستخدام مفككات متسلسلة فوريه للمعاملات الدورية، يستنتج المؤلفون مسألة قيم ذاتية غير خطية ذات أبعاد لانهائية. يؤدي هذا إلى تعريف "مصفوفة هيل الكسرية"، $H^{(\alpha)}_\infty(\lambda)$، حيث تتضمن الكتل القطرية حدوداً من الشكل $(\lambda + ik\omega)^\alpha I - J_k$.
5. تحليل الأنظمة الفرعية ذات الزمن الثابت (LTI): قبل معالجة الحالة الدورية، يحلل المؤلفون معادلات LTI-FODEs ذات الذاكرة اللانهائية. ويقارنون حلول نظام Liouville-Weyl بنظام Caputo (ذات الحد السفلي صفر) لفهم سلوك الحلول المتلاشية.

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

• إطار العمل للحلول الدورية: تؤسس الورقة مسألة قيمة أولية صارمة لصياغة معادلات Liouville-Weyl ذات الرتب الكسرية التي تسمح بوجود حلول دورية تماماً، متجاوزة بذلك القيد الذي تفرضه مشتقات Caputo.
• توسيع نظرية فلوكيت: يستنتج المؤلفون نهج مصفوفة هيل الكسرية لـ LTP-FODEs. ويثبتون أنه إذا وجد حل على شكل $y(t) = e^{\lambda t} p(t)$ مع $\text{Re}(\lambda) \ge 0$، فيجب أن يحقق مسألة القيم الذاتية غير الخطية $H^{(\alpha)}_\infty(\lambda)p_\infty = 0$.
• القصور في الحلول المتلاشية: النتيجة المركزية هي أن منهج هيل الكسري المقترح (وصيغة الأسية بشكل عام) يمكنه فقط تقييم النمو الأسي (عدم الاستقرار). بينما يفشل في رصد التلاشي الجبري.
  • في حالة LTI، يوضح المؤلفون أنه بينما ترتبط الحلول المتباعدة بقيم ذاتية ذات أجزاء حقيقية موجبة، فإن الحلول المتقاربة (التي تتلاشى جبرياً، $O(t^{-\alpha})$) لا ترتبط بقيم ذاتية ذات أجزاء حقيقية سالبة في المعادلة المميزة المستمدة من الصيغة الأسية.
  • يُعزى هذا القصور إلى حقيقة أن حلول معادلات Liouville-Weyl ذات المعاملات المستقرة تظهر عادةً تلاشياً جبرياً بدلاً من التلاشي الأسي. الصيغة الأسية غير صالحة لهذه الأنماط المتلاشية لأنها لا تنتمي إلى الفضاء الدالي المطلوب لتعريف مشتقة Liouville-Weyl بالطريقة القياسية.
• التحقق العددي: يبرهن المؤلفون على منهجيتهم عددياً على أنظمة قياسية وثنائية الأبعاد (بما في ذلك نظام من نوع Mathieu كسري). ويظهرون أن مصفوفة هيل الكسرية المبتورة تحدد بنجاح القيم الذاتية ذات الأجزاء الحقيقية الموجبة التي تقابل المسارات غير المستقرة، مما يؤكد فائدة الطريقة في اكتشاف عدم الاستقرار.
• خصائص مصفوفة هيل الكسرية: تثبت الورقة أن القيم الذاتية لمصفوفة هيل الكسرية المبتورة تأتي في أزواج مترافقة معقدة وتقع ضمن مناطق تشبه مناطق جيرشورين (Gershgorin-like regions) في المستوى المركب.

الأهمية والادعاءات
تدعي الورقة أن الميزة الأساسية لهذا العمل هي تقديم رؤية مفادها أن الحلول الدورية ذات التلاشي الأسي لا توجد في LTP-FODEs، في حين يمكن ملاحظة وتقييم الحلول ذات النمو الأسي عبر نهج من نوع هيل.

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

وتختتم الورقة بتواضع، مشيرة إلى أن غياب القيم الذاتية ذات الأجزاء الحقيقية غير السالبة في مسألة هيل الكسرية يعد مؤشراً قوياً على الاستقرار، ولكنه يتطلب مزيداً من البحث ليُعتبر حاسماً. علاوة على ذلك، حدد المؤلفون بناء مؤثر مونو درومي (monodromy operator) لـ LTP-FODEs كتحدٍ مفتوح، مقترحين أن طيف هذا المؤثر سيكون ضرورياً لتوصيف الاستقرار بشكل كامل، بما في ذلك أنماط التلاشي الجبري التي يفتقدها منهج هيل الحالي.