Luận văn thạc sĩ sự ổn định của phương pháp sắp xếp spline đối với phương trình vi tích phân

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

Đ IN H T H Ị T H U

S ự Ổ N Đ ỊN H C Ủ A P H Ư Ơ N G P H Á P S A P X E P
SPL IN E ĐỐ I VỚI P H Ư Ơ N G T R ÌN H
V I-T ÍC H P H Â N

L U Ậ N VĂN TH ẠC s ĩ TO Á N HỌC
C h u yên ngành: T oán giải tích
M ã số : 60 46 01 02

Người hướng dẫn khoa học
T S . N g u y ễn V ăn Tuấn

H À N Ộ I, 2015


Lời cảm ơn
Tôi xin bày tỏ lòng biết ơn chân thành và sâu sắc đến TS Nguyễn
Văn Tuấn, ngưòi thầy đã định hướng chọn đề tài và nhiệt tình hướng
dẫn để tôi có thể hoàn thành luận văn này.
Tôi cũng 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 giảng dạy 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ôi trong suốt quá trình học tập tại trường.
Nhân dịp này tôi cũng xin gửi lòi cảm ơn đến gia đình, bạn bè đồng
nghiệp, BGH và tổ KHTN trường THCS Xuân Hòa thị xã Phúc Yên
tỉnh Vĩnh Phúc đã cổ vũ, động viên, tạo điều kiện để tôi hoàn thành
luận văn này.
Hà Nội, tháng 6 năm 2015

Tác giả

Đ in h T h ị T h u


Lời cam đoan
Tôi xin cam đoan, dưới sự chỉ bảo và hướng dẫn của TS Nguyễn Văn
Tuấn, luận văn chuyên ngành Toán giải tích với đề tằi:((S ự ổ n đ ị n h
c ủ a p h ư ơ n g p h á p s ắ p x ế p s p l i n e đ ố i v ớ i p h ư ơ n g t r ì n h v i t íc h
p h â n ” được hoàn thành bởi sự nhận thức và tìm hiểu của bản thân
tác giả.
Trong quá trình nghiên cứu và thực hiện luận văn, tác giả đã kế thừa
những kết quả 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 2015
Tác giả

Đ in h T h ị T h u


Muc luc

M ỏ dầu

5

1

K iến th ứ c cơ bản

8

1.1

Không gian v e c t o ................................................

8

1.2

Không gian định c h u ẫ n ......................................

12

1.2.1

Khái niệm không gian định chuẫn

. .

12

1.2.2

Sự hội tụ trong không gian định chuẫn

16

1.3

Không gian H i l b e r t .............................................

16

1.4

Không gian các hàm s p l i n e ...............................

19

1.4.1

Spline đa thức bậc ba vòi mốc cách đều

19

1.4.2

Spline đa thức tổng quát

24

1.5

2

Sai số, tốc đô hôi t u .................

28

1.5.1

Sai số

28

1.5.2

Xấp xỉ tốt nhất

1.5.3

Tốc đô hôi tu của nghiêm xấp x ỉ ...........................

30

1.5.4

Ma trận đường chéo trội

31

1.5.5

Các khái niệm cờ bản của lý thuyết ổn định . . .

32

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

29

Sự ổn đ ịn h củ a phương pháp sắp xếp sp lin e đối với
phư ơng trìn h vi tíc h ph ân

36

2.1

36

Dịnh nghĩa phường pháp sắp xếp spline

2


2.2

Sự ốn định của phương pháp sắp xếp spline vdi phương
trình vi phân
2.2.1

38

Sử dụng ma trận đưòng chéo trội nghiên cứu tính
ốn định của phương trình vi p h â n ........................

2.2.2

Sự ốn định của phương pháp sắp xếp spline cho

phương trình vi phân bậc hai

2.3

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

45

Sự ốn định của phương pháp sắp xếp spline vdi phương
trình vi tích phân
2.3.1

48

Phương pháp sắp xếp spline cho phương trình vi
tích phân Volterra bậc hai

2.3.2

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

48

Sự ốn định của phương pháp sắp xếp spline cho
phương trình vi tích phân Volterra bậc hai

3

38

54

ứ n g dụng

58

3.1

ữ n g dụng với phương trình vi phân

58

3.2

Ung dụng với phương trình vi tích phân

60

K ế t lu ận

67
68

T ài liệu t h a m k h ả o

3


B Ả N G KÍ HIỆU
N

N*
R

Tập số tự nhiên
Tập số tự nhiên khác không
Tập số thực
c Tập số phức
C[o Ị,] Tập tấ t cả các hàm số thực liên tục trên \a, b]
Ss(ĩĩ) Tập tấ t cả các hàm spline đa thức bậc 3
II • II Chuẩn

4


Mở đầu
1. Lí do chọn đề tài
Trong khoa học tự nhiên, kĩ thuật, trong kinh tế, cũng như trong các
lĩnh vực khác của cuộc sống ta gặp rất nhiều bài toán đưa tới việc nghiên
cứu các phương trình vi phân, phương trình vi tích phân.... Giải đúng
phương trình vi tích phân rất khó vì vậy ngưòi ta thường áp dụng các
phương pháp xấp xỉ để giải. Có rất nhiều phương pháp giải gần đúng
khác nhau, phương pháp sắp xếp spline là một phương pháp thường
được lựa chọn.
Ưu điểm của phương pháp sắp xếp spline là sử dụng các hàm đa thức
tính toán để giải. Các hàm đa thức dễ dàng lập trình đưa lên máy tính,
tính toán thuận lợi, hiệu quả. Trong một số trường hợp phương pháp
sắp xếp spline thường đạt tốc độ hội tụ nhanh, độ chính xác của nghiệm
gần đúng tốt hơn các phương pháp khác. Có thể khái quát cho nghiệm
xấp xỉ bằng spline bậc cao hoặc các hàm B-spline.
Sự ổn định của nghiệm xấp xỉ luôn được các nhà Toán học trong

và ngoài nước quan tâm nghiên cứu. Với mong muốn tìm hiểu về sự ổn
định của nghiệm xấp xỉ và phương pháp sắp xếp spline nhằm nâng cao
kiến thức đã học trong chương trình đại học và cao học nên em chọn đề
tài này để nghiên cứu và làm luận văn tốt nghiệp cho mình.

5


2. M ục đích nghiên cứu
Tìm hiểu khái niệm và các định lý cơ bản của phương pháp sắp xếp
spline .
Nghiên cứu sự ổn định của phương pháp sắp xếp spline với phương
trình vi phân và phương trình vi tích phân.

3. N hiệm vụ nghiên cứu
Nghiên cứu ổn định của phương pháp sắp xếp spline với phương trình
vi phân, phương trình vi tích phân.
Nghiên cứu về lập trình Maple để ứng dụng.

4. Đ ối tượng và phạm vi nghiên cứu
Đối tượng nghiên cứu: “Sự ổn định của phương pháp sắp xếp spline ”.
Phạm vi nghiên cứu: các khái niệm, định lý các kết quả cơ bản của
phương pháp sắp xếp spline. Các phương trình vi phân, vi tích phân.
Lập trình Maple với phương pháp sắp xếp spline.

5. Phương pháp nghiên cứu
Sử dụng các phương pháp phân tích, tổng hợp, phương pháp lấy ý
kiến chuyên gia.

6. Đ óng góp mới

Sẽ nghiên cứu sự ổn định của một lớp phương trình vi tích phân
bằng phương pháp sắp xếp spline, có thể chứng minh được sự ổn định

6


của một lớp phương trình vi tích phân bằng phương pháp sắp xếp spline.

7


Chương 1

Kiến thức cơ bản
Chương này trình bày một số không gian thường dùng như: Không
gian vectơ, không gian định chuẩn, không gian Hilbert, không gian các
hàm spline, sai số, tốc độ hội tụ, sự ổn định của nghiệm để phục vụ
chứng minh ở chương sau.

1.1

K hông gian vectơ

Đ ị n h n g h ĩa 1.1.1. Cho V là một tập khác rỗng mà các phần tử được
kí hiệu : ã, /3, 7 ,. .. và I là một trường mà các phần tử được kí hiệu:
x ,y ,z,. . .
Giả sử trên V được trang bị hai phép toán, gồm:
1. Phép toán cộng, kí hiệu + :
Y X Y — >V.
(ã, ệ ) I— >ã + ệ .

2. Phép toán nhân, kí hiệu là • :
K X V — >V.
(X, ã) I— > X ■ã.

8


thỏa mãn các tiên đề sau:
• ã + ặ = ặ + ã ì Va, ¡3 G V;

• ( â + ¡3) + 7 = ã + 0 + 7 ), Va, ệ, 7 G V;

• tồn tại 9 G V sao cho ớ + a = a + ớ = a , Va G V ;

• Với mỗi a tồn tại a ' G V sao cho a' + ã = ã + a' = 0]

• (2; + y ) ã = x ã + y ã , Va G V và 2;, y G K;

• a:(a + ¡3) = a:a + xj3, Va, ¡3 G V và 2; G K;

• 2;(ya) = (a:y)a,Va G V và x , y G K;

• 1 • ã = ã, Va G V và 1 là phần tử đơn vị của trường K.
Khi đó V cùng với hai phép toán trên gọi là không gian vectơ
trên trường K, hay K-không gian vectơ,hay không gian tuyến tính.

Khi K = R thì V được gọi là không gian vectơ thực.
Khi K =

c thì

V được gọi là không gian vectơ phức.

Ví dụ: Dễ dàng kiểm tra Cịa, b] là một không gian vectơ.

9


Đ ịn h n g h ĩa 1.1.2. Trong không gian vectơ V
Hệ vectơ ( « 1 , .. . , ã n) được gọi là độc lập tuyến tính nếu hệ thức:
Aiõ:i + • • • + Ằnã n = 0
Chỉ xảy ra khi Ai = • • • = A„ = 0.
Hệ vectơ («!, .. . , ã n) được gọi là phụ thuộc tuyến tính nếu nó không
độc lập tuyến tính.
Ví dụ: Trong không gian vectơ thực R 2 cho hệ ba vectơ:
ã i = (2;0);ữ2 = (0;4);ữ3 = (4; 4)
Hệ

độc lập tuyến tính vì :
AiCüi + A202 = 0 => (2Ai; 0) + (0; 4 A2) = (0; 0)
^ ( 2 A 1;4A2) = ( 0 ; 0 ) ^ A 1 = A2 = 0

Hệ (0 1 , 0:2 , 0:3) phụ thuộc tuyến tính vì 2c?i + c?2 — Ổ3 = 0.
Đ ị n h n g h ĩ a 1.1.3. Giả sử V là một không gian vectơ
Một hệ vectơ của V được gọi là một hệ sinh của V nếu mọi vectơ của
V đều biểu thị tuyến tính qua hệ đó.
Khi V có một hệ sinh gồm hữu hạn phần tử thì V được gọi là không
gian vectơ hữu hạn sinh.
Một hệ vectơ của V được gọi là một cơ sở của V nếu nó là hệ sinh
độc lập tuyến tính.

Đ ị n h n g h ĩ a 1 .1 .4 . Cho V là không gian vectơ có cơ sở gồm hữu hạn
phần tử thì số phần tử trong cơ sở đó được gọi là số chiều của không
gian vectơ.

10


Khi Y là một K-không gian vectơ có số chiều n ta kí hiệu:
d im V = n( hay dim K V = n).
Nếu V = {0} ta quy ước dim V = 0.
Nếu V không có cơ sở nào gồm hữu hạn phần tử thì nó được gọi là
không gian vectơ vô hạn chiều.
Ví dụ: Trong K —không gian vectơ I n xét hệ vectơ:

(e) = {ẽỉ = (1, 0, . . . ,0),ẽỉ = (0, 1,... , 0 ) , . . . ,ẻ; = ( 0 , 0 , . . . , 1)}
n
thì với mọi ã = ( x i , X 2 , ■■■, x n) G I n ta đều có ã =

Xịẽị

nên (e) là

i=1

n
một hệ sinh của I n. Mặt khác, nếu có
i=1

AịéTị = 0 thì (Ai, ■■■, A„)

( 0 , . . . , 0). Suy ra Ai = • • • = A„ = 0.
Vậy hệ (e) còn là một hệ độc lập tuyến tính. Do đó hệ (e) là một cơ
sở của K". Cơ sở này được gọi là cơ sở chính tắc ( hay cơ sở tự nhiên )
của K n. Từ đó suy ra dirriKn = n.

...

Đ ịn h n gh ĩa 1 1 5 Tập con

w

Ỷ 0

một i f -không gian vectơ E

được gọi là không gian vectơ con của E nếu nó ổn định với hai phép
toán của E , nghĩa là thỏa mãn các điều kiện sau:
1. V a J e W , a + ß e W ,
2. Va G

w và Va: G K

thì x ã G

w.

( xem [6])

11

1.2

K hông gian định chuẩn

1.2.1

K h á i niệm k h ô n g gian đ ịn h chuẩn

Đ ị n h n g h ĩ a 1.2.1. Ta gọi không gian định chuẩn (hay không gian tuyến
tính định chuẩn) là không gian tuyến tính X trên trường p ( P = K
hoặc p = C) cùng với một ánh xạ từ X vào tập số thực K, kí hiệu là
||.|| và đọc là chuẩn,thỏa mãn các tiên đề sau đây:
1. (Va; G X ) ||a;|| > 0, ||a;|| = 0

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

2. (Va; G X ) (Va G p ) ||aa;|| = |a | ||a;|| ;
3. (Va;, y e X ) \\x + y\\ < ||a;|| + ||y|| .
Số ||a;|| gọi là chuẩn của vectơ X. Ta cũng kí hiệu không gian định
chuẩn là X . Các tiên đề trên gọi là hệ tiên đề chuẩn.
Ví dụ: Không gian K2 là không gian định chuẩn với chuẩn thường
chọn là chuẩn:

Ngoài ra còn có những chuẩn khác chẳng hạn như:
X 1

Xị \ + \ x 2

hay

m a x { | Xị I, I x 2 |}.
trong đó X = (a;i,a; 2) G M2. Ví dụ: Không gian C[a,b] = { / : [a,b] —>
R I / liên tục trên [ữ, ồ]} là không gian định chuẩn


Đ ị n h n g h ĩ a 1.2.2. Cho không gian tuyến tính X với hai chuẩn lllli và
II ||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à ra > 0 sao cho:
ra II Xi ||<|| x 2 ||< M II Xi II, Va; G X .
Đ ị n h lý 1.2.1. Cho không gian định chuẩn X . Đối với hai vectơ bất kì
X, y G X ta đặt
d { x , y ) = \\x - y\\ .

(1.1)

Khi đó d là một metric trên X .
Chứng minh của định lý trên dễ dàng suy ra từ hệ tiên đề chuẩn và
hệ tiên đề tuyến tính.
Mọi không gian định chuẩn đều có thể trở thành không gian metric
với metric (1.1). Do đó mọi khái niệm,mệnh đề đã đúng trong không
gian metric đều đúng trong không gian định chuẩn.
Đ ị n h n g h ĩ a 1.2.3. Dãy (x n) trong không gian định chuẩn X được gọi
là hội tụ đến Xq & X nếu
lim \\xn — XqII = 0.

n —>00

Khi đó ta kí hiệu
lim x n = Xq hoặc x n —> Xq khi n —>• oo.

n —>00

Đ ị n h n g h ĩ a 1.2.4. Dãy điểm (x n) trong không gian định chuẩn X gọi
là dãy cơ bản, nếu
lim

m ,n —>oo

\\xn — x m II = 0 .

Đ ị n h n g h ĩ a 1.2.5. 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ụ.
13


Đ ịn h n gh ĩa 1.2.6. Giả sử không gian định chuẩn X là một không gian
metric đầy đủ (với khoảng cách d ( x , y ) = 11rc — y II). Khi đó X được gọi
là một không gian định chuẩn đầy đủ, hay còn gọi là không gian Banach.
Đ ịn h n gh ĩa 1.2 .7 . Cho hai không gian tuyến tính X và Y trên trường
p . Ánh xạ A từ không gian định chuẩn X vào không gian định chuẩn Y
được gọi là ánh xạ tuyến tính hay toán tử tuyến tính nếu Ả thỏa mãn:
1. A(X + y) = A x + Ay\
2. A (a x ) = a A (a;).
-

A được gọi là toán tử cộng tính nếu A chỉ thỏa mãn

1.

- A được gọi là toán tử thuần nhất nếu A chỉ thỏa mãn 2.
- Khi Y = p thì toán tử tuyến tính A được gọi là phiếm hàm tuyến
tính.
Đ ịn h n gh ĩa 1.2 .8 . Cho hai không gian định chuẩn X và Y . Toán tử
tuyến tính A từ không gian X vào không gian Y gọi là bị chặn nếu tồn
tại hằng số c > 0 sao cho:
II A x ||< c II X II, với mọi

X

G X.

A được gọi là ánh xạ giới nội.
Đ ịn h lý 1.2 .2 . Ánh xạ A : X —> Y tuyến tính, liên tục khi và chỉ khi
A giới nội.
Chứng minh. Giả sử A giới nội. Lấy {x n} c X , x n —> X tương đương với
x n — X —> 0.
Ta có ||A(a;n) — A (s)|| = ||A(a;n — ar)II < k \\xn — s|| —> 0.
Suy ra d (A ( x n), A ( x ) ) = IIA ( x n) — ^4(rr) 11 —> 0 suy ra A ( x n) —> A {x)

14


do đó A liên tục.
Ngược lại, giả sử Ả liên tục nhưng Ả không giới nội.
Tức Vra > 0, 3 x n G X : ||A(æm)|| > ra ||æm||.
II ■ Ta

Ta đăt y m =
Suy ra {Vm}

đươc ||ym|| =

= — —> 0, ra —> 00.

0. Ta có \\A{ym)\\ =

\\A (ym)\\ -/* 0 o- A ( y m)

> 1. Suy ra

0 = A(0) (mâu thuẫn)

Vậy A giới nội.
Ta có điều phải chứng minh.
Đ ị n h n g h ĩ a 1.2.9. Cho hai không gian định chuẩn X và Y. Kí hiêu
L(J*r, Y ) là tập tấ t cả các toán tử tuyến tính bị chặn từ không gian X
vào không gian Y . Ta đưa vào L (X , Y ) hai phép toán:
1. Tổng của hai toán tử A, B G L (X , Y ) là toán tử,kí hiệu A + B , xác
định bởi biểu thức.
(A + B ) ( x ) = A x + B x , với mọi

2. Tích

vô hướng của a G p ( p =

R hoặc

P =

X

& X;

c)

với toán tử

AG

L (X , V) là toán tử, kí hiệu a A , được xác định bởi biểu thức
(aA)(x) = a (A x)
Dễ kiểm tra được rằng A + B G L (X , V), a.A

g

L (X , V) và hai

phép toán trên thỏa mãn tiên đề tuyến tính. Khi đó, tập L(Jf, V)
trở thành một không gian tuyến tính trên trường p .
Đ ị n h lý 1.2.3. Nếu Y là một không gian Banach thì L (X , Y ) là không
gian Banach.
( xem [5])
15


1 .2 .2

Sự h ộ i t ụ tr o n g k h ô n g g ian đ ịn h chu ẩn

Giả sử X là không gian định chuẩn và {æ„}“=1 c X , Xo £ X .
1. x n —> Xq ( dãy x n hội tu tới £o) có nghĩa là II x n — Xo II—> 0.
2. Nếu x n —> Xo thì II x n II—>\\

II, tức là chuẩn II x n II là một hàm

Xo

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ì 3 M £ R, M >
0, Vn, Il x n ||< M .
4. Nếu x n ->• Xo, y n -> Vo thì x n + y n ->• x 0 + y0.
5. Nếu x n —> X o ,a n —> CKo thì x na n —> X qOIo , V {a„}“=1 c l , ữ f l € R.
6. Một dãy cơ bản trong không gian định chuẩn X là một dãy { x n} c
X sao cho:
lim \\xn — x m II = 0.
m,n—>oo
Nếu trong không gian định chuẩn mọi dãy cơ bản đều hội tụ,tức là:
11x n — x m \\ —?• 0 kéo theo sự tồn tại

Xo

£ X sao cho x n —> Xo, thì không

gian đó được gọi là không gian đủ thường gọi là không gian Banach.

1.3

K hông gian H ilbert

Đ ịn h n gh ĩa 1.3.1. (Tích vô hướng) Cho không gian tuyến tính X trên
trường p (p là trường số thực R hoặc trường số phức c ). Ta gọi là tích
vô hướng trên không gian X mọi ánh xạ từ tích Descartes X
trường p ,k í hiệu (., .), thỏa mãn tiên đề:
1. (Væ, y £ X) ( y , x ) = {x,y)\
2. (Væ, y,

Z

£ X ) (x + y,

z)

= (x,

z)

16

+ (y,

z)

;

X

X vào

3. (Va;, y G X) (Va G p ) (a x , y) = a { x , y ) ]
4. (Va: G X ) (a;,a;) > 0, nếu X Ỷ 0 (ớ là kí hiệu phần tử không),
(a;, a;) = 0, nếu X = 6.
Các phần tử x , y , z , . . . gọi là các nhân tử của tích vô hướng, số (a;, y)
gọi là tích vô hướng của hai nhân tử X và y ,các tiên đề trên gọi là tiên
đề tích vô hướng.
Đ ị n h lý 1.3.1. Đối với mỗi X G X . Ta đặt
||z|| =

(x, x).

( 1 .2 )

Khi đó V x , y G X ta có bất đẳng thức Schwarz

\{x,y)\ < ||x|| \\y\\ .

(1.3)

Công thức (1^2) xác định một chuẩn trên không gian X .
Đ ị n h n g h ĩa 1.3.2. Không gian tuyến tính trên trưòng p cùng với một
tích vô hướng gọi là không gian tiền Hilbert.
Đ ị n h n g h ĩ a 1.3.3. Ta gọi một tập H Ỷ 0 sồm những phần tử X, y, z , . . .
nào đấy là không gian Hilbert,nếu tập H thỏa mãn các điều kiện:
1. H là không gian tuyến tính trên trường P\
2. H được trang bị một tích vô hướng (., .);
3. H là không gian Banach với chuẩn ||a;|| = - \ / ( x , x ) , x G H. Ta gọi
mọi không gian tuyến tính con đóng của không gian Hilbert H là
không gian Hilbert con của không gian H .
Đ ị n h n g h ĩa 1.3.4. (Trực giao) Cho không gian Hilbert H. Hai phần

tử x , y G H gọi là trực giao,ký hiệu X-Ly,nếu (x , y ) = 0.

17


Đ ị n h n g h ĩ a 1.3.5. Cho không gian Hilbert H và tập con Ả c H, Ả Ỷ
0. Phần tử

X

G H gọi là trực giao với tập A, nếu X-Ly (Vy E A ) và kí

hiệu X-LA.
Đ ị n h lý 1.3.2 (Định lý hình chiếu lên không gian con). Cho không gian
Hilbert H

và

Hq

là

không gian con của H . Khi

đó

phần tứ bất kì

X

GH

biểu diễn một cách duy nhất dưới dạng
X

= y + z,y

e

H 0, z

e

H 0.

(1.4)

Phần tử y trong biểu diễn (1.4) gọi là hình chiếu của phần tử X lên
không gian con H q.
Đ ị n h n g h ĩ a 1.3.6. (Hệ trực chuẩn) Cho không gian Hilbert H. Một
tập (còn gọi là hệ thống) gồm hữu hạn hay đếm được các phần tử
ien)n> 1 c H gọi là một hệ trực chuẩn, nếu
(eí ì ej) — õii
ỏij là kí hiệu Kroneckes,ổịj = 0 với

% Ỷ

(1-5)
=

1 v<3i í =

j,{hj

=

1, 2, . . . ) .
N h ậ n x é t: Không gian định chuẩn và không gian Hilbert có hai cấu
trúc tôpô và đại số.
Về cấu trúc tôpô: Họ lân cận của 0, IÁ = {Ua} aeI, Ua là lân cận của
0.
X

€

X, {x +

Ua} aeI là họ lân cận của

X.

Đ ị n h n g h ĩa 1.3.7. Mọi họ các lân cận của điểm 0 (Kí hiệu: ỈA ) được
gọi là họ cơ sở của lân cận nếu:
1. u bất kỳ là lân cận của điểm 0 thì tồn tại ƯQ c u sao cho ƯQ c U]

2.

Với U \, U2 E ỈA thì U\

n U2

E ỈA]
18


3. Với ưị G u , i = 1, • • • ,00 thì

u Uị G U;

i=1
4. Với

w &u,

tồn tại ƯQ G u sao cho ƯQ + ƯQ c

w.

1.4 K hông gian các hàm spline
1.4.1

S p lin e đa t h ứ c bậc b a vổi m ố c cách đ ề u

Xét phân hoạch 7ĩ trên đoạn [a, b] với các mốc nội suy
a = t0 < ti < t 2 < . . . < t n = b.
Kí hiệu hị = tị — t ị - 1 , nếu hi = h = const thì các mốc nội suy
to, t \ , ¿2 J ■■■Jt n gọi là các mốc nội suy cách đều.
Đ ị n h n g h ĩa 1.4.1. Một spline đa thức bậc ba trên đoạn [ữ, b] với phân
hoạch 7Ĩ là hàm số y = s ( t ) thỏa mãn hai điều kiện sau:
1. sịt) G c 2[a, 6];

2. Hạn chế của s(t) trên mỗi khoảng A ị = \tị]tị+1 ] là đa thức s(t) | a ,
với deg(s(t) I A) < 3, Vị = 0,1, 2, • • • ,n.
Không gian gồm tấ t cả các hàm số s(t) thỏa mãn hai điều kiện trên
kí hiệu là ^ (

tĩ).

Từ định nghĩa ta có không gian S 3(7ĩ) chứa tấ t cả các đa thức có bậc
nhỏ hơn hoặc bằng 3. Dễ dàng kiểm tra các tiên đề của không gian vec
tơ suy ra ^ ( tĩ) là không gian tuyến tính.
M ệ n h đ ề 1.4.1. Không gian

là không gian tuyến tính và không

gian đó chứa tất cả các đa thức có bậc nhỏ hơn bằng 3.
B à i t o á n 1.4.1. Tồn tại duy nhất hàm số s ( t ) G

19

thỏa mãn hệ


điều kiện:
V ( í o ) = /'(to),
< s{ti) =

0 < i < n,

(1.6)

ks'(ín) =
Khi đó s ( t ) được gọi là đa thức nội suy spline bậc ba của hàm số f ( t ) .
Xây dựng sự tồn tại của hàm s ( t ) với các mốc nội suy cách đều
i(b — a)
tị = to-ị------------ ,trong đó chúng ta bô sung thêm 4 môc nội suy t n- 2 <
n
t n- 1 < ÍQVầ t n+2 > t n+1 > t n đồng thòi xét hàm số B ị ( t ) được xác định
bởi công thức:
' { t - t ị - 2)3, í e [ti -2,ti- 1 ],
h 3 + 3h 2{t - t i- i ) + 3h{t - t ị - i ) 2 - 3 (í - í i - i ) 3,

Bị (í) = ^ < h3 + 3h2(tị+i - í) + 3/i(íi+i - í ) 2 - 3 (ti+i - í ) 3,
t G [ti,ti+1],
(tị+ 2 0 ) í € [íị+1; íị+ 2])
k0, í ^ [íị_2,íị + 2]'
(1.7)
Bằng cách thay vào (1.7) các hàm số B ị ( t ) liên tục khả vi hai lần
trên K, mà
' 4 ,3 = i,
Bi{tj) = < 1, j = ỉ — 1 hoặc j = ỉ + 1,
0 , j = ỉ — 2 hoặc j = ỉ + 2,

( 1 ,8 )

Đồng thời B ị ( t ) = 0 với t > tị+2 và t < t ị - 2 Các hàm ổ-spline khác không nhỏ gọn nhất với các mốc nội suy
t —2 < t - 1 < to < . . . < t n < t n+1 < t n + 2 đó là, bất kì spline đa thức bậc
3 s (í) đồng nhất triệt tiêu bằng 0 ngoài khoảng ( t j - 2 , t j +2 )Hơn nữa mỗi Bị(t) là bậc 3 trên [tj,tj+1 ] nên Bị(t) G So(ĩĩ).
Tính B j( t), B'-(t), B"{t) chúng ta có bảng sau:
20

t

t ị -2

tj _ 1

B iW

0

1

4

B '(i)

0

h

3

0

B'Ị{t)

0

h2

6

+1

12

1

3

h2

Bảng 1.1: Giá trị

tj +2

0

h

0

h2

0

6

B'ị (í), B'ị (í)

Giả sử B = -B- 1 ,£?0) . . . , B n+1 và £?3(7ĩ) = spanB.
Dễ thấy hệ B là độc lập tuyến tính và B ị (ĩ ĩ ) là không gian n + 3
chiều.
Đ ịn h lý 1.4.1. Có duy nhất hàm s ( t ) G £?3(tt) thỏa mãn bài toán (1.6) ■
Chứng minh. Giả sử s ( t ) G

thì

sịt) — X -\B _ i( t ) + XqB q{ì ) + • • • “h ^n + l -®71+ 1(í).

(1.9)

s{t) thỏa mãn bài toán (1.6ị) nên chúng ta có:

s'{t0) = X - ị B - i ( t ) + XqB qÍẠ,) + . . ■“b ^n + l-^n + l (í) = /'(to),
<

s(t») = I - l ổ - l ( í i ) + a:05o(íi) + . . . + In + iổ n + iíí,) = /( ti ) ,
0 < ĩ < n,
® (^n)

2C—\B —\{tn) ~\~XoBo[tn) “b ■■■“b XnjỊ-i BnJr\(
f (tn) ■
( 1.10 )

Đây là hệ phương trình tuyến tính gồm n + 3 phương trình dạng:

21

Ax = b, với
" X-1 '

-f'(to)-

x0

f(to )

X\

f(tl)

X=

,b =

Xn

f{tn )

.Xn +1.
và ma trận hệ số:
B'_ l{to)

B'o (to)

B [ ( k ) • ■ B'n+1 (to)

B _ l{to)

Bo (to)

Bi(to)

■ ■ B n-ị- (to)

l(íl)

Bo (ti)

B i( ti)

. ■ B n-ị- (ti)

B_ l{tn)

Bo (tn)

■ B n-ị- (tn)

B'_ 1{tn) B'o (tn)

3 0 3
h
h

■ B'n+1 (tn)_

0

0

0

1

4 1

0

0

0

0

1 4

1

0

0

0

0 1

4

1

0

0

0 0

0

1

3
hJ
A là ma trận có đường chéo trội nên không suy biến từ đó hệ (1.10 |có
0

0 0

0

nghiệm duy nhất.
22

0

Đ ị n h lý 1.4.2. Với các không gian

và B^(Tĩ)nêu trên chúng ta có:

5 3(tt) = s 3(7ĩ).

(1.11)

Chứng minh: Từ định nghĩa 5 3(7ĩ) ta có 5 3(7ĩ) c ¿^( tĩ).
Chúng ta đi chứng minh ¿^( tĩ) D 5 3(7ĩ ).
Lấy / ( í ) G Ssịĩĩ) khi đó f ' { t 0), f ' { t n) và f { t ị ) , 0 < i < n đều xác
định.
Giả sử s ( t ) G B ị (tĩ) là hàm spline duy nhất được xác định trong
định lý (|l.4.1 ),đặt / ( í ) — s ( t ) = g ( t) thì g(tị) = 0, 0 < i < n.
Vì f { t ) , g{t) G c 2[a, b] nên g(t) G c 2[a, h] theo định lý Rolle thì g'(t)
có n nghiệm Ui thỏa mãn tị < Ui < tị + 1 đồng thòi t o , t n là hai nghiệm
của g'(t). Như vậy g'(t) có ít nhất n + 2 nghiệm do đó g"{t) có ít nhất
n + 1 nghiệm Zị và Xi < Zị < ĩ/i, 0 < i < n .
Nhưng g"{t) là đa thức có bậc cao nhất bằng 1 trên t i , t i+ 1 với các
điểm lưới của phân hoạch 7Ĩ , vì g"(t) nhận Zị,i = 0,1, 2, . . . , n là các
nghiệm nên chúng ta suy ra g"(t) = 0 trên tị, tị + 1 .
Do đó g"(t) = 0 trên tị, tị + 1 từ đó chúng ta được: g(t) = a t + /3 . Mà
g{to) = g{tn) = 0 suy ra a = /3 = 0 hay g(t) = 0 trên t 0, t n .
Do đó
s{t) = f { t ) G S 3(7t).
Kết luận
ổ 3(7ĩ )

= 5 3(7r).

H ệ q u ả 1.4.1. ^ ( tĩ) là không gian tuyến tính n + 3 chiều với hệ cơ sở
B = j - 6 - i , B 0, . . . , B n+Ì I .

23