Giáo trình Cơ sở điều khiển tự động: Phần 1 cung cấp cho người đọc những kiến thức như: đại cương về hệ thống điều khiển; mô tả toán học hệ thống điều khiển liên tục; đặc tính động học; ổn định của hệ thống tự động tuyến tính; chất lượng hệ thống điều khiển tự động. Mời các bạn cùng tham khảo!
TS NGUYỄN ĐỨC THÀNH (CHỦ BIÊN) ThS NGUYỄN VĂN SƠN ThS VÕ THÀNH NHÂN CƠ SỞ ĐIỀU KHIỂN TỰ ĐỘNG NHÀ XUẤT BẢN ĐẠI HỌC THỦ DẦU MỘT 2015 TRƢỜNG ĐẠI HỌC THỦ DẦU MỘT Khoa Điện - Điện Tử Nguyễn Đức Thành Nguyễn Văn Sơn Võ Thành Nhân CƠ SỞ ĐIỀU KHIỂN TỰ ĐỘNG NHÀ XUẤT BẢN ĐẠI HỌC THỦ DẦU MỘT BÌNH DƢƠNG 2015 LỜI NĨI ĐẦU Quyển sách “Cơ sở điều khiển tự động” trình bày nội dung Lý thuết điều khiển tự động, nội dung thường giảng dạy trường Đại học Cao đẳng chuyên ngành Kỹ thuật điều khiển Tự động hoá, Điện tử- Tự động, Cơ điện tử số chuyên ngành khác Nội dung sách gồm chín chương Chương 1: Đại Cương Hệ Thống Điều Khiển Chương Mô Tả Toán Học Hệ Thống Điều Khiển Liên Tục Chương Đặc Tính Động Học Chương Ổn Định Của Hệ Thống Tự Động Tuyến Tính Chương Chất Lượng Hệ Thống Điều Khiển Tự Động Chương Thiết Kế Bộ Điều Khiển Liên Tục Chương Hệ Thống Điều Khiển Rời Rạc Chương Thiết Kế Bộ Điều hiển ời ạc Chương Hệ Thống Điều Khiển Phi Tuyến Cuối chương có câu hỏi tập Cơng cụ Matlab dùng xuyên suốt ví dụ nhằm giảm thời gian tính tốn giúp sinh viên tiếp cận phương pháp tính tốn đại Thạc sĩ Nguyễn Văn Sơn soạn chương 8, Thạc sĩ Võ Thành Nhân soạn chương 9, Tiến sĩ Nguyễn Đức Thành soạn phần lại sách chỉnh lý tồn nội dung Quyển sách dùng làm tài liệu giảng dạy/ tham khảo cho mơn học có tên liên quan Cơ sở điều khiển tự động, Cơ sở tự động học, Lý thuyết điều khiển tự động, Kỹ thuật điều khiển tự động… Tuỳ theo lượng thời gian phân bố cho môn học trình độ sinh viên, số chương giảng dạy, số chương dành cho sinh viên tự đọc Với khối lượng 30 tiết lý thuyết, nội dung giảng dạy lớp chương 2, 3, 4, 5, 6, với số đề mục nhỏ cho sinh viên tự đọc Với khối lượng 45 tiết lý thuyết/bài tập, nội dung giảng dạy lớp chương 1,2, 3, 4, 5, 6, Với khối lượng 60 tiết lý thuyết/bài tập, nội dung giảng dạy lớp chương 1, 2, 3, 4, 5, 6, Nhóm biên soạn có nhiều cố gắng nhiên cịn nhiều thiếu sót, mong lượng thứ Các đóng góp chân thành người đọc cho sách tốt xin gởi địa chỉ: hoa Điện- Điện tử, Trường Đại Học Thủ Dầu Một, Trần Văn Ơn, Thành phố Thủ Dầu Một, Bình Dương Bình Dương, 20-10-2015 Nhóm biên soạn MỤC LỤC Chương1 Đại Cương Hệ Thống Điều Khiển 1.1 Giới thiệu 1.1.1 Các định nghĩa 1.1.2 Mục đích lý thuyết điều khiển tự động 1.2 Phân Loại HTĐK 1.3 Lịch Sử Phát Triển ĐKTĐ 1.4 Ví Dụ HTĐKTĐ Chương Mơ Tả Tốn Học Hệ Thống Điều Khiển Liên Tục 2.1 Biến Đổi Laplace 16 2.1.1 Định nghĩa 16 2.1.2 Các tính chất biến đổi Laplace 16 2.1.3 Bảng biến đổi Laplace 18 2.1.4 Biến đổi Laplace ngược 19 2.1.5 Phân tích phân thức tối giản 19 2.2 Hàm Truyền 23 2.2.1 Hàm truyền hệ ngõ vào, ngõ SISO 23 2.2.2 Hàm truyền hệ nhiều ngõ vào nhiều ngõ MIMO 28 2.3 Phương Trình Trạng Thái 29 2.3.1 Vế phải phương trình vi phân khơng có đạo hàm tín hiệu vào 29 2.3.2 Vế phải phương trình vi phân có đạo hàm tín hiệu vào 30 2.3.3 Dạng tắc điều khiển 33 2.3.4 Dạng tắc quan sát 34 2.3.5 Chuyển đổi phương trình trạng thái sang hàm truyền 37 2.3.6 Định lý Cayley-Hamilton 39 2.4 Sơ Đồ Khối 41 2.4.1 Khái niệm 41 2.4.2 Đại số sơ đồ khối 42 2.5 Sơ Đồ Dịng Tín Hiệu 46 2.5.1 Định nghĩa 46 2.5.2 Công thức Mason 46 2.6 Tuyến Tính Hố 48 2.7 Mơ Hình Tốn Hệ Thống Cơ Điện 50 2.7.1 Hệ thống 50 2.7.2 Con lắc ngược 53 2.7.3 Động chiều nam châm vĩnh cửu 56 2.7.3 Động chiều không chổi than (Brushless DC Motor BLDC) 59 2.7.4 Động AC 59 2.8 Hàm Truyền Hệ Thống Nhiệt 61 2.8.1 Các đại lượng nhiệt 61 16 2.8.2 Hàm truyền lò nhiệt điện trở 63 2.8.3 Hàm truyền máy nước nóng 64 2.9 Hàm Truyền Hệ Thống Mực Chất Lỏng 65 2.9.1 Hệ bồn chứa nước 66 2.9.2 Hệ hai bồn nước 67 2.10 Tương Đồng Giữa Các Hệ Thống 68 2.11 Hệ Thống Có Trễ 69 Kết luận 69 Bài tập 70 Chương Đặc Tính Động Học 3.1 Đáp Ứng Thời Gian 74 3.2 Đáp Ứng Tần Số 77 3.2.1 Biểu đồ Bode 78 3.2.2 Biểu đồ Nyquist 81 3.2.3 Biểu đồ Nichols 81 3.3 Đặc Tính Động Học Khâu Bậc Nhất 82 3.3.1 Khâu tỷ lệ 82 3.3.2 Khâu quán tính 82 3.3.3 Khâu tích phân 84 3.3.4 Khâu vi phân 84 3.3.5 Khâu 1+sT 85 3.4 Đặc Tính Động Học Khâu Bậc Hai 86 3.4.1 Khâu bậc hai tiêu chuẩn 86 3.4.2 Khâu bậc hai có zero 91 3.4.3 Khâu tích phân quán tính 93 74 3.5 Khâu Trì Hỗn e sTD 96 3.5 Vẽ Biểu Đồ Bode Hàm Truyền Tổng Quát 96 Kết luận 100 Bài tập 100 Chương Ổn Định Của Hệ Thống Tự Động Tuyến Tính 4.1 Khái Niệm Về Ổn Định 103 4.2 Tiêu Chuẩn Routh-Hurwitz 105 4.2.1 Tiêu chuẩn Routh 107 4.2.2 Tiêu chuẩn Hurwitz 114 4.3 Quĩ Tích Nghiệm 114 4.3.1 Định nghĩa 114 4.3.2 Quy tắc vẽ quỹ tích nghiệm 115 4.4 Tiêu Chuẩn Ổn Định Nyquist 121 4.5 Tiêu Chuẩn Ổn Định Bode 127 103 Kết luận 128 Bài tập 129 Chương Chất Lượng Hệ Thống Điều Khiển Tự Động 132 5.1 Khái Niệm 132 5.2 Sai Số Xác Lập 133 5.2.1 Sai số xác lập với hàm nấc 133 5.2.2 Sai số xác lập với hàm dốc 135 5.2.3 Sai số xác lập với hàm parabole 136 5.2.4 Trường hợp H(s) ≠ 139 5.2.5 Điều kiện hàm truyền kín để triệt tiêu sai số xác lập 139 5.3 Triệt Tiêu Nhiễu Loạn 140 Độ Quá Điều Chỉnh Và Thời Gian Xác Lập 144 5.4.1 Thêm zero vào hàm truyền hệ bậc hai 145 5.4.2 Hàm truyền bậc ba 148 5.5 Các Hàm Truyền Tiêu Chuẩn 152 5.5.1 Các hàm tiêu 152 5.5.2 Hàm truyền tối ưu 153 5.6 Khảo Sát Chất Luợng Dựa Trên Biểu Đồ Bode 155 Kết luận 156 Bài tập 156 Chương Thiết Kế Bộ Điều Khiển Liên Tục 6.1 Khái Niệm 159 6.2 Điều Khiển PID 163 6.2.1 Điều khiển tỷ lệ P 163 6.2.2 Điều khiển PI 168 6.2.3 Thiết kế điều khiển PI hệ bậc hai dùng quỹ tích nghiệm 170 6.2.4 Điều khiển PID 172 6.2.5 Thiết kế PID dùng phương pháp giải tích 173 6.2.6 Thiết kế dùng tiêu chuẩn ITAE 175 6.2.7 Phương pháp Ziegler Nichols 176 6.2.8 Thiết kế PID dùng Matlab 180 6.2.9 Các điều khiển PID công nghiệp 182 6.3 Điều Khiển Sớm Trễ Pha 183 6.3.1 Thiết kế hiệu chỉnh sớm pha 186 6.3.2 Thiết kế hiệu chỉnh trễ pha 190 6.3.3 Thiết kế hiệu chỉnh sớm-trễ pha 194 6.4 Thiết Kế Dùng Matlab 199 6.5 Thiết Kế Bộ Điều Khiển Nhiều Vòng Hồi Tiếp 205 6.6 Thiết Kế Bộ Điều Khiển Dùng Phương Pháp Đặt Cực 208 6.6.1 Tính ma trận hồi tiếp 208 159 6.6.2 Điều khiển tín hiệu theo tín hiệu vào 213 6.6.3 Điều khiển hồi tiếp trạng thái ngõ 215 6.7 Hệ Thống Quan Sát- Điều Khiển Trạng Thái 220 6.7.1 Quan sát trạng thái đầy đủ 220 6.7.2 Kết hợp điều khiển quan sát trạng thái đầy đủ 222 6.7.3 Quan sát giảm bậc 231 6.7.4 Đặt cực cho hệ MIMO 232 6.8 Điều Khiển LQR 233 Kết luận 236 Bài tập 236 Chương Hệ Thống Điều Khiển Rời Rạc 7.1 iới Thiệu Hệ Thống Điều Khiển Rời Rạc 238 7.1.1 Lấy mẫu tín hiệu 238 7.1.2 ộ giữ ậc ero 239 7.1.3 Thực điều khiển số 240 7.2 iến Đổi Z 241 7.2.1 iến đổi v tính chất 241 7.2.2 iến đổi ngược 244 7.3 H m Truyền Hệ Thống Rời Rạc 249 7.3.1 Biến đổi phương trình sai phân h m truyển 249 7.3.2 Biến đổi hàm truyền liên tục sang rời rạc 249 7.3.3 H m truyền v ng kín rời rạc 251 7.4 Phân Tích Ổn Định Tr n Mặt Ph ng Z 253 7.4.1 Ti u chuẩn ổn định Routh-Hurwit cho hệ thống rời rạc 255 7.4.2 Ti u chuẩn ổn định Jury 256 7.5 Quỹ Tích Nghiệm Của Hệ Thống Điều Khiển Rời Rạc 257 7.6 Đáp Ứng Tần Số Hệ Thống Điều Khiển Rời Rạc 259 7.7 Hệ Phương Trình Trạng Thái Hệ Thống Điều Khiển Rời Rạc 261 7.7.1 Dạng tắc điều khiển 262 7.7.2 Dạng tắc quan sát 264 7.7.3 Dạng tắc đường chéo 265 7.7.4 Hệ thống điều khiển 267 ` 7.7.5 Hệ thống quan sát 268 7.7.6 Chuyển hệ phương trình trạng thái li n tục sang rời rạc 268 7.7.7 H m truyền rời rạc t phương trình trạng thái 270 Kết luận 271 Bài tập 271 Chương Thiết Kế Bộ Điề hiển ời ạc 8.1 Chất Lượng Hệ Thống Điều Khiển Rời Rạc 278 8.2 Hiệu Chỉnh PID 281 8.2.1 ộ hiệu chỉnh PID 281 238 278 8.2.2 Thực hiệu chỉnh PID ng vi điều khiển 288 8.3 ổ Chính Sớm Pha –Trễ Pha 290 8.3.1 Bổ sớm pha rời rạc 290 8.3.2 ổ trễ pha rời rạc 292 8.4 Thiết kế giải tích mặt ph ng z 297 8.5 Điều khiển đặt cực 300 8.5.1 Thiết kế đặt cực 301 8.5.2 Đặt cực triệt tiêu sai số xác lập 304 8.6 Đáp Ứng Dead eat 306 8.7 ộ Quan Sát Trạng Thái 310 8.7.1 Thiết kế quan sát 310 8.7.2 ộ quan sát v điều khiển đặt cực 313 8.8 ộ Quan Sát iảm ậc 314 8.9 Điều Khiển LQR 319 Kết luận 322 Bài tập 322 Chương Hệ Thống Điều Khiển Phi Tuyến 9.1 Khái Niệm ề Hệ Thống Phi Tuyến 325 9.1.1 Khái Niệm 325 9.1.2 Tính chất v đặc điểm riêng phi tuyến 325 9.1.3 Các phương pháp khảo sát hệ phi tuyến 326 9.2 Phương Pháp H m Mơ Tả Tuyến Tính Hóa Điều H a) 327 9.2.1 Khái Niệm 327 9.2.2 H m mô tả hay hệ số khuếch đại phức khâu phi tuyến 328 9.2.3 H m mô tả khâu phi tuyến điển hình 328 9.2.3 Chu trình giới hạn 331 9.2.4 Xét ổn định chu trình giới hạn 332 9.3 Phương Pháp Mặt Ph ng Pha 334 9.3.1 Điểm bất thường 335 9.3.2 Vẽ chân dung pha 336 9.3.3 Chân dung pha hệ tuyến tính 337 9.3.4 Chân dung pha hệ phi tuyến 338 9.4 Ti u Chuẩn Ổn Định Lyapunov 339 9.4.1 Ổn định theo nghĩa Lyapunov 339 9.4.2 Định lý Lyapunov thứ 341 9.4.3 Định lý trực tiếp Lyapunov 341 9.5 Thiết Kế Bộ Điều Khiển Phi Tuyến 346 9.5.1 Phương pháp hồi tiếp tuyến tính hóa 347 9.5.2 Phương pháp điều khiển thích nghi theo mơ hình mẫu (MRAC) 347 Kết luận 349 Bài tập 349 Tài liệu tham khảo 325 CHƯƠNG ĐẠI CƯƠNG VỀ HỆ THỐNG ĐIỀU KHIỂN 1.1 GIỚI THIỆU 1.1.1 Các định nghĩa Hệ thống (System): tập hợp phần tử có quan hệ hữu với nhau, tác động chi phối lẫn theo quy luật định, ví dụ hệ thống sưởi, hệ thống điều hòa, hệ thống giao thông, hệ thống điện, hệ thống ngân hàng, hệ thống giáo dục, hệ thống trị, hệ thống tuần hồn sinh vật… Điều khiển (Control): tác động lên hệ thống để đạt mục tiêu đó, ví dụ điều khiển vận tốc hướng xe chạy đường, điều khiển phanh chống bó (ABS AntiLock Braking System), điều khiển tần số điện áp máy phát, diều khiển giao thông Hệ thống điều khiển (Control System) (HTĐK): bao gồm Đối tượng điều khiển, gọi tắt Đối tượng (Object) Bộ điều khiển (Controller), Bộ điều khiển giám sát tác động đến đối tượng thông qua thay đổi lượng cung cấp cho đổi tượng, ví dụ tăng cơng suất điện lị nóng hơn, người lái xe thay đổi vận tốc xe cách thay đổi lượng hịa khí đốt Đối tượng đa dạng, hệ thống vật lý ta có đối tương cơ, nhiệt, điện…; kinh tế, sinh học… có nhiều loại đối tượng khác nhau.Thuật ngữ tiếng Anh thường dùng từ Plant, Process để Đối tượng Điều chỉnh (Regulate) trì giá trị ổn định cho trạng thái hệ thống, điều chỉnh có phạm vi hẹp điều khiển, ta có điều chỉnh, hệ thống điều chỉnh, ví dụ hệ thống điều chỉnh thân nhiệt sinh vật , ổn áp 5V ổn định tần số lưới điện Nhiễu loạn (Disturbance): tác động không mong muốn lên hệ thống gây giảm chất lượng điều khiển, ví dụ đóng mở cửa vào phịng có điều hồ, thay đổi tải động Bộ điều khiển cần phải lấy thông tin đối tượng môi trường xung quanh cảm biến để tác động phù hợp lên đối tượng Nếu việc điều khiển có tham gia người, ta gọi điều khiển tay (manual control), ví dụ lái xe có tài xế; khơng có tham gia người ta có điều khiển tự động, ví dụ xe ô tô không người lái Giáo trình khảo sát hệ thống điều khiển tự động Đối tượng hệ thống có ngõ vào tín hiệu điều khiển ngõ tín hiệu điều khiển, điều khiển nhận ngõ vào tín hiệu đặt xuất tín hiệu tín hiệu điều khiển Hệ thống điều khiển có ngõ vào tín hiệu đặt tín hiệu tín hiệu điều khiển, bám theo tín hiệu đặt (Hình 1.1) Các tín hiệu vào vơ hướng, ta gọi hệ thống SISO (Single Input Single Output) hệ thống có tín hiệu vào tín hiệu ra, trường hợp ngược CƠ SỞ ĐIỀU KHIỂN TỰ ĐỘNG CHƢƠNG lại tín hiệu vào vector, hệ thống gọi MIMO (Multiple Input Multiple Output) hay hệ thống đa biến (Multivariable system) Hình 1.1 Hệ thống điều khiển 1.1.2 Mục đích lý thuyết điều khiển tự động Lý thuyết điều khiển tự động (LTĐKTĐ) khảo sát hành vi đối tượng đưới tác động tín hiệu vào điều khiển, tính tốn thiết kế điều khiển để bảo đảm tín hiệu đạt mong muốn tín hiệu đặt LTĐKTĐ kết hợp tảng toán học kỹ thuật, áp dụng cho kỹ thuật phát triển thêm cho ngành khác xã hội, kinh tế, tài chính, trị, sinh lý học, y sinh… Điều khiển tự động quan trọng khoa học kỹ thuật, công nghiệp dân dụng Các toán thường gặp điều khiển điều khiển vị trí, vận tốc, nhiệt độ, áp suất, độ ẩm, lưu lượng… Nội dung LTĐKTĐ bao gồm ba vấn đề chính: - Nhận dạng hệ thống (Identification): xây dụng mơ hình tốn đối tượng cần điều khiển thơng qua tín hiệu kích thích tín hiệu - Phân tích hệ thống (Analyse): Khảo sát đáp ứng với tín hiệu vào khác dựa vào mơ hình tốn đối tượng - Thiết kế điểu khiển (Design Controller): Tinh toán điều khiển để tín hiệu hệ thống đạt yêu cầu tín hiệu đặt 1.2 PHÂN LOẠI HTĐK HTĐK chia làm HTĐK vòng hở HTĐK vòng kín - HTĐK vịng hở (Hình 1.2) tác động đến ngõ vào đối tượng thơng qua tín hiệu điều khiển, tín hiệu khơng phụ thuộc ngõ thực tế, dịng tín hiệu hệ thống theo chiều, chịu tác động nhiễu loạn bên ngoài, tín hiệu thay đổi khơng mong muốn, ví dụ người lái xe để tay ga cố định cho đường vận tốc 40km/g xuống dốc vận tốc tăng lên lên dốc vận tốc giảm xuống, lấy ví dụ khác điều khiển vận tốc động chiều cách đưa điện áp cố định vào động cơ, có nhiễu loạn tải tác động vận tốc quay động giảm CHƢƠNG Hình 5.16 Đáp ứng điều khiển vận tốc với nhiễu loạn Ví dụ 5.6 ĐỘ QUÁ ĐIỀU CHỈNH VÀ THỜI GIAN XÁC LẬP K Cho hệ thống với hàm truyền đối tượng G(s)= , dùng điều khiển tỷ lệ s( s a ) hàm truyền kín Gk ( s ) KKc Chương khảo sát đặc tính động học khâu s as KKc n2 bậc hai hàm truyền tùy theo KKc ta đặc tính động học s 2 n s n2 khác với n2 KKc , a=2ζωn Nhắc lại kết chương 3: - Phần trăm vọt lố : POT= 100 e - Thời gian xảy đỉnh vọt lố tmax - Thời gian tăng tr 1 % , n 0.8 2.5 - Thời gian xác lập 2% n ts n - Độ vọt lố thời gian xác lập ngược nhau, hệ thống có vọt lố nhỏ thời gian xác lập tăng lên - Hệ thống có băng thơng lớn tác động nhanh - Đối với hệ bậc hai có zero hay hệ bậc cao hơn, có nhiều thơng số tham gia nên khó tìm cơng thức tính tốn hai thơng số quan trọng này, sau khảo sát số trường hợp cụ thể 144 CHẤT LƢỢNG HỆ THỐNG ĐIỀU KHIỂN TỰ ĐỘNG Hình 5.17 Cực Đáp ứng nấc khâu bậc hai 5.4.1 Thêm zero vào hàm truyền hệ bậc hai a/ Thêm zero vào hàm truyền hở hệ bậc hai KKc ( s b) Xét hàm truyền vòng thuận Gc(s)G(s)= , trường hợp xảy s( s a) - điều khiển P với G(s)= K ( s b) hay s( s a) - điều khiển PD với G(s)= K hay s( s a) - điều khiển PI với G(s)= K sa hay - điều khiển sớm trễ pha với G(s)= K s Hàm truyền kín s n2 ( 1) KKc ( s b) b Gk(s)= s ( KKc a) s bKKc s 2 n s n2 Ta nhận thấy zero s=-b xen vào tử mẫu hàm truyền hệ kín đặc tính động học hệ kín phụ thuộc vị trí zero so với cực Dựa vào quỹ tích nghiệm Hình 5.18 ta nhận thấy b < a hệ kín có hai cực thật âm cịn b > a có trường hợp cực phức KKc có giá trị nhỏ Vậy b< a khơng có vọt lố, cực tiến zero –a, cực lại tiến - ∞ Khi b>a, hệ số đệm thay đổi theo KKc, cực thực hay phức, ta chọn giá trị phù hợp tùy theo tiêu chuẩn thiết kế Cực hệ kín là: 145 CHƢƠNG Hình 5.18 Qũy tích nghiệm hệ kín có hàm truyền vòng thuận K ( s b) s( s a) ( KKc a) ( KKc ) a 2(a 2b) KKc Giả sử chọn cực kép để đáp ứng khơng vọt lố nhanh (KKc+a)2=4b KKc , hai giá trị tương ứng nghiệm phương trình (KKc)2+(2a-4b)KKc+a2=0, KKc)=2b-a± (2b a)2 a 2b-a ± b(b a) , cực kép - b(2b a b(b a) ) Chọn dấu + đáp ứng nhanh có vọt lố Chương tính tốn POT hệ Tóm lại đưa thêm zero vào hàm truyền thuận (điều khiển PD) ảnh hưởng tốt đến đặc tính động học hệ kín, giảm thời gian xác lập, giảm vọt lố so với điều khiển P Hình 5.19 Đáp ứng nấc hệ kín có hàm truyền vịng thuận 146 K ( s b) s( s a) CHẤT LƢỢNG HỆ THỐNG ĐIỀU KHIỂN TỰ ĐỘNG , chọn điều khiển PD có hàm s( s 2) 2( K P K D s) K D ( s K P / K D ) , hàm truyền Gc K P K D s , hàm truyền vòng hở Gh s(s 2) s(s 2) Ví dụ 5.7 : Hàm truyền đối tượng G n2 (1 sK D / K p ) 2( K P K D s) KD truyền kín Gk , n2 K P , , zero s (2 K D )s K P s 2n s n 2KP KP , chọn zero ≥ -2 khơng có vọt lố, chọn KP lớn đáp ứng nhanh, giả sử chọn KD 20 s 20 Nếu chọn zero nhỏ -2 có đáp ứng có vọt lố, ví dụ KP=20, KD=5, 40 10 s Gk , ωn=6.246 rad/s, ζ=0.9487, dùng công thức ta POT=6.7% thời s 12 s 40 điểm 0.3927 sec Ta viết chương trình Matlab vẽ đáp ứng với trường hợp điều khiển P điều khiển PD, nhận thấy điều khiển PD tốt điều khiển P (Hình 5.20) KP=10, KD=5, triệt tiêu cực hệ hở lúc Gk close all kd=5; for kp= [0.5 16 ] ; t=0:0.01:5; g=tf([2],[1 0]); gc=tf([kd kp],[1]); gk=feedback(g*gc,1); step(gk,t) ; hold on; end title( 'Unit-step responses of PD control 2nd order system' grid on ) Hình 5.20 Đáp ứng nấc điều khiển P PD 147 CHƢƠNG b/ Thêm zero vào hàm truyền kín bậc hai Nếu so sánh hàm truyền hệ thống kín có zero khơng có zero ta nhận thấy zero tử số hàm truyền làm tăng vọt lố giảm thời tăng đáp ứng nấc (1 Tz s) Xét hàm truyền kín Gk ( s) n , dùng Matlab ta vẽ đáp ứng nấc s 2n s n2 Hình 5.21 với thông số ωn=1, ζ=0.5, Tz=0, 1, 2, 4, 6, 10 5.4.2 Hàm truyền bậc ba Đối với hệ thống có hàm truyền thuận GcG K (b2 s b1s b0 ) , hệ kín có ba a3s a2 s a1s a0 cực, ảnh hưởng đến đáp ứng hệ thống ổn định Nếu a0=0 hệ kín có ssxl với hàm nấc triệt tiêu Việc tính toán POT ts phức tạp, phụ thuộc vào vị trí cực zero Nhìn chung hệ bậc ba có thời gian xác lập lớn hệ bậc hai, nều hàm truyền có zero vọt lố tăng Hình 5.21 Đáp ứng nấc hàm truyền hệ kín bậc hai có zero a/ Thêm cực vào hàm truyền thuận bậc hai Xét hàm truyền hở bậc hai thêm khâu quán tính Gh n2 , ta s( s 2n )(1 Tp s ) khảo sát ảnh hưởng cực thêm vào đáp ứng nấc hệ kín điều khiển tỷ lệ K Pn2 Hàm truyền hệ kín Gk Tp s (2nTp 1) s 2n s K Pn2 Khâu qn tính thêm vào làm hệ thống ổn định, làm tăng vọt lố thời tăng tr Điều khơng có lợi cho chất lượng hệ thống Hệ ổn định 148 CHẤT LƢỢNG HỆ THỐNG ĐIỀU KHIỂN TỰ ĐỘNG 2 (2nTp 1) Tp KPn Dùng Matlab khảo sát đáp ứng nấc với ωn=1, ζ=1, KP=1, Tp=0,1,2, (Hình 5.22) b/ Thêm cực vào hàm truyền kín bậc hai n2 Hàm truyền kín có dạng Gk , cực thêm vào -1/Tp làm giảm ( s 2n s n2 )(1 Tp s) vọt lố tăng thời gian xác lập Với thông số ωn=1, ζ=0.5, Tp=0,0.5 1,2, ta Hình 5.23 Tăng Tp vọt lố giảm thời gian xác lập tăng, ngược với trường hợp mục a Hình 5.22 Đáp ứng nấc thêm cực vào hàm truyền thuận hệ bậc hai Hình 5.23 Đáp ứng nấc thêm cực vào hàm truyền kín bậc hai c/ Điều khiển tích phân hệ bậc hai khơng có zero Xét đối tượng có hàm truyền G , dùng điều khiển tích phân ta có s as b pωn2 K , hàm truyền kín có dạng Gk (s) , khơng GcG s( s as b) ( s 2ζωn s ωn2 )( s p) có zero, K tăng hệ kín ổn định (Hình 5.24), điều khơng có lợi cho hệ thống d/ Hệ bậc ba có zero, điều khiển PI Thêm zero vào điều khiển ta có điều khiển PI, GcG truyền kín Gk (s) K ( s b0 ) , hàm s( s a1s a0 ) K ( s b0 ) , lúc hệ kín có cực tiến -b0, , đáp ( s 2ζωn s ωn2 )( s p) ứng hệ kín tùy thuộc vị trí zero cực hệ hở giá trị K 149 CHƢƠNG Hình 5.24 Quỹ tích nghiệm hàm truyền vịng khơng có zero Khảo sát quỹ tích nghiệm Hình 5.25, 5.26 cho thấy zero tiến bên trái hệ thống ổn định, nên chọn zero nằm bên phải cực hệ hở, lúc hệ kín ln ln ổn định Điều chứng minh qua tiêu chuẩn Routh, phương trình đặc trưng s( s a1s a0 ) + K ( s b0 ) = s3 a1s (a0 K )s Kb0 =0 Điều kiện ổn định a1(a0+K)>Kb0 hay a1a0>K(b0- a1) Điều kiện luôn thỏa mãn b0< a1, lúc vế phải bất đẳng thức âm cịn vế trái dương Vì a1=-(tổng hai cực hệ hở), b0=-zero hệ hở nên chọn zero nằm bên phải cực hệ hở, hệ kín ln ln ổn định Nếu b0≥ a1thì tùy giá trị K hệ kín ổn định Dùng Matlab ta vẽ đáp ứng nấc cho hai hệ Hình 5.27 d/ Hệ bậc ba có hai zero Dùng điều khiển PID , Gc K ( s b1s b0 ) K ( s b1s b0 ) , GcG s s( s a1s a0 ) , hàm truyền hệ kín có phương trình đặc trưng s (a1 K )s (a0 Kb1 )s Kb0 Hệ thống kín có ba cực, cực tiến -∞, hai cực tiến hai zero, cực zero hệ hở thực âm, hệ kín ln ổn định, khơng có vọt lố Điều kiện ổn định (a1 K )( a0 Kb1 ) Kb0 K 2b1 K (a0 b1a1 b0 ) a1a0 Chọn b0 a1 s=tf(‘s’); g1=(s+1.5)/(s*(s+2)*(s+4)); gk1=feedback(8*g1,1); g2=(s+1.5)/(s*(s^2+4*s+8)); gk2=feedback(8*g2,1); subplot(1,2,1); rlocus(g1) subplot(1,2,2); rlocus(g2) figure subplot(1,2,1); step(gk1) subplot(1,2,2) step(gk2) Hình 5.27 Đáp ứng nấc hệ bậc ba 151 CHƢƠNG 5.5 CÁC HÀM TRUYỀN TIÊU CHUẨN Ở mục trước khảo sát chất lượng hệ thống điều khiển dựa sai số xác lập, độ vọt lố thời gian xác lập Mục có mục đích đánh giá chất lượng hệ thống điều khiển dựa cực tiểu hàm tiêu liên quan đến sai số 5.5.1 Các hàm tiêu Khi thiết kế hệ thống điều khiển ta thường phải chọn số thông số điều khiển để làm cực tiểu hàm tiêu có dạng T J f (r (t ) , y (t ) ,e(t ) ,t )dt Khi hệ thống điều khiển gọi hệ thống tối ưu theo tiêu J Sau số tiêu thông dụng dựa sai số i/ Chỉ tiêu ISE (Integral of square of error) T Cực tiểu tích phân bình phương sai số ISE e ( t )dt ii/IAE (Integral of absolute error) T IAE e( t )dt iii/ITAE (Integral of time multiplied by absolute error) T ITAE e( t )tdt iv/ITSE (Integral of time multiplied by squared error) T ITSE te2 ( t )dt Ví dụ 5.8: Cho hệ thống với hàm truyền kín Gk ( s ) 2.5 , tìm K cực tiểu s Ks 2.5 tiêu chuẩn ISE y( t ) e n t 1 sin(d t ) 10 10 K Sai số e( t ) 152 10 10 K e Kt 10 K sin( t ) e Kt 10 K sin( t ) CHẤT LƢỢNG HỆ THỐNG ĐIỀU KHIỂN TỰ ĐỘNG Bình phương sai số e2 ( t ) 10 10 K Kt e sin ( t ) 10 K 10 Chỉ tiêu ISE J= e ( t )dt 10 K 2 Kt J e (1 cos( 10 K t 2 ))dt 0.2 K 2K Đạo hàm J theo K cho ta K= 10 n tương ứng hệ số đệm ζ=0.5 5.5.2 Hàm truyền tối ưu Để thuận tiện cho việc sử dụng, có bảng hàm truyền kín tối ưu theo tiêu ITAE, hàm truyền khơng có zero sai số xác lập với hàm nấc không a0 Gk ( s ) n n 1 s an 1s an s n a0 Sau mẫu số hàm truyền đến bậc (D Graham and R.C Lathrop,1953) Bảng 5.2 s n s 1.4n s n2 s3 1.75n s 2.15n2 s n3 s 2.1n s 3.4n2 s 2.7n3s n4 s5 2.8n s 5n2 s3 5.5n3s 3.4n4 s n5 Bảng 5.2 cơng bố năm 1953, đáp ứng nấc có độ vọt lố cao, sau nhiều tác giả đưa bảng 5.3 (Y Cao, 1989) có đáp ứng tốt Bảng 5.3 s n s 1.505n s n2 s3 1.783n s 2.172n2 s n3 s 1.953n s3 3.347n2 s 2.648n3s n4 s5 2.068n s 4.499n2 s 4.675n3s 3.257n4 s n5 153 CHƢƠNG Hình 5.28 Đáp ứng nấc hàm truyền hệ kín cực tiểu tiêu chuẩn ITAE từ bậc đến Đối với tín hiệu vào hàm dốc ta có hàm truyền kín a1s a0 Gk ( s ) n n 1 s an 1s an s n a0 Đa thức mẫu số bảng 5.4: Bảng 5.4 s 3.2n s n2 s3 1.75n s 3.25n2 s n3 s 2.41n s3 4.93n2 s 5.14n3s n4 s5 2.19n s 6.5n2 s3 6.30n3s 5.24n4 s n5 Nguyên tằc thiết kế điều khiển theo ITAE ta chọn điều khiển cho hàm T (s) truyền kín có dạng Gk ( s) với M(s) đa thức mẫu số giống bảng hàm tiêu chuẩn M (s) sau dùng tiền lọc hàm truyền P( s ) nj T (s) ngõ vào để triệt tiêu zero (đối với đáp ứng nấc), ý dùng tiền lọc làm tăng thời gian xác lập hệ thống 154 CHẤT LƢỢNG HỆ THỐNG ĐIỀU KHIỂN TỰ ĐỘNG Hình 5.29 Thiết kế điều khiển với lọc ngõ vào 5.6 KHẢO SÁT CHẤT LUỢNG DỰA TRÊN BIỂU ĐỒ BODE Các mục khảo sát chất lượng hệ thống dựa vào hàm truyền đáp ứng thời gian Trong mục ta xét liên quan chất lượng với biểu đồ Bode Ta xét hệ kín hồi tiếp đơn vị có sơ đồ Hình 5.30, vẽ biểu đồ Bode vịng hở G(jω), biểu đồ Bode biên cắt trục hoành tần số ωcg | G(jωcg)|=1, với hệ bậc hai: Hình 5.30 n2 G( s) s( s 2n ) cg n 4 2 cg 90o tan 1 2n 2 Góc pha tương ứng G( jcg ) 90o tan 1 Dự trữ pha PM=180o+ G( jcg ) =tan-1 4 2 2 4 2 PM phụ thuộc hệ số đệm Dùng Matlab ta vẽ PM theo ζ (Hình 5.31) 155 CHƢƠNG Hình 5.31 Dự trữ pha theo hệ số đệm Ta nhận thấy ζ