ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA TOÁN - CƠ - TIN HỌC PHẠM KIM QUÝ BÁN KÍNH ỔN ĐỊNH CỦA HỆ PHƯƠNG TRÌNH VI PHÂN ĐẠI SỐ CÓ CHẬM LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội - 2015 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA TOÁN - CƠ - TIN HỌC PHẠM KIM QUÝ BÁN KÍNH ỔN ĐỊNH CỦA HỆ PHƯƠNG TRÌNH VI PHÂN ĐẠI SỐ CÓ CHẬM LUẬN VĂN THẠC SĨ KHOA HỌC Ngành: Toán ứng dụng Mã số: 60460112 Cán hướng dẫn: PGS TSKH Vũ Hoàng Linh Hà Nội - 2015 LỜI CẢM ƠN Trước trình bày nội dung luận văn, em xin bày tỏ lòng biết ơn sâu sắc tới PGS TSKH Vũ Hoàng Linh, người dành nhiều thời gian, công sức để hướng dẫn tận tình bảo suốt trình thực luận văn Nhân em xin gửi lời cảm ơn đến ban lãnh đạo thầy cô giáo, anh/chị cán trường ĐHKHTN - ĐHQGHN nói chung khoa Toán - Cơ - Tin học nói riêng tạo điều kiện thuận lợi nhất, giúp đỡ em thời gian em học tập, nghiên cứu trường Tôi xin cảm ơn anh chị bạn chuyên ngành Toán ứng dụng động viên ý kiến trao đổi quý báu thân thời gian qua Cuối xin bày tỏ lòng biết ơn gia đình, người thân chỗ dựa tinh thần vật chất sống học tập Hà Nội, ngày 01 tháng 05 năm 2015 Học viên Phạm Kim Quý DANH MỤC KÍ HIỆU VÀ CHỮ VIẾT TẮT • AC([0, ∞), Cn): Không gian hàm liên tục tuyệt đối từ [0, ∞) vào n C • C: Tập số phức • Ckw(I, Cn): Không gian hàm khả vi liên tục khúc cấp k từ I p vào Cn • diag(σ1, σ2): Ma trận đường chéo với phần tử chéo σ1, σ2 • K: K = R K = C • PTVP: Phương trình vi phân • PTVP ĐS: Phương trình vi phân đại số • R: Tập số thực • rank A: Hạng ma trận A Mục lục Lời cảm ơn Mở đầu Một số kiến thức chuẩn bị 1.1 Một số khái niệm ma trận 1 Ma trận Metzler, ma trận dương Nghịch đảo Drazin 61 Khai triển kì dị 1 Phổ số 10 1.2 Chuẩn véc-tơ chuẩn ma trận 12 1.3 Một số khái niệm phương trình vi phân 13 Một số kết bán kính ổn định 16 2.1 Bán kính ổn định hệ phương trình vi phân đại số 16 2.2 Bán kính ổn định phương trình vi phân thường có chậm 27 2 Bán kính ổn định PTVP thường có chậm 28 2 Hệ dương có chậm 30 Bán kính ổn định hệ phương trình vi phân đại số có chậm 33 3.1 Các khái niệm mở đầu 34 3.2 Tính ổn định mũ phương trình vi phân đại số có chậm 39 3.3 Tính ổn định mũ vững 44 Tài liệu tham khảo 58 Mở đầu Bài toán bán kính ổn định PTVP ĐS có chậm (Delay Differential Algebraic Equations) toán nghiên cứu tính ổn định, ổn định vững xây dựng công thức tính toán bán kính ổn định thực/phức cho PTVP ĐS có chậm, dạng: Ex(t) = Ax(t) + Dx(t − τ ), ˙ E, A, D ∈ Cn⋅n, x : I → Cn, I = [0, ∞), τ > độ trễ thời gian, det E = Trong tài liệu này, tính chất Π hệ gọi vững tính chất bảo toàn nhiễu tùy ý ε (đủ nhỏ) tác động lên hệ Ngoài việc quan tâm tới tính vững tính chất, người ta quan tâm tới độ vững tính chất mà đại lượng quan trọng để đánh giá khái niệm bán kính thuộc tính (được đo mê-tric tương thích) Trong khuôn khổ luận văn, tính chất Π xét tính ổn định, hệ xét hệ PTVP ĐS có chậm tuyến tính hệ số hằng, chịu tác động nhiễu có cấu trúc PTVP ĐS có chậm trường hợp tổng quát PTVP ĐS (Differential Algebraic Equations) PTVP thường có chậm (Delay Ordinary Differential Equations) Trong PTVP ĐS mô hình toán học cho nhiều hệ động lực nhiều lĩnh vực ứng dụng chẳng hạn mô mạch điện, lý thuyết điều khiển, động lực chất lỏng, kĩ thuật hóa học, PTVP ĐS có chậm cần thiết để mô hình hóa tác động không tức thời (có chậm) Không giống trường hợp PTVP thường có chậm PTVP ĐS, việc nghiên cứu tính ổn định PTVP ĐS có chậm gặp nhiều khó khăn bao gồm phần ràng buộc đại số độ trễ thời gian, chí lý thuyết tồn nghiệm thu nhiều kết Khó khăn rõ rệt phân tích tính ổn định Hầu hết kết biết tính ổn định PTVP ĐS có chậm trường hợp quy có hệ số vài trường hợp có dạng đặc biệt Nhiều kết biết PTVP thường có chậm PTVP ĐS chuyển sang PTVP ĐS có chậm Bài báo [5] sở thực luận văn Trong tài liệu này, tác giả nghiên cứu tính ổn định hệ thông qua mối quan hệ tập phổ với − tập C với số điều kiện kèm theo Và để thu công thức tính toán bán kính ổn định PTVP ĐS có chậm, việc phân tích phức tạp với việc sử dụng hàm truyền G(λ)(transf erf unctions) Trong luận văn, tác giả đề cập đến dạng PTVP tuyến tính hệ số Luận văn gồm 56 trang, phần mở đầu, kết luận tài liệu tham khảo, luận văn gồm có ba chương: ⋄ Chương Một số kiến thức chuẩn bị Trong chương này, tóm tắt số kiến thức sử dụng luận văn, chủ yếu kiến thức mở rộng ma trận, véc-tơ chuẩn ⋄ Chương Một số kết bán kính ổn định Nội dung chương giới thiệu số kết công thức bán kính ổn định PTVP ĐS PTVP thường có chậm tuyến tính hệ số trường hợp đặc biệt phương trình vi phân đại số có chậm ⋄ Chương Bán kính ổn định hệ phương trình vi phân đại số có chậm Chương nội dung luận văn Trong đó, phân tích chứng minh kết bán kính ổn định phức PTVP ĐS có chậm tuyến tính hệ số Và kết đưa công thức tính toán bán kính ổn định Quá trình tìm hiểu, nghiên cứu vấn đề không tránh khỏi sai sót, hạn chế Do đó, em mong nhận góp ý thầy cô để luận văn hoàn chỉnh Chương Một số kiến thức chuẩn bị Trong chương này, trình bày tóm tắt số kiến thức cần dùng cho phân tích, chứng minh luận văn vài ví dụ minh họa Cụ thể số kiến thức mở rộng ma trận, chuẩn vài kiến thức PTVP 1 1.1.1 Một số khái niệm ma trận Ma trận Metzler, ma trận dương Nghịch đảo Drazin Định nghĩa 1.1 Cho ma trận A = [aij] ∈ Rn⋅n, ≤ i, j ≤ n Khi đó: A gọi ma trận Metzler tất phần tử, ngoại trừ phần tử đường chéo chính, không âm, tức aij ≥ 0, ∀i = j A gọi ma trận không âm (nonnegative matrix) viết A ≥ aij ≥ 0, ∀i, j = 1, 2, , n A gọi ma trận dương (positive matrix) tất phần tử A dương, tức aij > 0, ∀i, j = 1, 2, , n, kí hiệu A > Trong đại số tuyến tính, biết đến khái niệm ma trận nghịch đảo ma trận vuông khả nghịch Mở rộng khái niệm có khái niệm nghịch đảo khác, chẳng hạn: nghịch đảo nhóm, nghịch đảo Drazin, nghịch đảo suy rộng Trong phần trình bày khái niệm nghịch đảo Drazin vài kết liên quan Định nghĩa 1.2 Cho ma trận A ∈ Cn⋅n Khi đó: Số tự nhiên k gọi số A kí hiệu ind(A) = k k số tự nhiên nhỏ thỏa mãn rank Ak = rank Ak+1 Ma trận X ∈ Cn⋅n gọi nghịch đảo Drazin A X thỏa mãn đồng thời biểu thức Ak X A = Ak , X AX = X , AX = X A Trong đó, k = ind(A) Nghịch đảo Drazin ma trận A kí hiệu AD Từ định nghĩa ta có rằng, khái niệm nghịch đảo thông thường trường hợp đặc biệt nghịch đảo Drazin, tức A khả nghịch − theo nghĩa thông thường AD = A Ta có số kết sau nghịch đảo Drazin Định lý 1.3 Trong định lý ta xét ma trận vuông Khi ta có khẳng định sau: (a) Nghịch đảo Drazin ma trận A tồn nhất, (b) Nghịch đảo Drazin ma trận lũy linh ma trận không, (c) Nếu P ma trận chiếu, P = P , có số ind P ≤ P D = P , ∗ ∗ (d) (A )D = (AD) , (e) (AT )D = (AD)T Ví dụ sau ma trận nghịch đảo Drazin ma trận suy biến Ví dụ 1.1 Xét ma trận: 100 A= 0 000 Ta có rank A = 2, rank A2 = rank A3 = nên ind(A) = − Vì det A = 0, nên không tồn A Tuy nhiên ta kiểm tra 100 X= 0 000 thỏa mãn điều kiện Định nghĩa 1.2, tức AD = X 1.1.2 Khai triển kì dị Khai triển kì dị (Singular Value Decomposition) công cụ đại số tuyến tính mạnh hữu dụng, sử dụng nhiều toán liên quan đến ma trận mà áp dụng phương pháp khử Gauss hay phân tích LU cho kết với sai số lớn Phân tích SVD dựa định lý sau, xem [6] Định lý 1.4 Cho A ∈ Cm⋅n Khi tồn ma trận trực giao U ∈ Cm⋅m, V ∈ Cn⋅n ma trận đường chéo D := diag(σ1, , σr) σi, ≤ i ≤ r, bậc hai dương (kể bội) giá trị riêng ∗ ma trận A A thỏa mãn ∗ A=UD0V 00 Ta thường ký hiệu Σ := D ∈ Rm⋅n khai triển 00 A = U ΣV ∗ gọi khai triển kỳ dị ma trận A Các véc-tơ cột ma trận U gọi véc-tơ kỳ dị trái, véc-tơ cột ma trận V gọi véc-tơ kỳ dị phải, σi gọi giá trị kỳ dị ma trận A Để tìm khai triển kì dị ma trận A ta tìm véc-tơ riêng ∗ ∗ ∗ ma trận A A AA Cụ thể véc-tơ riêng đơn vị A A ∗ véc-tơ cột V , véc-tơ riêng đơn vị AA véc-tơ cột U , giá trị kỳ dị A bậc hai giá trị riêng Re λ≥0 Khi nhiễu có cấu trúc chấp nhận được, tức phương trình bị nhiễu (3.24) tính chất lạ với khối cỡ d, a h 50 Chứng minh Ta chứng minh E1 −1 B1 −1 sup Re λ≥0 G(λ) ≤ C A2 Ta viết lại G sau − G(λ) = CH(λ) −λB1 B2 e λE1 − A1 − e =C  −A2 − e =: CF (λ) Do ta có B22 e e − λτ D − −e λτ D3  −A2 − e   F ( λ) =  33 31 B22 e − λτ B − λτ B −λτ B /λ 31 e −λτ B B2  −D3   −  F (λ) = .0 32 33 B −B / λ −e −E + A / λ + e   −λτ D /λ Vì 0 32 e  − λτ B − λτ B − λτ B −λB11 B21 e  Nếu λ = điều tương đương với −A − e −λτ D   0 − λτ D − λE1 − A1 − e λτ D1 (3.37) − − λτ D − −λB B e λτ B  11 21 31 − λτ D −e B2 0 B3 D3 − λτ B −1  011 + +e − λτ E A1 D1/λ /λ B E 3 32  lim Re λ→+∞ −A − e −λτ D −D3    = − A1 2D −e − λτ B 31 / B2 e −λτ B B3 32 li λ→ +∞ 0 m ∞   = 11 B 22 0 nên limRe λ→+ F (λ) tồn λ B 1 −B / λ Re B E1 −1 B1 0 B3 B2 0 B3 Re λ→+∞ li m F ( λ) = − A2 D3 Do ta có E1 l i + G(λ) = C Re l i + F (λ) = −C A2 λ→ ∞ λ→ ∞ m m Re 51 D3 −1 B1 0 B2 , 0B 3 (3.37) thỏa mãn Rõ ràng sup G(λ) −1 lim ≤ Re λ≥0 −1 G(λ) Re λ→+∞ Từ mệnh đề 3.11 ta có ∆< sup G(λ) −1 Re λ≥0 Khi đó, hệ bị nhiễu (3.24) tính chất lạ với khối cỡ d, a, h (3.1) (đpcm) Chúng ta kết hợp kết để mô tả bán kính ổn định PTVP ĐS có chậm tính chất lạ tác động nhiễu có cấu trúc phù hợp Định lý 3.13 Giả sử (3.1) ổn định mũ, tính chất lạ chịu tác động nhiễu có cấu trúc (3.24) với ma trận cấu trúc B1, B2, B3 thỏa mãn (3.36) Khi rC ( E , A, D ; B , C ) = sup G(λ ) −1 Re λ≥0 Hơn nữa, ∆ < rC(E, A, D; B, C) (3.24) tính chất lạ với khối cỡ d, a, h với (3.1) Chứng minh Từ (3.28), ta có sup rC ( E , A, D ; B , C ) ≤ G(λ ) −1 Re λ≥0 Để chứng minh bất đẳng thức ngược lại, gọi ∆ nhiễu tùy ý làm cho (3.1) tính ổn định mũ Giả sử ∆< sup G(λ ) −1 Re λ≥0 Vì (3.1) tính chất lạ ổn định mũ, ta có α(H) < 0; từ Mệnh đề 3.9, ta có α(H∆) < Khi đó, từ Mệnh đề 3.12 phương trình 52 (3.38) bị nhiễu (3.24) tính chất lạ, đó, theo Định lý 3.4, ta thu phương trình bị nhiễu (3.24) ổn định mũ, điều mâu thuẫn Vậy sup ∆≥ G(λ ) Re λ≥0 −1 −1 rC(E, A, D; B, C) ≥ supRe λ≥0 G(λ ) , tức ta có (3.38) Cuối cùng, từ Mệnh đề 3.12, ta có (3.24) tính chất lạ ∆ < rC(E, A, D; B, C) (đpcm) Nhận xét Do nguyên lý cực đại [5] nên cận G(λ) nửa phải mặt phẳng phức đạt điểm xác định trục ảo vô Đối với PTVP ĐS có chậm tính chất lạ supRe ≥ G(λ) = sup λ=0 G(λ) Re λ , tức rC ( E , A, D ; B , C ) = sup Re λ=0 G(λ) −1 Như hệ quả, ta có kết tương tự cho trường hợp đặc biệt hệ tính chất lạ mà cặp (E, A) quy ind(E, A) ≤ Hệ 3.14 Xét hệ (3.1) với (E, A) quy thỏa mãn ind(E, A) ≤ 1, giả sử hệ ổn định mũ có dạng (1.1) Nếu hệ chịu tác động nhiễu có cấu trúc (3.24), ma trận cấu trúc B1 thỏa mãn − W 1B1 = B011 , với B11 ∈ Cd⋅p1 Khi bán kính ổn định có cấu trúc cho công thức rC ( E , A, D ; B , C ) = sup Re λ=0 G(λ) −1 Với PTVP ĐS chậm, ta biết [4] mục 2.1 rằng, nhiễu bảo toàn cấu trúc lũy linh dạng Weierstrass-Kronecker, ta mô tả bán kính ổn định có cấu trúc trường hợp (E, A) quy ind(E, A) > Ở mục 3.2, ta thấy tính ổn định mũ mô tả phổ H ta giả sử N D21 = N D22 = Sau đây, ta giả sử tính chất ổn định mũ 53 bảo toàn phương trình bị nhiễu (3.24), ma trận cấu trúc B1, B2, B3 thỏa mãn − − − W 1B1 = B011 , W 1B2 = B021 , W 1B3 = B31 , N B32 = 0, (3.39) B3 Bj1 ∈ Cd⋅pj, j = 1, 2, 3, B32 ∈ C(n−d)⋅p3 W ∈ Cn⋅n, N ∈ C(n−d)⋅(n−d) dạng (1.1) Xét hệ bị nhiễu không làm thay đổi cấu trúc lũy linh dạng Weierstrass-Kronecker (E, A), tức ma trận lũy linh N không gian bất biến trái tương ứng liên kết với giá trị đặc trưng ∞ bảo toàn, xem [2] cho trường hợp ind(E, A) = D = Tương tự phương pháp tiếp cận [2], ta đưa khoảng cách tới cặp gần có cấu trúc lũy linh thay đổi dn (E, A, D; B, C) = inf{ ∆ : (3.24) không bảo toàn cấu trúc lũy linh} C Với giả thiết (3.39), ta có kết sau cho trường hợp PTVP ĐS chậm Mệnh đề 3.15 Xét (3.1) với (E, A) quy ind(E, A) > 1, chịu tác động nhiễu biến đổi thỏa mãn (3.39) Khi đó, khoảng cách đến hệ gần có cấu trúc lũy linh thay đổi cho công thức − , d n ( E , A , D ; B , C ) = C 1 B1 C C = [C11 C12] với C11 ∈ Cq⋅r, C12 ∈ Cq⋅(n−r) Chứng minh Từ (3.39), cấu trúc lũy linh phương trình bị nhiễu (3.24) bảo toàn ma trận bị nhiễu Ir + B11∆1C11 không suy biến Do đó, sử dụng lại khoảng cách ma trận không suy biến đến ma trận suy biến gần nhất, xem [13], ta thu d n ( E , A , D ; B , C ) = C 1 B1 −1 C (đpcm) Định lý 3.16 Xét (3.1) ổn định mũ với (E, A) quy, ind(E, A) > giả sử (3.1) chịu tác động nhiễu biến đổi thỏa mãn (3.39) 54 Khi đó, bán kính ổn định cho công thức rC ( E , A, D ; B , C ) = sup Re λ=0 G(λ) −1 Hơn nữa, ∆ < rC(E, A, D; B, C) (E+C∆1B1, A+C∆2B2) quy cấu trúc lũy linh dạng Weierstrass-Kronecker (1.1) hệ bị nhiễu ổn định mũ Chứng minh Với giả thiết (3.39), tính toán đơn giản, ta G1(λ) = C11B11 , Re l i + G (λ) = l i + G (λ) = Re λ → ∞ λ→ ∞ lim m m Re λ→+∞ Do đó, limRe λ→+ G(λ) = C11B11 Từ thực tế supRe λ≥0 G(λ) ≥ ∞ limRe λ→+ G(λ ) Mệnh đề 3.15, phần lại chứng minh tương ∞ tự với chứng minh Định lý 3.13 Sử dụng nguyên lý cực đại, ta nhận supRe λ≥0 G(λ) = supRe λ=0 G(λ) (đpcm) Để minh họa kết phần này, ta xét ví dụ sau Ví dụ 3.4 Xét PTVP ĐS tuyến tính tính chất lạ 100 −1 020 0 x (t) = x(t) + 1 x(t − 1), 000 00 001 với (E, A) suy biến, chịu tác động nhiễu có cấu trúc 100 E= 0 ❀E= 000 −1 A= ❀A= 00 020 + δ1 δ1 δ1 00, 00 −1 + δ δ 2 δ δ2 + δ2 δ2 , 0 δ3 + δ3 2 δ3 D = 1 ❀ D = δ3 1 + δ3 + δ3 001 δ3 δ3 + δ3 Và giới thiệu lại dạng (3.24) với δ1 δ1 δ1 B1 = , B2 = , B3 = , C = I , ∆ = δ δ 2 δ 0 δ3 δ3 δ3 (3.40) 55 −λ 1 + λ −2(2 + e )  Ta D có = , H (λ  )= λ E − A − e − λ y −λ − −2 − e −e λ dễ kiểm tra α(H) < (3.40) ổn định mũ d − e Bằng tính toán đơn ự với n chuẩn − g λ giản, ta có  1+ λ l − +  1/2 Hơn nữa, người ta 0 kiểm tra ∆ với λ 1+ ∆ này, = hệ bị G(λ) = 2−ee−−λλ  nhiễu 2+e1−λ không  − − với chuẩn max C3⋅3 nên sup G(iπ) ∞ ∈i G(λ) λ R = G(λ) / 0, có = Do đó, Định lý 3.13 ta có ∞ = Ta ý rằng, việc sử dụng (3.29), nhiễu gây bất ổn định xâ = sup chất lạ α(H∆) = ∞ rC(E, A, D; B, C) λ∈iR có tính nghĩa hệ bị n h i ễ , u k h ô / n g ổn định tiệm cận Kết luận Luận văn hoàn thành Khoa Toán - Cơ - Tin học trường Đại học KHTN, ĐHQGHN sở báo [5] tính ổn định ổn định vững PTVP ĐS có chậm Kết luận văn xây dựng lại điều kiện cần đủ cho tính ổn định vững PTVP ĐS có chậm tác động nhiễu có cấu trúc, đồng thời công thức tính toán bán kính ổn định phức, vài trường hợp đặc biệt xét Chương PTVP ĐS PTVP thường có chậm Các kết thu luận văn khiêm tốn, phạm vi nghiên cứu chưa rộng (chỉ hệ tuyến tính hệ số hằng) Tuy nhiên, số vấn đề đặt nghiên cứu phổ điều kiện tính ổn định mũ PTVP ĐS trường hợp tổng quát công thức bán kính ổn định thực toán mở 57 Tài liệu tham khảo [1] Ascher U.M., Petzold L.R (1998): Computer methods for Ordinary Differential equations and differential-algebraic equations, SIAM [2] Byers R., Nichols N.K (1993): "On the stability radius of a generalized state-space system", Linear Algebra and its Applications, Vol 188-189, 113134 [3] Nguyễn Hữu Dư, Vũ Hoàng Linh (2005): "Implicit-system approach to the robust stability for a class of singularly perturbed linear systems", Systems and Control Letters, Vol 54, 33-41 [4] Nguyễn Hữu Dư, Vũ Hoàng Linh, Mehrmann V (2013), "Robust Stability of Differential-Algebraic Equations", Surveys in DifferentialAlgebraic Equations I, 63-95, Springer [5] Nguyễn Hữu Dư, Vũ Hoàng Linh, Mehrmann V., Đỗ Đức Thuận (2013): "Stability and robust stability of linear time-invariant delay differential-algebraic equations", SIAM Journal on Matrix Analysis and Applications, Vol 34, 1-32 [6] Golub G.H., Van Loan C.F (1996): Matrix Computations, 3rd edn The Johns Hopkins University Press, Baltimore P [7] Ha P., Mehrmann (2012): "Analysis and reformulation of linear delay differential-algebraic equations", Electronic Journal of Linear Algebra, Vol 23, 703-730 [8] Hale J.K (1977): Theory of functional-differential equations, SpringerVelag, Newyork-Heidelberg-Berlin 58 [9] Kunkel and V Mehrmann (2006): Differential-algebraic equations Analysis and numerical solution, EMS Publishing House, Z¨rich, u Switzerland [10] Vũ Hoàng Linh, Đỗ Đức Thuận (2015), "Spectrum-Based Robust Stability Analysis of Linear Delay Differential-Algebraic Equations", Numerical Algebra, Matrix Theory, Differential-Algebraic Equations and Control Theory, Springer International Publishing Switzerland, 533557 [11] Michiels W (2011): Spectrum-based stability analysis and stabilisation of systems described by delay differential algebraic equations, IET Control theory Appl [12] Qiu L., Bernhardsson B., Rantzer A., Davison E.J., Young P.M., Doyle J.C (1995): "A formula for computation of the real stability radius", Automatica, Vol 31, 879-890 [13] Nguyễn Khoa Sơn, Đỗ Đức Thuận (2011): "On the radius of surjectivity for rectan-gular matrices and its application to measuring stabilizability of linear systems under structured perturbations", Journal of Nonlinear and Convex Analysis, Vol 12 , 441-453 59 ... phương trình vi phân đại số 16 2.2 Bán kính ổn định phương trình vi phân thường có chậm 27 2 Bán kính ổn định PTVP thường có chậm 28 2 Hệ dương có chậm 30 Bán kính ổn định. .. chậm ⋄ Chương Bán kính ổn định hệ phương trình vi phân đại số có chậm Chương nội dung luận văn Trong đó, phân tích chứng minh kết bán kính ổn định phức PTVP ĐS có chậm tuyến tính hệ số Và kết đưa... Chương Một số kết bán kính ổn định Nội dung chương giới thiệu số kết công thức bán kính ổn định PTVP ĐS PTVP thường có chậm tuyến tính hệ số trường hợp đặc biệt phương trình vi phân đại số có chậm