Phương pháp ổn định lyapunov nghiên cứu sự ổn định toàn cục của một số mô hình dịch tễ học

ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC SƯ PHẠM NGUYỄN HỮU THẮNG PHƯƠNG PHÁP ỔN ĐỊNH LYAPUNOV NGHIÊN CỨU SỰ ỔN ĐỊNH TOÀN CỤC CỦA MỘT SỐ MÔ HÌNH DỊCH TỄ HỌC LUẬN VĂN THẠC SĨ TOÁN HỌC THÁI

ĐẠI HỌC THÁI NGUYÊN

TRƯỜNG ĐẠI HỌC SƯ PHẠM

NGUYỄN HỮU THẮNG

PHƯƠNG PHÁP ỔN ĐỊNH LYAPUNOV NGHIÊN CỨU SỰ ỔN ĐỊNH TOÀN CỤC CỦA

MỘT SỐ MÔ HÌNH DỊCH TỄ HỌC

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

THÁI NGUYÊN-2020

ĐẠI HỌC THÁI NGUYÊN

TRƯỜNG ĐẠI HỌC SƯ PHẠM

NGUYỄN HỮU THẮNG

PHƯƠNG PHÁP ỔN ĐỊNH LYAPUNOV NGHIÊN CỨU SỰ ỔN ĐỊNH TOÀN CỤC CỦA

LỜI CAM ĐOAN

Tôi xin cam đoan rằng nội dung trình bày trong luận văn này là trung thực

và không trùng lặp với đề tài khác Tôi cũng xin cam đoan rằng mọi sự giúp đỡ cho việc thực hiện luận văn này đã được cảm ơn và các thông tin trích dẫn trong luận văn đã được chỉ rõ nguồn gốc

Thái Nguyên, tháng 11 năm 2020

Người viết luận văn

Nguyễn Hữu Thắng

Mặc dù đã rất cố gắng để luận văn đượ c hoàn thiện một cách tốt nhất nhưng

do điều kiện thời gian và năng lực của bản thân còn hạn chế, luận văn này không tránh khỏi những thiếu sót Tôi rất mong nhận được sự đóng góp ý kiến của thầy cô

và các bạn để luận văn này được hoàn thiện hơn

MỤC LỤC

Lời cảm ơn i

Lời cam đoan ii

Mục lục iii

Danh sách hình vẽ iv

Mở đầu 1

1 Kiến thức chu ẩn bị 4

1.1 Bài toán giá trị ban đầu 4

1.2 Định lý ổn định Lyapunov 7

1.3 Phương pháp Ru nge-Kutta bốn nấc kinh điển 9

1.4 Một số mô hình dịch tễ học cổ điển 10

2 Tính chất ổn định toàn cục của các mô hình dịch tễ SIR, SIRS, SIS và S I với các biến điều khiển phản hồi 13

2.1 Các mô hình SIR và SIRS 13

2.1.1 Phân tích ổn định 14

2.1.2 Các mô phỏng số 17

2.2 Mô hình SIS cổ điển 19

2.2.1 Phân tích ổn định 19

2.2.2 Các mô phỏng số 20

2.3 Mô hình SIS với tỷ lệ mắc bệnh chuẩn 22

2.3.1 Mô hình toán học 22

2.3.2 Phân tích ổn định 24

2.3.3 Các mô phỏng số 25

2.4 Các mô hình SI với các biến điều khiển phản hồi 27

2.4.1 Mô hình toán học 27

2.4.2 Phân tích ổn định 28

2.4.3 Các mô phỏng số 31

Kết luận chung 33

Tài liệu tham khảo 34

DANH SÁCH HÌNH VẼ

1.1 Sơ đồ lan truyền của mô hình SIS đơn giản 11

2.1 Sơ đồ lan truyền của mô hình SIS đơn giản 14

2.2 Nghiệm của mô hình trong Ví dụ 2.1 18

2.3 Nghiệm của mô hình trong Ví dụ 2.2 18

2.4 Sơ đồ lan truyền của mô hình SIS (2.14) 19

2.5 Nghiệm của mô hình trong Ví dụ 2.3 21

2.6 Nghiệm của mô hình trong Ví dụ 2.4 22

2.7 Nghiệm của mô hình (2.19) trong Ví dụ 2.5 26

2.8 Nghiệm của mô hình (2.19) trong Ví dụ 2.6 26

2.9 Nghiệm của mô hình (2.34) với các tham số cho bởi (2.42) 31

2.10 Nghiệm (S, I) của mô hình với điều khiển phản hồi 32

2.11 Nghiệm (u 1, u 2) của mô hình với điều khiển phản hồi 32

MỞ ĐẦU

Các phương trình vi phân đạo hàm thường có một vai trò nổi bật trong

cả lý thuyết lẫn ứng dụng Chúng thường được sử dụng để mô hình hóa mộtcách hiệu quả nhiều hiện tượng và quá trình quan trọng nảy sinh trong thếgiới thực Một cách tổng quát, một phương trình vi phân thường có thể viếtdưới dạng

hàm trơn, tức là khả vi liên tục Việc nghiên cứu các phương trình vi phânthường nảy sinh trong các lĩnh vực ứng dụng có một vai trò đặc biệt quantrọng Chủ đề này đã thu hút được sự quan tâm của nhiều nhà nghiên cứu

ở nhiều lĩnh vực khác nhau trong suốt nhiều năm qua, chẳng hạn, thiết lập

mô hình, nghiên cứu định tính và các lời giải số

Một trong những ứng dụng nổi bật và rất quan trọng của phương trình

vi phân thường là để mô hình toán học các bệnh truyền nhiễm trong dịch tễhọc (epidemic models) Như chúng ta đã biết, dịch tễ học là một môn khoahọc nghiên cứu tình trạng sức khỏe, bệnh tật và các yếu tố liên quan ở cấp

độ dân số Dịch tễ học là khoa học nền tảng của y tế công cộng với vai trò

cơ bản là nâng cao sức khỏe cộng đồng Một số thành tựu quan trọng củadịch tễ học có thể kể đến như thanh toán bệnh đậu mùa, điều trị nhiễm độcMethyl thủy ngân, điều trị bệnh sốt thấp tim và bệnh thấp tim, kiểm soátlây truyền và phòng ngừa các bệnh truyền nhiễm,

Trong toán học, sự lan truyền của nhiều bệnh truyền nhiễm, tiêu biểunhư sởi, quai bị, rubella, thủy đậu, có thể được mô hình thông qua các

hệ phương trình vi phân thường Việc thiết lập mô hình toán học và nghiêncứu các tính chất của chúng giúp chúng ta hiểu rõ các cơ chế lây lan củabệnh dịch, từ đó đề xuất các chính sách hiệu quả để phòng ngừa, kiểm soát

và điều trị bệnh tật Nói riêng, việc nghiên cứu tính chất ổn định tiệm cậntoàn cục của các mô hình dịch tễ học có vai trò đặc biệt quan trọng trongthực tế Bài toán này đã thu hút được sự quan tâm của nhiều nhà toán học,

kỹ thuật, sinh học, dịch tễ học trong suốt nhiều thập kỷ qua Một trongnhững cách tiếp cận thành công nhất tới bài toán này là phương pháp ổnđịnh Lyapunov Phương pháp này nghiên cứu sự ổn định của các hệ độnglực dựa trên việc xác định một hàm số phù hợp, được gọi là hàm Lyapunov.Cho tới nay, sự ổn định của rất nhiều các hệ phương trình vi phân quantrọng nảy sinh trong các lĩnh vực ứng dụng nói chung và các mô hình dịch

tễ học nói riêng đã được thiết lập thành công dựa trên phương pháp ổn địnhLyapunov Các kết quả thu được là rất quan trọng trong cả lý thuyết địnhtính của phương trình vi phân cũng như trong khía cạnh ứng dụng

Chính vì những lý do trên, đề tài "Phương pháp ổn định Lyapunov nghiêncứu sự ổn định toàn cục của một số mô hình dịch tễ học" được thực hiện vớimục tiêu tìm hiểu việc sử dụng phương pháp ổn định Lyapunov cho một số

mô hình dịch tễ học Các ý nghĩa thực tế thu được từ việc thiết lập ổn địnhcho các mô hình này cũng được nghiên cứu và tìm hiểu Ngoài việc tìm hiểu

về mặt lý thuyết, các mô phỏng số cũng được trình bày để minh họa cho cáckết quả lý thuyết

Ngoài phần "Mở đầu", "Kết luận" và "Tài liệu tham khảo", các kết quả

SIRS, SIS và SI với các biến điều khiển phản hồi

Chương 1

Kiến thức chuẩn bị

Nội dung của chương này trình bày một số kiến thức chuẩn bị quan trọngbao gồm bài toán giá trị ban đầu đối với phương trình vi phân cấp một, định

lý ổn định Lyapunov, phương pháp Runge-Kutta bốn nấc kinh điển và một

số mô hình dịch tễ học cổ điển Phần trình bày của chương này được dựatrên các tài liệu [1, 3, 7, 9, 10, 11]

(t, y)trong đótlà một vô hướng thực (real scalar) và giả sử rằng f : D → Rn

là một hàm trơn, tức là khả vi liên tục Một phương trình vi phân thường làmột phương trình có dạng

Khi một phương trình vi phân được sử dụng để mô tả sự tiến hóa củamột biến trạng thái trong một quá trình vật lý, một bài toán cơ bản và quantrọng là xác định các giá trị tương lai của biến trạng thái từ giá trị ban đầucủa nó Trong trường hợp này, mô hình toán học được cho bởi cặp phươngtrình

˙

y = f (t, y), y(t0) = y0, (1.3)trong đó phương trình thứ hai được gọi là điều kiện ban đầu Bài toán (1.3)còn được gọi là bài toán giá trị ban đầu, hoặc bài toán Cauchy

Các vấn đề cơ bản của lý thuyết chung về phương trình vi phân là sự tồntại, tính duy nhất và sự phụ thuộc liên tục vào dữ liệu đầu của bài toán giátrị ban đầu

Mệnh đề 1.1 Bài toán (1.3) là tương đương với phương trình tích phân

miền D = 0 ≤ t ≤ T, |y| < ∞ Giả sử thêm rằng f (t, y) là liên tục

|f (t, y) − f (t, ˆy| ≤ L|y − ˆy|, ∀(t, y), (t, ˆy) ∈ D

Khi đó ta có các khẳng định sau:

đoạn [0, T ] của bài toán (1.3) Hơn nữa, nghiệm này là khả vi

một nghiệm của bài toán nhưng khác giá trị ban đầu thì

|y(t) − ˆy(t)| ≤ eLt|y(0) − ˆy(0)|

3 Tổng quát hơn, nếu yˆ thỏa mãn phương trình với vế phải có nhiễu

dy

dt = f (y), y(t0) = y0

và ta gọi (1.3) là hệ ô-tô-nôm hoặc dừng

Ví dụ 1.1 Phương trình phân rã tuyến tính

˙

y = λy, y(0) = y0,

có nghiệm duy nhất xác định bởi

y(t) = y0eλt

Ví dụ 1.2 Xét mô hình hai loài bao gồm động vật ăn thịt và con mồi Giả

triển của quần thể được mô tả bởi hệ Lotka-Volterra sau đây

˙u = au − buv, ˙v = −cv + duv,

trong đó a, b, c và d là các hằng số dương Trong trường hợp này, chúng takhông thể tìm được nghiệm chính xác của mô hình Tuy nhiên, chúng ta cóthể tìm được nghiệm xấp xỉ cho mô hình nhờ các phương pháp xấp xỉ bao gồmcác phương pháp giải tích và các phương pháp số.

˙z = ˙y = f (z + y∗) = g(z)

Định nghĩa 1.4 ([11]) Một điểm cân bằng y∗ được gọi là ổn định tiệm cận

khả vi liên tục sao cho

Chứng minh Xem ([11]).

triệt để (radially unbounded)

Về mặt lý thuyết ta có thể chỉ ra sự tồn tại và duy nhất nghiệm của bàitoán giá trị ban đầu nhưng nói chung việc tìm nghiệm chính xác là rất khókhăn và phức tạp, thậm chí là không thể Nói chung việc giải gần đúng bàitoán giá trị ban đầu là hầu như không thể tránh khỏi Vì thế các phươngpháp gần đúng bao gồm phương pháp giải tích và các phương pháp số đượcnhiều nhà toán học quan tâm phát triển

Trong các mô phỏng số của luận văn, chúng tôi sử dụng phương phápRunge-Kutta bốn nấc kinh điểm để giải số bài toán giá trị ban đầu Để giảigần đúng bài toán (1.1), đầu tiên chúng tôi sử dụng lưới điểm

t0 < t1 < t2 < < tN = T,

pháp xác định như dưới đây được gọi là phương pháp Runge-Kutta bốn nấckinh điển [3]

y1, y2, Các phân tích toán học chỉ ra rằng phương pháp Runge-Kutta

Trong mục này chúng tôi giới thiệu một số mô hình dịch tễ học cổ điểnbao gồm các mô hình SI, SIR và SIRS

Các bệnh truyền nhiễm như sởi, quai bị, rubella, thủy đậu có thể được

mô hình hóa bằng cách phân loại các cá thể trong cộng đồng tùy theo trạngthái của họ đối với bệnh tật: khỏe mạnh, bị lây nhiễm và miễn dịch Cácbệnh do vi-rút hoặc vi khuẩn gây ra không được mô hình hóa ở cấp độ dân

số mà chỉ gián tiếp thông qua số cá thể bị nhiễm bệnh Các trạng thái bệnh

bệnh nhưng có thể bị nhiễm bệnh và trở thành người truyền bệnh

khả năng truyền bệnh cho người khác

hồi phục và những người có miễn dịch vĩnh viễn, hoặc bị cô lập (khôngthể truyền bệnh cho người khác) cho đến khi hồi phục và miễn dịch vĩnhviễn xảy ra

thể phục hồi và phát triển miễn dịch:

3 SIRS - Chỉ có miễn dịch tạm thời S → I → R → S.

Ví dụ 1.3 (Mô hình SIS đơn giản [1, 10])

tễ SIS và có thể nói đây là mô hình đơn giản nhất trong dịch tễ toán học

Sơ đồ lan truyền của mô hình (1.13) được mô tả như trong Hình 1.1 dướiđây

Hình 1.1: Sơ đồ lan truyền của mô hình SIS đơn giản.

Ví dụ 1.4 (Mô hình SI đơn giản [1, 10])

Chương 2

Tính chất ổn định toàn cục của các

mô hình dịch tễ SIR, SIRS, SIS và

SI với các biến điều khiển phản hồi

Trong chương này trình bày các hàm Lyapunov cho các mô hình dịch tễSIR, SIRS, SIS và SI, tính chất ổn định toàn cục của các điểm cân bằng.Các kết quả lý thuyết được hỗ trợ bởi các mô phỏng số Phần trình bày củachương này được dựa trên các tài liệu [2, 4, 5, 6, 8, 12]

Dựa trên các giả thiết cổ điển đề xuất trong [2, 4, 5], toàn bộ dân số kích

(do điều trị) sẽ kéo theo miễn dịch bền vững hoặc miễn dịch tạm thời Ở

hình dựa trên các giả thiết trên được biết đến là mô hình SIR (có sự miễndịch bền vững) hoặc mô hình SIRS (không có sự miễn dịch bền vững) Sơ đồchuyển trạng thái của mô hình được mô tả trong hình 2.1

Hình 2.1: Sơ đồ lan truyền của mô hình SIS đơn giản.

Dựa trên một số giả thiết cần thiết, chúng ta nhận được mô hình vi phândưới đây

Chứng minh Giả sử (S, I) là một điểm cân bằng của hệ (2.1) Khi đó (S, I)

Giải và biện luận hệ trên ta thu được điều phải chứng minh

Để xây dựng hàm Lyapunov cho hệ (2.1), trước hết chúng ta thực hiệnphép đổi biến

N

R∗0, I

ˆδ

Dễ dàng kiểm tra rằng P∗ và I∗ thỏa mãn

định tiệm cận toàn cục

Định lý 2.2 ([8]) Nếu R0 > 1 thì điểm cân bằng trong E∗ của mô hình(2.1) là ổn định tiệm cận toàn cục.

P N (1 −

P

P0)

2 − ˆγ(1 − R0)I ≤ 0,

là ổn định tiệm cận toàn cục Nghiệm số thu được từ phương pháp

0 200 400 600 800 1000 1200 1400

E*

Hình 2.2: Nghiệm của mô hình trong Ví dụ 2.1.

Ví dụ 2.2 Xét mô hình (2.1) với các tham số cho bởi

γ = 0.2, α = 0.25, β = 0.2, p = 0.25, σ = 0.1, δ = 0.16, N = 1000

(2.13)

ổn định tiệm cận toàn cục Nghiệm số thu được từ phương pháp Runge-Kuttavới h = 10−5 và t ∈ [0, 100] trong trường hợp này được biểu diễn trong Hình

E*

Hình 2.3: Nghiệm của mô hình trong Ví dụ 2.2.

2.2 Mô hình SIS cổ điển

Hình 2.4: Sơ đồ lan truyền của mô hình SIS (2.14).

Dễ dàng kiểm tra được rằng

βS

∗I∗

N = γN + (δ − pγ)I

∗ − σS∗ = (δ + σ +  − pγ)I∗ (2.15)

cận toàn cục

nên theo định lý ổn định toàn cục ta có điều phải chứng minh.

sinh học Tuy nhiên bằng cách tiếp cận tương tự như trong mục 2.1, bằng

Ví dụ 2.3 Xét mô hình (2.14) với các tham số cho bởi

γ = 0.5,  = 0.0001, β = 0.5, p = 0.25, σ = 0.1, δ = 0.1, N = 1000

(2.17)

ổn định tiệm cận toàn cục Nghiệm số thu được từ phương pháp Runge-Kuttavới h = 10−5 và t ∈ [0, 100] được biểu diễn trong trường hợp này được biểu

S

0 500 1000 1500 2000 2500 3000 3500 4000 4500

5000

E*

Hình 2.5: Nghiệm của mô hình trong Ví dụ 2.3.

Ví dụ 2.4 Xét mô hình (2.14) với các tham số cho bởi

γ = 0.6,  = 0.0001, β = 0.6, p = 0.25, σ = 0.1, δ = 0.1, N = 1000

(2.18)

ổn định tiệm cận toàn cục Nghiệm số thu được từ phương pháp Runge-Kuttavới h = 10−5 và t ∈ [0, 100] được biểu diễn trong trường hợp này được biểu

0 1000 2000 3000 4000 5000

E*

Hình 2.6: Nghiệm của mô hình trong Ví dụ 2.4.

Trong mục này, chúng tôi trình bày lớp hàm Lyapunov mới để nghiên cứutính chất ổn định toàn cục của một mô hình SIS với tỷ lệ mắc bệnh chuẩn.Phần trình bày của mục này được dựa trện tài liệu [12]

Để xác định các điểm cân bằng của (2.19), ta xét hệ

Λ − βSI

S + I − µS + φI = 0,βSI

cân bằng dương xác định bởi

2.3.2 Phân tích ổn định

∗

S∗ + I∗.

(2.30)

Tính toán đạo hàm của hàm V dọc quỹ đạo nghiệm của (2.19) ta nhận được

∗ + I∗)

βI∗

I − I∗I

V0(E∗) = 0 Do đó, theo định lý ổn định Lyapunov ta kết luận được E∗ là

ổn định tiệm cận toàn cục Nghiệm của mô hình được biểu diễn trong Hình

0 5 10 15 20 25

E 0

Hình 2.7: Nghiệm của mô hình (2.19) trong Ví dụ 2.5.

Ví dụ 2.6 Xét mô hình (2.19) với các tham số

Λ = 0.25, β = 0.9, α = 0.25, µ = 0.3, φ = 0.25

(2.4358, 0.3072) là ổn định tiệm cận toàn cục Nghiệm của mô hình được

S

0 5 10 15 20 25 30

E*

Hình 2.8: Nghiệm của mô hình (2.19) trong Ví dụ 2.6.

2.4 Các mô hình SI với các biến điều khiển phản hồi

Trong phần này, chúng tôi trình bày một dịch tễ SI với các biến điều khiểnphản hồi được đề xuất bởi Chen và Sun [6] và tính chất ổn định toàn cụccủa mô hình Cùng với đó, các mô phỏng số được trình bày để hỗ trợ chocác kết quả lý thuyết

2.4.1 Mô hình toán học

Xét mô hình dịch tế SI như sau:

˙S(t) = S(t) r − aS(t) − bI(t)
,

˙I(t) = I(t) bS(t) − µ − f I(t)
,

(2.34)

trong đó dân số được chia thành hai nhóm: nhóm những người khỏe mạnh

dương Chi tiết của mô hình được trình bày trong [6]

Để tìm điểm cân bằng dương của (2.34), ta giải hệ

Sử dụng (2.35) cùng một số phép biến đổi đại số ta thu được

˙

V (S, I) = (S − S∗)− a(S − S∗) − b(I − I∗)+ (I − I∗)b(S − S∗) − f (I − I∗)

= a(S − S∗)2 − f (I − I∗)2

định tiệm cận toàn cục Định lý được chứng minh

Do đó, Chen và Sun [6] đã đề xuất một cách tiếp cận hiệu quả để hạn chế

và ngăn chặn tác nhân gây bệnh, đó là mô hình với các điều khiển phản hồi

Trong mục tiếp theo chúng ta sẽ trình bày tính chất ổn định toàn cục của

mô hình (2.37) để chứng minh rằng mô hình là hiệu quả trong hạn chế và

ngăn chặn tác nhân gây bệnh

2.4.2 Phân tích ổn định

Trước tiên ta định nghĩa số

Υ0 = (br − aµ)e1

Để tìm các điểm cân bằng có thể của mô hình (2.37) ta xét hệ

Nhờ các hàm Lyapunov, chúng ta thiết lập tính chất ổn định toàn cục của

cục

Chứng minh (a) Chứng minh (i)