روش رانگ-کوتا

Runge-Kutta Method of Order 4
روش رانگ کوتا در متلب
معرفی

روش رانگ-کوتای مرتبه 4 یکی از روش های انتگرال گیری عددی از معادلات دیفرانسیل معمولی (ODE) است که برای حل معادلات به شکل

dy/dx = f(x,y)

مورد استفاده قرار می گیرد. در این روش داریم:

y_{i+1}=y_i+{1/6}delim{[}{k_1 + 2 k_2 + 2 k_3 + k_4}{]}

که در آن

k_1 =h f(x_i,y_i)
k_2 =h f(x_i + {1/2}h,y_i + {1/2}k_1)
k_3 =h f(x_i + {1/2}h,y_i + {1/2}k_2)
k_4 =h f(x_i + h,y_i + k_3)

 هنگامی که تابع f(x,y) تنها تابعی از x باشد روش رانگ-کوتای مرتبه چهار به روش انتگرال گیری سیمپسون تبدیل می شود.

    برنامه نویسی با MATLAB

    در اینجا با استفاده از نرم افزار متلب (MATLAB) برنامه ای جهت حل معادلات دیفرانسیل معمولی مرتبه اول به روش رانگ کوتای مرتبه 4 ارائه گردیده است. دو کد برای این منظور در یک فایل فشرده ارائه گردیده است:

    1. کد صریح (explicit) که در خروجی روند حل کامل مسئله را نمایش می دهد

    2. کد غیرصریح (معمولی) که تنها پاسخ نهایی را نمایش می دهد

    لازم به ذکر است که برنامه ی ارائه شده قادر به حل تمامی مثال های قابل حل با روش رانگ کوتا بوده و به صورت کاملا عمومی (general) کدنویسی شده است.

    ورودی ها و خروجی ها

    ورودی:

    1. تابع (f(x,y ورودی
    2. نقطه ی ابتدایی بازه ی حل
    3. نقطه ی انتهایی بازه ی حل
    4. مقدار اولیه ی تابع y در نقطه ی ابتدایی بازه
    5. تعداد گام های حل n

    خروجی:

    1. محاسبه و نمایش طول گام حل h
    2. نمایش روند حل به صورت کامل در هر گام
    3. نمایش جدول مقادیر تابع حاصل از حل به روش اولر
    4. رسم نمودار حاصل از حل

    تصاویر اجرای برنامه

    مشاهده ی ورودی و خروجی برنامه در یک مثال نمونه

    مشاهده ی ورودی و خروجی برنامه ی غیر صریح در یک مثال

    نمودار خروجی رسم شده برای مثال نمونه

    این برنامه با نسخه های متلب/MATLAB سال های 2010-2013 تست شده است. در صورت استفاده از نسخه ی سال های دیگر لایبریکا تضمین کننده ی اجرای صحیح برنامه نمی باشد.