Regula falsi

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

نوعان تاريخيان

يمكن التمييز تاريخياً بين نوعين أساسيين من طريقة الموقف الخاطئ، وهما الموقف الخاطئ البسيط والموقف الخاطئ المزدوج .

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

أx=ب،{\displaystyle ax=b,}

إذا كانت قيمتا a و b معلومتين، تبدأ الطريقة باستخدام قيمة إدخال اختبارية x ، ثم إيجاد قيمة الإخراج المقابلة b عن طريق الضرب: ax ′ = b ′ . بعد ذلك ، يتم إيجاد الإجابة الصحيحة عن طريق التعديل النسبي، x = b / b x .

كمثال، لننظر إلى المسألة رقم 26 في بردية رايند ، والتي تطلب حل المعادلة (المكتوبة بالصيغة الحديثة) x + x / 4 = 15. يتم حل هذه المسألة باستخدام طريقة الوضع الخاطئ. [ 1 ] أولًا ، نفترض أن x = 4 لنحصل على المعادلة 4 + 4 / 4 = 5. هذا الافتراض صحيح لأنه يُعطي قيمة صحيحة . مع ذلك، فإن 4 ليس حلًا للمعادلة الأصلية، لأنه يُعطي قيمة أصغر بثلاث مرات من القيمة الصحيحة. للتعويض، نضرب x ( المُعيّنة حاليًا على 4) في 3 ونُعوّض مرة أخرى لنحصل على 12 + 12 / 4 = 15 ، مما يُؤكد أن الحل هو x = 12 .

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

و(x)=أx+ج=0،{\displaystyle f(x)=ax+c=0,}

إذا كان معروفاً أن

و(x1)=ب1؛و(x2)=ب2.{\displaystyle {\begin{aligned}f(x_{1})&=b_{1};\\f(x_{2})&=b_{2}.\end{aligned}}}

يُعدّ وضع الخطأ المزدوج مكافئًا رياضيًا للاستيفاء الخطي . باستخدام زوج من مدخلات الاختبار وزوج المخرجات المقابل، تكون نتيجة هذه الخوارزمية كما يلي: [ 2 ]

x=ب1x2-ب2x1ب1-ب2،{\displaystyle x={\frac {b_{1}x_{2}-b_{2}x_{1}}{b_{1}-b_{2}}},}

سيتم حفظها وتنفيذها عن ظهر قلب. في الواقع، القاعدة كما وردت في كتاب روبرت ريكورد " أساسيات الفنون " (حوالي 1542) هي: [ 2 ]

اجْتَهْدَمْ فِي هَذَا الْعَمَلِ فَإِنَّ الْحَقَّ قَدْ تَقْدِمُ. بِالْحَقِّ قَدْ تَسْبِعُ الْحَقَّ. وَاعْمَلْ أَوْ افْعَلَ الْأَوَّلِ بِالْأَسئلِ، حَتَّى لَمْ يَكُونَ الْحَقَّ فِيهِ. كَانَ الْبَاطِلُ مَكَانًا جَيِّدًا، وَبِالْأَحْصَادِ مِنْهُ سَيُوجَدُ الْحَقَّ بِهِ سَيُوجَدُ الْحَقَّ بِهِ بَعْضًا. مِنَ الْبَاطِلِ يُوجَدُّ بَعْضًا، مِنَ الْبَاطِلِ يُوجَدُّ بَعْضًا. بِالْأَحْصَادِ تَتَعَالَجُ الْأَضْلِيَّاتِ، كُلُّ الحَقَّ بَاطِلٍ لِيُوجَدُ الْحَقَّ.

بالنسبة للدالة الخطية الأفينية ،

و(x)=أx+ج،{\displaystyle f(x)=ax+c,}

يوفر الوضع الخاطئ المزدوج الحل الدقيق، بينما بالنسبة للدالة غير الخطية f فإنه يوفر تقريبًا يمكن تحسينه بشكل متتابع عن طريق التكرار .

تاريخ

توجد تقنية الوضع الخاطئ البسيطة في الألواح المسمارية من الرياضيات البابلية القديمة ، وفي البرديات من الرياضيات المصرية القديمة . [ 3 ] [ 1 ]

ظهرت مغالطة الموقف المزدوج الخاطئ في أواخر العصور القديمة كخوارزمية حسابية بحتة. في النص الرياضي الصيني القديم المسمى " الفصول التسعة في الفن الرياضي" (九章算術)، [ 4 ] والذي يعود تاريخه إلى الفترة ما بين 200 قبل الميلاد و100 ميلادي، خُصص معظم الفصل السابع لهذه الخوارزمية. هناك، تم تبرير الإجراء بحجج حسابية ملموسة، ثم طُبق بشكل إبداعي على مجموعة واسعة من المسائل الكلامية، بما في ذلك مسألة تتضمن ما نسميه اليوم خطوط القاطع على القطع المخروطي . ومن الأمثلة الأكثر شيوعًا مسألة "الشراء المشترك" هذه التي تتضمن شرط "الفائض والعجز": [ 5 ]

الآن، يتم شراء سلعة بشكل جماعي؛ يساهم كل شخص بـ 8 عملات معدنية، والزيادة 3 عملات؛ يساهم كل شخص بـ 7 عملات معدنية، والنقص 4 عملات. ما هو عدد الأشخاص، وسعر السلعة، وما هو سعر كل منهم؟ الجواب: 7 أشخاص، سعر السلعة 53. [ 6 ]

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

خصّص ليوناردو البيزي ( فيبوناتشي ) الفصل الثالث عشر من كتابه "ليبر أباتشي " (1202م) لشرح وتوضيح استخدامات الوضع الخاطئ المزدوج، وأطلق على هذه الطريقة اسم "ريجوليس إلشاتين" نسبةً إلى طريقة "الخطائين" التي تعلمها من مصادر عربية . [ 7 ] وفي عام 1494م، استخدم باتشولي مصطلح "إل كاتايم" في كتابه "سوما دي أريثميتيكا" ، وربما استقاه من فيبوناتشي. وحذا حذوه كتّاب أوروبيون آخرون، وقدّم بعضهم ترجمةً إلى اللاتينية أو اللغة العامية. فعلى سبيل المثال، ترجم تارتاليا النسخة اللاتينية من مصطلح باتشولي إلى اللغة العامية "الأوضاع الخاطئة" عام 1556م. [ 8 ] وكاد مصطلح باتشولي أن يختفي في الأعمال الأوروبية في القرن السادس عشرم، وأصبحت هذه التقنية تُعرف بأسماء مختلفة مثل "قاعدة الخطأ" و"قاعدة الوضع" و"قاعدة الوضع الخاطئ". يظهر مصطلح Regula Falsi كنسخة لاتينية من Rule of False في وقت مبكر من عام 1690. [ 2 ]

شعر العديد من المؤلفين الأوروبيين في القرن السادس عشر بالحاجة إلى الاعتذار عن اسم المنهج في علم يسعى إلى إيجاد الحقيقة. على سبيل المثال، يقول همفري بيكر في عام 1568: [ 2 ]

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

التحليل العددي

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

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

بفرض معادلة، انقل جميع حدودها إلى طرف واحد بحيث تصبح على الصورة f ( x ) = 0 ، حيث f دالة ما للمتغير المجهول x . تُسمى القيمة c التي تحقق هذه المعادلة، أي f ( c ) = 0 ، جذرًا أو صفرًا للدالة f ، وهي حل للمعادلة الأصلية. إذا كانت f دالة متصلة ، ووجدت نقطتان a₀ و b₀ بحيث تكون f ( a₀ ) و f ( b₀ ) متعاكستين في الإشارة، فإنه وفقًا لنظرية القيمة المتوسطة ، يكون للدالة f جذر في الفترة ( a₀ , b₀ ) .

توجد العديد من خوارزميات إيجاد الجذور التي يمكن استخدامها للحصول على تقريبات للجذر. إحدى أكثرها شيوعًا هي طريقة نيوتن ، ولكنها قد تفشل في إيجاد الجذر في بعض الحالات، وقد تكون مكلفة حسابيًا لأنها تتطلب حساب مشتقة الدالة . لذا، ثمة حاجة إلى طرق أخرى، ومنها طرق التحديد بنقطتين . تعمل هذه الطرق عن طريق إنتاج سلسلة من الفترات المتناقصة [ a <sub>k</sub> , b<sub> k</sub> ] ، في الخطوة k ، بحيث تحتوي الفترة ( a <sub>k</sub> , b<sub> k</sub> ) على جذر للدالة f .

طرق تحديد المواقع بنقطتين

تبدأ هذه الطرق بقيمتين لـ x ، يتم تحديدهما مبدئيًا بالتجربة والخطأ، حيث تكون إشارة f ( x ) متعاكسة. بافتراض استمرارية الدالة، يُضمن وجود جذر للدالة f بين هاتين القيمتين، أي أن هاتين القيمتين تُحيطان بالجذر. ثم تُختار نقطة تقع تمامًا بين هاتين القيمتين وتُستخدم لإنشاء فترة أصغر تُحيط بجذر. إذا كانت النقطة المختارة هي c ، فإن الفترة الأصغر تمتد من c إلى نقطة النهاية حيث تكون إشارة f ( x ) معاكسة لإشارة f ( c ) . في الحالة النادرة التي تكون فيها f ( c ) = 0 ، يكون قد تم العثور على جذر وتتوقف الخوارزمية. وإلا، تُكرر العملية حسب الحاجة للحصول على تقريب للجذر بالدقة المطلوبة.

يمكن اعتبار النقطة المختارة في أي فترة زمنية حالية بمثابة تقدير للحل. وتتضمن الاختلافات في هذه الطريقة أساليب مختلفة لحساب هذا التقدير للحل.

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

أبسط طريقة، وتُسمى طريقة التنصيف ، تحسب تقدير الحل كنقطة منتصف فترة التحديد. أي، إذا كانت فترة التحديد الحالية في الخطوة k هي [ a k , b k ] ، فإن تقدير الحل الجديد c k يُحسب كما يلي:

جك=أك+بك2.{\displaystyle c_{k}={\frac {a_{k}+b_{k}}{2}}.}

وهذا يضمن أن يكون c k بين a k و b k ، مما يضمن التقارب نحو الحل.

بما أن طول فترة التحديد يتقلص إلى النصف في كل خطوة، فإن خطأ طريقة التنصيف يتقلص، في المتوسط، إلى النصف مع كل تكرار. وبالتالي، كل 3 تكرارات، تكتسب الطريقة دقةً تقارب 2³ ، أي ما يعادل منزلة عشرية تقريبًا. يُشار إلى هذا عادةً بالتقارب من الدرجة الأولى، أي أن عدد خانات الدقة يتناسب طرديًا مع عدد التكرارات المستخدمة.

طريقة regula falsi (الموضع الزائف).

أول دورتين من طريقة الوضع الخاطئ. يُظهر المنحنى الأحمر الدالة f ، بينما تُمثل الخطوط الزرقاء القواطع.

من الممكن تحسين معدل تقارب طريقة التنصيف باستخدام تقدير حل مختلف.

تحسب طريقة "الوضع الخاطئ المزدوج" تقدير الحل الجديد كنقطة تقاطع القطعة المستقيمة الواصلة بين طرفي الدالة على الفترة المحددة مع محور السينات. وباختصار، يتم تقريب الجذر باستبدال الدالة الأصلية بقطعة مستقيمة على الفترة المحددة، ثم تطبيق صيغة "الوضع الخاطئ المزدوج" الكلاسيكية على تلك القطعة المستقيمة. [ 9 ]

بتعبير أدق، لنفترض أن الفترة المحددة في التكرار رقم k هي ( a <sub>k</sub> , b<sub> k</sub> ) . ارسم الخط المار بالنقطتين ( a <sub>k</sub> , f ( a<sub> k</sub> )) و ( b <sub>k</sub> , f ( b<sub> k</sub> )) ، كما هو موضح. هذا الخط هو قاطع أو وتر لمنحنى الدالة f . معادلته، بصيغة الميل والنقطة ، هي:

y-و(بك)=و(بك)-و(أك)بك-أك(x-بك).{\displaystyle yf(b_{k})={\frac {f(b_{k})-f(a_{k})}{b_{k}-a_{k}}}(x-b_{k}).}

اختر الآن c k لتكون نقطة تقاطع هذا الخط مع المحور السيني، أي قيمة x التي عندها y = 0 ، ثم استبدل هذه القيم للحصول على

و(بك)+و(بك)-و(أك)بك-أك(جك-بك)=0.{\displaystyle f(b_{k})+{\frac {f(b_{k})-f(a_{k})}{b_{k}-a_{k}}}(c_{k}-b_{k})=0.}

حل هذه المعادلة لإيجاد قيمة c k يعطينا:

جك=بك-و(بك)بك-أكو(بك)-و(أك)=أكو(بك)-بكو(أك)و(بك)-و(أك).{\displaystyle c_{k}=b_{k}-f(b_{k}){\frac {b_{k}-a_{k}}{f(b_{k})-f(a_{k})}}={\frac {a_{k}f(b_{k})-b_{k}f(a_{k})}{f(b_{k})-f(a_{k})}}.}

يتميز هذا الشكل المتناظر الأخير بميزة حسابية عند استخدام العمليات الحسابية ذات الفاصلة العائمة : فمع اقتراب الحل، تتقارب قيمتا a <sub>k</sub> و b<sub> k</sub> بشكل كبير، وتكونان دائمًا تقريبًا بنفس الإشارة. قد يؤدي هذا النوع من الطرح إلى فقدان الدقة بسبب الاختزال . ولأن f ( b<sub> k</sub> ) و f ( a<sub> k</sub> ) دائمًا بإشارتين متعاكستين، فإن "الطرح" في بسط الصيغة المحسّنة هو في الواقع عملية جمع (وكذلك الطرح في المقام).

في التكرار رقم k ، تُحسب قيمة c<sub> k </sub> كما سبق، ثم إذا كانت f ( ak ) و f ( ck ) لهما نفس الإشارة، نُعيّن a <sub>k +1</sub> = c<sub> k</sub> و b <sub>k +1 </sub> = b <sub> k </sub> ، وإلا نُعيّن a <sub>k +1</sub> = a <sub>k</sub> و b <sub> k +1 </sub> = c <sub> k</sub> . تُكرر هذه العملية حتى يتم تقريب الجذر بدقة كافية. تُستخدم الصيغة المذكورة أعلاه أيضًا في طريقة القاطع .

بالنسبة للدوال غير الخطية، عندما يتقلص نطاق البحث إلى حدٍّ كافٍ بحيث تصبح إشارة المشتقة الثانية ثابتة في جميع أنحاء النطاق، تصبح إحدى نقطتي نهاية البحث ثابتة، بينما تتقارب الأخرى نحو الجذر. وبالتالي، فإن أفضل تقدير للحل هو آخر قيمة محسوبة لـجك{\displaystyle c_{k}}مع ذلك، ولأن الفترة تتوقف عن التقلص، لا تستطيع طريقة "النمط الخاطئ" أن تضاهي ضمان دقة طريقة التنصيف. في بعض الحالات، قد ينخفض ​​معدل التقارب عن معدل طريقة التنصيف. يُفضّل عمومًا استخدام نسخ مُعدّلة من " النمط الخاطئ" لأنها تُعالج هذه العيوب بأقل تكلفة.

تحليل

بما أن نقطتي النهاية الابتدائيتين a₀ و b₀ تُختاران بحيث تكون f ( a₀ ) و f ( b₀ ) متعاكستين في الإشارة ، ففي كل خطوة، تقترب إحدى نقطتي النهاية من جذر الدالة f . إذا كانت المشتقة الثانية للدالة f ذات إشارة ثابتة (أي لا توجد نقطة انعطاف ) في الفترة، فإن إحدى نقطتي النهاية (التي تكون عندها f بنفس الإشارة) ستظل ثابتة في جميع التكرارات اللاحقة، بينما يتم تحديث نقطة النهاية المتقاربة. ونتيجة لذلك، وعلى عكس طريقة التنصيف ، لا يؤول عرض القوس إلى الصفر (إلا إذا كان الصفر عند نقطة انعطاف يكون عندها sign( f ) = −sign( f '' ) ). وبالتالي، فإن التقريب الخطي للدالة f ( x ) ، المستخدم لاختيار الموضع الخاطئ، لا يتحسن بالسرعة المثلى.

أحد الأمثلة على هذه الظاهرة هو الدالة

و(x)=2x3-4x2+3x{\displaystyle f(x)=2x^{3}-4x^{2}+3x}

على القوس الأولي [ −1,1 ]. الطرف الأيسر، −1 ، لا يتم استبداله أبدًا (لا يتغير في البداية وبعد التكرارات الثلاثة الأولى، تكون f " سالبة على الفترة) وبالتالي فإن عرض القوس لا يقل أبدًا عن 1. ومن ثم، تقترب نقطة النهاية اليمنى من 0 بمعدل خطي (يزداد عدد الأرقام الصحيحة خطيًا، بمعدل تقارب 2/3).

بالنسبة للدوال غير المتصلة، لا يُتوقع من هذه الطريقة سوى إيجاد نقطة تغير عندها الدالة إشارتها (على سبيل المثال، عند x = 0 للدالة 1/ x أو دالة الإشارة ). بالإضافة إلى تغيرات الإشارة، من الممكن أيضًا أن تتقارب الطريقة إلى نقطة تكون عندها نهاية الدالة صفرًا، حتى لو كانت الدالة غير مُعرَّفة (أو لها قيمة أخرى) عند تلك النقطة (على سبيل المثال، عند x = 0 للدالة المعطاة بالصيغة f ( x ) = abs( x ) − عندما x ≠ 0، و f (0) = 5 ، بدءًا من الفترة [-0.5، 3.0]). من الممكن رياضيًا، مع الدوال غير المتصلة ، أن تفشل الطريقة في التقارب إلى نهاية صفرية أو تغيير في الإشارة، لكن هذا لا يُشكِّل مشكلة عملية، إذ يتطلب ذلك سلسلة لانهائية من التطابقات لكلا نقطتي النهاية حتى تتعثر في التقارب إلى نقاط عدم اتصال حيث لا تتغير الإشارة، على سبيل المثال عند x = ±1 .

و(x)=1(x-1)2+1(x+1)2.{\displaystyle f(x)={\frac {1}{(x-1)^{2}}}+{\frac {1}{(x+1)^{2}}}.}

تتجنب طريقة التنصيف مشكلة التقارب الافتراضية هذه.

تحسينات في regula falsi

على الرغم من أن طريقة "النمط الكاذب" تتقارب دائمًا، وعادةً ما يكون التقارب أسرع بكثير من طريقة "التنصيف"، إلا أن هناك حالات قد تُبطئ تقاربها، وأحيانًا إلى درجةٍ تُعيق ذلك. هذه المشكلة ليست حكرًا على طريقة " النمط الكاذب" : فباستثناء طريقة "التنصيف"، قد تُعاني جميع طرق حل المعادلات العددية من مشكلة بطء التقارب أو عدم التقارب في ظل ظروفٍ معينة. في بعض الأحيان، تتباعد طريقة نيوتن وطريقة القاطع بدلًا من التقارب، وغالبًا ما يحدث ذلك في ظل نفس الظروف التي تُبطئ تقارب طريقة "النمط الكاذب" .

لكن، على الرغم من أن regula falsi هي واحدة من أفضل الطرق، وحتى في نسختها الأصلية غير المحسنة، فإنها غالبًا ما تكون الخيار الأفضل؛ على سبيل المثال، عندما لا يتم استخدام طريقة نيوتن لأن حساب المشتقة يستغرق وقتًا طويلاً للغاية، أو عندما تفشل طريقة نيوتن والاستبدالات المتتالية في التقارب.

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

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

تقوم خوارزمية إلينوي بتقسيم قيمة y لنقطة النهاية المحفوظة إلى النصف في حساب التقدير التالي عندما تكون قيمة y الجديدة (أي f ( ck ) ) لها نفس إشارة القيمة السابقة ( f ( ck - 1 ))، مما يعني أنه سيتم الاحتفاظ بنقطة نهاية الخطوة السابقة. وبالتالي:

جك=12و(بك)أك-و(أك)بك12و(بك)-و(أك){\displaystyle c_{k}={\frac {{\frac {1}{2}}f(b_{k})a_{k}-f(a_{k})b_{k}}{{\frac {1}{2}}f(b_{k})-f(a_{k})}}}

أو

جك=و(بك)أك-12و(أك)بكو(بك)-12و(أك)،{\displaystyle c_{k}={\frac {f(b_{k})a_{k}-{\frac {1}{2}}f(a_{k})b_{k}}{f(b_{k})-{\frac {1}{2}}f(a_{k})}},}

يتم تقليل وزن إحدى قيم نقطة النهاية لإجبار الدالة التالية على الظهور على ذلك الجانب من الدالة. [ 10 ] قد يبدو العامل 1/2 المستخدم أعلاه عشوائيًا، ولكنه يضمن تقاربًا فائق الخطية (بشكل تقاربي، ستنفذ الخوارزمية خطوتين عاديتين بعد أي خطوة معدلة، ولها رتبة تقارب 1.442 ) . توجد طرق أخرى لاختيار إعادة التحجيم تُعطي معدلات تقارب فائق الخطية أفضل. [ 11 ]

يُطلق بعض الباحثين على التعديل المذكور أعلاه لـ regula falsi اسم خوارزمية إلينوي . [ 10 ] [ 12 ] وقد لخص فورد (1995) هذه الخوارزمية وغيرها من المتغيرات فوق الخطية المماثلة لطريقة الوضع الخاطئ وحللها. [ 11 ]

خوارزمية أندرسون-بيورك

لنفترض أنه في التكرار رقم تكون فترة التحديد [ ak , bk ] ، وأن القيمة الوظيفية للتقدير المحسوب الجديد ck لها نفس إشارة f ( bk ) . في هذه الحالة، تكون فترة التحديد الجديدة [ ak +1 , bk +1 ] = [ ak , ck ] ، وقد تم الاحتفاظ بنقطة النهاية اليسرى. (حتى الآن، هذا يُشابه خوارزمية Regula Falsi العادية وخوارزمية إلينوي).

لكن، في حين أن خوارزمية إلينوي تضرب f ( a k ) في 1 / 2 ، فإن خوارزمية أندرسون-بيورك تضربها في m ، حيث m لها إحدى القيمتين التاليتين: [ 13 ]

م=1-و(جك)و(بك)،م={ملو م>0،12خلاف ذلك.{\displaystyle {\begin{aligned}m'&=1-{\frac {f(c_{k})}{f(b_{k})}},\\m&={\begin{cases}m'&{\text{if }}m'>0,\\{\frac {1}{2}}&{\text{otherwise.}}\end{cases}}\end{aligned}}}

بالنسبة للجذور البسيطة، فإن طريقة أندرسون-بيورك تؤدي أداءً جيداً للغاية من الناحية العملية. [ 14 ]

طريقة ITP

منحκ1(0،)،κ2[1،1+ϕ){\displaystyle \kappa _{1}\in (0,\infty ),\kappa _{2}\in \left[1,1+\phi \right)}،ن1/2(ب0-أ0)/2ϵ{\displaystyle n_{1/2}\equiv \lceil (b_{0}-a_{0})/2\epsilon \rceil } ون0[0،){\displaystyle n_{0}\in [0,\infty )} أينϕ{\displaystyle \phi }النسبة الذهبية12(1+5){\displaystyle {\tfrac {1}{2}}(1+{\sqrt {5}})}، في كل تكرارج=0،1،2...{\displaystyle j=0,1,2...} تقوم طريقة ITP بحساب النقطةxنقص الصفيحات المناعي{\displaystyle x_{\text{ITP}}}الخطوات الثلاث التالية:

  1. [خطوة الاستيفاء] حساب نقاط التنصيف ونقاط الانحراف الكاذب: x1/2أ+ب2{\displaystyle x_{1/2}\equiv {\frac {a+b}{2}}} و xوبو(أ)-أو(ب)و(أ)-و(ب){\displaystyle x_{f}\equiv {\frac {bf(a)-af(b)}{f(a)-f(b)}}} ؛
  2. [خطوة الاقتطاع] قم بتحريك المُقدِّر نحو المركز: xتxو+σدلتا{\displaystyle x_{t}\equiv x_{f}+\sigma \delta } أين σلافتة(x1/2-xو){\displaystyle \sigma \equiv {\text{sign}}(x_{1/2}-x_{f})}ودلتامين{κ1|ب-أ|κ2،|x1/2-xو|}{\displaystyle \delta \equiv \min\{\kappa _{1}|b-a|^{\kappa _{2}},|x_{1/2}-x_{f}|\}} ؛
  3. [خطوة الإسقاط] إسقاط المُقدِّر على فاصل minmax:xنقص الصفيحات المناعيx1/2-σρك{\displaystyle x_{\text{ITP}}\equiv x_{1/2}-\sigma \rho _{k}}أينρكمين{ϵ2ن1/2+ن0-ج-ب-أ2،|xت-x1/2|}{\displaystyle \rho _{k}\equiv \min \left\{\epsilon 2^{n_{1/2}+n_{0}-j}-{\frac {b-a}{2}},|x_{t}-x_{1/2}|\right\}}.

قيمة الدالةو(xنقص الصفيحات المناعي){\displaystyle f(x_{\text{ITP}})}يتم الاستفسار عن هذه النقطة، ثم تُختزل الفترة لتشمل الجذر مع الاحتفاظ بالفترة الفرعية التي تحتوي على قيم دالة ذات إشارات معاكسة على كل طرف. يضمن هذا الإجراء المكون من ثلاث خطوات أن يتمتع التقدير بخصائص الحد الأدنى والحد الأقصى لطريقة التنصيف، بالإضافة إلى التقارب فوق الخطي لطريقة القاطع. وقد لوحظ تفوقها على كل من طرق التنصيف وطرق الاستيفاء في حالة الدوال الملساء وغير الملساء. [ 15 ]

الاعتبارات العملية

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

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

ويُستثنى من ذلك حالة اضطرار برنامج الحاسوب إلى حل المعادلات مرات عديدة أثناء تشغيله. فحينها قد يكون الوقت المُوفّر باستخدام الطرق الأسرع كبيراً.

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

عندما يصبح التغير في قيمة y ضئيلاً للغاية، ويتغير x أيضاً بشكل طفيف جداً، فمن المرجح ألا تواجه طريقة نيوتن أي مشكلة، وستتقارب. لذا، في ظل هذه الظروف المواتية، يمكن اللجوء إلى طريقة نيوتن إذا رغبنا في تقليل الخطأ إلى أدنى حد ممكن وتحقيق تقارب سريع جداً.

مثال: نمو نبات البردي

في الفصل  السابع من كتاب "الفصول التسعة" ، يمكن ترجمة مسألة إيجاد الجذر إلى اللغة الحديثة على النحو التالي:

مشكلة الفائض والعجز رقم 11:

  • نما نبات البردي بمقدار 3 وحدات في  يومه الأول. وفي نهاية كل يوم، لوحظ أن النبات قد نما بمقدار نصف نمو اليوم السابق.
  • نما نبات السعد بمقدار وحدة واحدة في يومه الأول. وفي نهاية كل يوم، يكون النبات قد نما بمقدار ضعف  نمو اليوم السابق.
  • أوجد الوقت [بالأيام الجزئية] الذي يصبح فيه نبات الهراوة بنفس طول نبات البردي.

إجابة:(2+613){\displaystyle \left(2+{\frac {6}{13}}\right)}أيام؛ الارتفاع هو(4+810+6130){\displaystyle \left(4+{\frac {8}{10}}+{\frac {6}{130}}\right)}وحدات.

توضيح:

  • لنفترض أنه اليوم  الثاني. يبلغ طول نبات الهراوة أقصر من نبات البردي بمقدار 1.5  وحدة.
  • لنفترض أنه اليوم  الثالث. نبات الهراوة أطول من نبات البردي بمقدار 1.75  وحدة. ∎
رسم بياني للدالة F ، وجذرها الدقيق (النقطة K )، والجذر التقريبي

لفهم هذا، سنقوم بنمذجة ارتفاعات النباتات في اليوم n ( n = 1، 2، 3...) وفقًا لسلسلة هندسية .

ب(ن)=أنا=1ن312أنا-1{\displaystyle B(n)=\sum _{i=1}^{n}3\cdot {\frac {1}{2^{i-1}}}\quad }نبات البردي
ج(ن)=أنا=1ن12أنا-1{\displaystyle C(n)=\sum _{i=1}^{n}1\cdot 2^{i-1}\quad }اندفاع النادي

من أجل تدوين أفضل، دع ك=أنا-1 .{\displaystyle \ k=i-1~.}أعد كتابة سلسلة ارتفاع النبات ب(ن)، ج(ن) {\displaystyle \ B(n),\ C(n)\ }بدلالة k واستدعاء صيغة المتسلسلة الهندسية.

 ب(ن)=ك=0ن-1312ك=3(1-(12)ن-1+11-12)=6(1-12ن){\displaystyle \ B(n)=\sum _{k=0}^{n-1}3\cdot {\frac {1}{2^{k}}}=3\left({\frac {1-({\tfrac {1}{2}})^{n-1+1}}{1-{\tfrac {1}{2}}}}\right)=6\left(1-{\frac {1}{2^{n}}}\right)}
 ج(ن)=ك=0ن-12ك=  1-2ن 1-2 =2ن-1 {\displaystyle \ C(n)=\sum _{k=0}^{n-1}2^{k}={\frac {~~1-2^{n}}{\ 1-2\ }}=2^{n}-1\ }

الآن، استخدم القاعدة الكاذبة لإيجاد جذر (ج(ن)-ب(ن)) {\displaystyle \ (C(n)-B(n))\ }

 F(ن):=ج(ن)-ب(ن)=62ن+2ن-7 {\displaystyle \ F(n):=C(n)-B(n)={\frac {6}{2^{n}}}+2^{n}-7\ }
  • تعيين x1=2 {\displaystyle \ x_{1}=2\ }واحسب F(x1)=F(2) {\displaystyle \ F(x_{1})=F(2)\ }وهو ما يساوي -1.5 (العجز).
  • تعيين x2=3 {\displaystyle \ x_{2}=3\ }واحسب F(x2)=F(3) {\displaystyle \ F(x_{2})=F(3)\ }وهو ما يساوي 1.75 (الزيادة).

الجذر المقدر (التكرار الأول):

 x^ =  x1F(x2)-x2F(x1) F(x2)-F(x1) =  2×1.75+3×1.5 1.75+1.5 = 3213  2.4615 {\displaystyle \ {\hat {x}}~=~{\frac {~x_{1}F(x_{2})-x_{2}F(x_{1})~}{F(x_{2})-F(x_{1})}}~=~{\frac {~2\times 1.75+3\times 1.5~}{1.75+1.5}}~=~{\frac {32}{13}}~\approx ~2.4615\ }

لإيجاد الجذر الدقيق، دعy=2ن{\displaystyle y=2^{n}}لذلك نسعى إلى حل6/y+y-7 = 0{\displaystyle 6/y+y-7~=~0}اضرب فيy{\displaystyle y}للحصول على المعادلة التربيعيةy2-7y+6=0{\displaystyle y^{2}-7y+6=0}والتي لها جذورy=0{\displaystyle y=0}(وهو أمر زائف هنا) وy=6{\displaystyle y=6}(الجذر الذي نريده). وبالتالين = سجل2(6)  2.5849625{\displaystyle n~=~\log _{2}(6)~\approx ~2.5849625}ويبلغ هامش الخطأ في تقديرنا 4.78%.

مثال على التعليمات البرمجية

هذا البرنامج النموذجي، المكتوب بلغة البرمجة C ، هو مثال على خوارزمية إلينوي. لإيجاد العدد الموجب x حيث cos( x ) = ، يتم تحويل المعادلة إلى صيغة إيجاد الجذر f ( x ) = cos( x ) - = 0 .

#include <stdio.h> #include <math.h>double f ( double x ) { return cos ( x ) - x * x * x ; } /* a,b: نهايتا فترة البحث  e: نصف الحد الأعلى للخطأ النسبي  m: الحد الأقصى لعدد التكرارات */ double falsi_method ( double ( * f )( double ), double a , double b , double e , int m ) { double c , fc ; int n , side = 0 ; /* القيم الابتدائية عند نهايتي الفترة */ double fa = f ( a ); double fb = f ( b );for ( n = 0 ; n < m ; n ++ ) { c = ( fa * b - fb * a ) / ( fa - fb ); if ( fabs ( b - a ) < e * fabs ( b + a )) break ; fc = f ( c );إذا كان ( fc * fb > 0 ) { /* fc و fb لهما نفس الإشارة، انسخ c إلى b */ b = c ; fb = fc ; إذا كان ( side == -1 ) fa /= 2 ; side = -1 ; } وإلا إذا كان ( fa * fc > 0 ) { /* fc و fa لهما نفس الإشارة، انسخ c إلى a */ a = c ; fa = fc ; إذا كان ( side == +1 ) fb /= 2 ; side = +1 ; } وإلا { /* fc * f_ صغير جدًا (يبدو كصفر) * / break ; } } return c ; }int main ( void ) { printf ( "%0.15f \n " , falsi_method ( & f , 0 , 1 , 5E-15 , 100 )); return 0 ; }

بعد تشغيل هذا الكود، تكون الإجابة النهائية تقريبًا 0.865474033101614.

انظر أيضاً

  • طريقة ITP ، وهي طريقة معدلة تضمن التقارب الأدنى والأقصى والتقارب فوق الخطي.
  • طريقة ريدرز ، وهي طريقة أخرى لإيجاد الجذور تعتمد على طريقة الوضع الخاطئ
  • طريقة برنت

مراجع

  1. 1 2 كاتز، فيكتور ج. (1998)، تاريخ الرياضيات (  الطبعة الثانية)، أديسون ويسلي لونجمان، ص 15 ، ISBN  978-0-321-01618-8
  2. 1 2 3 4 سميث، دي إي (1958) [1925]، تاريخ الرياضيات ، المجلد الثاني، دوفر، الصفحات 437-441 ، ISBN   978-0-486-20430-7{{citation}}عدم توافق رقم ISBN / التاريخ ( مساعدة )
  3. شابير، جان لوك، محرر. (2012) [1999]. "3. طرق الموقف الخاطئ" . تاريخ الخوارزميات: من الحصاة إلى الشريحة الإلكترونية . سبرينغر. ص 86-91 . ISBN  978-3-642-18192-4.
  4. نيدهام، جوزيف (1959). الرياضيات وعلوم السماء والأرض . العلم والحضارة في الصين. المجلد 3. مطبعة جامعة كامبريدج. ص 147 وما بعدها. ISBN   978-0-521-05801-8.{{cite book}}عدم توافق رقم ISBN / التاريخ ( مساعدة )
  5. "تسعة فصول" . www-groups.dcs.st-and.ac.uk . تم الاطلاع عليه بتاريخ 16 فبراير 2019 .
  6. شين، كانغشين؛ كروسلي، جون ن.؛ لون، أنتوني واه-تشيونغ (1999). الفصول التسعة في الفن الرياضي: دليل وتعليق . مطبعة جامعة أكسفورد. ص 358. ISBN  978-7-03-006101-0.
  7. 1 2 شوارتز، آر كيه (2004). قضايا في أصل وتطور حساب الخطائين (الحساب بالوضع الخاطئ المزدوج) . الاجتماع الثامن لشمال إفريقيا حول تاريخ الرياضيات العربية. رادس، تونس.متاح على الإنترنت على الرابط التالي: http://facstaff.uindy.edu/~oaks/Biblio/COMHISMA8paper.doc مؤرشف بتاريخ 15 سبتمبر 2011 في أرشيف الإنترنت ( Wayback Machine ) ونسخة مؤرشفة (PDF) . مؤرشف من النسخة الأصلية (PDF) بتاريخ 16 مايو 2014. تم الاطلاع عليه بتاريخ 8 يونيو 2012 .{{cite web}}: CS1 maint: archived copy as title ( link )
  8. ^ تراتاتو العامة ، المجلد. ط، البندقية، 1556، ص. فول. 238، v، Regola Helcataym (vocabulo Arabo) che in nostra lingua vuol dire delle false Positioni  
  9. ^ كونتي، إس دي؛ بور، كارل دي (1965). التحليل العددي الأولي: نهج خوارزمي ( الطبعة الثانية). ماكجرو هيل. ص. 40. أو سي إل سي 1088854304 .   
  10. 1 2 دالكويست, جيرموند ; بيورك، آكي (2003) [1974]. الطرق العددية . دوفر. ص 231 – 232. ISBN  978-0486428079.
  11. 1 2 فورد، جيه إيه (1995). "خوارزميات محسّنة من نوع إلينوي للحل العددي للمعادلات غير الخطية" . معاملات ACM في البرمجيات الرياضية . 30 : 64-85 . تم الاسترجاع في 1 يوليو 2025 .
  12. داول، م.؛ جاريت، ب. (1971). "طريقة معدلة لقاعدة الخطأ لحساب جذر المعادلة". BIT . 11 (2): 168–174 . doi : 10.1007/BF01934364 . S2CID 50473598 . 
  13. كينغ، ريتشارد ف. (أكتوبر 1983). "خوارزمية أندرسون-بيورك للمتتاليات الخطية". رياضيات الحساب . 41 (164): 591-596 . doi : 10.2307/2007695 . JSTOR 2007695 . 
  14. غالدينو، سيرجيو (2011). "مجموعة من طرق إيجاد الجذور باستخدام قواعد خاطئة" . وقائع المؤتمر العالمي للهندسة والتكنولوجيا لعام 2011. 1. تم الاطلاع عليه بتاريخ 9 سبتمبر 2016 .
  15. أوليفيرا، IFD؛ تاكاهاشي، RHC (2020-12-06). "تحسين طريقة التنصيف مع الحفاظ على متوسط ​​الأداء الأمثل لـ Minmax" . معاملات ACM في البرمجيات الرياضية . 47 (1): 5:1–5:24. doi : 10.1145/3423597 . ISSN 0098-3500 . S2CID 230586635 .  

للمزيد من القراءة