Improvement of conformal maps combined with the Sinc approximation for derivatives over infinite intervals

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

تتناول الورقة البحثية الاشتقاق العددي للدوال المعرفة على فترات غير منتهية، وتحديداً $(0, \infty)$ و $(-\infty, \infty)$. يعتمد النهج القياسي على تقريب سينك (Sinc approximation)، الذي يجمع بين أساس دالة سينك مع خريطة مطابقة (conformal map) لتحويل المجال غير المنتهي إلى مجال منتهٍ أو أكثر سهولة في التعامل معه.

لقد وضع ف. ستينجر (F. Stenger) سابقاً صيغ تقريب فعالة للمشتقات من خلال الجمع بين تقريب سينك وخرائط مطابقة محددة ($\psi_i$). تحقق هذه الصيغ تقارباً ذا أس جذري (root-exponential convergence) من الرتبة $O(\exp(-c\sqrt{n}))$، حيث $n$ هو عدد نقاط التجزئة. ومع ذلك، فإن معدل التقارب يعتمد بشكل كبير على المعاملات $d$ (عرض شريط التحليل في المستوى المركب) و $\mu$ (المتعلقة بمعدل اضمحلال الدالة).

المشكلة المحددة التي تستهدفها هذه الدراسة هي معدل التقارب دون الأمثل في حالتين محددتين حددهما ستينجر:

• الحالة 2: الفترة $(0, \infty)$ حيث تضمحل الدالة أسياً عندما $t \to \infty$. خريطة ستينجر: $\psi_2(x) = \text{arsinh}(e^x)$.
• الحالة 4: الفترة $(-\infty, \infty)$ حيث تضمحل الدالة جبرياً عندما $t \to -\infty$ وأسياً عندما $t \to \infty$. خريطة ستينجر: $\psi_4(x) = \sinh(\log(\text{arsinh}(e^x)))$.

في هاتين الحالتين، تحد الخرائط المطابقة الموجودة من قيم $d$ و $\mu$ التي يمكن تحقيقها، مما يضع سقفاً لسرعة التقارب.

2. المنهجية

يقترح المؤلفون نهجاً جديداً يتمثل في استبدال خرائط ستينجر المطابقة بخرائط جديدة محسنة مع الاحتفاظ بإطار عمل تقريب سينك.

الخرائط المطابقة المقترحة

• للحالة 2 ($t \in (0, \infty)$):
  استبدال $\psi_2(x) = \text{arsinh}(e^x)$ بـ:
  $$\phi_2(x) = \log(1 + e^x)$$
• للحالة 4 ($t \in (-\infty, \infty)$):
  استبدال $\psi_4(x) = \sinh(\log(\text{arsinh}(e^x)))$ بـ:
  $$\phi_4(x) = 2 \sinh(\log(\log(1 + e^x)))$$

الإطار النظري

يشتق المؤلفون صيغ تقريب جديدة للمشتقة رقم $m$ للدالة $f^{(m)}(t)$:
$$f^{(m)}(t) \approx \sum_{k=-M}^{N} \frac{f(\phi_i(kh))}{g_i(\phi_i(kh))} \left( \frac{d}{dt} \right)^m \left\{ g_i(t) S(k, h)(\phi_i^{-1}(t)) \right\}$$
حيث $g_i(t)$ هي دوال وزن يتم اختيارها لقمع التباعد عند النهايات، و $S(k, h)$ هي دالة سينك.

يتضمن جوهر المنهجية تحليلاً نظرياً دقيقاً للخطأ لإثبات أن هذه الخرائط الجديدة تسمح بقيم أكبر للمعاملات $d$ و $\mu$ في حد التقارب. تتضمن استراتيجية الإثبات ما يلي:

1. خطأ التجزئة: استخدام فضاء الدالة $H_1(D_d)$ وصيغة فآ د برونو (Faà di Bruno) لتحديد حدود مشتقات الدوال المركبة التي تتضمن الخرائط العكسية الجديدة.
2. خطأ القطع: تحليل اضمحلال الدالة $F(\zeta) = f(\phi(\zeta))/g(\phi(\zeta))$ في الشريط المركب لإنشاء حدود لذيل المتسلسلة.
3. التركيب: دمج هذه الحدود لاستخلاص تقدير الخطأ النهائي.

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

1. خرائط مطابقة جديدة: تقديم $\phi_2$ و $\phi_4$ كبدائل متفوقة لخرائط ستينجر $\psi_2$ و $\psi_4$ لتقريب المشتقة.
2. نظريات تقارب محسنة:
   • النظرية 3.1 (الحالة 2): تثبت التقارب مع $d < \pi$ (مقارنة بـ $d < \pi/2$ لستينجر).
   • النظرية 3.2 (الحالة 4): تثبت التقارب مع $d < \pi$ وتسمح بقيمة أكبر لـ $\mu$ (تحديداً، يمكن مضاعفة معامل الاضمحلال $\beta$ فعلياً مقارنة بصياغة ستينجر بسبب بنية الخريطة الجديدة).
   • يظل حد الخطأ $O(n^{(m+1)/2} \exp(-\sqrt{\pi d \mu n}))$, ولكن زيادة $d$ و $\mu$ تؤدي إلى تصغير الثابت في الأس بشكل كبير، مما يسرع التقارب.
3. إثباتات دقيقة: إثباتات مفصلة باستخدام التحليل المركب، وتحديداً التعامل مع التفردات وخصائص الاضمحلال للخرائط العكسية $\phi_i^{-1}(t)$ ومشتقاتها.

4. النتائج

أجرى المؤلفون تجارب عددية باستخدام حسابات الدقة المزدوجة (ودقة 20 رقماً لحسابات مشتقة محددة لتجنب الضجيج) للتحقق من النظرية.

• حالة الاختبار 1 (الحالة 2): $f(t) = \sqrt{\frac{t}{1+t}} e^{-t}(1-e^{-t})^2$ على $(0, \infty)$.
  • المعاملات: ستينجر ($d=1.57, \mu=0.5$) مقابل المحسنة ($d=3.14, \mu=0.5$).
  • النتيجة: تقاربت الصيغة المحسنة بشكل أسرع بكثير. انخفض الخطأ بسرعة أكبر مع زيادة عدد تقييمات الدالة ($M+N+1$).
• حالة الاختبار 2 (الحالة 4): $f(t) = \frac{1}{(4+t^2)(1+e^{\pi t/2})}$ على $(-\infty, \infty)$.
  • المعاملات: ستينجر ($d=1.57, \mu \approx 0.785$) مقابل المحسنة ($d=2.07, \mu \approx 1.57$).
  • النتيجة: أظهرت الصيغة المحسنة تقارباً أسرع. بالإضافة إلى ذلك، بدا ملف الخطأ لصيغة ستينجر "متعرجاً" (أقل سلاسة)، بينما أظهرت الصيغة المحسنة ملف تقارب سلس ورتيب.

في كلتا الحالتين، تفوقت الصيغة المحسنة (1.6) باستمرار على صيغة ستينجر (1.5)، مما يؤكد التنبؤ النظري بأن قيم $d$ و $\mu$ الأكبر تؤدي إلى أداء أفضل.

5. الأهمية

• الكفاءة: يحقق النهج المقترح نفس مستوى الدقة بعدد أقل من تقييمات الدالة مقارنة بطريقة ستينجر الأصلية. وهذا أمر بالغ الأهمية لحل المعادلات التفاضلية حيث تكون المشتقات من الرتب العليا مطلوبة بشكل متكرر.
• التقدم النظري: توضح الورقة أن اختيار الخريطة المطابقة أمر بالغ الأهمية حتى ضمن فئة تحويل "الأس المنفرد" (single-exponential). وهي توفر مخططاً لتحسين هذه الخرائط من خلال تحليل عرض شريط التحليل ($d$) ومعدلات الاضمحلال ($\mu$).
• الاتجاهات المستقبلية: يشير المؤلفون إلى أنه بينما تم تحسين تحويلات الأس المنفرد، فقد يكون من الممكن تحقيق مكاسب إضافية من خلال الانتقال إلى تحويلات الأس المزدوج (double-exponential transformations) (المعروفة بتسريع التقارب في تقريبات التكامل). ويشيرون إلى أن هذا هو موضوع بحثهم المستمر.

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