Matematikte
Runge-Kutta-Fehlberg yöntemi (RKF yöntemi),
adi diferansiyel denklemlerin sayısal çözümü için kullanılan etkili bir yöntemdir.
Bu yöntem, klasik Runge-Kutta yaklaşımına dayanır ve özellikle
adım uzunluğunun otomatik olarak kontrol edilebilmesi avantajı ile öne çıkar.
Yöntem, Alman matematikçi
Erwin Fehlberg tarafından geliştirilmiştir.
Temel fikir, aynı adım içinde hem
4. dereceden hem de
5. dereceden iki farklı yaklaşık çözüm elde ederek
bunların farkından yerel hatayı tahmin etmektir. Böylece çözümün istenen doğrulukta elde edilmesi için
adım uzunluğu \(h\) büyütülüp küçültülebilir.
Aşağıdaki başlangıç değer problemini ele alalım:
\(y' = f(t,y)\), \(y(t_0)=y_0\)
Burada amaç, verilen başlangıç koşulunu sağlayan çözümün sayısal olarak yaklaşık değerlerini bulmaktır.
Runge-Kutta-Fehlberg yönteminde her adımda aşağıdaki altı ara değer hesaplanır:
\(k_1 = h f(t_n, y_n)\)
\(k_2 = h f\!\left(t_n+\displaystyle \frac{1}{4}h,\; y_n+\displaystyle \frac{1}{4}k_1\right)\)
\(k_3 = h f\!\left(t_n+\displaystyle \frac{3}{8}h,\; y_n+\displaystyle \frac{3}{32}k_1+\displaystyle \frac{9}{32}k_2\right)\)
\[k_4 = h f\!\left(t_n+\displaystyle \frac{12}{13}h,\; y_n+\displaystyle \frac{1932}{2197}k_1-\displaystyle \frac{7200}{2197}k_2+\displaystyle \frac{7296}{2197}k_3\right)\]
\(k_5 = h f\!\left(t_n+h,\; y_n+\displaystyle \frac{439}{216}k_1-8k_2+\displaystyle \frac{3680}{513}k_3-\displaystyle \frac{845}{4104}k_4\right)\)
\(k_6 = h f\!\left(t_n+\displaystyle \frac{1}{2}h,\; y_n-\displaystyle \frac{8}{27}k_1+2k_2-\displaystyle \frac{3544}{2565}k_3+\displaystyle \frac{1859}{4104}k_4-\displaystyle \frac{11}{40}k_5\right)\)
Bu ara değerler kullanılarak aynı adım için iki farklı yaklaşık çözüm hesaplanır.
Bunlardan biri 4. dereceden, diğeri ise 5. dereceden yaklaşık çözümdür.
4. dereceden yaklaşık çözüm:
\(y_{n+1}^{(4)} = y_n + \displaystyle \frac{25}{216}k_1 + \displaystyle \frac{1408}{2565}k_3 + \displaystyle \frac{2197}{4104}k_4 - \displaystyle \frac{1}{5}k_5\)
5. dereceden yaklaşık çözüm:
\(y_{n+1}^{(5)} = y_n + \displaystyle \frac{16}{135}k_1 + \displaystyle \frac{6656}{12825}k_3 + \displaystyle \frac{28561}{56430}k_4 - \displaystyle \frac{9}{50}k_5 + \displaystyle \frac{2}{55}k_6\)
Burada \(h\), adım uzunluğunu ifade eder.
İki çözüm arasındaki fark, yerel kesme hatasının yaklaşık bir göstergesi olarak kullanılır:
\(\text{Hata} \approx \left| y_{n+1}^{(5)} - y_{n+1}^{(4)} \right|\)
Eğer bu hata kabul edilebilir sınırlar içindeyse adım kabul edilir.
Hata büyükse adım küçültülerek hesap tekrar yapılır.
Hata çok küçükse sonraki adım için daha büyük bir \(h\) seçilerek hesaplama daha verimli hale getirilebilir.
Runge-Kutta-Fehlberg yönteminin başlıca avantajları şunlardır:
- Aynı anda iki farklı mertebeden çözüm üretir.
- Yerel hata tahmini yapmaya imkân verir.
- Adım uzunluğu otomatik olarak ayarlanabilir.
- Hem doğruluk hem de hesap verimliliği açısından güçlü bir yöntemdir.
Bu nedenle RKF yöntemi, mühendislik, fizik ve uygulamalı matematikte
başlangıç değer problemlerinin çözümünde yaygın olarak kullanılmaktadır.