1.4. RungËes ir Kuto metodas - techmat.vgtu.lt

More documents

Recommendations

Info

24 1 SKYRIUS. PDL PRADINIS UŽDAVINYS(1.38) algoritmas0, 5 0, 50 1⎧⎪⎨⎪⎩K 1 = F (t n , Y n ),K 2 = F ( t n + 1 2 τ n+1, Y n + 1 2 τ n+1K 1),Y n+1 − Y nτ n+1= K 2 .Pateiksime m-pakopio Rungės ir Kuto metodo algoritma.˛procedure Rungės ir Kuto metodo algoritmasbegin1. Pradinės sąlygos skaičiavimast = 0, n = 0for i = 1 to MY N[i] = U 0 [i]end for2. Sprendinio skaičiavimas nauju laiku t n+1while ( (t + τ n+1 ) ≤ T )for i = 1 to MY S[i] = Y N[i]end for3. Tarpinių vektorių K j skaičiavimasfor i = 1 to MK 1 [i] = f i (t, Y S)end forfor j = 2 to mfor i = 1 to MY N[i] = Y S[i]for k = 1 to j − 1Y N[i] = Y N[i] + τ n+1 b jk K k [i]end forend for
1.4. RUNGĖS IR KUTO METODAS 25endfor i = 1 to MK j [i] = f i (t + τ n+1 a j , Y N)end forend for4. Vektoriaus Y n+1 skaičiavimasfor i = 1 to MY N[i] = Y S[i]for j = 1 to mY N[i] = Y N[i] + τ n+1 σ j K j (i)end forend fort = t + τ n+1n = n + 1end while1.4.3. Rungės ir Kuto metodo koeficientu˛radimasAlgoritmo koeficientai a i , b ij , σ i parenkami taip, kad metodo aproksimacijostikslumas būtų didžiausias. Išnagrinėsime dvipakopį Rungės ir Kuto metoda. ˛ Iššio pavyzdžio bus aišku, kaip koeficientai randami ir bendruoju atveju.Imsime vienos diferencialinės lygties atvejį, kai M = 1. Lygčių sistemu˛analizėatliekama visiškai taip pat, be to, tik Rungės ir Kuto metodams, kuriu˛tikslumoeilė didesnė už ketvirtaj ˛ a, ˛ atsiranda papildomu˛lygčiu, ˛ kurias turi tenkinti metodokoeficientai. Kadangi paprastai randame ne konkretu˛vienintelį nurodyto tikslumoRungės ir Kuto metoda, ˛ o metodu˛šeima, ˛ priklausančia˛nuo keliu˛laisvu˛parametru,˛tai dažniausiai pavyksta patenkinti ir šias papildomas lygtis.Dvipakopio Rungės ir Kuto metodo formulėse yra keturi laisvi parametrai a 2 ,b 21 , σ 1 , σ 2 :K 1 = f(t n , y n ),K 2 = f(t n + a 2 τ, y n + b 21 τK 1 ),y n+1 = y n + τ(σ 1 K 1 + σ 2 K 2 ) .Perrašykime paskutiniaj ˛ a˛lygybę standartine forma:y n+1 − y n= σ 1 f(t n , y n ) + σ 2 f ( t n + a 2 τ, y n + b 21 τf(t n , y n ) ) .τ
Page 2 and 3: 22 1 SKYRIUS. PDL PRADINIS UŽDAVIN

1.4. RungËes ir Kuto metodas - techmat.vgtu.lt

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?

1.4. RungËes ir Kuto metodas - techmat.vgtu.lt