1 / 49

BÖLÜM 7

BÖLÜM 7. ADİ DİFERANSİYEL DENKLEMLER. BÖLÜM 7.1. Giriş. 7 .1 Giriş. İlk derslerde Newton’un 2. yasasına dayanılarak türetilen ve düşen bir paraşütçünün hızının zamana göre değişimini ifade eden aşağıdaki bağıntıdan söz edilmişti;

teal
Download Presentation

BÖLÜM 7

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. BÖLÜM7 ADİ DİFERANSİYEL DENKLEMLER

  2. BÖLÜM 7.1 Giriş

  3. 7.1 Giriş İlk derslerde Newton’un 2. yasasına dayanılarak türetilen ve düşen bir paraşütçünün hızının zamana göre değişimini ifade eden aşağıdaki bağıntıdan söz edilmişti; Burada v paraşütçünün hızı (bağımlı değişken), t de zaman (bağımsızdeğişken) idi. Denklemdeki diğer büyüklükler ise: g yerçekimi ivmesi, c direnç katsayısı, m de düşen paraşütçünün toplam kütlesi ollarak belirtilmişti. Bu şekilde bilinmeyen bir fonksiyon ve onun türevleri cinsinden ifade edilen denklemlere diferansiyel denklemlerdenir. Diferansiyel denklemler bazen oran denklemleri olarak da adlandırılır. Çünkü, denklemden de görüldüğü üzere diferansiyel denklemler, değişkenlerin değişim oranlarını değişkenlerin kendilerinin ve çeşitli parametrelerin fonksiyonları cinsinden ifade eden bir yapıya sahiptirler. Bir çok fiziksel olay, en iyi, değişim oranları şeklinde matematiksel olarak formüle edilebildiğinden, bu tür denklemlere mühendisliğin her alanında sıklıkla rastlanır ve bu nedenle diferansiyel denklemler mühendislikte önemli bir yere sahiptirler.

  4. 7.1 Giriş Eğer diferansiyel denklem tek bir bağımsız değişkeniçeriyorsa, bu tür diferansiyel denklemlere adi diferansiyel denklemler(ADD) denir. Eğer birden fazla bağımsız değişkensöz konusu ise, bu tür diferansiyel denklemlere de kısmi türevli diferansiyel denklemler(KDD) denir. Diferansiyel denklemler mertebelerine göre de sınıflandırılırlar. Diferansiyel denklemdeki en yüksek türev mertebesi, diferansiyel denklemin mertebesi olarak adlandırılır. Örneğin, bir kütle yay sönüm elemanı sistemindeki kütlenin konumunu ifade eden aşağıdaki diferansiyel denklem ikinci mertebeden adi diferansiyel denklemdir. Bu denklemde m kütle, csönüm katsayısı, k da yay katsayısıdır. n. mertebeden bir diferansiyel denklem, n tane 1. dereceden diferansiyel şeklinde ifade edilebilir. Örneğin yukarıdaki denklemde aşağıdaki gibi yeni bir y değişkeni tanımlanırsa;

  5. 7.1 Giriş ve bu değişkenin t ye göre diferansiyeli alınırsa; elde edilir. (PT7.3) ve (PT7.4) ifadeleri (PT7.2) de yerine yazılırsa; veya elde edilir. (PT7.3) ve (PT7.5) (veya (PT7.6)) denklemleri, (PT7.2) denklemine eşdeğeriki adet birinci mertebeden diferansiyel denklemden oluşan bir diferansiyel denklem takımıdır.

  6. 7.1 Giriş Bir diferansiyel denklemin çözümü, bağımsız değişken ve parametrelerin, orijinal diferansiyel denklemi sağlayan özel bir fonksiyonudur. Bu kavramı aşağıdaki fonksiyon üzerinde inceleyelim; Şekilde grafiği görülen yukarıdaki fonksiyon dördüncü dereceden bir polinomdur. FIGURE PT7.3 Plots of (a) y versus x and (b) dy/dx versus x for the function

  7. 7.1 Giriş (PT7.12) eşitliğinin x ‘e göre diferansiyeli alınırsa; elde edilir. Elde edilen bu eşitlik orjinal fonksiyonun davranışını farklı bir şekilde ifade etmektedir. Orijinal fonksiyon her bir x değerine karşılık gelen y değerini gösterirken, elde edilen eşitlik her bir x değerine karşılık, bu noktada orijinal fonksiyonun sahip olduğu eğimi göstermektedir. Yandaki şekilde (Şekil PT7.3) orjinal fonksiyonun ve diferansiyelinin (eğiminin) x ‘e göre değişimi görülmektedir. FIGURE PT7.3: Plots of (a) y versus x and (b) dy/dx versus x for the function

  8. 7.1 Giriş Dikkat edilirse, orijinal fonksiyonun yatay olduğu (eğimin sıfır olduğu) noktalarda diferansiyel sıfır değeri almaktadır. Yine diferansiyelin ekstremum yaptığı noktalarda orijinal fonksiyon mutlak değerce en büyük eğimlere sahip olmaktadır. Yukarıda bir fonksiyondan hareket ederek buna karşılık gelen diferansiyel denklemin elde edilmesi incelendi. Halbuki bu konunun amacı tam tersine olarak, verilen bir diferansiyel denklemden hareketle, bu diferansiyel denkleme karşılık gelen orijinal fonksiyonun belirlenmesidir. Belirlenmek istenen orijinal fonksiyona diferansiyel denklemin çözümü adı verilir. FIGURE PT7.3

  9. 7.1 Giriş Şimdi yukarıda elde ettiğimiz diferansiyel denklemden hareketle orijinal denklemi belirlemeye çalışalım. Bunun için (PT7.13) eşitliğinin integralini almamız gerekir; İntegrasyon gerçekleştirilirse; elde edilir. Elde edilen çözüm, orijinal denklemle bir fark dışında özdeştir. Orijinal denklemin diferansiyeli alınırken yok olan 1 sabiti yerine integrasyon sonunda C sabiti gelmektedir. C sabitine integrasyonsabiti adı verilir. Böyle bir sabitin ortaya çıkması çözümün tek olmaması anlamına gelmektedir. Diğer bir ifadeyle sonsuz sayıda C değerine karşılık gelecek sonsuz sayıda çözüm söz konusudur.

  10. 7.1 Giriş Yandaki şekilde diferansiyel denklemi sağlayan olası altı çözüm görülmektedir. Dolayısıyla tek bir çözüm elde edilmek isteniyorsa, diferansiyel denklemle birlikte bazı ilave şartların da bilinmesi gerekir. Birinci mertebeden bir diferansiyel denklem için bu şart da bir tanedir ve başlangıç şartı(ya da değeri) adı verilir. Örneğin x=0 için y=1 şeklindeki bir başlangıç şartı bilinirse, diferansiyel denklemi sağlayan tek çözümü belirlemek mümkün olur. FIGURE PT7.4 Six possible solutions for the integral of dy. Each conforms to adifferent value of the constant of integration C.

  11. 7.1 Giriş Bu değerler (PT7.14) denkleminde yerine yazılırsa; elde edilir. Buradan C =1 olarak belirlenir.Belirlenen değer (PT7.14)’de yerine yazılırsa, diferansiyel denklemi ve başlangıç şartını sağlayan tek çözüm; şeklinde belirlenmiş olur. Başlangıç şartının bilinmesiyle çözüm fonksiyonu sanki başlangıç noktasından sabitlenerek tek bir çözümün elde edilmesini sağlamaktadır.

  12. 7.1 Giriş Fiziksel problemlere karşılık gelen diferansiyel denklemler için başlangıç değerleri önemli bir anlama sahiptir. Örneğin düşen paraşütçü probleminde başlangıç şartının sıfır olması, paraşütçünün düşey doğrultudaki başlangıç hızının sıfır olması anlamına gelmektedir. n. mertebeden bir diferansiyel denklem n tane birinci mertebeden diferansiyel denklemden oluşan denklem takımına karşılık geldiğine göre, n. mertebeden diferansiyel denklemin çözümü için n tane şarta ihtiyaç duyulacağı açıktır. Eğer bu şartların tümü başlangıç değeri için veriliyorsa, bu tür problemlere başlangıç değer problemleri, diğer noktalar için de verilmiş ise sınırdeğerproblemleri adını alır.

  13. 7.2Euler Yöntemi Bu bölümde; yapısındaki adidiferansiyeldenklemlerin sayısal yöntemlerle çözümü incelenecektir. Bu yöntemler şeklinde formüle edilebilecek genel bir yapıya sahiptir. Bu eşitliğe göre  eğimi ve yi eski değerinden yararlanılarak, h adımı sonundaki yi+1 yeni değeri belirlenebilir (Şekil 25.1); FIGURE 25.1 Graphical depiction of a onestep method.

  14. 7.2Euler Yöntemi Bu formül adım adım uygulanarak belirli bir aralıktaki çözüm adım adım belirlenmiş olur. Bütün tek adımlı yöntemlerde aynı temel düşünce kullanılmaktadır. Farklı yöntemler, eğimin belirlenmesindeki farklılıklardan kaynaklanmaktadır. Eğimin belirlenmesinde en basit yol, adımın başındaki eğimi doğrudan kullanmaktır. Adımın başındaki eğimin adım boyunca geçerli olduğu varsayımına dayalı olan bu yaklaşım Euler yöntemi olarak adlandırılır. FIGURE 25.1 Graphical depiction of a onestep method.

  15. 7.2Euler Yöntemi Birinci türev xi noktasında fonksiyonun eğimi için doğrudan bir tahmin oluşturur; Burada f(xi,yi), diferansiyel denklemin xi, yinoktasındaki değerinden başka bir şey değildir. Bu eğim değeri (25.1) eşitliğinde yerine yazılırsa; FIGURE 25.2 Euler’s method. Bu formül, Euler (Euler-Cauchy veya nokta-eğim) yöntemi olarak adlandırılır.

  16. 7.2Euler Yöntemi Example 25.1Aşağıdaki diferansiyel denklemi x=0 ilex=4 aralığında,h=0.5 adımı ile, x=0 day=1 başlangıç koşulu altında Euler yöntemiyle çözünüz. Not: Analitik çözüm Çözüm:

  17. 7.2Euler Yöntemi x=0 ve y(0)=1 için (0,1) noktasında çözüm fonksiyonunun eğimi; şeklinde belirlenir. Buna göre adımın sonundaki y değeri, Euler yönteminden; bulunur. Halbuki x=0.5 deki analitik çözümise şeklindedir. İki çözüm arasındaki hata oranı; kadar olmaktadır.

  18. 7.2Euler Yöntemi İkinci adım (2h=2*0,5=1) için çözüm; şeklinde elde edilir. x=1.0 deki analitik çözüm 3.0 olup, hata oranı; olmaktadır.

  19. 25.1 EULER YÖNTEMİ Tablo 25.1 başlangıç şartı x=0’da y=1 için integralinin kesin ve yaklaşık değerlerinin karşılaştırılması. Yaklaşık değerler, 0.5’lik bir adım büyüklüğünde Euler yöntemi kullanılarak hesaplanmıştır. Yerel hatalar, tek bir adım boyunca oluşan hataya göre olup Taylor serisi açılımı ile hesaplanmıştır. Genel hata, önceki ve şimdiki adımlar arasındaki farktır.

  20. 7.2Euler Yöntemi Dikkat edileceği üzere sayısal çözüm analitik çözümün genel karakteristiğini yakalamakta fakat hatanın mertebesi oldukça büyük olmaktadır. Hatanın büyüklüğü adım miktarının büyük olmasından kaynaklanmaktadır. Adım miktarı küçültüldükçe hata mertebesi de azalacaktır. FIGURE 25.3 Comparison of the true solution with a numerical solution using Euler’s method for the integralof dy from x 0 to x 4 with a step size of 0.5. The initialcondition at x 0 is y 1.

  21. 7.2Euler Yöntemi Örneğin adım miktarı yarıya düşürüldüğünde aşağıdaki çözüm elde edilmektedir. FIGURE 25.4 (a) Comparison of two numerical solutions with Euler’s method using step sizes of 0.5 and 0.25. 

  22. 7.3Heun Yöntemi Euler yönteminde hatanın temel bir kaynağı aralığın başlangıcındaki eğimin, tüm aralık boyunca geçerli olduğunun varsayılmasıdır. Bu yetersizliğin olumsuz etkisini azaltmak için iki basit düzeltme yapılabilir. Bunlardan biri, aralığın sadece başlangıcındaki eğimi kullanmak yerine, biri aralığın başlangıcındaki, diğeri de aralığın sonundaki olmak üzere iki noktadaki eğimin ortalamasının aralık boyunca geçerli olduğunun varsayılmasıdır. Bu yaklaşım Heun yöntemi olarak adlandırılmaktadır (Şek. 25.9). Euler yönteminde aralığın başındaki eğim; şeklinde belirlenmekteydi. Bu değer Heun yönteminde yi+1 değerinin lineer ektrapolasyonunda kullanılmaktadır;

  23. 7.3Heun Yöntemi Standart Euler yönteminden farklı olarak (25.13) eşitliği ile belirlenen y0i+1 değeri son değer değil, bir ara tahmin değeridir ve bu nedenle ‘0’ üst indisi ile gösterilmiştir. (25.13) eşitliği bu nedenle “tahmin denklemi (predictorequation)” olarak adlandırılır ve aralığın sonundaki eğim değerinin; şeklinde hesaplanabilmesi için bir tahmin değeri oluşturmaktadır. FIGURE 25.9 Heun yönteminin grafik olarak açıklanması. (a) Deneme (b) düzeltme.

  24. 7.3Heun Yöntemi Böylece (25.12) ve (25.14) eğimleri, aralık için ortalama eğimin belirlenmesinde kullanılabilir; Bu ortalama eğim artık yi den yararlanılarak lineer bir ekstrapolasyonlayi+1 in belirlenmesinde kullanılabilir; Bu son denklem “düzeltme denklemi (correctorequation)” olarak adlandırılır.

  25. 7.3Heun Yöntemi Böylece Heun yöntemi bir “tahmin etme-düzeltme (predictor-corrector approach)” yöntemi olarak tanımlanabilir; Tahmin etme (Şek. 25.9a): Düzeltme (Şek. 25.9b); Dikkat edileceği gibi (25.16) eşitliğinin her iki tarafında da yi+1 bulunduğundan, bu eşitlik iteratif bir yapıdadır ve eski değerler kullanılarak daha iyileştirilmiş yeni değerlerin belirlenmesinde kullanılabilir.

  26. 7.3Heun Yöntemi Tabii ki bu iteratif yaklaşımın analitik çözümün bulunmasını sağlamasını bekleyemeyiz. İyileştirme ancak yaklaşımın sahip olduğu belirli bir kesme hatası ile meydana gelecektir. Bu durum aşağıda bir örnekle incelenecektir. İterasyonun sonlandırılmasında daha önceden açıklandığı üzere aşağıdaki yaklaşık hata oranı hesabı kullanılacaktır; burada yj-1i+1veyji+1 değerleri sırasıyla bir önceki ve son iterasyon adımında belirlenen değerlerdir.

  27. 7.3Heun Yöntemi Example 25.5y’ =4 e0.8x – 0.5ydiferansiyel denklemini x=0 danx=4 e kadarh=1 adım miktarı ile integre ediniz. Başlangıç şartı: x=0 day=2. Çözüm:Sayısal çözümü gerçekleştirmeden önce analitik çözümün şeklinde olduğunu hatırlatalım ve bu çözümü Tablo 25.2 deki analitik değerlerin hesaplanmasında kullanıldığını belirtelim. İlk olarak (x0,y0) daki eğim değeri bulunabilir; Şimdi tahmin denklemini [(25.15) eşitliği] kullanılarak y’nin 1’deki ilk tahmini çözümü belirleyebiliriz;

  28. 7.3Heun Yöntemi Belirlenen bu değer normal olarak standart Euler yöntemiyle belirlenecek olan değerdir. Tablo 25.2 deki analitik değer ile karşılaştırılırsa bulunan yaklaşık çözümün %25.3 gerçek hata oranına sahip olduğunu görebiliriz. Tablo 25.2 İterasyonlar 1 15 (bu sütunu dikkate almayınız) x ytrue yHeunt(%) yHeun t(%) 0 2.0000000 2.0000000 0.00 2.0000000 0.00 1 6.1946314 6.7010819 8.18 6.3608655 2.68 2 14.8439219 16.3197819 9.94 15.3022367 2.68 3 33.6771718 37.1992489 10.46 34.7432761 3.17 4 75.3389626 83.3377674 10.62 77.7350962 3.18

  29. 7.3Heun Yöntemi Şimdi yi+1 deki çözümü iyileştirmek için y01 değeri aralığın sonundaki eğimin belirlenmesinde kullanılabilir; Belirlenen bu değer ile aralığın başındaki eğim kullanılarak x=0 dan 1 e kadar olan aralık boyunca ortalama eğim belirlenebilir; Bu eğim değeri, gerçek ortalama eğim değeri olan 4.1946 a daha yakın bir değerdir. Şimdi bu eğim değeri düzeltme denkleminde [(25.16) eşitliği] yerine yazılarak x=1 deki iyileştirilmiş yaklaşık çözüm belirlenebilir; Bu yaklaşık çözüm %8.18 gerçek hata oranına sahip olup, iterasyon yapılmaksızın belirlenen bu değer Euler yöntemi ile elde edilenden daha yüksek bir doğruluğa sahiptir.

  30. 7.3Heun Yöntemi Şimdi bulunan değerin (25.16) eşitliğinin sağ tarafına yeniden konularak hesabın tekrarlanmasıyla y1 çözümü iteratif olarak iyileştirilebilir; Elde edilen son çözüm % 1.31 yaklaşık hata oranına sahiptir. Bu çözüm kullanılarak iterasyona devam edilirse; Bu kez yaklaşık hata oranı %3.03 olmaktadır. Bu da iteratif yaklaşımlarda bazen hatanın büyüyebileceğini bize göstermektedir. Not: yukarıdaki hesaplamalar iteratif hesaplamalardır. Aşağıdaki tabloda Heun yöntemine göre hesaplanan 2,3 ve 4. değerleri hesaplayınız. Bu ödevdir ve sorumlusunuz.

  31. 7.3Heun Yöntemi Bu tür hata artışları özellikle büyük adım değerlerin kullanılması durumlarında meydana gelebilir ve bu gibi durumlar, ard arda gelen iterasyonlar sonucunda iteratif yaklaşımlarda mutlaka yakınsamanın meydana geleceği şeklinde genel bir sonuca ulaşmamızı engeller. Bununla birlikte yeterince küçük adımlar kullanıldığında, ard arda gelen iterasyonlar sonucu çözüm belirli bir değere yakınsayabilir. Yukarıdaki örnekte 15 iterasyon sonucunda 6.360865 değerine %2.68 yaklaşık hata oranı ile ulaşılmaktadır. Tablo 25.2 de her bir adım için tek iterasyon ve 15 iterasyon sonucu elde edilen değerlerin karşılaştırılması; Tablo 25.2İterasyonlar 1 15 (bu sütunu dikkate almayınız) x ytrue yHeunt(%) yHeun t(%) 0 2.0000000 2.0000000 0.00 2.0000000 0.00 1 6.1946314 6.7010819 8.18 6.3608655 2.68 2 14.8439219 16.3197819 9.94 15.3022367 2.68 3 33.6771718 37.1992489 10.46 34.7432761 3.17 4 75.3389626 83.3377674 10.62 77.7350962 3.18

  32. Runga-Kutta Yöntemleri 7.4Orta-Nokta Runga-Kutta Yöntemi Euler yönteminin iyileştirilmesi için yapılabilecek bir diğer iyileştirme Orta-nokta R-K yöntemi olarak adlandırılan, önce aralığın orta noktasının belirlenmesi için Euler yöntemi kullanılmakta, sonra da bu noktanın eğiminin aralık boyunca geçerli olduğu varsayımıyla aralık sonundaki değer belirlenmektedir. FIGURE 25.12 : Orta nokta yönteminin grafiksel açıklanması. (a) Eq. (25.25) and (b) Eq. (25.27).

  33. 7.4Orta-Nokta Runga-Kutta Yöntemi Eğim olarak aralığın orta noktasının eğimi kullanıldığından, yöntemin Euler yöntemine göre daha iyi sonuçlar vereceği tahmin edilebilir. FIGURE 25.12 Graphical depiction of the midpoint method. (a) Eq. (25.25) and (b) Eq. (25.27).

  34. 7.4.4.Mertebeden Klasik Runga-Kutta Yöntemi Çeşitli farklı Runge-Kutta yöntemleri vardır. Ancak bunların tümü aşağıdaki gibi ortak bir genel formda yazılabilir; burada (xi,yi,h) artımfonksiyonuolarak adlandırılır ve temsilieğimolarak düşünülebilir. Artım fonksiyonu aşağıdaki genel yapıda yazılabilir; burada a ‘lar katsayılar olup, k‘lar da çözüm fonksiyonunun çeşitli noktalarındaki eğimleri olarak düşünülebilir. Bunların belirlenmesi konusuna burada girilmeyecektir. Ancak uygulamada en çok kullanılan R-K yöntemi olan 4. mertebeden R-K yöntemi için katsayıların verilmesi ile yetinilecektir.

  35. 7.4.4.Mertebeden Klasik Runga-Kutta Yöntemi En çok kullanılan R-K yöntemi 4. derecedendir. Aşağıdaki form en yaygın olarak kullanılan 4. mertebeden klasik R-K yöntemidir ve artım fonksiyonu ile birlikte şeklinde yazılabilir. Buradaki k katsayıları aşağıdaki gibi belirlenmiştir; FIGURE 25.15 4. Dereceden R-K yöntemleriyle karşılaştırıldığında eğim tahminlerinin grafiksel olarak açıklanması.

  36. 7.4.4.Mertebeden Klasik Runga-Kutta Yöntemi denklemleri yukarıdaki gibi verilir. Yalnızca x’in fonksiyonu olan ADD’ler için klasik 4. dereceden R-K yöntemi Simpson’un 1/3 kuralına benzerdir. Ayrıca,, integral aralığı için iyileştirilmiş bir eğim elde etmek amacıyla geliştirilmiş olan çoklu tahmin açısından, 4. dereceden R-K yöntemi Heun yöntemiyle benzerdir. Aşağıdaki şekilde gösterildiği gibi, k’ların her biri bir eğimi göstermektedir. Bu durumda aşağıdaki eşitlik iyileştirilmiş bir eğim elde etmek için bunların bir ağırlıklı ortalamasıdır. FIGURE 25.15 Graphical depiction of the slopeestimates comprising the fourth-order RK method.

  37. 7.4.4.Mertebeden Klasik Runga-Kutta Yöntemi Example 25.7Aşağıdaki diferansiyel denklemleri x=0 ile 0.5 aralığında integre etmek için klasik 4.mertebeden R-K yöntemini kullanınız. (a) h =0.5 adımı ile y(0)=1 başlangıç şartı altında. (b) h =0.5 adımı ile y(0)=2 başlangıç şartı altında. Çözüm: (a) (25.40a) ila (25.40d) eşitlikleri kullanılarak k katsayıları hesaplanırsa: k1=8.5, k2=4.21875, k3=4.21875 ve k4=1.25 elde edilir (slayt 38). Bu değerler (25.40) eşitliğinde yerlerine yazılırsa; şeklinde analitik değerle çakışan sonuç bulunur.

  38. 7.4.4.Mertebeden Klasik Runga-Kutta Yöntemi (b)Bu kez aralığın başındaki eğim; olur. Bu değer kullanılarak orta noktadaki (0.5*h=0.25) y değeri ve eğim aşağıdaki gibi hesaplanabilir; Bu eğim değeri kullanılarak, orta noktada yeni bir eğim değeri hesaplanabilir;

  39. 7.4.4.Mertebeden Klasik Runga-Kutta Yöntemi Bu eğim kullanılarak aralığın sonundaki y değeri ve eğimi; şeklinde belirlenir. Bu dört eğimden yararlanılarak, ortalama bir eğim değeri ve bu eğim değeri yoluyla aralığın sonundaki son y değeri aşağıdaki gibi hesaplanır; Bu değer analitik değer olan 3.751152 ‘e çok yakındır.

  40. 7.5Diferansiyel Denklem Sistemleri (25.4) Mühendislik problemlerinin bir çoğunda tek bir diferansiyel denklemden ziyade bir diferansiyel denklem takımının çözümü gerekir. Böyle bir denklem takımı aşağıdaki gibi gösterilebilir; Bu şekildeki bir denklem sisteminin çözülebilmesi için n adet başlangıç şartına ihtiyaç vardır. Bu tür denklem takımlarının çözümü, tek bir denklemi çözmek için kullanılan yöntemlerin n denklem için genişletilmesiyle elde edilir. Burada her bir adım için tüm denklemler sıra ile çözüldükten sonra bir sonraki adıma geçilerek çözüm sürdürülür.

  41. 7.5.1Euler Yöntemi Example 25.9Aşağıdaki diferansiyel denklem takımını x=0ilex=2 aralığında verilen şartlar altında Euler yöntemi ile çözünüz. x=0 için y1=4, ve y2=6. Adım büyüklüğüh=0.5, Çözüm: Euler yöntemi her bir değişken için formüle edilirse; Not: Aşağıda y2 hesaplanırken, y1‘in yeni değeri kullanılmayacaktır.

  42. 7.5.1Euler Yöntemi İşlemler benzer tarzda devam ettirilirse (lütfen devam ettirelim); x y1y2 0 4 6 0.5 3 6.9 1.0 2.25 7.715 1.5 1.6875 8.44525 2.0 1.265625 9.094087 elde edilir.

  43. 7.5.2Runge-Kutta Yöntemi Herhangi bir R-K yöntemi Euler yöntemindekine benzer mantıkla kolaylıkla uygulanabilir. Burada dikkat edilmesi gereken durum, belirli bir noktada önce bütün denklemler için eğimlerin belirlenip, sonra y değerlerinin belirlenmesidir. Örneğin dördüncü mertebeden R-K yöntemi için, önce bütün denklemler için aralığın başındaki k değerleri (k1 ‘ler) bulunur, sonra aralığın orta noktasında bağımlı değişkenler hesaplanır. Benzer işlemler diğer k ‘lar için de tekrarlanır. Example 25.10Yukarıdaki örnekteki denklemleri dördüncü mertebeden R-K yöntemi ile çözünüz. Çözüm:Her bir denklem için aralığın başındaki eğimler belirlenirse;

  44. 7.5.2Runge-Kutta Yöntemi Bu değerlerden yararlanılarak aralığın orta noktasındaki y1vey2değeri; şeklinde belirlenir. Şimdi de bu değerler kullanılarak orta noktadaki eğimlerin ilk değerini hesaplamak için kullanılabilir;

  45. 7.5.2Runge-Kutta Yöntemi Şimdi bu değerler kullanılarak, aralığın orta noktasındaki ikinci y değerleri belirlenebilir; Son değerlerden yararlanılarak, aralık ortasındaki eğimlerin ikinci değeri belirlenir; Bu değerler yardımıyla da aralık sonundaki y değerleri belirlenir;

  46. 7.5.2Runge-Kutta Yöntemi ve bu değerler son noktadaki eğimleri hesaplamak için kullanılabilir; Her bir denklem için belirlenen tüm eğim değerleri ve aralık başındaki y değerleri kullanılarak aralık sonunda bulunmak istene y değerleri hesaplanabilir;

  47. 7.5.2Runge-Kutta Yöntemi İşlemlere benzer şekilde devam edilirse, elde edilen sonuçlar aşağıdaki tabloda görülmektedir.

  48. 7.5.2Runge-Kutta Yöntemi Lütfen bu değerleri işlemleri yaparak elde edelim; x y1y2 0 4 6 0.5 3.115234 6.857670 1.0 2.426171 7.632106 1.5 1.889523 8.326886 2.0 1.471577 8.946865 

More Related