Tính ổn định của phương trình sai phân dạng tuyến tính, tựa tuyến tính và ứng dụng

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN ————– * ————— NGUYỄN HỒNG QN TÍNH ỔN ĐỊNH CỦA PHƯƠNG TRÌNH SAI PHÂN DẠNG TUYẾN TÍNH, TỰA TUYẾN TÍNH VÀ ỨNG DỤNG LUẬN VĂN THẠC SĨ TOÁN HỌC Hà Nội - Năm 2012 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN ————— * —————- NGUYỄN HỒNG QN TÍNH ỔN ĐỊNH CỦA PHƯƠNG TRÌNH SAI PHÂN DẠNG TUYẾN TÍNH, TỰA TUYẾN TÍNH VÀ ỨNG DỤNG Chun ngành: TỐN GIẢI TÍCH Mã số : 60 46 01 LUẬN VĂN THẠC SĨ TOÁN HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC PGS.TS NGUYỄN SINH BẢY Hà Nội - Năm 2012 Mục lục Mở đầu phương trình sai phân 1.1 Các khái niệm 1.2 Phương trình sai phân vơ hướng 1.2.1 1.2.2 1.3 Phương trình sai phân tuyến tính tr 1.3.1 1.3.2 1.3.3 Tính ổn định phương trình sai phân 2.1Khái niệm ổn định nghiệm phương t 2.2Phương pháp bất đẳng thức ngh Định tính vài mơ hình dạng sai phân 3.1 Mơ hình Cobweb (về thị trườ 3.2 Tăng trưởng GDP (Gross 3.3 Mô hình quần thể cạnh tranh Kết luận Tài liệu tham khảo i Mở đầu Chúng ta làm việc nhiều với trình thời gian liên tục Trên loại thời gian có lý thuyết phương trình vi phân Trong luận văn này, kiến thức trình bày nghiên cứu thực loại thang thời gian khác, gọi thang thời gian rời rạc Thực tế cho thấy phần lớn liệu thường lưu giữ xử lý với trình thời gian Tập thời điểm rời rạc đơn giản, phổ biến tiện lợi sử dụng tập thời điểm cách độ dài h > 0; thời điểm t0 I := ft0 + kh : k = 0; 1; 2; :::g Tập thời điểm quy tập số nguyên: Z := f0; 1; 2; :::g, ta gọi đơn giản lưới số nguyên Tương tự khái niệm đạo hàm trình thời gian liên tục ta có khái niệm sai phân cấp khái niệm phương trình sai phân Mục tiêu luận văn tìm cách giải số lớp phương trình đơn giản, nghiên cứu định tính phương trình sai phân cuối tìm vài ứng dụng thơng qua mơ hình cụ thể thực tiễn Luận văn cấu trúc thành ba chương sau: Chương trình bày kiến thức tổng quan phương trình sai phân vài ứng dụng trực tiếp Chương trình bày tính ổn định phương trình sai phân, phương pháp nghiên cứu tính ổn định Chương Trình bày ứng dụng lý thuyết định tính phương trình sai phân để nghiên cứu định tính vài mơ hình dạng sai phân Luận văn thực khoa Toán - Tin - Cơ học, trường Đại học Khoa học Tự nhiên Hà Nội hướng dẫn khoa học PGS TS Nguyễn Sinh Bảy Nhân dịp tác giả muốn bày tỏ lòng biết ơn chân thành kính trọng tới PGS hướng dẫn tận tình cho tác giả suốt q trình hồn thành luận văn này, từ việc định hình luận văn, hướng dẫn đọc tài liệu, đầu ví dụ, kiểm tra kiến thức khuyến khích động viên tác giả gặp khó khăn nghiên cứu Tác giả xin chân thành cảm ơn Ban Giám hiệu, phịng Sau Đại học, khoa Tốn - Cơ - Tin học trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội tạo điều kiện thuận lợi cho tác giả suốt trình học tập trường Tác giả xin chân thành cảm ơn lãnh đạo thầy cô, đồng nghiệp trường MỞ ĐẦU THPT Mai Châu, Huyện Mai Châu, Tỉnh Hịa Bình - nơi tác giả cơng tác gia đình, người thân bạn bè ln tạo điều kiện, động viên, khuyến khích tác giả q trình học tập hồn thành luận văn Mặc dù tác giả có nhiều cố gắng luận văn khơng tránh khỏi nhiều thiếu sót Tác giả kính mong rộng lượng tha thứ xin tiếp thu ý kiến góp ý từ Thầy, Cô Bạn Xin chân thành cảm ơn! Hà Nội, tháng 12 năm 2012 Nguyễn Hoàng Quân Bảng ký hiệu r(A) - tập giải toán tử tuyến tính A s (A) - tập phổ tốn tử tuyến tính A F(n; m) - ma trận hệ K - lớp hàm + Hahn R := [0; +¥) + Z := f0; 1; 2; 3; :::g Z := f0; 1; 2; 3; :::g d R Không gian véc tơ d- chiều Z(n0) := fn0; n0 + 1; n0 + 2; :::g Z(m; n) := fm; m + 1; m + 2; :: : ng k D x(n) - sai phân bậc k hàm x(:) n Chương Mở đầu phương trình sai phân 1.1 Các khái niệm Lưới Z sai phân Cho điểm t0 R khoảng cách h : < h < +¥ Tập I = ft0 + nh : n = 0; 1; 2; ::::g: gọi lưới thời gian rời rạc cách với bước lưới h > 0, thời điểm t0 R Trường hợp đặc biệt: Nếu lấy t0 = coi h = đơn vị thời gian tập I trở thành tập số nguyên Z I = f0 + n : n = 0; 1; 2; ::::g := Z: Trường hợp riêng: với n = 0; 1; 2; ::: ta có tập số nguyên không âm: I + = f0; 1; 2; 3; :::g := Z : Kí hiệu: + R = [0;+¥): Z(n0) = fn0; n0 + 1; n0 + 2; :::; g (n0 Z): Z(m; n) = fm; m + 1; m + 2; :::; n 1; ng (m < n): Chương Mở đầu phương trình sai phân k Giả sử f ánh xạ từ Z vào R (hoặc vào không gian tổng quát X): f :Z ! R k k Z 3n ! f (n) R : k Định nghĩa 1.1.1 Giả sử f ( ) hàm số xác định tập Z, nhận giá trị R Khi đó, hiệu sau gọi sai phân cấp hàm f ( ) n Z: Sai phân cấp hai D Sai phân cấp k D k Tính chất thường dùng sai phân cấp: (xem [3,6]): DC = (C số) k m D x = : < đa thứ k N k k D [ax(n) + b y(n)] = aD x(n) + b D y(n) (a; b R) k å D x(n) = D n=M k x(N + 1) D k x(M) Khái niệm phương trình sai phân Trong luận văn, ký hiệu x(n) xn hiểu Định nghĩa 1.1.2 Giả sử x(n); n Z hàm số chưa biết cần tìm từ đẳng thức n; x(n) sai phân cấp đến cấp k x(n), thiết phải k có mặt D x(n), đẳng thức gọi phương trình sai phân cấp k Nói cách khác, phương trình sai phân cấp k với hàm cần tìm x(n) đẳng thức có dạng sau: k F(n; D x(n); D k x(n); :::; Dx(n); x(n)) = 0: Chương Mở đầu phương trình sai phân Từ Định nghĩa 1:1:1 sai phân cấp, ta thấy phương trình sai phân cấp k đưa dạng tương đương sau (không khuyết x(n) x(n + k)): G(n; x(n + k); x(n + k 1); :::; x(n + 1); x(n)) = 0: Trường hợp riêng sau (1:5) gọi phương trình sai phân cấp k dạng tắc x(n + k) = f (n; x(n + k 1); x(n + k 2); :::; x(n + 1); x(n)): Phương trình sai phân cấp k có dạng sau gọi phương trình tuyến tính cấp k x(n + k) + ak 1(n)x(n + k Nếu f (n) ta có phương trình sai phân tuyến tính x(n + k) + ak 1(n)x(n + k Nếu hệ số ai(n) khơng phụ thuộc vào n ta có phương trình sai phân tuyến tính hệ số x(n + k) + ak 1x(n + k 1) + + a1x(n + 1) + a0x(n) = 0: Trong trường hợp phương trình tuyến tính hệ số x R , phương trình nghiệm phức sau gọi phương trình đặc trưng phương trình trên: k l + ak 1l k + + a1l + a0 = (l C): Một vài tính chất phương trình sai phân tuyến tính Với phương trình sai phân (các cấp) ta có khái niệm nghiệm tổng quát, nghiệm riêng tương tự với phương trình vi phân Nghiệm tổng quát phương trình sai phân cấp k chứa k số tuỳ ý C 1;C2; :::;Ck: x(n) = f (n;C1;C2; :::;Ck): 0 Với giá trị cụ thể C1 ;C2 ; :::;Ck , ta có nghiệm riêng 0 x(n) = f (n;C1 ;C2 ; :::;Ck ): 0 Thông thường nghiệm riêng xác định theo điều kiện ban đầu: Cho trước (n 0; x1 ; x2 ; :::; xk ), nói nghiệm x(:) phương trình sai phân cấp k nói thỏa mãn điều kiện 0 ban đầu (n0; x1 ; x2 ; :::; xk ) nếu: 0 x(n0) = x1 ; x(n0 1) = x2 ; :::; x(n0 k + 1) = xk0: Một vài tính chất tập nghiệm: Chương Mở đầu phương trình sai phân Nếu x1(n) x2(n) nghiệm riêng (1:8) với số a; b có x(n) = ax1(n) + b x2(n) nghiệm riêng (1:8) i) Nếu x1(n); x2(n), , xk(n) nghiệm riêng độc lập tuyến tính (1:8) nghiệm tổng quát (1:8) ii) x(n) = C1x1(n) +C2x2(n) + +Ckxk(n) với C1;C2; :::;Ck số tùy ý Nếu x(n) nghiệm tổng quát (1:8) xˆ(n) nghiệm riêng (1:7) x(n) = x(n) + xˆ(n) nghiệm tổng quát (1:7) iii) iv) Nguyên lý chồng chất nghiệm: Giả sử x1(n); x2(n) tương ứng nghiệm riêng hai phương trình: x(n + k) + ak x(n + k) + ak x(n) = x1(n) + x2(n) nghiệm riêng phương trình x(n + k) + ak 1(n)x(n + k 1) + + a1(n)x(n + 1) + a0(n)x(n) = f1(n) + f2(n): Phương trình sai phân phi tuyến dạng tắc Một phương trình sai phân tuyến tính cấp k dạng tổng qt khơng đưa dạng tắc dạng tuyến tính nói chung ta chưa có cách giải Mọi phương trình dạng tắc giải cách truy hồi Phương trình sai phân tắc cấp k (1:5) (trong khơng gian X đó) thường viết theo cách sau x(n + 1) = f (n; x(n); x(n 1); :::; x(n k + 1)): Định lý tồn nghiệm phương trình khơng địi hỏi tính liên tục, tính Lipschitz hàm f Điều đơn giản so với trường hợp 0 phương trình vi phân Với điều kiện ban đầu x(n 0) = x1 ; x(n0 1) = x2 ; :::; x(n0 k + 1) = xk , việc tìm cơng thức nghiệm riêng, thỏa mãn điều kiện ban đầu thực cách truy hồi liên tiếp n0 sau: x(n0 + 1) = f (n0; x(n0); x(n0 1); :::; x(n0 k + 1)) = 0 f (n0; x1 ; x2 ; :::; xk ): x(n0 + 2) = f (n0 + 1; x(n0 + 1); x(n0); :::; x(n0 = 0 k + 2)) x(n0 + 3) = f (n0 + 2; x(n0 + 2); x(n0 + 1); :::; x(n0 = 0 f (n0 + 1; f (n0; x1 ; x2 ; :::; xk ); x2 ; :::; xk ::: k + 3)) 1): Chương Định tính vài mơ hình dạng sai phân Ký hiệu lại A = tắc: Giá cân P xác định Ta xét định tính nghiệm cân ba trường hợp: mS (VìA= a) Khi m < A < , ta có: Vậy điểm cân P ổn định tiệm cận (Hình H.2) b) Khi A = D 1, ta có phương trình P(n + 1) = P(n) + B Lấy điều kiện ban đầu: P(0) = P0; ta có: P(1) =P(0)+B P(2) =P(1)+B = ( P0+B)+B = P0 33 Chương Định tính vài mơ hình dạng sai phân P(n) = ( B Giá dao động quanh giá trị cân P = (hình H.3) c) Khi A < 1, ta có n P(n) = A P0 + Điểm cân P không ổn định (H.4) 34 Chương Định tính vài mơ hình dạng sai phân 3.2 Tăng trưởng GDP (Gross domestic product) + Lưu ý ta xét trình với biến thời gian n Z Đơn vị thời gian tính theo năm Ta cần số ký hiệu: Y(n) - tổng thu nhập quốc nội (của quốc gia) năm thứ n G (n) - tổng chi dùng phủ năm thứ n C(n) - tổng chi dùng dân sinh (các hộ gia đình, doanh nghiệp tư nhân) năm thứ n I(n) - tổng đầu tư năm thứ n Trong năm tổng GDP quốc gia lấy tổng ba đại lượng: tổng chi dùng Chính phủ, tổng chi dùng dân sinh tổng đầu tư: Y (n) = C(n) + I(n) + G(n) Giả sử hàm chi dùng dân sinh dự kiến sở từ tổng thu nhập năm trước liền kề: đó, a (0; 1), gọi xu hướng chi dùng cận biên Giả sử lượng đầu tư tỷ lệ với độ tăng lên hàm chi dùng dân sinh: số b > gọi số quan hệ Giả sử lượng chi dùng Chính phủ g > 0; không thay đổi theo năm Từ (3.5), (3.6), (3.7) ta có phương trình Y (n + 2) Phương trình đặc trưng (3.8) là: Điểm cân (3.8) là: Gọi y (n) nghiệm tổng quát phương trình Y (n + 2) 35 Chương Định tính vài mơ hình dạng sai phân Do y nghiệm riêng phương trình khơng (3.8) nên nghiệm tổng quát (3.8) là: y(n) = Mọi nghiệm phương trình khơng (3.8) dần tới y nghiệm phương trình đặc trưng (3:9) nằm phía hình trịn đơn vị phức, nghĩa maxfjl1j; jl2jg < 1; li s (A) g Định lý 3.2.1 Điểm cân y = y(n + 2) a phương trình (3.8) (hay mơ hình GDP) a(1 + b )y(n + 1) + ab y(n) = g ổn định tiệm cận < > > a Chứng minh Điều kiện đủ Xét phương trình: y(n +2) a(1 +b )y(n +1) +ab y(n) = g : Phương trình nhất: y(n + 2) Phương trình đặc trưng: P(l ) = l P( Ở a0 = + a + 2ab ; a1 = 2(1 Ma trận Hurwitz Điều kiện đủ để nghiệm y = hệ (3.10) ổn định tiệm cận D01 >a < > : 36 D2 > > Chương Định tính vài mơ hình dạng sai phân Điều kiện cần Để hệ số đơn giản, trước tiên ta xét cho phương trình tổng qt, sau thay hệ số cụ thể mơ hình Xét phương trình tổng quát: x(n + 2) + p1x(n + 1) + p2x(n) = 0: Phương trình đặc trưng : l + p1l + p2 = 0: Giả sử nghiệm x (3.12) ổn định tiệm cận 2 Ta xét hai khả sau: D = p1 4p2 D = p1 a) Khi D = p1 4p2 0; phương trình đặc trưng có hai nghiệm thực l 1;2 = Do hệ ổn định tiệm cận nên ta có: jl1j < jl0j < Ta thấy, kl1k < , p1 Đầu tiên, ta xét bất đẳng thức bên phải: ( q p1 + > p12 4p2 < p12 kl2k < , p1 Ta xét bất đẳng thức bên trái (3.16): p1 p1 < ( + p1 4p1 > p1 4p2 Do p1 < nên bất đẳng thức bên trái (3.14) hiển nhiên thỏa mãn Từ (3.15), p1 > Kết hợp (3.15)và (3.17), ta được: Lại p1 4p2 nên p2 Tóm lại, hệ (3.12) ổn định tiệm cận ta có: 37 Chương Định tính vài mơ hình dạng sai phân < > > p2 < 1 + p1 + p2 > : (3.17) p + p2 > Áp dụng (3.17) cho mơ hình: y(n + 2) a(1 + b )y(n + 1) + ab y(n) = 0, ta được: Nghiệm y(n) ổn định tiệm cận < > ab < 1 a>0 : + a + 2ab > > Trường hợp a) chứng minh xong b) Khi D = p hợp 4p2 < 0; phương trình đặc trưng có hai nghiệm phức (thực sự) liên l 1;2 = Do hệ ổn định tiệm cận nên ta có: Từ (3.19) (3.20) ta có hệ phương trình sau: ( p12 < p2 < Tiếp theo, ta p1 + p2 > 0: p Thật vậy, từ < p2 < p1 < Xét hàm f (x) = + x Vậy f (x) > 0; 38 Chương Định tính vài mơ hình dạng sai phân Tóm lại, từ giả thiết nghiệm tầm thường y = ổn định tiệm cận ta được: >0 < p2 < 1; > > < > > + p1 + p2 > 0; > :1 p1 + p2 > 0: Điều kiện cần chứng minh tương tự cho trường hợp trước 3.3 Mơ hình quần thể cạnh tranh lồi Giả sử ta có mơi trường khép kín ta quan tâm đến loài Gọi u(k) lượng cá thể lồi thời kỳ k (k Z ) Xét phương trình u(k + 1) = Trong đó: a số gọi số tăng trưởng riêng loài (a > 0), b > gọi hệ số cạnh tranh lồi Khi b = phương trình tăng trưởng khơng có cạnh tranh Đó phương trình tăng trưởng Malthus u(k + 1) = au(k): Với b > phương trình (3.21), bu(k) nhỏ (so với 1) mẫu số gần với 1, phương trình tăng trưởng gần với phương trình Mathus (có hình tựa mũ) Khi bu(k) gần so với 1thì thành phần đáng kể, làm cho đà tăng trưởng bị chậm lại (mẫu số lớn 1) Bây ta xét tính chất mơ hình cho phương trình (3.21): ) Điểm cân bằng: Gọi u điểm cân (3.21), ta xác định u từ phương trình: u = + bu , " Khi a = hai điểm cân 0, lồi diệt vong nên ta không xét tiếp trạng thái tầm thường Ta quan tâm đến điểm cân u = a b với a 6= 1; b > a Để điểm cân u = b > 0, ta giả thiết a > chứng minh a trường hợp này, nghiệm u = b > ổn định tiệm cận 39 Chương Định tính vài mơ hình dạng sai phân ) Tiếp theo, cách đổi biến, ta đưa điểm cân u điểm cân 0: Đặt z = u u , ta đưa phương trình (3.21) phương trình Khi đó, điểm cân u = ) Để xét tính ổn định tiệm cận nghiệm cân z = phương trình (3.22), ta chọn hàm Lyapunov (ở lân cận đủ nhỏ z = 0): V (z) = ln (1 + z) Đầu tiên, ta tồn hàm a(:) K ; b(:) K , cho: a(kzk) V (z) b(kzk); 8z thuộc lân cận đủ nhỏ z = 0: Ta thấy, kzk đủ bé cho: e Vậy z = điểm cực tiểu hàm f (x) Do tồn lân cận đủ nhỏ z = cho ln có f (z) > (trừ điểm z = 0), hay ln(1 + z) < z; 8z 6= 0; đủ nhỏ: ) Với z < (đủ gần 0), đặt t = z Khi đó, kln(1 + z)k < kzk trở thành kln(1 t)k < t(t > đủ nhỏ) Do t < đủ gần (sao cho kln(1 t) < 0k), bất phương trình trở thành: ln(1 t) < t (t > 0): Xét f (t) = t + ln(1 t), ta có: f (0) = 0; f (0) = 0; f 00 (0) = > Vậy t = điểm cực tiểu địa phương ) f (t) > hay t > kln(1 t)k hay kzk > kln(1 + z)k z < đủ gần 0: 40 Chương Định tính vài mơ hình dạng sai phân Tóm lại, ta có hàm b(kzk) = kzk Tiếp theo, ta tìm hàm a(:) K cho a(kzk) V (z) Xét lân cận U(0) = f1 < z < e 1g Khi đó, z U(0) kln(1 + z)k < 4 Nếu chọn a(kzk) = kzk = z , ta có: a(kzk) < V (z) Thật vậy, Xét hàm f (z) = z ln (1 + z): Khi ta có : f (0) = 0; f (0) = 0; f 00 (0) = < Nên z = điểm cực đại Vậy tồn lân cận z = cho f (z) < f (0) = hay z V (z) < , z < V (z) Vậy, chọn hàm a(kzk) = kzk ; a(:) K Cuối cùng, ta tồn lân cận z = cho DV (z) < (tìm a; b để DV (z) < lân cận z = ) Nhắc lại phương trình: z(k + 1) = DV (z) = ln (1 + Lưu ý ta xét cho a > Khi đó, tồn kzk đủ nhỏ để a + bz > Vì vâỵ, với z > ta có: a + bz > ) < + z < + z ) ln2(1 + z ) < ln2(1 + z ) hay DV < a+bz a+bz z z z z a > ) 9zja+bz > ) > 1+ a+ bz > 1+ ) ln (1+ a+ bz ) < ln (1+ ) hay DV < 0: Vậy V (z) hàm Lyapunov thỏa mãn điều kiện Định lý 2.13 Nghiệm tầm thường z = ổn định tiệm cận Từ đó, đổi biến ngược lại, ta có: a u(k) u = b ổn định tiệm cận Tóm tắt chương Chương trình bày ứng dụng lý thuyết định tính để khảo sát tính chất vài mơ hình sai phân thực tế Các mơ hình thầy hướng dẫn lập mơ hình tác giả thực việc xét tính ổn định, tính dao động 41 Kết luận Luận văn cho tổng quan ngắn gọn khái niệm phương trình sai phân, cách k giải số lớp phương trình đơn giản R R Tiếp theo, luận văn trình bày vắn tắt kiến thức tính ổn định nghiệm phương trình sai phân sâu vào phương pháp bất đẳng thức nghiên cứu tính ổn định Cuối cùng, luận văn trình bày cách ứng dụng lý thuyết định tính (tính ổn định, tính dao động) vào việc phân tích số mơ hình thực tiễn 42 Tài liệu tham khảo [1] Nguyễn Thế Hồn Phạm Phu, Cơ sở phương trình vi phân lý thuyết ổn định, NXB Giáo dục (2000) [2] Vũ Ngọc Phát, Nhập môn Lý thuyết điều khiển Toán học, NXB Đại học Quốc gia Hà nội (2001) [3] Lê Đình Thịnh, Đặng Đình Châu Lê Đình Định, Phương trình sai phân số ứng dụng NXB ĐHQG Hà Nội (2001) [4] E Liz and J.B Ferreiro, A Note on the Global Stability of Generalized Differ-ence Equations, Applied Mathematics Letters 15 (2002) 655-659 [5] Eduardo Liz, Stability of non-autonomous difference equations: simple ideas lead-ing to useful results, Journal of Difference Equations and Applications , Vol 17, (2011), 203–220 [6] R.P Agarwal, Diference Equations and Inequalities, Theory Methods and Appli-cations, Marcel Dekker New York (2000) [7] S Udpin and P Niamsup, New discrete type inequalities and global stability of nonlinear difference equations, Appl Math Letters (2009), 22, 856859 [8] M.I Gil, Diference Equations in Normed Spaces, Stability and Oscilations, North Holland (2006) [9]N S Bay and V N Phat, Stability analysis of nonlinear retarded difference equations in Banach spaces, J Comp and Math with Appl., 45 (2003), 951-960 [10] N S Bay, Stability and stabilization of nonlinear time-varying delay systems with non-autonomous kernels, Adv in Nonl Var Ineq., Volume 13, 2, (2010), 59-69 43 ... tính ổn định phương trình sai phân, phương pháp nghiên cứu tính ổn định Chương Trình bày ứng dụng lý thuyết định tính phương trình sai phân để nghiên cứu định tính vài mơ hình dạng sai phân Luận... đầu phương trình sai phân 1.1 Các khái niệm 1.2 Phương trình sai phân vơ hướng 1.2.1 1.2.2 1.3 Phương trình sai phân tuyến tính tr 1.3.1 1.3.2 1.3.3 Tính ổn định phương trình sai phân. .. Chương trình bày khái niệm sai phân, phương trình sai phân, cơng thức nghiệm tổng qt phương trình sai phân tuyến p tính R ; R mối liên hệ qua lại chúng 18 Chương Tính ổn định phương trình sai phân