Luận văn tính ổn định của các phương pháp runge kutta ẩn
Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (953.37 KB, 83 trang )
B ộ GIÁO DỤC VÀ ĐÀO TẠO
TRƯỜNG ĐAI
• HOC
• s ư PHAM
• HÀ NÔI
• 2
= = m = =
Đ ổ VĂN LẼ
TÍNH ỔN ĐỊNH CỦA CÁC PHƯƠNG PHÁP
RUNGE- KUTTA ẨN
Chuyên ngành
: Toán giải tích
Mã số
: 60460102
LUẬN
VĂN THẠC
SĨ TOÁN HỌC
•
•
•
Người hướng dẫn khoa học: TS.Nguyễn Văn K hải
Hà Nội -2015
L Ờ I CẢM ƠN
Tác giả xin bày tỏ lòng biết ơn sâu sắc tới
T.s Nguyễn Văn Khải, người
đã định hướng chọn đề tài và tận tình hướng dẫn để tôi hoàn thành luận vãn.
Tôi cũng xin bày tỏ lòng biết ơn chân thành tới phòng Sau đại học, các
giảng viên dạy cao học chuyên ngành Toán giải tích, trường Đại học Sư phạm
Hà Nội 2 đã giúp đỡ tôi trong suốt quá trình học tập.
Nhân dịp này tôi cũng xin được gửi lời cảm ơn chân thành tới gia đình
bạn bè đã luôn động viên, cổ vũ, tạo mọi điều kiện thuận lợi cho tôi trong quá
trình học tập và hoàn thành luận văn.
Hà Nội, tháng 12 năm 2015
Tác giả
Đỗ Văn Lễ
LỜ I CAM ĐOAN
Luận văn được hoàn thành tại trường Đại học Sư phạm Hà Nội 2 dưới sự
hướng dẫn của TS. Nguyễn Văn Khải.
Tôi xin cam đoan luận văn “Tính ổn định của các phương pháp
R unge-K utta ẩn” không trùng với bất kỳ luận văn nào khác. Nếu sai tôi xin
hoàn toàn chịu trách nhiệm.
Hà Nội, tháng 12 năm 2015
Tác giả
ĐỖ Văn Lễ
MỤC LỤC
T rang
M Ở ĐẦU
1
CHƯƠNG 1: PHƯƠNG PHÁP RUNGE - KUTTA
3
1.1 Bài toán giá trị ban đầu
3
1.2 Phương pháp Euler
4
1.3 Phương pháp một bước
6
1.4 Phương pháp Runge - Kutta
16
1.4.1 Khái niệm và phân loại
16
1.4.2 Tính phù hợp của phương pháp Runge - Kutta
18
1.4.3 Bậc của phương pháp Runge - Kutta
19
1.4.4 Sự hội tụ của phương pháp Runge - Kutta
20
1.4.5 Tính ổn định của phương pháp Runge - Kutta
22
CHƯƠNG 2:TÍNH ỔN ĐỊNH CỦA CÁC PHƯƠNG PHÁP
RUNGE - KUTTA ẨN
27
2.1 Kiến thức chuẩn bị
27
2.2 Các khái niệm ổn định
52
2.3 Điều kiện đủ B-, BS- và BS- ổn định
57
2.3.1 BS- và BSI- ổn định
57
2.3.2 B -ổn định
62
2.4 Các kết quả ổn định YỚi một số lược đồ Runge - Kutta
65
2.4.1 Các phương pháp Gauss và Radau IA
65
2.4.2 Công thức Radau IIA
74
2.4.3 Các phương pháp Lobatto IIIC
77
KẾT LUẬN
79
TÀI LIỆU THAM KHẢO
80
1
M Ở ĐẦU
1. Lý do chọn đề tài
Phương trình vi phân là mô hình mô tả tốt các quá trình chuyển động
trong tự nhiên và kĩ thuật. Đe nghiên cứu phương trình vi phân, người ta
thường tiếp cận theo hai hướng: nghiên cứu định tính và giải số.
Mặc dù đã có lịch sử phát triển hàng trăm năm, do còn nhiều bài toán cần
giải quyết, giải số phương trình vi phân thường vẫn thu hút sự quan tâm mạnh
mẽ của các nhà toán học và các nhà nghiên cứu ứng dụng.
Trong giải số phương trình vi phân, người ta thường cố gắng tìm ra
những phương pháp hữu hiệu bảo đảm sự hội tụ, tính ổn định và tính chính
xác cao. Một trong các phương pháp đó là phương pháp Runge-Kutta ẩn, có
nhiều ứng dụng vì vậy tôi đã chọn đề tài: “Tính ổn định của các phương
pháp Runge - K utta ẩn”.
2. Mục đích nghiên cứu
Luận vãn sẽ nghiên cứu tính ổn định của các phương pháp Runge - Kutta ẩn.
3. Nhiệm vụ nghiên cứu
Nghiên cứu về phương pháp Runge-Kutta, tính phù họp, tính ổn định,
tính hội tụ của phương pháp.
4. Đối tượng và phạm vỉ nghiền cứu
Phương trình vi phân;
Phương pháp Runge-Kutta ẩn.
5. Phương pháp nghiên cứu
2
Phương pháp phân tích và tổng hợp tài liệu đã có từ đó hệ thống lại các
vấn đề liên quan tới đề tài.
6. D ự kiến đóng góp của đề tài
- Trình bày một cách hệ thống phương pháp Runge-Kutta;
- Trình bày tính ổn định của các phương pháp Runge-Kutta ẩn.
3
CHƯƠNG1
PHƯƠNG PHÁP RUNGE - KUTTA
1.1 Bài toán giá trị ban đầu
Định nghĩa 1.1.1. Giả sử G là một miền trong R 2 và hàm f :G —»R , hàm
khả vi liên tục y : \a,b\ —> R được gọi là nghiệm của phương trình vi phân
thường phân bậc nhất y = / ( * , 7 ) nếu:
■{x: A x ) \ e G
„
y 0 ) = /0 ,;K * ) )
v « k i].
Định nghĩa 1.1.2. Giả sử G là một miền trong R n+l và f : G ^ > R n, hàm khả
vi liên tục y : ịa,b]—>Rn được gọi là nghiệm của phương trình vi phân
thường y ' = f { x , y ) nếu:
\ (x: y ( x ) ] * G
„
V ,6 k 4 ].
Bài toán giá trị ban đầu: Cho hàm số y :[a,b] —> R n và hàm / : R x R ”
Giải bài toán giá trị ban đầu:
ị y' = f ( vx , y );
\y(a) = y 0.
Định lí 1.1.1. ( Định lỉ về sự tồn tại và duy nhất nghiệm)
Cho hàm số f : R x R n —»R n xác định liên tục trên miền:
(1.1)
R".
4
D = |(x , y ) : a < x < ố ,-0 0 < y t < +00, Vi = l,w j vớz'
Aữí/ hạn . Giả sử sự tồn
tại hằng số L sao cho:
||/ ( * > j) - / ( * > / ) ! ^
V ( x ,.y ) ,( x ,/) e D ,
khỉ đó với y ữ e R" luôn tồn tại duy nhất nghiệm của bài toán giá trị ban đầu
(1.1) sao cho y liên tục, khả vi v (x ,
e D.
1.2 Phương pháp Euler
Định nghĩa 1.2.1. Phương pháp Euler đế tìm nghiệm của bài toán Cauchy
( 1.2)
'y' = f ( x , y )
y { x°) = y ° ’
là việc tỉnh xấp xỉ y . thay cho giá trị của nghiệm ,y(* ) tại những mốc cách
đều với bước lưới h là:
X. = x a + h ,
yj+1= y j +hf ( x j , yj ),
j = 1,2,3...
j = 1,2,3...
Có 3 cách giải thích khác nhau cho công thức xấp xỉ của phương pháp Euler
là:
+) Thay đạo hàm bởi tỉ số sai phân:
h
+) Từ phương trình tích phân tương đương có:
5
^ 0 0 ® ;y(*o) +1
*0
xấp xỉ bởi quy tắc hình chữ nhật:
] f ( ỉ , y ( ỉ ) d ệ ~ h f ( x 0, y0)*0
+) Sử dụng công thức Taylor:
yCO=y ( x0+ h ) = .y(*o) +
h2
+2
y' (*0 +ỡh) ’
với 0 < 0 < 1 và bỏ qua những số hạng còn lạ i, nghĩa là xấp xỉ:
>>(*„)+ Ạy'(*„).
Mỗi cách giải thích đó mở ra một khả năng cải tiến phương pháp Euler.
Chẳng hạn thay vì sử dụng công thức hình chữ nhật, ta sử dụng công thức
hình thang chính xác hơn:
Ì f ( ậ , y ( ậ ) ) d ậ * 7 ( / ( w ( * o ) ) + / ( * 1 .:Kxi)))
*0
^
J 'i = J 'o + |( / C W o ) + / ( W i ) ) -
(1-3)
Quá trình được lặp lại và ta có định nghĩa sau đây:
Định nghĩa 1.2.2. Phương pháp Euler ẩn đổi với việc tìm nghiệm bằng sổ của
bài toán Cauchy (1.2) sẽ xấp xỉ y . tới nghiệm đúng y ( x ) tại những mốc
cách đều X. = x0 + j h , j > 0 bởi công thức:
6
y j+í= y j + 1^ { f { x . , y ]) + f { x . +ỉ,y.+ỉ)), j > 0 . (1.4)
Chú ỷ: Phương trình không tuyến tính (1.4) của phương pháp Euler ẩn có thể
được giải bởi dãy xấp xỉ liên tiếp với giả thiết hẳng số Lipschitz L của f và
bước h thỏa mãn Lh < 2.
Chứng minh. Ta phải giải phương trình (1.3) đối với y ì .
Đặt g(y) = y 0 + | ( / ( W o ) ) + / ( w ) ) > suy ra (1-3) được viết lại là:
= g ( j 1) ,h à m g là co YÌ
1^00 - g ( z ) I= 2 1 / ( w ) - f ( xr>z )\ ^ ^ - \ y - z \ < \ y - z \>
theo nguyên lí ánh xạ co ta có yì là tồn tại và duy nhất.
Định nghĩa 1.2.3. Phương pháp Euler cải tiến để tìm nghiệm của bài toán
Cauchy (1.2) là việc xấp xỉ y . tới nghiệm đúng
) tại những mốc cách
đều:X —xữ+ j h , j = 1,2,3... bởi công thức sau:
y J+1 =yj
+ 2 Ừ ( xJ’y J)+
+ hf(X],yj))\ J > 0.
1.3 Phương pháp một bước
Định nghĩa 1.3.1. Phương pháp một bước để tìm nghiệm bằng sổ của bài
toán giá trị ban đầu:
ịy ' =f(x,y)
1 ^ 0
) = y0
7
là việc xây dựng xấp xỉ y. tới nghiệm đủng y ị x .) tại các điểm cách đều:
x . = x 0+jh, j = 1,2,3...
với y . xác định bởi:
y J+1 = yj + h(p(x.,y.;h) với j > 0,
trong đó hàm:
ẹ?:G x(0,+oo)—»iỉ,
được xác định bởi hàm f { x , ỳ ) .
N hận xét:
+ Phương pháp Euler là phương pháp một bước với hàm ọ{x, y\ h) = f ( x , y ) .
+ Phương pháp Euler cải tiến là phương pháp một bước với hàm:
e G , giả sử rj = ĩj(ệ) là nghiệm duy nhất
của bài toán giá trị ban đầu:
77' = /(£77)>
Jì = y
,
với dữ liệu ban đầu (x,y), khỉ đó:
A (x,y,h) =
+ h) - 1 (*)] -
8
được gọi là sai số chặt cụt địa phương.
Phương pháp 1 bước được gọi là phù hợp nểu lim A (x,y\h) = 0 hội tụ đều
V ( x , ỵ ) e G và được gọi là phù hợp bậc p nếu: |л(дг,_у,/г)| < khp với
e ơ , VA > Q,k là hằng sổ.
Không giảm tính tổng quát ta giả sử / l à liên tục đều và bị chặn trên G.
Định lí 1.3.3. Phương pháp một bước là phù hợp khỉ và chỉ khi:
lim A (x,y;h) = f ( x , y ) ,
hội tụ đều với v ( x ,y ) e G .
Chứng minh. Giả sử / bị chặn, ta có:
t
ri(x + t ) - r i ( x) = |
t
7
'(^ + 5)í/5 = |/ ( л : + £,77(л; + £))^5 —>0 khi í —>0 hội tụ
0
0
đều với v ( j c , j ) e ơ .
Do / liên tục đều nên:
ịỵj]'{x + t)-r i'{ x Ỵ ịd t < max ìj'(x + t ) - ì ] ' ( x ) I
0
= m ax|/(jc + í,7(jc + í ) ^ - / ( j c , 7 ( x ) |^ 0
hội tụ đều với v (jc ,j) e G .
Tù đó ta có:
9
A (x,y,h) +
= —j[ 7 ( * + 1) -
h0
hội tụ đều với
T ]{x Ỵ ịd t
—»0
v(jc,>>) e G.
Vậy:
Định lí 1.3.4. Phương pháp Euler là phù hợp. Neu hàm f có đạo hàm riêng
ỉiên tục trên G thì phương pháp Euler là phù hợp bậc nhất.
Chứng minh:
Theo phương pháp Euler thì: ọ{x, y, h) = / ộ t,y) suy ra phương pháp Euler là
phù hợp.
Do / có đạo hàm riêng liên tục nên từ rj' = f(ệ,7]) suy ra TỊ có đạo hàm
riêng liên tục cấp hai và
TJ" = f x{Z,Tl)+ f
’Tì) ■
(1.5)
Theo công thức khai triển Taylor ta có:
^J]{x +h )-ĩ](x)\-r]\x)
h
T]"ịx + 6hì)I < Kh, với 0 < 9 < 1
và K là hằng số sao cho:
f , ( ỉ , n ) +f A ĩ > n ) f ( ỉ , v h 2 K .
10
Vậy phương pháp Euler là phù hợp bậc nhất.
Định lí 1.3.5. Phương pháp Euler cải tiến là phù hợp. Nếu hàm f có đạo
hàm liên tục đến cấp hai trên G thì phương pháp Euler cải tiến là phù hợp
bậc 2.
Chứng minh.
Ta có:
+ f ( x + h,y + hf(x,ỳ))~\h^ > f ( x , y ) (hội tụ đều).
Do / có đạo hàm liên tục đến cấp hai theo (1.5) ta suy ra TỊ có đạo hàm liên
tục đến cấp ba:
V = f„ (Ẹ,rị) +24 (ệ,rj
+4
{ Ẹ , v ) f 2+
+ f y( ệ, ĩ ỉ ) f x(ệ,ĩỉ)+ / ' ( ệ ^ í ị ệ ^ Theo công thức khai triển Taylor ta có:
ĩj(x + h ) ~
7 0
h2
) - hĩ]\x) - y (*) = Ỵ\ĩj"(x + 0h)\
^
với 0 < 9 < 1 và K J được chọn sao cho:
If ( x + h, y + k ) ~ f ( x , y ) - hfy (*,
trong đó K 2 là một hằng số thỏa mãn:
Lấy k = hf( x, y) và theo (1.5) ta có:
- hfy (*, j ) | < ^ K 2(|ã| + |ả:|)2,
11
If ( x + h, y + h f 0,;>0) - f ( x , y ) - hrj"(x)\ < ị x 2(1 + K ữỴ h 2,
với K ữ là hằng số sao cho 1
/ <
K 0, do đó:
ọ(x,y;h)-f(x,y)-^rj\x)
(1.7)
Từ (1.6) và (1.7) ta có:
K ^ K ^ + K ,)' h1
Vậy phương pháp Euler cải tiến là phù hợp bậc hai.
Định nghĩa 1.3.6. Giả sử ta chia đoạn [a,b] bởi các mốc cách đều:
x . = x ữ+ j h j = 0,1,2,...,«
xa= a , x n =b,
xấp xỉ tới giá trị y . thay cho giá trị của nghiệm >>(■*.) của bài toán giá trị
ban đầu:
\y ' = f ( x , y )
^ (X ) = y0
bằng phương pháp một bước, khỉ đó:
e. = e.Ọì) = y.
j = 0,1,2,...,«
được gọi là sai số toàn cục và
E = EỌÌ) = max e i h )
}=0....n 1
12
được gọi là sai số toàn cục lớn nhất.
Định nghĩa 1.3.7. Phương pháp một bước được gọi là hội tụ đều nếu
XìmEỢí) = 0 và được gọi là hội tụ bậc p nếu E(h) < Hhp,\fh > 0 và H là
0
một hằng số nào đó.
Bổ đề 1.3.8. Cho 1^1 là dãy sổ thực thỏa mãn:
(với hằng số A>0 và B > 0 ), khi đó ta luôn có:
ị \ < \ ị y - ì )
j = 0,1,2.
Chứng minh. Ta chứng minh bằng quy nạp.
Đánh giá là đúng với j > 0.
Áp dụng bất đẳng thức 1+ A < eA ta có:
< (1 + A ) \ ị y + (1 + A ) * ( e " - 1) + B
<\Ịt \eu" u + - { é j'n* - ì ) .
Ả
Vậy đánh giá đúng với j +1.
Định lí 1.3.9. Giả sử hàm (p được mô tả trong phương pháp một bước là liên
tục đối với biến h và thỏa mãn điều kiện Lipschitz:
\(pịx,y,h) -
< M \y - z\
13
với
v(jtj) ,( í ,z ) e G ,V Ã
đủ nhỏ, khi đó phương pháp một bước hội tụ nếu
và chỉ nếu nó phù hợp.
Chứng minh.
*) Phù họp suy ra hội tụ:
Giả sử phương pháp một bước là phù hợp. Xét 2 sai số liên tiếp ta có:
e» - ej = [ y » - y j ] - [ y ( x» - * M ) ]
= h
Do đó:
e -4.1
7+1
( 1.8)
~ ej < h ị M ị y . - y [ x . ) ị + c(h)ị
với c(Ã) = max à (^ ,j(x );/ỉ) thỏa mãn c (/ỉ)—>0 khi h —>0
a
I
'
' '
'
'
Do phương pháp 1 bước phù hợp nên suy ra:
e,-+i < ịl + hM) e . +hc(h) , j = 0,1,2...,«
Áp dụng bổ đề 1.3.8 với A - hM và B = hcịh) và sử dụng e0 = 0 ta có:
< C^ - [ e M{XrX° v ớ i ý = 0,l,...,w
suy ra:
E(h) <
Vậy phương pháp một bước là hội tụ.
*) Hội tụ suy ra phù hợp:
_ !^4°o
(1.9)
14
Giả sử phương pháp 1 bước là hội tụ nghĩa là khi h —» 0 thì xâp xỉ:
y j+í= y J +h
(1.10)
hội tụ đến nghiệm của bài toán giá trị ban đầu.
Đ ặ t:g (x ,j) = ọ ị x , j;0 ) . Từ định lý 1.3.3 ta thấy phương pháp một bước là
phù hợp đối với bài toán Cauchy:
Từ đó ta thấy rằng tính phù hợp dẫn đến sự hội tụ , công thức gần đúng (1.10)
cũng hội tụ đều đến nghiệm của (1.11) và nghiệm của bài toán Cauchy là
trùng nhau.
Do đ ó ta c ó : /( * 0,.y0) = g(x0,.y0) với V(*0, j 0) e G
suy ra:
hợp.
Định lí 1.3.10. Giả sử rằng phương pháp một bước thỏa mãn các giả thiết
của định lỉ 1.3.9 và nó là phù hợp bậc p nghĩa là:
< Khp
khi đó:
nghĩa là sự hội tụ cũng có bậc p .
Chứng minh.
j = 0,l,2...,n
15
Do
cỌì) = m.dĩĩỏịk{x,y{xy,h^ị
nên từ (1.9) suy ra:
o) _ A <
<
M v
'
<*,-*.) _ ỵ.ỵp
m
Hệ quả 1.3.11. Phương pháp Euler và Euler cải tiến là hội tụ.
Neu f khả vỉ thì phương pháp Euler hội tụ bậc một.
Neu f khả vi liên tục đến cấp hai thì phương pháp Euler hội tụ bậc hai.
Ta có thể xây dựng phương pháp một bước với bậc cao hơn như sau:
Cho một tập hợp các số thực:
s, ; / = 2,3,...,m
c(i;i = l,2 ,...,l-l;l = 2,3,...,w
a , ; ỉ = ỉ,2,...,m
Các hệ thức:
= /(W j)
K
=
f ( XJ
+
K = f ( Xj +
S 2k
y
j
+
c 21k l h )
S 3h > J ; +
c 31k lh +
c 32k 2 h )
k m = Jf iVx .j + s m h,v.
+ h ^Ỵ c m i kí)S
7y J
i=1
được tính toán 1 cách đệ quy và dãy xấp xỉ được xác định bởi:
m
y j+i = y j + h l l a ik r
i=1
16
Phương pháp Euler được mô tả bởi m =l, «! = 1.
Phương pháp Euler được mô tả bởi m =l, s2 = 1,c21 = 1,a 1= a 2 = —.
Mục đích cơ bản trong việc xây dựng phương pháp bậc cao hơn là cho 1
số m, xác định các hệ số cao cho bậc của tính phù hợp và hội tụ là lớn nhất có
thể.
Năm 1895, Runge giới thiệu phương pháp 1 bước để giải phương trình vi
phân và được mở rộng hệ phương trình vi phân bởi Kutta năm 1901. Phương
pháp này được gọi là phương pháp Runge-Kutta. Đây là phương pháp thành
công nhất của lớp các phương pháp một bước và ngày nay được sử dụng rộng
rãi.
Định nghĩa 1.3.12. Phương pháp Runge-Kutta để tỉnh nghiệm bằng số của
bài toán Cauchy (1.2) xây dựng bởi sự xấp xỉ y. tới nghiệm chỉnh xác ,y(* )
tại các mốc cách đều:X = x0 + jh với j = 1,2,3,... với bước h bằng cách sử
dụng phương pháp 1 bước với bậc cao hơn:
K =/ ( w , )
K=f(xj +ị,yj +ịK)
K = f{x.+ ị,y.+ ị\)
K = f ( xj + K y Ị +hk3)
y ì« = 6 (*1 + 2K + 2K + k, ) + y . .
1.4 Phương pháp Runge - K utta
1.4.1 K hái niệm và phân loại
17
Định nghĩa 1.4.1. Dạng tổng quát của phương pháp Runge —Kutta s nấc cho
bài toán giá trị ban đầu (1.1) là:
y,„ = y .+ hĨLb,K
“
k . = f
X + c.h ,y
^
.
>
+ h ỵ^ a..k.
J=1
-
(L 1 2 )
,i = l,s
'J
ta giả sử điều kiệnc. =ỵ^a ,i = l,s luôn được thỏa mãn.
.7=1
Để thuận tiện cho việc trình bày ta ghi các hệ số xuất hiện trong các công thức
(1.12) vào bảng được gọi là bảng Butcher:
cs
an
a2l
al2 - a u
... a2s
«22
a n
a.i
b i
...
ass
b2 ... bs
Ký hiệu: A = [ ạ .] là ma trận cấp s x s , c = ịcl,c2,...,c Ỵ ,b = ịb1,b2,...,b Ỵ là
các vec tơ s chiều. Khi đó phương pháp Runge-Kutta (1.12) được xác định
duy nhất bởi bảng Butcher:
Nhân
xét: Nếu ta đăt:
•
•
18
s
Yi =y»+hH aíìf ( x»+cjh’Yj)’
i = 1,5
k, = f ( x. + c№ ) >
i = ỉ,s .
Khi đó phương pháp Runge-Kutta có thể viết dưới dạng:
r
Phân loai:
+ Nếu a = 0, V/ > ỉ,ỉ = 1,5 hay A là tam giác dưới chặt thì phương pháp
Runge-Kutta (1.12) gọi là phương pháp Runge-Kutta hiển (hay phương pháp
Runge-Kutta cổ điển).
+ Nếu a = 0, V/ > ỉ,ỉ = 1,5 hay A là tam giác dưới chặt thì phương pháp
Rune-Kutta (1.12) gọi là phương pháp Runge-Kutta nửa ẩn.
+ Trong trường hợp còn lại thì phương pháp Runge - Kutta (1.12) được
gọi là phương pháp Runge-Kutta ẩn. A không là ma trận tam giác dưới.
1.4.2 Tính phù hợp của phương pháp Runge-K utta
Xét bài toán giá trị ban đầu:
y' = f ( x, y ) ,
y( a) = y„>
Xét lớp phương pháp tổng quát:
Ẹ a ^ = h ệ f{ y n+k,yn+k- „ - ,y n,\-,h )
19
khi:
s
k = l , al = l , a0= - l , ệ f ( y ri,xri;h) = ỵ b iki,
i=1
k.i = -ỉf
(
V
X n + c .i h 3, yy n + h y a ij. k .j
,
i = 1,5,
j= '
thì ta có phương pháp Runge-Kutta dạng (1.12). Khi đó điều kiện cần và đủ
để phương pháp Runge-Kutta (1.12) phù hợp là:
2 > ;= °
j= 0
o
0/ừ(* ),* ;0) = f { \ , y { x j )
° ẳ i=*l ' =1
Vậy với 2^ ồ. = 1 thì phương pháp Runge-Kutta là phù hợp.
j=l
1.4.3 Bậc của phương pháp Runge-K utta
Định nghĩa 1.4.2. Sai số chặt cụt địa phương của phương pháp Runge-Kutta
(1.12) tại
X
J
tại T J được xác định bởi công thức:
T +l ^ y ( x a+l) - y a+l
với giả thiết y n = y ( x n),Y.= y ( x n + cJBh ) j = 1,5
s
hay là Tn+l := y m+l - y n- hỴ^bk.
i=1
20
Định nghĩa 1.4.3. Bậc của phương pháp Runge-Kutta (1.12) là số nguyên p
lớn nhất sao cho T J =0(Aì>+1^.
Gọi y +1 là giá trị tính theo phương pháp Runge-Kutta tại X +ì, y = y ( x ) là
giá trị đúng, khi đó y n+l = y n+ ệ f ( y n,xn; h ) .
Do vậy: y ( x n+l) - y n+1=Tn+1 .
Khai triển Taylor tại X ta có:
Tn+l = y(xn+1) - y { x M) ~ ệ f ( y ( x n),xn;h)
= hy \ x n) - h f (xn,y(xn)) + 0( h1)
= 0 (h2),
như vậy mọi phương pháp Runge-Kutta (1.12) phù hợp đều có bậc p > 1 .
1.4.4 Sự hội tụ của phương pháp Runge-K utta
Định nghĩa 1.4.4. Phương pháp Runge-Kutta (1.12) hội tụ khỉ và chỉ khỉ
phương pháp Runge —Kutta (1.12) là phù hợp và thỏa mãn điều kiện nghiệm.
Đa thức đặc trưng của phương pháp Runge-Kutta (1.12) là:
p(ệ) = z - h
luôn thỏa mãn điều kiện nghiệm.
Vì vậy phương pháp Runge-Kutta (1.12) hội tụ khi và chỉ khi phương pháp
Runge-Kutta (1.12) thỏa mãn điều kiện phù hợp. Nghĩa là phương pháp
Runge-Kutta (1.12) hội tụ khi và chỉ khi Ỳ^b - l .
i
Ví du
=1
21
Ví dụ 1: Phương pháp Runge - Ku tta hiển 2 nấc cho bởi bảng Buttcher:
0
0
0
1
2
ỉ
2
0
0
1
(Phương pháp Euler cải tiến: ỏ, = 0 A =1>C2 = | ) .
Ví dụ 2: Phương pháp Runge - Ku tta nửa ẩn với bảng Butcher:
0
1
2
1
2
0
1
2
0
1
2
1
2
(Nhận xét: s = 2; p = 2 ta có phươELg pháp Runge - Kutta 2 nấc cấp 2).
Ví dụ 3: Xét sự hội tụ của phương pháp Runge - Kutta cho bởi:
K = f { xn ơ n)
f l
1 A
h
y„+l-y„ = - ị ( K + ?>K) >với: K = f x n + ị ^ y n + ị hk1
V
k
T
1
J
_
_
•
7
•
L òi giải:
Ta có bảng Butcher là:
3
=
f
./
f
X
n
-3
2
+ —h,y
3
i
2
+ —hk,
3
