Một số phương pháp lặp giải phương trình và hệ phương trình phi tuyến

1



























2
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 TS.Khuất Văn Ninh.
Tôi xin chân thành cảm ơn sâu sắc TS.Khuất Văn Ninh, người đã luôn
quan tâm, động viên, tận tình hướng dẫn tôi trong quá trình làm luận văn.
Tôi xin gửi lời cảm ơn chân thành tới Ban giám hiệu trường Đại học sư
phạm Hà Nội 2, phòng sau đại học, các thầy cô giáo trong nhà trường, các
thầ
y cô giáo dạy cao học chuyên ngành Toán giải tích đã tạo điều kiện thuận
lợi cho tôi trong quá trình học tập và nghiên cứu.
Tôi bày tỏ lòng biết ơn tới gia đình, người thân, bạn bè cùng học, đội ngũ
bảo vệ an ninh khu vực đã động viên và tạo mọi điều kiện thuận lợi để tôi
hoàn thành bản luận văn này.

Hà Nội, tháng 6 năm 2011



Lê Th
ị Thu Phương









3
LỜI CAM ĐOAN

Tôi xin cam đoan bản luận văn này là công trình nghiên cứu của riêng
tôi dưới sự hướng dẫn của TS.Khuất Văn Ninh.

Hà Nội, tháng 6 năm 2011



Lê Thị Thu Phương


















4
MỤC LỤC

Trang
Trang phụ bìa…………………………………………………………… 1
Lời cảm ơn ……………………………………………………………… 2
Lời cam đoan …………………………………………………………… 3
Mục lục…………………………………………………………………… 4
Bảng các kí hiệu ………………………………………………………… 5
Mở đầu…………………………………………………………………… 6
Nội dung ………………………………………………………………… 8
Chương 1: Kiến thức cơ sở …………………………………………….8
1.1 Không gian véc tơ ………………………………………………… 8
1.2 Các không gian quan trọng ……………………………………… 11
1.3 Đạo hàm Gateaux và đạo hàm Frechet ……………………………15
Chương 2: Một số phương pháp lặp giải phương trình và hệ
phương trình phi tuyến………………………………… 19
2.1 Phương pháp Newton và một số
biến thể của nó………………… 19
2.2 Phương pháp cát tuyến…………………………………………… 27
2.3 Một số biến thể …………………………………………………….44
2.4 Phương pháp sử dụng tính liên tục của ánh xạ …………………….52
2.5 Các phương pháp đặc biệt đối với hàm một biến ………………….57
2.6 Bàn thảo về các phương pháp lặp giải phương trình và hệ
phương trình phi tuyến …………………………………………….62
Chương 3: Bài tập…………………………………………………… 67
Kết luận và kiến nghị…………………………………………………….71
Tài liệu tham khảo……………………………………………………….72



5

BẢNG CÁC KÝ HIỆU

R
Tập hợp số thực.
n
R
Không gian véc tơ Ơclit thực
n
chiều.
),(
#
YXL Không gian các ánh xạ tuyến tính từ
X
vào Y .
),( YXL Không gian các ánh xạ tuyến tính liên tục từ
X
vào Y .
),(:
##
KXLX = Không gian các phiếm hàm tuyến tính trên
X
.
),(:
*
KXLX = Không gian các phiếm hàm tuyến tính liên tục trên
X
.
),( KnmMat × Tập các ma trận m dòng n cột.
),( KnMat Tập các ma trận vuông n dòng n cột.
),( KnGL Tập các ma trận vuông cấp n không suy biến.

)(XL Không gian các toán tử tuyến tính từ
X
vào
X
.
C
Tập hợp số phức.
K
Tập hợp số thực hoặc phức.
*
K
Tập hợp số thực hoặc phức khác không.











6
MỞ ĐẦU

1. Lý do chọn đề tài
Vào những năm 70 của thế kỷ 20, một số nhà toán học đã nghiên cứu về
giải các phương trình và hệ phương trình dạng:

yFx
=
(1)
trong đó
F
là một toán tử từ một tập
X
đến một tập
Y
,
YyXx ∈
∈
,

Trường hợp đặc biệt của (1) là:

0
=
Fx (2)
với
.1,:), ,(
1
≥→= nRRffF
nn
n

hay:
⎪
⎩
⎪

⎨
⎧
=
=
0), ,(
0), ,(
1
11
nn
n
xxf
xxf
M

Với các hiểu biết ban đầu và qua tham khảo một số tài liệu liên quan, tôi
thấy:
Phạm vi ứng dụng của lý thuyết toán tử là rất rộng lớn. Phạm vi ứng dụng này
càng rộng và càng có hiệu lực thực tiễn trước sự phát triển nhanh chóng của
máy tính điện tử với sự phát triển mạnh mẽ các công trình nghiên cứu xấp xỉ.
Việc giải xấp xỉ các phương trình, hệ ph
ương trình dạng (2) là phù hợp với
năng lực của tôi. Có nhiều phương pháp giải, song phương pháp lặp là
phương pháp có thể lập trình trên máy tính điện tử. Vì vậy tôi đã chọn đề tài:
“Một số phương pháp lặp giải phương trình và hệ phương trình phi tuyến
n ẩn số”
2. Mục đích nghiên cứu
Nghiên cứu một số phương pháp lặp giải phương trình và hệ phương
trình phi tuyến n
ẩn số, ứng dụng vào các bài tập cụ thể có sử dụng máy tính

7
điện tử để giải. Đánh giá về những nghiên cứu khoa học của mình. Nêu ra
những đóng góp của đề tài. Đề xuất các kiến nghị.
3. Nhiệm vụ nghiên cứu
Nghiên cứu một số phương pháp lặp giải phương trình và hệ phương
trình phi tuyến n ẩn số.
4. Đối tượng và phạm vi nghiên cứu
Nghiên cứu một số phương pháp lặp giải phương trình và hệ phương
trình phi tuyến n ẩn số. Nghiên cứu các ứng dụng của lý thuyết vào các bài
toán cụ thể có sử dụng ngôn ngữ lập trình Maple hoặc Pascal.
5. Phương pháp nghiên cứu
Nghiên cứu lý thuyết, áp dụng lý thuyết vào giải toán số.
6. Dự kiến đóng góp mới của đề tài
Hệ thống hoá các phương pháp lặp giải phương trình và hệ phương trình
phi tuyến n ẩn số. Lập trình các bài toán trên máy tính điện tử bằng ngôn ngữ
lập trình Maple hoặc Pascal.













8

NỘI DUNG
CHƯƠNG 1: KIẾN THỨC CƠ SỞ
1.1. Không gian véc tơ.
1.1.1. Định nghĩa: Cho
X
là một tập hợp,
K
là một trường số. Một không
gian véc tơ
X
trên trường
K
là một tập hợp
X
cùng với hai ánh xạ:
(phép cộng):
yxyx
XXX
+
→
×
a),(

(phép nhân với một vô hướng):
xx
XXK
αα
a),(
→
×

thoả mãn các tiên đề sau:
KXzyx
∈
∀
∈
∀
β
α
,;,, .
1.
x
yy
x
+=+
(phép cộng có tính giao hoán)
2.
)()( zyxzyxzyx
+
+
=++=++ (phép cộng có tính kết hợp)
3.
xxxX
=
+=+∈∃
θ
θ
θ
:
(

θ
gọi là phần tử không của
X
)
4.
θ
=
+
−
=−+∈−∃ xxxxXx )()(: (
x
−
gọi là phần tử đối của
x
)
5.
)()( xx
β
α
β
α
=
6.
xxx
β
α
β
α
+=+ )(
7.

yxyx
α
α
α
+=+ )(
8.
xxx == 11 (1 là phần tử đơn vị của trường
K
)
Các phần tử:
θ
và
x
−
là duy nhất, ngoài ra xx 1)(
−
=
−

1.1.2. Cơ sở Hamel và số chiều của không gian véc tơ.
Cho
X
là không gian véc tơ trên trường
K
.
1.1.2.1. Định nghĩa: Giả sử
{
}
Xxx
n

⊂, ,
1
là một hệ véctơ tuỳ ý,
{}
K
n
⊂
α
α
, ,
1
là một bộ số tuỳ ý. Ta gọi
∑
=
=
n
i
ii
xy
1
α
là tổ hợp tuyến tính của
các véc tơ
n
xx , ,
1
trong
X
(hay còn nói là
y

biểu thị tuyến tính qua các véc
tơ
n
xx , ,
1
).

9
Hệ véc tơ
{}
Xxx
n
⊂, ,
1
được gọi là độc lập tuyến tính nếu
K
n
∈∀
α
α
, ,
1
, từ đẳng thức: nix
i
n
i
ii
,1,0
1
=∀=⇒=

∑
=
αθα
.
Hệ véc tơ
{}
Xxx
n
⊂, ,
1
được gọi là phụ thuộc tuyến tính nếu nó không
độc lập tuyến tính.
1.1.2.2. Định nghĩa: Tập
φ
≠
⊂ MXM , gọi là tập độc lập tuyến tính nếu mọi
hệ hữu hạn
{}
{
}
KMxx
nn
⊂⊂
α
α
, ,,, ,
11
từ đẳng thức
nix
i

n
i
ii
,1,0
1
=∀=⇒=
∑
=
αθα
.
Tập
X
M
⊂ là tập độc lập tuyến tính tối đại của
X
nếu nó là tập độc lập
tuyến tính và nếu có
B
là tập độc lập tuyến tính của
X
mà BM ⊆ thì
B
M
=

1.1.2.3. Định lí và định nghĩa: Cho hệ hữu hạn vectơ
{}
XxxX
m
⊂= , ,

11
khi đó số phần tử của mọi hệ con độc lập tuyến tính
tối đại của hệ
1
X
là bằng nhau. Số đó được gọi là hạng của hệ véc tơ
1
X
.
1.1.2.4. Định nghĩa: Tập
X
M
⊂ gọi là cơ sở Hamel của
X
nếu
M
là tập
độc lập tuyến tính và mọi phần tử của
X
đều biểu thị tuyến tính duy nhất qua
các phần tử của
M
(hay
M
là tập độc lập tuyến tính tối đại của
X
). Khi đó số
chiều của
X
, kí hiệu là Xmid , chính là số phần tử của

M
.
1.1.3. Không gian véc tơ con.
Cho
X
là không gian véc tơ trên trường
K
.
1.1.3.1. Định nghĩa: Tập
X
M
⊂ gọi là không gian véc tơ con của
X
nếu
M

đóng kín đối với các phép toán của
X
, nghĩa là:
MxKMx
MyxMyx
∈∀∈∀∈
∈
∀
∈
+
,;
,;
αα

(hay
M
là một không gian véc tơ với các phép toán cảm sinh từ
X
).
1.1.3.2. Định nghĩa: Giả sử
M
là một tập con của
X
. Khi đó tập tất cả các tổ

10
hợp tuyến tính của các phần tử thuộc
M
lập thành một không gian con của
X
, ký hiệu là
M
hay
Mnasp
.
Nhận xét:
M là giao của mọi không gian véc tơ con của
X
chứa
M
,
đó là không gian véc tơ con nhỏ nhất của
X
chứa

M
và gọi là không gian con
sinh bởi
M
(hay còn gọi là bao tuyến tính của
M
)
1.1.4. Ánh xạ và ma trận.
Cho
YX , là hai không gian véc tơ trên cùng trường
K
. Ánh xạ

YXF →: .
1.1.4.1. Định nghĩa:
F
gọi là ánh xạ tuyến tính nếu và chỉ nếu:
.2.,1,,;)()(
11
≥=∈∈∀=
∑∑
==
nniKXxxFxF
ii
n
i
ii
n
i
ii

ααα

1.1.4.2. Ma trận: Giả sử
X
và Y là các không gian hữu hạn chiều,

{}
n
uuu , ,
1
= là một cơ sở của
X
,
{
}
m
vvv , ,
1
=
là một cơ sở của Y ,
F
là ánh
xạ tuyến tính. Khi đó
F hoàn toàn được xác định bởi:
njvabuF
m
i
ijijj
,1;)(
1

===
∑
=

Nếu
YyyyXxxxxFy
mn
∈
=
∈
== ), ,(,), ,(),(
11
thì
.,1;
1
mixay
n
j
jjii
==
∑
=

Hay
⎪
⎪
⎩
⎪
⎪
⎨

⎧
+++=
+++=
+++=
nnmmmm
nn
nn
xaxaxay
xaxaxay
xaxaxay




2211
22221212
12121111


11
Ta có ma trận:
⎟
⎟
⎟
⎟
⎟
⎠
⎞
⎜
⎜

⎜
⎜
⎜
⎝
⎛
==
nmmm
n
n
ji
aaa
aaa
aaa
aA
L
MMMM
L
L
21
22221
11211
)(
và gọi là ma trận của ánh
xạ tuyến tính
F
đối với các cơ sở nói trên. Từ đó
F
hoàn toàn được xác định
bởi ma trận
)(

ji
aA =
này.
Định nghĩa:
{
}
)(:)( xFyXxYyFmIXF =∈∃∈== gọi là ảnh của
F
.
{}
YY
xFXxFKerF
θθ
=∈==
−
)()(
1
gọi là hạt nhân (hạch) của F .
hạng
=F
hạng
{
}
==
=
FmImiduF
n
jj 1
)(
hạng

A
.
Định nghĩa và tính chất: Cho ma trận
),( KnMatA
∈
thì
A
là ma trận không
suy biến nếu và chỉ nếu
A
là ma trận khả nghịch (hoặc nếu và chỉ nếu
0≠Ated ), nghĩa là tồn tại ma trận nghịch đảo
1−
A
của
A
để IAAAA ==
−− 11

(
I
là ma trận đơn vị).
1.1.5. Định nghĩa không gian afin: Cho
X
là không gian véc tơ thực.
Tập
X
A
⊂ gọi là tập afin nếu AytxtRtAyx ∈−+⇒
∈

∀
∈
∀
)1(,,
Tập
X
A
⊂ gọi là không gian con afin của
X
(hay không gian véc tơ
con của
X
) nếu và chỉ nếu
A
là tập afin chứa phần tử X
∈
θ

1.2. Các không gian quan trọng.
1.2.1. Không gian Metric.
1.2.1.1. Định nghĩa
Định nghĩa 1: Cho tập hợp
φ
≠
X , ánh xạ RXXd →
×
:
thoả mãn với
Xzyx ∈∀ ,, ta có:
1.

yxyxdyxd
=
⇔
=
≥ 0),(;0),( (tiên đề đồng nhất)
2.
),(),( xydyxd = (tiên đề đối xứng)
3.
),(),(),( yzdzxdyxd
+
≤ (tiên đề tam giác)

12
Ánh xạ
d
gọi là metric (hay khoảng cách) trên
X
. Các tiên đề 1,2,3 gọi
là hệ tiên đề metric. Tập
X
cùng với ánh xạ
d
như thế gọi là một không gian
metric, kí hiệu là
),( dX
Định nghĩa 2: Dãy điểm
{
}
n
x trong không gian metric

),( dX
được gọi là hội
tụ đến
Xx ∈
0
nếu
εε
<≥∀∈∃>∀ ),(,:,0
00
*
0
xxdnnNn
n
. Ta viết:

0
xxmil
n
n
=
∞→
hoặc
0
xx
n
n
⎯⎯→⎯
∞→
.
Dãy điểm

{}
),( dXx
n
⊂ gọi là dãy cơ bản (hay dãy Cauchy) nếu:
εε
<≥∀∈∃>∀ ),(,,:,0
0
*
0 mn
xxdnnmNn (hay 0),(
,
=
∞→
mn
nm
xxdmil )
Không gian metric
),( dX
là không gian đầy nếu mọi dãy cơ bản trong
nó đều hội tụ.
1.2.1.2. Nguyên lí Cantor về dãy hình cầu thắt dần.
Không gian metric
),( dX là không gian đầy khi và chỉ khi mọi dãy
hình cầu đóng thắt dần đều có một điểm chung duy nhất.
1.2.1.3. Nguyên lí Banach về ánh xạ co.
Định nghĩa: Cho hai không gian metric
),(
1
dX và ),(
2

dY . Ánh xạ: YXF →:
gọi là ánh xạ co nếu tồn tại
)1;0[
∈
α
sao cho:
XxxxxdxFxFd
∈
∀
≤
21211212
,;),(),(
α
.
Nguyên lí Banach về ánh xạ co: Mọi ánh xạ co
),(),(: dXdXF →

không gian metric đầy vào chính nó đều có một điểm bất động duy nhất
(nghĩa là tồn tại duy nhất một điểm
***
: xxFXx =∈ ).
1.2.2. Định nghĩa không gian định chuẩn và không gian Banach.
Định nghĩa 1: Ta gọi không gian định chuẩn
X
là không gian véc tơ
X
trên
trường
K
cùng với một ánh xạ

RX →:.
(gọi là chuẩn) thoả mãn các tiên đề
sau:

13
Với
KXyx
∈
∀
∈∀
α
,,
ta có:
1.
θ
=⇔=≥ xxx 0;0
2.
xx
αα
=
3.
yxyx +≤+

Nếu ta đặt
yxyxd −=),(
thì d là một mêtric trên
X
.
Định nghĩa 2: Dãy điểm
{

}
n
x trong không gian định chuẩn
X
gọi là dãy cơ
bản (hay dãy Cauchy) nếu:
0
,
=−
∞→
mn
nm
xxmil .
Định nghĩa 3: Dãy điểm
{
}
n
x
trong không gian định chuẩn
X
gọi là hội tụ
tới
Xx ∈ nếu 0=−
∞→
xxmil
n
n
.
Ta viết:
xxmil

n
n
=
∞→
hay xx
n
n
⎯⎯→⎯
∞→
.
Định nghĩa 4: 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.2.3. Không gian véc tơ tôpô.
1.2.3.1. Định nghĩa: Cho
X
là một tập bất kỳ. Ta nói họ
τ
các tập con của
X
là một tôpô trên
X
nếu:
(i)
τ
φ
∈X, .

(ii) Giao của một số hữu hạn các tập thuộc
τ
thì cũng thuộc
τ
(
τ
kín
đối với phép giao hữu hạn).
(iii) Hợp của một số tuỳ ý các tập thuộc
τ
thì cũng thuộc
τ
(
τ
kín
đối với phép hợp bất kỳ).
Tập
X
cùng với một tôpô
τ
trên
X
gọi là không gian tôpô
),(
τ
X
. Các
tập thuộc họ
τ
đều là các tập mở.

14
1.2.3.2. Định nghĩa và định lí.
Định nghĩa 1: Tập
V
trong không gian tôpô
),(
τ
X
gọi là một lân cận của
điểm
Xx ∈ nếu VDxD ⊆
∈
∈∃ :
τ
.
Định lí: Ký hiệu
x
ϑ
là họ tất cả các lân cận của Xx
∈
. Ta có
φ
ϑ
≠
x
và:
i)
VxV

x
∈
⇒∈
ϑ
.
(ii)
xx
UUVV
ϑ
ϑ
∈
⇒⊆∈ , .
(iii)
xx
VUVU
ϑ
ϑ
∈
∩
⇒∈, .
(4i)
WyWWxVWV
yx
∈
∀
∈
∈
⊆
∃
⇒∈ ,,:

ϑ
ϑ
.
Định nghĩa 2: Họ con
xx
U
ϑ
⊂
gọi là hệ cơ sở các lân cận của
x
nếu:
VUUUV
xx
⊆
∈
∃
∈
∀
:,
ϑ

Định nghĩa 3: Không gian tôpô
),(
τ
X gọi là tách Hausdorff nếu với hai điểm
bất kỳ
yxXyx ≠∈ ,, tồn tại các lân cận tương ứng
yx
VV , của
x

và y sao cho
φ
=∩
yx
VV
.
1.2.3.3. Định nghĩa.
Định nghĩa 1: Giả sử
X
là không gian véc tơ trên trường
K
. Ta nói tôpô
τ

trên không gian véc tơ
X
là tương thích với cấu trúc đại số nếu các phép toán
đại số trong
X
liên tục trong tôpô đó, nghĩa là:
(i)
y
x
+ là hàm liên tục của hai biến
y
x
,
(tức là mọi lân cận V của
điểm
y

x
+ đều có một lân cận
x
U của điểm
x
và một lân cận
y
U của điểm y
sao cho nếu
VyxUyUx
yx
∈
+
⇒∈∈
1111
, ).
(ii)
x
α
là hàm liên tục của hai biến
x
,
α
(tức là mọi lân cận V của
x
α
đều có một số 0>
ε
và một lân cận U của
x

sao cho
VxUx ∈⇒∈<−
1111
,
αεαα
).
Định nghĩa 2: Không gian véc tơ
X
trên đó có một tôpô tương thích với cấu

15
trúc đại số gọi là không gian véc tơ tôpô.
Định nghĩa 3: Tập con
A
trong không gian véc tơ
X
gọi là tập lồi nếu:
AyxAyx
∈
−
+
⇒
≤
≤
∈
∀ )1(10;,
α
α
α
.

Định nghĩa 4: Không gian véc tơ tôpô
X
gọi là không gian lồi địa phương
nếu mọi phần tử của
X
đều có cơ sở lân cận thành lập từ các tập lồi (hay nếu
phần tử
X∈
θ
có cơ sở lân cận thành lập từ các tập lồi).
1.2.3.4. Ánh xạ afin.
Định nghĩa 1: Cho
YX ,
là hai không gian véc tơ, ánh xạ RYX →×:.,. gọi
là một dạng song tuyến tính tách được theo từng biến nếu:

RYyyXxyxyxyyx ∈∀∈∀∈∀+=+
µλµλµλ
,,,,;,,,
212121


RYyXxxyxyxyxx ∈∀∈∀∈∀+=+
µλµλµλ
,,,,;,,,
212121


{
}

0,:,
00
≠∈∃−∈∀ yxYyXx
X
θ


{
}
0,:,
00
≠∈∃−∈∀ yxXxYy
Y
θ

Định nghĩa 2: Cho
X
là không gian véc tơ tô pô lồi địa phương tách
Hausdorff. Ta viết:
####
,;)(, XxXxxxxx ∈∈=
(đây là một dạng song tuyến tính trên
#
X
X
×
).
Ánh xạ affine trên
X
là hàm có dạng:

αϕ
+=
#
,)( xxx với RXxXx ∈∈∈
α
,,
##


ϕ
liên tục khi và chỉ khi
*#
Xx ∈ .
1.3. Đạo hàm Gateaux và đạo hàm Frechet.
Định nghĩa 1: Cho YX , là hai không gian định chuẩn, U là một tập con mở
của
X
, ánh xạ YUF →: . Khi đó:
(i) Ánh xạ tuyến tính liên tục
YXT →: là đạo hàm Gateaux của
F
tại
Ux ∈
0
nếu và chỉ nếu:

16
0)()(
1

:
00
0
=−−+∈∀
→
TuxFuxFmilXu
ττ
τ
τ

(ii) Ánh xạ tuyến tính liên tục
YXT →:
là đạo hàm Frechet của
F

tại
Ux ∈
0
nếu và chỉ nếu:

0)()(
1
:
00
0
=−−+∈∀
→
TuxFuxF
u
milXu

u

(iii) Ánh xạ tuyến tính liên tục
YXT →: là đạo hàm Gateaux yếu
của
F tại Ux ∈
0
nếu và chỉ nếu:
[
]
0)()(
1
:,
00*
0
**
=−−+∈∀∈∀
→
TuxFuxFxmilXxXu
ττ
τ
τ

(4i) Ánh xạ tuyến tính liên tục
YXT →: là đạo hàm Frechet yếu của
F
tại Ux ∈
0
nếu và chỉ nếu:

[
]
0)()(
1
:
00*
0
**
=−−+∈∀
→
TuxFuxFx
u
milXx
u

Các đạo hàm được định nghĩa trên đây đều được kí hiệu là
)('
0
xFT = .
Định nghĩa 2: Cho
,2;,, ,
1
≥nYXX
n
là các không gian định chuẩn, ánh xạ:
YXXF
n
→×× :
1
, với mọi

nn
XXxxx ××∈= ), ,(
11
, ta cố định
n
n
XXxxx ××∈
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⎟
⎠
⎞
⎜
⎝
⎛
⎟
⎠
⎞
⎜
⎝
⎛
= , ,
1
0

1
00

xét ánh xạ:
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⎟
⎠
⎞
⎜
⎝
⎛
⎟
⎠
⎞
⎜
⎝
⎛
⎟
⎠
⎞
⎜
⎝
⎛

⎟
⎠
⎞
⎜
⎝
⎛
==→
+− ni
i
i
iiii
xxxxxFxFniYXF
0
1
0
1
0
1
0
, ,,,, ,)(,,1;:
Nếu
i
F
có đạo hàm Frechet tại điểm
i
x
⎟
⎠
⎞
⎜

⎝
⎛
0
thì đạo hàm đó gọi là đạo
hàm riêng Frechet của
F
theo
i
x tại điểm
0
x . Ta viết:
i
i
i
xFx
x
F
⎟
⎠
⎞
⎜
⎝
⎛
=
⎟
⎠
⎞
⎜
⎝
⎛

∂
∂
00
')(

Khi
RYXX
n
====
1
thì đạo hàm riêng Frechet này trùng với đạo
hàm riêng thông thường.

17
Ví dụ 1: Cho
YX ,
là hai không gian định chuẩn,
YXF →:
là ánh xạ hằng
(tức là
YbconstbFxXx
∈
=
=∈∀ ,:
). Khi đó
∧
=⇒∈∀
θ
xFXx ' (ánh xạ
không).

Ví dụ 2:
RRf →: là hàm số thực, với Rx ∈
0
ta có )('
0
xf là đạo hàm của f
tại
0
x theo nghĩa thông thường.
Ví dụ 3: Giả sử ánh xạ:
nxxFxFxxx
RRF
nn
n
2), ,(), ,(
:
11
≥==
→
a

có các đạo hàm riêng theo
n
xx , ,
1
liên tục. Khi đó mọi
n
n
Rxxx ∈= ), ,(
1

ta có:
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
∂
∂
∂
∂
= )(), ,()('
1
x
x
F
x
x
F
xF
n

và
∑
=
∂
∂
=⇒∈=∀

n
j
j
j
n
n
h
x
F
hxFRhhh
1
1
)('), ,(

Ví dụ 4: Nếu mỗi ánh xạ
)2(.,1,2;:), ,()(
1
≥=≥⎯→⎯= mminRRxxfxf
n
nii

có đạo hàm tại
x
thì ánh xạ:
()
)(), ,(
:), ,(
1
1
xfxfxF

RRffF
m
mn
m
=
⎯→⎯=

là hàm véc tơ nhiều biến có đạo hàm tại
x
và:

⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎜
⎜
⎜

⎜
⎜
⎜
⎜
⎜
⎝
⎛
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
∂
=
⎟
⎟
⎠
⎞
⎜

⎜
⎝
⎛
∂
∂
=
⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
=
×
n
mmm
n
n
nm
j
i
m
x
f
x
f

x
f
x
f
x
f
x
f
x
f
x
f
x
f
x
x
f
xf
xf
xF
K
MMMM
K
L
M
21
2
2
2
1

2
1
2
1
1
1
1
)(
)('
)('
)('

(ma trận Jacobi của ánh xạ
F
).

18

⎟
⎟
⎟
⎠
⎞
⎜
⎜
⎜
⎝
⎛
=⇒∈=∀
hxf

hxf
hxFRhhh
m
n
n
)('
)('
)('), ,(
1
1
M
Định lí hàm số ngược: Cho
YX , là hai không gian định chuẩn, U là tập con
mở trong
X
, ánh xạ
YUF →:
là phép đồng phôi từ
U
vào tập con mở
YUFV ⊂= )( . Giả sử
F
có đạo hàm Frechet tại điểm Ux ∈
0
và YXxF →:)('
0

là phép đồng phôi tuyến tính. Khi đó ánh xạ ngược
XVF →
−

:
1
có đạo hàm
Frechet tại điểm
VxFy ∈=
00
và
100101
)(')(')()(')(
−−−
== xFFxFyF


















19

CHƯƠNG 2: MỘT SỐ PHƯƠNG PHÁP LẶP GIẢI PHƯƠNG TRÌNH
VÀ HỆ PHƯƠNG TRÌNH PHI TUYẾN

2.1. Phương pháp Newton và một số biến thể của nó.
2.1.1. Phương pháp dây cung song song.
Cho
RRUf →⊂: là hàm số thực một biến có một nghiệm
*
x
, ta thay
giá trị của
f
tại một xấp xỉ
0
x
của
*
x
bởi một hàm tuyến tính:
)()()(
00
xfxxxl +−=
α

với độ lệch
0≠
α
thích hợp, sau đó lấy nghiệm
1
x của l làm một xấp xỉ mới

của
*
x
. Lặp lại cách làm này với
α
cố định, ta có phương pháp lặp:

, 1,0),(
11
=−=
−+
kxfxx
kkk
α
(2.1.1)
gọi là phương pháp dây cung song song trong không gian một chiều và
được biểu diễn hình học như hình 2.1 dưới đây:

Hình 2.1
Ta mở rộng từ hàm
f đến hàm
nn
RRDF →⊂: trong không gian n chiều,
2≥n , giả sử
F
có nghiệm
*
x , có ma trận Jacobi ),( RnGLA
∈
không đổi, ta

thay giá trị của
F
tại một xấp xỉ
k
x của
*
x bởi hàm afin:
kk
k
xFxxAxL +−= )(

20
, đặt
1+k
x
là nghiệm duy nhất của
θ
=
xL
k
(với
, 1,0;)0, ,0( =∈= kR
n
θ
).
Ta có phương pháp dây cung song song trong không gian
n chiều

, 1,0,
11

=−=
−+
kxFAxx
kkk
(2.1.2)
Về mặt hình học,
1+k
x chính là giao điểm của n siêu phẳng:

., ,1;0)()(
1
nixfxxa
k
i
n
j
k
jjji
==+−
∑
=
với siêu phẳng
θ
=x
trong
1+n
R
.
Để áp dụng phương pháp dây cung song song thì điều quan trọng là
phải chọn ma trận

A
thích hợp. Ta giới thiệu một vài cách chọn sau đây:
Cách 1:
IA
λ
= với
*
R∈
λ
nào đó. Từ đó ta thấy phương pháp lặp trong
không gian một chiều (2.1.1) được gắn vào mỗi thành phần riêng
i
f
của
F
.
Cách 2: Đối với quá trình lặp trong không gian một chiều (2.1.1), ta có
)('
0
xf là hệ số góc của tiếp tuyến của f tại
0
x , ta chọn )('
0
xf=
α
, cách chọn
này hoàn toàn hợp lý.
Đối với quá trình lặp trong không gian
n
chiều (2.1.2), ta chọn

)('
0
xFA = , với )(' xF là đạo hàm Gato của
F
tại
x
, khi đó (2.1.2) là phương
pháp Newton cải biên:

, 1,0,)('
101
=−=
−+
kxFxFxx
kkk
(2.1.3)
Cách 3: Chọn
A
một cách ngẫu nhiên tuỳ thuộc vào dạng đặc biệt nào đó của
F
.
Ví dụ:

xGxAxF
−
=
(2.1.4)
với
),( RnGLA∈ nào đó, G là một hàm phi tuyến, ta có phương pháp Picard:

, 1,0,
11
==
−+
kxGAx
kk
(2.1.5)
Vì
xFAxxGA
11 −−
−= nên rõ ràng (2.1.5) có thể được viết lại dưới dạng (2.1.2)
Chú ý: Có nhiều cách chọn ma trận
A
khác nữa trong phương pháp dây

21
cung song song . Tuy nhiên, yêu cầu cơ bản là phải chọn
A
sao cho quá trình
lặp (2.1.2) ít nhất phải là hội tụ địa phương. Nghĩa là khi
0
x đủ gần tới
nghiệm
*
x của 0=xF thì chắc chắn ta có:
*
xxmil
k
k
=

∞→
.
Ta có thể chứng minh được rằng khi
)('
*
xF tồn tại thì điều kiện cần và
đủ để (2.1.2) hội tụ địa phương là:

(
)
1)('
*1
<−=
−
xFAI
ρσ
(2.1.6)
Trong đó
ρ
là bán kính phổ của ma trận. Hơn nữa, ta còn có thể chứng minh
được rằng
σ
càng nhỏ thì hội tụ càng nhanh.
Do
*
x
chưa biết nên khó chọn trước được ma trận
A
để có (2.1.6),
trong đó cách chọn lý tưởng vẫn là

)('
*
xFA = . Từ đó buộc ta phải đi xét các
quá trình lặp:
, 1,0,
11
=−=
−+
kxFAxx
k
k
kk

với
)('
k
k
xFA =
được thay đổi từ bước nọ đến bước kia sao cho:
)('
*
xFAmil
k
k
=
∞→

Nguyên mẫu của các phương pháp như thế gọi là phương pháp
Newton. Phương pháp Newton đóng vai trò trọng tâm trong tất cả các phương
pháp lặp giải phương trình và hệ phương trình phi tuyến n ẩn số.

2.1.2. Phương pháp Newton.
Cho
RRUf →⊂: là hàm số thực một biến có một nghiệm
*
x , quá
trình lặp:

, 1,0),()('
11
=−=
−+
kxfxfxx
kkkk
(2.1.7)
là phương pháp Newton trong không gian một chiều và được biểu diễn hình
học như hình 2.2 sau đây:

22

Hình 2.2
Ta mở rộng từ hàm
f
đến hàm
nn
RRDF →⊂: trong không gian n
chiều,
2≥n , có nghiệm
*
x , bằng cách thay )('
k

xf trong (2.1.7) bằng đạo hàm
Gato
)('
k
xF của
F
tại
k
x , ta có phương pháp Newton trong không gian n
chiều:

, 1,0,)('
11
=−=
−+
kxFxFxx
kkkk
(2.1.8)
Bước lặp từ
k
x đến
1+k
x ở (2.1.8) được mô tả hình học là: Mỗi thành
phần
i
f của
F
được xấp xỉ bởi hàm afin:

)())((')(

k
i
kk
i
xfxxxfxL +−=
(2.1.9)
mà hàm afin này mô tả siêu phẳng tiếp xúc với
i
f tại
k
x , sau đó lấy
1+k
x là
giao điểm của
n siêu phẳng (2.1.9) với siêu phẳng
θ
=
x trong
1+n
R
.
Khi
1=n
thì (2.1.8) chính là (2.1.7). Vậy (2.1.8) là công thức tổng quát
của phương pháp Newton từ một chiều đến
n chiều.
Chú ý: Có nhiều phương pháp lặp trong không gian
n
chiều khác nữa
mà cũng có điều kiện khi cho

1
=
n thì nó trở thành (2.1.7), chẳng hạn ta xét
quá trình lặp:

, 1,0;)1()('
11
=−+−=
−+
kxGnxFxFxx
kkkkk
(2.1.10)

23
với
nn
RRG →: là ánh xạ tuỳ ý nào đó.
Tầm quan trọng của phép lặp Newton (2.1.8) dựa trên cơ sở là: từ
những điều kiện đã cho đối với
F , bất đẳng thức dạng:

2
**1
xxcxx
kk
−≤−
+
(2.1.11)
luôn thoả mãn, miễn là các
k

x đủ gần với nghiệm
*
x . Điều này chứng tỏ rằng
sai số thứ
1
+
k
tỉ lệ thuận với bình phương của sai số thứ
k
, bởi vậy sự hội tụ
là rất nhanh mỗi khi các sai số là rất nhỏ. Đây còn được gọi là tính chất của sự
hội tụ bậc hai; nói chung tính chất này không đúng với các quá trình lặp dạng
(2.1.10); bất đẳng thức (2.1.11) kết hợp với (2.1.8) tạo ra phương pháp
Newton - một tâm điểm trong việc nghiên cứu các phương pháp lặp giải
phương trình và hệ phương trình phi tuyến n ẩn số.
Ta hoàn toàn có thể ch
ứng minh được (2.1.11). Hơn nữa, ta có thể thấy
rằng bất đẳng thức (2.1.11) có được là do phép lấy đạo hàm sau đây của
phương pháp Newton:
Nếu
F
có đạo hàm Frechet tại
k
x
thì:

)())(('
*** kkkk
xxRxxxFxFxF −+−+==
θ

(2.1.12)
ở đây:
0
)(
0
=
→
h
hR
mil
h

Từ đó nếu
k
x gần với nghiệm
*
x thì hiển nhiên số hạng )(
* k
xxR −
giảm đi và xấp xỉ hiệu
k
xx −
*
bằng nghiệm h của hệ tuyến tính:

kk
xFhxF −=)(' (2.1.13)
Hay nói cách khác,
kkkkk
xFxFxhxx

11
)('
−+
−=+=
được lấy là một xấp xỉ mới của
*
x .
Nếu
''F bị chặn trong một lân cận của
*
x thì ta có:
2
**
)( xxxxR
kk
−≤−
α


24
và giả sử )('
*
xF không suy biến từ đó dẫn tới đánh giá (2.1.11).
Mọi phép lấy đạo hàm có liên quan của phương pháp Newton hoàn
toàn có thể thực hiện được dựa vào định lí hàm số ngược.
Thực chất là nếu
F khả vi liên tục và )(' xF không suy biến trong một
lân cận U của
*
x

thì theo định lý hàm số ngược ta có: thu hẹp
U
F
của
F
trên
U có một nghịch đảo
1−
=
U
FG
khả vi và
1
)(')('
−
= xFxFG
.
Do đó nếu ta khai triển
G
theo
k
xF , ta có:
)()(')()(')()(
1* kkkkkkkk
xFRxFxFxxFRxFxFGxFGGx
∧
−
∧
+−=+−==
θ

(2.1.14)
Ta thấy (2.1.14) lại cho ta phương pháp Newton khi bỏ đi số hạng cuối.
2.1.3. Một số biến thể của phương pháp Newton.
Mặc dù phương pháp Newton là lý thuyết có sức thu hút người đọc,
nhưng nó rất khó khăn trong thực hành. Thực tế tại mỗi bước ta cần phải có
nghiệm của hệ tuyến tính (2.1.13), hiếm khi ta tính chính xác được giá trị
1
)('
−k
xF , đặc biệt là các bài toán nảy sinh từ các phương trình đạo hàm riêng
trong đó kích cỡ của hệ có thể là vài nghìn; đây là một việc hết sức khó khăn.
Hơn nữa, tại mỗi bước không chỉ cần tới
n thành phần của
k
xF mà còn cần
có cả
2
n dạng )('
k
xF , và nếu các đạo hàm riêng )(
k
ij
xf∂ có dạng hàm số phức
tạp thì ta hoàn toàn không cần tính toán các đạo hàm riêng này. Hầu hết các
phương pháp của chương này đều tránh việc tính toán một cách chính xác các
đạo hàm. Hai phép xấp xỉ khác nhau điển hình thường được sử dụng để tránh
tính
)(' xF đơn giản nhất là xấp xỉ các đạo hàm riêng .,1,);( njixf
ij
=∂ bằng

các tỷ sai phân:

⎥
⎦
⎤
⎢
⎣
⎡
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
+−
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
+=∂
∑∑
−
==
1

11
1
)(
j
k
k
kii
j
k
k
kii
ji
ij
ehxfehxf
h
xf
(2.1.15)
và

25

[
]
)()(
1
)( xfehxf
h
xf
i
j

jii
ji
ij
−+=∂
(2.1.16)
với
Rh
ji
∈ là các tham số phân biệt, )0, ,0,1,0, ,0(=
j
e là véc tơ toạ độ thứ
j
.
Biến thể 1: Tổng quát, cho
p
Rh ∈
là một véc tơ tham số,
),( hx
ji
∆
là các xấp
xỉ khác nhau của
)(xf
ij
∂
với tính chất: mỗi khi
)(xf
ij
∂
tồn tại thì với

.,1, nji =

ta có:

)(),(
0
xfhxmil
ijji
h
∂
=
∆
→
(2.1.17)
Khi đó, với ma trận sai số:

(
)
),(),( hxhxJ
ji
∆
=
(2.1.18)
thì quá trình lặp:

, 1,0;),(
11
=−=
−+
kxFhxJxx

kkkkk
(2.1.19)
được gọi là phương pháp Newton rời rạc với các véc tơ tham số
pk
Rh ∈

được phép thay đổi theo chỉ số lần lặp.
Chọn
k
h đơn giản nhất là
pk
Rhh ∈= (các tham số khác nhau được giữ
cố định cùng bằng một hằng số trong suốt quá trình lặp).
Ta dễ thấy phép lặp chỉ có một tốc độ hội tụ tuyến tính, và để tiến dần
đến sự hội tụ nhanh của phương pháp Newton thì cần phải có
0=
∞→
k
k
hmil .
Muốn vậy ta có nhiều cách chọn
k
h
; chẳng hạn chọn
hh
k
k
γ
=
, với

p
Rh ∈
cố
định và dãy
{}
R
k
⊂
γ
sao cho 0
=
∞→
k
k
mil
γ

Đặc biệt là các phương pháp thú vị khác sinh ra khi lấy
k
h là các hàm
số.
Thực chất mỗi phương pháp được nói tới trong chương này là một biến
thể của phương pháp Newton theo một nghĩa nào đó.
Một trong những quy tắc được gắn với bất kỳ phương pháp lặp nào là
quy tắc “làm giảm chuẩn” theo nghĩa: