الفرق المحدود

الفرق المحدود هو تعبير رياضي على الصورة f ( x + b ) − f ( x + a ) . تُستخدم الفروق المحدودة (أو نسب الفروق المرتبطة بها) غالبًا كتقريبات للمشتقات، كما هو الحال في التفاضل العددي .

عامل الفرق ، والذي يُشار إليه عادةً بـΔ{\displaystyle \Delta }( دلتا الكبيرة )، هو العامل الذي يحول الدالة f إلى الدالةΔ[و]{\displaystyle \Delta [f]}محدد بواسطة Δ[و](x)=و(x+1)-و(x).{\displaystyle \Delta [f](x)=f(x+1)-f(x).} معادلة الفرق هي معادلة دالية تتضمن عامل الفرق المحدود، تمامًا كما تتضمن المعادلة التفاضلية المشتقات . ثمة أوجه تشابه عديدة بين معادلات الفرق والمعادلات التفاضلية. ويمكن كتابة بعض العلاقات التكرارية على شكل معادلات فرق باستبدال رمز التكرار بالفروق المحدودة.

في التحليل العددي ، تُستخدم الفروق المحدودة على نطاق واسع لتقريب المشتقات ، ويُستخدم مصطلح "الفروق المحدودة" غالبًا كاختصار لعبارة "تقريب الفروق المحدودة للمشتقات". [ 1 ] [ 2 ] [ 3 ]

طُرحت الفروق المحدودة لأول مرة على يد بروك تايلور عام ١٧١٥، ودُرست أيضًا ككائنات رياضية مجردة قائمة بذاتها في أعمال جورج بول (١٨٦٠)، وإل إم ميلن-تومسون (١٩٣٣)، وكارولي جوردان (١٩٣٩). تعود أصول الفروق المحدودة إلى إحدى خوارزميات جوست بورغي ( حوالي ١٥٩٢ ) وأعمال أخرى، من بينها أعمال إسحاق نيوتن . ويمكن اعتبار حساب الفروق المحدودة الرسمي بديلًا لحساب المتناهيات في الصغر . [ ٤ ]

الأنواع الأساسية

الأنواع الثلاثة للفروق المحدودة. يعطي الفرق المركزي حول x أفضل تقريب لمشتقة الدالة عند x .

تُعتبر ثلاثة أنواع أساسية شائعة الاستخدام: الفروق المحدودة الأمامية ، والفروق المحدودة الخلفية ، والفروق المحدودة المركزية . [ 1 ] [ 2 ] [ 3 ]

أالفرق الأمامي ، المشار إليه بـΔح[و]{\displaystyle \Delta _{h}[f]}الدالة f هي دالة معرفة على النحو التاليΔح[و](x)=و(x+ح)-و(x).{\displaystyle \Delta _{h}[f](x)=f(x+h)-f(x).}

بحسب التطبيق، قد تكون المسافة h متغيرة أو ثابتة. وعندما لا يتم تحديدها، تكون القيمة الافتراضية لـ h هي 1؛ أي Δ[و](x)=Δ1[و](x)=و(x+1)-و(x).{\displaystyle \Delta [f](x)=\Delta _{1}[f](x)=f(x+1)-f(x).}

أتستخدم طريقة الفرق العكسي قيم الدالة عند x و x h ، بدلاً من القيم عند x + h و x :ح[و](x)=و(x)-و(x-ح)=Δح[و](x-ح).{\displaystyle \nabla _{h}[f](x)=f(x)-f(xh)=\Delta _{h}[f](xh).}

وأخيراً،يُعطى الفرق المركزي بواسطةدلتاح[و](x)=و(x+ح2)-و(x-ح2)=Δح2[و](x)+ح2[و](x).{\displaystyle \delta _{h}[f](x)=f(x+{\tfrac {h}{2}})-f(x-{\tfrac {h}{2}})=\Delta _{\tfrac {h}{2}}[f](x)+\nabla _{\tfrac {h}{2}}[f](x).}

العلاقة بالمشتقات

يلعب تقريب المشتقات بواسطة الفروق المحدودة دورًا محوريًا في طرق الفروق المحدودة للحل العددي للمعادلات التفاضلية ، وخاصة مسائل القيم الحدية .

مشتقة الدالة f عند النقطة x تُعرَّف بالنهايةو(x)=ليمح0و(x+ح)-و(x)ح.{\displaystyle f'(x)=\lim _{h\to 0}{\frac {f(x+h)-f(x)}{h}}.}

إذا كانت قيمة h ثابتة (غير صفرية) بدلاً من أن تقترب من الصفر، فسيتم كتابة الطرف الأيمن من المعادلة أعلاه على النحو التالي و(x+ح)-و(x)ح=Δح[و](x)ح.{\displaystyle {\frac {f(x+h)-f(x)}{h}}={\frac {\Delta _{h}[f](x)}{h}}.}

لذا، فإن الفرق الأمامي مقسومًا على h يُقارب المشتقة عندما تكون h صغيرة. ويمكن استنتاج الخطأ في هذا التقريب من نظرية تايلور . بافتراض أن f قابلة للتفاضل مرتين، لدينا Δح[و](x)ح-و(x)=o(ح)0مثل ح0.{\displaystyle {\frac {\Delta _{h}[f](x)}{h}}-f'(x)=o(h)\to 0\quad {\text{عندما }}h\to 0.}

ينطبق نفس القانون على الفرق العكسي: ح[و](x)ح-و(x)=o(ح)0مثل ح0.{\displaystyle {\frac {\nabla _{h}[f](x)}{h}}-f'(x)=o(h)\to 0\quad {\text{عندما }}h\to 0.}

مع ذلك، فإن الفرق المركزي (أو ما يُسمى بالفرق المُمركز) يُعطي تقريبًا أكثر دقة. إذا كانت الدالة f قابلة للتفاضل ثلاث مرات، دلتاح[و](x)ح-و(x)=o(ح2).{\displaystyle {\frac {\delta _{h}[f](x)}{h}}-f'(x)=o\left(h^{2}\right).}

تكمن المشكلة الرئيسية في طريقة الفروق المركزية في أن الدوال المتذبذبة قد تُعطي مشتقة صفرية. فإذا كانت f ( nh ) = 1 عندما يكون n فرديًا، و f ( nh ) = 2 عندما يكون n زوجيًا، فإن f ′( nh ) = 0 إذا حُسبت باستخدام طريقة الفروق المركزية . ويُصبح هذا الأمر إشكاليًا بشكل خاص إذا كان مجال f منفصلًا. انظر أيضًا: المشتقة المتناظرة .

يُعرّف المؤلفون الذين يعتبرون الفروق المحدودة بمثابة تقريبات للفروق المحدودة، الفروق الأمامية/الخلفية/المركزية بأنها نواتج القسمة الواردة في هذا القسم (بدلاً من استخدام التعريفات الواردة في القسم السابق). [ 1 ] [ 2 ] [ 3 ]

الفروق من الرتبة العليا

وبطريقة مماثلة، يمكن الحصول على تقريبات الفروق المحدودة للمشتقات والمؤثرات التفاضلية من الرتب العليا. على سبيل المثال، باستخدام صيغة الفروق المركزية المذكورة أعلاه لـ f ′( x + h / 2 ) و f ′( xh / 2 ) وتطبيق صيغة الفروق المركزية لمشتقة f عند x ، نحصل على تقريب الفروق المركزية للمشتقة الثانية لـ f : 

مركزي من الدرجة الثانية
و"(x)دلتاح2[و](x)ح2=و(x+ح)-و(x)ح-و(x)-و(x-ح)حح=و(x+ح)-2و(x)+و(x-ح)ح2 .{\displaystyle f''(x)\approx {\frac {\delta _{h}^{2}[f](x)}{h^{2}}}={\frac {{\frac {f(x+h)-f(x)}{h}}-{\frac {f(x)-f(xh)}{h}}}{h}}={\frac {f(x+h)-2f(x)+f(xh)}{h^{2}}}~.}

وبالمثل، يمكننا تطبيق صيغ التفاضل الأخرى بطريقة تكرارية.

مهاجم من الدرجة الثانية
و"(x)Δح2[و](x)ح2=و(x+2ح)-و(x+ح)ح-و(x+ح)-و(x)حح=و(x+2ح)-2و(x+ح)+و(x)ح2 .{\displaystyle f''(x)\approx {\frac {\Delta _{h}^{2}[f](x)}{h^{2}}}={\frac {{\frac {f(x+2h)-f(x+h)}{h}}-{\frac {f(x+h)-f(x)}{h}}}{h}}={\frac {f(x+2h)-2f(x+h)+f(x)}{h^{2}}}~.}
المعادلة العكسية من الدرجة الثانية
و"(x)ح2[و](x)ح2=و(x)-و(x-ح)ح-و(x-ح)-و(x-2ح)حح=و(x)-2و(x-ح)+و(x-2ح)ح2 .{\displaystyle f''(x)\approx {\frac {\nabla _{h}^{2}[f](x)}{h^{2}}}={\frac {{\frac {f(x)-f(xh)}{h}}-{\frac {f(xh)-f(x-2h)}{h}}}{h}}={\frac {f(x)-2f(xh)+f(x-2h)}{h^{2}}}~.}

وبشكل أعم، تُعطى الفروق الأمامية والخلفية والمركزية من الرتبة n على التوالي،

إلى الأمام
Δحن[و](x)=ج=0ن(-1)ن-ج(نج)و(x+جح)،{\displaystyle \Delta _{h}^{n}[f](x)=\sum _{j=0}^{n}(-1)^{n-j}{\binom {n}{j}}f{\bigl (}x+jh{\bigr )},}
إلى الخلف
حن[و](x)=ج=0ن(-1)ج(نج)و(x-جح)،{\displaystyle \nabla _{h}^{n}[f](x)=\sum _{j=0}^{n}(-1)^{j}{\binom {n}{j}}f(x-jh),}
وسط المدينة
دلتاحن[و](x)=ج=0ن(-1)ج(نج)و(x+(ن2-ج)ح) .{\displaystyle \delta _{h}^{n}[f](x)=\sum _{j=0}^{n}(-1)^{j}{\binom {n}{j}}f\left(x+\left({\frac {n}{2}}-j\right)h\right)~.}

تستخدم هذه المعادلات معاملات ذات الحدين بعد علامة الجمع الموضحة كما يلي: (نج) .{\textstyle \ {\binom {n}{j}}~.}يوفر كل صف من صفوف مثلث باسكال معاملًا لكل قيمة من قيم j .

لاحظ أن الفرق المركزي، بالنسبة لقيم n الفردية ، سيحتوي على h مضروبة في أعداد غير صحيحة. غالبًا ما يمثل هذا مشكلة لأنه يؤدي إلى تغيير فترة التقطيع. يمكن معالجة هذه المشكلة باستبدال متوسطدلتان[و](x-ح2){\textstyle \delta ^{n}[f](x-{\tfrac {h}{2}})} و دلتان[و](x+ح2).{\textstyle \delta ^{n}[f](x+{\tfrac {h}{2}}).}

تُسمى الفروق الأمامية المطبقة على متتالية ما أحيانًا بالتحويل الثنائي للمتتالية، ولها عدد من الخصائص التوافقية المهمة. يمكن حساب الفروق الأمامية باستخدام تكامل نورلوند-رايس . يُعد التمثيل التكاملي لهذه الأنواع من المتسلسلات مثيرًا للاهتمام، لأنه غالبًا ما يمكن حساب التكامل باستخدام التوسع التقاربي أو تقنيات نقطة السرج ؛ في المقابل، قد يكون حساب متسلسلة الفروق الأمامية عدديًا صعبًا للغاية، لأن معاملات ذات الحدين تنمو بسرعة مع ازدياد قيمة n .

إن العلاقة بين هذه الفروق ذات الرتبة الأعلى والمشتقات المقابلة لها علاقة مباشرة، دنو(دx)ن(x)=Δحن[و](x)حن+o(ح)=حن[و](x)حن+o(ح)=دلتاحن[و](x)حن+o(ح2).{\displaystyle {\frac {\mathrm {d} ^{n}f}{(\mathrm {d} x)^{n}}}(x)={\frac {\Delta _{h}^{n}[f](x)}{h^{n}}}+{\mathcal {o}}(h)={\frac {\nabla _{h}^{n}[f](x)}{h^{n}}}+{\mathcal {o}}(h)={\frac {\delta _{h}^{n}[f](x)}{h^{n}}}+{\mathcal {o}}\!\left(h^{2}\right).}

يمكن أيضًا استخدام الفروق ذات الرتب الأعلى لإنشاء تقريبات أفضل. كما ذُكر أعلاه، فإن الفرق من الرتبة الأولى يُقارب المشتقة من الرتبة الأولى حتى حد من الرتبة h . ومع ذلك، فإن الجمع Δح[و](x)-12Δح2[و](x)ح=-و(x+2ح)-4و(x+ح)+3و(x)2ح{\displaystyle {\frac {\Delta _{h}[f](x)-{\frac {1}{2}}\Delta _{h}^{2}[f](x)}{h}}=-{\frac {f(x+2h)-4f(x+h)+3f(x)}{2h}}} تقارب f ′( x ) حتى حد من الرتبة h 2. ويمكن إثبات ذلك عن طريق توسيع التعبير أعلاه في متسلسلة تايلور ، أو باستخدام حساب الفروق المحدودة، الموضح أدناه.

إذا لزم الأمر، يمكن توسيط الفرق المحدود حول أي نقطة عن طريق مزج الفروق الأمامية والخلفية والمركزية.

أحيانًا، قد تكون المشتقات ذات الرتب المنخفضة لدالة ما معروفة تحليليًا، بينما لا تكون المشتقات ذات الرتب العالية كذلك. في هذه الحالات، يمكن تقريب المشتقات ذات الرتب العالية باستخدام الفروق المحدودة للمشتقات ذات الرتب المنخفضة، وهو ما يكون غالبًا أكثر دقة واستقرارًا عدديًا من الفروق المحدودة للدالة f ( x ) نفسها. يُطلق على هذا أحيانًا اسم التفاضل شبه العددي. [ 5 ] على سبيل المثال، عندما تكون المشتقة من الرتبة الأولى f ′( x ) متاحة، ولكن المشتقة من الرتبة الثانية f ″( x ) غير متاحة، يمكن تقريب الأخيرة باستخدام الفروق المركزية من الرتبة الثانية لـ f ′( x ) . و"(x)و(x+ح)-و(x-ح)2ح .{\displaystyle f''(x)\approx {\frac {f'(x+h)-f'(x-h)}{2h}}~.}

كثيرات الحدود

بالنسبة لكثير حدود معين من الدرجة n 1 ، معبر عنه بالدالة P ( x ) ، مع الأعداد الحقيقية a 0 و b والحدود ذات الرتبة الأدنى (إن وجدت) المميزة بـ lot : P(x)=أxن+بxن-1+ ل.o.ت.{\displaystyle P(x)=ax^{n}+\;bx^{n-1}+~l.o.t.}

بعد إجراء n من الفروق الزوجية، يمكن تحقيق النتيجة التالية، حيث h 0 هو عدد حقيقي يشير إلى الفرق الحسابي: [ 6 ]Δحن[P](x)=أحنن!{\displaystyle \Delta _{h}^{n}[P](x)=ah^{n}n!}

يبقى فقط معامل الحد ذي الرتبة الأعلى. وبما أن هذه النتيجة ثابتة بالنسبة لـ x ، فإن أي فروق ثنائية أخرى ستكون قيمتها صفرًا .

البرهان الاستقرائي

الحالة الأساسية

ليكن Q ( x ) متعدد الحدود من الدرجة 1 : Δح[سؤال](x)=سؤال(x+ح)-سؤال(x)=[أ(x+ح)+ب]-[أx+ب]=أح=أح11!{\displaystyle \Delta _{h}[Q](x)=Q(x+h)-Q(x)=[a(x+h)+b]-[ax+b]=ah=ah^{1}1!}

وهذا يثبت ذلك بالنسبة للحالة الأساسية.

خطوات استقرائية

ليكن R ( x ) متعدد حدود من الدرجة m  1 حيث m 2 ومعامل الحد ذي الرتبة الأعلى a 0. بافتراض صحة ما يلي لجميع متعددات الحدود من الدرجة m  1 : Δحم-1[R](x)=أحم-1(م-1)!{\displaystyle \Delta _{h}^{m-1}[R](x)=ah^{m-1}(m-1)!}

ليكن S ( x ) متعدد حدود من الدرجة m . مع فرق واحد بين كل زوج من الحدود: Δح[S](x)=[أ(x+ح)م+ب(x+ح)م-1+كثير]-[أxم+بxم-1+كثير]=أحمxم-1+كثير=تي(x){\displaystyle \Delta _{h}[S](x)=[a(x+h)^{m}+b(x+h)^{m-1}+{\text{l.o.t.}}]-[ax^{m}+bx^{m-1}+{\text{l.o.t.}}]=ahmx^{m-1}+{\text{l.o.t.}}=T(x)}

بما أن ahm 0 ، ينتج عن ذلك متعددة حدود T ( x ) من الدرجة m − 1 ، حيث ahm هو معامل الحد ذي الرتبة الأعلى. وبالنظر إلى الافتراض أعلاه و m − 1 فرقًا ثنائيًا (مما ينتج عنه إجمالي m فرقًا ثنائيًا لـ S ( x ) )، يمكن إيجاد ما يلي: Δحم-1[تي](x)=أحمحم-1(م-1)!=أحمم!{\displaystyle \Delta _{h}^{m-1}[T](x)=ahm\cdot h^{m-1}(m-1)!=ah^{m}m!}

وبهذا يكتمل البرهان.

طلب

يمكن استخدام هذه المتطابقة لإيجاد متعددة الحدود ذات أدنى درجة والتي تقطع عددًا من النقاط ( x , y ) حيث يكون الفرق على المحور السيني من نقطة إلى أخرى ثابتًا h 0. على سبيل المثال، بالنظر إلى النقاط التالية:

xy
14
4109
7772
102641
136364

يمكننا استخدام جدول الفروقات، حيث توجد العلاقة التالية بين جميع الخلايا الموجودة على يمين أول قيمة لـ y والخلايا الموجودة في العمود المجاور لها مباشرةً على اليسار، وذلك بالنسبة للخلية ( a + 1, b + 1) ، مع كون الخلية الموجودة في أعلى اليسار عند الإحداثيات (0, 0) : (أ+1،ب+1)=(أ،ب+1)-(أ،ب){\displaystyle (a+1,b+1)=(a,b+1)-(a,b)}

لإيجاد الحد الأول، يمكن استخدام الجدول التالي:

xyΔ yΔ 2 yΔ 3 y
14
4109105
7772663558
10264118691206648
13636437231854648

وهذا يُؤدي إلى قيمة ثابتة 648. والفرق الحسابي هو h = 3 ، كما هو مُبين أعلاه. ونظرًا لعدد الفروق الزوجية اللازمة للوصول إلى هذه القيمة الثابتة، يُمكن استنتاج أن هذه كثيرة حدود من الدرجة الثالثة . وبالتالي، باستخدام المتطابقة أعلاه: 648=أ333!=أ276=أ162{\displaystyle 648=a\cdot 3^{3}\cdot 3!=a\cdot 27\cdot 6=a\cdot 162}

بحل المعادلة لإيجاد قيمة a ، نجد أنها تساوي 4. وبالتالي، فإن الحد الأول من كثيرة الحدود هو 4 × 3 .

ثم، بطرح الحد الأول، مما يقلل من درجة متعددة الحدود، وإيجاد الفرق المحدود مرة أخرى:

xyΔ yΔ 2 y
14 − 4(1) 3 = 4 − 4 = 0
4109 − 4(4) 3 = 109 − 256 = −147-147
7772 − 4(7) 3 = 772 − 1372 = −600-453-306
102641 − 4(10) 3 = 2641 − 4000 = −1359-759-306
136364 − 4(13) 3 = 6364 − 8788 = −2424-1065-306

هنا، يتم الوصول إلى الثابت بعد إجراء فرقين زوجيين فقط، وبالتالي تكون النتيجة كالتالي: -306=أ322!=أ18{\displaystyle -306=a\cdot 3^{2}\cdot 2!=a\cdot 18}

وبحل المعادلة لإيجاد قيمة a ، وهي −17 ، فإن الحد الثاني لكثير الحدود هو −17 × 2 .

بالانتقال إلى الحد التالي، عن طريق طرح الحد الثاني:

xyΔ y
10 − (−17(1) 2 ) = 0 + 17 = 17
4-147 - (-17(4) 2 ) = -147 + 272 = 125108
7-600 - (-17(7) 2 ) = -600 + 833 = 233108
10-1359 - (-17(10) 2 ) = -1359 + 1700 = 341108
13-2424 - (-17(13) 2 ) = -2424 + 2873 = 449108

وبالتالي، يتم تحقيق الثابت بعد فرق ثنائي واحد فقط: 108=أ311!=أ3{\displaystyle 108=a\cdot 3^{1}\cdot 1!=a\cdot 3}

يمكن إيجاد أن a = 36، وبالتالي فإن الحد الثالث من كثيرة الحدود هو 36x . بطرح الحد الثالث:

xy
117 − 36(1) = 17 − 36 = −19
4125 − 36(4) = 125 − 144 = −19
7233 − 36(7) = 233 − 252 = −19
10341 − 36(10) = 341 − 360 = −19
13449 - 36(13) = 449 - 468 = -19

بدون أي فروق ثنائية، وُجد أن الحد الرابع والأخير من كثيرة الحدود هو الثابت -19 . وبالتالي، تم إيجاد كثيرة الحدود ذات أدنى درجة والتي تقطع جميع النقاط في الجدول الأول: 4x3-17x2+36x-19{\displaystyle 4x^{3}-17x^{2}+36x-19}

حبات ذات أحجام عشوائية

باستخدام الجبر الخطي، يمكن بناء تقريبات الفروق المحدودة التي تستخدم عددًا عشوائيًا من النقاط على يسار نقطة التقييم وعددًا (قد يكون مختلفًا) من النقاط على يمينها، لأي رتبة من المشتقات. يتضمن ذلك حل نظام خطي بحيث يكون توسيع تايلور لمجموع تلك النقاط حول نقطة التقييم هو الأقرب إلى توسيع تايلور للمشتقة المطلوبة. يمكن تمثيل هذه الصيغ بيانيًا على شبكة سداسية أو معينية الشكل. [ 7 ] يُفيد هذا في اشتقاق دالة على شبكة، حيث كلما اقتربنا من حافة الشبكة، يجب أخذ عينات من عدد أقل فأقل من النقاط على أحد الجانبين. [ 8 ] يمكن بناء تقريبات الفروق المحدودة لقوالب غير قياسية (وحتى غير صحيحة) بالنظر إلى قالب عشوائي ورتبة مشتقة مطلوبة. [ 9 ]

ملكيات

  • لجميع قيم k و n الموجبةΔكحن(و،x)=ج1=0ك-1ج2=0ك-1جن=0ك-1Δحن(و،x+ج1ح+ج2ح++جنح).{\displaystyle \Delta _{kh}^{n}\left(f,x\right)=\sum \limits _{j_{1}=0}^{k-1}\sum \limits _{j_{2}=0}^{k-1}\cdots \sum \limits _{j_{n}=0}^{k-1}\Delta _{h}^{n}\left(f,x+j_{1}h+j_{2}h+\cdots +j_{n}h\right).}
  • قاعدة لايبنتز :Δحن(وز،x)=ك=0ن(نك)Δحك(و،x)Δحن-ك(ز،x+كح).{\displaystyle \Delta _{h}^{n}\left(fg,x\right)=\sum \limits _{k=0}^{n}{\binom {n}{k}}\Delta _{h}^{k}\left(f,x\right)\cdot \Delta _{h}^{n-k}\left(g,x+kh\right).}

في المعادلات التفاضلية

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

تتواجد التطبيقات الشائعة لطريقة الفروق المحدودة في مجالات العلوم والهندسة الحاسوبية، مثل الهندسة الحرارية وميكانيكا الموائع وما إلى ذلك.

متسلسلة نيوتن

تتكون متسلسلة نيوتن من حدود معادلة نيوتن التفاضلية الأمامية ، التي سميت على اسم إسحاق نيوتن ؛ في جوهرها، هي صيغة استيفاء غريغوري-نيوتن [ 10 ] (التي سميت على اسم إسحاق نيوتن وجيمس غريغوري )، والتي نُشرت لأول مرة في كتابه Principia Mathematica عام 1687، [ 11 ] [ 12 ] أي النظير المنفصل لتوسيع تايلور المستمر،

و(x)=ك=0Δك[و](أ)ك!(x-أ)ك=ك=0(x-أك)Δك[و](أ)،{\displaystyle f(x)=\sum _{k=0}^{\infty }{\frac {\Delta ^{k}[f](a)}{k!}}\,(x-a)_{k}=\sum _{k=0}^{\infty }{\binom {x-a}{k}}\,\Delta ^{k}[f](a),}

وهذا ينطبق على أي دالة متعددة الحدود f وعلى العديد من الدوال التحليلية (وليس جميعها) . (لا ينطبق هذا عندما تكون f من النوع الأسي).π{\displaystyle \pi }وهذا واضحٌ بسهولة، حيث أن دالة الجيب تتلاشى عند المضاعفات الصحيحة لـπ{\displaystyle \pi }(إن متسلسلة نيوتن المقابلة تساوي صفرًا تمامًا، لأن جميع الفروق المحدودة تساوي صفرًا في هذه الحالة. ومع ذلك، من الواضح أن دالة الجيب لا تساوي صفرًا.) هنا، التعبير (xك)=(x)كك!{\displaystyle {\binom {x}{k}}={\frac {(x)_{k}}{k!}}} هو معامل ذات الحدين ، و (x)ك=x(x-1)(x-2)(x-ك+1){\displaystyle (x)_{k}=x(x-1)(x-2)\cdots (x-k+1)} هو " المضروب الهابط " أو "المضروب السفلي"، بينما يُعرَّف الناتج الفارغ ( x ) 0  بأنه 1. في هذه الحالة بالذات، هناك افتراض لخطوات الوحدة للتغيرات في قيم x ، h = 1 للتعميم أدناه.

لاحظ التطابق الرسمي لهذه النتيجة مع نظرية تايلور . تاريخيًا، هذه النتيجة، بالإضافة إلى متطابقة تشو-فاندرموند ، (x+y)ن=ك=0ن(نك)(x)ن-ك(y)ك،{\displaystyle (x+y)_{n}=\sum _{k=0}^{n}{\binom {n}{k}}(x)_{n-k}\,(y)_{k},} (النتائج المترتبة على ذلك، والمتوافقة مع نظرية ذات الحدين )، يتم تضمينها في الملاحظات التي تطورت إلى نظام حساب الظلال .

قد تكون متسلسلات نيوتن أفضل من متسلسلات تايلور عند تطبيقها على كميات منفصلة مثل اللف المغزلي الكمومي (انظر تحويل هولشتاين-بريماكوف )، أو دوال المؤثرات البوزونية ، أو إحصاءات العد المنفصلة. [ 13 ]

لتوضيح كيفية استخدام صيغة نيوتن في الممارسة العملية، لننظر إلى الحدود القليلة الأولى لمضاعفة متتالية فيبوناتشي f = 2، 2، 4، ... يمكن إيجاد متعددة حدود تعيد إنتاج هذه القيم، عن طريق حساب جدول الفروق أولاً، ثم استبدال الفروق التي تتوافق مع x 0 (المسطرة) في الصيغة كما يلي، xو=Δ0Δ1Δ212_0_222_234و(x)=Δ01+Δ1(x-x0)11!+Δ2(x-x0)22!(x0=1)=21+0x-11+2(x-1)(x-2)2=2+(x-1)(x-2){\displaystyle {\begin{matrix}{\begin{array}{|c||c|c|c|}\hline x&f=\Delta ^{0}&\Delta ^{1}&\Delta ^{2}\\\hline 1&{\underline {2}}&&\\&&{\underline {0}}&\\2&2&&{\underline {2}}\\&&2&\\3&4&&\\\hline \end{array}}&\quad {\begin{aligned}f(x)&=\Delta ^{0}\cdot 1+\Delta ^{1}\cdot {\dfrac {(x-x_{0})_{1}}{1!}}+\Delta ^{2}\cdot {\dfrac {(x-x_{0})_{2}}{2!}}\quad (x_{0}=1)\\\\&=2\cdot 1+0\cdot {\dfrac {x-1}{1}}+2\cdot {\dfrac {(x-1)(x-2)}{2}}\\\\&=2+(x-1)(x-2)\\\end{aligned}}\end{matrix}}}

في حالة وجود خطوات غير منتظمة في قيم x ، يحسب نيوتن الفروق المقسمة . Δج،0=yج،Δج،ك=Δج+1،ك-1-Δج،ك-1xج+ك-xج{ك>0،جالأعلى(ج)-ك}،Δ0ك=Δ0،ك{\displaystyle \Delta _{j,0}=y_{j},\qquad \Delta _{j,k}={\frac {\Delta _{j+1,k-1}-\Delta _{j,k-1}}{x_{j+k}-x_{j}}}\quad \ni \quad \left\{k>0,\;j\leq \max \left(j\right)-k\right\},\qquad \Delta 0_{k}=\Delta _{0,k}} سلسلة المنتجات، P0=1،Pك+1=Pك(ξ-xك)،{\displaystyle {P_{0}}=1,\quad \quad P_{k+1}=P_{k}\cdot \left(\xi -x_{k}\right),} والمتعددة الحدود الناتجة هي الضرب القياسي ، [ 14 ]و(ξ)=Δ0P(ξ).{\displaystyle f(\xi )=\Delta 0\cdot P\left(\xi \right).}

في التحليل باستخدام الأعداد p -adic ، تنص نظرية ماهلر على أنه يمكن إضعاف الافتراض بأن f دالة متعددة الحدود وصولاً إلى الافتراض بأن f متصلة فقط.

تُقدّم نظرية كارلسون شروطًا ضرورية وكافية لكي تكون متسلسلة نيوتن فريدة، إن وُجدت. مع ذلك، فإن متسلسلة نيوتن غير موجودة، بشكل عام.

تُعد سلسلة نيوتن، إلى جانب سلسلة ستيرلنغ وسلسلة سيلبرغ ، حالة خاصة من سلسلة الفروق العامة ، والتي يتم تعريفها جميعًا من حيث الفروق الأمامية ذات المقياس المناسب.

في صيغة مختصرة وأكثر عمومية قليلاً، ومع وجود عقد متساوية البعد، تصبح الصيغة كالتالي: و(x)=ك=0(x-أحك)ج=0ك(-1)ك-ج(كج)و(أ+جح).{\displaystyle f(x)=\sum _{k=0}{\binom {\frac {x-a}{h}}{k}}\sum _{j=0}^{k}(-1)^{k-j}{\binom {k}{j}}f(a+jh).}

حساب التفاضل والتكامل للفروق المحدودة

يمكن اعتبار الفرق الأمامي بمثابة مؤثر ، يُسمى مؤثر الفرق ، والذي يحول الدالة f إلى Δh [ f ] . [ 15 ] [ 16 ] هذا المؤثر يُعادلΔح=تيح-أنا،{\displaystyle \Delta _{h}=\operatorname {T} _{h}-\operatorname {I} ,}حيث T h هو عامل الإزاحة بخطوة h ، المعرف بواسطة T h [ f ]( x ) = f ( x + h ) ، و I هو عامل الهوية .

يمكن تعريف الفرق المحدود من الرتب العليا بطريقة تكرارية على النحو التالي: Δ n h ≡ Δ hn − 1 h ) . وهناك تعريف مكافئ آخر هو: Δ n h ≡ [T h − I] n .

إن عامل الفرق Δ h هو عامل خطي ، وعلى هذا النحو فإنه يحقق Δ h [ α f + β g ]( x ) = α Δ h [ f ]( x ) + β Δ h [ g ]( x ) .

كما أنه يفي بقاعدة خاصة من قواعد لايبنتز :

Δح(وز)(x)=(Δحو(x))ز(x+ح)+و(x)(Δحز(x)) .{\displaystyle \operatorname {\Delta } _{h}{\bigl (}fg{\bigr )}(x)={\bigl (}\operatorname {\Delta } _{h}f(x){\bigr )}g(x+h)+f(x){\bigl (}\operatorname {\Delta } _{h}g(x){\bigr )}~.}

تنطبق قواعد لايبنتز المماثلة على الفروق الخلفية والمركزية.

بتطبيق متسلسلة تايلور رسميًا بالنسبة إلى h ، نحصل على معادلة المؤثر Δح=حد+12!ح2د2+13!ح3د3+=هـحد-أنا،{\displaystyle \operatorname {\Delta } _{h}=h\operatorname {D} +{\frac {1}{2!}}h^{2}\operatorname {D} ^{2}+{\frac {1}{3!}}h^{3}\operatorname {D} ^{3}+\cdots =e^{h\operatorname {D} }-\operatorname {I} ,} حيث يرمز D إلى مُعامل الاشتقاق التقليدي المتصل، الذي يُحول f إلى مشتقتها f . يكون هذا التوسع صحيحًا عندما يؤثر كلا الطرفين على دوال تحليلية ، وذلك لقيم h صغيرة بما فيه الكفاية ؛ في الحالة الخاصة التي تنتهي فيها سلسلة المشتقات (عندما تكون الدالة التي يُعمل عليها متعددة حدود منتهية )، يكون التعبير دقيقًا، لجميع أحجام الخطوات المنتهية h . وبالتالي، T h = e h D ، وبقلب الدالة الأسية رسميًا نحصل على حد=ln(1+Δح)=Δح-12Δح2+13Δح3- .{\displaystyle h\operatorname {D} =\ln(1+\Delta _{h})=\Delta _{h}-{\tfrac {1}{2}}\,\Delta _{h}^{2}+{\tfrac {1}{3}}\,\Delta _{h}^{3}-\cdots ~.} هذه الصيغة صحيحة بمعنى أن كلا العاملين يعطيان نفس النتيجة عند تطبيقهما على متعدد الحدود.

حتى بالنسبة للدوال التحليلية، لا يُضمن تقارب المتسلسلة على اليمين؛ فقد تكون متسلسلة تقاربية . مع ذلك، يمكن استخدامها للحصول على تقريبات أكثر دقة للمشتقة. على سبيل المثال، يؤدي الاحتفاظ بالحدين الأولين من المتسلسلة إلى الحصول على التقريب من الرتبة الثانية لـ f  ′( x ) المذكور في نهاية القسم §  الفروق من الرتب العليا .

الصيغ المماثلة لمؤثرات الفروق الخلفية والمركزية هي حد=-ln(1-ح) و حد=2أرسينه(12دلتاح) .{\displaystyle h\operatorname {D} =-\ln(1-\nabla _{h})\quad {\text{ and }}\quad h\operatorname {D} =2\operatorname {arsinh} \left({\tfrac {1}{2}}\,\delta _{h}\right)~.}

يرتبط حساب الفروق المحدودة بحساب الظلال في التوافقية. ويعود هذا التطابق المنهجي الملحوظ إلى تطابق مبدلات كميات الظلال مع نظائرها في الوسط المتصل ( حدود h → 0 ).

[Δحح،xتيح-1]=[د،x]=أنا.{\displaystyle \left[{\frac {\Delta _{h}}{h}},x\,\operatorname {T} _{h}^{-1}\right]=[\operatorname {D} ,x]=I.}

وبالتالي فإن عددًا كبيرًا من العلاقات التفاضلية الرسمية لحساب التفاضل والتكامل القياسي التي تتضمن الدوال f ( x ) يتم ربطها بشكل منهجي بنظائر الفروق المحدودة الظلية التي تتضمن f ( x T −1 h ) .

على سبيل المثال، فإن النظير الظلي للحد الأحادي x n هو تعميم للعامل التنازلي المذكور أعلاه ( رمز بوخامر k ). (x)ن=(xتيح-1)ن=x(x-ح)(x-2ح)(x-(ن-1)ح)،{\displaystyle (x)_{n}=\left(x\operatorname {T} _{h}^{-1}\right)^{n}=x\left(x-h\right)\left(x-2h\right)\cdots {\bigl (}x-\left(n-1\right)h{\bigr )},} لهذا السبب. Δحح(x)ن=ن(x)ن-1،{\displaystyle {\frac {\Delta _{h}}{h}}(x)_{n}=n(x)_{n-1},} ومن ثم صيغة نيوتن للاستيفاء المذكورة أعلاه (عن طريق مطابقة المعاملات في توسيع دالة عشوائية f ( x ) في مثل هذه الرموز)، وهكذا.

على سبيل المثال، الجيب الظلي هو الخطيئة(xتيح-1)=x-(x)33!+(x)55!-(x)77!+{\displaystyle \sin \left(x\operatorname {T} _{h}^{-1}\right)=x-{\frac {(x)_{3}}{3!}}+{\frac {(x)_{5}}{5!}}-{\frac {(x)_{7}}{7!}}+\cdots }

كما هو الحال في حد الاستمرارية ، فإن الدالة الذاتية لـ Δh / h هي أيضًا دالة أسية.

Δحح(1+λح)xح=Δححهـln(1+λح)xح=λهـln(1+λح)xح،{\displaystyle {\frac {\Delta _{h}}{h}}(1+\lambda h)^{\frac {x}{h}}={\frac {\Delta _{h}}{h}}e^{\ln(1+\lambda h){\frac {x}{h}}}=\lambda e^{\ln(1+\lambda h){\frac {x}{h}}},}

وبالتالي ، يمكن بسهولة وبدقة تحويل مجاميع فورييه لدوال الاستمرارية إلى مجاميع فورييه الظلية ، أي التي تتضمن نفس معاملات فورييه مضروبة في هذه الدوال الأسية الأساسية الظلية. [ 17 ] وبالتالي، فإن هذه الدالة الأسية الظلية تُعادل الدالة المولدة الأسية لرموز بوخامر .

وهكذا، على سبيل المثال، تُقابل دالة ديراك دلتا نظيرتها الظلية، وهي دالة الجيب الأساسية.دلتا(x)الخطيئة[π2(1+xح)]π(x+ح)،{\displaystyle \delta (x)\mapsto {\frac {\sin \left[{\frac {\pi }{2}}\left(1+{\frac {x}{h}}\right)\right]}{\pi (x+h)}},} وهكذا دواليك. [ 18 ] يمكن في كثير من الأحيان حل المعادلات الفرقية بتقنيات مشابهة جدًا لتلك المستخدمة في حل المعادلات التفاضلية .

إن المؤثر العكسي لمؤثر الفرق الأمامي، وبالتالي التكامل الظلي، هو المجموع غير المحدد أو مؤثر الفرق المضاد.

قواعد حساب التفاضل والتكامل لمؤثرات الفروق المحدودة

على غرار قواعد إيجاد المشتقة ، لدينا:

  • قاعدة الثوابت : إذا كان c ثابتًا ، فإنΔج=0{\displaystyle \Delta c=0}
  • الخطية : إذاكان a و b ثابتين ،Δ(أو+بز)=أΔو+بΔز{\displaystyle \Delta (af+bg)=a\Delta f+b\Delta g}

تنطبق جميع القواعد المذكورة أعلاه بنفس القدر على أي عامل فرق كما تنطبق على Δ ، بما في ذلك δ و .

  • قاعدة المنتج :Δ(وز)=وΔز+زΔو+ΔوΔز(وز)=وز+زو-وز{\displaystyle {\begin{aligned}\Delta (fg)&=f\,\Delta g+g\,\Delta f+\Delta f\Delta g\\[4pt]\nabla (fg)&=f\,\nabla g+g\,\nabla f-\nabla f\nabla g\end{aligned}}}
  • قاعدة القسمة :(وز)=(المحقق[وزوز])/(زالمحقق[زز11]){\displaystyle \nabla \left({\frac {f}{g}}\right)=\left.\left(\det {\begin{bmatrix}\nabla f&\nabla g\\f&g\end{bmatrix}}\right)\right/\left(g\cdot \det {\begin{bmatrix}g&\nabla g\\1&1\end{bmatrix}}\right)}أو(وز)=زو-وزز(ز-ز){\displaystyle \nabla \left({\frac {f}{g}}\right)={\frac {g\,\nabla f-f\,\nabla g}{g\cdot (g-\nabla g)}}}
  • قواعد الجمع :ن=أبΔو(ن)=و(ب+1)-و(أ)ن=أبو(ن)=و(ب)-و(أ-1){\displaystyle {\begin{aligned}\sum _{n=a}^{b}\Delta f(n)&=f(b+1)-f(a)\\\sum _{n=a}^{b}\nabla f(n)&=f(b)-f(a-1)\end{aligned}}}

انظر المراجع. [ 19 ] [ 20 ] [ 21 ] [ 22 ]

التعميمات

  • يُعرَّف الفرق المحدود المعمم عادةً على النحو التالي :Δحμ[و](x)=ك=0شمالμكو(x+كح)،{\displaystyle \Delta _{h}^{\mu }[f](x)=\sum _{k=0}^{N}\mu _{k}f(x+kh),}حيث μ = ( μ₀, …, μₙ) هو متجه معاملاته. يُعدّ الفرق اللانهائي تعميمًا إضافيًا، حيث يُستبدل المجموع المحدود أعلاه بمتسلسلة لانهائية. وهناك طريقة أخرى للتعميم تتمثل في جعل المعاملات μₖ تعتمد على النقطة x : μₖ = μₖ ( x ) ، وبالتالي اعتبار الفرق المحدود الموزون . كما يمكن جعل الخطوة h تعتمد على النقطة x : h = h ( x ) . تُفيد هذه التعميمات في بناء معاملات استمرارية مختلفة .
  • يمكن اعتبار الفرق المعمم بمثابة حلقات متعددة الحدود R [ T h ] . ويؤدي ذلك إلى جبر الفرق.
  • يتم تعميم عامل الفرق إلى عكس موبيوس على مجموعة مرتبة جزئياً .
  • كعامل التفاف: من خلال شكلية جبر الوقوع ، يمكن تمثيل عوامل الفرق وانعكاس موبيوس الأخرى عن طريق الالتفاف مع دالة على المجموعة المرتبة جزئيًا، تسمى دالة موبيوس μ ؛ بالنسبة لعامل الفرق، فإن μ هي المتتالية (1، -1، 0، 0، 0، ...) .

الفروق المحدودة متعددة المتغيرات

يمكن دراسة الفروق المحدودة في أكثر من متغير واحد. وهي مماثلة للمشتقات الجزئية في عدة متغيرات.

بعض تقريبات المشتقات الجزئية هي: وx(x،y)و(x+ح،y)-و(x-ح،y)2حوy(x،y)و(x،y+ك)-و(x،y-ك)2كوxx(x،y)و(x+ح،y)-2و(x،y)+و(x-ح،y)ح2وyy(x،y)و(x،y+ك)-2و(x،y)+و(x،y-ك)ك2وxy(x،y)و(x+ح،y+ك)-و(x+ح،y-ك)-و(x-ح،y+ك)+و(x-ح،y-ك)4حك.{\displaystyle {\begin{aligned}f_{x}(x,y)&\approx {\frac {f(x+h,y)-f(x-h,y)}{2h}}\\f_{y}(x,y)&\approx {\frac {f(x,y+k)-f(x,y-k)}{2k}}\\f_{xx}(x,y)&\approx {\frac {f(x+h,y)-2f(x,y)+f(x-h,y)}{h^{2}}}\\f_{yy}(x,y)&\approx {\frac {f(x,y+k)-2f(x,y)+f(x,y-k)}{k^{2}}}\\f_{xy}(x,y)&\approx {\frac {f(x+h,y+k)-f(x+h,y-k)-f(x-h,y+k)+f(x-h,y-k)}{4hk}}.\end{aligned}}}

بدلاً من ذلك، بالنسبة للتطبيقات التي يكون فيها حساب f هو الخطوة الأكثر تكلفة، ويجب حساب كل من المشتقة الأولى والثانية، فإن الصيغة الأكثر كفاءة للحالة الأخيرة هي وxy(x،y)و(x+ح،y+ك)-و(x+ح،y)-و(x،y+ك)+2و(x،y)-و(x-ح،y)-و(x،y-ك)+و(x-ح،y-ك)2حك،{\displaystyle f_{xy}(x,y)\approx {\frac {f(x+h,y+k)-f(x+h,y)-f(x,y+k)+2f(x,y)-f(x-h,y)-f(x,y-k)+f(x-h,y-k)}{2hk}},} بما أن القيم الوحيدة التي يجب حسابها والتي ليست مطلوبة بالفعل للمعادلات الأربع السابقة هي f ( x + h , y + k ) و f ( xh , yk ) .

للوظائف ذاتشمال{\displaystyle N}المتغيراتو(x1،x2،...،xشمال){\displaystyle f(x_{1},x_{2},\ldots ,x_{N})}، وتقييم كاملم{\displaystyle m}يتطلب حساب موتر المشتقة من الرتبة n باستخدام الفروق المحدودةيا(شمالم){\displaystyle O(N^{m})}استدعاءات الدالةو{\displaystyle f}(حيث استخدمنا رمز Big O للدلالة على سلوك القياس التقاربي)، أويا(شمالم-م){\displaystyle O(N^{m-m'})}نداءاتم{\displaystyle m'}المشتقة من الرتبة n للدالةو{\displaystyle f}(أينم<م{\displaystyle m'<m}ومع ذلك، بالنسبة للعديد من فئات الدوال، فإنم{\displaystyle m}يكون موتر المشتقة من الرتبة n متفرقًا، أو قد تكون رتبته منخفضة في الكتل غير القطرية. في هذه الحالات، قد توجد خوارزميات قادرة على تقدير ذلك عدديًا.م{\displaystyle m}مشتقة من الرتبة n باستخدام أقل منيا(شمالم-م){\displaystyle O(N^{m-m'})}نداءاتم{\displaystyle m'}المشتقة من الرتبة n، على سبيل المثال عندمام=2{\displaystyle m=2}وم=1{\displaystyle m'=1}في الحالة الأخيرة، من الممكن تقدير مصفوفة هيسيان باستخدام فقطيا(1){\displaystyle O(1)}التدرجات ، بدلاً منيا(شمال){\displaystyle O(N)}التدرجات كما هو مطلوب بواسطة خوارزمية الفروق المحدودة التقليدية. [ 23 ]

انظر أيضاً

مراجع

  1. 1 2 3 بول ويلموت؛ سام هاويسون؛ جيف ديوين (1995). رياضيات المشتقات المالية: مقدمة للطلاب . مطبعة جامعة كامبريدج . ص 137. ISBN  978-0-521-49789-3.
  2. 1 2 3 بيتر أولفر (2013). مقدمة في المعادلات التفاضلية الجزئية . سبرينغر ساينس آند بيزنس ميديا. ص 182. ISBN  978-3-319-02099-0.
  3. 1 2 3 م حنيف تشودري (2007). تدفق القنوات المفتوحة . سبرينغر. ص 369. ISBN  978-0-387-68648-6.
  4. جوردان، المرجع السابق، ص 1 وميلن-تومسون، ص 21. ميلن-تومسون، لويس ميلفيل (2000): حساب الفروق المحدودة (دار نشر تشيلسي، 2000) ISBN 978-0821821077
  5. وانغ، بو؛ لو، شاوهانغ؛ وانغ، زيكوان؛ ليو، وينجيان (2025). "O1NumHess: خوارزمية هيسيان شبه عددية سريعة ودقيقة باستخدام O (1) تدرج فقط". arXiv : 2508.07544 [ physics.chem-ph ].
  6. "الفروق المحدودة لكثيرات الحدود" . divisbyzero.com . 13 فبراير 2018.
  7. فريزر، دنكان سي. (1 يناير 1909). "حول التمثيل البياني لصيغ الاستيفاء" . مجلة معهد الاكتواريين . 43 (2): 235-241 . doi : 10.1017/S002026810002494X . تاريخ الاسترجاع: 17 أبريل 2017 .
  8. ملاحظات
  9. حاسبة معاملات الفروق المحدودة
  10. بوركارد بولستر / عالم رياضيات (2021). "لماذا لا يُدرَّس حساب نيوتن للتفاضل والتكامل 'ماذا يأتي بعد ذلك؟'" على يوتيوب
  11. ^ نيوتن، إسحاق، (1687). المبادئ ، الكتاب الثالث، ليما الخامس، الحالة 1
  12. ياروسلاف ف. بلاغوشين (2018). "ثلاث ملاحظات حول تمثيلات سير وهاس لدوال زيتا" (ملف PDF) . مجلة الأعداد الصحيحة (المجلة الإلكترونية لنظرية الأعداد التوافقية) . 18A : 1-45 . arXiv : 1606.02044 . doi : 10.5281/zenodo.10581385 .
  13. كونيغ، يورغن؛ هوخت، فريد (2021). "توسيع متسلسلة نيوتن لدوال المؤثرات البوزونية" . SciPost Physics . 10 (1): 007. arXiv : 2008.11139 . Bibcode : 2021ScPP...10....7K . doi : 10.21468/SciPostPhys.10.1.007 . S2CID 221293056 . 
  14. ريختماير، د. ومورتون، ك.و.، (1967). طرق الفرق لمسائل القيمة الأولية ، الطبعة الثانية، وايلي، نيويورك.
  15. بول، جورج (1872). رسالة في حساب الفروق المحدودة ( الطبعة الثانية). ماكميلان وشركاه عبر أرشيف الإنترنت . كما صدرت طبعة معاد طباعتها من دار نشر دوفر عام 1960.
  16. جوردان، تشارلز (1965) [1939]. حساب التفاضل والتكامل للفروق المحدودة . دار نشر تشيلسي. ISBN 978-0-8284-0033-6 عبر كتب جوجل.
  17. زاكوس، سي. (2008). "التشوهات الظلية في الزمكان المنفصل". المجلة الدولية للفيزياء الحديثة أ . 23 (13): 200-214 . arXiv : 0710.2306 . Bibcode : 2008IJMPA..23.2005Z . doi : 10.1142/S0217751X08040548 . S2CID 16797959 . 
  18. كورترايت، تي إل؛ زاكوس، سي كيه (2013). "دليل الظل" . حدود في الفيزياء . 1 : 15. arXiv : 1304.0429 . Bibcode : 2013FrP.....1...15C . doi : 10.3389/fphy.2013.00015 . S2CID 14106142 . 
  19. ليفي، هـ.؛ ليسمان، ف. (1992). معادلات الفروق المحدودة . دوفر. ISBN 0-486-67260-3.
  20. أميس، دبليو إف (1977). الطرق العددية للمعادلات التفاضلية الجزئية . نيويورك، نيويورك: أكاديميك برس. القسم 1.6. ISBN  0-12-056760-1.
  21. هيلدبراند، إف بي (1968). معادلات الفروق المحدودة والمحاكاة . إنجلوود كليفس، نيوجيرسي: برنتيس هول. القسم 2.2. 
  22. فلاجويه، فيليب؛ سيدجويك، روبرت (1995). "تحويلات ميلين والتقارب: الفروق المحدودة وتكاملات رايس" (ملف PDF) . علوم الحاسوب النظرية . 144 ( 1-2 ): 101-124 . doi : 10.1016/0304-3975(94)00281-M .
  23. وانغ، بو؛ لو، شاوهانغ؛ وانغ، زيكوان؛ ليو، وينجيان (2025). "O1NumHess: خوارزمية هيسيان شبه عددية سريعة ودقيقة باستخدام تدرجات O(1) فقط". arXiv : 2508.07544 [ physics.chem-ph ].
  • ريتشاردسون، سي إتش (1954): مقدمة في حساب الفروق المحدودة ( نسخة فان نوستراند الإلكترونية (1954))
  • ميكنز، ر. إي. (1991): المعادلات التفاضلية: النظرية والتطبيقات (تشابمان آند هول/سي آر سي) ISBN 978-0442001360