خوارزمية نيفيل

في الرياضيات ، تُعدّ خوارزمية نيفيل خوارزميةً تُستخدم لاستيفاء كثيرات الحدود، وقد اشتقّها عالم الرياضيات إريك هارولد نيفيل عام 1934. بمعلومية n + 1 نقطة، يوجد كثير حدود وحيد من الدرجة ≤ n يمرّ بهذه النقاط. تقوم خوارزمية نيفيل بتقييم هذا كثير الحدود.

تعتمد خوارزمية نيفيل على صيغة نيوتن لكثير الحدود الاستيفائي وعلاقة التكرار للفروق المقسمة . وهي مشابهة لخوارزمية أيتكن (نسبةً إلى ألكسندر أيتكن )، والتي لم تعد مستخدمة في الوقت الحاضر.

الخوارزمية

بفرض وجود مجموعة من n + 1 نقطة بيانات ( xᵢ , yᵢ ) حيث لا توجد نقطتان xᵢ متطابقتان، فإن متعددة الحدود الاستيفائية هي متعددة الحدود p من الدرجة n على الأكثر والتي تتمتع بالخاصية التالية :

p ( xi ) = yi لجميع قيم i من 0 إلى n .

هذه المعادلة متعددة الحدود موجودة وفريدة من نوعها. تقوم خوارزمية نيفيل بتقييم المعادلة متعددة الحدود عند نقطة ما x .

لنفترض أن pᵢ, يمثل متعدد الحدود من الدرجة j i الذي يمر بالنقاط ( xⱼ , yⱼ ) حيث k = i , i + 1, ..., j . ويحقق pᵢ , العلاقة التكرارية التالية

صأنا،أنا(x)=yأنا،{\displaystyle p_{i,i}(x)=y_{i},\,}0أنان،{\displaystyle 0\leq i\leq n,\,}
صأنا،ج(x)=(x-xأنا)صأنا+1،ج(x)-(x-xج)صأنا،ج-1(x)xج-xأنا،{\displaystyle p_{i,j}(x)={\frac {(x-x_{i})p_{i+1,j}(x)-(x-x_{j})p_{i,j-1}(x)}{x_{j}-x_{i}}},\,}0أنا<جن.{\displaystyle 0\leq i<j\leq n.\,}

يمكن لهذه العلاقة التكرارية حساب p 0, n ( x )، وهي القيمة المطلوبة. هذه هي خوارزمية نيفيل.

على سبيل المثال، بالنسبة لـ n = 4، يمكن للمرء استخدام التكرار لملء الجدول المثلثي أدناه من اليسار إلى اليمين.

ص0،0(x)=y0{\displaystyle p_{0,0}(x)=y_{0}\,}
ص0،1(x){\displaystyle p_{0,1}(x)\,}
ص1،1(x)=y1{\displaystyle p_{1,1}(x)=y_{1}\,}ص0،2(x){\displaystyle p_{0,2}(x)\,}
ص1،2(x){\displaystyle p_{1,2}(x)\,}ص0،3(x){\displaystyle p_{0,3}(x)\,}
ص2،2(x)=y2{\displaystyle p_{2,2}(x)=y_{2}\,}ص1،3(x){\displaystyle p_{1,3}(x)\,}ص0،4(x){\displaystyle p_{0,4}(x)\,}
ص2،3(x){\displaystyle p_{2,3}(x)\,}ص1،4(x){\displaystyle p_{1,4}(x)\,}
ص3،3(x)=y3{\displaystyle p_{3,3}(x)=y_{3}\,}ص2،4(x){\displaystyle p_{2,4}(x)\,}
ص3،4(x){\displaystyle p_{3,4}(x)\,}
ص4،4(x)=y4{\displaystyle p_{4,4}(x)=y_{4}\,}

ينتج عن هذه العملية p 0,4 ( x )، قيمة متعددة الحدود التي تمر عبر نقاط البيانات n + 1 ( x i , y i ) عند النقطة x .

تحتاج هذه الخوارزمية إلى O ( n 2 ) عملية حسابية للفاصلة العائمة لاستكمال نقطة واحدة، و O ( n 3 ) عملية حسابية للفاصلة العائمة لاستكمال متعدد الحدود من الدرجة n.

يمكن الحصول على مشتقة متعددة الحدود بنفس الطريقة، أي:

صأنا،أنا(x)=0،{\displaystyle p'_{i,i}(x)=0,\,}0أنان،{\displaystyle 0\leq i\leq n,\,}
صأنا،ج(x)=(x-xأنا)صأنا+1،ج(x)+صأنا+1،ج(x)-(x-xج)صأنا،ج-1(x)-صأنا،ج-1(x)xج-xأنا،{\displaystyle p'_{i,j}(x)={\frac {(x-x_{i})p'_{i+1,j}(x)+p_{i+1,j}(x)-(x-x_{j})p'_{i,j-1}(x)-p_{i,j-1}(x)}{x_{j}-x_{i}}},\,}0أنا<جن.{\displaystyle 0\leq i<j\leq n.\,}

تدوين بديل أسهل للتطبيق الحاسوبي

في الصيغ أعلاه، إذا أخذنا درجة كثيرات الحدود الاستيفائية المتتالية d = j i وغيرنا الترميز إلى p d , i ,

ص0،أنا(x)=yأنا،{\displaystyle p_{0,i}(x)=y_{i},}د=0{\displaystyle d=0}
صد،أنا(x)=(x-xأنا)صد-1،أنا+1(x)-(x-xأنا+د)صد-1،أنا(x)xأنا+د-xأنا،{\displaystyle p_{d,i}(x)={\frac {(x-x_{i})p_{d-1,i+1}(x)-(x-x_{i+d})p_{d-1,i}(x)}{x_{i+d}-x_{i}}},}1دن،0أنان-د{\displaystyle 1\leq d\leq n,0\leq i\leq nd}

القيمة النهائية p n ,0 (في هذه الصيغة) هي القيمة المطلوبة التي تم استكمالها.

بما أن عدد العناصر المحسوبة، أي نطاق يتناقص مع كل قيمة d متتالية ، يمكن استخدام مصفوفة خطية لتحسين كفاءة الذاكرة، حيث يتم استبدال p i وتجاهل d . (على سبيل المثال:)

يمكن حساب المشتقة (باستخدام قاعدة الضرب ) بنفس الطريقة كما يلي:

ص0،أنا(x)=0،{\displaystyle p'_{0,i}(x)=0,}د=0{\displaystyle d=0}
صد،أنا(x)=(x-xأنا)صد-1،أنا+1(x)+صد-1،أنا+1(x)-(x-xأنا+د)صد-1،أنا(x)-صد-1،أنا(x)xأنا+د-xأنا،{\displaystyle p'_{d,i}(x)={\frac {(x-x_{i})p'_{d-1,i+1}(x)+p_{d-1,i+1}(x)-(x-x_{i+d})p'_{d-1,i}(x)-p_{d-1,i}(x)}{x_{i+d}-x_{i}}},}1دن،0أنان-د{\displaystyle 1\leq d\leq n,0\leq i\leq nd}

كما في السابق، فإن pn ,0 (في هذه الصيغة) هو المشتق.

بما أن هذا يعتمد على القيم المحسوبة المتتالية لـ p لكل d ، فإنه يمكن حسابه ضمن نفس الحلقة. إذا تم استخدام مصفوفات خطية لـ p و p′ لتحسين الكفاءة، فيجب حساب قيم p ′ قبل استبدال قيم p .

تطبيق على التفاضل العددي

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

مراجع