Kỹ thuật điều khiển nâng cao - TS.Nguyễn Viễn Quốc

Chia sẻ: Le Tan Kiep | Ngày: | Loại File: PDF | Số trang:29

Thêm vào BST

Báo xấu

149
lượt xem 22
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Mời các bạn cùng tham khảo tài liệu kỹ thuật điều khiển nâng cao để nắm được các nội dung kiến thức cần thiết và hỗ trợ trong việc học môn này được tốt hơn.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Kỹ thuật điều khiển nâng cao - TS.Nguyễn Viễn Quốc

KỸ THUẬT ĐIỀU KHIỂN NÂNG CAO Giảng viên: TS. Nguyễn Viễn Quốc Email: vienquoc@gmail.com Nội dung Chương 1: Điều khiển dựa trên mô hình trạng thái Chương 2: Điều khiển tối ưu Chương 3: Điều khiển mờ Chương 4: Mạng nơron nhân tạo Tài liệu tham khảo 1. Rolands S. Burns, Advanced Control Engineering, 2001 (Chapter 8, 9, 10) 2. Nguyễn Thị Phương Hà, Lý thuyết Điều khiển Hiện đại, NXB ĐHQG, 2012. 1
Chương 1: Điều khiển dựa trên mô hình trạng trái 1.1) Mô hình trạng thái 1.1.1) Định nghĩa - Xét hệ thống bậc n, m ngõ vào, p ngõ ra. m ngõ vào Hệ thống p ngõ ra … … - Mô hình không gian trạng thái (state-space model) gọi tắt là “mô hình trạng thái” có dạng: & x = Ax + Bu (1.1)a y = Cx (1.2)b trong đó: x : vectơ trạng thái, (n ´ 1) u : vectơ ngõ vào, (m ´ 1) y : vectơ ngõ ra, (p ´ 1) A : ma trận hệ thống (hay ma trận trạng thái), (n ´ n) B : ma trận vào, (n ´ m) C: ma trận ra, (p ´ n) D : ma trận truyền thẳng (feedthrough matrix), (p ´ n) - Phương trình 1.1a gọi là phương trình trạng thái; phương trình 1.1b gọi là phương trình ngõ ra. - Sơ đồ khối chi tiết hệ thống: 2
- Trong trường hợp ngõ ra hệ thống không chịu tác động trực tiếp từ ngõ vào thì D = 0 (ma trận không) - Hệ thống bậc n có n biến trạng thái. Các biến trạng thái thay đổi - Trường hợp hệ bậc 2 (n = 2), 2 trạng thái ( ), ( ) vẽ nên quỹ theo thời gian “vẽ” nên một quỹ đạo trong không gian n chiều. đạo gọi là chân dung pha (phase portrait) trong không gian 2 chiều gọi là mặt phẳng pha (phase plane). 1.1.2) Thành lập mô hình trạng thái - 2 cách thiết lập mô hình trạng thái: o Từ phương trình vi phân o Từ hàm truyền. VD: Tìm mô hình trạng thái của hệ thống trong hình dưới đây: … 3
VD: … 1.1.3) Xác định hàm truyền, zero, cực từ mô hình trạng thái: ̇ = + - Xét hệ thống có mô hình trạng thái: = det( − )=0 - Cực của hệ thống là nghiệm phương trình: Phương trình trên được gọi là phương trình đặc tính của hệ thống. - Zero của hệ thống là nghiệm phương trình - Hàm truyền được xác định theo công thức: 4
( ) ( )= = ( − ) ( ) 0 1 1 ̇ = + VD: Xác định cực - zero và hàm truyền của hệ thống: −2 −1 1 = [1 0] Giải:… 1.1.4) Tính không duy nhất của mô hình trạng thái: - Xét hệ thống: Ta cũng có thể biểu diễn hệ thống trên theo vectơ trạng thái z qua = phép đổi biến: ̇= + trong đó, là ma trận n x n không suy biến. = ̇= + Suy ra = Hệ thống trên có thể được biểu diễn bởi mô hình trạng thái theo ̇= + biến trạng thái z như sau: = = trong đó: = 5
= - Như vậy, tùy vào cách đặt biến trạng thái mà một hệ thống có thể được biểu diễn bởi nhiều mô hình trạng thái khác nhau. 0 1 0 ̇ = + VD: Cho hệ thống: −6 −5 1 = [6 0] 1 1 = 1 2 Đổi biến trạng thái với: 1 1 = Giải 1 2 1 1 Đổi biến trạng thái với: = 1 2 13 20 = = −12 −18 −1 = = 1 = = [6 6] 15 22 −1 ̇= + Phương trình trạng thái theo biến trạng thái z: −13 −19 1 = [6 6] 1.2) Lời giải phương trình trạng thái ̇( ) = ( )+ ( ) - Xét hệ bậc nhất: ( ) − (0) = ( ) + ( ) Biến đổi Laplace 2 vế: (0) ⇒ ( )= + ( ) − − 6
Biến đổi Laplace ngược, ta có: ( )= (0) + ( ) ( ) có thể được xác định theo công thức: = 1+ + + ⋯+ 2! ! ̇( ) = ( )+ ( ) - Tương tự, xét hệ bậc > 1: ( ) − (0) = ( )+ ( ) Biến đổi Laplace 2 vế: ⇒( − ) ( ) = (0) + ( ) ⇒ ( )=( − ) (0) + ( − ) ( ) Biến đổi Laplace ngược, ta có: ( )= (0) + ( ) ( ) Đặt ( )= , phương trình trên trở thành: ( )= ( ) (0) + ( − ) ( ) ( ) gọi là ma trận chuyển trạng thái. Nó có thể được tính theo các công thức sau: = + + + ⋯+ 2! ! hoặc =ℒ {( − ) } VD: Xét hệ thống mô tả bởi phương trình trạng thái: 7
Cho m = 1 kg, C = 3 Ns/m, K = 2 N/m, u(t) = 0. Hãy xác định: ( ) và ( ) a) Phương trình đặc tính, cực hệ thống, và b) (0) = 1.0; c) Đáp ứng quá độ của các biến trạng thái từ điều kiện đầu: ̇ (0) = 0 Giải a) PTĐT: b) c) u(t) = 0 8
1.3) Điều khiển được – Quan sát được 1.3.1) Điều khiển được ̇ = + - Xét hệ thống bậc n có mô hình trạng thái: = - Hệ thống được gọi là điều khiển được (controllable) nếu tồn tại ( ) đến trạng thái ( ) bất kỳ trong khoảng thời gian hữu hạn. luật điều khiển u(t) có thể đưa hệ thống từ trạng thái ban đầu - Điều kiện để hệ thống điều khiển được (toàn phần) là ma trận =[ … ] điều khiển được , Lưu ý: rank( ) = Û det ( ) ≠ 0. có hạng là n. VD: Cho hệ thống: 9
Hệ thống có điều khiển được không? −1 2 4 ̇= + VD: Cho hệ thống: −3 4 6 = [1 −2] Hệ thống có điều khiển được không? =[ ]= 4 8 Giải 6 12 det =0⇒ ( )
VD: Cho hệ thống: Hệ thống có quan sát được không? = [1 −1] Giải: −2 0 = [1 −1] = [−5 5] 3 −5 1 −1 = −5 5 det = 0 Þ Hệ thống không quan sát được 1.4) Điều khiển hồi tiếp trạng thái ̇ = + - Xét hệ thống có mô hình trạng thái: = với luật điều khiển hồi tiếp trạng thái: = − . ̇ = + ( − ) - Phương trình trạng thái của hệ thống kín: 11
Þ ̇=( − )+ det( − + )=0 - Phương trình đặc trưng của hệ kín: Nghiệm của phương trình này chính là cực hệ kín. Bài toán Phân bố cực (Pole Placement): - Nếu hệ thống điều khiển được, ta có thể chọn K để hệ kín có cực ở bất kỳ vị trí nào mong muốn. - Các phương pháp xác định K: o So sánh trực tiếp o Công thức Ackerman 1.4.1) Xác định K bằng phương pháp so sánh trực tiếp: - Nếu vị trí cực mong muốn là , , …, , ta có phương trình đặc ( )=( − )( − )⋯( − )= trưng mong muốn là: + +⋯+ + =0 - Đồng nhất các hệ số của phương trình trên với phương trình đặc det( − + )=0 trưng hệ kín: Ta sẽ có hệ n phương trình để xác định n phần tử của ma trận . 0 1 0 ̇= + VD: Cho hệ thống hở: 0 −4 1 = [1 0] = = −2. Xác định ma trận hồi tiếp trạng thái K để hệ thống kín có cực tại Kiểm tra tính điều khiển được: 12
=[ ]= 0 1 1 −4 det ≠ 0 Þ rank( ) = 2 Þ hệ thống điều khiển được. ( ) = ( + 2)( + 2) = +4 +4 =0 Phương trình đặc trưng mong muốn: (1) 1 0 0 1 0 [ ] =0 det − + Phương trình đặc trưng hệ thống hồi tiếp trạng thái: 0 1 0 −4 1 Þ + (4 + ) + = 0 (2) =4 So sánh (1) và (2): 4+ =4Þ =0 Ma trận hồi tiếp trạng thái cần tìm: = [4 0] 1.4.2) Xác định K bằng công thức Ackermann: = [0 0 ⋯ 0 1] ( ) - Công thức Ackermann: ( )= + +⋯+ + trong đó: là ma trận điều khiển được, ( : là hệ số phương trình đặc trưng mong muốn) VD: Như VD trên, nhưng xác định bằng công thức Ackermann. = [0 1] ( ) Giải: =[ ]= 0 1 Þ = 4 1 1 −4 1 0 ( )= + + 0 1 0 1 1 0 4 0 = +4 +4 = 0 −4 0 −4 0 1 0 4 13
4 1 4 0 = [0 1] = [4 0] Suy ra: 1 0 0 4 0 1 0 0 VD: Cho hệ thống hở: ̇= 0 0 1 + 3 −4 −7 −3 1 = [0 0 1] cặp cực phức với = 0.6, = 10 và cực thứ ba là cực thực tại Hãy xác định bộ điều khiển hồi tiếp trạng thái sao cho hệ kín có −20. 1.5) Bộ quan sát (ước lượng) trạng thái - Luật điều khiển hồi tiếp trạng thái đòi hỏi phải đo được tất cả các trạng thái. Trong thực tế, thường là không đo được tất cả các trạng thái do thiếu cảm biến hoặc có những trạng thái không phải là những đại lượng vật lý nên không thể đo được ® cần phải ước lượng trạng thái. - Việc ước lượng trạng thái được thực hiện bởi bộ ước lượng trạng thái, hay còn gọi là bộ quan sát. o Bộ quan sát đủ bậc: ước lượng tất cả các biến trạng thái. o Bộ quan sát giảm bậc: ước lượng một vài trạng thái trong số các trạng thái của hệ thống. - Bộ quan sát Luenberger: 14
̇ = + Xét hệ thống: = (1.5.1) Gọi là trạng thái ước lượng. Phương trình trạng thái bộ quan ̇ = + + ( − ) sát: (1.5.2) trong đó là ma trận độ lợi của bộ quan sát. ̇ − ̇ = ( − )− ( − ) Trừ 2 phương trình (1.5.1) và (1.5.2), ta có: = − , phương trình (3) viết lại như sau: (1.5.3) ̇ =( − ) Đặt (1.5.4) ̇ =( − ) + + Từ (2), phương trình trạng thái bộ quan sát được viết lại như sau: (1.5.5) 15
định bởi trị riêng của − Nhận xét: Từ pt (1.5.4), ta thấy đặc tính động học của quyết , cũng chính là cực của bộ quan sát (1.5.5). - Việc thiết kế bộ quan sát chính là xác định để cực của bộ quan sát ở vị trí mong muốn. Thông thường cực của bộ quan sát cần phải được thiết kế sao cho nó nhanh hơn cực của bộ điều khiển ít nhất là 5 lần. Bài toán Phân bố cực (Pole Placement) cho bộ quan sát: - Nếu hệ thống quan sát được, ta có thể chọn để bộ quan sát có cực ở bất kỳ vị trí nào mong muốn. - Phương pháp xác định tương tự như pp xác định , bao gồm: o So sánh trực tiếp o Công thức Ackerman 1.5.1) Xác định Ke bằng phương pháp so sánh trực tiếp ( )=( − )( − )⋯( − )= - Phương trình đặc trưng mong muốn + +⋯+ + =0 - Đồng nhất các hệ số của phương trình trên với phương trình đặc det( − + )=0 trưng của bộ quan sát: Ta sẽ có hệ n phương trình để xác định n phần tử của ma trận . 0 1 0 ̇= + VD: Cho hệ thống hở: −2 −3 1 = [1 0] Thiết kế bộ quan sát có = 10 rad/s và = 0,5. 16
Giải 1 0 = = Kiểm tra tính quan sát được của hệ thống: 0 1 rank( ) = 2 Þ Hệ thống quan sát được. ( )= +2 + = + 10 + 100 = 0 Phương trình đặc trưng mong muốn: det( − + )=0 Phương trình đặc trưng của bộ quan sát: 1 0 0 1 [1 0] det( − + ) = det − + 0 1 −2 −3 = + (3 + ) + (3 +2+ )=0 3+ = 10 Đồng nhất hệ số 2 phương trình, ta có: 3 +2+ = 100 = 7, = 77 7 = Giải ra: 77 Ma trận độ lợi bộ quan sát cần tìm: 1.5.2) Xác định Ke bằng công thức Ackermann = ( ) [0 0 ⋯ 0 1] - Công thức Ackermann: ( )= + +⋯+ + trong đó: là ma trận quan sát được, ( : là hệ số phương trình đặc trưng mong muốn) VD: Như ví dụ trên, nhưng xác định bằng công thức 1 0 = = Ackermann. 0 1 17
( )= + 10 + 100 = 0 Phương trình đặc trưng mong muốn: ( )= + + Suy ra: 0 1 0 1 1 0 = + 10 + 100 −2 −3 −2 −3 0 1 0 = ( ) Ma trận Ke được xác định theo công thức: 1 1.6) Tính đối ngẫu giữa bài toán thiết kế bộ điều khiển và thiết kế bộ quan sát: - Bài toán phân bố cực thiết kế bộ điều khiển: Đồng nhất hệ số 2 det( − + )= ( ) đa thức - Bài toán phân bố cực thiết kế bộ quan sát: Đồng nhất hệ số 2 đa det( − + )= ( ) thức det( − + ) = det( − + ) - Ta có: = det( − +( ) ) = det( − + ) - Nhận xét: 2 bài toán này tương tự nhau: Thiết kế bộ điều khiển Thiết kế bộ quan sát 18
1.7) Hệ thống điều khiển hồi tiếp trạng thái ước lượng - Xét hệ thống điều khiển như hình: ̇ = + Hệ thống hở: = (1.6.1) =− Luật điều khiển hồi tiếp trạng thái ước lượng: (1.6.2) thống là = − , thì: Gọi sai lệch giữa trạng thái ước lượng và trạng thái thực của hệ = − (1.6.3) ̇ = − ( − )=( − ) + Kết hợp các phương trình (1.6.1), (1.6.2), (1.6.3), ta có: (1.6.4) ̇ =( − ) Từ (1.5.4), phương trình sai lệch: (1.6.5) 19
̇ − = Gộp các phương trình (1.6.4) và (1.6.5): ̇ 0 − (1.6.6) - Phương trình (1.6.6) mô tả đặc tính động học của hệ thống điều khiển hồi tiếp trạng thái ước lượng. Từ đó ta có thể suy ra − det − =0 phương trình đặc trưng của hệ thống: 0 − Þ det( − + ) det( − + )=0 (1.6.7) Nhận xét: o Cực của hệ kín bao gồm cực của bộ quan sát và cực của hệ thống điều khiển hồi tiếp trạng thái. o Việc thêm bộ quan sát không làm thay đổi vị trí của cực hệ thống điều khiển hồi tiếp trạng thái. - Từ nhận xét trên, ta có kết luận: Bài toán thiết kế hệ thống hồi tiếp trạng thái ước lượng có thể phân tách thành 2 bài toán: o Thiết kế bộ điều khiển hồi tiếp trạng thái Þ o Thiết kế bộ quan sát bộ quan sát Þ - Cực của bộ quan sát cần phải được chọn sao cho nó ít ảnh hưởng đến cực của bộ điều khiển hồi tiếp trạng thái Þ cực bộ quan sát cần được chọn ở xa trục ảo hơn (nhanh hơn) cực bộ điều khiển. 1.8) Bộ quan sát giảm bậc: - Trong trường hợp hệ thống bậc n có m trạng thái có thể đo chính xác, chỉ cần ước lượng n-m trạng thái còn lại Þ sử dụng bộ quan - Để đơn giản, ta xét trường hợp = 1. Giả sử trạng thái sát giảm bậc. đo được chính xác. Phương trình ngõ ra: 20