zunia.vn

Tuyển sinh 2024 dành cho Gen-Z

zunia.vn

» Khoa Học Tự Nhiên

Ứng dụng phương trình vi phân có chậm vào mô hình bài toán dân số

Chia sẻ: _ _ | Ngày: | Loại File: PDF | Số trang:5

Báo xấu

47
lượt xem 3
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Bài viết tập trung về phương trình vi phân có chậm thông qua việc nghiên cứu tính ổn định, ổn định tiệm cận, ổn định tiệm cận toàn cục dựa trên định lý hàm Lyapunov bằng phương pháp tính để ứng dụng vào mô hình dân số. Trong bài viết có sử dụng phần mềm Maple để tính các thông số cần thiết cho mô hình, sau đó đưa thông số cần tính vào phần mềm mô phỏng Matlab thông qua việc lập code.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Ứng dụng phương trình vi phân có chậm vào mô hình bài toán dân số

TẠP CHÍ KHOA HỌC ĐẠI HỌC VĂN LANG Lê Nguyễn Hạnh Vy và tgk ỨNG DỤNG PHƯƠNG TRÌNH VI PHÂN CÓ CHẬM VÀO MÔ HÌNH BÀI TOÁN DÂN SỐ APPLYING DELAY DIFFERENTIAL EQUATION IN POPULATION PROBLEM MODEL LÊ NGUYỄN HẠNH VY và TRẦN LƯU CƯỜNG TÓM TẮT: Bài viết tập trung về phương trình vi phân có chậm thông qua việc nghiên cứu tính ổn định, ổn định tiệm cận, ổn định tiệm cận toàn cục dựa trên định lý hàm Lyapunov bằng phương pháp tính để ứng dụng vào mô hình dân số. Trong bài viết có sử dụng phần mềm Maple để tính các thông số cần thiết cho mô hình, sau đó đưa thông số cần tính vào phần mềm mô phỏng Matlab thông qua việc lập code. Kết quả thu được là quỹ đạo nghiệm và biểu đồ Phase của mô hình. Nó cho ta thấy mức độ tăng trưởng dân số phụ thuộc vào thời gian có chậm. Và được đưa vào thực tiễn, tạo ra nhiều bước phát triển mới cho các ngành khoa học khác. Từ khóa: hàm Lyapunov; phương trình vi phân có chậm; mô hình Lotka-Volterra; bài toán phát triển dân số. ABSTRACT: The paper focuses on the application of delay differential equation by digging deeply into stabilization, asymptotic stability, and globally asymptotically stable based on Lyapunov's theorem by calculating method in order to apply in population model. In this paper, we used Maple software to calculate the necessary parameters for the model, and then put the parameters to be calculated in Matlab simulation software through coding. The results are the root locus and the Phase diagram of the model. It indicates that the rate of population growth depends on delay time. And put into practice, creating many new developments for other sciences. Key words: Lyapunov function; delay differential equation; Lotka-Volterra model; population growth problem. 1. ĐẶT VẤN ĐỀ Lý thuyết phương trình vi phân có chậm Nhiều hiện tượng thực tế cuộc sống trong được phát triển rộng rãi bởi Bellman và Cooke vật lý, kỹ thuật, sinh học, y học… có thể được [4], Hale [5], Dirver [6], El’sgol’ts và Norkin mô hình hóa bởi giá trị ban đầu của phương [7] và hiện nay có một cuốn sách mới nói về  . vấn đề này của Hale và Verduyn Lunel [8],  trình vi phân thường:  t )  g (t , x(t )), t  t0 x ( Kolmanowskii và Myshkis [9]… Việc nghiên  x(t 0 )  x0  cứu này yêu cầu đòi hỏi không chỉ về mặt lý Tuy nhiên, để mô hình phù hợp với thực tế thuyết mà cả tính ứng dụng rộng rãi, thu hút hơn, người ta đã sử dụng mô hình hóa bởi được sự quan tâm của nhiều nhà toán học và đã phương trình vi phân có chậm như sau: đưa ra nhiều kết quả quan trọng. Nó đã góp . phần xây dựng lý thuyết chung cho ngành toán x(t )  f (t , xt ), t  t0 học và các ngành khoa học khác. Nó có mặt và góp phần nâng cao tính hấp dẫn, lý thú, tính  ThS. Trường Đại học Văn Lang, lenguyenhanhvy1991@gmail.com  TS. Trường Đại học Văn Lang, cuong.tl@vlu.edu.vn, Mã số: TCKH24-06-2020 86
TẠP CHÍ KHOA HỌC ĐẠI HỌC VĂN LANG Số 24, Tháng 11 – 2020 đầy đủ sâu sắc, tính hiệu quả, giá trị của nhiều Trong đó     C biểu diễn trạng thái ngành như tối ưu, điều khiểu tối ưu, giải tích ban đầu hoặc trạng thái dữ liệu. số, tính toán khoa học... Vì vậy, lý thuyết này 2.2. Tính ổn định Lyapunov trong phương đã trở thành một trong các lĩnh vực toán học trình vi phân có chậm hiện đại nhất, có khả năng ứng dụng trong 2.2.1. Định nghĩa cơ bản về tính ổn định nhiều lĩnh vực khác như: Vật lý học, Cơ học, Lyapunov Kinh tế học, Sinh thái học, Hóa học… Định nghĩa 1: Nghiệm tầm thường của 2. NỘI DUNG phương trình (2) được gọi là ổn định Lyapunov Giới thiệu hàm Lyapunov trong phương khi t   nếu   0, t0  0,     t0 ,    0, a trình vi phân có chậm. Ứng dụng vào mô hình phát triển ổn định dân số. sao cho     xt  t ,     , t  t0 . Bảng ký hiệu: Định nghĩa 2: Nghiệm tầm thường của phương trình (2) được gọi là ổn định đều khi khi t   nếu số  trong định nghĩa không phụ thuộc vào t0 . Định nghĩa 3: Nghiệm tầm thường của phương trình (2) được gọi là ổn định tiệm cận đều theo Lyapunov khi t   nếu i ) Nghiệm tầm thường ổn định đều ii )   0 (không phụ thuộc vào t0 ) 2.1. Giới thiệu phương trình vi phân có   C,     lim x  t0 ,   0 chậm rời rạc t  Dạng tổng quát của phương trình vi phân Định nghĩa 4: Xét phương trình (2), hàm có chậm rời rạc: khả vi liên tục V : R   C  R được gọi là x  t   f  t , x  t  1  , ..., x  t   n   , t  t0 (1) . hàm Lyapunov nếu tồn tại các hằng cố a, b, c >0 thỏa mãn: Với  i  i t  x t    0, t  t0 , i  1, ..., n i) a x  t   V  t , xt   b xt , 2 2 . được gọi là các chậm rời rạc   x t  f t , xt  (2) . ii) V  t , xt   c x  t  2 với mọi nghiệm Trong đó, xt  x  t    ,  r ,0 , là   x(t) của (2) một hàm thuộc không gian các hàm liên tục từ  Định nghĩa 5: Nếu ánh xạ V : R  C  R  r,0 vào R n  . Ký hiệu: C  C0  r ,0 ,R n  liên tục từ x  t0 ,  là nghiệm của phương trình f :   R n là hàm cho trước, với (2) thỏa điều kiện ban đầu  t0 ,   , ta có:   R  C. . . Bài toán giá trị ban đầu:   V  V t ,   lim h0 h   1   V t  h, xt h t0 ,    V t0 ,   . .  x  t   f  t , xt  , t  t0 Hàm V  t0 ,   gọi là đạo hàm trên, bên phải   x0  x  t0        theo t của hàm V  t0 ,  dọc theo nghiệm của hệ (2). 2.2.2. Định lý cơ bản về tính ổn định Lyapunov [1, tr.27] 87
TẠP CHÍ KHOA HỌC ĐẠI HỌC VĂN LANG Lê Nguyễn Hạnh Vy và tgk Cho hàm u  s  , v  s  , w  s  : R   R  liên i (0)  0,i  1,2, 2 (r)  0 với r > 0, tục và không giảm, u  s   0, v  s   0,  s  0 lim 1 (r)   và V : C  R là hàm vi r  và u  0     0     0   0 . Những phát biểu phân vô hướng liên tục. Với tập S của nghiệm . sau là đúng: x (t )  F  xt  thỏa mãn: i) Nếu hàm V : R  C  R thỏa mãn . V    1 (  (0) ),V    2 (  (0) )      V t ,    v    , V t ,       0  u  0 thì Khi đó x = 0 là ổn định tiệm cậm đối với nghiệm x=0 là ổn định đều. tập S, nghĩa là nhứng nghiệm bị chặn trong S ii) Nếu thêm vào (i), lim u  s    thì s  hội tụ đến x = 0 khi t   nghiệm của (2) là bị chặn đều (tức là với bất kỳ Hàm Lyapunov V:   0, tồn tại       0 sao cho, với Hàm Lyapunov V được xác định trong tập hợp   R,   C ,    , ta có  x  , t    , t   . 1 , 2  : 1  C  , 0 , R  , 2     (0),    , 0 : iii) Nếu thêm vào (i), w( s )  0 với s > 0     1 2 V  V 1 , 2  V0 1 (0), 2 (0) thì nghiệm x=0 là ổn định tiệm cận đều. 2 2.3. Ứng dụng phương trình vi phân có   V1 1 (0), 2 (0)  V2 1 , 2   chậm giải các bài toán mô hình phát triển Trong đó, hàm vô hướng V0, V1, V2 được dân số xác định như sau: 2.3.1. Mô hình Lotka – Volterra có chậm đơn V0 (t )  V0 (u1 (t ), u2 (t ))  ln(1  u1 (t ))   ln(1  u2 (t )) Mô hình Lotka – Volterra có chậm đơn: V1  V1 (u1 (t ), u 2 (t ))  u1  z1   (u 2  z 2 ), Xét hệ động vật ăn thịt con mồi Lotka – t t * 2 Volterra có chậm có dạng [3] V2  V2 (1 , 2 )  P  ds  [x (1 1 ( ))1 ( )] d t  s . x  x ( t )( r  ax ( t )  by (t )), t t * 2  Q  ds  [x (1 1 ( ))2 ( )] d . . t  s y  y (t )(  d  cx (t   )), Mô hình: Trong đó: x(t), y(t) là mật độ dân số của con mồi và kẻ săn mồi tại thời điểm t tương ứng. r,a, b,c,d, là các hằng số dương  x0        0 Điều kiện ban đầu:       , 0   x 0 0   y 0  0 Trong đó   C   ,0 , R  , R  x : x  0 và   max     :     ,0 Hình 1. Quỹ đạo nghiệm của mô hình Lotka – Volterra có chậm đơn Bổ đề: Để xây dựng Lyapunov, ta sẽ sử dụng bổ đề sau: Cho 1 (.) và 2 (.) là những hàm vô hướng liên tục không âm sao cho 88
TẠP CHÍ KHOA HỌC ĐẠI HỌC VĂN LANG Số 24, Tháng 11 – 2020 Khi đó x = 0 là ổn định tiệm cận đối với tập S, nghĩa là những nghiệm trong S hội tụ đến x = 0 khi t   . Hàm Lyapunov V: Hàm Lyapunov V được xác định trong: (1 , 2 )  C ( ,0), R 2 ), i (0)  0,  [ ,  max(1 , 2 )], i  1, 2 1 V  V (1 , 2 )  V0 2 (1 (0), 2 (0)) , 2 Hình 2. Biểu đồ phase mô hình Lotka – V1 (1 (0), 2 (0))  V2 (1 , 2 ) Volterra có chậm đơn Trong đó các hàm vô hướng V0 , V1 , V2 2.3.2. Mô hình Lotka – Volterra có chậm kép được xác định như sau: Mô hình Lotka – Volterra có chậm kép [2] V0 (t )  ln(1  u(t )   ln(1  u2 (t )) , Ta xét hệ động vật ăn thịt con mồi Lotka – Volterra với có chậm 1 , 2  R : [0, ) rời V1 (t )  u1  z1   (u2  z2 ), t t rạc riêng biệt V2 (t )  Cu  ds  [1  u2 ( )]2u12 (   2 )d  . t 1 s  x(t )  x(t )[r1  ax(t )  by (t   1 )] . t  A  [1  u (s)] u 2 2  y (t )  y (t )[  r2  cx(t   2 )  dy (t )] 1 1 ( s)ds t 2 Trong đó: r1 , r2 , b, c là hằng số dương và t t . .  B  ds  [1  u1 ( )]2u2 2 (  1 )d a, d là hằng số âm, x(t ) , y (t ) là mật độ dân số t  2 s t 1 của con mồi và kẻ săn mồi tương ứng Cu1  (u  u42 )ds 2 1 Điều kiện ban đầu: t 1 2 ( x( ), y ( ))  (1 ( ), 2 ( ))  0 1 t  ( Du  BD )  [1  u2 ( s)]2u22 ( s)ds   [  max( 1 , 2 ), 0] 2 t 1  (0)  0,  (0)  0  1 t 1  (u 2  B 2 2  u42 )ds . Trong đó t 2 2 2   (1 , 2 )  C([  max(1, 2 ),0], R ) 2  Mô hình: Bổ đề: [2] Để xây dựng hàm Lyapunov, ta sẽ sử dụng bổ đề sau: Cho 1 (.) và  2 (.) là những hàm vô hướng liên tục không âm sao cho  i (0)  0 , i  1, 2; lim  1 ( r )   ,  2 (r )  0 với r  r  0 . Cho V : C  R là hàm vi phân vô hướng liên tục và S là tập con rỗng của C thỏa mãn: V ( )   1 (  (0) , V ( )   1 (  (0) . . Hình 3. Quỹ đạo nghiệm của mô hình Lotka – Volterra có hai sự chậm trễ 89
TẠP CHÍ KHOA HỌC ĐẠI HỌC VĂN LANG Lê Nguyễn Hạnh Vy và tgk 3. KẾT LUẬN Bài viết đã giải quyết được vấn đề được ra về tính ổn định hàm Lyapunov, mô phỏng được quỹ đạo nghiệm bài toán Lotka-Vollterra trong mặt phẳng Phase, sử dụng phần mềm Maple và Malab. Ngoài ra, bằng công cụ này ta sẽ có hướng nghiên cứu tiếp là dùng công cụ giải tích hàm để tìm điều kiện tường minh thay thế phương pháp hàm Lyapunov cho tính ổn định của hệ phương trình vi phân có chậm, như: độ Hình 4. Biểu đồ Phase của mô hình Lotka - Volterra đo ma trận, ma trận Metzler. có hai sự chậm trễ TÀI LIỆU THAM KHẢO [1] Y. Kuang (1992), Delay Differential Equations with Applications in Population Dynamics. [2] Y. Tang, E. Beretta and F. Solimano (2001), Stability analysis of a volterra predator - prey system with two delays, Volume 9, number 1. [3] Y. Kuang, E. Berrtta (1995), Convergence Results in a Well-Known Delayed Predator-Prey Sytem. [4] R. Bellman, K. L. Cooke (1963), Differential – Difference Equation, Academic Press, New York-London. [5] J. K. Hale (1977), Theory of Functiona Differential Equations, Springer-Verlag, New York. [6] R. D. Driver (1963), Existence theory for a delay-differential system, \textit{Contributions to Differential Equations}. [7] L. E. El'sgol'ts và S. B. Norkin (1973), Introduction to the Theory and Application of Differential Equations with Deviating Arguments, Academic Press, New York. [8] J. K. Hale và S. M. Verduyn Lunel (1977), Introduction to Functional Differential Equations, \textit{Applied Mathematical Sciences 99}, Spring-Verlag, New York. [9] V. Kolmanovskii và A. Myshkis (1992), Applied Theory of Functional Diffenrential Equations, Kluwer, Dordrecht. Ngày nhận bài: 23-6-2020. Ngày biên tập xong: 02-11-2020. Duyệt đăng: 27-11-2020 90

CÓ THỂ BẠN MUỐN DOWNLOAD

THÔNG TIN

TRỢ GIÚP

HỖ TRỢ KHÁCH HÀNG

Theo dõi chúng tôi

Chịu trách nhiệm nội dung:

Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA

LIÊN HỆ

Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM

Hotline: 093 303 0098

Email: support@tailieu.vn

Giấy phép Mạng Xã Hội số: 670/GP-BTTTT cấp ngày 30/11/2015 Copyright © 2022-2032 TaiLieu.VN. All rights reserved.