1 LỜI CAM ĐOAN Tôi xin cam đoan đây là công trình nghiên cứu của riêng tôi, được hoàn thành dưới sự hướng dẫn của GS TSKH Vũ Ngọc Phát Các kết quả viết chung với tác giả khác đã được sự nhất trí của đ[.]
1 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi, hồn thành hướng dẫn GS TSKH Vũ Ngọc Phát Các kết viết chung với tác giả khác trí đồng tác giả đưa vào luận án Các kết phát biểu luận án hoàn toàn trung thực chưa cơng bố cơng trình khác Nghiên cứu sinh Lê Anh Tuấn LỜI CẢM ƠN Luận án thực hoàn thành hướng dẫn nghiêm khắc, tận tình, chu đáo GS TSKH Vũ Ngọc Phát Tôi xin bày tỏ kính trọng lịng biết ơn sâu sắc đến Thầy, người dìu dắt tơi bước vào nghiên cứu khoa học gần mười năm qua Ngoài dẫn mặt khoa học, động viên lòng tin tưởng mà Thầy dành cho nguồn động lực lớn thúc đẩy tơi tiến trình nghiên cứu Tôi xin trân trọng gửi lời cảm ơn đến Ban Giám hiệu, Phòng Sau Đại học, Ban Chủ nhiệm Khoa Toán - Tin, Trường Đại học Sư phạm Hà Nội, đặc biệt thầy cô giáo Bộ mơn Giải tích, Khoa Tốn - Tin, Trường Đại học Sư phạm Hà Nội giúp đỡ, động viên, tạo môi trường học tập nghiên cứu thuận lợi cho suốt quãng thời gian làm nghiên cứu sinh Tơi xin bày tỏ lịng biết ơn đến Ban Giám hiệu trường Đại học Khoa học Huế, thầy cô anh chị em đồng nghiệp công tác Khoa Toán, Trường Đại học Khoa học Huế tạo điều kiện thuận lợi để hỗ trợ suốt q trình học tập nghiên cứu Tơi xin dành lời cảm ơn sau cho đại gia đình tơi, người ln u thương, chia sẻ, động viên tơi vượt qua khó khăn, thử thách khoa học sống để hoàn thành luận án Mục lục Lời cam đoan Lời cảm ơn Mục lục Danh mục ký hiệu MỞ ĐẦU LÝ DO CHỌN ĐỀ TÀI TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU MỤC ĐÍCH, ĐỐI TƯỢNG VÀ PHẠM VI NGHIÊN CỨU 16 PHƯƠNG PHÁP NGHIÊN CỨU 16 KẾT QUẢ CỦA LUẬN ÁN 17 BỐ CỤC CỦA LUẬN ÁN 18 Chương MỘT SỐ KIẾN THỨC CHUẨN BỊ 19 1.1 BÀI TOÁN ỔN ĐỊNH VÀ ỔN ĐỊNH HĨA HỆ PHƯƠNG TRÌNH CĨ TRỄ 19 1.1.1 Bài toán ổn định 19 1.1.2 Bài toán ổn định hóa 26 1.2 BÀI TOÁN ĐIỀU KHIỂN H∞ 27 1.2.1 Không gian H∞ 1.2.2 Bài toán điều khiển H∞ 27 29 1.3 BẤT ĐẲNG THỨC MA TRẬN TUYẾN TÍNH 31 Chương BÀI TOÁN ĐIỀU KHIỂN H∞ CHO LỚP HỆ NƠ-RON CÓ TRỄ BIẾN THIÊN HỖN HỢP 34 2.1 PHÁT BIỂU BÀI TOÁN 34 2.2 KẾT QUẢ CHÍNH 37 2.3 VÍ DỤ MINH HỌA 48 Chương BÀI TOÁN ĐIỀU KHIỂN H∞ TRONG THỜI GIAN HỮU HẠN CHO LỚP HỆ RỜI RẠC TUYẾN TÍNH CĨ TRỄ BIẾN THIÊN THEO THỜI GIAN DẠNG KHOẢNG 51 3.1 KHÁI NIỆM ỔN ĐỊNH TRONG THỜI GIAN HỮU HẠN 51 3.2 PHÁT BIỂU BÀI TOÁN 53 3.3 CÁC KẾT QUẢ CHÍNH 55 3.4 CÁC VÍ DỤ MINH HỌA 65 Chương BÀI TOÁN ĐIỀU KHIỂN H∞ TRONG THỜI GIAN HỮU HẠN CHO LỚP HỆ NƠ-RON RỜI RẠC SUY BIẾN CÓ TRỄ BIẾN THIÊN THEO THỜI GIAN DẠNG KHOẢNG 70 4.1 SƠ LƯỢC VỀ HỆ RỜI RẠC SUY BIẾN TUYẾN TÍNH 70 4.2 PHÁT BIỂU BÀI TOÁN 74 4.3 CÁC KẾT QUẢ CHÍNH 77 4.4 CÁC VÍ DỤ MINH HỌA 97 KẾT LUẬN VÀ KIẾN NGHỊ 101 DANH MỤC CÁC CƠNG TRÌNH KHOA HỌC CỦA TÁC GIẢ ĐÃ CÔNG BỐ LIÊN QUAN ĐẾN LUẬN ÁN 103 TÀI LIỆU THAM KHẢO 104 DANH MỤC CÁC KÝ HIỆU R, R+ tập số thực tập số thực không âm tương ứng Z+ tập số nguyên không âm N tập số tự nhiên C tập số phức Re(s) phần thực số phức s Rn khơng gian Euclide thực n chiều Rn×r khơng gian ma trận thực có kích thước (n × r) n P tích vơ hướng hai véc tơ x, y Rn : xT y = xi yi i=1 P 1/2 n n chuẩn Euclide véc tơ x ∈ R , kxk = xi hx, yi = xT y kxk I i=1 ma trận vuông đơn vị với số chiều phù hợp ∗ phần tử đường chéo ma trận đối AT ma trận chuyển vị ma trận A A−1 ma trận nghịch đảo ma trận A A−T viết tắt (A−1 )T λ(A) tập hợp tất giá trị riêng ma trận A λmax (A) := max{Reλ : λ ∈ λ(A)} λmin (A) := min{Reλ : λ ∈ λ(A)} xứng σmax (A) giá trị suy biến (singular value) lớn ma trận A A>0 A ma trận nửa xác định dương, tức xT Ax > ∀x ∈ Rn A>0 A ma trận xác định dương, tức xT Ax > ∀x ∈ Rn \ {0} diag{A1 , , An } ma trận đường chéo với Ai phần tử thứ i đường chéo K C([a, b], Rn ) tập hàm liên tục không giảm u : R+ −→ R+ , u(0) = 0, u(s) > ∀s > không gian hàm liên tục [a, b], nhận giá trị Rn với chuẩn kxkC = max kx(t)k a6t6b C ([a, b], Rn ) không gian hàm khả vi liên tục [a, b], nhận giá trị Rn với chuẩn kxkC = max {kx(t)k, kx(t)k} ˙ a6t6b L2 ([0, ∞), Rn ) không gian hàm ω : [0, ∞) −→ Rn bình phương khả R∞ tích [0, ∞), nghĩa kω(t)k2 dt < ∞ LMI bất đẳng thức ma trận tuyến tính (viết tắt cụm từ tiếng Anh “linear matrix inequality”) FTS tính ổn định thời gian hữu hạn (viết tắt cụm từ tiếng Anh “finite-time stability”) LS tính ổn định Lyapunov (viết tắt cụm từ tiếng Anh “Lyapunov stability”) RFDE phương trình vi phân hàm có trễ (viết tắt cụm từ tiếng Anh “retarded functional differential equation”) MỞ ĐẦU LÝ DO CHỌN ĐỀ TÀI Lý thuyết ổn định nhánh quan trọng lý thuyết định tính hệ phương trình vi phân mà nhà toán học người Nga A.M Lyapunov khởi xướng từ năm cuối kỷ XIX Với bề dày lịch sử kỷ đến thời điểm lý thuyết ổn định Lyapunov cịn lĩnh vực nghiên cứu có sức lơi lớn toán học với ngày nhiều ứng dụng quan trọng tìm thấy học, vật lý, hóa học, cơng nghệ thơng tin, sinh thái, mơi trường, v.v trở thành nhánh nghiên cứu thiếu lý thuyết hệ thống ứng dụng [18, 20, 24, 29, 30] Cùng với tính ổn định nghiệm, người ta cịn quan tâm tới việc thiết kế điều khiển cho tác động vào hệ điều khiển, hệ trở nên ổn định Bài toán gọi tốn ổn định hóa hệ điều khiển người ta bắt đầu nghiên cứu tính ổn định hóa hệ điều khiển từ năm 1960 Mặt khác, mơ hình tốn học (được xây dựng từ toán kỹ thuật thực tiễn) thường xuất độ trễ thời gian Các đại lượng trễ hình thành cách tự nhiên, khơng thể tránh khỏi trình truyền tải, xử lý liệu người ta diện nhiều ảnh hưởng đến dáng điệu tính chất hệ, có tính ổn định, tính chất thiết yếu hệ kỹ thuật [18, 28, 43] Chính vậy, việc nghiên cứu tính ổn định điều khiển cho hệ có trễ tốn có ý nghĩa thực tế, nhiều học giả quan tâm năm gần [2, 8, 12, 14, 41, 57] Các hướng nghiên cứu quan trọng bao gồm việc đánh giá định tính phụ thuộc độ trễ tính ổn định xây dựng tiêu chuẩn mới, tân tiến để áp dụng cho nhiều mơ hình tổng quát phức tạp hơn, phù hợp với mơ hình kỹ thuật đại Bên cạnh đó, trình thực tiễn thường xảy cách khơng chắn (nghĩa là, có xuất đại lượng “nhiễu” hệ thống) Các nhiễu xuất sai số vận hành, ảnh hưởng lẫn thành tố hệ thống hệ thống khác Vì vậy, việc địi hỏi phải biết xác tất tham số hệ mơ hình điều khơng tưởng khó vận dụng thực tế Do đó, việc đánh giá tối ưu mức ảnh hưởng nhiễu đầu hệ thống (bài toán H∞ ) tốn có tính thời sự, nhiều nhà toán học kỹ sư quan tâm nghiên cứu Các cách tiếp cận khác phát triển số lượng lớn kết quan trọng điều khiển H∞ cho nhiều lớp hệ có trễ công bố thời gian qua [4, 8, 13, 44, 51, 53, 57, 59, 64] Tuy nhiều vấn đề mở thú vị quan trọng lý thuyết lẫn ứng dụng chưa giải quyết, đặc biệt kết có toán H∞ cho lớp hệ điều khiển có trễ tổng qt cịn khiêm tốn cần tiếp tục nghiên cứu sâu Đó động lực để thực đề tài TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU Trong cách tiếp cận theo miền thời gian (time-domain approach), phương pháp Lyapunov trực tiếp công cụ hữu hiệu để nghiên cứu toán ổn định điều khiển H∞ cho hệ có trễ như: hệ tuyến tính, hệ phi tuyến, hệ nơ-ron, hệ suy biến, v.v Qua đó, điều kiện giải toán điều khiển H∞ cho hệ ô-tô-nôm thiết lập dạng bất đẳng thức ma trận tuyến tính phương trình Riccati đại số; cịn với hệ khơng ơ-tơ-nơm điều kiện giải tốn thiết lập thơng qua phương trình Riccati vi phân Hệ nơ-ron có trễ vừa đề cập đến lớp hệ phương trình vi phân hàm đặc biệt, nghiên cứu cách rộng rãi hai thập kỷ qua ứng dụng thành công nhiều lĩnh vực như: nhớ kết hợp (associative memory), nhận dạng phân loại mẫu, xử lý tín hiệu, xử lý ảnh, giải tốn tối ưu, v.v Mặc dù có số cơng trình đề cập đến tốn điều khiển H∞ cho hệ nơ-ron có trễ [35, 40, 46, 47, 48] chủ đề lâu đạt trọn vẹn điều thúc đẩy quan tâm đáng kể luận án Vì lý đó, lớp hệ đề cập luận án toán điều khiển H∞ hệ nơ-ron có trễ biến thiên hỗn hợp (nghĩa yếu tố trễ gồm hai loại: trễ dạng rời rạc trễ dạng tích phân): x(t) ˙ = −Ax(t) + W0 f (x(t)) + W1 g(x(t − h(t))) + W2 Z t c(x(s))ds t−k(t) + Bu(t) + Cω(t) z(t) = Ex(t) + M x(t − h(t)) + N u(t), x(t) = ϕ(t), t ∈ [−d, 0], (1) t > 0, d = max{h2 , k}, x(t) = [x1 (t), x2 (t), , xn (t)]T ∈ Rn véc tơ trạng thái hệ nơ-ron; u(t) ∈ Rm véc tơ điều khiển; ω(t) ∈ Rr biến nhiễu/không chắn; z(t) ∈ Rs hàm quan sát đầu hệ nơ-ron; A = diag{a1 , a2 , , an } ma trận đường chéo dương; W0 , W1 , W2 , B, C, E, M, N ma trận thực cho trước có số chiều thích hợp; f (·), g(·), c(·) hàm kích hoạt hệ; h(t), k(t) hàm trễ hệ thỏa mãn điều kiện h1 h(t) h2 , k(t) k Năm 2009, toán ổn định mũ cho hệ nơ-ron x(t) ˙ = −(A+∆A(t))x(t)+(W0 +∆W0 (t))f (x(t))+(W1 +∆W1 (t))f (x(t−h(t))) với hàm trễ h(t) biến thiên liên tục dạng khoảng có đạo hàm bị chặn xét Kwon Park [32] Cịn tốn ổn định hóa dạng mũ tác giả Phat, Trinh [45] đề xuất vào năm 2010 cho hệ nơ-ron 10 với trễ hỗn hợp x(t) ˙ = −Ax(t) + W0 f (x(t)) + W1 g(x(t − h(t))) + W2 Z t c(x(s))ds + Bu(t), t−k(t) hàm trễ h(t), k(t) giả thiết thỏa mãn điều kiện: h(t) ˙ h, h(t) δ < 1, k(t) k ∀t > Khơng lâu sau đó, kết mở rộng sang trường hợp trễ rời rạc h(t) hàm liên tục, nhận giá trị khoảng hai tác giả Thuan, Phat [52] Năm 2012, Sakthivel cộng [47] xét tốn điều khiển H∞ cho hệ nơ-ron có trễ hỗn hợp (và khơng có trễ hàm quan sát) x(t) ˙ = −(A + ∆A)x(t) + (W0 + ∆W0 )f (x(t)) + (W1 + ∆W1 )g(x(t − h(t))) Z t c(x(s))ds + u(t) + (C + ∆C)ω(t), + (W2 + ∆W2 ) t−k(t) z(t) = Ex(t), ˙ với hàm trễ h(t), k(t) thỏa mãn: h(t) h, h(t) δ, k(t) k ∀t > Trong cơng trình này, tác giả thu tính ổn định hóa dạng tiệm cận điều kiện H∞ Sang năm 2013, tác giả Phat, Trinh [46] tiếp tục nghiên cứu toán điều khiển H∞ cho hệ nơ-ron có trễ x(t) ˙ = −Ax(t) + W0 f (x(t)) + W1 g(x(t − τ1 (t))) + Bu(t) + Cω(t), z(t) = Ex(t) + M h(x(t − τ2 (t))) + N u(t), với hai trường hợp xét: hàm trễ τ1 (t), τ2 (t) khả vi có đạo hàm bị chặn số thực dương bé hàm trễ bị chặn không thiết khả vi Từ đó, tác giả thu tính ổn định hóa dạng mũ điều kiện H∞ ứng với trường hợp Như vậy, kết nêu tính ổn định điều khiển H∞ phần lớn bị hạn chế giả thiết độ trễ hàm khả vi có đạo hàm bị chặn đơn giản hàm bị chặn Hiện việc nghiên cứu toán điều khiển H∞ cho lớp hệ phương trình (1) với độ trễ h(t) liên tục, khơng địi hỏi ... dù có số cơng trình đề cập đến tốn điều khiển H∞ cho hệ nơ-ron có trễ [35, 40, 46, 47, 48] chủ đề lâu đạt trọn vẹn điều thúc đẩy quan tâm đáng kể luận án Vì lý đó, lớp hệ đề cập luận án toán điều. .. biến hay hệ thơng thường • Phạm vi nghiên cứu ◦ Nội dung 1: Bài toán điều khiển H∞ cho lớp hệ nơ-ron có trễ biến thiên hỗn hợp ◦ Nội dung 2: Bài toán điều khiển H∞ thời gian hữu hạn cho lớp hệ rời... cứu luận án ? ?Bài toán điều khiển H∞ cho số lớp hệ phương trình có trễ? ?? Cụ thể hơn, yếu tố trễ quan tâm hàm biến thiên theo thời gian, có giá trị thuộc khoảng R N tùy trường hợp mà hệ xét hệ suy