Các phương pháp hàm Spline và một số ứng dụng

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



LÃ THỊ NGỌ


CÁC PHƯƠNG PHÁP HÀM SPLINE VÀ
MỘT SỐ ỨNG DỤNG

Chuyên ngành: Toán giải tích
Mã số: 60 46 01 02

LUẬN VĂN THẠC SĨ TOÁN HỌC

Người hướng dẫn khoa học: TS. NGUYỄN VĂN TUẤN




HÀ NỘI, 2015

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. Nguyễn Văn Tuấn
.
Tác giả xin bày tỏ sự kính trọng, lòng biết ơn sâu sắc tới TS. Nguyễn Văn Tuấn,
người thầy đã luôn tận tình hướng dẫn, chỉ bảo, động viên và khuyến khích tác giả
trong học tập, nghiên cứu và hoàn thành luận văn. Đồng thời, tác giả cũng xin được

gửi lời cảm ơn chân thành tới các thầy cô giáo trong Khoa Toán, trường Đại hoc sư
phạm Hà Nội 2 cùng với các thầy cô tham gia giảng dạy cao học chuyên ngành Toán
giải tích đã giúp đỡ tác giả trong suốt thời gian họ c tập. Nhân dịp này tác giả cũng
xin được gửi lời cảm ơn chân thành tới gia đình, bạn bè và đồng nghiệp trong trường
THPT Đa Phúc - Sóc Sơn - Hà Nội ( nơi tác giả đang công tác) đã luôn cổ vũ, động
viên, tạo điều kiện giúp đỡ tác giả trong suốt quá trình học tập và thực hiện luận văn
này.
Hà Nội, ngày 05 tháng 12 năm 2014
Tác giả
Lã Thị Ngọ
LỜI CAM ĐOAN
Tôi xin cam đoan luận văn này là công trình nghiên cứu của riêng tôi, được
hoàn thành dưới sự hướng dẫn của
TS. Nguyễn Văn Tuấn
. Trong khi nghiên cứu
luận văn tôi đã kế thừa thành quả khoa học của các nhà khoa học với sự trân trọng
biết ơn.
Hà Nội, ngày 05 tháng 12 năm 2014
Tác giả
Lã Thị Ngọ
Mục lục
Lời cảm ơn i
Lời cam đoan ii
Bảng ký hiệu v
Mở đầu 1
1 KIẾN THỨC CHUẨN BỊ 3
1.1 Không gian tuyến tính . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.1.1 Định nghĩa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.1.2 Vectơ độc lập tuyến tính, không gian con . . . . . . . . . . . . . 4
1.2 Không gian định chuẩn . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

1.2.1 Khái niệm không gian định chuẩn . . . . . . . . . . . . . . . . . 6
1.2.2 Sự hội tụ trong không gian định chuẩn . . . . . . . . . . . . . . 7
1.2.3 Sai số . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.2.4 Xấp xỉ tốt nhất . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2.5 Tốc độ hội tụ của nghiệm xấp xỉ . . . . . . . . . . . . . . . . . 9
1.2.6 Ma trận đường chéo trội . . . . . . . . . . . . . . . . . . . . . . 10
1.3 Không gian Hilbert . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.3.1 Mở đầu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.3.2 Định nghĩa không gian Hilbert . . . . . . . . . . . . . . . . . . . 11
1.3.3 Định nghĩa hệ trực chuẩn . . . . . . . . . . . . . . . . . . . . . 11
1.3.4 Một số kết quả cơ bản . . . . . . . . . . . . . . . . . . . . . . . 13
iii
2 MỘT SỐ VẤN ĐỀ VỀ SPLINE VÀ B-SPLINE 15
2.1 Mở đầu về hàm spline và B-spline . . . . . . . . . . . . . . . . . . . . . 15
2.1.1 Tổ hợp lồi và bao lồi . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.2 Các khái niệm cơ bản . . . . . . . . . . . . . . . . . . . . . . . . 17
2.1.3 Nội suy các đường cong đa thức . . . . . . . . . . . . . . . . . . 18
2.1.4 Xây dựng đường cong spline . . . . . . . . . . . . . . . . . . . . 19
2.1.5 Giới thiệu đường cong spline theo các hàm số cơ bản . . . . . . 23
2.2 Tính chất cơ bản của spline và B-spline . . . . . . . . . . . . . . . . . . 25
2.2.1 Một vài hệ quả đơn giản của hệ thức truy hồi . . . . . . . . . . 25
2.2.2 Các tổ hợp tuyến tính của B - spline . . . . . . . . . . . . . . . 29
2.2.3 Ma trận biểu diễn của các B - spline . . . . . . . . . . . . . . . 32
2.2.4 Thuật toán để ước lượng một spline . . . . . . . . . . . . . . . 35
3 NỘI SUY HÀM SPLINE 36
3.1 Các phương pháp xấp xỉ địa phương . . . . . . . . . . . . . . . . . . . 36
3.1.1 Nội suy tuyến tính . . . . . . . . . . . . . . . . . . . . . . . . . 36
3.1.2 Nội suy bậc ba Hermite . . . . . . . . . . . . . . . . . . . . . . 37
3.1.3 Ước lượng các đạo hàm . . . . . . . . . . . . . . . . . . . . . . . 39
3.2 Nội suy spline bậc ba . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40

3.3 Phép tính xấp xỉ spline tổng quát . . . . . . . . . . . . . . . . . . . . . 42
3.4 Áp dụng . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3.4.1 Nội suy bằng hàm spline bậc 1 . . . . . . . . . . . . . . . . . . 43
3.4.2 Nội suy bằng hàm spline bậc 2 . . . . . . . . . . . . . . . . . . 53
Kết luận 57
Tài liệu tham khảo 58
iv
BẢNG KÝ HIỆU
N Tập số tự nhiên
Z Tập các số nguyên
R Tập số thực
C Tập số phức
S
3
Không gian các hàm Spline bậc 3
C
[a;b]
Không gian các hàm số thực liên tục trên đoạn [a; b]
MỞ ĐẦU
1. Lý do chọn đề tài
Trong khoa học, kĩ thuật chúng ta thường gặp rất nhiều bài toán cần phải tính
được giá trị của hàm số tại một điểm nhưng để tính đúng giá trị của hàm số tại một
điểm của một số hàm gặp rất nhiều khó khăn ví dụ như hàm số mũ, hàm lượng giác,
hàm số logarit, Để giải quyết vấn đề này người ta đã nghiên cứu nhiều phương pháp
khác nhau. Trong đó phương pháp hàm Spline đang được nhiều nhà toán học trong
và ngoài nước quan tâm nghiên cứu. Áp dụng hàm spline và phương pháp nội suy để
xấp xỉ hàm số người ta chia khoảng xác định thành nhiều đoạn, trên mỗi đoạn ta xấp
xỉ bằng một hàm spline, từ đó ta xấp xỉ được hàm số đã cho.Tính xấp xỉ giá trị của
hàm số tại một điểm bằng phương pháp hàm Spline rất thuận lợi vì nó là những hàm
đa thức nên việc tính toán, lập trình với hàm đa thức rất thuận tiện và dễ dàng.

Do vậy, với sự hướng dẫn của TS. Nguyễn Văn Tuấn, tôi đã nghiên cứu luận văn: “
Các phương pháp hàm Spline và một số ứng dụng”.
Bố cục của luận văn gồm 3 chương:
Chương 1 của luận văn trình bày một số khái niệm cơ bản để sử dụng cho các chương
sau.
Chương 2 của luận văn trình bày về các khái niệm và tính chất của hàm spline và
B-spline.
Chương 3 của luận văn trình bày về nội suy hàm spline và ứng dụng phần mềm Maple
vào xấp xỉ hàm số bằng hàm spline.
2. Mục đích nghiên cứu
Nghiên cứu để nắm được một số phương pháp hàm spline.
Ứng dụng hàm Spline để tính giá trị của hàm số tại một điểm và một số ứng dụng
khác.
3. Nhiệm vụ nghiên cứu
Nghiên cứu các khái niệm, các tính chất của hàm Spline và B- spline.
Xấp xỉ hàm số bằng hàm spline.
Ứng dụng phần mềm Maple vào xấp xỉ hàm số bằng hàm spline.
2
4. Đối tượng và phạm vi nghiên cứu
Đối tượng nghiên cứu: Các hàm spline, phương pháp spline
Phạm vi nghiên cứu: Các khái niệm, tính chất, ứng dụng vào xấp xỉ hàm số.
5. Phương pháp nghiên cứu
Nghiên cứu tài liệu, phân tích, tổng hợp, tham khảo ý kiến chuyên gia.
6. Giả thuyết khoa họ c
Áp dụng phương pháp spline để xấp xỉ một lớp hàm có nhiều ứng dụng trong
thực tế.
Chương 1
KIẾN THỨC CHUẨN BỊ
1.1 Không gian tuyến tính
1.1.1 Định nghĩa

Định nghĩa 1.1.1. Cho tập hợp X khác rỗng cùng với một phép toán hai ngôi viết
theo lối cộng (+) và một ánh xạ ψ : K × X → X. Với mỗi α ∈ K và mỗi x ∈ X thì
phần tử ψ(α, x) được gọi là tích của số α với phần tử x và được kí hiệu αx. Giả sử các
điều kiện sau được thỏa mãn:
1) x + y = y + x, ∀x, y ∈ X;
2) (x + y) + z = x + (y + z), ∀x, y, z ∈ X;
3) Tồn tại duy nhất phần tử θ sao cho x + θ = θ + x, ∀x ∈ X (phần tử này gọi là
phần tử không);
4) Ứng với mỗi phần tử x ∈ X, tồn tại duy nhất phần tử −x ∈ X sao cho x+(−x) = θ
phần tử (−x) được gọi là phần tử đối của x;
5) 1 · x = x, ∀x ∈ X;
6) α(βx) = (αβ)x, ∀α, β ∈ K, ∀x ∈ X;
7) (α + β)x = αx + βx, ∀x ∈ X, α, β ∈ K;
8) α(x + y) = αx + αy, ∀x, y ∈ X, α ∈ K;
Khi đó ta nói rằng X là một không gian tuyến tính trên trường K, K là trường số thực
hoặc trường số phức và mỗi phần tử ∀x ∈ X được gọi là một vectơ; còn các điều kiện
trên gọi là các tiên đề về không gian tuyến tính.
3
4
Ví dụ 1.1.1. Trong mặt phẳng thực R
2
Tập X = R
2
là tập
R
2
= {(x
1
, x
2

) : x
1
và x
2
là các số thực}
Với mỗi số thực α và các vectơ x = (x
1
, x
2
), y = (y
1
, y
2
) ∈ X, phép cộng và
nhân vô hướng được định nghĩa:
x + y = (x
1
+ y
1
, x
2
+ y
2
)
αx = (αx
1
, αx
2
)
là không gian tuyến tính.

Ví dụ 1.1.2. Không gian C[a, b]
Không gian C[a, b] là không gian các hàm liên tục trên [a, b].
Với mỗi số thực α và f(t), g(t) ∈ C[a, b], phép cộng và nhân vô hướng được định
nghĩa:
(f + g)(t) = f(t) + g(t), a ≤ t ≤ b
(αf)(t) = αf(t)
là không gian tuyến tính.
Ví dụ 1.1.3. P
n
[a, b], không gian các đa thức bậc n trên đoạn [a, b] là không gian tuyến
tính.
1.1.2 Vectơ độc lập tuyến tính, không gian con
Định nghĩa 1.1.2. Cho X là một không gian tuyến tính.
Một tổ hợp tuyến tính của các vectơ x
1
, x
2
, . . . , x
n
∈ X là một tổng có dạng:
α
1
x
1
+ α
2
x
2
+ . . . + α
n

x
n
Các vectơ x
1
, x
2
, . . . , x
n
được gọi là độc lập tuyến tính nếu
α
1
x
1
+ α
2
x
2
+ . . . + α
n
x
n
= 0 ⇒ α
1
= α
2
= . . . = α
n
== θ
Các vectơ x
1

, x
2
, . . . , x
n
được gọi là phụ thuộc tuyến tính nếu chúng không độc
lập tuyến tính, tức là tồn tại những số α
1
, α
2
, . . . , α
n
trong đó có ít nhất một số khác
0, sao cho:
α
1
x
1
+ α
2
x
2
+ . . . + α
n
x
n
= 0
5
Ví dụ, hai vectơ x và (−x) là phụ thuộc tuyến tính vì :
1 · x + 1 · (−x) = 0
Nếu trong các vectơ x

1
, x
2
, . . . , x
n
có một vectơ bằng 0 thì chúng là phụ thuộc
tuyến tính.
Một không gian tuyến tính X được gọi là không gian k chiều nếu trong X có k
vectơ độc lập tuyến tính và không có k + 1 vectơ độc lập tuyến tính.
Trong trường hợp này một tập k vectơ độc lập tuyến tính của X gọi là một cơ
sở của nó.
Các không gian k chiều, với k là một số nguyên không âm bất kì gọi là không
gian hữu hạn chiều.
Một không gian vô hạn chiều tức là sao cho với mọi k đều tìm được k vectơ độc
lập tuyến tính của nó.
Ví dụ 1.1.4. R
k
= {(a
1
, a
2
, ··· , a
k
) |a
i
∈ R} là không gian k chiều, với cơ sở là:
x
1
= (1, 0, . . . , 0), x
2

= (0, 1, 0, . . . , 0), . . . , x
n
= (0, 0, . . . , 1)
Không gian R[a, b] là vô hạn chiều vì với mọi k ta luôn có k phần tử của nó độc lập
tuyến tính, đó là:
x
1
(t) = t, x
2
(t) = t
2
, . . . , x
n
(t) = t
n
Nếu X là không gian k chiều và x
1
, x
2
, . . . , x
k
là một cơ sở của nó thì mọi x
thuộc X đều được biểu diễn duy nhất dưới dạng
x = α
1
x
1
+ α
2
x

2
+ . . . + α
k
x
k
Các số α
1
, α
2
, . . . , α
k
là tọa độ của vectơ x đối với cơ sở x
1
, x
2
, . . . , x
k
. Nếu ta
làm phép ánh xạ 1 − 1 : x ↔ (α
1
, α
2
, . . . , α
k
) thì đó là một phép đẳng cấu giữa X và
R
k
. Như vậy không gian tuyến tính k chiều bao giờ cũng đẳng cấu với không gian R
k
.

Định nghĩa 1.1.3. Một tập con không rỗng M của một không gian X gọi là một
không gian con, nếu nó kín đối với phép cộng phần tử với phần tử và phép nhân phần
tử với một số, nghĩa là:
1) ∀x, y ∈ M ⇒ x + y ∈ M.
2) ∀x ∈ M, α ∈ R ⇒ αx ∈ M.
Cho A là một tập con bất kì, khác rỗng của X. Khi đó, bao giờ cũng có ít nhất
một không gian con bao hàm A, đó chính là X. Vậy họ các không gian con bao hàm
6
A khác rỗng, giao của họ các không gian ấy cũng là một không gian con và là không
gian con nhỏ nhất bao hàm A. Không gian này gọi là không gian con sinh bởi tập A.
Như vậy, không gian con sinh bởi A là tập tất cả các tổ hợp tuyến tính hữu
hạn:
α
1
x
1
+ α
2
x
2
+ . . . + α
k
x
k
của những phần tử của A.
1.2 Không gian định chuẩn
1.2.1 Khái niệm không gian định chuẩn
Định nghĩa 1.2.1. Một không gian định chuẩn là một không gian tuyến tính X, trong
đó ứng với mỗi phần tử x ∈ X, ta có một số ∥x∥, gọi là chuẩn của nó, sao cho với mọi
x, y ∈ X, và mọi số α thỏa mãn 3 điều kiện sau:

1) ∥x∥ > 0 nếu x ̸= 0; ∥x∥ = 0 nếu x = 0.
2) ∥αx∥ = |α| · ∥x∥.
3) ∥x + y∥ ≤ ∥x∥ + ∥y∥ (bất đẳng thức tam giác).
Ví dụ 1.2.1. Không gian R
2
là không gian định chuẩn với chuẩn thường chọn là chuẩn:
∥x∥
2
=

x
2
1
+ x
2
2
Ngoài ra còn có những chuẩn khác chẳng hạn như:
∥x∥
1
= |x
1
| + |x
2
|
hay
∥x∥
∞
= max{|x
1
|, |x

2
|}
trong đó x = (x
1
, x
2
) ∈ R
2
.
Ví dụ 1.2.2. Không gian C[a, b] = {f : [a, b] → R |f liên tục trên [a, b]} là không gian
định chuẩn với chuẩn
∥f(t)∥ = max
a≤t≤b
|f(t)|
7
Định nghĩa 1.2.2. Cho không gian tuyến tính X với hai chuẩn ∥∥
1
và ∥∥
2
. Hai chuẩn
này được gọi là tương đương nếu tồn tại hằng số M > 0 và m > 0 sao cho:
m∥x
1
∥ ≤ ∥x
2
∥ ≤ M∥x
1
∥, ∀x ∈ X.
Trong ví dụ 1.3.1 thì cả 3 chuẩn đôi một tương đương. Chẳng hạn:
∥x

2
∥ ≤

2∥x∥
2
∞

1
2
=
√
2∥x∥
∞
.
Mặt khác:
∥x∥
∞
= max{|x
1
|, |x
2
|} ≤ (x
2
1
+ x
2
2
)
1
2

= ∥x∥
2
.
Do đó chọn M =
√
2, m = 1, ta có:
∥x∥
∞
≤ ∥x∥
2
≤
√
2∥x∥
∞
.
1.2.2 Sự hội tụ trong không gian định chuẩn
Giả sử X là không gian định chuẩn và {x
n
}
∞
n=1
⊂ X, x
0
∈ X.
1) x
n
−→ x
0
(dãy x
n

hội tụ tới x
0
) có nghĩa là ∥x
n
− x
0
∥ −→ 0.
2) Nếu x
n
−→ x
0
thì ∥x
n
∥ −→ ∥x
0
∥, tức là chuẩn ∥x
n
∥ là một hàm liên tục của x.
3) Mọi dãy hội tụ đều bị chặn, tức là: nếu x
n
hội tụ thì ∃M ∈ R, M > 0, ∀n, ∥x
n
∥ ≤
M.
4) Nếu x
n
−→ x
0
, y
n

−→ y
0
thì x
n
+ y
n
−→ x
0
+ y
0
.
5) Nếu x
n
−→ x
0
, α
n
−→ α
0
thì x
n
α
n
−→ x
0
α
0
, ∀{α
n
}

∞
n=1
⊂ R, α
0
∈ R.
6) Một dãy cơ bản trong không gian định chuẩn X là một dãy {x
n
} ⊂ X sao cho:
lim
m,n→∞
∥x
n
− x
m
∥ = 0.
Nếu trong không gian định chuẩn mọi dãy cơ bản đều hội tụ, tức là: ∥x
n
−x
m
∥ →
0 kéo theo sự tồn tại x
0
∈ X sao cho x
n
→ x
0
, thì không gian đó được gọi là không
gian đủ thường gọi là không gian Banach.
8
1.2.3 Sai số

Trong thực tế khi giải quyết các bài toán về kĩ thuật và vật lý ta thường không
biết chính xác giá trị của một đại lượng nào đó. Số liệu ban đầu mà ta có trong các
bài toán trên được gọi là số gần đúng. Nếu đánh giá được độ lệch của số gần đúng với
số đúng thì ta đánh giá được chất lượng của việc giải quyết bài toán. Do đó đi nghiên
cứu và đánh giá sự sai khác giữa số gần đúng và số đúng là yêu cầu bắt buộc trong
việc giải bài toán.
Định nghĩa 1.2.3. Số a được gọi là số gần đúng của số a
∗
nếu a sai khác với a
∗
không
nhiều.
Kí hiệu a ≈ a
∗
.
Định nghĩa 1.2.4. Đại lượng ∆ = |a − a
∗
| được gọi là sai số thực sự của a.
Nói chung ta không biết a
∗
nên ta không biết ∆. Tuy nhiên ta có thể ước lượng
sai số thực sự của a bằng số dương ∆a ≥ 0 sao cho:
|a − a
∗
| △
a
. (1.1)
Định nghĩa 1.2.5. Số ∆a nhỏ nhất thỏa mãn điều kiện (1.2.1) gọi là sai số tuyệt đối
của số gần đúng a.
Khi đó a

∗
= a ± ∆a.
Định nghĩa 1.2.6. Số δa =
∆a
|a|
được gọi là sai số tương đối của a.
Ví dụ 1.2.3. Giả sử a
∗
= π và a = 3, 14.
Do 3, 14 < π < 3, 15 = 3, 14 + 0, 01 nên ta có thể lấy ∆a = 0, 01.
Do 3, 14 < π < 3, 142 = 3, 14 + 0, 002 nên ta có thể lấy ∆a = 0, 002.
Ví dụ 1.2.4. Đo độ dài đoạn thẳng AB và CD ta thu được a = 10m ± 0, 01m; b =
1m ± 0, 01m.
Khi đó, ta có: δa =
0, 01
10
= 0, 1%; δb =
0, 01
1
= 1%.
Vậy phép đo đoạn thẳng AB chính xác hơn đoạn thẳng CD tuy chúng có cùng
sai số tuyệt đối ∆a = ∆b = 0, 01m.
• Sai số tuyệt đối cũng như sai số tương đối của một số gần đúng a của số đúng a
∗
là không duy nhất.
• Độ chính xác của một phép đo phản ánh qua sai số tương đối.
9
1.2.4 Xấp xỉ tốt nhất
Định nghĩa 1.2.7. Cho X là không gian định chuẩn với chuẩn ∥∥, M ⊂ X và p ∈ X.
Điểm y

0
∈ M được gọi là xấp xỉ tốt nhất tới p từ M nếu:
∥p − y
0
∥ ≤ ∥p − y∥, ∀y ∈ M
Xấp xỉ tốt nhất có thể tồn tại, có thể không và sự tồn tại nếu có cũng không
phải là duy nhất.
Ví dụ về sự không tồn tại xấp xỉ tốt nhất.
Cho X = R
2
, M = (x, 0) : x ̸= 0, p = (0, 1) thì trong trường hợp này không có
xấp xỉ tốt nhất từ M tới p.
Ví dụ về sự không duy nhất xấp xỉ tốt nhất:
Cho X = R
2
, M = (1, y) : y ∈ R ∪ (−1, y) : y ∈ R, p = (0; 0). Trong trường hợp
này tồn tại hai xấp xỉ tốt nhất z
1
= (−1, 0), z
2
= (1, 0).
Định lý 1.2.1. Nếu X là không gian định chuẩn với chuẩn ∥∥ và X
N
là không gian
con hữu hạn chiều của X thì với mỗi x ∈ X tồn tại xấp xỉ tốt nhất x
N
∈ X
N
, tức là:
∥x − x

N
∥ = min
y∈X
N
∥x − y∥
Chứng minh. Lấy z ∈ X
N
và đặt d = ∥x − z∥
K = z ∈ X
N
: ∥x − z∥ ≤ d
Từ ∥x∥ là hàm liên tục của x nên K là tập đóng và bị chặn. Mà K là không
gian hữu hạn chiều nên K compact.
Đặt g(z) = ∥x − z∥, z ∈ K. Khi đó, g là hàm liên tục của z.
Do K compact nên g đạt giá trị nhỏ nhất tại một số điểm x
N
∈ K.
Vậy ∥x − x
N
∥ = min
y∈K
∥x − y∥.
1.2.5 Tốc độ hội tụ của nghiệm xấp xỉ
Cho đoạn [a, b], chia đoạn [a, b] thành n phần bằng nhau bởi các điểm chia
x
i
, i = 0, n thỏa mãn:
x
0
= a < x

1
< x
2
< . . . < x
n
= b
10
Đặt h =
b − a
n
.
Giả sử x là nghiệm đúng và x
∗
là nghiệm xấp xỉ của phương trình đã cho (theo
phương pháp gần đúng nào đó). Nếu có:
∥x − x
∗
∥ = 0(h
k
)
thì x
∗
được gọi là hội tụ bậc k về nghiệm x.
1.2.6 Ma trận đường chéo trội
Định nghĩa 1.2.8. Cho ma trận vuông A = (a
ij
)
n
i,j=1
.

Ta nói ma trận A có tính chất đường chéo trội nếu nó thỏa mãn một trong hai
tính chất sau:
•

n
j=1,j̸=i
|a
ij
| < |a
ii
|, ∀i = 1, 2, . . . , n;
•

n
i=1,i̸=j
|a
ij
| < |a
jj
|, ∀j = 1, 2, . . . , n.
Định lý 1.2.2. Ma trận A có tính chất đường chéo trội thì A không suy biến.
Khi đó hệ phương trình












a
11
x
1
+ a
12
x
2
+ ···+ a
1n
x
n
= y
1
a
21
x
1
+ a
22
x
2
+ ···+ a
2n
x
n
= y

2
···
a
n1
x
1
+ a
n2
x
2
+ ···+ a
nn
x
n
= y
n
luôn có nghiệm.
Đặt ma trận A = (a
ij
)
n
i,j=1
, y = (y
1
, y
2
, ··· , y
n
)
T

, ta được phương trình Ax = y
luôn có nghiệm.
1.3 Không gian Hilbert
1.3.1 Mở đầu
Định nghĩa 1.3.1. Cho X là một không gian tuyến tính.
Ánh xạ ψ : X × X → R thỏa mãn các điều kiện:
1. ψ(x, x) ≥ 0, ∀x ∈ X;
11
2. ψ(x, x) = 0 ⇔ x = θ;
3. ψ(x, y) = ψ(y, x), ∀x, y ∈ X ;
4. ψ(αx
1
+ βx
2
, y) = αψ(x
1
, y) + βψ(x
2
, y), ∀x
1
, x
2
, y ∈ X và ∀α, β ∈ R.
được gọi là môt tích vô hướng trên X, còn ψ(x, y) được gọi là tích vô hướng của
hai phần tử x,y và thường được kí hiệu là (x, y).
Nhận xét. Nếu X là một không gian tuyến tính trên đó có xác định một tích vô
hướng (·),khi đó ánh xạ ∥· ∥ : X → R xác định bởi ∥x∥ =

(x, x) là một chuẩn trên
X và X cùng với chuẩn đó là một không gian tuyến tính định chuẩn.

Chuẩn xác định như trên được gọi là chuẩn cảm sinh bởi tích vô hướng.
Từ đó có ánh xạ d : X × X → R xác định bởi:
d(x, y) = ∥x − y∥ =

(x − y, x − y)
là một hàm khoảng cách trên X và (X, d) là một không gian metric. Khoảng
cách d vừa xác định được gọi là khoảng cách cảm sinh bởi tích vô hướng.
1.3.2 Định nghĩa không gian Hilbert
Định nghĩa 1.3.2. Cho không gian tuyến tính X cùng với tích vô hướng (·) . Nếu
cùng với khoảng cách d cảm sinh bởi tích vô hướng mà (X, d) trở thành một không gian
metric đủ thì lúc đó X cùng với tích vô hướng (·) được gọi là một không gian Hilbert.
1.3.3 Định nghĩa hệ trực chuẩn
Định nghĩa 1.3.3. Hệ vô hạn các phần tử {x
i
}
i∈I
thuộc không gian tuyến tính X được
gọi là độc lập tuyến tính nếu như mọi hệ con hữu hạn các phần tử của nó là độc lập
tuyến tính.
Định nghĩa 1.3.4. Cho X là một không gian Hilbert. Hệ các phần tử {e
i
}
i∈I
của X
được gọi là trực chuẩn nếu:
(e
i
, e
j
) = δ

ij
=

1 nếu i = j
0 nếu i ̸= j
Định lý 1.3.1. Giả sử {x
i
}
i∈I
là một hệ độc lập tuyến tính trong không gian Hilbert
X. Khi đó có thể xây dựng được một hệ {e
i
}
i∈I
trực chuẩn.
Thật vậy: Đặt e
1
=
x
1
∥x
1
∥
, y
2
= x
2
− (x
2
, e

1
)e
1
và e
2
=
y
2
∥y
2
∥
Giả sử đã có e
1
, e
2
, ··· , e
k−1
.
12
Ta đặt y
k
= x
k
−


k−1
i=1
(x
k

, e
i
)

e
i
và e
k
=
y
k
∥y
k
∥
. Khi đó hệ {e
i
}
i∈I
là hoàn
toàn xác định (vì nếu tồn tại một tỉ số k sao cho ∥y
k
∥ = 0 ⇔ y
k
= θ thì dẫn đến hệ
{x
1
, x
2
, ··· , x
k−1

} là phụ thuộc tuyến tính, trái với giả thiết). Dễ thấy: (e
1
, e
1
) = 1.
Xét (e
2
, e
1
) =

y
2
∥y
2
∥
, e
1

=
1
∥y
2
∥
(y
2
, e
1
) =
1

∥y
2
∥
(x
2
− (x
2
, e
1
).e
1
, e
1
) =
1
∥y
2
∥
[(x
2
, e
1
) −
(x
2
, e
1
)]
Vậy có (e
2

, e
1
) = 0. Dễ thấy (e
2
, e
2
) = 1.
Bằng quy nạp toán học ta thấy với k > h:
(e
k
, e
h
) =

y
k
∥y
k
∥
, e
h

=
1
∥y
k
∥
(y
k
, e

h
) =
1
∥y
k
∥
(x
k
−

k−1
i=1
(x
k
, e
i
)e
i
, e
h
)
=
1
∥y
k
∥
[(x
k
, e
h

) − (x
k
, e
h
)].
Từ đó (e
k
, e
h
) = 0, với k > h, ngoài ra (e
k
, e
k
) = 1 là rõ ràng. Như vậy hệ {e
i
}
i∈I
là hệ
trực chuẩn.
Quá trình xây dựng hệ {e
i
}
i∈I
từ hệ {x
i
}
i∈I
độc lập tuyến tính như trên được gọi là
quá trình trực chuẩn hóa Hilbert - Schmidt.
Ví dụ 1.3.1. Xét X = R

n
, với x = (x
1
, , x
n
) ∈ R
n
, y = (y
1
, , y
n
) ∈ R
n
, đặt
(x, y) =

n
i=1
x
i
y
i
. Có thể thấy R
n
cùng với tích vô hướng được xác định như trên là
một không gian Hilbert.
Ví dụ 1.3.2. Xét X = L
2
[a, b] là không gian các hàm bình phương khả tích trên đoạn
[a, b] bao gồm các hàm thực x(t) xác định, bình phương khả tích trên [a, b] sao cho:

b

a
p(t)x
2
(t)dt < +∞
trong đó p(t) là hàm trọng p(t) thường được chọn thỏa mãn các điều kiện: xác định và
khả tích trên [a, b], p(t) ≥ 0 trên [a, b] và p(t) = 0 chỉ trên một tập có độ đo 0). Ta
trang bị trên L
2
[a, b] một tích vô hướng bằng cách đặt: với x(t), y(t) ∈ L
2
[a, b] thì:
(x, y) =
b

a
p(t)x(t)y(t)dt
(có thể thấy tích phân này tồn tại hữu hạn ∀x(t), y(t) ∈ L
2
[a, b] do bất đẳng thức
Bunhiacopski dạng tích phân).
Không gian L
2
[a, b] với tích vô hướng vừa xác định là một không gian Hilbert.
Ví dụ 1.3.3. Xét trường hợp cụ thể của L
2
[a, b] ở trên với a = −1, b = 1, p(t) ≡ 1, và
xét hệ đa thức x
1

(t) ≡ 1, x
2
(t) ≡ t, , x
k
(t) = t
k−1
, k  2. Hãy trực giao hóa hệ {x
k
(t)}
nói trên tương tự quá trình trực chuẩn hóa Hilbert - Schmidt(chính xác đến một hằng
số nhân).
13
Nhận thấy: ∥x
1
∥ = 2, e
1
=
x
1
∥x
1
∥
, thay số ta có e
1
=
1
2
. Dễ thấy (x
2
, e

1
) =
1

−1
tdt = 0
nên y
2
= x
2
= t, ∀t ∈ [−1, 1].
Vậy ∥y
2
∥ =

1

−1
t.tdt

1
2
=

1
3
t
3





1
−1

1
2
=

2
3
. Vì e
2
=
y
2
∥y
2
∥
,thay số có e
2
=

3
2
.t, ∀t ∈
[−1, 1].
Ta có: (x
3
, e

1
) =
1

−1
t
2
.dt =
1
3
t
3




1
−1
=
2
3
(x
3
, e
2
) =
1

−1
t

2
.

3
2
.tdt =

3
2
.
1

−1
t
3
dt = 0
y
3
= x
3
− [(x
3
, e
1
)e
1
+ (x
3
, e
2

)e
2
], thay số được:
y
3
= t
2
−
1
3
⇒ ∥y
3
∥ =

1

−1

t
2
−
1
3

2
dt

1
2
,

rút gọn có: ∥y
3
∥ =
2
3
.

2
5
, từ đó có e
3
=
3
2
.

5
2
.

t
2
−
1
3

.v.v
Quá trình cứ tiếp tục như vậy, ta sẽ được một hệ trực chuẩn{e
i
}. Tuy nhiên do ta chỉ

quan tâm đến tính trực giao của hệ nên có thể nhân mỗi e
i
với một hằng số thích hợp
để được một vecto mới, vẫn kí hiệu là e
i
nhưng với dạng đơn giản hơn, như sau:
e
1
(t) ≡ 1, e
2
(t) ≡ t, e
3
=
1
2
(3.t
2
− 1)
e
4
(t) ≡
5t
3
− 3t
2
, e
5
(t) =
35t
4

− 30t
2
+ 3
8
e
6
(t) =
1
8
(63t
5
− 70t
3
+ 15t)v.v
Hệ đa thức {e
i
(t)}
i∈N
∗
trực giao thu được như trên gọi là hệ đa thức trực giao Legendre.
1.3.4 Một số kết quả cơ bản
Định nghĩa 1.3.5. Giả sử X là không gian Hilbert, còn {e
i
}
∝
i=1
là hệ trực chuẩn trong
X . Với mỗi x ∈ X, ta xét S
n
=


k
i=1
c
i
e
i
, với c
i
= (x, e
i
), thì S
n
được gọi là tổng
Fourier của x và các c
i
là các hệ số Fourier của x đối với hệ {e
i
}
∝
i=1
Nếu lim
n→∝
∥S
n
−x∥ =
0 thì người ta nói tổng Fourier S
n
hội tụ đến x và viết: x =


∝
i=1
(x, e
i
)e
i
.
Định lý 1.3.2. Cho không gian Hilbert X và {e
i
}
∝
i=1
là một hệ trực chuẩn của nó. Khi
đó các phát biểu sau tương đương:
1. ∀x ∈ X, x =

∝
i=1
(x, e
i
)e
i
.
14
2. ∀x ∈ X, ∥x∥
2
=

∝
i=1

(x, e
i
)
2
3. Nếu z là một phần tử trong X sao cho (z, e
i
) = 0, ∀i ∈ N
∗
thì z = θ.
4. Bao đóng của không gian con sinh bởi {e
i
}
∝
i=1
trùng với X.
Chương 2
MỘT SỐ VẤN ĐỀ VỀ SPLINE VÀ
B-SPLINE
2.1 Mở đầu về hàm spline và B-spline
2.1.1 Tổ hợp lồi và bao lồi
a) Tổ hợp lồi
Định nghĩa 2.1.1. Cho c
1
, c
2
là hai số, λ ∈ [0; 1]. Khi đó, số
c = (1 − λ).c
1
+ λ.c
2

(2.1)
gọi là tổ hợp lồi của hai số c
1
và c
2
. Trường hợp đặc biệt, có thể chọn trung bình cộng
của c
1
và c
2
, tức λ = 1 /2.
Trong R
2
, cho c
1
= (x
1
, y
1
), c
2
= (x
2
, y
2
) với x
i
, y
i
∈ R, i = 1, 2. Tổ hợp lồi của hai

điểm c
1
, c
2
là các điểm c(x; y) được xác định như sau:
(C) = {c(x, y) |x = (1 − λ)x
1
+ λx
2
, y = (1 − λ)y
1
+ λy
2
} (2.2)
0 ≤ λ ≤ 1
Để thuận tiện ta quy ước kí hiệu
c = (1 − λ)c
1
+ λc
2
(2.3)
thay cho giải thích (2.2)
b) Bao lồi của một tập hợp các điểm
Định nghĩa 2.1.2. • Bao lồi của hai điểm:
Trong không gian R
2
, cho c
1
= (x
1

, y
1
), c
2
= (x
2
, y
2
), x
i
, y
i
∈ R, i = 1, 2. Tập hợp
15
16
Hình 2.1: Tổ hợp lồi của c
1
và c
2
các điểm c thỏa mãn: c = (1 −λ).c
1
+ λ.c
2
, (với 0 ≤ λ ≤ 1) gọi là bao lồi của hai
điểm c
1
và c
2
.
• Bao lồi của n điểm:

Giả sử (c
i
)
n
i=1
là n điểm trong R
2
. Bao lồi của n điểm là tập hợp các điểm c
thỏa mãn: c = λ
1
.c
1
+ λ
2
.c
2
+ + λ
n
.c
n
, với n số λ
i
, thỏa mãn
n

i=1
λ
i
= 1 và
0 ≤ λ

i
≤ 1, i = 1, 2. , n.
Hình 2.2: Bao lồi của 3 điểm c
1
, c
2
, c
3
17
a,hai ñiểm
b,ba ñiểm
c,bốn ñiểm
d,năm ñiểm
Hình 2.3: Bao lồi (phần được bôi đen) của các điểm (các chấm đen)
2.1.2 Các khái niệm cơ bản
Cho hai điểm c
0
= (x
0
, y
0
) và c
1
= (x
1
, y
1
), x
i
, y

i
∈ R, i = 0, 1. Gọi AB là đoạn
thẳng đi qua c
0
và c
1
thì đoạn thẳng AB là bao lồi của hai điểm trên. Phương trình
của đoạn thẳng AB là:

x = x
0
+ (x
1
− x
0
).λ
y = y
0
+ (y
1
− y
0
)λ
⇔

x = (1 − λ)x
0
+ x
1
λ

y = (1 − λ)y
0
+ y
1
λ.
, λ ∈ [0; 1] (2.4)
Đặt λ =
t − t
0
t
1
− t
0
, t ∈ [t
0
; t
1
] ⇒ 1 − λ = 1 −
t − t
0
t
1
− t
0
=
t
1
− t
t
1

− t
0
Từ (2.4) ta suy ra





x =
t
1
− t
t
1
− t
0
x
0
+
t − t
0
t
1
− t
0
x
1
y =
t
1

− t
t
1
− t
0
y
0
+
t − t
0
t
1
− t
0
y
1
18
Kí hiệu:
q (t|c
0
, c
1
; t
0
, t
1
) =
t
1
− t

t
1
− t
0
c
0
+
t − t
0
t
1
− t
0
c
1
(2.5)
với t ∈ [t
0
, t
1
] , ∀t
0
, t
1
∈ R. Biểu thức trong 2.5 là một tổ hợp lồi của c
0
và c
1
. Nếu đặt
λ =

t − t
0
t
1
− t
0
thì 2.5 sẽ trở thành 2.3. Từ 2.5 có thể biểu thị phương trình đường thẳng
AB
y = f(x) =
x
1
− x
x
1
− x
0
.y
0
+
x − x
0
x
1
− x
0
.y
1
với (x; y) ∈ AB, A(x
0
; y

0
); B(x
1
; y
1
).
2.1.3 Nội suy các đường cong đa thức
a) Nội suy bậc hai của ba điểm
Giả sử c
0
, c
1
; c
2
∈ R là ba điểm đã cho, (t
j
)
2
j=0
là các số thực cho trước.
Đặt q
0;1
(t) = q(t|c
0
, c
1
, t
0
, t
1

) =
t
1
− t
t
1
− t
0
c
0
+
t − t
0
t
1
− t
0
c
1
và q
1,1
(t) = q (t|c
1
, c
2
; t
1
, t
2
) =

t
2
− t
t
2
− t
1
c
1
+
t − t
1
t
2
− t
1
c
2
thì
q
0,2
(t) = q (t|c
0
, c
1
, c
2
; t
0
, t

1
, t
2
) =
t
2
− t
t
2
− t
0
.q
0,1
(t) +
t − t
0
t
2
− t
0
.q
1,1
(t)
là đường cong bậc hai ẩn t.
q
0,2
(t) đi qua ba điểm c
0
, c
1

, c
2
đã cho vì
q
0,2
(t
0
) = q
0,1
(t
0
) = c
0
q
0,2
(t
1
) =
t
2
− t
1
t
2
− t
0
.q
0,1
(t
1

) +
t
1
− t
0
t
2
− t
0
.q
1,1
(t
1
) =
t
2
− t
1
t
2
− t
0
.c
1
+
t
1
− t
0
t

2
− t
0
.c
1
= c
1
q
0,2
(t
2
) = q
1,1
(t
2
) = c
2
Ta nói q
0,2
(t) là đường cong bậc hai 2 nội suy qua ba điểm c
0
, c
1
, c
2
b) Nội suy đa thức bậc cao
• Nội suy đa thức bậc ba.
Giả sử (c
i
)

3
i=0
là các điểm đã cho, gọi t = (t
i
)
3
i=0
, t
i
∈ R là các tham số cho trước.
Khi đó q
0,2
(t) là đường cong bậc hai nội suy cho các điểm c
0
, c
1
, c
2
và q
1,2
(t) là
đường cong bậc hai nội suy cho các điểm c
1
, c
2
, c
3
thì:
q
0,3

(t) =
t
3
− t
t
3
− t
0
.q
0,2
(t) +
t − t
0
t
3
− t
0
.q
1,2
(t)
là đường cong nội suy bậc ba cho bốn điểm c
0
, c
1
, c
2
, c
3
.
19

• Nội suy đa thức bậc d
Định nghĩa 2.1.3. Giả sử d ∈ N,(c
i
)
d
i=0
là các điểm đã cho, (t
i
)
d
i=0
, (t
i
∈ R, i =
0, d), là các tham số cho trước.
Từ định nghĩa nội suy bậc hai, bậc ba ta có thể xây dựng định nghĩa nội suy bậc
d bằng quy nạp. Cụ thể đặt:
q
0,d
(t) =
t
d
− t
t
d
− t
0
.q
0,d−1
(t) +

t − t
0
t
d
− t
0
.q
1,d−1
(t) (2.6)
thì q
0,d
(t) là đường cong bậc d. (Đồ thị của đa thức bậc d đi qua c
i
, i = 0, d)
Ta xây dựng thuật toán tìm q
o,d
(t) như sau:
Thuật toán 2.1.1 (Phương pháp Neville-Aitken).
Bước 1: Tính q
i,0
(t)
Bước 2: Sử dụng công thức quy nạp lần lượt tìm q
i,1
; q
i,2
q
i,d
như sau:
q
i,r

(t) =
t
i+r
− t
t
i+r
− t
i
q
i,r−1
(t) +
t − t
i
t
i+r
− t
i
q
i+1,r−1
(t)
Với i = 0, 1, d − r và r = 1, 2, d.
Các phép tính liên quan tới thuật toán tính q
0,d
(t) được tóm tắt như sau:
2.1.4 Xây dựng đường cong spline
a) Spline tuyến tính
Định nghĩa 2.1.4. 1. Cho hai điểm c
1
và c
2

và các tham số (t
i
)
3
i=2
với t
1
tùy ý,
t
2
< t
3
, đặt:
p(t|c
1
, c
2
; t
2
, t
3
) =
t
3
− t
t
3
− t
2
c

1
+
t − t
2
t
3
− t
2
c
2
, t ∈ [t
2
, t
3
] ,
thì p(t|c
1
, c
2
; t
2
, t
3
) là đoạn thẳng qua c
1
, c
2
.
2. Cho (c
i

)
n
i=1
là n điểm cho trước. Chọn (t
i
)
n+1
i=2
với t
i
< t
i+1
, i = 2, n
Xác định:
f (t) =











p (t|c
1
, c
2

; t
2
, t
3
) , t ∈ [t
2
, t
3
)
p (t|c
2
, c
3
; t
3
, t
4
) , t ∈ [t
3
, t
4
)
. . .
p (t|c
n−1
, c
n
; t
n
, t

n+1
) , t ∈ [t
n
, t
n+1
]
(2.7)
Thì f(t) là spline tuyến tính nội suy các (c
i
) đã cho. Các điểm (c
i
)
n
i=1
là các điểm
điều khiển, (t
i
)
n+1
i=2
gọi là các điểm nút.