Determining the Rolle function in Hermite interpolatory approximation by solving an appropriate differential equation

ملخص تقني: تحديد دالة رول في التقريب الاستيفائي لـ هيرميت

بيان المشكلة
تتناول الورقة البحثية التحدي المتمثل في تحسين دقة التقريبات الاستيفائية لهيرميت (Hermite). بينما يقوم كثير حدود هيرميت القياسي $H_{2n+1}(x)$ باستيفاء دالة $f(x)$ ومشتقتها $f'(x)$ عند $n+1$ من العقد، فإن الخطأ النقطي $\Delta(x | H_{2n+1})$ محكوم بحد يتضمن المشتقة رقم $(2n+2)$ للدالة $f$ عند نقطة مجهولة $\xi(x)$، تُعرف باسم دالة رول (Rolle function). تكمن المشكلة الجوهرية في أن $\xi(x)$ غير معروفة بشكل عام، مما يمنع الحساب المباشر أو التصحيح لترم الخطأ. يسعى المؤلف إلى إيجاد طريقة لتحديد $\xi(x)$ عددياً لتسهيل الحصول على تقريب أكثر دقة لـ $f(x)$.

المنهجية
توسع المنهجية المقترحة تقنية طُبقت سابقاً على استيفاء لاجرانج (Lagrange) لتشمل حالة هيرميت. يتضمن هذا النهج الخطوات التالية:

1. اشتقاق معادلة تفاضلية: انطلاقاً من صيغة الخطأ القياسية لاستيفاء هيرميت، قام المؤلف باشتقاق حد الخطأ بالنسبة لـ $x$. ومن خلال تطبيق قاعدة السلسلة وعزل مشتق دالة رول، $\frac{d\xi}{dx}$، تم اشتقاق معادلة تفاضلية عادية (ODE) من الدرجة الأولى. تربط هذه المعادلة التفاضلية بين $\frac{d\xi}{dx}$ والدوال المعروفة $f(x)$ و $H_{2n+1}(x)$ ومشتقاتها، وبين $\xi(x)$ المجهولة.
2. صياغة مسألة القيمة الابتدائية: لحل المعادلة التفاضلية، يلزم وجود قيمة ابتدائية $\xi_z = \xi(x_z)$ عند نقطة $x_z$ قريبة من (ولكن ليست عند) عقدة استيفاء. وبما أن الخطأ يكون صفراً عند العقد نفسها، فلا يمكن تحديد $\xi$ هناك. اختار المؤلف نقطة $x_z$ بالقرب من عقدة وقام بحل معادلة الخطأ الجبرية عددياً لإيجاد قيمة ابتدائية صالحة لـ $\xi_z$.
3. الحل العددي: تم حل المعادلة التفاضلية المستمدة عددياً باستخدام طريقة رونج-كوتا (Runge-Kutta) من الدرجة السابعة. تشير الورقة إلى إمكانية وجود حلول رياضية متعددة للقيمة الابتدائية؛ ويتم اختيار الحل الصحيح بناءً على القيد الذي يقضي بأن يجب أن تقع $\xi(x)$ تماماً داخل الفترة المعرفة بعقد الاستيفاء ($x_0 < \xi(x) < x_n$).
4. تصحيح الخطأ: بمجرد تحديد $\xi(x)$، يصبح حد $f^{(2n+2)}(\xi(x))$ في صيغة الخطأ معلوماً. قام المؤلف بتقريب هذا الحد باستخدام كثير حدود $H_\xi(x)$ (عبر مطابقة المربعات الصغرى أو الـ splines التكعيبية). يُضرب هذا الكثيرة الحدود في العامل العقدي $\prod (x-x_k)^2$ ويُقسم على $(2n+2)!$ لتشكيل كثير حدود الخطأ $E(x)$.
5. التقريب المنقح: يتم بناء التقريب المحسن النهائي عن طريق إضافة كثير حدود الخطأ إلى كثير حدود هيرميت الأصلي: $f(x) \approx H_{2n+1}(x) + E(x)$.

النتائج الرئيسية
تم اختبار المنههجية على الدالة $f(x) = e^x \sin x$ فوق العقد $\{0, 3\pi/2\}$ (حيث $n=1$).

• تحديد دالة رول: نجح الحل العددي في توليد دالة رول $\xi(x)$ استوفت قيود الفترة. وقد تم التحقق من الحل العددي للمعادلة التفاضلية عبر مقارنة الطرفين الأيسر والأيمن لمعادلة الخطأ، حيث أظهرت النتائج أنهما متطابقان جوهرياً.
• تحسين الدقة:
  • أظهر كثير حدود هيرميت الأصلي $H_3(x)$ خطأ تقريب أقصى قدره 7.04.
  • من خلال إضافة كثير حدود الخطأ $E(x)$ المشتق من مطابقة المربعات الصغرى من الدرجة 11، انخفض الخطأ الأقصى إلى $7.3 \times 10^{-8}$.
  • استخدام الـ spline التكعيبي لتوليد تقريب لحد المشتقة قلل الخطأ الأقصى بشكل أكبر إلى حوالي $10^{-12}$.
• الحفاظ على خصائص الاستيفاء: تتحقق الورقة من أن التقريب المصحح $H_{2n+1}(x) + E(x)$ يحتفظ بخصائص استيفاء هيرميت الأصلية، حيث يطابق كلاً من قيم الدالة وقيم المشتقة عند العقد.
• التكامل العددي: توضح الورقة بإيجاز أن تكامل كثير الحدود المصحح عبر الفترة يعطي نتيجة أكثر دقة بكثير ($\sim 3 \times 10^{-12}$ خطأ) مقارنة بتكامل كثير حدود هيرميت الأصلي ($\sim 0.7$ خطأ).

الأهمية والادعاءات
تدعي الورقة أن تحديد دالة رول عبر مسألة القيمة الابتدائية يسمح بالتقريب الصريح للخطأ النقطي في استيفاء هيرميت. تكمن الأهمية الرئيسية في القدرة على بناء تقريب جديد لكثير حدود يكون أكثر دقة بعدة مراتب من كثير حدود هيرميت القياسي، مع الحفاظ بدقة على شروط الاستيفاء الأصلية.

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