Các phương pháp giải bài toán cauchy đối với phương trình vi phân thường

Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

1

LỜI CẢM ƠN.

Để hoàn thành khoá luận tốt nghiệp em đã nhận được sự dìu dắt, chỉ bảo
và tạo điều kiện giúp đỡ của các thầy cô trong khoa Toán nói chung và trong
tổ Giải tích nói riêng, đặc biệt là sự hướng dẫn, chỉ bảo và giúp đỡ hết sức tận
tình của thầy giáo TS.Khuất Văn Ninh.
Qua đây, em xin bày tỏ lời cảm ơn chân thành tới thầy giáo TS.Khuất
Văn Ninh. Em xin gửi lời cảm ơn sâu sắc tới các thầy cô giáo trong tổ Giải
tích, các thầy cô giáo trong khoa Toán, cảm ơn gia đình, bạn bè và các bạn
sinh viên quan tâm và đóng góp ý kiến cho đề tài của em.

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

2

LỜI CAM ĐOAN.

Kết quả của đề tài này là do sự nỗ lực cố gắng tìm tòi của bản thân. Em
xin cam đoan kết quả nghiên cứu của em không trùng với kết quả của các tác
giả khác.

Hà Nội ngày 18/5/2009.
Sinh viên:
Phạm Thị Hoa.

Phạm Thị Hoa K31E-SPToán


3

Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

MỤC LỤC
LỜI CẢM ƠN

1

LỜI CAM ĐOAN

2

LỜI NÓI ĐẦU

5

NỘI DUNG KHOÁ LUẬN

7

Chƣơng 1: Kiến thức cơ sở.

7

I. Các khái niệm

7

1. Số gần đúng

7

2. Sai số

7

3. Sai phân

8

II. Khái quát về phương trình vi phân

8

1. Một số khái niệm

8

2. Bài toán Cauchy đối với phương trình vi phân thường cấp 1

8

3. Bài toán Cauchy đối với hệ phương trình vi phân thường cấp 1

10

4. Điều kiện Lipschitz

10

Chƣong 2: Các phƣơng pháp giải gần đúng bài toán Cauchy đối với
phƣơng trình vi phân thƣờng.

11

I. Các phương pháp giải tích.

11

1. Phương pháp lặp đơn.

11

2. Phương pháp xấp xỉ liên tiếp Picard.

13

3. Phương pháp chuỗi số nguyên.

15

II. Các phương pháp số.

16

1. Phương pháp Euler.

16

2. Phương pháp Euler Cauchy

18

3. Phương pháp Runge Kutta.

20

4. Phương pháp Adams.

25

III.Ứng dụng của tin học để giải phương trình vi phân thường.
Phạm Thị Hoa K31E-SPToán

29


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

4

1. Ứng dụng của chương trình MapleV

29

2. Ứng dụng của ngôn ngữ lập trình Pascal.

30

Chƣơng 3: Các bài tập ứng dụng.

38

I. Các bài tập ứng dụng của các phương pháp giải tích.

38

II. Các bài tập ứng dụng của các phương pháp số.

45

KẾT LUẬN

53

TÀI LIỆU THAM KHẢO

54

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

5

LỜI NÓI ĐẦU
Thế kỷ XXI là thế kỷ bùng nổ của công nghệ thông tin, ứng dụng của
công nghệ thông tin có đóng góp to lớn và hiệu quả trong mọi mặt của đời
sống. Và cũng từ rất lâu tin học đã được ứng dụng vào trong môn Toán. Có
những số liệu tính toán quá cồng kềnh và những bài toán phức tạp chúng ta
không thể giải bằng tay được nhưng nếu dùng các lập trình trên máy vi tính
thì chúng ta có kết quả rất nhanh gọn và chính xác.
Các bạn sinh viên đã được học môn phương trình vi phân từ kì II năm
thứ ba, vì thế các bạn đã rất quen thuộc với dạng toán tìm nghiệm đúng của
bài toán Cauchy đối với phương trình vi phân thường. Thế nhưng có nhiều
trường hợp nghiệm đúng của các phương trình vi phân không thể tìm được.
Bởi vậy để tìm nghiệm của chúng, ta phải áp dụng các phương pháp gần đúng
khác nhau. Ở mỗi phương pháp chúng ta có thể dùng lập trình Pascal hay sử
dụng thuật toán MapleV để giải các bài toán này.
Với mong muốn học hỏi tích luỹ thêm cho mình những kỹ năng và kinh
nghiệm khi tiếp cận với ứng dụng của công nghệ thông tin váo việc giải toán
đồng thời để hiểu sâu hơn về phương trình vi phân em mạnh dạn chọn đề tài
là: “Các phương pháp giải bài toán Cauchy đối với phương trình vi phân
thường”.
Nội dung của khoá luận gồm 3 chương:
Chƣơng 1: Kiến thức cơ sở.
Chương này nhằm trình bày các khái niệm và định lý cơ bản nhất về
các vấn đề có liên quan đến nội dung trong chương 2 sẽ trình bày.

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

6

Chƣơng 2: Các phƣơng pháp giải gần đúng bài toán Cauchy đối
với phƣơng trình vi phân thƣờng.
I. Các phương pháp giải tích.
1. Phương pháp lặp đơn.
2. Phương pháp xấp xỉ liên tiếp Picard.
3. Phương pháp chuỗi số nguyên.
II. Các phương pháp số.
1. Phương pháp Euler.
2. Phương pháp Euler-Cauchy.
3. Phương pháp Runge-Kutta.
4. Phương pháp Adams.
III. Ứng dụng của tin học để giải phương trình vi phân thường.
1. Ứng dụng của chương trình MapleV
2. Ứng dụng của ngôn ngữ lập trình Pascal.
Chƣơng 3: Các bài tập ứng dụng.
I. Các bài tập ứng dụng của các phương pháp giải tích.
II. Các bài tập ứng dụng của các phương pháp số.
Mặc dù đã có nhiều cố gắng song do thời gian có hạn và điều kiện
nghiên cứu còn hạn chế đồng thời kiến thức của bản thân người làm khoá luận
còn chưa vững nên khoá luận không tránh khỏi những thiếu sót, em rất mong
nhận được sự quan tâm góp ý của các thầy cô giáo, các bạn sinh viên cũng
như các bạn đọc quan tâm đến vấn đề này để khoá luận được hoàn thiện hơn.
Hà Nội 18/5/2009.
Sinh viên:

Phạm Thị Hoa.

Phạm Thị Hoa K31E-SPToán

Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

7

CHƢƠNG 1
KIẾN THỨC CƠ SỞ.

I. Các khái niệm:
1. Số gần đúng
Trong tính toán, ta thường phải là việc với các giá trị gần đúng của các
đại lượng. Ta nói a là số gần đúng của a * , nếu a không sai khác a * nhiều.
2. Sai số
a) Sai số tuyệt đối, sai số tương đối.
+) Sai số tuyệt đối: Đại lượng  : a  a* gọi là sai số thật sự của a .
Do không biết a * nên ta cũng không biết  . Tuy nhiên ta có thể tìm được

a  0 , gọi là sai số tuyệt đối của a , thoả mãn điều kiện:
a  a*  a hay a  a  a*  a  a.
+) Sai số tương đối: đại lượng  a :

a
.
a

b) Sai số thu gọn:
Một số thập phân a có dạng tổng quát như sau:

a     p10 p   p110 p1  ...   ps10 ps 

trong đó 0 ≤  i ≤ 9  i  p  s,..., p  1 ;  p  0 là những số nguyên. Nếu
s  , a là số thập phân vô hạn. Thu gọn một số a là vứt bỏ một số các chữ

số bên phải a để được một số a ngắn gọn hơn và gần đúng nhất với a .
c) Sai số tính toán: Các số vốn đã có sai số, còn thêm sai số thu gọn nên
khi tính toán sẽ xuất hiện sai số tính toán.

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

8

3. Sai phân
+) Sai phân: Giả sử

f : R  R là một hàm số cho trước và

h  const  0 . Ta gọi sai phân cấp 1 của

f  x

là đại lượng

f  x   f  x  h   f ( x) .
+) Tỷ sai phân cấp 1 của f  x  là

f(x)
.

h

Một cách tổng quát n f ( x)    n1 f ( x)  ,  n  1 , 0 f ( x) : f ( x) .
II. Khái quát về phƣơng trình vi phân:
1. Một số khái niệm:
a) Phương trình vi phân cấp 1
Phương trình vi phân cấp 1 có dạng tổng quát: F ( x, y, y ')  0 (a) trong
đó hàm F xác định trong miền D   3 . Nếu trong miền D , từ phương trình
(a) ta có thể giải được y ' : y '  f ( x, y ) thì ta được phương trình vi phân cấp 1
đã giải ra đạo hàm.
b) Phương trình vi phân cấp n:
Phương trình vi phân cấp n có dạng tổng quát F ( x, y, y ',..., y ( n) )  0 (b).
Hàm F xác định trong một miền G nào đấy của không gian  n2 . Trong
phương trình (b) có thể vắng một số biến x, y, y ',..., y ( n1) nhưng y ( n ) nhất thiết
phải có mặt. Nếu từ (b) ta giải ra được đạo hàm cấp cao nhất, tức là phương
trình (b) có dạng: y ( n)  f ( x, y, y ',..., y( n1) ) thf ta được phương trình vi phân
cấp n đã giải ra đối với đạo hàm cấp cao nhất.
2. Bài toán Cauchy đối với phƣơng trình vi phân thƣờng cấp 1
Xét bài toán

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

9

t, x   R  0,T    x0  r, x0  r  .
Trong đó x (t ) là hàm một biến xác định trên  0,T  với  0,T  cho
trước, hàm f (t , x) và x0 cho trước được gọi là bài toán Cauchy đối với

phương trình vi phân thường cấp 1, điều kiện (2) được gọi là điều kiện
Cauchy hay điều kiện ban đầu.
a) Định lý 1(định lý tồn tại nghiệm)
Xét bài toán (1-2),  t , x   R  0,T    x0  r , x0  r  .
Nếu f (t , x) là hàm liên tục trên hình chữ nhật R(r  0 cố định) thì tồn
tại ít nhất một nghiệm x (t ) của phương trình (1) thỏa mãn điều kiện (2) tức là
x (t ) là nghiệm của bài toán (1-2).

b) Định lý 2(định lý duy nhất nghiệm)
Xét bài toán (1-2). Nếu f (t , x) là hàm liên tục trên hình chữ nhật
R(r  0 cố định) và f (t , x) thoả mãn điều kiện Lipschitz tho biến x trên hình

chữ nhật R tức là:
.
Trong đó N là hằng số (gọi là hằng số Lipschitz) thì nghiệm của bài
toán (1-2) xác định là duy nhất.
c) Định lý (định lý tồn tại và duy nhất nghiệm)
Xét bài toán (1-2),  t , x   R  0,T    x0  r , x0  r  .
Hàm f (t , x) xác định trong R(r  0 cố định) thoả mãn 2 điều kiện :
a. f (t , x) liên tục trên R và do R đóng
ới M  max f (t , x) .
 t , x R

Phạm Thị Hoa K31E-SPToán

và bị chặn cho nên


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

10

b. f (t , x) thoả mãn điều kiện Lipchitz
là hằng số thì tồn tại duy nhất nghiệm
x (t ) của bài toán (1-2) xác định trên [0,T].

3. Bài toán Cauchy với hệ hai phƣơng trình vi phân:

4. Điều kiện Lipschitz.
Ta nói rằng trong miền G hàm f ( x, y ) thoả mãn điều kiện Lipschitz

 

theo biến y nếu tồn tại hằng số L > 0 sao cho đối với hai diểm x, y  G,

 x, y   G bất kì, ta có bất đẳng thức
.

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

11

CHƢƠNG 2
CÁC PHƢƠNG PHÁP GIẢI GẦN ĐÚNG BÀI TOÁN CAUCHY ĐỐI
VỚI PHƢƠNG TRÌNH VI PHÂN THƢỜNG
Các phương pháp giải phương trình vi phân chia làm 2 nhóm:
*Nhóm các phương pháp giải tích cho phép tìm nghiệm gần đúng dưới

dạng biểu thức giải tích.
*Nhóm các phương pháp số cho phép tìm nghiệm dưới dạng bảng.
Sau đây, ở mỗi nhóm chúng ta sẽ xét một vài phương pháp cụ thể.
I. Các phƣơng pháp giải tích.
1. Phƣơng pháp lặp đơn.
a) Nội dung phương pháp.
Xét bài toán giá trị ban đầu sau đây:
y' 

dy
 f ( x, y )
dx

(1)

y '( ) =

(2)

Giả sử f ( x, y ) là hàm liên tục trên
R=

].

Khi đó bài toán (1), (2) tương đương với phương trình tích phân
y ( x) =

+

(3)

Nghiệm của phương trình (3) được xác định bằng dãy
x

 yn ( x)  y0    f ( s, yn1 ( s )  ds

x0

 y0 ( x)  y0

Phạm Thị Hoa K31E-SPToán

(4)


12

Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

Giả sử hàm f ( x, y ) liên tục trong R và trên đó thoả mãn điều kiện
Lipschitz theo biến thứ hai.
,
Ngoài ra, giả sử M  max f ( x, y ) , h  min(a,
R

b
).
M

Khi đó dãy các xấp xỉ (4) sẽ hội tụ tới nghiệm duy nhất của phương

trình (3) trên [

.

Tốc độ hội tụ được xác định bởi công thức
, n= 0, 1, 2, …
Với y ( x ) là nghiệm đúng của (3),

(5)

được tính theo (4).

b) Ví dụ.
Giải bài toán sau bằng phương pháp lặp đơn.
, với y (0)=0.

y '=

Giải
Hàm f ( x, y ) =

liên tục trên toàn mặt phẳng nên a, b có thể

chọn tuỳ ý.
Chọn

=0 thì xấp xỉ đầu tiên được xác định:
.

Tương tự:

Chúng ta ước lượng sai số của

=

.

=

.+

(x) bằng công thức (5)

Chọn a =1, b =0,5 thì ta có:

Phạm Thị Hoa K31E-SPToán

.


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

13

M  max f ( x, y )  max x 2  y 2  1,25.
G

G

L  max f ' ( x, y )  max 2 y  1.

G

G

h  min(a,

b
0,5
)  min(1,
)  0,4 .
M
1,25

Do vậy trên đoạn [0; 0,4 ] ta có
,
Từ đó

= 0,00133.

2. Phƣơng pháp xấp xỉ liên tiếp Picard.
a) Nội dung phương pháp.
*) Tìm nghiệm của bài toán Cauchy (1-2) là phương pháp xấp
xỉ liên tiếp. Gỉả sử các điều kiện tồn tại duy nhất nghiệm được thoả mãn. Việc
giải bài toán tương đương với việc tìm nghiệm của phương trình tích phân
sau:
y ( x) =

+

Vì ẩn y ( x ) tham gia vào dưới dấu tích phân nên gọi (1.1) là phương trình

tích phân.
*) Nội dung của phương pháp:Thay cho việc tìm nghiệm đúng
của (1.1) ta tìm nghiệm gần đúng thứ n theo công thức quy nạp:

(Để đơn giản) gọi nghiệm gần đúng đầu tiên là

Phạm Thị Hoa K31E-SPToán

,


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

14

(Theo lý thuyết phương trình vi phân): Nếu f ( x, y ) xác định liên tục
đồng thời thoả mãn điều kiện Lipschitz

trong G:

theo biến y thì quá trinh xấp xỉ liên tiếp (1.2) hội tụ đều về nghiệm duy nhất
của nghiệm (1.1) trong đoạn

h = min (a, ) trong đó:
M  max f ( x, y ) .
( x , y )G

Tốc độ hội tụ của phương pháp đặc trưng bởi hàm:
trong đó y ( x ) là nghiệm đúng của (1.1),
được tính từ công thức (1.2).

Người ta chứng minh được rằng
, n= 0, 1, 2,…
Từ đó suy ra khi n   thì
trên [

 0 và

(1.4).
(x) hội tụ đều về y ( x )

].
b) Ví dụ: bằng phương pháp Picard giải bài toán sau:

y '  y 2  x, y(0)  0 .
Bài làm.
Dễ thấy hàm số y '  y 2  x thoả mãn định lý về sự tồn tại và duy nhất
nghiệm trên toàn mặt phẳng.

y0 ( x) =0;

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng
x

15

x

1
y1 ( x)  y0 ( x) +  f (t , y0 (t ))dt   tdt  x 2
2
0
0
x

x

1 4
1 2 x5
y1 ( x)  y0 ( x) +  f (t , y1 (t ))dt   (t  t )dt  x 
4
2
20
0
0

y1 ( x)  y0 ( x) +
x


0

x

x11
1 4 1 7
1 10
1 2 x5 x8


+
.
f (t , y2 (t ))dt   (t  t  t 
t )dt  x 
160
4400
4
20
400
2
20
0

3. Phƣơng pháp chuỗi hàm nguyên.
a) Nội dung phương pháp.
Giả sử hàm f ( x, y ) ở vế phải của (1) là giải tích trong lân cận điểm
(

) nghĩa là f ( x, y ) khai triển được thành chuỗi nguyên:


f ( x, y ) =

a

k , m 0

km

( x  x0 )k ( y  y0 )m (1.6)

hội tụ trong lân cận đó. Với giả thiết đó bài toán Cauchy (1-2 ) là giải tích
trong lân cận đủ bé của điểm


y ( x) = 
k 0

nếu tìm được các giá trị

và có thể biểu diễn dưới dạng chuỗi Taylor:

k
y ( k ) ( x0 )
(
x

x
)
0
k!

(1.7)

( ) thì nghiệm gần đúng có thể lấy chẳng hạn

k
y ( k ) ( x0 )
=

( x  x0 ) .
k
!
k 0
N

Trong đó N là 1 số tự nhiên nào đó, nếu ta tính được các đại lượng
( ), (k=0, 1, ,…,N). Để tính

y ' ( )= f (

( ) ta xuất phát từ bài toán (1-2), ta có:
)

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

Lấy đạo hàm các cấp cả 2 vế (1) rồi thế x =
giá trị

vào ta được lần lượt các

( ), (k=0, 1, ,…
b) Ví dụ.

Bằng phương pháp chuỗi hàm nguyên giải bài toán sau:
x(

)=2; y (1)=0.
Giải

2
y'=   ey ;
x

y (1)=0; y ' (1)  1 ;
y '' =

2
 y 'e y ;
2
x

y '' (1)=1;
y '"  

4
 y " e y  ( y ') 2 e y ;
3
x

y '" (1)  2 ;

12
 y ''' e y  3 y ' y "e y  ( y ')3 e y ;
4
x
y (4) (1)  6;

y (4) 

y k ( x0 )
1
1
1
y ( x)  
( x  x0 ) k  ( x  1)  ( x  1) 2  ( x  1)3  ( x  1) 4 ;
k!
2
3
4
k 0
4

(4)

II. Các phƣơng pháp số.
1. Phƣơng pháp Euler.
a) Nội dung phương pháp.
Xét bài toán Cauchy:
y '  f ( x, y ), a  x  b (1)

y(a)  y0

(2)

Phạm Thị Hoa K31E-SPToán

16

Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

17

Chia đoạn  a, b thành các đoạn nhỏ bởi các điểm chia xi (i=0, 1, …, N)
sao cho: a  x0  x1  ...  xN  b .
Giả sử hàm f ( x, y ) có các đạo hàm riêng bậc m liên tục trên

R= a,b   y0 , Y , m  .
+) Với xi đủ gần x0 để tính được giá trị gần đúng của nghiệm bài toán
(1-2) tại điểm xi ta sử dụng công thức tính của phương pháp chuỗi hàm
nguyên trong nhóm các phương pháp giải tích đã trình bày ở mục trước như
sau:
j
y ( j ) ( x0 )
y ( x)  
( x  x0 ) .
j!
j 1
m

+) Với xi ở xa x0 để tính giá trị gần đúng y( xi )  yi ta sử dụng công
thức:

yi1  yi  hi f ( xi , yi ), i  0,1,..., N  1 (3)
trong đó hi  xi 1  xi . Trong phần bài tập ta sẽ chỉ xét phép phân hoạch đoạn

 a, b là phép phân hoạch đều tức là:

h  hi  xi1  xi , i  0,1,..., N  1. Khi đó

công thức (3) trở thành:

yi1  yi  hf ( xi , yi ), i  0,1,..., N  1 (3’).
Giả sử:

f  M1,

f
f
 M2,
 M 3 với x  a, b, y  y0 ,Y  thì sai số
x
y

của phương pháp Euler được đành giá theo công thức:

y ( xi )  yi 

M 4 h* M 3 ( xi  x0 )
e
2M 3

trong đó M 4  M 2  M1M 3 , h* = max hi
0i N-1

b) Ví dụ: Giải bài toán sau bằng phương pháp Euler

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

y'  y 

18

2x
, y (0)  1, x   0,1, h  0,2.
y

Giải
Áp dụng công thức:

 yi 1  yi  hf ( xi , yi )

 h  xi 1  xi , i  0,1,..., N  1
Ta có bảng sau:
Nghiệm
đúng

Δ
y=
0

0

1,0000

0

1,0000

0,2000

1,0000

1

0,2

1,2000

0,3333

0,8667

0,1733

1,1832

2

0,4

1,3733

0,5928

0,7805

0,1561

1,3416

3

0,6

1,5294

0,7846

0,7458

0,1492

1,4832

4

0,8

1,6786

0,9532

0,7254

0,1451

1,6124

5

1,0

1,8237

1,7320

2. Phƣơng pháp Euler-Cauchy.
a) Nội dung của phương pháp.
Xét bài toán Cauchy (1-2). Ta cũng phân hoạch  a, b bởi các điểm
N

chia {xi }i 0 ( a  x0  x1  ...  xN  b ). Phương pháp Euler đã trình bày ở trên
tính giá trị gần đúng yi  y( xi ) tuy đơn giản nhưng độ chính xác chưa cao. Để
khắc phục nhược điểm này chúng ta sử dụng phương pháp Euler-Cauchy. Sơ
đồ tính toán của phương pháp này như sau:
Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

 yi*1  yi  hf ( xi , yi )

h


*
 yi 1  yi   f ( xi , yi )  f ( xi 1 , yi 1 ) 
2

 i  0,1,..., N  1

19

(4)

Sai số được xác định bởi công thức:

Rr (h) 

r(l 1) ( )
(l  1)!

hl 1, 0    h.

Để làm tăng độ chính xác cho phương pháp Euler-Cauchy chúng ta áp
dụng phương pháp lặp đơn cho tính các giá trị yi như sau:

 yi(0)
1  yi  hf ( xi , yi )

h
 ( k 1)
( k 1)
 yi 1  yi   f ( xi , yi )  f ( xi 1 , yi 1 ) 

2

(k )
 yi 1  yi 1

(5)

b)Ví dụ: Tìm nghiệm gần đúng của bài toán sau bằng phương pháp
Euler-Cauchy:
y'  y 

2x
, y (0)  1, x   0,1, h  0,2.
y

Bài làm
Áp dụng công thức (4) ta có bảng sau:
i

xi

yi

h
fi
2

xi+1

yi+1*

h *
f i+1
2

h
y  ( f 
i
i
2

0

0

1

0,1

0,2

1,2

0,0867

0,1867

Phạm Thị Hoa K31E-SPToán

f

*

)

i 1

Nghiệm
đúng
1,0000


20

Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

1 0,2 1,1867 0,0850

0,4

1,3566 0,0767

0,1617

1,1832

2 0,4 1,3484 0,0755

0,6

1,4993 0,0699

0,1454

1,3416

3 0,6 1,4938 0,0690

0,8

1,6180 0,0651

0,1341

1,4832

4 0,8 1,6272 0,0645

1,0

1,7569 0,0618

0,1263

1,6125

5 1,0 1,7542

1,7325

Nhận xét: Nhược điểm của phương pháp Euler là ở chỗ trong

fi  yi1  yi chỉ tính đến đạo hàm tại điểm ( xi , yi ) ,

fi  hf ( xi , yi ) mà

không chú ý tới sự thay đổi của đạo hàm trong khoảng  xi , xi1  từ đó ta thấy
rằng nếu hàm f ( x, y ) thay đổi nhiều và không tuyến tính thì sai số mắc phải
sẽ lớn. Phương pháp áp dụng công thức Euler-Cauchy sẽ khắc phục nhược
điểm này.
3. Phƣơng pháp Runge―Kutta.
a) Nội dung phương pháp.
Xét bài toán Cauchy (1 - 2) ta chia đoạn  a, b thành các đoạn nhỏ bởi
các điểm chia xi (i=0, 1, …, N) sao cho: a  x0  x1  ...  xN  b trong đó

b  a 
.
xi  x0  ih , i  1,2,..., N , h  0, N  
 h 
Xuất phát từ giá trị ban đầu y0 tìm được giá trị gần đúng y1 tại

x1  x0  h theo công thức:
y( x0  h)  y1  y0  y0  y0   pr1K1 (h)  ...  prr Kr (h) 
trong đó:

i  x0  i h,i  y0  i1K1 (h)  ...  ii1Ki1 (h),i , ij , là các hằng số; 1 = 0;
Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

21

Ki (h)  hf (i ,i ) Các i , ij , pri được chọn sao cho với l đủ lớn, hàm

r (h)  y( x0  h)  y0  y0   pr1K1 (h)  ...  prr Kr (h) thoả mãn:

r (0)  r' (0)  ...  r(l ) (0)  0,r(l 1) (0)  0 .
Như vậy sai số mắc phải trong mỗi bước là:

Rr (h) 

r(l 1) ( )
(l  1)!

hl 1,0    h.

Từ: r( s ) (0)  0, s  0,1,..., l

 y0(i )  pr1K1(i ) (0)  pr 2 K 2(i ) (0)  ...  prr K r(i ) (0)
ta rút ra: 
i  1,2,..., l.

Sau đây chúng ta sẽ xét một vài trường hợp riêng thường dùng.
1. Trường hợp r=1: Chúng ta có phương pháp Euler.
2.Trường hợp r=2: Chúng ta có phương pháp Euler-Cauchy.
3. Trường hợp r=3.
Trong trường hợp này chúng ta có thể xây dựng được một trong các công
thức thông dụng sau:
1

y ( x0  h)  y1  y0  ( K1  4 K 2  K 3 ) .
6

Với :

h
K
K1  hf ( x0 , y0 ), K 2  hf ( x0  , y0  1 )
2
2

K3  hf ( x0  h, y0  K1  2K2 ) .
Sai số được xác định bởi công thức:

R3 (h) 

3(4) ( )
24

h4 .

4. Trường hợp r=4.

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

22

Một trong các công thức thông dụng nhất ứng với trường hợp này
là:

yi1  yi  yi ;

(6)

1
yi  ( K1(i )  2 K 2( i )  2 K 3( i )  K 4( i ) ) ;
6

(7)

K1(i )  hf ( xi , yi ) ;

(8)

h
K (i )
K 2(i )  hf ( xi  , yi  1 ) ;
2
2

(9)

h
K 2(i )
 hf ( xi  , yi 
);
2

2

(10)

K

(i )
3

K4(i )  hf ( xi  h, yi  K3(i ) ) ;

(11)

Sơ đồ tính cho phương pháp Runge-Kutta:
i

x

y

K  hf ( x, y )

y

x0

y0

K1(0)

K1(0)

K2(0)

2K 2(0)

K1(0)
2

h
x0 
2

y0 

h
x0 
2

K 2(0)
y0 
2

K3(0)

2K3(0)

x0  h

y0  K3(0)

K4(0)

K4(0)

0

y

1

x1

y1

Sai số được xác định bởi công thức:

Phạm Thị Hoa K31E-SPToán


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

R3 (h) 

 4(5) ( )
120

23

h5 .

5. Trường hợp r  6 : Chúng ta có phương pháp Runge-Kutta-Fehlberg.
Các công thức của phương pháp này như sau:

yn*1  yn 

16
6656
28561
9
2
K1 
K3 
K 4  K5  K6 ;
135
12825
56430
50
55

Trong đó:

K1  hf ( xi , yi ) ;
h
1
K 2  hf ( xi  , yi  K1 ) ;
4
4
K 3  hf ( xi 

3h
3K 9 K 2
, yi  1 
);
8
32
32

K 4  hf ( xi 

12h
1932 K1 7200 K 2 7296 K 3
, yi 


);
13
2197
2197
2197

K 5  hf ( xi  h, yi 

439 K1
3680 K 3 845K 4
 8K 2 

);
216
513

4104

h
8K
3544 K 3 1859 K 4 11K 5
K 6  hf ( xi  , yi  1  2 K 2 


).
2
27
2565
4104
40

Ở phần bài tập chúng ta chỉ xét với trường hợp r  4 (RK4)
b) Ví dụ: Giải bài toán sau bằng phương pháp Euler
y'  y 

2x
, y (0)  1, x   0,1, h  0,2.
y

Ta có bảng sau:
i

xi

yi

f ( xi , yi )

K  hf ( xi , yi )

yi

0

1

0

1

0,1

0,1

1,05

0,05

1,145238

0,114524

0,229048

Phạm Thị Hoa K31E-SPToán

Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

24

1,05

0,057262

1,159071

0,115907

0,231814

1,1

0,115907

1,310740

0,131074

0,131074
0,115323

1

1,1

0,115323

1,309678

0,130968

0,130968

1,15

0,180807

1,464447

0,146445

0,292889

1,15

0,188546

1,477805

0,147791

0,295581

1,20

0,263114

1,638523

0,163852

0,163852
0,147215

2

1,2

0,262538

1,637563

0,163756

0,163756

1,25

0,344416

1,801066

0,180107

0,360213

1,25

0,352591

1,814146

0,181415

0,362829

1,30

0,443953

1,983005

0,198301

0,198301
0,180805

3

1,3

0,443388

1,982135

0,198214

0,198214

1,35

0,524495

2,153696

0,215370

0,430739

1,35

0,551073

2,166404

0,216604

0,443281

1,4

0,660028

2,42897

0,234290

0,234290
0,216087

4

1,4

0,659475

2,342107

Phạm Thị Hoa K31E-SPToán

0,234211

0,234211


Khoá luận tốt nghiệp: Các phƣơng pháp giải bài toán Cauchy đối với ptvp thƣờng

25

1,45

0,776580

2,521146

0,252115

0,504229

1,45

0,785532

2,533493

0,253349

0,506700

1,50

0,912824

2,717099

0,271710

0,271711
0,252808

5

1,5

0,912283

2,716377

4. Phƣơng pháp Adams.
a) Nội dung phương pháp.
Năm 1855, nhà toán học người Anh Adams đề xuất một phương pháp đa
bước giải bài toán Cauchy theo yêu cầu của ông Bashforth, một chuyên gia kỹ
thật pháo binh Anh. Kết quả của Adams sau này bị quên lãng. Mãi đến đầu
thế kỷ XX, nhà toán học Na Uy Stermer trong khi tính quỹ đạo các hạt tích
điện rời xa mặt trời với vận tốc lớn, đã phát minh lại công thức Adams. Sau
này viện sỹ Krylov (Nga) đã có công hoàn thiện phương pháp Adams.
Tư tưởng chung của phương pháp Adams như sau:
Xét bài toán Cauchy (1-2), ta chia đoạn ta chia đoạn  a, b thành các đoạn
nhỏ bởi các điểm chia xi (i=0, 1, …, N) sao cho: a  x0  x1  ...  xN  b
trong đó h  xi1  xi .
Tích phân hai vế từ xn đến xn1 ta được:

yn1  yn 

x n 1



y '( x)dx .

xn
Đặt t 

x  xn
, ta có:

h
1

yn1  yn  h  y '( xn  th)dx . (12)
0

Phạm Thị Hoa K31E-SPToán