طريقة القاطع

أول دورتين من طريقة القاطع. يُمثل المنحنى الأحمر الدالة f ، بينما تُمثل الخطوط الزرقاء القواطع. في هذه الحالة تحديدًا، لن تتقارب طريقة القاطع إلى الجذر الظاهر.

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

الطريقة

طريقة القاطع هي طريقة عددية تكرارية لإيجاد جذر دالة f . بمعلومية قيمتين ابتدائيتين x₀ و x₁ ، تعمل الطريقة وفقًا للعلاقة التكرارية .

xن=xن-1-و(xن-1)xن-1-xن-2و(xن-1)-و(xن-2)=xن-2و(xن-1)-xن-1و(xن-2)و(xن-1)-و(xن-2).{\displaystyle x_{n}=x_{n-1}-f(x_{n-1}){\frac {x_{n-1}-x_{n-2}}{f(x_{n-1})-f(x_{n-2})}}={\frac {x_{n-2}f(x_{n-1})-x_{n-1}f(x_{n-2})}{f(x_{n-1})-f(x_{n-2})}}.}

هذه علاقة تكرارية غير خطية من الدرجة الثانية، وهي معرفة جيدًا بمعلومية الدالة f والقيمتين الابتدائيتين x₀ و x₁ . من الناحية المثالية، ينبغي اختيار القيم الابتدائية قريبة من الصفر المطلوب.

اشتقاق الطريقة

انطلاقًا من القيمتين الابتدائيتين x₀ و x₁، نرسم خطًا مستقيمًا يمر بالنقطتين (x₀, f( x₀ ) ) و ( x₁ , f ( x₁ ) ) ، كما هو موضح في الصورة أعلاه . معادلة هذا الخط ، بصيغة نقطة-نقطة ، [ 2 ] هي

y=و(x1)-و(x0)x1-x0(x-x1)+و(x1).{\displaystyle y={\frac {f(x_{1})-f(x_{0})}{x_{1}-x_{0}}}(x-x_{1})+f(x_{1}).}

جذر هذه الدالة الخطية ، أي قيمة x التي تحقق y = 0، هو

x=x1-و(x1)x1-x0و(x1)-و(x0).{\displaystyle x=x_{1}-f(x_{1}){\frac {x_{1}-x_{0}}{f(x_{1})-f(x_{0})}}.}

ثم نستخدم هذه القيمة الجديدة لـ x كـ x2 ونكرر العملية، باستخدام x1 و x2 بدلاً من x0 و x1 . نستمر في هذه العملية، ونحل لإيجاد x3 و x4 ، وهكذا ، حتى نصل إلى مستوى عالٍ من الدقة (فرق صغير بما فيه الكفاية بين xn و xn - 1 ) .

x2=x1-و(x1)x1-x0و(x1)-و(x0)،x3=x2-و(x2)x2-x1و(x2)-و(x1)،xن=xن-1-و(xن-1)xن-1-xن-2و(xن-1)-و(xن-2).\begin{aligned}x_{2}=x_{1}-f(x_{1}){\frac {x_{1}-x_{0}}{f(x_{1})-f(x_{0})}},\\[6pt]x_{3}=x_{2}-f(x_{2}){\frac {x_{2}-x_{1}}{f(x_{2})-f(x_{1})}},\\[6pt]&\,\,\,\vdots \\[6pt]x_{n}=x_{n-1}-f(x_{n-1}){\frac {x_{n-1}-x_{n-2}}{f(x_{n-1})-f(x_{n-2})}}.\end{aligned}}}

التقارب

التكرارxن{\displaystyle x_{n}}تتقارب نتائج طريقة القاطع إلى جذرو{\displaystyle f}إذا كانت القيم الأوليةx0{\displaystyle x_{0}}وx1{\displaystyle x_{1}}قريبة بما يكفي من الجذر وو{\displaystyle f}حسن السلوك. عندماو{\displaystyle f}الدالة قابلة للتفاضل مرتين بشكل مستمر، والجذر المعني هو جذر بسيط، أي أن له رتبة 1، ورتبة التقارب هي النسبة الذهبية.φ=(1+5)/21.618.{\displaystyle \varphi =(1+{\sqrt {5}})/2\approx 1.618.}[ 3 ] هذا التقارب هو فوق الخطي ولكنه دون التربيعي.

إذا لم تكن القيم الأولية قريبة بما يكفي من الجذر أوو{\displaystyle f}إذا لم يكن سلوك الدالة منتظمًا، فلا يوجد ضمان لتقارب طريقة القاطع على الإطلاق. لا يوجد تعريف عام لعبارة "قريب بما فيه الكفاية"، ولكن معيار التقارب يتعلق بمدى "تذبذب" الدالة على الفترة بين القيم الابتدائية. على سبيل المثال، إذاو{\displaystyle f}تكون قابلة للتفاضل على تلك الفترة، وهناك نقطة حيثو=0{\displaystyle f'=0}في هذه الفترة، قد لا تتقارب الخوارزمية.

مقارنة مع طرق أخرى لإيجاد الجذور

لا تتطلب طريقة القاطع، ولا تضمن، بقاء الجذر محصورًا بين تكرارات متسلسلة، كما هو الحال في طريقة التنصيف ، ولذلك فهي لا تتقارب دائمًا. تستخدم طريقة الوضع الخاطئ (أو القاعدة الخاطئة ) نفس صيغة طريقة القاطع، ولكنها لا تطبق الصيغة على الجذر.xن-1{\displaystyle x_{n-1}}وxن-2{\displaystyle x_{n-2}}، مثل طريقة القاطع، ولكن علىxن-1{\displaystyle x_{n-1}}وفي التكرار الأخيرxك{\displaystyle x_{k}}بحيثو(xك){\displaystyle f(x_{k})}وو(xن-1){\displaystyle f(x_{n-1})}لها إشارة مختلفة. هذا يعني أن طريقة الوضع الخاطئ تتقارب دائمًا؛ ولكن فقط بترتيب تقارب خطي. يمكن تحقيق التقارب بترتيب تقارب فوق الخطي، كما هو الحال في طريقة القاطع، من خلال تحسينات على طريقة الوضع الخاطئ (انظر Regula falsi § تحسينات في Regula falsi ) مثل طريقة ITP أو طريقة إلينوي .

يمكن اشتقاق صيغة التكرار لطريقة القاطع من صيغة طريقة نيوتن

xن=xن-1-و(xن-1)و(xن-1){\displaystyle x_{n}=x_{n-1}-{\frac {f(x_{n-1})}{f'(x_{n-1})}}}

باستخدام تقريب الفروق المحدودة ، لقيمة صغيرةϵ=xن-1-xن-2{\displaystyle \epsilon =x_{n-1}-x_{n-2}}:

و(xن-1)=ليمϵ0و(xن-1)-و(xن-1-ϵ)ϵو(xن-1)-و(xن-2)xن-1-xن-2{\displaystyle f'(x_{n-1})=\lim _{\epsilon \rightarrow 0}{\frac {f(x_{n-1})-f(x_{n-1}-\epsilon )}{\epsilon }}\approx {\frac {f(x_{n-1})-f(x_{n-2})}{x_{n-1}-x_{n-2}}}}

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

إذا قارنا طريقة نيوتن بطريقة القاطع، نلاحظ أن طريقة نيوتن تتقارب أسرع (من الرتبة 2 مقابل الرتبة φ ≈ 1.6 للنسبة الذهبية ). [ 3 ] مع ذلك، تتطلب طريقة نيوتن حساب كليهما  و{\displaystyle f}ومشتقاتهو{\displaystyle f'}في كل خطوة، بينما تتطلب طريقة القاطع فقط تقييمو{\displaystyle f}لذلك، قد تكون طريقة القاطع أسرع في بعض الأحيان من الناحية العملية. على سبيل المثال، إذا افترضنا أن تقييمو{\displaystyle f}إذا استغرق حساب المشتقة نفس الوقت اللازم لحساب المشتقة ، وتجاهلنا جميع التكاليف الأخرى، فيمكننا تنفيذ خطوتين من طريقة القاطع (تقليل لوغاريتم الخطأ بمعامل φ² ≈ 2.6) بنفس تكلفة خطوة واحدة من طريقة نيوتن (تقليل لوغاريتم الخطأ بمعامل 2)، لذا فإن طريقة القاطع أسرع. في الأبعاد الأعلى، قد يصبح حساب المجموعة الكاملة من المشتقات الجزئية المطلوبة لطريقة نيوتن، أي مصفوفة جاكوبي ، أكثر تكلفة بكثير من حساب الدالة نفسها. مع ذلك، إذا أخذنا في الاعتبار المعالجة المتوازية لحساب المشتقة أو المشتقات، فقد تكون طريقة نيوتن أسرع من حيث وقت المعالجة، على الرغم من أنها لا تزال تتطلب عمليات حسابية أكثر إجمالاً.   

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

رياضياً، فإن الشكلين التاليين لطريقة القاطع متكافئان: xن+1=xن-و(xن)[xن-xن-1و(xن)-و(xن-1)]{\displaystyle x_{n+1}=x_{n}-f(x_{n})\left[{\frac {x_{n}-x_{n-1}}{f(x_{n})-f(x_{n-1})}}\right]} و xن+1=xن-1و(xن)-xنو(xن-1)و(xن)-و(xن-1).{\displaystyle x_{n+1}={\frac {x_{n-1}f(x_{n})-x_{n}f(x_{n-1})}{f(x_{n})-f(x_{n-1})}}.}

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

  • وهو على شكلxن+1=xن-و(xن)Δ{\displaystyle x_{n+1}=x_{n}-f(x_{n})\Delta }، لعدد ماΔ{\displaystyle \Delta }حتى لوΔ{\displaystyle \Delta }إذا كانت القيمة غير دقيقة، فسيكون التغيير في تقدير الجذر صغيرًا عند التقارب لأنو(xن){\displaystyle f(x_{n})}سيكون صغيرًا أيضًا.
  • أما الشكل الثاني فهو عرضة للإلغاء الكارثي ، لأنو(xن)و(xن-1){\displaystyle f(x_{n})\to f(x_{n-1})}مع اقتراب التقارب، يمكن أن يكون خطأ الإلغاء في المقام كبيرًا.

تعميم

تُعد طريقة برودن تعميمًا لطريقة القاطع لأكثر من بُعد واحد.

يُظهر الرسم البياني التالي الدالة f باللون الأحمر وخط القاطع الأخير باللون الأزرق الغامق. في الرسم البياني، يبدو أن نقطة تقاطع خط القاطع مع المحور السيني تُقارب جذر الدالة f بشكل جيد .

مثال

لنفترض أننا نبحث عن حل لكثير الحدود

و(x)=2x3-x-3.{\displaystyle f(x)=2x^{3}-x-3\,.}

تطبيق طريقة القاطع مع التخمينات الأوليةx0=1{\displaystyle x_{0}=1}وx1=2{\displaystyle x_{1}=2}يعطي التقريبات التالية:

التكرارxن{\displaystyle x_{n}}و(xن){\displaystyle f(x_{n})}
21.153846153846153{\displaystyle 1.153846153846153}-1.0814747382794727{\displaystyle -1.0814747382794727}
31.334985133795837{\displaystyle 1.334985133795837}0.4233966484501184{\displaystyle 0.4233966484501184}
41.284021551332461{\displaystyle 1.284021551332461}-5.00597×10-2{\displaystyle -5.00597\times 10^{-2}}
51.289410061128775{\displaystyle 1.289410061128775}-1.91967×10-3{\displaystyle -1.91967\times 10^{-3}}
61.289624937539790{\displaystyle 1.289624937539790}9.30251×10-6{\displaystyle 9.30251\times 10^{-6}}
71.289623901294108{\displaystyle 1.289623901294108}-1.77635×10-15{\displaystyle -1.77635\times 10^{-15}}

الخوارزمية

يصف الكود الزائف التالي طريقة القاطع لتقريب جذر المعادلةو(x)=0{\displaystyle f(x)=0}. يُعيد تقديرًاx2{\displaystyle x_{2}}ذلك يرضي|x2-x1|<ε{\displaystyle |x_{2}-x_{1}|<\varepsilon }، أينε{\displaystyle \varepsilon }هو هامش سماحية يحدده المستخدم.

المدخلات: الدالة f ، التقريبات الأولية x0 و x1 ، التسامح ε ، الحد الأقصى للتكرارات Nmax. من أجل n = 0 إلى Nmax، نفّذ ما يلي: x2x1 - f ( x1 )·( x1 - x0 ) / ( f ( x1 ) - f ( x0 )). إذا كان | x2x1 | < ε ، فأرجع x2 . نهاية الحلقة. x0 ، x1x1 ، x2. نهاية الحلقة. أرجع x2.

ملحوظات

  1. باباكستانتينو، جوانا؛ تابيا، ريتشارد (2013). "أصل وتطور طريقة القاطع في بُعد واحد" . المجلة الرياضية الأمريكية الشهرية . 120 (6): 500-518 . doi : 10.4169/amer.math.monthly.120.06.500 . JSTOR 10.4169/amer.math.monthly.120.06.500 . S2CID 17645996 .  
  2. مارسدن، جيرولد (1985). حساب التفاضل والتكامل 1. سبرينغر-فيرلاغ نيويورك، ص 31. ISBN  978-1-4612-5024-1.
  3. 1 2 شانسون، جيفري ر. (3 أكتوبر 2024). "رتبة التقارب" . ليبرتيكستس ماثيماتيكس . تم الاسترجاع في 3 أكتوبر 2024 .

انظر أيضاً

مراجع