CHƯƠNG II Số hóa bởi Trung tâm Học liệu – Đại học Thái Nguyên http //www lrc tnu edu vn 1 MỤC LỤC Trang Mở đầu 2 Chƣơng I Một số khái niệm về hệ phƣơng trình vi phân đại số 5 1 1 Phép chiếu Chỉ số của[.]
MỤC LỤC Trang Mở đầu Chƣơng I Một số khái niệm hệ phƣơng trình vi phân đại số 1.1 Phép chiếu - Chỉ số cặp ma trận 1.2 Hệ phương trình vi phân đại số tuyến tính với hệ số 1.3 Phân rã hệ phương trình vi phân đại số thành hệ phương trình vi phân thường hệ phương trình đại số 10 1.4 Sự ổn định (Lyapunov) hệ phương trình vi phân đại số 13 Chƣơng II Bán kinh ổn định hệ phƣơng trình vi phân đại số tuyến tính với ma trận hệ số 15 2.1 Bán kính ổn định phức hệ phương trình vi phân đại số 15 2.2 Liên hệ bán kính ổn định thực bán kính ổn định phức hệ phương trình vi phân đại số 24 Chƣơng III Bán kính ổn định hệ phƣơng trình vi phân đại số tuyến tính với nhiễu động 34 3.1 Hệ phương trình vi phân đại số tuyến tính với hệ số biến thiên 35 3.2 Nghiệm yếu khái niệm ổn định 37 3.3 Công thức bán kính ổn định 44 3.4 Các trường hợp đặc biệt 55 Kết luận 59 Tài liệu tham khảo 60 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên Tai ngay!!! Ban co the xoa dong chu nay!!! http://www.lrc-tnu.edu.vn MỞ ĐẦU Từ cuối kỷ XIX nhiều nhà khoa học quan tâm tìm lời giải cho tốn ổn định chuyển động Ở thời điểm đó, người ta đưa nhiều định nghĩa khác khái niệm này, chẳng hạn định nghĩa A.Poincaré, V.Rumyantsev, Chỉ từ A.M Lyapunov (1857-1918) công bố công trình “Bài tốn tổng qt tính ổn định chuyển động” vào năm 1892 Nga dịch sang tiếng Pháp (Problème général de la stabilité du mouvement) năm 1907, lý thuyết ổn định nghiên cứu cách có hệ thống trở thành phận quan trọng lý thuyết định tính phương trình vi phân Kể từ đó, lý thuyết ổn định nhiều nhà khoa học khắp giới quan tâm nghiên cứu Đến nay, kỷ trôi qua, lý thuyết ổn định lĩnh vực tốn học nghiên cứu sơi thu nhiều thành tựu rực rỡ, sâu sắc, như: vật lý, khoa học kỹ thuật công nghệ, sinh thái học, Lyapunov giải toán ổn định hai phương pháp, phương pháp số mũ đặc trưng Lyapunov (còn gọi phương pháp phổ hay phương pháp thứ Lyapunov) phương pháp hàm Lyapunov (còn gọi phương pháp thứ hai Lyapunov) Vào năm 70 kỷ trước, số tốn có liên quan đến phương trình vi phân dạng: A t x '(t ) +B t x(t ) đó, A det A t ,B C I , L Rn , x : I Rn , I a, , a số, t I Đây dạng đặc biệt phương trình vi phân đại số (differential algebraic equation-DAE) Ngay sau đó, loại phương trình vi phân nhiều nhà tốn học sâu nghiên cứu Để nghiên cứu DAE người ta thường làm sau: phân rã chúng nhờ phép chiếu để hệ phương trình vi phân thường hệ phương trình đại số Ngồi ra, Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn vài phương pháp khác Đến người ta tìm nhiều kết cho phương trình vi phân đại số tương tự phương trình vi phân thường chẳng hạn lý thuyết Floquet, tính ổn định tiệm cận nghiệm phương trình với ma trận hệ số Trong hai thập kỷ qua, từ khái niệm bán kính ổn định mà D.Hinrichsen A.J.Pritchard đưa ra, hai ơng hình thành hướng nghiên cứu nghiên cứu tính ổn định vững hệ động lực dựa khái niệm bán kính ổn định Hướng nghiên cứu thu hút ý tâm huyết nhiều nhà tốn học tính hiệu tính thời ứng dụng tốn kỹ thuật Nhóm tác giả Nguyễn Hữu Dư, Vũ Hoàng Linh nghiên cứu ổn định hệ phương trình vi phân đại số với ma trận hệ số phụ thuộc tham số thời gian đưa cơng thức bán kính ổn định báo “Stability radii for linear time - varying differential - algebraic equations with respect to dynamic perturbations” đăng tải JOURNAL OF DIFFERENTIAL EQUATIONS, June 2006 Đây báo sở để thực luận văn Luận văn gồm 61 trang, phần mở đầu, kết luận tài liệu tham khảo, gồm có ba chương: Chương I: Một số khái niệm hệ phương trình vi phân đại số Chương trình bày kiến thức sở để sử dụng chương sau Chương II: Bán kính ổn định hệ phương trình vi phân đại số tuyến tính với ma trận hệ số Chương trình bày tốn tính bán kính ổn định cho hệ phương trình vi phân đại số tuyến tính dạng Ax '(t ) - Bx(t ) A, B ma trận thực, det A Chương III: Bán kính ổn định hệ phương trình vi phân đại số tuyến tính với nhiễu động Chương nghiên cứu hệ phương trình vi phân đại số tuyến tính biến đổi theo thời gian có dạng: Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn A t x' t A Lloc 0, ; K n n B t x t ,t Lloc 0, ; K n , B n , cơng thức bán kính ổn định đưa Luận văn hoàn thành khoa Toán, trường Đại học Sư phạm - Đại học Thái Nguyên hướng dẫn ân cần, tỉ mỉ khoa học Cô giáo - Tiến sĩ Đào Thị Liên Qua tơi xin bày tỏ lịng biết ơn sâu sắc công lao vô bờ cô không quản thời gian cơng sức hướng dẫn tơi hồn thành luận văn Tôi xin cảm ơn thầy giáo khoa Tốn, khoa Sau Đại học, trường Đại học Sư phạm - Đại học Thái Nguyên đào tạo tạo điều kiện tốt để luận văn hồn thành Sau tơi xin bày tỏ tình cảm tha thiết dành cho gia đình tơi, quan nơi công tác (Trường PT Vùng Cao - Việt Bắc) động viên, tạo điều kiện cho yên tâm học tập, nghiên cứu Mặc dù cố gắng, song luận văn khó tránh khỏi hạn chế thiếu sót Chúng tơi mong nhận ý kiến đóng góp để luận văn hoàn thiện Thái Nguyên, tháng năm 2008 Học viên cao học Lƣu Thị Thu Hoài Số hóa Trung tâm Học liệu – Đại học Thái Nguyên http://www.lrc-tnu.edu.vn CHƢƠNG I MỘT SỐ KHÁI NIỆM VỀ HỆ PHƢƠNG TRÌNH VI PHÂN ĐẠI SỐ 1.1 Phép chiếu - Chỉ số cặp ma trận Định nghĩa 1.1.1 Cho P L P gọi phép chiếu P2 P Nhận xét 1.1.2 i) Cho P phép chiếu Khi đó, ta có: KerP ii) Mỗi phân tích n U Im P n V tồn phép chiếu P cho imP = U KerP = V, P gọi phép chiếu lên U dọc theo V Đặt Q:=I – P Q phép chiếu phép chiếu lên V dọc theo U Định nghĩa 1.1.3 (Chỉ số ma trận) Cho A L n Số tự nhiên k gọi số ma trận A, ký hiệu indA, số nhỏ mà KerAk indA k KerAk : KerAk KerAk Định lý 1.1.4 Với A L n ta ln có: imAk KerAk n với k thoả mãn 0 toán tử nhân L Lp 0, T ; K q , Lp 0, T ; K m I 0 L cho tốn tử FQG 1E khơng ổn định Chứng minh Để chứng minh mệnh đề trên, ta chọn T > cho 1 1 esssup FQG E esssup FQG E t T /3 0 L /3 (3.3.9) t sau ta xây dựng dãy đơn điệu ngặt Tn n 0, T cho 1 esssup FQG E t Tn ,Tn 0 L /3 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 53 http://www.lrc-tnu.edu.vn với n = 0, 1, … Theo định nghĩa cận đúng, tồn tập 0, T dương đo cho X 0 L FQG 1E t Cho 1 /3 X , t độ đo X cho a inf t,t X , b sup t,t X Rõ X ràng a b T Đặt T0 a Với n = 0, 1, … chọn Tn đo Tn , Tn X / 2n Dễ thấy dãy Tn n X T cho độ thoả mãn điều kiện Và ta dễ thấy tồn dãy f n suppfn FQG 1Ef n fn Tn ,Tn , 0 L 1 Lp 0, T ; K m , với Thay đổi điều kiện bổ đề 3.3.7 yếu chút, ln tồn tốn tử nhân L Lp 0, T ; K q , Lp 0, T ; K m , n = 0, 1, … cho n FQG f n n n 0 L supp Δ n h fn 1, , Tn , Tn , h Lp 0, T ; K q , h L p 0, T ; K q với supp h Ta đặt f : fn , h: n Tn ,Tn n n h h n với h L p 0, T ; K q g : f0 Dễ thấy L Lp 0, T ; K q , Lp 0, T ; K m nhân quả, 0 L FQG 1E f I Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 54 f0 g http://www.lrc-tnu.edu.vn Xuất phát từ f L p 0, T ; K m , g Lp 0, T ; K m mà toán tử I FQG 1E có nghịch đảo khơng bị chặn L Lp 0, T ; K m , Lp 0, T ; K m Từ mệnh đề 3.3.3 – 3.3.11 ta xây dựng cơng thức cho bán kính ổn định Định lý 3.3.12 Cho giả thiết A1, A2 Khi d K A, B; E , F sup L t0 0 ,L t0 Hệ 3.3.13 Cho A, B, E, F thực giả thiết A1, A2 Khi dC A, B; E, F dR A, B; E, F Chú ý 3.3.14 Ta ý tính đơn điệu L t coi hàm t (xem bổ đề 3.2.8 (b)), ta có sup L t0 lim L t0 t0 t0 0 So sánh với (3.4), ta thấy số hạng đặc biệt L độ đo trì tính vững số Đó khác biệt DAEs ODEs 3.4 Các trƣờng hợp đặc biệt 3.4.1 Các hệ nửa hiển Hệ (3.1) gọi nửa hiển, nghĩa A I n1 0 B t , B11 t B21 t B12 t B22 t (3.4.1) I n1 ma trận đơn vị, Bij i,j ma trận có số chiều tương ứng Giả thiết số có nghĩa B22 t khả nghịch hầu khắp nơi 0, Trong nhiều ứng dụng, hệ DAEs xuất dạng nửa hiển Cho tập hợp Q diag 0, I n2 dễ dàng có Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 55 http://www.lrc-tnu.edu.vn G I n1 B12 B22 , 0 B221B21 I Qs Ngồi ra, ta có t, s t, s B221B21 t, s t , s toán tử sinh khai triển phương trình vi phân thường y' B12 B221B21 y, B11 (3.4.2) coi ổn định mũ Giả thiết A2 tương đương với giả thiết tính bị chặn thực B221 , B221B21 B12 B221 Ta xây dựng hai ma trận E,F viết lại dạng suy biến sau E1 , E2 E F F1 F2 ma trận có số chiều tương ứng Bằng vài tính tốn với ma trận ta có L t0 u t t F1 22 F2 B B21 t, B12 B221E2 E1 u F2 B221E2u t , d t0 t u) t (L F2 B221E2u t , t (3.4.3) t0 Theo định lý 3.3.12, ta có d K A, B; E , F lim L t0 t0 1 22 , esssup F2 B E2 t t 3.4.2 Hệ quy ẩn hệ đại số Trước hết, ta xem xét hệ (3.1) với số hạng đầu A không suy biến hầu khắp nơi Và ta giả sử A bị chặn hầu khắp nơi 0, Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 56 Nhân hai http://www.lrc-tnu.edu.vn vế hệ (3.1) với A ta hệ hiển hoàn toàn nghiên cứu [12] Áp dụng định lý 3.3.12 với việc chọn P I, Q tầm thường nhất, công thức bán kính ổn định có giống [12] Trường hợp A = 0, hệ (3.1) thật trở thành hệ đại số B t x t , Ta đặt Q Điều tương đương với B I, P nơi 0, t khả nghịch hầu khắp nghịch đảo thực bị chặn Dễ thấy hệ có nghiệm tầm thường Hệ đại số ổn định tất hàm q Lp 0, ; K n Hệ không B t x t x Lp 0, ; K n q t có nghiệm hệ phụ thuộc vào tính liên tục vế phải (Theo chuẩn Lp ) Theo ý nghĩa kết mục 3.3, toán tử nhiễu với 0 L 1 esssup FB E t (3.4.4) , t hệ nhiễu lại ổn định (3.4.4) cho cận tốt nhất, nghĩa với tồn nhân với 0 L làm tính ổn định hệ đại số 3.4.3 Hệ bất biến theo thời gian Ta giả sử tất ma trận A, B, E,F bất biến theo thời gian Dễ thấy giả thiết A2 không cần thiết Bằng phép biến đổi Fourier – Plancherel [11], phát biểu sau, mà thực chất mở rộng định lý 2.1 [11] hệ phương trình vi phân đại số có số 1, chứng minh Mệnh đề 3.4.3.1 Cho A, B, E,F bất biến theo thời gian, hệ (3.1) có số ổn định mũ Nếu chọn p = (xét tính ổn định L 2) Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 57 http://www.lrc-tnu.edu.vn L t0 Hàm F L0 sup F A B E iR A B E gọi hàm truyền nhân tạo liên kết tới (3.1) Ta ý ổn định mũ hệ bất biến theo thời gian (3.1) có nghĩa xác tất giá trị riêng cặp ma trận (A, B) có phần thực âm Như hàm truyền xác định trục ảo iR mặt phẳng phức Kết là, định lý 3.3.12 viết lại sau Định lý 3.4.3.2 Cho A, B, E, F bất biến theo thời gian, hệ (3.1) có số ổn định mũ Nếu chọn p = d C A, B; E , F L0 sup F A B E iR Chứng minh 0 Ta cần L L Theo bổ đề 3.2.8 (c), ta dễ dàng xác định giới hạn lim F A B E FQG 1E cách biến đổi cặp ma trận hệ số (A, B) theo chuẩn Weierstrass – Kronecker (xem [5], [9]) theo hệ nửa hiển (3.4.1) làm rõ ta có: sup F A B E FQG 1E 0 L iR Từ định lý 3.4.3.2, việc tính bán kính ổn định hệ bất biến theo thời gian giải cách tối ưu giải số, ví dụ xem [4], [18] Cuối lưu ý đẳng thức định lý 3.4.3.2 mang nghĩa xác bán kính ổn định phức nhiễu động nghiên cứu trùng với bán kính ổn định phức nhiễu tĩnh (nghĩa là ma trận bất biến theo thời gian) xét [6], [7], [18] Vì định lý 3.4.3.2 tổng quát kết trước hệ tuyến tính ODEs [10] Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 58 http://www.lrc-tnu.edu.vn KẾT LUẬN Bản luận văn trình bày sở luận án [2] báo [8] bán kính ổn định hệ phương trình vi phân đại số tuyến tính với nhiễu động thu số kết sau: Đưa công thức cho tốn tìm bán kính ổn định phức hệ phương trình vi phân đại số tuyến tính hệ số hằng, bán kính ổn định phức bán kính ổn định thực hệ phương trình vi phân đại số Một định nghĩa xây dựng bán kính ổn định cho phương trình vi phân đại số tuyến tính có hệ số phụ thuộc tham số thời gian đưa Ở mục cuối cùng, vài trường hợp đặc biệt xem xét, phân tích Những kết thu khiêm tốn, điều kiện đưa tương đối ngặt song kết bước đầu mà thu mang ý nghĩa định Những kết mở số vấn đề để tiếp tục nghiên cứu: Bán kính ổn định thực hệ phương trình vi phân đại số Nếu điều kiện hệ dương đảm bảo d d R khơng Nới lỏng giả thiết A1, A2 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 59 http://www.lrc-tnu.edu.vn TÀI LIỆU THAM KHẢO Tiếng Việt [1] Lê Thị Lan (2000), Một số tính chất hệ phương trình vi phân tuyến tính đại số 1, số 2, Luận văn thạc sĩ Toán học, Đại học Sư phạm Hà Nội, Hà Nội [2] Đào Thị Liên (2004), Về ổn định hệ phương trình vi phân hệ phương trình vi phân đại số, Luận án Tiến sĩ Toán học, Đại học Sư phạm Hà Nội, Hà Nội [3] Phạm Văn Việt (2005), Về ổn định nghiệm số phương trình vi phân đại số, Luận án Tiến sĩ Toán học, Đại học Sư phạm Hà Nội, Hà Nội Tiếng Anh [4] Bracke M (2000), On stability radii of parametrized linear diffirential – algebraic systems, PhD thesis, University of Kaiserslautern [5] Brenan K.E., Campbell S.L., Petzold L.R (1996), Numerical Solution of Initial Value Problems in Diffirential - Algebraic Equations, SIAM, Philadelphia, PA [6] Nguyễn Hữu Dư (1999), “Stability radii for diffirential - algebraic equations”, Vietnam J Math (27), 379-382 [7] Nguyễn Hữu Dư, Vũ Hoàng Linh (2006), “Robust stability of implicit linear systems containing a small parameter in ther leading tern”, IMA J math Control Inform, (23), 67-84 [8] Nguyễn Hữu Dư, Vũ Hoàng Linh (2006) , “Stability radii for linear timevarying diffirential - algebraic equations with respect to dynamic perturbations”, Journal of Differential Equations, (230), 579-599 [9] Griepentrog E., März R (1986), Diffirential - Algebraic Equations and Their Numerical Treatment, Teubner - Texte Math., Teubner, Leipzig Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 60 http://www.lrc-tnu.edu.vn [10] Hinrichsen D., Pritchard A.J (1992), “Destabilization by output feedback, Differential Integral Equations” (2), 357-386 [11] Hinrichsen D., Pritchard A.J (1986), “Stability radius for structured perturbations and the algebraic Riccati equation”, Systems Control Lett (8), 105-113 [12] Jacob B (1998), “A formula for the stability radius of time-varying systems”, J Diffirential Equations (142), 167-187 [13] Kolmogorov A.N., Fomin S.V (1970), Introductory Real Analysis, revised English ed., Dover, New York [14] Lamour R., März R., Mattheij R.M.M (1992), On the stability behavious of systems obtained by index - reduction, Fachbereich Mathematic, O - 1086 Berlin, PF1297, BRD, Berlin [15] März R (1991), On quasilinear index diffirential algebraic equations, Fachbereich Mathematic, O - 1086 Berlin, PF1297, BRD, Berlin [16] März R (1992), Numerical methods for diffirential - algebraic equations, in: Acta Numer., Cambridge Univ Press, Cambridge, pp 141-198 [17] März R (1997), Extra-ordinary Diffirential Equation Attempts to an analysis of Algebraic System, Preprint N0 97-8, Humboldt Universitatzu Berlin [18] Qiu L., Davision E.J (1992), “The stability robustness of generalized eigenvalues”, IEEE trans Automat Control, (37), 886-891 Số hóa Trung tâm Học liệu – Đại học Thái Nguyên 61 http://www.lrc-tnu.edu.vn