روش تکرار نقطه ثابت

در بحث محاسبه غلظتهای تعادلی در واکنشهای تعادلی، همچنین خوب است با بحث حل معادلات با استفاده از روش­های عددی آشنا باشیم. همانطور که تاکنون دیده­ایم، در اغلب موارد برای بدست آوردن غلظت­های تعادلی نیازمند حل معادلات و بدست آوردن مجهولات درنظر گرفته شده هستیم. برای حل معادلات بطور کلی دو دسته روش درنظر گرفته می­شود: روش­های تحلیلی و روش­های عددی یا معادلاً روش­های مستقیم (Direct Methods) و روش­های تکرار شونده (Iterative Methods). روش­های مستقیم جواب دقیق یک معادله را براساس یک الگوریتم معین و با طی تعدادی مراحل مشخص نتیجه می­دهند. بعنوان مثال روش دلتا یک روش مستقیم برای یافتن جواب­ها یا ریشه­های معادلات درجه دو است. متأسفانه برای حل بسیاری از معادلات روش­های مستقیم وجود ندارد و ما برای حل اغلب معادلات مجبوریم به روش­های تکرار شونده یا عددی متوسل شویم. روش های عددی متنوعی وجود دارند و بسیاری از نرم­افزارها در کامپیوترها یا ماشین حساب­ها از این روش­ها به یافتن جواب استفاده می­کنند. در اینجا ما یک روش عددی برای حل معادلات به نام تکرار نقطه ثابت (Fixed-point Iteration) را ارائه می­دهیم.

نقاط ثابت یک تابع نقاطی هستند که برای آنها خروجی تابع با ورودی برابر است. بعبارتی برای تابع ، نقاطی که در  صدق می­کنند، نقاط ثابت تابع  محسوب می­شوند. صفرهای تابع  نقاطی هستند که بازای آنها  است یا بعبارتی صفرهای تابع ، جواب­های معادله  هستند. در روش تکراری نقطه ثابت برای یافتن جواب­های معادله  (یا همان صفرهای تابع )، ابتدا  را به صورت  یا  درمی­آوریم (بعنوان مثال  است). سپس با استفاده از روش زیر نقاط ثابت تابع  را که همان صفرهای تابع  یا جواب­های معادله  است، می­یابیم. روش بدین صورت است که یک حدس اولیه برای  بعنوان جواب بصورت دلخواه درنظر می­گیریم و برای  درنظر گرفته شده  را محاسبه می­کنیم. سپس  محاسبه شده را بعنوان  جدید درنظر می­گیریم و مجدداً  را اینبار برای  جدید محاسبه می­کنیم و به همین صورت ادامه می­دهیم. (آیا می­توانید بگوئید چرا برای چنین روش­هایی واژه تکرار شونده را (Iterative Methods) بکار می­برند؟!!). اگر خوش شانس باشیم، خواهیم دید که  و  به همدیگر میل می­کنند یا بعبارتی مقدار آنها به همدیگر نزدیک­تر می شود. بنابراین یکی از نقاط ثابت تابع  یا درواقع جواب معادله  برای ما معلوم می­شود. این روش به صورت الگوریتم زیر می­تواند نشان داده شود:


 

 


 


دقت داشته باشید که در نمایش الگوریتم، از بیضی برای نمایش شروع و پایان، از مستطیل برای بیان دستورالعمل، از لوزی برای بیان شرط و از دایره برای بیان نتیجه استفاده می­کنند. برای این روش مثال زیر را درنظر بگیرید:

مثال: یکی از جواب­های معادله زیر را به روش تکرار نقطه ثابت بیابید.

 

حل:

ابتدا  را بصورت  داریم. بعنوان مثال داریم:

 

 

بنابراین  می­تواند برابر  درنظر گرفته شود. بعنوان حدس اولیه برای ،  را درنظر می­گیریم. سپس محاسباتی به صورت زیر انجام می­دهیم:

 

  

  

  

 

 

 

 

 

 

همانطور که دیده می­شود، هم مقدار  و هم مقدار  به  نزدیک می­شود و بنابراین  یکی از نقاط ثابت تابع  و یکی از جواب­های معادله  است (این مطلب را چگونه می­توانید چک کنید؟!)

دقت داشته باشید که چه برای این مسئله و چه برای مسائل دیگر می­توان ها و حدسهای اولیه گوناگونی درنظر گرفت. بعنوان مثال می­توانستیم یک  دیگر به صورت زیر بدست آوریم:

 

بنابراین  می­تواند برابر  درنظر گرفته شود. بعنوان حدس اولیه برای ،  را درنظر می­گیریم. سپس محاسباتی به صورت زیر انجام می­دهیم:

 

 

 

 

 

 

 

 

                                             

 

بنابراین با تغییر  و حدس اولیه برای  یکی دیگر از جواب­های معادله  را پیدا نمودیم.

*****

دقت داشته باشید که در مثال قبل خوش شانس بودیم که به جواب رسیدیم!!! در حالت کلی ممکن است در یک حلقه بدام بیفتیم، نتوانیم  را حساب کنیم و...، یا بعبارتی نتوانیم به جواب برسیم. در اینصورت بایستی یا حدس اولیه یا فرم تابع  را عوض نمائیم و امیدوار باشیم که با مقدار جدید برای حدس اولیه یا فرم جدید تابع  به جواب برسیم. همچنین باتوجه به مثال قبل می­بینیم که مقدار  به دو صورت زیر به جواب نهایی  میل می­کند (آیا می­توانید دو صورت زیر را توصیف نمائید؟!!):

 

 

 

 


حال کاربرد این روش را در یک مسئله تعادل نشان می­دهیم.

مثال: واکنش تعادلی زیر را درنظر بگیرید:

 

در شروع محلول محتوی  از  و مواد بی­اثر بر واکنش تعادلی بالا است. به محلول مقدار زیادی  می­افزائیم. بعد از برقراری تعادل، غلظت  و  را بدست آورید.

حل:

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

 

                                                :شروع

                                     :تعادلی

باتوجه به غیرکمّی بودن واکنش و ناچیز بودن ، از  در مقابل  صرف­نظر می­کنیم و می­نویسیم:

 

 

 

دقت کنید که  کوچک است، اما در مقابل  قابل صرف­نظر کردن نیست و اجازه نداریم  را بطور تقریبی برابر  فرض نمائیم. بنابراین بایستی به صورت زیر عمل کنیم:

 

 

معادله حاصله برای ، یک معادله از درجه  نسبت به  است که فرمول معینی برای حل آن وجود ندارد و بنابراین بروش تکرار نقطه ثابت به صورت زیر متوسل می­شویم:

 

باتوجه به عملیات بالا،  به فرم  است. با شروع از حدس اولیه  داریم:

 

 

 

 

 

                          

 

بنابراین  است و غلظت­های تعادلی  و  بصورت زیر هستند:

 

*****

خوشبختانه پیاده­سازی روش تکرار نقطه ثابت با ماشین حساب­های دارای دکمه  ساده است. کافیست حدس اولیه را در ماشین حساب وارد کنید (مثلاً حدس اولیه را نوشته و دکمه = یا  یا   یا... را بزنید). سپس  را بصورت  وارد نمائید (مثلاً اگر  است، در ماشین حساب  را وارد نمائید) و سپس آنقدر دکمه = یا  یا    یا... را بزنید که عدد نمایش داده شده در صفحه نمایش ماشین حساب تقریباً ثابت شود. با اینکار خود ماشین حساب  را حساب کرده و بعنوان  جدید درنظر می­گیرد و اگر خوش­شانس باشیم، با تکرار اینکار به جواب خواهیم رسید (و اگر خوش­شانس نبودیم، چه کار باید کنیم؟!!!).

 

 

 

||||||||||||||||||||||||||||||||||||||||||||||||||

**************************************************

سایت: www.youngchemist.com

مولف: محمد شاهی

نظرات، پیشنهادات و انتقادات: chemistry.shahi@gmail.com

**************************************************

||||||||||||||||||||||||||||||||||||||||||||||||||