روش نیوتن

در آنالیز عددی روش نیوتن ، که همچنین به عنوان روش نیوتن-رافسون  نیز شناخته میشود الگوریتم ریشه یابی است که تقریب های خوبی در نزدیکی ریشه یک تابع (صفرهای یک تابع) میزند.در  پایه ای ترین حالت، الگوریتم نیوتن برای یک تابعی چون   با متغیر  و با مشتق  به همراه حدس اولیه  بکار میرود. اگر تابع حدس کافی و دقیقی را برآورد سازد و همچنین حدس اولیه نزدیک به ریشه تابع مفروض باشد (که با همگرایی تقریب ها این موضوع روشن می شود) آنگاه تقریب بهتری نسبت به  به حساب می آید.چرا که با احتساب همگرایی جواب ها، هر تقریب نسبت به تقریب قبل از خودش از دقت بالاتری برخوردار بوده و به ریشه تابع نزدیک تر است.به لحاظ هندسینقطه ای است که محور  و خط مماس تابع   در نقطهٔ  یکدیگر را قطع میکنند. شکل عمومی الگوریتم نیوتن به شرح زیر میباشد:

که در اصل از رابطه  بدست امده است. میدانیم که در نقطهٔ برخورد تابع با محور مقدار تابع صفر خواهد بود  لذا  که در آخر با تقسیم بر میتوان رابطه را به فرم رو به رو بازنویسی کرد:

همانطور که مشهود است روش نیوتن-رافسون از سری تیلور ناقص تابع مفروض به عنوان یک تقریب خطی حول نقطهٔ حدس اولیه بهره میبرد و از این جهت تقریب را ناقص میگویند که نیازی به نوشتن سری تابع تا مراتب بالاتر نبوده و به همان دو جمله ابتدایی بسنده میکند که این موضوع نیز دلیلی بر تقریب خطی بودن روش نیوتن میباشد. همچنین چون این روش معادلهٔ یک تابع را تا معادلهٔ یک تابع درجه یک تقیل میدهد، لذا صرف نظر از اینکه تابع چند ریشه دارد،  در نهایت الگوریتم  تنها یک جواب بدست می آورد.

این روش همچنین میتواند در توابع مختلط و دستگاه معادلات بکار رود.

شرح 
ساز و کار روش نیوتن شروع مراحل تخمین ریشه با انتخاب یک حدس اولیه است که به اندازه لازم به مقدار واقعی ریشه نزدیک باشد ، سپس میتوان با استفاده از خط مماس تابع را تقریب زد و نهایتاً با توجه به خط مماس طول از مبدا (ریشه تقریب خطی) را معین ساخت. ریشه تقریب معمولاً تقریب بهتری از ریشه واقعی تابع نسبت به حدس اولیه است و الگوریتم به همین منوال میتواند تکرار شود.

چگونگی نامگذاری الگوریتم 
اسم " Newton's method " از شرح آیزاک نیوتن در حالت خاصی از روش مذکور در "آنالیز بوسیله سری های بیکران" (نوشته شده در 1669،منتشر شده در سال 1711 توسط ویلیام جونز  که در اصل کاری ریاضیاتی از نیوتن بود) و در De metodis fluxionum et serierum infinitarum (نوشته شده در 1671،که در سال 1736 توسط جان کولسون ترجمه و منتشر شد)مشتق شده است.

عدم موفقیت الگوریتم 
روش نیوتن تنها زمانی یافتن ریشه تقریبی را تضمین میکند که شرایطی برقرار باشد.
نقاطی با شروع بد
در برخی مواقع شرایطی که برای همگرایی لازم است برقرار است،اما نقطهٔ که به عنوان حدس اولیه انتخاب شده دربازه ای که روش نیوتن همگرایی دارد قرار نمی‌گیرد.این شرایط زمانی میتواند اتفاق بیفتد که برای مثال همچنان که  به  یا  میل میکند ریشه  تابع  به طور متناوب و به اندازه کافی و مطلوبی به صفر نزدیک شود. در اینگونه موارد روش بهتری مانند روش دوبخشی (تنصیف) باید بکار گرفته شود تا تخمین بهتری بدست آورد.
زمانی که نقطهٔ تکرار ثابت است
به تابع : دقت کنید:

این تابع در نقطهٔ  ماکسیمم دارد و ریشه های آن   و  میباشد. اگر تکرار را از  شروع کنیم (نقطه ای که مشتق آن صفر و شیب خط تابع افقی است) تعریف نشده خواهد بود:

حتی اگر مشتق تابع در آن نقطه نزدیک به صفر باشد آنگاه تکرار مرحله، تقریب به مراتب بدتری را به همراه خواهد داشت.
حدس اولیه به یک دایره منتهی میشود
در برخی توابع، بعضی از نقاطی که به عنوان حدس اولیه انتخاب میشوند ممکن است که به یک دایره ختم شوند که این میتواند مانع از همگرایی شود.
برای مثال تابع  را در نظر بگیرید و  را به عنوان حدس اولیه اختیار کنید. نخستین مرحله جواب را  بدست می آوردو دومین مرحله دوباره به  میرسد،بنابراین جواب ها مدام بین دو مقدار مذکور تناوب میکند.
در حقیقت این دو نقطه ثابت هستند.در حقیقت رفتار در این حالت میتواند بسیار پیچیده باشد. جواب واقعی معادله  می باشد.

مثالها 
جذر یک عدد
با توجه به یافتن جذر  یک عدد، روش نیوتن یکی از چندین راهی است که میتوان از آن بهره گرفت.
برای مثال اگر هدف یافتن جذر 612 باشد ، این موضوع را میتوان معادل جواب  دانست. سپس تابعی را که قصد داریم تا روش نیوتن را برای آن بکار ببندیم میتوان به فرم  در نظر گرفت که مشتق آن نیز  خواهد بود.اگر حدس اولیه را  در نظر بگیریم آنگاه الگوریتم نیوتن به صورت زیر محاسبه را انجام خواهد داد:

جاییکه که زیر رقم صحیح خط کشیده شده است ،تنها با چند تکرار میتوان به یک جواب دقیق با ارقام اعشاری دست یافت. در صورت هگرایی هرچه روش را تکرار کنیم جواب ها دقیق تر و ارقام بعد اعشار نیز به مراتب از دقت بیش تری برخوردار خواهند بود. همانطور که مشخص است برخی از ارقام بعد ممیز در هر مرحله تکرار میشوند که این نشان از دقت و قطعیت آنها دارد. اما به صورت کلی یافتن چنین جواب هایی با بهره گیری از روش های عددیابی چون روش نیوتن،
روش تکرار ساده ( Iteration method)، دو بخشی (Bisection method)، وتری (Secant method)، نابجایی (False position)، روش مولر (Muller's method) و ... نمیتوان به کامل ترین جواب دست یافت، چرا که ارقام بعد ممیز همواره ادامه دارند و خاتمه نمی‌یابند از این رو تنها دقیق ترین و کامل ترین جواب تنها در صورتی بدست می آید که تکرار را تا بینهایت ادامه دهیم،با فرض اینکه ریشهٔ تابع  باشد، به زبان ریاضیاتی این موضوع را میتوان در قالب زیر بیان کرد:

یافتن جواب 
میتوان معادله را به فرم  بازنویسی کرد و همچنین  بدین صورت میباشد. از آنجایی که که برای هر   همواره  میباشد لذا برای هر  نیز  خواهد بود. میدانیم که پاسخ معادله بین  و  قرار دارد لذا حدس اولیه را با  شروع میکنیم ( توجه داشته باشید که اگر حدس اولیه را با  شروع میکردیم آنگاه به یک نتیجه تعریف نشده می‌رسیدیم،این موضوع اهمیت انتخاب حدس اولیه را به خوبی روشن می‌سازد که به اندازه کافی باید به جواب واقعی نزدیک باشد).

زیر ارقام صحیح و قطعی در مثالهای بالا خط کشیده شده است. به خصوص   که تا 12 رقم بعد ممیز صحیح میباشد. میبینیم که تعداد ارقام صحیح از 2 رقم برای   تا 10 رقم برای   در حال افزایش است که نشان دهنده همگرایی است.

مقایسه سرعت همگرایی در روش های مختلف 

در این قسمت اهمیت انتخاب روش مناسب به خوبی مشخص خواهد شد، چرا که سرعت همگرایی الگوریتم های مذکور با یکدیگر متفاوت است.

برای مثال تابع را در نظر بگیرید:

 

در این مثال سرعت همگرایی دو روش وتری (خط قاطع) و روش نیوتن با یکدیگر تقریباً برابر است، همچنین مشاهده میشود که کند ترین روش در این مثال الگوریتم نابجایی میباشد، اما به‌طور کل نمیتوان این مثال خاص را به تمامی مثال ها تعمیم داد چرا که همواره سرعت همگرایی متغیر بوده و به شرایط مختلفی بستگی دارد، لذا انتخاب بهترین روش میتواند سهم به سزایی در سرعت عملیات تقریب زنی داشته باشد.همچنین مسئله ی دیگر که مطرح است این موضوع میباشد که هر کدام از روش ها مزایا و معایب خود را دارا هستند، برای مثال از معایب روش نیوتن عدم جواب در صورتی ست که مشتق تابع در نقطه ی مورد نظر صفر و یا نزدیک به صفر باشد و یا ایجاد سختی در هنگامی ست که مشتق گیری تابع عملیاتی دشوار باشد و از مزایای آن همگرایی سریع در صورت انتخاب حدس اولیه مناسب است ،اما باید توجه داشت همگرایی روش نیوتن-رافسون برخلاف برخی روش ها مانند روش دوبخشی تضمینی نیست.

نقش انتخاب بازه در تعیین ریشه 
یکی از موارد کلیدی در اینکه روش مورد استفاده به کدام یک از ریشه های تابع (در صورت وجود بیش از یک ریشه) همگرا شود انتخاب بازه میباشد.
برای مثال تابع  را در نظر بگیرید و نقطه ی  را به عنوان حدس اولیه اختیار کنید:

در معادله ی فوق بازه ی انتخاب شده برای هر چهار روش یکسان بوده اما ریشه ی  تقریبی حاصل از روش نیوتن متفاوت از سایرین است.

منابع 

 

الگوریتم‌های پیدا کردن ریشه
روش‌های بهینه‌سازی