Luận văn một số phương pháp giải xấp xỉ phương trình vi tích phân tuyến tính fredholm

BỘ GIÁO D Ụ C VÀ ĐÀO TẠO
TR Ư Ờ NG ĐẠI HỌC s ư PH Ạ M H À NỘI 2

N G U Y Ễ N TH Ị T H U HÀ

M ỘT SỐ PH Ư Ơ N G P H Á P GIẢI X A P x ỉ
PH Ư Ơ N G T R ÌN H VI - TÍCH P H Â N
T U Y Ế N TÍN H FREDHOLM

LU Ậ N VĂN TH ẠC s ĩ TO Á N HỌC

H à N ội, 2016


BỘ GIÁO D Ụ C VÀ ĐÀO TẠO
T R Ư Ờ N G Đ Ạ I HỌC s ư P H Ạ M H À N Ộ I 2

N G U Y Ễ N TH Ị T H U H À

M ỘT SỐ PH Ư Ơ N G P H Á P GIẢI X A P x ỉ
PH Ư Ơ N G T R ÌN H VI - TÍCH P H Â N
T U Y Ế N TÍN H FREDHOLM

L U Ậ N V Ă N T H Ạ C s ĩ T O Á N HỌC
C huyên ngành: Toán giải tích
M ã số : 60 46 01 02

Người hướng dẫn khoa học
P G S .T S . K huất Văn N in h

H À N Ộ I, 2016

LỜI CẢM ƠN

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 PGS.TS. Khuất Văn Ninh
Tác giả xin bày tỏ lòng biết ơn sâu sắc nhất tới PGS.TS. Khuất Văn
Ninh, người đã định hướng chọn đề tài và tận tình hướng dẫn để tác giả
hoàn thành luận văn này.
Tác giả xin bày tỏ lòng biết ơn chân thành tới Phòng Sau đại học, các
thầy cô giáo 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ác giả trong suốt quá trình học tập và
hoàn thành luận văn tốt nghiệp.
Tác giả xin được gửi lời cảm ơn chân thành tới gia đình, bạn bè, người
thân đã luôn động viên, cổ vũ, tạo mọi điều kiện thuận lợi cho tác giả
trong quá trình học tập và hoàn thành luận văn.
Hà Nội, tháng 6 năm 2016
Tác giả

Nguyễn Thị Thu Hà


LỜI CAM Đ O A N

Tôi xin cam đoan, dưới sự hướng dẫn của PGS.TS. Khuất Văn Ninh,
luận văn Thạc sỹ chuyên ngành Toán Giải tích với đề tài “Một số phương
phấp giải xấp xỉ phương trình vi-tích phân tuyến tính Fredholm” do tôi
tự làm. Các kết quả và tài liệu trích dẫn được chỉ rõ nguồn gốc.
Trong quá trình nghiên cứu thực hiện luận văn, tôi đã kế thừa những
thành tựu của các nhà khoa học với sự trân trọng và biết ơn.

Hà Nội, tháng 6 năm 2016
Tác giả

Nguyễn Thị Thu Hà


M ục lục

M ỏ đầu

3

1

5

K IẾ N TH Ứ C C H U Ẩ N BỊ
1.1

Một số kiến thức về Giải tích hàm

5

1.1.1

5

1.1.2
1.1.3
1.2

1.3

2

3

Khống gian metric
Không gian định chuẩn
Không gian

7

và các tính c h ấ t ..........................

8

Môt, số kiến th ứ c v ề eũải t í c h ............................................................

9

1.2.1

C h u ỗi líív t h ừ a , ................

9

1.2.2

Tích phân phụ thuộc tham số và các tính chất . .

C ị a b -ị

10

11

Một số kiến thức về giải tích số
1.3.1

Phương pháp cầu phươnị

1.3.2

Sai phân và các tính chấ1

11

T

12

P H Ư Ờ N G P H Ấ P G ĨẲ Ĩ T ÍC H G IẢ I X Á P x ì P H Ư Ờ N G
T R ÌN H V Ĩ-T ÍC H P H Â N T U Y Ế N t í n h FR E D H O L M

14

2.1

Định lý về sự tồn tại nghiệm của phương trình

.............

14

2.2

Phương pháp tính toán trực tiếp

.......................................

15

2.3

Phương pháp phân tích A d o m ia n .......................................

20

2.4

Phương pháp c h u ỗ i................................................................

28

P H Ư Ờ N G P H Ấ P G ĨẲ Ĩ SỐ P H Ư Ờ N G T R ÌN H V Ĩ-T ÍC H
1


P H Ầ N T U Y Ế N T ÍN H FR ED H O L M

3.1

3.2

32

Phương pháp giải số phương trình vi-tích phân tuyến tính
F re d h o lm l................................................................................

32

Các ví dụ minh họa và ứng dụng Maple trong tính toán .

34

K ết luận

47

Tài liệu tham khảo

48

2


Mở đầu
1. Lí do chọn đ ề tà i
Toán học là một môn khoa học gắn liền với thực tiễn. Cùng với sự
phát triển của nội tại toán học và các ngành khoa học khác, toán học

chia thành toán lý thuyết và toán ứng dụng.
Trong lĩnh vực toán ứng dụng, thường gặp rất nhiều bài toán có liên
quan đến việc giải phương trình vi-tích phân. Phương trình vi-tích phân
tuyến tính Fredholm là loại phương trình xuất hiện trong toán học và
các ngành khoa học ứng dụng và từ lâu đã được các nhà toán học quan
tâm nghiên cứu.Việc tìm nghiệm chính xác của phương trình nói trên
gặp nhiều khó khăn. Vì vậy người ta nghiên cứu việc giải xấp xỉ phương
trình đó.
Phương trình vi-tích phân tuyến tính Fredholm có thể giải bằng các
phương pháp khác nhau. Trong đó, phương pháp giải tích cho nghiệm
dưới dạng biểu thức giải tích và phương pháp số cho nghiệm thu được
dưới dạng bảng số. Trong quá trình giải, ta có thể kết hợp sử dụng phần
mềm Maple trong tính toán.
Với mong muốn tìm hiểu và nghiên cứu sâu hơn về vấn đề này, dưới
sự hướng dẫn của thầy giáo PGS.TS. Khuất Văn Ninh tôi đã nghiên cứu
đề tầỉ“Một số phương phấp giải xấp xỉ phương trình vi-tích phẫn tuyến
tính Fredholm” để thực hiện luận văn của mình.

3


2. M ụ c đích n gh iên cứu
Luận văn sẽ nghiên cứu một số phương pháp giải phương trình vi-tích
phân tuyến tính Fredholm và ứng dụng Maple trong tính toán.

3. N h iệ m v ụ n gh iên cứu
Nghiên cứu một số phương pháp giải xấp xỉ phương trình vi-tích phân
tuyến tính Fredholm.

4. Đ ố i tư ợ n g và p h ạm v i n gh iên cứu

- Phương trình vi-tích phân tuyến tính Fredholm
- Các phương pháp giải xấp xỉ phương trình vi-tích phân tuyến tính
Fredholm

5. P h ư ơ n g p h áp n gh iên cứu
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 đ ón g góp
Hệ thống lại một số phương pháp giải phương trình vi-tích phân tuyến
tính Predholm và ứng dụng của phương pháp đó vào giải các phương
trình cụ thể. Áp dụng phần mềm Maple trong tính toán.

4


Chương 1
K IẾN THỨC C H U Ẩ N BỊ
1.1

M ộ t số kiến th ứ c về G iải tích hàm

Mục này nhắc lại một số kết quả về giải tích hàm, được trích dẫn chủ
yếu từ tài liệu [5].
1.1.1

K hông gian m etric

Cho X là một tập tùy ý
Đ ịnh nghĩa 1.1.1. Một metric trong X là một ánh xạ

d X

X

X —ỳ K,

thỏa mãn các điều kiện sau đây
(i) d ( x ,y ) ^ 0,Vx,y e X;
(ii) d(x,y) = 0

X

= y\

(hi) d(x,y) = d ( y ,x ) ,V x ,y e X ;
(iv) d(x, y) < d(x, z) + d(z, y) , Vs , í / GX;
Một không gian metric là một tập hợp cùng với một metric trong
tập hợp ấy. Các phần tử của một không gian metric được gọi là
điểm của không gian ấy. số d(x, y) được gọi là khoảng cách giữa
các điểm

X

và

y.

5

Đ ịnh nghĩa 1.1.2. Một dãy các điểm (xn) , n = 1,2,... trong không
gian metric X được gọi là hội tụ đến điểm a € X nếu
lim d(a, x n) = 0.

ra—
>oo

Khi đó ta kí hiệu lim x n = a hoặc x n !-»■ a khi Ĩ Ỉ 4 oo
n—
>00
Đ ịn h nghĩa 1.1.3. Dãy điểm

được gọi là dãy cơ bản trong không

gian metric X nếu với mọi £ > 0 cho trước, tồn tại một số n 0 sao cho
với mọi n > n0 và m > n 0 ta đều có
d{xn, x m) < £.
Nói cách khác ta có
lim d(xn, x m) = 0.

n,m—>00

Dễ thấy mọi dãy điểm hội tụ trong không gian metric đều là dãy cơ
bản.
Đ ịnh nghĩa 1.1.4. Một không gian metric X được gọi là đầy đủ nếu
mọi dãy cơ bản trong X đều hội tụ tới một phần tử trong X .
Đ ịnh nghĩa 1.1.5. Cho X , Y là hai không gian metric tùy ý. Ánh xạ
f x —>• Y được gọi là một ánh xạ co nếu tồn tại một số a với 0 < a < 1
sao cho với mọi x,x ' € X ta đều có
d {f{x ),f{x ')) < a d ( x , x ’),

và a được gọi là hệ số co của / .
Hiển nhiên một ánh xạ co là ánh xạ liên tục đều.
Đ ịnh lý 1.1.1. (Nguyên lý ánh xạ co) Giả sử X là một metric đầy đủ
và f : X —>• X là một ánh xạ co của X vào chính nó. Khi đó tồn tại một
6


và chỉ một điểm

X*

e X sao cho f(x*) =

X*.

Hơn nữa

X*

ỉà giới hạn của

dãy (x n) được xây dựng như sau
x ữ tùy ý thuộc X ,

Xn + 1

— f ( x n) , n > 0 và tốc độ hội tụ được đánh giá

theo cônq thức
d (x n, X*) < ------ d (x1: x 0) ,

1 —a
trong đó a là hệ số co của f .
1.1.2

K hông gian định chuẩn

Cho X là một không gian vectơ trên trường p (p = K hoặc

c)

Đ ịnh nghĩa 1.1.6. Một chuẩn, kí hiệu ||.|| trong X là một ánh xạ từ X
vào R thỏa mãn các điều kiện
(i) ||x|| > 0 với mọi l ẽ X ;
(ii) ||x|| = 0 khi và chỉ khi

X

= 0 [9 là kí hiệu phần tử không);

(iii) \\Ằx\\ = ỊAỊ ỊỊa^ll với mọi số À € p và với mọi X € X;
(iv) IIX + y II < ||x|| + \\yII với mọi x , y E X .
Số ||x|| được gọi là chuẩn (hay độ dài) của vectơ

X

e X.

Đ ịn h nghĩa 1.1.7. Một không gian vectơ X cùng với một chuẩn xác
định trong không gian ấy gọi là một không gian định chuẩn (thực hoặc
phức, tùy theo p là thực hoặc phức).

Đ ịn h lý 1.1.2. Giả

sử

X

ỉà

một không gian định chuẩn. Với mọi

đặt
d( x, y) = \\x - y II .
Khi đó d là một metric trên X .

7

X

GX


Đ ịnh nghĩa 1.1.8. Dãy (x n) trong không gian định chuẩn X được gọi
là hội tụ đến x 0 e X nếu lim II
n —>00

—íEoll = 0. Khi đó ta kí hiệu

lim x n — x 0 hay x n —>x 0(n —> oo).

n —>OQ

Đ ịnh nghĩa 1.1.9. Dãy (:rn) trong không gian định chuẩn X được gọi
là một dãy cơ bản nếu
lim I I —x m II = 0.

n .m —> oc

Đ ịnh nghĩa 1.1.10. Không gian định chuẩn X gọi là không gian Banach, nếu mọi dãy cơ bản trong X đều hội tụ.
1.1.3

K hông gian Cịab] và các tín h chất

Đ ịn h nghĩa 1.1.11. C[a 6] là tập tấ t cả các hàm số giá trị thực xác định
và liên tục trên đoạn [a, b], —oo < a < b < +oo.
Các tín h chất
(i) Không gian Cịa 6] là không gian metric.
Vz, y € CịaM, d (x, y) = max \x (t ) - y { t ) I.
a(ii) Không gian C[a 6] là không gian định chuẩn.
||x|| = max \x ( t ) |.
ã<ỉ(iii) Không gian Cịa 6] là không gian Banach.
(iv) Tập tấ t cả các đa thức với hệ số hữu tỷ trù m ật trong Cịab-ị. Cho
nên Cịa 6] là không gian tách được.
Đ ịnh nghĩa 1.1.12. Không gian C ^ 6Ị gồm tấ t cả các hàm X (t) xác
định trên đoạn [ữ, 6] và có đạo hàm liên tục đến cấp n, với chuẩn được
xác định bởi
INI = max { |z ( í) |, |íc,( t ) |,..., \xn(t)\}
a8

1.2

M ộ t số kiến th ứ c về giải tích

1.2.1

Chuỗi lũy thừa
+ 00

Đ ịnh nghĩa 1.2.1. Chuỗi lũy thừa là một hàm dạng Xì an{x — x 0)n
n=0
trong đó x ữ, a0, dị, a2, ... là những số thực.
Điểm x ữ được gọi là tâm của chuỗi lũy thừa. Để ý rằng chuỗi lũy thừa
luôn luôn hội tụ tại điểm

X

= x ữ.
+ 00

Nếu đặt y =

x ữ thì ta có thể đưa chuỗi lũy thừa về dạng Xì anyn,
n=0
chuỗi có tâm tại y — 0
X —

Các tín h chất của tổn g chuỗi lũy thừa

+ 00

Đ ịnh lý 1.2.1. Giả sử chuỗi lũy thừa

anx n có bán kính hội tụ R > 0,
n=0
khi đó tổng S ( x ) của nó là một hàm liên tục trong khoảng hội tụ ( - R , R).
+ 00

Đ ịnh lý 1.2.2. Giả sử chuỗi lũy thừa
Khi đó tổng

s

anx n có bán kính hội tụ R > 0.
71=0

của nó là hàm khả tích trên mọi đoạn [a, b] nằm trong

khoảng hội tụ ( - R , R) và

/

b

pb
an / x ndx

+°°

S{x)dx — ^
n= 0

a

Dặc biệt nếux € (—R , R ) thì Jq S(t)dt = X) “ n+1
n=0
+ 00

Đ ịn h lý 1.2.3. Giả sử chuỗi lũy thừa
vaV
S(x ) —

n=0

anx n có bán kính hội tụ R > 0

anx7lì x ẽ ( —R, R)n= 0

Khi đó
+ 00

+ Chuỗi lũy thừa

nanx n~l nhận được bằng cách đạo hàm từng số
71=0

hạng của chuỗi lũy thừa đã cho, cũng có bán kính hội tụ là R.
9

+ Tổng

s là hàm khả vi trong khoảng hội tụ ( - R , R )

và

nanx n 1
s ' ( x ) = ^ 2 nơnx
n=0
1.2.2

Tích phân phụ thu ộc tham số và các tín h chất

Đ ịnh nghĩa 1.2.2. Giả sử f ( x , y ) là một hàm số xác định với X thuộc
đoạn [a, 6] và y thuộc một tập hợp số thực Y nào đó, sao cho với mỗi y
cố định thuộc Y hàm f ( x , y ) khả tích trong đoạn [a,h].
Đặt I(y) = Ị bf( x ,y ) d x .
Khi đó I(y) là một hàm số xác định trên tập Y và được gọi là tích
phân phụ thuộc tham số của hàm f ( x , y ) trong đoạn [a, b].
Các tín h chất của tích phân xác định phụ thu ộc tham số
Giả sử f ( x , y) là hàm số xác định trong hình chữ nhật
D = [a, 6; c, dị = [a, 6] X [c,d] = { ( x ,y ) ,a < X < b;c < y < d}
Đ ịnh lý 1.2.4. Nếu hàm f ( x , y ) xác định và liên tục trong hình chữ
nhật D thì tích phân phụ thuộc tham số I ( y ) = f b f ( x ì y)dx là một hàm
liên tục trong đoạn [c, d].
Đ ịnh lý 1.2.5. Giả sử f ( x , y ) là hàm số xác định trong hình chữ nhật
D liên tục theo X € [a,b] với mỗi y cố định thuộc đoạn [c,d]. Hơn nữa
f ( x , y ) có đạo hàm riêng I~(x,y) là một hàm liên tục trong hình chữ
nhật D. Khi đó tích phân phụ thuộc tham số

b
a
là một hàm khả vi và

10


Đ ịnh lý 1.2.6. Nếu f ( x , y ) là hàm liên tục trong hình chữ nhật D =
[a, 6] X [c, d] thì ta có công thức
j \ y ) d y = Ị* {
Hay là

/

d

pb
dy

u

1.3
1.3.1

y)đx} dy =

y

rb
f (x, y)dx =

Ja

‘ f ( x , y ) d y } dx.

pd
dx

Ja

f ( x :y)dy.
Jc

M ộ t số kiến th ứ c về giải tích số
Phương pháp cầu phương

Cho hàm / xác định trên đoạn [a,b], f là hàm số liên tục trên đoạn
[ữ, 6] do đó / khả tích trên đoạn [a, 6].
Ta chia đoạn [a, b] thành n phần
a = x 0 < Xị < x 2 < ... < x n = b
Công thức sau được gọi là công thức cầu phương

/

b

n

(a;) dx —

E ẢkV 0xk) + R n

( 1. 1)

k= 1

trong đó, Ak và Xỵ- tương ứng là hệ số và nút của công thức cầu phương,
R n (tp)~ phần dư của công thức cầu phương.
Nếu như chọn công thức hình thang, thì chúng ta có
Aị

An

- h , Aỵ

h, k

2,..., n

1ị

b —a
Xk = a + (k — 1) h, k = 1,..., n , h = -----n —1
Rn {Kv)

{b-aý
12 (n — lj
11

[ ô2 ( ^ ) '

dy2

y=z,a

Nếu xét công thức Sympson thì sau khi đặt n = 2m + 1, ta sẽ có
A\ = A 2m+1 = —, A 2 = Aị = ... = A 2m = —h,
A 3 — A 5 — ... — A 2m_\ — —h,
b —a
Xỵ — a + (k — 1) h : k — 1 , 2 m + 1, h — ------- ;
2m
1 (b -a f \d U K J
n(
90 (2m)4 L
J„=í,a<{<„
Nếu một quy tắc nào đó được chọn thì các đại lượng A k,X]¡;, Rỵ có thể
được viết một cách tương tự.
1.3.2

Sai phân và các tín h chất

Đ ịnh nghĩa 1.3.1. Giả sử y = f (x) là hàm số xác định trên tập X , h
là hằng số lớn hơn 0. Biểu th ứ c A / (x) = f (x + h) — f (x) được gọi là
sai phân cấp 1 của / (x) tại điểm X. Biểu thức
A 7 = A [A/ (*)]
= [ f ( x + 2h) - f (x + h)] - [f {x + h) - f (x)]
= A f (x + h) - A f (x)
được gọi là sai phân cấp 2 của / (x) tại

X.

Tương tự, ta có A kf — A [Afe-1/] được gọi là sai phân cấp k của /
tại X.
Các tín h chất của sai phân
(i) A* [ ĩ ± g ] = A

± A kg.

(ii) A* [A/ (x)] = AA‘ [/ (x )].
(iii) A n \pn (a:)] = const, A m \pn (x)] = 0, khi m > n, trong đó pn (x) là
đa thức cấp n của X.
12


(iv) / (x + nh) = £ Q A 7 ( x ) .
i= 0

(v) A”/ ( i ) = Ẻ (—l ) ‘C i/ [z + (n - t)].
¿=0

13


Chương 2
PH Ư Ơ N G P H Á P GIẢI TÍCH
GIẢI X Ấ P XỈ PH Ư Ơ N G T R ÌN H
VI-TÍCH P H Â N T U Y Ế N TÍNH
FREDHOLM
Xét phương trình vi-tích phân tuyến tính Fredholm được cho bởi công
thức

(X) = / (x) + í K (x, t) u (t ) dt, u ^ (a) =
Ja

0 < k < n — 1 (2.1)

ở đó vSn\ x ) là đạo hàm bậc n của hàm u ( x ) với biến số là X và bk là cho
trước. Giả sử K ( x , t ) liên tục trên D = [a, b] X [a, 6], M = max \K(x, t)\
(.x , t ) e D

2.1

Đ ịn h lý v ề sự tồ n tạ i n gh iệm củ a phư ơng trìn h

Đ ịnh lý 2.1.1. Giả sử hàm K ( x , t ) liên tục trên tập D và Mịb —à) < 1.
Khi đó phương trình (2.1) có nghiệm duy nhất u = u(a^) thỏa mãn điều
kiện ban đầu u

(a) =

0 < k < n — 1.

Nghiệm này xác định trên đoạn [a, 6].
Chứng minh
Đặt
A{x, u) = f{ x ) + / K{x,t)u{t)dt
Ja
14


Khi đó theo tính chất của tích phân phụ thuộc tham số A(x, u ) là hàm

số liên tục trên đoạn [a, 6].
Ta lại có
\A(x,u) — A(x,v)\ = ị ị K (x ,t)u (t)d t — ị K ( x ,t)v (t)d t
ị**a
Ja
— ị í K ( x , t ) [u(t) —v(t)] dt
IJ a
< ị \K(x, t)\ \u(t) —v(t)\ dt
Ja
< M í \u(t) —v(t)\dt
Ja
< M(b — a)\\u — v\\ = q IIu — 1>|| ,
với (q = M(b — a))
Khi đó theo định lý về sự tồn tại nghiệm của bài toán Côsi suy ra
phương trình (2.1) có nghiệm duy nhất u = u(a^) thỏa mãn điều kiện
ban đầu.

2.2

P h ư ơ n g p h áp tín h to á n trự c tiế p

Xét phương trình vi-tích phân Predholm được cho bởi công thức(|2.ip
Thay thế K ( x , t ) = g(x).h(t) vào công thức (2.1 )ta được
u ^n\ x ) = f{ x ) + g{x).

Ị h(t).u(t)dt,u^k\ a ) = 6jfc,0 < k < n — 1 (2.2)

Ja

Chúng ta có thể dễ dàng thấy rằng tích phân xác định trong phương

trình vi-tích phân (2.2) liên quan đến một tích phân hoàn toàn phụ thuộc
vào biến t. Điều này có nghĩa là tích phân xác định ở vế phải của (2.2)
là bằng một hằng số a. Nói cách khác, chúng ta đặt
a = í h(t)u(t)dt
**a
15

(2.3)


Do đó, phương trình (2.2) trở thành
u(n\ x ) = f ( x ) +

)

(2.4)

Tích phân cả hai vế của (2.4) n lần từ a tới X và sử dụng các điều kiện
ban đầu, chúng ta có thể tìm được một biểu thức của u ( x ) liên quan đến
hằng số a và biến X. Nghĩa là chúng ta có thể viết
u (x) = V(x, a )

(2.5)

Thay thế (2.5) vào vế phải của (2.3). Đánh giá tích phân và nghiệm
của phương trình, chúng ta xác định được hằng số a. Với cách làm này ta
sẽ thu được nghiệm chính xác u(x) bằng cách thay giá trị a vào phương
trình (2.5).
Phương pháp sẽ được minh họa bằng các ví dụ sau
V í dụ 2.2.1. Tìm nghiệm của phương trình vi- tích phân sau

u (x) = —1 + cosx + / tu(t)dt,u(0) = 0.
■'o
Bài giải
Đặt
a=

(2.6)

tu(t)dt

Khi đó phương trình đã cho trở thành
u (x)

= —1 +

cosx

Oi.

+

Tích phân 2 vế từ 0 đến X và sử dụng điều kiện ban đầu ta có
PX

u(x) = / (—1 + cost + at)dt
/

.

Oi

z
a

QV ,

= (—t + sỉnt + ^-t2) lo

= —X + sinx H----- X
2

16

2


Thay vào (2.6) ta được

<
^a=
w ư

T ' 16

3 (1 2 8 —7T4)

Vậy nghiệm chính xác u{x) = —X + sinx + 3( li*128)x2•
V í d ụ 2.2.2. Tìm nghiệm của phương trình vi-tích phân sau
„
fi
,

u (x) = —1 —sỉnx + / tu(t)dt, u{0) = 0, u (0) = 1.

Bài giải
Đặt
a —

tu(t)dt

(2.7)

Phương trình được viết dưới dạng
u (x) = —1 — sinx + a.
Tích phân hai vế từ 0 tới

X

hai lần và sử dụng điều kiện ban đầu ta

được
px

u (x) = / (—1 —sint + à)dt
•'O
= (—t + cost + at) lỗ
= x(a — 1) + cosx

17


Khi đó

rx

u (X) = / [t(a — 1) + cost] dt
do
t2
— {a — 11)) + sỉntị
sintịl
= —(a
= — (a — 1 ) + sinx
2

Thay vào (2.7) ta được
r+2
a = ự t —{a — 1) + sint dt
2
do
a —1
-t3 + tsint dt
(a ~ 11,417

8

u

[

tsintdt

(a - 1)
Õ2 4 + 1

8
'

& a= l
Vậy nghiệm chính xác là Ií(:r) = sinx.
V í d ụ 2.2.3. Tìm nghiệm của phương trình vi tích phân sau
u (à?) = e—2—x + ea:(3+ x)+ í (x—t)u(t)dt,u(0) = 0,w (0) = 1,u (0) = 2.
“'0
Bài giải
Ta có
u (x) = e —2 —X + ex(3 + x) + X ị u(t)dt — ị tu(t)dt
do
dữ
Đặt
= í u(t)dt
do

( 2.8)

OL2 = ị tu(t)dt
do

(2.9)

18


Khi đó phương trình đã cho trở thành
u (X) = e —2 —X + ex (3 + x) + XOLI —«2Tích phân cả hai vế của phương trình trên từ 0 tới X ba lần và sử dụng
điều kiện ban đầu ta được

/£»4

£.3

u ( x ) = ^—{a.1 - 1) + -—{e - 2 - a 2) + 2x2 + X + xex
Zi.E

o

Thay u ( x ) vào , ta được
r 1 t4 ,
t3 ,
Qíi — Ị ——(ơi — 1) + —(e — 2 — QÍ2) T 2t2 T t T tetdt
J0 z 4
0
= ^Ĩ2Õ^a i _ ^ + 24^e ~~ 2 _ “ 2) + g í3 + ~2 + Íeí - ei\ lo
1
_
ĩ
_
7
= ---- (ơi — 1) T — (e —2 —QÍ2) T —
120v
J 2ẩy
} 6

Oi2

= ị 7rr{a i — 1) + ^-(e —2 —a 2) + 2í3 + t 2 T t 2ê à t
J0 ¿4

0
= [jỊ Ị Í a i - 1) +

- 2 - a 2) + - i 4 + — + í V - 2íe‘ + 2e‘| lị

“ Ĩ 4 4 ( a i _ 1 ) + ằ ( e _ 2 _ “ 2:) + 6 + e
Giải hệ phương trình ta thu được
ca1
0-2

98711 + 7440e
88511
65738 + 88561e
88511

Vậy nghiệm chính xác

“(*) =

t4 ,10260 T 7440e x t3 , 242760 T 50e
“) 4" 2t2 T Í T íe4.
88511
^ 6 ^ 88511

19


2.3

P h ư ơ n g p h áp p h ân tích A d om ian

Trong phương pháp này chúng ta thực hiện bằng cách chuyển đổi
một phương trình vi-tích phân Fredholm về một phương trình tích phân
Fredholm. Sau đó ta giải phương trình tích phân đó bằng phương pháp
phân tích Adomian.
Không làm mất tính tổng quát chúng ta có thể xét một phương trình
vi-tích phân Fredholm bậc hai được cho bởi
u (:r) = f ( x ) + / K ( x ,t)u ( t) d t,u (a ) = d0,u (a) = di
•'a
Tích phân cả hai vế của (2.10) từ a đến X hai lần được
u{x) = a0 + diX + L~1(f(x)) + L~l { í K ( x ,t )u ( t) d t)
''a

(2-10)

(2.11)

ở đây điều kiện ban đầu u(a) = a0 và ú(a) = «1 được sử dụng, L~l là
một toán tử tích phân hai lần. Sau đó ta sử dụng chuỗi khai triển
00

u{x) = ^ 2 un(x).

(2.12)

n= 0

Thay vào cả hai vế của phương trình (2.11) thu được
oo

«5

oc

^ 2 un{x) = a0 + ũịX + L~1Ự(x)) + L~x{ j K ( x , t ) ^ 2 u n(t)dt) (2.13)
n=0
°
n=0
hoặc tương đương
u0(x) + Ui(x) + u2(x) + u3(x) + ... = a0 + a,i(x) + L~1(f(x))+
+ L~1( ị K ( x , t ) u 0(t)dt) + L~1( ị K ( x , t)ui{t)dt) + ...
Ja
Ja
Do đó để xác định các thành phần u0(x), Ui(x), u2(x)... của nghiệm
u(x). Chúng ta đặt một quan hệ lặp
u0(x) = a0 + dịX + L~1{f{x))

Ị

uk+i{x) = L~l { K ( x , t ) u k(t)dt), k > 0
Ja
20


ở đó th à n h ph ần u0(x ) bằn g tấ t cả các số hạn g kh ông bao gồm phần

bên trong dấu tích phân của (2.13). Ta phải xác định các thành phần
Uị(x)J > 0, nghiệm u(x) của (2.10) có được sau khi thu được Ui(x) . Sử
dụng (2.12) ta thu được chuỗi hội tụ tới nghiệm chính xác. Để sử dụng
tốt phương pháp phân tích Adomian chúng ta kết hợp với phương pháp

phân tích cải biên.
Phương pháp phân tích A dom ian cải biên
Theo phương pháp phân tích Adomian nghiệm được tìm dưới dạng
chuỗi vô hạn u(x)
oc
u(x) =

un(x)
n=0

(2.14)

Thay thế vào cả hai vế của phương trình tích phân Fredholm
u(x) = f ( x ) + X í K (x,t)u(t)dt
Ja
Ta có công thức truy toán
u0(x) = f ( x ),
Uk+i(x) = À /f b K (x, t)uk(t)dt, k > 0.
Ja
Nhìn vào công thức trên các thành phần un(x), n > 0 sẽ dễ dàng tìm
được. Nếu hàm số f ( x ) là tổng của hai hàm thành phần, là fi(x)vầ Ỉ 2{x).
Nói cách khác chúng ta có thể đặt
f ( x ) = fị(x) + f 2{x).

(2.15)

Khi đó ta sử dụng phương pháp đồng nhất hóa thành phần u0{x) bằng
một bộ phận của f ( x ) , là fi(x) hoặc f 2{x). Và sử dụng phép truy toán

21