Nội suy newton và bài toán biên hỗn hợp thứ nhất của phương trình vi phân

ĐẠI HỌC THÁI NGUYÊN
TRƯỜNG ĐẠI HỌC KHOA HỌC
- - - - - - - - - - - - - - - - - -
Trần Thu Hiền
NỘI SUY NEWTON
VÀ BÀI TOÁN BIÊN HỖN HỢP THỨ NHẤT
CỦA PHƯƠNG TRÌNH VI PHÂN
LUẬN VĂN THẠC SỸ TOÁN HỌC
Chuyên ngành: TOÁN ỨNG DỤNG
Mã số: 60.46.01.12
Người hướng dẫn khoa học
GS. TSKH. NGUYỄN VĂN MẬU
THÁI NGUYÊN - NĂM 2014
Mục lục
Mở đầu 3
1 Một số kiến thức chuẩn bị 6
1.1 Công thức nội suy Taylor . . . . . . . . . . . . . . . . . . . . 8
1.2 Công thức nội suy Newton . . . . . . . . . . . . . . . . . . . 9
1.3 Khai triển Taylor - Gontcharov . . . . . . . . . . . . . . . . . 12
1.3.1 Khai triển Taylor-Goncharov với phần dư dạng Lagrange 18
1.3.2 Khai triển Taylor-Goncharov với phần dư dạng Cauchy 19
2 Bài toán biên hỗn hợp thứ nhất của phương trình vi phân 21
2.1 Bài toán Cauchy . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.1.1 Bài toán Cauchy trừu tượng . . . . . . . . . . . . . . 23
2.1.2 Bài toán Cauchy của phương trình vi phân . . . . . . 29
2.2 Bài toán biên hỗn hợp thứ nhất . . . . . . . . . . . . . . . . 32
2.2.1 Bài toán biên hỗn hợp thứ nhất trừu tượng . . . . . . 32
2.2.2 Bài toán biên hỗn hợp thứ nhất của phương trình vi
phân . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
3 Ví dụ áp dụng 44
Kết luận 49

Tài liệu tham khảo 51
1
Một số kí hiệu dùng trong luận văn
Trong luận văn có sử dụng một số kí hiệu sau đây:
+ L(X) : tập tất cả các toán tử tuyến tính tác động trong X
+ domA: miền xác định của tập A
+ L
0
(X) = A ∈ L(X) : domA = X
+ R(X): tập tất cả các toán tử khả nghịch phải trong domA
+ R
D
: tập tất cả các toán tử khả nghịch phải của D ∈ R(X).
2
Mở đầu
Lý thuyết nội suy và các vấn đề liên quan là một lĩnh vực tuy không mới
nhưng luôn là chuyên đề được nhiều người quan tâm nghiên cứu. Nó không
chỉ quan trọng đối với toán học, đặc biệt là đối với chuyên ngành Đại số và
Giải tích mà nó còn làm công cụ để giải quyết nhiều vấn đề của thực tiễn
cuộc sống.
Trong toán học, có rất nhiều cách để xác định được một hàm số nhưng
để xác định được chính xác một hàm số mà chỉ nhờ vào một số giá trị rời
rạc của hàm số đó và đạo hàm của nó thì phương pháp nội suy là hữu hiệu
nhất mà trong đó phải kể đến bài toán nội suy Taylor và nội suy Newton:
Bài toán nội suy Taylor. Hãy xác định đa thức N(x) có bậc không quá
n (deg N(x) ≤ n) và thỏa mãn các điều kiện
N
(i)
(x
0

) = a
i
, i = 0, 1, . . . , n.
Bài toán nội suy Newton. Hãy xác định đa thức N(x) có bậc không quá
n (deg N(x) ≤ n) và thỏa mãn các điều kiện
N
(i)
(x
i
) = a
i
, i = 0, 1, . . . , n.
Từ bài toán nội suy Taylor và nội suy Newton, ta có thể phát triển để
nghiên cứu về phương trình vi phân - một vấn đề được đề cập rất nhiều trong
chương trình toán phổ thông cũng như chương trình toán ở các trường cao
đẳng, đại học với bài toán Cauchy (bài toán ban đầu) và biên hỗn hợp thứ
nhất như sau
Bài toán Cauchy. Tìm tất cả các nghiệm của phương trình
n

i=0
a
i
(t)x
(n−i)
(t) = y(t),
thỏa mãn điều kiện ban đầu
x
(i)
(t

0
) = b
i
, i = 0, 1, . . . , n.
3
Bài toán biên hỗn hợp thứ nhất của phương trình vi phân. Tìm tất cả
các nghiệm của phương trình
n

i=0
a
i
(t)x
(n−i)
(t) = y(t)
thỏa mãn điều kiện sau
x
(j)
(t
j
) = b
j
, j = 0, 1, . . . , n.
Trong thực tiễn, phương trình vi phân có rất nhiều ứng dụng. Chúng ta
có thể nghiên cứu về sự gia tăng dân số, về sự phân rã phóng xạ, về sự nóng
lên hoặc nguội đi của vật thể
Với những lí do trên đây nên tác giả chọn: "Nội suy Newton và bài toán
biên hỗn hợp thứ nhất của phương trình vi phân" làm đề tài nghiên cứu luận
văn tốt nghiệp.
Khi thực hiện luận văn này tác giả thấy được việc nghiên cứu là hữu ích

đối với bản thân, trước hết là hình thành và rèn luyện cho người viết kĩ năng
cũng như kinh nghiệm nghiên cứu một vấn đề khoa học. Đồng thời người
viết luôn được trau dồi và cập nhật những kiến thức mới. Ngoài ra luận văn
còn là tài liệu tham khảo, nghiên cứu cho các học viên cao học, giảng viên
cũng như sinh viên các trường trong cả nước.
Ngoài mục lục, lời nói đầu, kết luận và tài liệu tham khảo, luận văn bao
gồm 3 chương như sau
Chương 1. Một số kiến thức chuẩn bị
Chương này ngoài phần đầu nhắc lại khái niệm và các tính chất liên quan
đến đa thức đại số là ba phần: phần thứ nhất trình bày về công thức nội
suy Taylor, phần thứ hai trình bày về bài toán và công thức nội suy Newton,
phần thứ ba đưa ra khai triển Taylor-Gontcharov và chứng minh các định lý
làm nền tảng kiến thức cho các chương sau.
Chương 2. Bài toán biên hỗn hợp thứ nhất của phương trình vi phân
Chương này là phần chính của luận văn, tác giả trình bày bài toán Cauchy
và bài toán biên hỗn hợp thứ nhất trừu tượng. Sau đó áp dụng khảo sát bài
toán Cauchy và bài toán biên hỗn hợp thứ nhất của phương trình vi phân.
Chương 3. Ví dụ áp dụng
Trong chương này, tác giả đưa ra các ví dụ điển hình để làm minh họa cụ thể
cho bài toán Cauchy và bài toán biên hỗn hợp thứ nhất của phương trình vi
phân đã đề cập trong Chương 2.
4
Tác giả đã tiến hành nghiên cứu, phân tích khá kỹ các tài liệu có liên
quan rồi từ đó tập trung đánh giá, tổng hợp lại để làm nền tảng kiến thức
cho luận văn. Nhưng do lần đầu tiên làm quen với việc nghiên cứu khoa học,
mặc dù đã có rất nhiều cố gắng, song do hạn chế nhiều về mặt thời gian, kiến
thức và kinh nghiệm nghiên cứu nên luận văn không thể tránh khỏi những
thiếu sót. Tác giả rất mong nhận được sự góp ý, chỉ bảo, bổ sung của các
thầy cô giáo cũng như sự quan tâm của các bạn bè đồng nghiệp để luận văn
ngày càng hoàn thiện hơn.

Hoàn thành được luận văn này, tác giả xin được chân thành cảm ơn thầy
giáo GS. TSKH. NGUYỄN VĂN MẬU là người trực tiếp hướng dẫn và
đã tận tình giúp đỡ tác giả trong suốt thời gian thực hiện luận văn. Tác giả
xin chân thành cảm ơn Ban chủ nhiệm khoa Toán, phòng Đào tạo trường
Đại học Khoa học - Đại học Thái Nguyên cùng các thầy, cô giáo đã tham
gia giảng dạy khóa học.
Xin cảm ơn các anh chị đồng nghiệp, các anh chị đồng môn là những
người luôn ủng hộ, động viên và tạo điều kiện cho tôi trong thời gian hoàn
thành luận văn.
5
Chương 1
Một số kiến thức chuẩn bị
Ở phần này, ta nhắc lại các kiến thức về đa thức đại số.
*) Đa thức đại số và các tính chất liên quan.
Định nghĩa 1.1 (xem [1]). Cho vành A là một vành giao hoán có đơn vị.
Ta gọi đa thức bậc n biến x là một biểu thức có dạng
P
n
(x) = a
n
x
n
+ a
n−1
x
n−1
+ · · · + a
1
x + a
0

(a
n
= 0)
trong đó các a
i
∈ A được gọi là hệ số, a
n
là hệ số bậc cao nhất và a
0
là hệ
số tự do của đa thức.
Nếu a
i
= 0, ∀i = 1, 2, . . . , n và a
0
= 0 thì ta có bậc của đa thức là 0.
Nếu a
i
= 0, ∀i = 0, 1, . . . , n thì ta coi bậc của đa thức là −∞ và gọi là
đa thức không (nói chung thì người ta không định nghĩa bậc của đa thức
không). Tập hợp tất cả các đa thức với hệ số lấy trong vành A được kí hiệu
là A [x].
Khi A = K với K là một trường thì vành K [x] là một vành giao hoán có
đơn vị. Ta thường xét A = Z hoặc A = Q hoặc A = R hoặc A = C. Khi đó
ta có các vành đa thức tương ứng là Z [x] , Q [x] , R [x] , C [x].
Định nghĩa 1.2 (xem [1]). Cho hai đa thức
f(x) = a
n
x
n

+ a
n−1
x
n−1
+ · · · + a
1
x + a
0
,
g(x) = b
n
x
n
+ b
n−1
x
n−1
+ · · · + b
1
x + b
0
.
Ta định nghĩa các phép tính số học
f(x) + g(x) = (a
n
+ b
n
)x
n
+ (a

n−1
+ b
n−1
)x
n−1
+ · · · + (a
1
+ b
1
)x + a
0
+ b
0
,
f(x) − g(x) = (a
n
− b
n
)x
n
+ (a
n−1
− b
n−1
)x
n−1
+ · · · + (a
1
− b
1

)x + a
0
− b
0
,
f(x)g(x) = c
2n
x
2n
+ c
2n−1
x
2n−1
+ · · · + c
1
x + c
0
,
trong đó
c
k
= a
0
b
k
+ a
1
b
k−1
+ · · · + a

k
b
0
, k = 0, 1, . . . , n.
6
Định lý 1.1 (xem [1]). Giả sử A là một trường, f(x) và g(x) = 0 là hai đa
thức của vành A [x], thế thì bao giờ cũng có cặp đa thức duy nhất g(x) và
r(x) thuộc A [x] sao cho
f(x) = g(x)q(x) + r(x) với deg r(x) < deg g(x).
Nếu r(x) = 0 ta nói f(x) chia hết cho g(x).
Giả sử f(x) = a
n
x
n
+ a
n−1
x
n−1
+ · · · + a
1
x + a
0
là đa thức tùy ý của vành
A [x], a là phần tử tùy ý của vành A, phần tử f(a) = a
n
a
n
+ a
n−1
a

n−1
+
· · · + a
1
a + a
0
có được bằng cách thay x bởi a được gọi là giá trị của f(x)
tại a.
Nếu f(a) = 0 thì ta gọi a là nghiệm của f(x). Bài toán tìm các nghiệm của
f(x) trong A gọi là giải phương trình đại số bậc n
a
n
x
n
+ a
n−1
x
n−1
+ · · · + a
1
x + a
0
= 0 (a
n
= 0)
trong A.
Định lý 1.2 (xem [1]). Giả sử A là một trường, a ∈ A và f(x) ∈ A [x]. Dư
số của phép chia f(x) cho x − a chính là f(a).
Định lý 1.3 (xem [1]). Số a là nghiệm của f(x) khi và chỉ khi f(x) chia hết
cho x − a.

Giả sử A là một trường, a ∈ A và f(x) ∈ A [x] và m là một số tự nhiên
lớn hơn hoặc bằng 1. Khi đó a là nghiệm bội cấp m của f(x) khi và chỉ khi
f(x) chia hết cho (x − a)
m
và f(x) không chia hết cho (x − a)
m+1
.
Trong trường hợp m = 1 thì ta gọi a là nghiệm đơn còn khi m = 2 thì a
được gọi là nghiệm kép. Số nghiệm của một đa thức là tổng số nghiệm của
đa thức đó kể cả bội của các nghiệm (nếu có). Vì vậy người ta coi một đa
thức có một nghiệm bội cấp m như một đa thức có m nghiệm trùng nhau.
Định lý 1.4 (xem [1]). Mỗi đa thức thực bậc n đều có không quá n nghiệm
thực.
Hệ quả 1.1. Đa thức có vô số nghiệm là đa thức không.
Hệ quả 1.2. Nếu đa thức f(x) có deg ≤ n mà nhận cùng một giá trị như
nhau tại n + 1 điểm phân biệt của đối số thì đó là đa thức hằng.
Hệ quả 1.3. Hai đa thức có deg ≤ n mà nhận n + 1 giá trị trùng nhau tại
n + 1 điểm phân biệt của đối số thì chúng đồng nhất bằng nhau.
7
Định lý 1.5 (xem [1]). Mọi đa thức f(x) ∈ C [x] bậc n có đúng n nghiệm
(tính cả bội của nghiệm).
Định lý 1.6 (xem [1]). Mọi đa thức f(x) ∈ R [x] có bậc n và có hệ số chính
(hệ số bậc cao nhất) a
n
= 0 đều có thể phân tích (duy nhất) thành nhân tử
dạng
f(x) = a
n
m


i=1
(x − d
i
)
s

k=1

x
2
+ b
k
x + c
k

với d
i
, b
k
, c
k
∈ R, 2s + m = n, b
2
k
− 4c
k
< 0, l, m, n ∈ N
∗
.
1.1 Công thức nội suy Taylor

Bây giờ, ta chuyển sang xét bài toán nội suy Taylor.
Bài toán 1.1 (Bài toán nội suy Taylor, xem [2]). Cho x
0
, a
k
∈ R với k =
0, 1, . . . , N − 1. Hãy xác định đa thức T (x) bậc không quá N − 1 (tức là
deg T (x)  N − 1) và thỏa mãn các điều kiện
T
(k)
(x
0
) = a
k
, ∀k = 0, 1, . . . , N − 1.
Giải. Ta dễ thấy rằng mọi đa thức đều viết được dưới dạng
T (x) =
N−1

k=0
α
k
(x − x
0
)
k
có deg T (x)  N − 1.
Ta cần đi xác định các hệ số α
k
∈ R sao cho T (x) thỏa mãn điều kiện

T
(k)
(x
0
) = a
k
, ∀k = 0, 1, . . . , N − 1.
Lần lượt lấy đạo hàm T (x) đến cấp thứ k, k = 0, 1, . . . , N − 1, tại x = x
0
và sử dụng giả thiết
T
(k)
(x
0
) = a
k
, ∀k = 0, 1, . . . , N − 1,
ta suy ra
α
k
=
a
k
k!
, ∀k = 0, 1, . . . , N − 1.
Thay giá trị của α
k
vào biểu thức của T (x), ta thu được
T (x) =
N−1


k=0
a
k
k!
(x − x
0
)
k
. (1.1)
8
Với mỗi k = 0, 1, . . . , N − 1, ta có
T
(k)
(x) = a
k
+
N−1

j=k+1
a
j
(j − k)!
(x − x
0
)
j−k
.
Do vậy đa thức T (x) thỏa mãn điều kiện
T

(k)
(x
0
) = a
k
, ∀k = 0, 1, . . . , N − 1.
Cuối cùng, ta phải chứng minh đa thức T (x) nhận được từ (1.1) là đa thức
duy nhất thỏa mãn điều kiện của bài toán nội suy Taylor (1.1), và ta gọi đa
thức này là đa thức nội suy Taylor.
Thật vậy, nếu có đa thức T
∗
(x), có bậc deg T
∗
(x)  N − 1 cũng thỏa
mãn điều kiện của bài toán (1.1) thì khi đó, đa thức
P (x) = T (x) − T
∗
(x)
cũng có bậc deg P (x)  N − 1 và đồng thời thỏa mãn điều kiện
P
(k)
(x
0
) = 0, ∀k = 0, 1, . . . , N − 1,
tức là, đa thức P (x) là đa thức có bậc không quá N − 1 (deg P(x)  N − 1)
mà lại nhận x
0
làm nghiệm với bội không nhỏ thua N, nên P (x) ≡ 0, và do
đó T (x) = T
∗

(x).
Nhận xét 1.1. Chú ý rằng đa thức nội suy Taylor T (x) được xác định từ
(1.1) chính là khai triển Taylor đến cấp thứ N −1 của đa thức T (x) tại điểm
x = x
0
.
1.2 Công thức nội suy Newton
Bài toán 1.2 (Bài toán nội suy Newton, xem [2]). Cho x
i
, a
i
∈ R với i =
0, 1, . . . n. Hãy xác định đa thức N(x) có bậc không quá n (deg N(x) ≤ n)
và thỏa mãn các điều kiện
N
(i)
(x
i
) = a
i
, i = 0, 1, . . . , n. (1.2)
Giải. Dễ dàng chứng minh các đẳng thức sau đây
N(x) = N(x
0
) +
x

x
0
N


(t
1
) dt
1
,
9
N

(t
1
) = N

(x
1
) +
t
1

x
1
N

(t
2
) dt
2
,
N


(t
2
) = N

(x
2
) +
t
2

x
2
N

(t
3
) dt
3
,

Từ đó ta thu được đa thức cần tìm có dạng
N(x) = a
n
x

x
0
t
1


x
1
t
2

x
2
. . .
t
n−1

x
n−1
dt
n
dt
n−1
. . . dt
1
+
+ a
n−1
x

x
0
t
1

x

1
t
2

x
2
. . .
t
n−2

x
n−2
dt
n−1
. . . dt
1
+ · · · + a
1
x

x
0
dt
1
+a
0
Ký hiệu
R
i
(x

0
, x
1
, x
2
, . . . x
i−1
, x) =
x

x
0
t
1

x
1
t
2

x
2
. . .
t
n−1

x
i−1
dt
1

. . .dt
i
.
Khi đó
N(x) = a
0
+ a
1
R(x
0
, x) + · · · + a
n−1
R
n−1
(x
0
, x
1
, . . . , x
n−2
, x)
+ a
n
R
n
(x
0
, x
1
, . . . , x

n−1
, x).
Đây là đa thức duy nhất thỏa mãn điều kiện của bài toán nội suy Newton
và ta gọi đa thức này là đa thức nội suy Newton.
Thật vậy, dễ thấy rằng deg N(x) ≤ n. Ngoài ra, ứng với mỗi i = 0, 1, . . . , n,
ta có
N
(i)
(x) = a
i
+ a
i+1
R (x
i
, x) + · · · + a
n
R
n−i
(x
i
, x
i+1
, . . . , x
n−1
, x) .
Từ đó suy ra
N
(i)
(x
i

) = a
i
, i = 0, 1, . . . , n.
Cuối cùng ta đi chứng minh tính duy nhất nghiệm của bài toán nội suy
Newton. Giả sử tồn tại đa thức N
1
(x), có deg N
1
(x) ≤ n cũng thỏa mãn điều
kiện bài toán. Khi đó, đa thức P (x) = N(x) − N
1
(x) cũng có deg P (x) ≤ n
và thỏa mãn điều kiện
P
(i)
(x
i
) = 0, ∀i = 0, 1, . . . , n.
10
Vậy, theo cách xây dựng đa thức N(x) bậc cao nhất bằng 0, ứng với trường
hợp a
i
= 0; i = 0, 1, . . . , n, ta suy ra P (x) ≡ 0 và do đó N(x) = N
1
(x).
Nhận xét 1.2. Trong trường hợp n cụ thể ta nhận được dạng của đa thức
nội suy Newton tương ứng.
*) Với n = 1 thì N(x) = a
0
+ a

1
R(x
0
, x) = a
0
+ a
1
(x − x
0
).
*) Với n = 2 thì
N(x) = a
0
+ a
1
R(x
0
, x) + a
2
R
2
(x
0
, x
1
, x)
= a
0
+ a
1

(x − x
1
) +
a
2
2!

(x − x
1
)
2
− (x
0
− x
1
)
2

.
*) Với n = 3 thì
N(x) = a
0
+ a
1
R(x
0
, x) + a
2
R
2

(x
0
, x
1
, x) + a
3
R
3
(x
0
, x
1
, x
2
, x)
= a
0
+ a
1
(x − x
0
) +
a
2
2!

(x − x
1
)
2

− (x
0
− x
1
)
2

+
+
a
3
3!

(x − x
2
)
3
− 3x(x
2
− x
1
)
2
+ (x
2
− x
0
)
2
+ 3x

0
(x
1
− x
2
)
2

.
Ví dụ 1.1. Xác định các tam thức bậc hai f(x) thỏa mãn các điều kiện
f
(n)
(2n + 1) = (−1)
n
(2n
2
− n − 1), n = 0, 1, 2.
Giải. Từ giả thiết bài toán, suy ra tam thức bậc hai f(x) cần xác định thỏa
mãn các điều kiện
f(0) = −1; f

(3) = 0; f

(5) = 5.
Theo công thức nội suy Newton trong trường hợp n = 2, ta nhận được tam
thức bậc hai
f(x) = −1 + 5

(x − 3)
2

2
−
4
2

hay
f(x) =
5
2
x
2
− 15x +
23
2
.
Ví dụ 1.2. Phát biểu một bài toán nội suy Newton ứng với đa thức nội suy
có dạng
P (x) = x(x + 1)(x + 2).
Giải. Ta chọn các nút nội suy x
0
= 0, x
1
= 1, x
2
= 2, x
3
= 3. Khi đó
P (x) = x
3
+ 3x

2
+ 2x
P

(x) = 3x
2
+ 6x + 2
P

(x) = 6x + 6
P

(x) = 6
11
Suy ra
P (0) = 0, P

(1) = 11, P

(2) = 18, P

(3) = 6.
Do vậy, bài toán được đặt ra là Xác định đa thức P (x) có bậc không quá 3
(deg P (x) ≤ 3) sao cho
P (0) = 0, P

(1) = 11, P

(2) = 18, P


(3) = 6.
1.3 Khai triển Taylor - Gontcharov
Bài toán 1.3. Giả sử hàm f xác định trên tập Ω ⊂ R trong đó Ω là
hợp của các khoảng mở trên trục thực. Giả sử f khả vi cấp i tại điểm
x
i
∈ Ω, i = 0, 1, 2, . . . , n. Hãy xác định các đa thức P
n
(x) có deg P (X) ≤ n
sao cho P
(i)
n
(x
i
) = f
(i)
(x
i
), i = 0, 1, . . . , n.
Giải. Đặt f
(i)
n
(x
i
) = a
i
với i = 0, 1, . . . , n. Khi đó đa thức
P
n
(x) = a

n
R
n
(x
0
, . . . , x
n−1
, x) + a
n−1
R
n−1
(x
0
, . . . , x
n−2
, x)
+ · · · + a
1
R(x
0
, x) + a
0
là đa thức duy nhất thỏa mãn yêu cầu của bài toán.
Ví dụ 1.3. Xét hàm số f(x) = ln x với nút nội suy x
0
= 2; x
1
= 3; x
2
= 4.

Giải. Ta có f
(k+1)
(x) =
(−1)
k
k!
x
k+1
, ∀k = 0, 1, . . . , n.
Khi đó đa thức P
2
(x) thỏa mãn P
(i)
2
(x
i
) = f
(i)
(x
i
), (i = 0, 1, 2) là duy nhất
và có dạng
P
2
(x) = f

(4)
x

2

t
1

3
dt
1
dt
2
+ f

(3)
x

2
dt
1
+ f(2)
= −
1
16

x
2
2
− 3x + 4

+
1
3
(x − 2) + ln 2.

Bài toán 1.4. Cho x
i
∈ R với i = 0, 1, . . . , n. Chứng minh rằng đa thức
w
n+1
(x) có bậc n + 1(deg w
n+1
(x) = n + 1), hệ số cao nhất bằng 1 và thỏa
mãn điều kiện w
(k)
n+1
(x
k
) = 0, k = 0, 1, . . . , n là tồn tại duy nhất và được xác
định bởi
w
n+1
(x) = (n + 1)!R
n+1
(x
0
, x
1
, . . . x
n
, x),
trong đó
R
n+1
(x

0
, x
1
, . . . x
n
, x) =
x

x
0
t
1

x
1
. . .
t
n

x
n
dt
n+1
dt
n
. . . .dt
1
.
12
Chứng minh.

*) Chứng minh tính tồn tại nghiệm
Đặt
w
n+1
(x) = (n + 1)!R
n+1
(x
0
, x
1
, . . . , x
n
, x)
= (n + 1)!
x

x
0
t
1

x
1
. . .
t
n

x
n
dt

n+1
.dt
n
. . . dt
2
dt
1
Dễ thấy deg w
n+1
(x) = n + 1 và hệ số của x
n+1
bằng 1. Ngoài ra ứng với mỗi
k = 0, 1, . . . , n ta có
w
(k)
n+1
(x) = (n + 1)!R
n−k+1
(x
k
, x
k+1
, . . . x
n
, x)
= (n + 1)!
x

x
k

t
k+1

x
k+1
t
k+2

x
k+2
. . .
t
n

x
n
dt
n+1
dt
n
. . . .dt
k
.
Suy ra
w
(k)
n+1
(x) = 0.
Do đó, đa thức w
n+1

(x) thỏa mãn w
(k)
n+1
(x
k
) = 0, k = 0, 1, . . . , n.
*) Chứng minh tính duy nhất nghiệm
Giả sử tồn tại đa thức w
∗
n+1
(x) có deg w
∗
n+1
(x) = n + 1 cũng thỏa mãn điều
kiện bài toán. Khi đó, đa thức P (x) = w
n+1
(x) − w
∗
n+1
(x) có deg P (x) ≤ n
và thỏa mãn điều kiện
P
(k)
(x
k
) = 0, k = 0, 1, . . . , n.
Suy ra P (x) ≡ 0 và do đó w
n+1
(x) ≡ w
∗

n+1
(x). 
Định nghĩa 1.3 (xem[2]). Đa thức
P
n
(f; x) =f
(n)
(x
n
)R
n
(x
0
, x
1
, x
2
, . . . x
n−1
, x)+
+f
(n−1)
(x
n−1
)R
n−1
(x
0
, x
1

, x
2
, . . . x
n−2
, x)+
+ · · · + f

(x
1
)R(x
0
, x) + f(x
0
),
được gọi là đa thức Newton theo bộ nội suy x
0
, x
1
, . . . , x
n
của hàm f, (f khả
vi cấp i tại điểm x
i
), i = 0, 1, . . . , n.
Ta đặt
f(x) =P
n
(f; x) + P
n+1
(f; x)

=f
(n)
(x
n
)R
n
(x
0
, x
1
, x
2
, . . . x
n−1
, x)+
+f
(n−1)
(x
n−1
)R
n−1
(x
0
, x
1
, x
2
, . . . x
n−2
, x)+

+ · · · + f

(x
1
)R(x
0
, x) + f(x
0
) + R
n+1
(f, x). (1.3)
13
Công thức (1.3) được gọi là công thức khai triển Taylor - Gontcharov.
Biểu thức R
n+1
(f; x) được gọi là phần dư của công thức khai triển Taylor
- Gontcharov, lời giải của bài toán ước lượng hiệu f(x) −P
n
(f; x) cũng chính
là ước lượng các biểu thức phần dư này.
Định lý 1.7 (xem [2]). Giả sử f : (a; b) → R khả vi và liên tục đến cấp
n + 1 trên khoảng (a; b); x
i
∈ (a; b), i = 0, 1, . . . , n. Khi đó
f(x) = P
n
(f; x) + R
n+1
(f; x),
trong đó R

n+1
(f; x) được xác định bởi
R
n+1
(f; x) =
x

x
0
t
1

x
1
t
2

x
2
. . .
t
n

x
n
f
(n+1)
(t
n+1
).dt

n+1
.dt
n
. . . dt
2
dt
1
. (1.4)
Chứng minh. Ta có
f(x) = f(x
0
) +
x

x
0
f

(t
1
)dt
1
= f(x
0
) +
x

x
0



f

(x
1
) +
t
1

x
1
f

(t
2
)dt
2


dt
1
= f(x
0
) +
x

x
0
f


(x
1
)dt
1
+
x

x
0
t
1

x
1
f

(t
2
)dt
2
dt
1
= f(x
0
) + f

(x
1
)R
1

(x
0
, x) +
x

x
0
t
1

x
1
f

(t
2
)dt
2
dt
1
= f(x
0
) + f

(x
1
)R
1
(x
0

, x) +
x

x
0
t
1

x
1


f

(t
2
) +
t
2

x
2
f

(t
3
)dt
3



dt
2
dt
1
= f(x
0
)+f

(x
1
)R
1
(x
0
, x) + f

(x
2
)R
2
(x
0
, x
1
, x)+
+
x

x
0

t
1

x
1
t
2

x
2
f

(t
3
)dt
3
dt
2
dt
1
= . . .
= f(x
0
)+f

(x
1
)R
1
(x

0
, x) + f

(x
2
)R
2
(x
0
, x
1
, x) + · · · +
+f
(n)
(x
n
)R
n
(x
0
, x
1
, . . . , x
n
) +
x

x
0
t

1

x
1
t
2

x
2
. . .
t
n

x
n
f
(n+1)
(t
n+1
)dt
n+1
.dt
n
. . . dt
2
dt
1
.
14
Vậy f(x) = P

n
(x) + R
n+1
(f; x) với
R
n+1
(f; x) =
x

x
0
t
1

x
1
t
2

x
2
. . .
t
n

x
n
f
(n+1)
(t

n+1
).dt
n+1
.dt
n
. . . dt
2
dt
1
. 
Nhận xét 1.3. Biểu thức trong (1.4) cho ta công thức xác định phần dư
R
n+1
(f, x) trong khai triển Taylor - Gontcharov của hàm f(x). Câu hỏi tự
nhiên đặt ra là liệu ta có thể đưa ra công thức đánh giá phần dư này giống
như trong khai triển Taylor ở dạng Lagrange và Cauchy được không?
Trong phần tiếp theo ta sẽ tìm cách đưa ra các đánh giá công thức phần dư
(1.4) dưới dạng Lagrange và Cauchy.
Trước hết, để đơn giản ta định nghĩa phép toán
(R
i
f)(x) =
x

x
1
f(s)ds, i = 0, 1, 2, . . . n.
Dễ dàng nhận thấy
(R
0

R
1
R
2
. . . R
n
f)(x) =
x

x
0
t
1

x
1
t
2

x
2
. . .
t
n

x
n
f(t
n+1
).dt

n+1
dt
n
. . .dt
1
,
(R
i
)(x) =
x

x
i
ds.
Và
(R
0
R
1
R
2
. . . R
n
)(x) =
x

x
0
t
1


x
1
t
2

x
2
. . .
t
n

x
n
dt
n+1
dt
n
. . . dt
1
=R
n+1
(x
0
, x
1
, . . . , x
n−1
, x).
Mệnh đề 1.1 (xem [4]). Với (R

i
f)(x) =
x

x
i
f(s)ds, ta luôn có
(R
k
f)(x) = (R
0
f)(x) − (R
0
f)(x
k
). (1.5)
(R
i
R
k
f)(x) = (R
i
R
0
f)(x) − (R
i
)(x)(R
0
f)(x
k

). (1.6)
Chứng minh. Ta đi chứng minh đẳng thức (1.5), ta có
(R
k
f)(x) =
x

x
k
f(s)ds =
x

x
0
f(s)ds −
x
k

x
0
f(s)ds = (R
0
f)(x) − (R
0
f)(x
k
).
15
Chứng minh đẳng thức (1.6), ta có
(R

i
R
k
f)(x) =
x

x
i
(R
k
f)(t)dt =
x

x
i


t

x
0
f(s)ds −
x
k

x
0
f(s)ds



dt
=
x

x
i
t

x
0
f(s)dsdt −
x

x
i
dt
x
k

x
0
f(s)ds
=(R
i
R
0
f)(x) − (R
i
)(x)(R
0

f)(x
k
) 
Mệnh đề 1.2 (xem [4]). Ta luôn có đẳng thức sau
(R
0
R
1
. . . R
n
f)(x) = (R
n+1
0
f)(x) −
n

k=1
(R
0
R
1
. . . R
k−1
)(x)(R
n−k+1
0
f)(x
k
),
(1.7)

trong đó
(R
0
R
1
. . . R
n
f)(x) =
x

x
0
t
1

x
1
t
2

x
2
. . .
t
n

x
n
f(t
n+1

).dt
n+1
dt
n
. . .dt
1
.
Chứng minh. Ta đi chứng minh (1.7) bằng phương pháp quy nạp.
Thật vậy,
*) Với n = 1 thì từ (1.6) ta suy ra
(R
0
R
1
f)(x) = (R
2
0
f)(x) − R
0
(x)(R
0
f)(x)
nên (1.7) đúng.
*) Giả sử (1.7) đúng với n, tức là
(R
0
R
1
. . . R
n

f)(x) = (R
n+1
0
f)(x) −
n

k=1
(R
0
R
1
. . . R
k−1
)(x)(R
n−k+1
0
f)(x
k
).
*) Ta đi chứng minh (1.7) đúng với n + 1. Ta có
(R
0
R
1
. . . R
n
R
n+1
f)(x) = R
0

R
1
. . . R
n−1
[(R
n
R
0
f)(x) − R
n
(x)(R
0
f)(x
n+1
)]
= R
0
R
1
. . . R
n
(R
0
f)(x) − R
0
R
1
. . . R
n
(R

0
f)(x
n+1
)
= R
n+1
0
(R
0
f)(x) −
n

k=1
(R
0
R
1
. . . R
k−1
)(x)(R
n−k+1
0
R
0
f)(x
k
)−
− R
0
R

1
. . . R
n
(R
0
f)(x
n+1
)
= (R
n+2
0
f)(x) −
n+1

k=1
(R
0
R
1
. . . R
k−1
)(x)(R
n−k+1
0
f)(x
k
)
16
Suy ra (1.7) đúng với n + 1.
Vậy,

(R
0
R
1
. . . R
n
f)(x) = (R
n+1
0
f)(x) −
n

k=1
(R
0
R
1
. . . R
k−1
)(x)(R
n−k+1
0
f)(x
k
).
Mệnh đề 1.3 (xem [4]). Ta luôn có đẳng thức sau
(R
k
0
f)(x) =

x

x
0
(x − s)
k−1
(k − 1)!
f(s)ds, k = 1, 2, . . . . (1.8)
Chứng minh. Ta đi chứng minh (1.8) bằng phương pháp qui nạp theo k.
Thật vậy,
*) Với k = 1 thì (1.8) luôn đúng.
*) Giả sử (1.8) đúng với n, tức là
(R
n
0
f)(x) =
x

x
0
(x − s)
n−1
(n − 1)!
f(s)ds.
*) Ta đi chứng minh (1.8) đúng vơi n + 1. Theo công thức tính tích phân
từng phần, ta có
x

x
0

(x − s)
n
n!
f(s)ds =
(x − s)
n
n!
s

x
0
f(r)dr


x
x
0
+
x

x
0
(x − s)
n−1
(n − 1)!


s

x

0
f(r)dr


ds
= 0 +
x

x
0
(x − s)
n−1
(n − 1)!


s

x
0
f(r)dr


ds
= (R
n
0


s


x
0
f(r)dr


)(x) = (R
n
0
R
0
f)(x)
= (R
n+1
0
f)(x).
Vậy (1.8) đúng với n + 1. Từ đó ta có điều phải chứng minh. 
Hệ quả 1.4. Từ các hệ thức (1.7) và (1.8) ta nhận được kết quả quan trọng
sau
(R
0
R
1
. . . R
n
f)(x) =
=
x

x
0

(x − s)
n
n!
f(s)ds −
n

k=1
R
k
(x
0
, . . . , x
k−1
, x)
x
k

x
0
(x
k
− s)
n−k
(n − k)!
f(s)ds.
17
Định lý 1.8 (xem [2]). Cho hàm số f(x) liên tục trên [a; b]. Khi đó, với mọi
x
0
; x ∈ [a; b] luôn tồn tại ξ = x

0
+ θ(x − x
0
) với 0 < θ < 1 sao cho
(R
k+1
0
f)(x) = f(ξ)
(x − x
0
)
k+1
(k + 1)!
, k = 0, 1, . . .
Chứng minh. Từ (1.8) ta có
(R
k+1
0
f)(x) =
x

x
0
(x − s)
k
k!
f(s)ds
Vì g(s) =
(x − s)
k

k!
không đổi dấu trên [x
0
; x]. Theo định lý về giá trị trung
bình của tích phân ∃ ξ : ξ = x
0
+ θ(x − x
0
) với 0 < θ < 1 sao cho
(R
k+1
0
f)(x) = f(ξ)
x

x
0
(x − s)
k
k!
ds = f(ξ)
x

x
0
(x − x
0
)
k+1
(k + 1)!

ds. 
Hệ quả 1.5. Hàm số f(x) liên tục trên [a; b] , x
i
∈ [a; b] , i = 0, 1, 2, . . . n
và x là điểm bất kỳ [a; b]. Khi đó ta có
(R
0
R
1
. . . R
n
f)(x) =
f(ξ)
(x − x
0
)
n+1
(n + 1)!
−
n

k=1
R
k
(x
0
, . . . , x
k−1
, x)
f(ξ

k
)
(n − k + 1)!
(x
k
− x
0
)
n−k+1
,
(1.9)
trong đó ξ nằm giữa x
0
; x, còn ξ
k
nằm giữa x
0
; x
k
.
Thật vậy, từ hệ quả 1.4 ta có
(R
0
R
1
. . . R
n
f)(x) =
=
x


x
0
(x − s)
n
n!
f(s)ds −
n

k1
R
k
(x
0
, . . . , x
k−1
, x)
x
k

x
0
(x
k
− s)
n−k
(n − k)!
f(s)ds.
Làm tương tự như định lý trên ta nhận được điều phải chứng minh. 
1.3.1 Khai triển Taylor-Goncharov với phần dư dạng Lagrange

Định lý 1.9 (xem [2]). Giả sử hàm số f(x) có đạo hàm liên tục cấp n + 1
trên [a; b] và x
i
∈ [a; b] , i = 0, 1, 2, . . . n; x là điểm bất kỳ thuộc [a; b] . Khi
đó ∃ ξ; ξ
i
: ξ = x
0
+ θ(x − x
0
) và ξ
i
= x
0
+ θ
i
(x
i
− x
0
) với 0 < θ; θ
i
< 1 để
phần dư của hàm f(x) có dạng
18
R
n+1
(f; x) =
f
(n+1)

(ξ)
(n + 1)!
(x − x
0
)
(n+1)
−
−
n

k=1
f
(n+1)
(ξ
k
)
(n − k + 1)!
R
k
(x
0
, . . . , x
k−1
, x)(x
k
− x
0
)
(n−k+1)
.

Vậy ta nhận được phần dư của công thức khai triển Taylor-Gontcharov của
hàm f(x) dưới dạng Lagrange như sau
R
n+1
(f; x) =
f
(n+1)
(ξ)
(n + 1)!
(x − x
0
)
(n+1)
−
−
n

k=1
f
(n+1)
(ξ
k
)
(n − k + 1)!
R
k
(x
0
, . . . , x
k−1

, x)(x
k
− x
0
)
(n−k+1)
.
(Phần dư dạng Lagrange)
Nhận xét 1.4. Trong trường hợp đặc biệt, khi x
i
= x
0
, ∀
i
= 1, 2, . . . , n thì
công thức phần dư dạng Lagrange trong khai triển Taylor-Gontcharov của
hàm f(x) là
R
n+1
(f; x) =
f
(n+1)
(ξ)
(n + 1)!
(x − x
0
)
(n+1)
sẽ trùng với công thức phần dư dạng Lagrange trong khai triển Taylor của
hàm f(x) tại điểm x

0
mà ta đã biết ở trong mục trước.
1.3.2 Khai triển Taylor-Goncharov với phần dư dạng Cauchy
Định lý 1.10 (xem [2]). Cho hàm số f(x) liên tục trên [a; b] , ∀x, x
0
∈ [a; b].
Khi đó, ∃θ ∈ (0, 1) sao cho
(R
n+1
0
f)(x) = f (x
0
+ θ(x − x
0
))
(1 − θ)
k
(x − x
0
)
k+1
k!
, k = 0, 1, 2 . . .
Hệ quả 1.6. Giả sử hàm số f(x) liên tục trên [a; b], x
i
∈ [a; b], i = 0, 1, 2 . . .
và x là điểm bất kỳ thuộc [a; b]. Khi đó, ta có
(R
0
R

1
. . . R
n
f)(x) = f (x
0
+ θ(x − x
0
))
(1 − θ)
n
(x − x
0
)
n+1
n!
−
−
n

k=1
R
k
(x
0
, . . . , x
k−1
, x)f (x
0
+ θ(x
k

− x
0
))
(1 − θ
k
)
n−k
(n − k)!
(x
k
− x
0
)
n−k+1
,
(1.10)
trong đó θ, θ
k
∈ (0, 1), k = 0, 1, 2 . . .
19
Định lý 1.11 (xem [4]). Giả sử hàm số f(x) có đạo hàm liên tục đến cấp
n + 1 trên [a; b] và x
i
∈ [a; b], i = 0, 1, 2 . . . , n; x là điểm bất kỳ thuộc [a; b].
Khi đó, phần dư dạng Cauchy trong khai triển Taylor-Gontcharov của hàm
f(x) có dạng
R
n+1
(f; x) = f
n+1

(x
0
+ θ(x − x
0
))
(1 − θ)
n
(x − x
0
)
n+1
n!
−
−
n

k=1
R
k
(x
0
, . . . , x
k−1
, x)f
n+1
(x
0
+ θ(x
k
− x

0
))
(1 − θ
k
)
n−k
(n − k)!
(x
k
− x
0
)
n−k+1
)
trong đó θ, θ
k
∈ (0, 1), k = 0, 1, 2 . . .
(Phần dư dạng Cauchy)
Nhận xét 1.5. Trong trường hợp đặc biệt, khi x
i
= x
0
, ∀i = 1, 2, . . . , n
thì công thức phần dư dạng Cauchy trong khai triển Taylor-Gontcharov của
hàm f(x) là
R
n+1
(f; x) = f
n+1
(x

0
+ θ(x − x
0
))
(1 − θ)
n
(x − x
0
)
n+1
n!
; 0 < θ < 1
sẽ trùng với công thức phần dư Cauchy trong khai triển Taylor của hàm số
f(x) tại điểm x
0
mà ta đã biết trong mục trước.
20
Chương 2
Bài toán biên hỗn hợp thứ nhất của
phương trình vi phân
Như ta đã biết ở chương 1, luận văn đã xét bài toán nội suy Taylor và
bài toán nội suy Newton. Ta có thể xét bài toán nội suy như trường hợp đặc
biệt của phương trình vi phân. Và trong chương này ta sẽ đi xét bài toán
Cauchy và bài toán biên hỗn hợp thứ nhất dạng trừu tượng cho lớp toán tử
khả nghịch phải đồng thời trình bày các bài toán Cauchy và bài toán biên
hỗn hợp thứ nhất cho phương trình vi phân thường.
2.1 Bài toán Cauchy
Định nghĩa 2.1 (xem [4]). Toán tử D ∈ L(X) được gọi là khả nghịch phải
nếu tồn tại toán tử R ∈ L
0

(X) sao cho ImR ⊂ dom D và DR = I, trong
đó I là toán tử đồng nhất.
Mệnh đề 2.1 (xem [4]). Nếu D ∈ R(X) và R ∈ R
D
thì D
k
R
k
= I, với
∀k ∈ N
+
.
Mệnh đề 2.2 (xem [4]). Nếu D ∈ R(X), R
1
, R
2
∈ R
D
và y
1
= R
1
x,
y
2
= R
2
x với x ∈ X. Khi đó, y
1
− y

2
∈ ker D.
Mệnh đề 2.3 (xem [4]). Nếu D ∈ R(X) và R ∈ R
D
thì mỗi nguyên hàm
của x ∈ X có dạng
R
D
x = {Rx + z; z ∈ ker D} = Rx + ker D.
Mệnh đề 2.4 (xem [4]). Nếu D ∈ R(X) và R ∈ R
D
thì
dom D = RX ⊕ ker D.
Định lý 2.1 (xem [4]-[6]). Giả sử D ∈ R(X) và R
1
∈ R
D
. Khi đó, mỗi
nghịch đảo phải của D có dạng R = R
1
+(I −R
1
D)A với A ∈ L
0
(X); AX ⊂
dom D.
21
Ví dụ 2.1. Cho X là không gian tất cả các dãy số {x
n
} ⊂ R và x ∈

X,(x = {x
n
}) , n = 0, 1, . . . . Ta định nghĩa toán tử D như sau
D
x
= {x
n+1
− x
n
} = {x
1
− x
0
, x
2
− x
1
, . . .}
Khi đó, dom D = X và D là toán tử khả nghịch phải. Thật vậy, ta giả sử
Rx = {0, x
0
, x
0
+ x
1
, x
0
+ x
1
+ x

2
, }
=

y
n
: y
0
= 0, y
n
=
n−1

k=0
x
k
; n = 1, 2,

Ta có DRx = Dy = {y
n+1
− y
n
} = {x
n
}.
Do đó, DR = I nên D là toán tử khả nghịch phải.
Ví dụ 2.2. Cho
X = C [0, 1] ; D =
d
dt

+ I; (Rx)(t) =
t

t
0
e
s−t
x(s)ds; t
0
∈ [0, 1]
với t
0
cố định và x
n
∈ C [0, 1]. Khi đó, D là toán tử khả nghịch phải.
Thật vậy, ta có
(DRx) (t) =


t

t
0
e
s−t
x(s)ds


/
+

t

t
0
e
s−t
x(s)ds
=


e
−t
t

t
0
e
s
x(s)ds


/
+ e
−t
t

t
0
e
s

x(s)ds
= − e
−t
t

t
0
e
s
x(s)ds + e
−t
e
t
x(t) + e
−t
t

t
0
e
s
x(s)ds = x(t).
Do đó, DR = I nên D là toán tử khả nghịch phải.
Định nghĩa 2.2 (Toán tử ban đầu, xem[4]). Toán tử F ∈ L(X) được gọi là
toán tử ban đầu của toán tử D ∈ R(X) ứng với nghịch đảo phải R của D
nếu
i) F là một phép chiếu lên không gian các hằng số, nghĩa là
F
2
= F, F X = ker D

ii) F R = 0
22
Hệ quả 2.1. F z = z với mỗi z ∈ ker D.
Hệ quả 2.2. DF = 0 trên X.
Định lý 2.2 (xem [4]). Giả sử D ∈ R(X). Điều kiện cần và đủ để toán tử
F ∈ L
0
(X) là toán tử ban đầu của D tương ứng với R ∈ R
D
là
F = I − RD trên dom D
Hệ quả 2.3. Nếu A ∈ L(X) khả nghịch thì không tồn tại toán tử ban đầu
khác 0.
Hệ quả 2.4. Họ R
D
= {R
γ
}
γ∈Γ
tất cả các nghịch đảo phải của toán tử
D ∈ R(X) xác định duy nhất một họ F
D
= {F
γ
}
γ∈Γ
các toán tử ban đầu của
D.
F
γ

= I − R
γ
D trên dom D ; ∀γ ∈ Γ.
2.1.1 Bài toán Cauchy trừu tượng
Giả sử D ∈ R(X), và F là toán tử ban đầu của D ứng với nghịch đảo
phải R. Bài toán giá trị ban đầu của toán tử Q [D] như sau:
Tìm tất cả các nghiệm của phương trình
Q [D] x =
M

m=0
N

n=0
D
m
A
mn
D
n
x = y, y ∈
X
, (2.1)
thỏa mãn điều kiện ban đầu
F
j
D
j
x = y
j

, y
j
∈ ker D (j =0, 1, . . . , M + N − 1). (2.2)
Định nghĩa 2.3 (xem [4]). Bài toán giá trị ban đầu (2.1)-(2.2) được gọi là
"đặt chỉnh" nếu có nghiệm duy nhất với mỗi y ∈ C; y
0
, y
1
, . . . , y
M+N−1
∈ C.
Bài toán (2.1)-(2.2) được gọi là "đặt không chỉnh" nếu tồn tại y ∈ Q [D] X
M+N
và y
0
, . . . , y
M+N−1
∈ ker D sao cho bài toán này không có nghiệm hoặc bài
toán thuần nhất (tức là y = y
i
= 0; j = 0, . . . ., M + N − 1) có nghiệm
không tầm thường.
Định nghĩa 2.4 (xem [4]-[6]). Giả sử Q [D] có dạng (2.1). Đặt
K =
M

m=0
N

n=0

R
M−m
B
mn
R
N−n
, (2.3)
23
trong đó
B
mn
=





∧
A
0n
nếu m = 0
∧
A
mn
−
M

k=m
F D
k−m

∧
A
kn
các trường hợp khác
∧
A
mn
=

0 nếu m = M, n = N
A
mn
các trường hợp khác
(m = 0, 1, . . . , M; n = 0, 1, . . . , N). Khi đó, I + K được gọi là "toán tử
giải" của bài toán (2.1)- (2.2) .
Bổ đề 2.1 (xem[4]). Đặt
Q =
M

m=0
N

n=0
R
M+N−m
B
mn
D
n
. (2.4)

Khi đó
QR
N
= R
N
K, (2.5)
trong đó K xác định bởi (2.3) và
D
M+N
(I + Q) = Q [D] , (2.6)
F D
j
(I + Q) = F D
j
(j = 0, 1, . . . , M + N − 1). (2.7)
Bổ đề 2.2 (xem[4]). Cho Q xác định bởi (2.4). Khi đó, bài toán giá trị ban
đầu (2.1)-(2.2) "đặt chỉnh" khi và chỉ khi I + Q khả nghịch trên X
M+N
.
Chứng minh. Theo bổ đề 2.1 ta có thể viết phương trình (2.2) dưới dạng
D
M+N
(I + Q)x = y.
Phương trình này tương đương với phương trình
(I + Q)x = R
M+N
y +
M+N−1

j=0

R
j
z
j
trong đó z
0
, z
1
, . . . , z
M+N−1
∈ kerD tùy ý.
Từ công thức (2.7) và phương trình cuối, ta suy ra bài toán giá trị ban đầu
(2.1)-(2.2) tương đương với phương trình
(I + Q)x = R
M+N
y +
M+N−1

j=0
R
j
y
j
. (2.8)
24