Thông tin tài liệu
ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC SƯ PHẠM –––––––––––––––––––––––– ĐỖ THỊ HƯƠNG lu an n va p ie gh tn to PHƯƠNG TRÌNH VI PHÂN ĐẠI SỐ CĨ TRỄ TRONG LÝ THUYẾT ĐIỀU KHIỂN d oa nl w nf va an lu LUẬN VĂN THẠC SĨ TOÁN HỌC z at nh oi lm ul z an Lu va http://www.lrc.tnu.edu.vn n Số hóa Trung tâm Học liệu – ĐHTN m co l gm @ THÁI NGUYÊN - 2016 ac th si ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC SƯ PHẠM –––––––––––––––––––––––– ĐỖ THỊ HƯƠNG lu an n va PHƯƠNG TRÌNH VI PHÂN ĐẠI SỐ CĨ TRỄ TRONG LÝ THUYẾT ĐIỀU KHIỂN gh tn to p ie Chuyên ngành: Tốn giải tích d oa nl w Mã số: 60 46 01 02 nf va an lu LUẬN VĂN THẠC SĨ TOÁN HỌC z at nh oi lm ul Người hướng dẫn khoa học: TS Đào Thị Liên z m co l gm @ an Lu THÁI NGUYÊN - 2016 n va http://www.lrc.tnu.edu.vn ac th Số hóa Trung tâm Học liệu – ĐHTN si LỜI CAM ĐOAN Tơi xin cam đoan cơng trình nghiên cứu riêng tôi, kết nghiên cứu trung thực chưa công bố cơng trình khác Thái Ngun, tháng năm 2016 Tác giả luận văn lu an Đỗ Thị Hương n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu Số hóa Trung tâm Học liệu – ĐHTNihttp://www.lrc.tnu.edu.vn n va ac th si LỜI CẢM ƠN Luận văn hồn thành khoa Tốn, trường Đại học sư phạm - Đại học Thái Nguyên hướng dẫn tận tình khoa học giáo - Tiến sĩ Đào Thị Liên Qua em xin bày tỏ lời cám ơn chân thành lòng biết ơn sâu sắc công ơn vô bờ bến cô không quản thời gian cơng sức hướng dẫn em hồn thành luận văn Em xin cảm ơn Ban Giám hiệu, phòng ban chức trường Đại học sư phạm - Đại học Thái Nguyên đặc biệt thầy cô giáo khoa tốn Bộ mơn Giải tích tạo điều kiện tốt để em hoàn thành luận văn Sau cùng, em xin bày tỏ tình cảm tha thiết dành cho gia lu đình, bạn bè đồng nghiệp động viên, tạo điều kiện cho em yên an n va tâm học tập nghiên cứu tn to Mặc dù cố gắng, song luận văn khó tránh khỏi hạn chế gh thiếu sót Em mong nhận ý kiến đóng góp q thầy p ie bạn để luận văn hoàn thiện w Thái Nguyên, tháng năm 2016 d oa nl Học viên cao học nf va an lu Đỗ Thị Hương z at nh oi lm ul z m co l gm @ an Lu Số hóa Trung tâm Học liệu – ĐHTNiihttp://www.lrc.tnu.edu.vn n va ac th si MỤC LỤC Trang Lời cam đoan i Lời cảm ơn ii Mục lục iii Danh mục kí hiệu viết tắt iv MỞ ĐẦU Chương 1: KIẾN THỨC CƠ SỞ 1.1 Phương trình vi phân đại số (DAE) lu 1.1.1 Một số khái niệm DAE an 1.1.2 Chỉ số DAE va n 1.1.3 Phương pháp giải số DAE 10 gh tn to 1.2 Phương trình vi phân thường có trễ (DODE) 11 ie 1.2.1 Một số khái niệm kết DODE 11 p 1.2.2 Phương pháp số giải DODE 16 oa nl w 1.3 Phương pháp Radau IIA tìm nghiệm cho ODE cương 16 1.3.1 Bậc hội tụ 18 d an lu 1.3.2 Các gián đoạn nghiệm 18 nf va 1.3.3 Giải phương trình phi tuyến 20 lm ul 1.4 Các hệ điều khiển 22 1.4.1 Hệ điều khiển tuyến tính liên tục 22 z at nh oi 1.4.2 Các hệ điều khiển rời rạc 28 Chương 2: PHƯƠNG TRÌNH VI PHÂN ĐẠI SỐ CÓ TRỄ TRONG LÝ z THUYẾT ĐIỀU KHIỂN 35 gm @ 2.1 Bài toán sisofeed5 35 l 2.2 Tính chất nghiệm 38 co m 2.3 Các hàm đầu vào 42 an Lu Số hóa Trung tâm Học liệu –iiiĐHTN http://www.lrc.tnu.edu.vn n va ac th si 2.4 Lựa chọn công thức 45 2.5 Ước lượng điều khiển sai số 47 2.6 Các bước “dài” 50 2.7 Các gián đoạn 52 2.7.1 Sự lan truyền 52 2.7.2 Theo dõi gián đoạn 55 KẾT LUẬN 59 TÀI LIỆU THAM KHẢO 60 lu an n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu Số hóa Trung tâm Học liệu – ĐHTNivhttp://www.lrc.tnu.edu.vn n va ac th si DANH MỤC CÁC KÍ HIỆU VIẾT TẮT lu an Code: Mã DAE: Phương trình vi phân đại số DDAE: Phương trình vi phân đại số có trễ DODE: Phương trình vi phân thường có trễ FOH: Mẫu bậc GENLTI: Hệ tuyến tính bất biến tổng quát LTI: Hệ tuyến tính bất biến ODE: Phương trình vi phân thường ZOH: Mẫu bậc không n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu Số hóa Trung tâm Học liệu – ĐHTNvhttp://www.lrc.tnu.edu.vn n va ac th iv si MỞ ĐẦU Các hệ tuyến tính bất biến (LTI) mơ hình lý thuyết điều khiển Các mơ hình với trạng thái có trễ xuất cách tự nhiên xây dựng mối liên hệ ngược hàm truyền sơ cấp với trễ đầu vào đầu Các nghiên cứu cho thấy đối số trễ thường có ảnh hưởng lớn đến dáng điệu nghiệm hệ động lực nói chung Những kết tương tự cho hệ phương trình vi phân đại số (DDAE) cịn ít, chúng có nhiều ứng dụng tự nhiên, hóa học, sinh học,… Do đó, luận văn tập trung trình bày nghiệm dạng số DDAE hai tác giả L F Shampine, P Gahinet đưa lu an báo “Delay differential algebraic equations in control theory” có va n dạng tn to x(t ) Ax(t ) B1u (t ) B2 w(t ) (0.1) (0.2) p ie gh z (t ) C2 x(t ) D21u (t ) D22 w(t ) w ma trận A, B1 , B2 , D21 , D22 không đổi hàm đầu vào u (t ) oa nl trơn khúc với t Hàm w xác định N số hạng vectơ (cột) d z với N trễ cố định , , N an lu mơ hình hóa bắt lm ul Việc nf va w(t ) [ z1 (t 1 ), , z N (t N )]T đầu t 0 từ điểm dừng z at nh oi u(t ), x(t ), z(t ) 0, 0, với t Trong suốt trình tính tốn giá trị x(t ), z (t ) 0, T f , phải ước lượng hàm đầu z y (t ) C1x(t ) D11u (t ) D12 w(t ) @ l gm (0.3) co xác định ma trận C1 , D11 , D12 Thường thì, hàm đầu vào m u(t ) có bước nhảy t Theo định nghĩa, nghiệm x(t ) liên tục an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th si bước nhảy u (t ) kéo theo bước nhảy z (t ), x(t ), y (t ) Khi giá trị trễ, nguyên nhân làm cho bước nhảy lan truyền khắp 0, T f Chúng ta thấy (0.1), (0.2) hệ DDAE Chúng ta cố gắng nghiên cứu sâu để tìm cách giải số cho phương trình vi phân đại số (DAE) phương trình vi phân thường có trễ (DODE) Có kết nghiên cứu DDAE nghiên cứu có phần lớn trọng việc ứng dụng DDAE Vì mơ hình hệ tuyến tính bất biến tổng quát (GENLTI) thường sử dụng thiết kế nên việc giải chúng cách nhanh chóng điều quan trọng Để làm điều đó, cố gắng tận dụng triệt để đặc tính lu an chúng Chẳng hạn như, xử lý trước mơ hình để đảm bảo n va tất giá trị trễ dương Những kết hai tác giả L F Shampine, P tn to Gahinet dựa kết trước tác U Ascher L ie gh Petzold trình bày nghiệm dạng số DDAE DDAE trung hòa (xem [5]); p phương pháp sử dụng máy tính để giải số phương trình vi phân thường (ODE) nl w DAE (xem [6]); W.H Enright H Hayashi trình bày cách giải DDAE d oa dựa phương pháp Runge−Kutta liên tục với điều khiển khuyết (xem [8]) an lu nhiều tác giả khác nf va Luận văn gồm 59 trang, phần mở đầu, kết luận tài liệu tham khảo, lm ul nội dung gồm có hai chương Chương Kiến thức sở z at nh oi Chương trình bày số khái niệm DAE, DODE hệ điều khiển tuyến tính liên tục rời rạc lý thuyết điều khiển Đây kiến z thức sở cần thiết để sử dụng chương @ gm Chương Phương trình vi phân đại số có trễ lý thuyết điều khiển co l Nội dung chương trình bày số nghiên cứu phương m trình vi phân đại số có trễ lý thuyết điều khiển mà hai tác giả L.F Shampine, an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th si P Gahinet đưa báo: “Delay differential algebraic equations in control theory” năm 2006 Ở thể rõ hệ bất biến tuyến tính mơ hình lý thuyết điều khiển, hệ tuyến tính tổng quát hóa bao gồm trạng thái trễ Kết mơ hình mơ tả hệ DDAE, người ta tận dụng triệt để tính chất đặc biệt mơ hình sinh từ ứng dụng lý thuyết điều khiển, để giải tốn khó với độ xác cao tốc độ đủ nhanh nhằm đạt mục đích nghiên cứu, ứng dụng lu an n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th si bước nhảy đặc trưng điểm D Sau đó, lấy tích phân hai thành phần hàm thu liên tục xấp xỉ, có gián đoạn đạo hàm cấp thấp Trong trường hợp thứ nhất, nguyên nhân gián đoạn lan truyền bước nhảy, có gián đoạn bắt nguồn từ hàm đầu vào tuyến tính khúc mẫu FOH Chúng ta biểu diễn thành phần nghiệm hàm đa thức bậc khúc Ngoại trừ điểm D , xấp xỉ x(t ) có đạo hàm cấp liên tục xấp xỉ z (t ) liên tục tn , tn1 Xấp xỉ x(t ) xác định hàm nội suy lu an Hermite bậc với giá trị độ nghiêng hai đầu mút khoảng n va Xấp xỉ z (t ) hàm nội suy bậc với giá trị đầu mút khoảng Các phần sau cung cấp thêm thông tin chi tiết Nếu tất ie gh tn to hai giá trị điểm p thuận lợi, hàm nội suy bảo đảm độ xác công thức Runge- nl w Kutta ẩn điểm lưới trường hợp mở rộng định nghĩa oa nghiệm dạng số cho hàm trơn khúc Chúng ta thay d nghiệm vào phương trình đánh giá thỏa mãn chúng Nếu an lu lm ul xác định nf va sử dụng trường hợp để nghiệm dạng số, hàm số dư R(t ) z at nh oi X (t ) AX (t ) B1u(t ) B2 W(t ) R t W (t ) Z1 (t 1 ), , Z N (t N ) T z gm @ Theo quan điểm phép phân tích sai số lùi, mơ tả nghiệm l dạng số gần nghiệm toán nhiễu đo chất lượng m co nghiệm theo cỡ nhiễu R(t ) Đây cách đo sai số tự nhiên suốt khoảng an Lu lấy tích phân định nghĩa rõ ràng trước cơng thức RungeSố hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 48 si Kutta đánh giá Đây điểm quan trọng cỡ số dư biện pháp có ý nghĩa độ xác nghiệm đánh giá nghiệm với độ tin cậy cỡ bước lớn nghiệm có nhiều gián đoạn đạo hàm cấp thấp Trong tình khó để tính tốn xác phép đo cỡ số dư sử dụng chuẩn tích phân Cụ thể, sử dụng chuẩn RMS có lấy trọng số tính giá trị gần tích phân với cơng thức cầu phương Lobatto điểm Theo cách xây dựng, số dư hai đầu mút bước, địi hỏi phải đánh giá số dư lần bước Để cho thuận tiện, giá trị lu gọi t L1 , t L2 , t L3 Mỗi kết đánh giá giá trị dư liên quan an n va đến việc đánh giá hàm nội suy x(t ) , z (t ) , đạo hàm hàm nội suy tn to x(t ) (như xấp xỉ với x(t ) phương trình vi phân) Đối với gh trọng số thành phần, sử dụng biện pháp để đo cỡ p ie thừa số tương ứng nghiệm, sử dụng năm giá trị x(t ) khoảng nl w bước Tất nhiên vào cuối bước giá trị ln có sẵn giá trị oa thấp giá trị tới hạn có sẵn, giá trị có từ tính tốn d số dư Một chuẩn tích phân tính xấp xỉ cơng thức cầu phương Lobatto an lu nhiều điểm đem lại chương trình đặc biệt mạnh mẽ đáng tin cậy để đánh nf va giá sai số nghiệm lm ul Các xấp xỉ cho giá trị độ nghiêng x(t ) hai đầu mút bước z at nh oi xác định hàm nội suy Hermite bậc 3, mà xấp xỉ x(t ) suốt bước Điều cho phép tính xấp xỉ, mà có đạo hàm liên tục z gm @ trừ trường hợp có tác động đặc biệt xem xét để tạo bước nhảy Tương tự vậy, xấp xỉ z (t ) thành lập điểm Lobatto để đánh giá l m co điều khiển sai số sử dụng để xác định xấp xỉ z (t ) suốt bước Số hóa Trung tâm Học liệu – ĐHTN an Lu Cụ thể, hàm nội suy Hermite bậc xác định phép nội suy điểm http://www.lrc.tnu.edu.vn n va ac th 49 si tn , tL , t L , tn 1 Kết gần liên tục trừ trường hợp có tác động đặc biệt để tạo bước nhảy Tại bước điều khiển kết cỡ bước chuẩn RMS có trọng số số dư Khi tất diễn tốt đẹp cỡ bước nói đến liên hệ cỡ số dư với cỡ sai số địa phương (xem [16]) Sau thực phép tính, tính tốn nghiệm cho thỏa mãn phương trình cách có nghĩa cơng thức khơng theo trật tự thơng thường cỡ bước lớn hay nghiệm không đủ trơn Một sai số tương đối cho phép 10 3 mã hóa cương để phù hợp với lu ứng dụng Người dùng cụ thể hóa vectơ có sai số cho phép theo giá trị mặc an 6 n va định 10 Các tính tốn thu với sai số tương đối gh tn to 2.6 Các bước “dài” Bây giả sử cỡ bước không lớn để công thức p ie vừa xem xét ẩn tuyến tính Có tốn phát sinh xử lý thông qua cỡ bước thật khơng thực tế chút độ trơn w oa nl nghiệm không địi hỏi điều Với trường hợp mục 2.7.1 lan truyền d thảo luận ví dụ choppy, có 0.0001 T f Nếu giới lu nf va an hạn cỡ bước đến , cần 50.000 bước để giải toán này, lm ul code ddaeresp giải 343 bước Khi lấy bước dài trễ ngắn nhất, thành phần z j (t j ) có đối số z at nh oi tn , tn1 Chúng ta chưa tính z j (t ) khoảng này, nghiệm dạng số xác định ẩn Điều thân phương pháp Runge- z gm @ Kutta tường minh l Để thực cách đơn giản sử dụng bước lặp dài m co Đầu tiên dự đoán z (t ) tn , tn1 cách ngoại suy xấp xỉ an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 50 si thuộc tn1 , tn Trong đoạn này, hàm sử dụng để tính giá trị z (t ) ba điểm mới, sử dụng để xác định xấp xỉ đa thức bậc cho z (t ) bước tại, cụ thể tL , t L , tn 1 Sử dụng giá trị này, hàm khác đánh giá w(t ) đối số cần để tính bậc công thức Runge-Kutta Với giai đoạn này, tính xn 1 , xn 1 sau tính zn1 Một hàm nội suy Hermite bậc cho x(t ) đánh giá để nhận xấp xỉ x(t ) x(t ) nút công thức cầu phương Lobatto nhằm đánh giá sai số Một xấp xỉ cho w(t ) đánh giá điểm tương tự Những giá trị giá trị lu an u(t ) x(t ) sau sử dụng để làm tròn xấp xỉ tới z (t ) tL , t L Nếu n va cỡ bước không lớn bước hồn thành chương trình tiếp tn to tục để đánh giá sai số Nếu không, chương trình kiểm tra thay đổi tương 3 việc tính tốn lặp lặp lại lên đến lần Rõ p ie gh đối zn1 Nếu lớn 10 ràng thuận lợi đánh giá u (t ) bên ngồi vịng lặp Các khoản tiết kiệm oa nl w hữu ích hồn tồn với cỡ bước tương đương với ước lượng d gián đoạn công thức cách kiểm tra xem gián đoạn ngầm nf va lặp an lu định hay chưa, liệu gián đoạn phải ước lượng lần lm ul Bởi thời gian chạy quan trọng ứng dụng, phải z at nh oi thường xuyên sử dụng phần mềm MATLAB để xác định nơi mà tích phân tự động ưu tiên tham số chi phí thời gian chúng Thường nơi khơng phải nơi dễ thấy Với trường hợp hàm giá trị w(t ) cổ chai Bởi z gm @ đối số w j (t ) z j (t j ) phụ thuộc vào j , giá trị w(t ) không vectơ co l Các biến đại số z (t ) xếp dạng hàm đa thức bậc khúc, m cần phải tìm vị trí khúc đại diện cho thành phần quan an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 51 si tâm sau nội suy Lagrange để đánh giá hàm Hermite bậc Bài toán trở nên phức tạp xuất bước nhảy gián đoạn Mặc dù mã hóa hàm cách cẩn thận, ứng viên sáng giá cho việc mã hóa C 2.7 Các gián đoạn Mơ hình GENLTI thường xuất nhiều bước nhảy gián đoạn z (t ) x(t ) gốc lan truyền suốt khoảng xét Chắc chắn khơng cịn điều nên làm theo dõi gián đoạn lời giải số cho toán phần mà nghiệm liên tục Việc tốn lu chí đến mức nhiều code khơng theo dõi gián đoạn có an va lựa chọn cho việc Chúng ta có thỏa hiệp Có lý vật lý n để nghĩ mơ hình GENLTI quan tâm, cỡ bước nhảy gh tn to giảm mơ hình hóa xuất Đây giả thuyết làm việc Nếu p ie giả thiết khơng có giá trị, tức mơ hình gặp tượng số bước nhảy tăng hay nghiệm không ổn định, hy vọng người dùng sớm dừng nl w việc chạy chương trình Đầu tiên xét làm mà bước nhảy d oa lan truyền sau nghiên cứu cách mà chương trình xử lý nf va 2.7.1 Sự lan truyền an lu bước nhảy lm ul Những nghiên cứu từ hàm đầu vào, thấy z at nh oi giả sử u (t ) x(t ) liên tục t Để thuận tiện, sử dụng ký hiệu chuẩn z (t 0) z (t 0) z (t ) cho cỡ bước nhảy t Sử z dụng phương trình (2.2) để thấy điểm t D gm @ (2.5) l z (t ) D22 w(t ) m co - Nếu D22 Trong tình (2.5) cho thấy bước nhảy an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 52 si z (t ) không lan truyền tất điểm Có hai tình khác z (t ) khơng có bước nhảy Nhìn chung hàm đầu vào, u(0) cho khơng có bước nhảy u (t ), tương ứng khơng có bước nhảy z (t ) biến đại số liên tục 0, T f Một tình phổ biến nhiều D21 Bước nhảy tượng trưng hàm đầu vào t bị triệt tiêu (2.2) hệ biến đại số liên tục 0, T f Trong thuật giải, nhận ba trường hợp đưa biến đổi đặc biệt cho chúng Trong trường hợp D22 vô hướng, từ (2.5) suy gián đoạn lu an biến đại số giảm nhanh theo tốc độ mũ D22 va n - Nếu D22 1, bước nhảy tăng nhanh theo tốc độ mũ mơ tn to gh hình khơng có ý nghĩa mặt vật lý Trường hợp tổng quát D22 phức p ie tạp nhiều có kết hạn chế Thường hầu hết oa nl w thành phần w(t ) (2.5) Theo định nghĩa, d điểm t 1 , , t N D , điểm khác thể có thuộc khơng lu nf va an thuộc tập Nếu t k D số hạng zk t k z liên tục lm ul Chúng ta chưa tìm thấy cách khai thác triệt để trường hợp tổng quát nhận số kết hay cách xét tình đặc biệt, z at nh oi có thành phần khác Cho D , có bước nhảy z j ( ) Bước nhảy z gm @ lan truyền tiếp đến j , 2 j , , theo định nghĩa điểm co l thuộc D Sự lan truyền gián đoạn phức tạp gián đoạn khác m lan truyền đến điểm khác trễ khác Bây an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 53 si nghiên cứu tình đơn giản có vài điểm khơng có cách khác để đạt lan truyền gián đoạn Đối với điểm có thành phần thứ j w(t ) Khi có z j (t ) D22 j j z j (t j ) - Nếu D22 j j điều cho thấy ảnh hưởng trễ j giảm nhanh theo tốc độ mũ Chú ý phân rã bảo toàn với mức lan truyền khơng phải khoảng lấy tích phân lu Kết cuối giúp hiểu tình khó an mặt tính toán, cụ thể j nhỏ nhiều so với trễ khác Một gián va n đoạn lan truyền đến , 2 , tăng lên khoảng gh tn to độ trễ nhỏ thứ hai Nếu D22 j j , thấy ảnh p ie hưởng nhỏ trễ j giảm nhanh theo tốc độ mũ Một ví dụ cụ thể kiểm nl w tra toán choppy Nhiệm vụ phải tính bước phản hồi hệ với d oa biến đầu vào, năm biến vi phân, ba biến đại số, biến đầu Các trễ lần 0.6047 0 D22 0 0.0216 0 0.5457 0.7887 nf va an lu lượt tương ứng 0.3, 0.0001, 0.0230 z at nh oi lm ul Trễ ngắn 0.0001 nhỏ so với T f 5, phải theo dõi tất gián đoạn khoảng việc lấy tích phân khơng z gm @ thực tế Vì phần tử 2, D22 nên trễ không ảnh hưởng nhiều đến co l việc lấy tích phân lời giải khơng gặp khó khăn mơ hình m này, giải 1.25s Tình hồn tồn khác trễ ngắn an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 54 si trễ thứ ba phần tử 3, 3 D22 đủ lớn cho ảnh hưởng trễ tương đối dài Với nhỏ so với T f , nên chưa thể xác định rõ ràng liệu theo dõi tất gián đoạn có thực tế không Ngay với việc theo dõi hạn chế thuật giải chúng ta, phải 12s để lấy tích phân tốn mà t 2.7.2 Theo dõi gián đoạn Như thấy, việc trở lại nghiệm cho thuật giải bước nhảy gián đoạn, lan truyền suốt khoảng tích phân quan tâm đặc biệt Điều làm phức tạp cho việc mã hóa, thuật giải theo dõi lu an gián đoạn, nghiệm dạng số có bước nhảy điểm gián đoạn t Biểu va n diễn nghiệm dạng số không thuận lợi áp dụng, ta phải trở lại tn to nghiệm, chia lân cận điểm gián đoạn thành hai khoảng kề ie gh t * , t * bên trái bên phải t sau lấy giới hạn trái phải điểm p t cách tương ứng Như thấy hình kết đồ thị liên w oa nl tục có góc nhọn điểm gián đoạn d Trong dự kiến thuật giải, thực lu nf va an tính tốn D với D t0 t1 tM tập hợp điểm gián đoạn tất mức sau thực phép xây dựng nghiệm D theo lm ul dõi tất gián đoạn, thấy tˆ tr j điểm gián đoạn z at nh oi * * tk , tk 1 kết luận mơ hình GENLTI có z nghiệm Biến đại số z (t ) thường gián đoạn điểm thuộc D , @ gm lại có đạo hàm liên tục cấp d Biến vi phân x(t ) liên tục lại có co l điểm gián đoạn đạo hàm cấp điểm thuộc D m Mọi thứ diễn thuận lợi gặp phải vấn đề cạn kiệt an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 55 si nhớ Sau đó, mã hóa tính tốn để tạo điểm gián đoạn khối theo yêu cầu Không may ứng dụng có nhiều điểm gián đoạn việc theo dõi tất điểm gián đoạn không thực tế Chúng ta tin ứng dụng này, tốn ổn định có nhiều bước nhảy, mà bước nhảy làm giảm độ lớn tiến hành tích phân Vì vậy, theo dõi bước nhảy dừng theo dõi chúng trở nên đủ nhỏ Theo quan điểm phép phân tích sai số lùi, mô tả nghiệm dạng số nghiệm toán nhiễu đo chất lượng nghiệm theo cỡ nhiễu R(t ) Đây cách lu đo sai số tự nhiên suốt khoảng lấy tích phân định nghĩa rõ ràng an trước cơng thức Runge-Kutta đánh giá Cách tiếp cận va n điểm quan trọng cỡ số dư giải pháp có ý nghĩa độ gh tn to xác nghiệm đánh giá nghiệm với độ tin cậy cỡ p ie bước lớn nghiệm gián đoạn lấy đạo hàm cấp thấp Trong tình w khó khăn sử dụng chuẩn tích phân để tính tốn phép đo cỡ số dư oa nl cách xác Cụ thể, sử dụng chuẩn RMS có lấy trọng số tính d giá trị gần tích phân với cơng thức cầu phương Lobatto điểm sau lu nf va an dừng theo dõi nghiệm trơn khúc vị trí điểm nối với nhau, nghiệm liên tục xấp xỉ Hình ảnh cho thấy chí lm ul chương trình dành lượng thời gian không nhiều để theo dõi gián z at nh oi đoạn Để thuật giải hiệu công cụ nghiên cứu, thấy phải điều chỉnh việc theo dõi gián đoạn z gm @ Chương trình phải thực để lan truyền gián đoạn tới số tối đa mức định sẵn Maxlevel = Hầu hết hy vọng l co theo dõi gián đoạn thành phần nghiệm liên tục xấp xỉ Như m giải thích mục 2.7.1 lan truyền, ba tình nhận biết an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 56 si cách dễ dàng theo nghiệm liên tục sau có gián đoạn lan truyền lần Trong tình này, lấy mức Maxlevel = Ngay với giới hạn mức 4, có nhiều gián đoạn, tiếp tục hạn chế số điểm gián đoạn mà theo dõi Tập hợp điểm gián đoạn D xây dựng nhờ gián đoạn lan truyền liên tục đến mức khác cách xóa lần lặp Nếu mức nào, thấy có nhiều MaxPoints = 100 điểm, phải hủy công thức Trong phần tiếp theo, D tập bị cắt ngắn điểm gián đoạn Khi theo dõi gián đoạn, bước thuật giải bước đến điểm lu thuộc D tính tốn bước nhảy x(t ) , z (t ) , y (t ) Giả sử p điểm cuối an n va D t max( , , N ) Tại điểm t p phương trình khơng phụ thuộc tn to vào giá trị biến điểm t p t Để chuẩn bị cho trình theo ie gh dõi gián đoạn t p thuật giải bước đến điểm gián đoạn , p khơng tính bước nhảy điểm Điều làm p p t w oa nl cho thuật giải nhận xuất gián đoạn lựa chọn cỡ bước d bước khử bước nhảy nghiệm dạng số Từ t p thuật giải lu nf va an tự sử dụng cỡ bước xuất để cung cấp độ xác độ ổn định Chúng ta sử dụng hình 2.4 để kiểm tra hiển thị xảy lm ul hạn chế việc theo dõi gián đoạn 0.3 0.2 0.1 z at nh oi y(t) z @ -0.1 -0.2 gm -0.3 -0.4 l -0.5 -0.6 co -0.7 m -0.8 -0.9 an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 57 si 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 t Hình 2.4 Ví dụ artefact3 Các đường cong thấp tính tốn với việc theo dõi gián đoạn thường lệ đường cong cao tính tốn sử dụng theo dõi kéo dài Đường cong xấp xỉ với y (t ) phía xấp xỉ với y (t ) 0, Chúng ta sử dụng toán thử artefact3 để xuất giới hạn việc theo dõi gián đoạn Việc mô tả chi tiết tốn khơng cần thiết khơng có ý nghĩa cần dùng hình ảnh để mơ tả Mục đích phải tính tốn bước tương ứng hệ xác định 0,1 lu an mà có biến đầu vào, năm biến vi phân, hai biến đại số biến đầu Trễ n va 0.3000, 0.0513 trễ hai biến vào Các đường tn to cong nằm hình 2.4 hiển thị biến đầu y (t ) tính với giá trị gh thường Maxlevel MaxPoints Trong hình 2.4 đường cong nằm cho thấy p ie xảy thuật giải theo dõi đến 50 lần nhiều mức cho phép w 50 lần nhiều gián đoạn Để dễ dàng so sánh mật độ gián đoạn oa nl hai lần lấy tích phân, thêm 0.2 vào nghiệm tốt trước vẽ d Như kỳ vọng, gián đoạn sắc - nhọn (sharp) bị chặn nhiều sau an lu ngừng theo dõi Tuy nhiên, xử lý gián đoạn phần thứ nf va tích phân gián đoạn cịn tương đối lớn thu lm ul nghiệm theo nghĩa số dư không đáng kể Khi điều xảy ra, khơng tốn tính nghiệm tốn thực tế với việc tiếp tục theo dõi, z at nh oi xuất nhu cầu cần thiết phải hạn chế theo dõi thuật giải để chúng trở nên hữu ích cơng cụ thiết kế z m co l gm @ an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 58 si KẾT LUẬN Qua việc nghiên cứu báo: “Delay differential algebraic equations in control theory” tác giả L.F Shampine P Gahinet tìm hiểu số tài liệu liên quan để trình bày luận văn Em trình bày số kết kiến thức sở nghiệm dạng số hệ tuyến tính bất biến tổng quát dạng x(t ) Ax(t ) B1u (t ) B2 w(t ) z (t ) C2 x(t ) D21u (t ) D22 w(t ) lý thuyết điều khiển sau trình bày phương pháp số cho hệ tuyến tính bất lu biến tổng quát giải chúng cách nhanh việc tích phân trực tiếp an hệ sử dụng code đặc biệt, cho hàm đầu vào tương ứng va n phương trình vi phân đại số có trễ lý thuyết điều khiển to tn Tiếc nhiều vấn đề chưa làm sáng tỏ phương ie gh trình vi phân đại số có trễ chẳng hạn: Bài tốn cương, tính A - ổn định, tính ổn p định tính xác phương pháp số, thuật giải công cụ nl w MATLAB Vì chưa có điều kiện thời gian em hy vọng vấn đề d oa tiếp tục nghiên cứu thời gian tới nf va an lu z at nh oi lm ul z m co l gm @ an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 59 si TÀI LIỆU THAM KHẢO Tiếng Việt Nguyễn Thị Thu Hường (2001), “Phương pháp số giải phương trình vi phân đại số”, Luận văn thạc sĩ Nguyễn Văn Mạnh (2008), “Phương pháp bước giải phương trình vi phân cấp 2”, Tạp chí Khoa học Công nghệ ĐHTN, (45) Nguyễn Văn Minh, Trần Thanh Tùng (2016), “Phương pháp bước giải phương trình vi phân y f ( x, y ) ’’, Tạp chí Khoa học Cơng nghệ ĐHTN lu an Vũ Ngọc Phát (2001), Nhập môn lý thuyết điều khiển Toán học, Nhà xuất n va Đại học Quốc Gia Hà Nội, 41-70 gh tn to Tiếng Anh p ie Ascher U., Petzold L (1995), “The numerical solution of delay-differential- Algebraic Equations of retarded anh neutral type”, SIAM J Numer Anal, 32 oa nl w 1635- 1657 Ascher U., Petzold L (1998), Computer Methods for Ordinary Differential d an lu Equations and Differential-Algebraic Equations, SAIM, Philadelphia, PA, nf va ISBN 978- 1- 61197- 139- lm ul Corwin S P., Sarafyan D., Thompson S., DKLAG6 (1997), “A code based on continuously imbedded sixth order Runge-Kutta methods for the solution z at nh oi of state dependent functional differential equations”, Appl, Numer, Math, 24 319-333 z Enright W H., Hayashi H (1997), “A delay differential equations solver @ gm based on a continuous Runge-Kutta method with defect control”, Numer, m co l Algorithms, 16 349-364 an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 60 si Gahinet P., Shampine L F (2004), “Software for Modeling and Analysis of Linear Systems with Delays”, Proc American Control conf Boston, 5600-5605 10 Guglielmi N., Hairer E (2001), “Implementing Radau IIA methods for stiff delay differential equations”, Computing 67 1-12 11 Ixaru L Gr., Rizea M (1997), “Four step methods for y f ( x, y ) ’’ Journal of Computational and Applied Mathematics, 79 87-99 12 Lambert J D (1991), “Nummerical Methods in Ordinary Differential Equations”, Wiley Chichester, 978-0-470-72335-7 13 Luzyanina T., Roose D (2004), Numerical local stability analysis of lu an differential algebraic equations with time delays, in: Fifth IFAC Workshop va n on Time-Delay Systems, K.U Leuven, Belgium, submitted for publication vol 283, Mathematics department, University of Manchester, UK, ISSN ie gh tn to 14 Paul C A H (1995), “A User-Guide to ARCHI”, Numer, Anal Rept., p 1360-1725 w 15 Roussel R Marc (2005), “Delay-differential equations”, p1-12 oa nl 16 Shampine L F (2005), “Solving ODEs and DDEs with residual control”, d Appl Numer Math, 52 113-127 lu nf va an 17 Shampine L F., Gladwell I., Thompson S (2003), “Solving ODEs with MATLAB”, Cambridge University Press, New York, 1-21 lm ul 18 Shampine L F., Gahinet P (2006), “Delay-differential-algebraic equations p574-588 z at nh oi in control theory”, Applied Numerical Mathematics, vol56, ISSUE 3-4, 19 Stephen L (2008), “Differential-algebraic equations”, Campbell et al z m co l gm @ Scholarpedia, 3(8):2849 an Lu Số hóa Trung tâm Học liệu – ĐHTN http://www.lrc.tnu.edu.vn n va ac th 61 si lu an n va p ie gh tn to d oa nl w nf va an lu z at nh oi lm ul z m co l gm @ an Lu n va ac th 62 http://www.lrc.tnu.edu.vn Số hóa Trung tâm Học liệu – ĐHTN si
Ngày đăng: 24/07/2023, 09:27
Xem thêm:
