عرض مصدر الاستكمال الرياضي بفروق نيوتن

{{لا مصدر|تاريخ=ديسمبر 2018}}
{{يتيمة|تاريخ=ديسمبر 2014}}

اكتشف طريقة الاستكمال العالم الهندي القديم ''آرياباتا'' حينما كان يقيس أطوال الأقواس المناظرة للزوايا في [[دائرة|الدائرة]] حيث وصل الي الاستكمال من الدرجة الثانية .

ويعتبر '''الاستكمال''' الرياضي الموضوع الأبرز في [[تحليل عددي|التحليل العددي]] ، إذ يشكل قلب ونواة التحليل العددي الكلاسيكي. وذلك لسبيبن رئيسيين، السبب الأول يعود لحاجتنا المستمرة في البحث عن قيمة لدالة من بيانات مجدولة أثناء معظم المسائل الحسابية، أما في تلك المسائل والنقاشات الغير مجدولة فلكي نجد قيمة للدالة عند واحدة أو أكثر من النقاط الغير مدرجة في جدول البيانات، فلابد لنا من أن نستكمل تلك الدالة ونستخدم طرق '''الاستكمال'''. والأكثر من ذلك أن الحاجة '''للاستكمال''' تكمن في كون أن البيانات المجدولة التي تُعطى إلينا في معظم المسائل تكون لها من الدقة العالية الشيء الكثير، حتى وإن كانت بيانات محدودة. لذلك قدم التحليل العددي الكلاسيكي مجموعة متطورة جدا من الطرق المختلفة للاستكمال الرياضي. أما بالنسبة للسبب الثاني لأهمية '''الاستكمال''' الرياضي فيعود لكون أن معظم الطرق العددية الكلاسيكية في شتى القطاعات قد تم استناتجها واشتقاقها من طرق '''الاستكمال'''({{وإو|عر=Interpolation%20|تر=Interpolation%20|لغ=en|نص=Interpolation Formulas}})، فتلك الطرق العددية المستخدمة في [[مشتق (رياضيات)|الاشتقاقات]]، [[تكامل]]ات [[معادلة تفاضلية عادية|المعادلات التفاضلية العادية]]، المعادلات التربيعية، وغيرها من قطاعات التحليل العددي الكلاسيكي قد طُورت واشتُـقّـت مباشرةً انطلاقًا من طرق '''الاستكمال''' الرياضي. بالرغم من أن الطرق المستخدمة في التحليل العددي الحديث لا تعتمد ذاك الاعتماد الكبير على طرق '''الاستكمال'''؛ لوجود طرق أخرى اشتُقت منها إلا أن هذا لا يتعارض مع الدور الكبير والفائدة الجمة '''للاستكمال''' وطرق '''الاستكمال'''.

== طرق الاستكمال ==
* '''الاستكمال''' الرياضي له العديد من الطرق، منها: 
* 	الاستكمال بطريقة لاجرانج ([[متعدد حدود لاغرانج|Lagrangian Interpolation]]).
* 	الاستكمال بطريقة لاجرانج على الفترات المتساوية (Lagrangian Interpolation at equal intervals).
* 	الاستكمال بالفروق المنتهية ([[فرق محدود|Finite Differences Formulas]]).
* 	'''الاستكمال بطريقة نيوتن الأمامية والخلفية''' (Newton’s Forward/Backward Formulas) .
* 	الاستكمال بطريقة جاوس الأمامية والخلفية ([https://web.archive.org/web/20180614123330/http://www.nptel.ac.in/courses/122104018/node114.html Gauss Forward/Backward Formulas]).
* 	الاستكمال بطريقة سترلينغ ([https://web.archive.org/web/20180614123330/http://www.nptel.ac.in/courses/122104018/node114.html Stirling’s Formula]).
* 	الاستكمال بطريقة بيزييل ([http://www.abecedarical.com/zenosamples/zs_bessel_interpolation.html Bessel’s Formula]).
ونظرًا لأهمية طرق '''الاستكمال''' باستخدام الفروق في الكثير من المسائل والتطبيقات الحسابية، خصصنا هذا الموضوع لمناقشة إحدى أهم هذه الطرق، ألا وهي '''الاستكمال بفروق نيوتن الأمامية والخلفية'''.

== تقديم رياضي للاستكمال ==
نفترض أن لدينا [[دالة]] <math>f(x) </math> معرفة وقابلة للاشتقاق، لها قيم معينة عند مجموعة محددة من النقاط، حيث أن هذه المجموعة من النقاط تسمى بالنقاط المجدولة (tabular points) التي يتم الاعتماد عليها كثيرًا أثناء الاستكمال. الآن الهدف من '''الاستكمال''' يكمن في الحصول على قيمة تقريبة للدالة عند واحدة أو أكثر من النقاط الأخرى الغير مدرجة في الجدول، بحيث أن هذا التقريب سيحد من الخطأ الناتج بين القيمة الفعلية والصحيحة للدالة عند النقطة والقيمة التقريبية لها بعد '''الاستكمال'''.                            وصلتنا الرياضية التالية هي الحصول على تقريب [[دالة|للدالة]] <math>f(x) </math> من خلال اعتبار أن النقاط والقيم المجدولة المعطاة تمثل الدالة <math>y(x) </math> التي تعطي قيمًا متساوية مع قيم [[دالة|الدالة]] <math>f(x) </math> ( وكذلك بالنسبة لقيم [[مشتق (رياضيات)|مشتقتها]] إن وجدت ).

 الدالة '''المستكملة''' تُعطى بالعلاقة : 
  (1)	 <math>f(x)= \sum_{j=1}^n l_j(x)f(a_j)+E(x) = y(x) + E(x) </math>

والعلاقة العامة أكثر '''للاستكمال''' هي: 
  (2)     <math>f(x)= \sum_{j=1}^n\sum_{i=0}^{m_j} A_{ij}(x)f^{(i)}(a_j)+E(x) </math>
الهدف هو تحديد <math>l_j(x) </math> بحيث أن                                                               
    (3)                    <math>E(a_j)=0,\qquad    j=1, ...,n</math>
تكون مستقلة عن [[دالة|الدالة]]   <math>f(x) </math> . و بصورة عامة، لأي نقطة غير مجدولة   <math>E(x)\ne 0 </math>  نحن نسعى بطرق '''الاستكمال''' إلى تحديد  <math>l_j(x) </math>   بحيث أن  (3) تتحقق، ولإيجاد تمثيل لـ  <math>E(x) </math>  بحيث يمكننا تقدير الخطأ أو على الأقل حده للقيم  <math>x\ne a_j,\;j=1, ...,n </math>

== الاستكمال بطريقة نيوتن ==
=== صيغة نيوتن الأمامية ===

'''تُعطى صيغة نيوتن الأمامية بالعلاقة التالية:'''

 (4) <math>y(m)=f_o + (m)_1 \triangle f_o + (m)_2 \triangle^2 f_o +...+ (m)_n \triangle^n f_o =\sum_{j=0}^n (m)_j \triangle^j f_o   </math>
حيث أن هذه الصيغة تم توليدها كغيرها من الصيغ باستخدام ما يعرف بمخطط لوزيينغ ([https://books.google.com.sa/books?id=Z2owE_0LQukC&pg=PA304&lpg=PA304&dq=Lozenge+diagram&source=bl&ots=2y3KXuSKwU&sig=8hotvaHUwX2SoO2Mjvx5SXWw3Wo&hl=en&sa=X&ei=TbaSVKChOMqr7AbB0YGoDg&ved=0CB4Q6AEwAA#v=onepage&q=Lozenge%20diagram&f=false%20Lozenge%20diagram The lozenge diagram])،حيث عمليات التوليد هذه من المخطط تؤكد أن أي صيغة للفروق تحوي n من الحدود، هي مكافئة جبريًا لصيغة لاجرانج على الفترات المتساوية. وسنقدم برهانًا يوضح كيف أن '''صيغة نيوتن الأمامية''' تحقق هذه النتيجة. 
 
كما أن برهان النتيجة السابقة يتطلب توضيح أن:
#  على الأقل واحدة من صيغ الفروق المنتهية '''للاستكمال''' لها هذه الخاصية، وسنعنى بإثبات تحقق ذلك على '''صيغة نيوتن الأمامية''' في الفقرة التالية.
#  جميع الصيغ التي تنتهي بنفس الفروق تكون متكافئة جبريًا، بغض النظر عن المسارات التي تم المرور عليه اأثناء العمل على مخطط لوزيينغ.

=== البــرهان ===
نريد أن نثبت أن العلاقة (4) تكافئ جبريًا [[متعدد حدود لاغرانج|صيغة لاجرانج '''الاستكمالية''']] على الفترات المتساوية للنقاط <math>a_o, a_1, ..., a_n </math> .
من كون <math>(m)_n </math> [[متعددة الحدود|كثيرة حدود]] من الدرجة  n في m، فيكفي إثبات أن <math>y(i) </math>  في العلاقة (4) تساوي <math>f_i ,i=0, ...,n </math> 
ومن ثم فإن <math>y(m) </math> ستكون كثيرة حدود وحيدة من الدرجة n، والمارة بال n+1 نقطة للدالة <math>f_i </math>.

باستخدام (4) والعلاقة <math>\triangle^i f(x)=\sum_{k=1}^j (-1)^{j-k} {j \choose K} f(x+kh) </math>
 
نجد أن

<math> y(i)=\sum_{j=0}^n (i)_j \triangle^j f_o = \sum_{j=0}^n\sum_{k=0}^j (-1)^{j-k} (i)_j {j \choose k} f_k </math>
 (5) <math>= \sum_{k=0}^n\sum_{j=k}^n (-1)^{j-k} (i)_j {j \choose k}f_k \quad i=0, ...,n </math>

عوامل <math>f_r </math> في <math>y(i) </math> تُعطى بالعلاقة:

 (6) <math>\sum_{j=r}^i (-1)^{j-r} (i)_j {j \choose r} </math>

لكل <math> r>i </math> فإن هذا المعامل مساوي للصفر؛ من كون  <math>(i)_j=0 </math> عندما <math>i<j </math>. عندما <math>r=i </math> فإن الحد الغير صفري الوحيد في (6) لكل  <math>j=r </math> يساوي 1. وعندما  <math>r<i </math> فإن العلاقة (6) يمكن كتابتها على الشكل:

 (7)   <math>\sum_{j=r}^i (-1)^{j-r} {i \choose j} {j \choose r} </math>

ومن ثم نجد أن الطرف الأيمن من العلاقة (5) ما هو إلا <math>f_i </math> ،وبذلك يتم البرهان.

باستخدام مخطط لوزيينغ ([https://books.google.com.sa/books?id=Z2owE_0LQukC&pg=PA304&lpg=PA304&dq=Lozenge+diagram&source=bl&ots=2y3KXuSKwU&sig=8hotvaHUwX2SoO2Mjvx5SXWw3Wo&hl=en&sa=X&ei=TbaSVKChOMqr7AbB0YGoDg&ved=0CB4Q6AEwAA#v=onepage&q=Lozenge%20diagram&f=false Lozenge diagram]) يمكن توليد العديد من صيغ '''الاستكمال'''، ومنها '''صيغة نيوتن الخلفية'''.

=== صيغة نيوتن الخلفية ===

في مخطط لوزيينغ بالابتداء من <math>f_o </math> والتحرك على طول المستقيمات سالبة الميل وباتجاه اليمين، نتوصل إلى أن :

 (8) <math>y(m)= f_o +(m)_1 \triangle f_{-1} + (m+1)_2 \triangle^2 f_{-2} + ... + (m+n-1)_n \triangle^n f_{-n}</math>
حيث هذه العلاقة تكافئ صيغة لاجرانج ([[متعدد حدود لاغرانج|Lagrangian formula]]) باستخدام النقاط  <math> a_o, a_{-1}, ..., a_{-n} </math> .

== الاستكمال بطريقة نيوتن وتطبيقات الحياة ==

تعتبر [[رياضيات|الرياضيات]] من أكثر العلوم صلة بالتطبيقات العملية الحياتية، لذلك فإنه ليس من الغريب أن يكون '''للاستكمال''' و'''الاستكمال بطريقة نيوتن''' النصيب الكبير من هذه التطبيقات. 
من أهم القطاعات التي يستخدم فيها '''الاستكمال''' بكثرة، نذكر منها ما يلي:
* قطاعات [[إحصاء|الإحصاء]] والتنبؤ الرياضي.
* القطاعات الهندسية، وخصوصَا [[فيزياء هندسية|الهندسة الفيزيائية]] و[[هندسة كيميائية|الكيميائية]].
* معالجة الصور وعملية إعادة تكوين [[بكسل|الكثافة النقطية للصورة]] ([https://web.archive.org/web/20150503180109/http://www.aols.com/colorite/interpolation.html decimation])
* مسائل [[ديناميكا حرارية|الديناميكا الحرارية]].
* قطاعات [[تصميم بمساعدة حاسوب|التصميم]] و[[رسوميات|الجرافيك]] والتي تستخدم بكثرة طرق '''الاستكمال''' الخطي([[استيفاء خطي]]) أو ما بعرف بـ"[http://graphics.wikia.com/wiki/Linear_interpolation LERP]".

=== أمثلة محلولة ===
'''''مثال(1):'''''                                                                                                                                                       فيما يلي جدول يحوي بيانات [[درجة حرارة|لدرجة الحرارة]] و [[ضغط|الضغط]] [[غاز|لغاز]]  . '''باستخدام صيغة نيوتن أوجد ضغط الغاز عند درجة حرارة 142°م.'''
{| class="wikitable sortable"
{| class="wikitable"
!درجة الحرارة ( [[درجة حرارة مئوية|م°]])  !! 140 !! 150 !! 160 !! 170 !! 180
|-
| الضغط ([[باسكال (وحدة)|باسكال]]) || 3.685 || 4.854 ||6.302 || 8.076 || 10.225
|}				

'''''الحل:'''''

نلاحظ أن 142°م تقترب من بداية الجدول لذلك نستخدم '''صيغة نيوتن الأمامية'''.
{| class="wikitable"
|-
! x (درجة الحرارة) !! y (الضغط) !!  <math>\triangle </math>  !!  <math>\triangle^2 </math>  !!  <math>\triangle^3 </math>!!  <math>\triangle^4 </math> !!
|-
| 140 || 3.685 || 1.169  ||0.279 || 0.047 || 0.002
|-
| 150 || 4.854 || 1.448 || 0.326 || 0.049 ||                                          
|-
|160 || 6.302 || 1.774 || 0.375 ||   || 
|-
| 170 || 8.076 || 2.149 ||   ||   || 
|-
| 180 || 10.225 ||   ||   ||   || 
|-
|}
لدينا :            <math>x_o=140, y_o=3.685, h=10 </math>    بالتالي :   <math>t=2\frac{x-x_o}{h}=\frac{142-140}{10}=0.2 </math>

لدينا '''صيغة نيوتن''': 
 <math>y_t=y_o+ t\triangle y_o+{t(t-1)\over2!}\triangle^2y_o+{t(t-1)(t-2)\over3!}\triangle^3y_o+ ...</math>
<math>=3.685+0.2(1.199) +{(0.2)(0.2-1)\over2!}0.279+{(0.2)(0.2-1)(0.2-2)\over3!}0.047+{0.2(0.2-1)(0.2-2)( 0.2-3)\over4!}0.002 </math>

<math>\simeq3.899 </math>

'''إذن ضغط الغاز عندما تبلغ درجة حرارته 142 م° هو 3.899 باسكال.'''

__________________________________________________

'''''مثال(2):'''''
'''استخدم صيغة نيوتن للاستكمال في إيجاد عدد الطلاب الذين حصلوا على علامات ليست أكثر من 45 من خلال البيانات التالية :'''
{| class="wikitable"
|-
! العلامات !! 30-40 !! 40-50 !! 50-60!! 60-70 !! 70-80
|-
| عدد الطلاب|| 31 || 42 || 51 || 35 || 31
|-
|}

'''''الحل:'''''
{| class="sortable"
|-
! العلامات !! عدد الطلاب !!  <math>\triangle </math>  !!  <math>\triangle^2 </math>  !!  <math>\triangle^3 </math>!!  <math>\triangle^4 </math> !!
|-
| أقل من40 || 31 ||  ||  ||  || 
|-
| أقل من50 || 42 || 42 || 9 ||  || 
|-
| أقل من 60|| 51|| 51 || 16-|| 25-||
|-
| أقل من 70|| 35|| 35|| 4-|| 17|| 37 ||
|-
| أقل من 80|| 31|| 31|| ||  ||
|}

بالتالي فإن:  <math>t={45-40\over10}={5\over10} </math>

وباستخدام '''صيغة نيوتن الأمامية''' نجد أن:

<math>y= 31+0.5(42)+\frac{0.5(0.5-1)}{2!}(9)+\frac{0.5(0.5-1)(0.5-2)}{3!}(-25)+\frac{0.5(0.5-1)(0.5-2)(0.5-3)}{4!}(37)</math>

<math>= 47.868\simeq 48 </math>

'''إذن عدد الطلاب الذين حصلوا على علامات ليست أكثر من 45 هو 48 طالب.'''

== مراجع ==
{{مراجع}}
# A First Course In Numerical Analysis, second edition, by: Anthony Ralston and Philip Rabinowitz 	
# Numerical Methods, by: S Kalavathy
{{شريط بوابات|رياضيات}}

[[تصنيف:رياضيات]]