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

الفرق المحدود هو تعبير رياضي على الصورة 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).}

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

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

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

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

The Newton series consists of the terms of the Newton forward difference equation, named after Isaac Newton; in essence, it is the Gregory–Newton interpolation formula[10] (named after Isaac Newton and James Gregory), first published in his Principia Mathematica in 1687,[11][12] namely the discrete analog of the continuous Taylor expansion,

f(x)=k=0Δk[f](a)k!(xa)k=k=0(xak)Δk[f](a),{\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),}

which holds for any polynomial function f and for many (but not all) analytic functions. (It does not hold when f is exponential typeπ{\displaystyle \pi }. This is easily seen, as the sine function vanishes at integer multiples of π{\displaystyle \pi }; the corresponding Newton series is identically zero, as all finite differences are zero in this case. Yet clearly, the sine function is not zero.) Here, the expression (xk)=(x)kk!{\displaystyle {\binom {x}{k}}={\frac {(x)_{k}}{k!}}} is the binomial coefficient, and (x)k=x(x1)(x2)(xk+1){\displaystyle (x)_{k}=x(x-1)(x-2)\cdots (x-k+1)} is the "falling factorial" or "lower factorial", while the empty product(x)0 is defined to be 1. In this particular case, there is an assumption of unit steps for the changes in the values of x, h = 1 of the generalization below.

Note the formal correspondence of this result to Taylor's theorem. Historically, this, as well as the Chu–Vandermonde identity, (x+y)n=k=0n(nk)(x)nk(y)k,{\displaystyle (x+y)_{n}=\sum _{k=0}^{n}{\binom {n}{k}}(x)_{n-k}\,(y)_{k},} (following from it, and corresponding to the binomial theorem), are included in the observations that matured to the system of umbral calculus.

Newton series expansions can be superior to Taylor series expansions when applied to discrete quantities like quantum spins (see Holstein–Primakoff transformation), bosonic operator functions or discrete counting statistics.[13]

To illustrate how one may use Newton's formula in actual practice, consider the first few terms of doubling the Fibonacci sequencef = 2, 2, 4, ... One can find a polynomial that reproduces these values, by first computing a difference table, and then substituting the differences that correspond to x0 (underlined) into the formula as follows, xf=Δ0Δ1Δ212_0_222_234f(x)=Δ01+Δ1(xx0)11!+Δ2(xx0)22!(x0=1)=21+0x11+2(x1)(x2)2=2+(x1)(x2){\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}}}

For the case of nonuniform steps in the values of x, Newton computes the divided differences, Δj,0=yj,Δj,k=Δj+1,k1Δj,k1xj+kxj{k>0,jmax(j)k},Δ0k=Δ0,k{\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}} the series of products, P0=1,Pk+1=Pk(ξxk),{\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