Relation between two Sinc-collocation methods for Volterra integral equations of the second kind and further improvement

فيما يلي ملخص تقني مفصل للورقة البحثية بعنوان: "العلاقة بين طريقتين من طرق التوفيق (Sinc-collocation) لمعادلات فولتيرا التكاملية من النوع الثاني وتحسينهما".

1. بيان المشكلة

تتناول الورقة البحثية الحل العددي لـ معادلات فولتيرا التكاملية من النوع الثاني، والتي تُعطى بالصيغة:
$$u(t) - \int_a^t k(t, s)u(s) \, ds = g(t), \quad a \leq t \leq b$$
حيث $k(t, s)$ هو النواة (kernel) و $g(t)$ هي دالة معلومة. والهدف هو تحديد الحل المجهول $u(t)$ بدقة وكفاءة عاليتين.

تركز الدراسة على طرق Sinc العددية، المعروفة بمعدلات تقارب أسية للدوال التحليلية. وتحديداً، تبحث الدراسة في العلاقة وفروق الأداء بين ثلاثة نهج قائمة تعتمد على Sinc:

1. طريقة Sinc-Nyström: تعتمد على التكامل غير المحدد لـ Sinc.
2. طريقة Stenger للتوفيق بـ Sinc (Stenger's Sinc-collocation Method): مشتقة أصلاً للنوى ذات المتغير الواحد ولكنها قابلة للتطبيق في الحالات العامة.
3. طريقة Rashidinia–Zarebnia للتوفيق بـ Sinc: نهج توفيق قياسي مع معالجة محددة للشروط الحدية.

2. المنهجية

تستخدم الورقة تحليلاً نظرياً مقارناً وتقترح طريقة هجينة جديدة. وتشمل المنهجيات الأساسية ما يلي:

• تحويلات المتغيرات: لتعيين الفترة المحدودة $[a, b]$ إلى الخط الحقيقي اللانهائي $\mathbb{R}$ (حيث تُعرف دالات Sinc)، تستخدم الورقة تحويلين:
  • تحويل الأس المنفرد (SE): $t = \psi_{SE}(x) = \frac{b-a}{2}\tanh(\frac{x}{2}) + \frac{b+a}{2}$.
  • تحويل الأس المزدوج (DE): $t = \psi_{DE}(x) = \frac{b-a}{2}\tanh(\frac{\pi}{2}\sinh x) + \frac{b+a}{2}$.
• تقريب Sinc والتكامل غير المحدد: تعتمد الطرق على تقريب الدوال وتكاملاتها باستخدام دالات أساس Sinc $S(j, h)(x)$ ودالة تكامل الجيب (Sine Integral function) $Si(x)$.
• تقريب Sinc المعمم: للتعامل مع الشروط الحدية غير الصفرية ($u(a) \neq 0$ أو $u(b) \neq 0$)، يستخدم المؤلفون مؤثر تقريب معمّم ($P_N$) يدمج دالات حدية خطية $\omega_a(t)$ و $\omega_b(t)$.

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

أ. العلاقة النظرية بين الطرق القائمة

تحلل الورقة بدقة العلاقة بين طريقة Stenger ($v_N^{SE}$) و طريقة Rashidinia–Zarebnia ($w_N^{RZ}$).

• التطابق عند نقاط التوفيق: تثبت الدراسة أنه بينما تولد الطريقتان دالتين مستمرتين مختلفتين عالمياً، إلا أنهما تعطيان قِيماً متطابقة عند نقاط التوفيق ($t_i$).
• عدم التكافؤ: رغم تطابقهما عند العقد، فإن الطريقتين ليستا متكافئتين بشكل عام لأنهما تستخدمان مؤثرات تقريب مختلفة ($P_N^{SE}$ مقابل $P_N^{RZ}$) للقيام بعملية الاستكمال (interpolation) بين العقد.
• معدلات التقارب: ثبت أن كلتا الطريقتين تحققان نفس معدل التقارب ذو الجذر الأسي:
  $$O(\sqrt{N} e^{-\sqrt{\pi d \alpha N}})$$
  حيث $d$ هو مجال التحليل و $\alpha$ هو رتبة استمرارية هولدر (Hölder continuity).

ب. تحليل التنفيذ العملي

• التعقيد: تتطلب طريقة Rashidinia–Zarebnia حل نظام خطي أكثر تعقيداً يتضمن الحدود الحدية صراحة في بنية المصفوفة. أما طريقة Stenger، فتُشتق عن طريق حل نظام Sinc-Nyström أولاً ثم تطبيق تقريب Sinc المعمم، مما يجعلها أبسط في التنفيذ.
• الكفاءة مقابل Nyström: بينما تمتلك طريقة Sinc-Nyström معدل تقارب نظري أفضل قليلاً (بدون عامل $\sqrt{N}$)، إلا أنها تتطلب تقييم دالة تكامل الجيب (Si) ضمن دالات الأساس، مما يسبب تكلفة حسابية عالية. وتتجنب طريقة Stenger للتوفيق بـ Sinc ذلك، مما يوفر توازناً أفضل بين الدقة ووقت الحساب.

ج. التحسين المقترح: طريقة DE-Sinc-collocation

تقترح الورقة طريقة جديدة: طريقة DE-Scollocation.

• الآلية: تجمع بين إجراء التوفيق الخاص بـ Stenger وتحويل الأس المزدوج (DE) بدلاً من تحويل tanh (SE) القياسي.
• التبرير النظري: يثبت المؤلفون أن هذه الطريقة الجديدة تحقق معدل تقارب أعلى بكثير:
  $$O(e^{-\pi d N / \log(2dN/\alpha)})$$
  هذا المعدل هو "شبه أسي" وهو يتفوق بشكل كبير على معدل التقارب ذو الجذر الأسي للطرق القائمة على SE.

4. النتائج

• الإثباتات النظرية:
  • أثبتت أن طريقتي Stenger و Rashidinia–Zarebnia تتطابقان عند نقاط التوفيق ولكن تختلفان عالمياً.
  • وضعت حدود خطأ صارمة لكلتا الطريقتين القائمتين على SE تحت فرضيات استمرارية هولدر.
  • اشتققت حد الخطأ لطريقة التوفيق الجديدة القائمة على DE، مما يظهر تفوقها على الطرق القائمة على SE.
• التجارب العددية:
  • تم الاختبار على معادلتين مرجعيتين من معادلات فولتيرا التكاملية ذات حلول معروفة ($u(t) = e^{-t^2}$ و $u(t) = 1+\sqrt{t}$).
  • الدقة: أظهرت طريقة DE-Sinc-collocation أفضل أداء، حيث حققت أعلى دقة لعدد معين من العقد $N$.
  • الكفاءة: عند رسم الخطأ مقابل وقت الحساب، تفوقت طريقة DE-Sinc-collocation بشكل كبير على طرق Sinc-Nystrom. كانت طرق Nyström أبطأ بسبب التكلفة الحسابية لدالة تكامل الجيب، رغم مزايا التقارب النظرية لها.
  • قدمت طريقة Stenger للتوفيق بـ SE أداءً مماثلاً تقريباً لطريقة Rashidinia–Zarebnia، لكنها كانت أسهل في التنفيذ.

5. الأهمية

1. توضيح الأدبيات: تحل الورقة البحثية الغموض المتعلق بالعلاقة بين اثنتين من أبرز طرق التوفيق بـ Sinc، حيث أثبتت أنهما مختلفتان ولكن تشتركان في رتبة التقارب والقيم عند العقد.
2. توجيهات التنفيذ: تقدم توصية واضحة للممارسين: طريقة Stenger هي المفضلة على طريقة Rashidinia–Zarebnia نظراً لبساطتها، وكلتاهما مفضلة على طريقة Sinc-Nyström للتطبيقات العامة بسبب انخفاض العبء الحسابي.
3. التقدم المنهجي: يمثل تقديم طريقة DE-Sinc-collocation تحسيناً كبيراً. فمن خلال الجمع بين البساطة الهيكلية لإجراء التوفيق الخاص بـ Stenger وخصائص التقارب المتفوقة لتحويل الأس المزدوج (DE)، حقق المؤلفون طريقة تتسم بأنها فعالة حسابياً وعالية الدقة.
4. التعميم: تبرر الدراسة نظرياً استخدام طريقة Stager لمعادلات $k(t, s)$ العامة (دوال بمتغيرين)، مما يوسع نطاق تطبيقها خارج صياغتها الأصلية للمتغير الواحد.

ختاماً، تثبت الورقة أن طريقة DE-Sinc-collocation هي الخيار الأمثل من بين التقنيات التي تمت مناقشتها، حيث توفر حلاً قوياً وسريعاً وعالي الدقة لمعادلات فولتيرا التكاملية من النوع الثاني.