intTypePromotion=3

Bài giảng Lý thuyết điều khiển hiện đại - Chương 3: Điều khiển bền vững

Chia sẻ: Lão Lão | Ngày: | Loại File: PDF | Số trang:0

0
1
lượt xem
0
download

Bài giảng Lý thuyết điều khiển hiện đại - Chương 3: Điều khiển bền vững

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

Chương này cung cấp cho người học những nội dung cơ bản như: Khái niệm điều khiển bền vững, chuẩn của tín hiệu, đại số ma trận, trị suy biến của ma trận – độ lợi chính (principal gain), ổn định nội, định lý độ lợi nhỏ (small gain theorem), phương pháp lqg (linear quadratic gaussian),... Mời các bạn cùng tham khảo để biết thêm các nội dung chi tiết.

Chủ đề:
Lưu

Nội dung Text: Bài giảng Lý thuyết điều khiển hiện đại - Chương 3: Điều khiển bền vững

  1. Chương 3 : Điều khiển bền vững Chương 3 ĐIỀU KHIỂN BỀN VỮNG 3.1 Giới thiệu 3.1.1 Khái niệm điều khiển bền vững Hệ thống điều khiển bền vững làm cho chất lượng của sản phẩm ổn định, không phụ thuộc vào sự thay đổi của đối tượng cũng như của nhiễu tác động lên hệ thống.Mục đích của điều khiển bền vững là chất lượng vòng kín được duy trì mặc dù có những sự thay đổi trong đối tượng. P0 :Mô hình chuẩn (mô hình danh định) PΔ :Mô hình thực tế với sai lệch Δ so với mô hình chuẩn Hình 3.1 : Mô hình điều khiển bền vững Cho tập mô hình có sai số PΔ và một tập các chỉ tiêu chất lượng, giả sử P0 ∈ PΔ là mô hình danh định dùng để thiết kế bộ điều khiển K.Hệ thống hồi tiếp vòng kín được gọi là có tính : - Ổn định danh định: nếu K ổn định nội với mô hình danh định P0 - Ổn định bền vững: nếu K ổn định nội với mọi mô hình thuộc PΔ - Chất lượng danh định: nếu các mục tiêu chất lượng được thỏa đối với mô hình danh định P0 Học kì 1 năm học 2005-2006
  2. PGS.TS Nguyễn Thị Phương Hà Trang 2 - Chất lượng bền vững: nếu các mục tiêu chất lượng được thỏa đối với mọi mô hình thuộc PΔ Mục tiêu bài toán ổn định bền vững là tìm bộ điều khiển không chỉ ổn định mô hình danh định P0 mà còn ổn định một tập các mô hình có sai số PΔ 3.1.2 Chuẩn của tín hiệu 3.1.2.1 Khái niệm chuẩn Trong điều khiển nói riêng cũng như trong các công việc có liên quan đến tín hiệu nói chung,thông thường ta không làm việc chỉ riêng với một tín hiệu hoặc một vài tín hiệu điển hình mà ngược lại phải làm việc với một tập gồm rất nhiều các tín hiệu khác nhau. Khi phải làm việc với nhiều tín hiệu khác nhau như vậy chắc chắn ta sẽ gặp bài toán so sánh các tín hiệu để chọn lọc ra được những tín hiệu phù hợp cho công việc. Các khái niệm như tín hiệu x1(t) tốt hơn tín hiệu x2(t) chỉ thực sự có nghĩa nếu như chúng cùng được chiếu theo một tiêu chuẩn so sánh nào đó. Cũng như vậy nếu ta khẳng định rằng x1(t) lớn hơn x2(t) thì phải chỉ rõ phép so sánh lớn hơn đó được hiểu theo nghĩa nào, x1(t) có giá trị cực đại lớn hơn , có năng lượng lớn hơn hay x1(t) chứa nhiều thông tin hơn x2(t)…..Nói một cách khác ,trước khi so sánh x1(t) với x2(t) chúng ta phải gắn cho mỗi một tín hiệu một giá trị đánh giá tín hiệu theo tiêu chuẩn so sánh được lựa chọn . Định nghĩa: Cho một tín hiệu x(t) và một ánh xạ x(t) →||x(t)|| ∈ R+ chuyển x(t) thành một số thực dương ||x(t)||.Số thực dương này sẽ được gọi là chuẩn của x(t) nếu nó thỏa mãn: a. ||x(t)|| ≥ 0 và ||x(t)|| = 0 khi và chỉ khi x(t) =0 (3.1) b. ||x(t)+y(t)|| ≤ ||x(t)|| + ||y(t)|| ∀ x(t), y(t) (3.2) c. ||ax(t)|| = |a|.||x(t)|| ∀ x(t) và ∀a ∈ R . (3.3) 3.1.2.2 Một số chuẩn thường dùng trong điều khiển cho một tín hiệu x(t): ∞ - Chuẩn bậc 1: || x(t ) ||1 = ∫ | x(t ) |dt (3.4) −∞ ∞ || x(t ) || 2 = ∫ | x(t ) | 2 - Chuẩn bậc 2: dt . (3.5) −∞ http://www.khvt.com
  3. Chương 3 : Điều khiển bền vững Bình phương chuẩn bậc hai chính là giá trị đo năng lượng của tín hiệu x(t). ∞ || x(t ) || p = ∫ | x(t ) | với p ∈ N p -Chuẩn bậc p: p dt (3.6) −∞ - Chuẩn vô cùng: || x(t ) ||∞ = sup | x(t ) | (3.7) t đây là biên độ hay đỉnh của tín hiệu Khái niệm chuẩn trong định nghĩa trên không bị giới hạn là chỉ cho một tín hiệu x(t) mà còn được áp dụng được cho cả vector tín hiệu gồm nhiều phần tử và mỗi phần tử lại là một tín hiệu. Xét một vector tín hiệu: ⎛ x1 (t ) ⎞ ⎜ ⎟ x(t) = ⎜# ⎟ ⎜ x (t ) ⎟ ⎝ n ⎠ - Chuẩn 1 của vector x: n x 1 = ∑ xi (3.8) i =1 - Chuẩn 2 của vector x: n ∑x 2 x 2 = i (3.9) i =1 - Chuẩn vô cùng của vector x: x ∞ = max xi (3.10) i =1, 2 ,..., n 3.1.2.3 Quan hệ của chuẩn với ảnh Fourier và ảnh Laplace: Để phục vụ mục đích sử dụng khái niệm chuẩn vào điều khiển ,ta cần quan tâm tới mối liên quan giữa chuẩn tín hiệu x(t) là ||x(t)|| với ảnh Fourier X(j ω ) cũng như ảnh Laplace X(s) của nó. Trang 3
  4. PGS.TS Nguyễn Thị Phương Hà Trang 4 Định lí 3.1: (Parseval) Chuẩn bậc hai của một tín hiệu x(t) và ảnh Fourier X(j ω ) của nó có quan hệ : ∞ ∞ 1 || x(t ) || 22 = ∫ | x(t ) | 2 dt = ∫ | X ( jω ) | 2 dω (3.11) −∞ 2π −∞ Cho tín hiệu nhân quả causal x(t). Gọi X(s) là ảnh Laplace của nó .Giả sử rằng X(s) có dạng thực -hữu tỷ với bậc của đa thức tử số không lớn hơn bậc đa thức mẫu số ,tức là: B( s ) b0 + b1 s + ..... + bm s m X ( s) = = với m < n (3.12) A( s ) a0 + a1 s + ..... + an s n Định lí 3.2: Xét tín hiệu nhân quả causal x(t) có X(s) dạng (3.12) .Để chuẩn bậc 1 của x(t) là một số hữu hạn ||x(t)||1= K < ∞ thì điều kiện cần và đủ là tất cả các điểm cực của X(s) phải nằm bên trái trục ảo (có phần thực âm) . 3.1.3 Đại số ma trận 3.1.3.1 Một số ma trận thường gặp: - Một ma trận A=(aij) có số hàng bằng số cột được gọi là ma trận vuông. Đường chéo nối các phần tử aii trong ma trận vuông được gọi là đường chéo chính .Đường chéo còn lại được gọi là đường chéo phụ. ⎡ a11 a12 " a1n ⎤ ⎢a a 22 " a 2 n ⎥⎥ A = ⎢ 21 (3.13) ⎢ # # # # ⎥ ⎢ ⎥ ⎣a n1 an 2 " a nn ⎦ - Một ma trận vuông A=(aij) có aij = 0 khi i ≠ j ,tức là các phần tử không nằm trên đường chéo chính đều bằng 0, được gọi là ma trận đường chéo. Ma trận đường chéo được ký hiệu bởi: ⎡a11 0 " 0 ⎤ ⎢0 a 22 " 0 ⎥⎥ A= ⎢ = diag(aij) (3.14) ⎢ # # # # ⎥ ⎢ ⎥ ⎣0 0 " a nn ⎦ http://www.khvt.com
  5. Chương 3 : Điều khiển bền vững ⎡1 0 " 0⎤ ⎢0 1 " 0⎥⎥ - Ma trận đường chéo I = diag(1) = ⎢ gọi là ma trận đơn vị. ⎢# # # #⎥ ⎢ ⎥ ⎣0 0 " 1⎦ - Ma trận vuông A=(aij) có aij = 0 khi i > j (hoặc i < j) được gọi là ma trận tam giác + Ma trận tam giác dưới ⎡ a11 0 " 0 ⎤ ⎢a a 22 " 0 ⎥⎥ A= ⎢ 21 (3.15) ⎢ # # # # ⎥ ⎢ ⎥ ⎣ a n1 an 2 " a nn ⎦ + Ma trận tam giác trên ⎡a11 a12 " a1n ⎤ ⎢0 a 22 " a 2 n ⎥⎥ A= ⎢ (3.16) ⎢ # # # # ⎥ ⎢ ⎥ ⎣0 0 " a nn ⎦ 3.1.3.2 Các phép tính về ma trận: - Phép cộng / trừ: Cho hai ma trận A=(aij) và B=(bij) cùng có m hàng và n cột .Tổng hay hiệu A ± B = C =(cij) của chúng được định nghĩa là một ma trận cũng có m hàng và n cột với các phần tử cij = aij + bij i=1,2,…..,m và j=1,2,…..,n. - Phép nhân với số thực: Cho ma trận A=(aij) có m hàng và n cột và một số vô hướng thực(phức) x tùy ý .Tích B = xA = Ax = (bij) được hiểu là ma trận cũng có m hàng và n cột với các phần tử Bij = x.aij i=1,2,….m và j=1,2,…..,n - Phép chuyển vị: Ma trận chuyển vị của ma trận A=(aij) với m hàng và n cột là ma trận AT = (aji) có n hàng và m cột được tạo từ ma trận A qua việc hoán chuyển hàng thành cột và ngược lại cột thành hàng. - Phép nhân ma trận: Cho ma trận A=(aik) có m hàng và p cột và ma trận B=(bkj) có p hàng và n cột ,tức là : Trang 5
  6. PGS.TS Nguyễn Thị Phương Hà Trang 6 + A=(aik) i=1,2,....,m và k=1,2,….,p + B=(bkj) k=1,2,….,p và j=1,2,…..,n Tích AB = C =(cij) của chúng là một ma trận có m hàng và n cột với các phần tử p Cij = ∑a k =1 b ik kj Một ma trận vuông A ∈ R n×n được gọi là ma trận trực giao nếu ATA=AAT=I 3.1.3.3 Hạng của ma trận: Cho n vector vi i=1,2,…,n Chúng sẽ được gọi là độc lập tuyến tính nếu đẳng thức a1v1+a2v2+…….+anvn=0 trong đó ai là những số thực (hoặc phức) sẽ đúng khi và chỉ khi a1 = a2 = …..=an = 0 Xét một ma trận A=(aij) bất kì có m hàng và n cột .Nếu trong số m vector hàng có nhiều nhất p ≤ m vector độc lập tuyến tính và trong số n vector cột có nhiều nhất q ≤ n vector độc lập tuyến tính thì hạng ma trận đươc hiểu là: Rank(A) = min{p,q} Một ma trận vuông A kiểu (n× n) sẽ được gọi là không suy biến nếu Rank(A)=n .Ngược lại nếu Rank(A)
  7. Chương 3 : Điều khiển bền vững Do phải tồn tại cả hai phép nhân AA-1 và A-1A cho ra kết quả có cùng kiểu nên ma trận A phải là một ma trận vuông,tức là phải có m = n.Hơn nữa do det(I) = 1 ≠ 0 nên: det(A)det(A-1) ≠ 0 => det(A) ≠ 0 và det(A-1) ≠ 0. (3.22) Vậy A phải là ma trận không suy biến. Ma trận nghịch đảo A-1 của A có tính chất sau: - Ma trận nghịch đảo A-1 của A là duy nhất (3.23) - Tập hợp tất cả các ma trận vuông cùng kiểu và không suy biến cùng với phép nhân ma trận tạo thành một nhóm (không giao hoán). (3.24) ⎡a b ⎤ 1 ⎡ d − b⎤ - Nghịch đảo ma trận kiểu (2× 2): A −1 = ⎢ ⎥ = ⎢ ⎥ (3.25) ⎣ c d ⎦ det( A) ⎣− c a ⎦ - (AB)-1 = B-1A-1 (3.26) - (A-1)T = (AT)-1 (3.27) ⎡1⎤ - Nếu A = diag(ai) và không suy biến thì A-1 = diag ⎢ ⎥ (3.28) ⎣ ai ⎦ Aadj - A-1 = (3.29) det( A) trong đó Aadj là ma trận có các phần tử a i j = (-1)i+jdet(Aij) với Aij là ma trận thu được từ A bằng cách bỏ đi hàng thứ j và như cột thứ i. - Cho ma trận A ∈ Rn × n không suy biến . Nếu U ∈ Rn × m và V ∈ Rn × m là hai ma trận làm cho (I+VTA-1U) cũng không suy biến thì (A+UVT)-1 = A-1 – A-1U(I+VTA-1U)-1VTA-1 (3.30) ⎡A A2 ⎤ - Cho ma trận vuông A = ⎢ 1 không suy biến,trong đó A1,A2,A3,A4 ⎣ A3 A4 ⎥⎦ cũng là các ma trận. Nếu A1 không suy biến và B = A4 – A3A1-1A2 cũng không suy biến thì −1 ⎡A −1 A2 ⎤ ⎡ A1 −1 + A −11 A2 B −1 A3 A1 −1 − A1 A2 B −1 ⎤ −1 A =⎢ 1 =⎢ ⎥ (3.31) ⎣ A3 A4 ⎥⎦ ⎣ − B −1 A3 A1 −1 B −1 ⎦ Trang 7
  8. PGS.TS Nguyễn Thị Phương Hà Trang 8 Nếu A4 không suy biến và C = A1 – A2A4-1A3 cũng không suy biến thì −1 ⎡A −1 A2 ⎤ ⎡ C −1 − C −1 A2 A4 −1 ⎤ A =⎢ 1 =⎢ −1 ⎥ (3.32) ⎣ A3 A4 ⎥⎦ −1 ⎣− A4 A3 AC −1 −1 −1 A4 + A4 A3C −1 A2 A3 ⎦ 3.1.3.5 Vết của ma trận: Cho ma trận vuông A=(aij) ,i,j=1,2,……,n kiểu (nxn).Vết của A được hiểu là tổng giá trị các phần tử trên đường chéo chính của A và được ký hiệu bằng trace(A): m trace= ∑ aii (3.33) i =1 Vết của ma trận có các tính chất: a. trace(AB) = trace(BA) (3.34) -1 b. trace(S AS) = trace(A) với S là ma trận không suy biến bất kì (3.35) 3.1.3.6 Giá trị riêng và vector riêng: Số thực λ được gọi là giá trị riêng và vector x được gọi là vector riêng bên phải ứng với giá trị riêng λ của A thỏa mãn: Ax = λ x ∀ x (3.36) ⇔ (A - λ I)x = 0 ∀ x (3.37) Giá trị riêng và vector riêng của ma trận A có những tính chất sau: a. Hai ma trận tương đương A và S-1AS luôn cùng giá trị riêng, nói cách khác giá trị riêng của ma trận bất biến với phép biến đổi tương đương: det(A- λ I)=det(S-1AS- λ I) (3.38) b. Các giá trị riêng của ma trận bất biến với phép chuyển vị, tức là: det(A- λ I)=det(AT- λ I) (3.39) c. Nếu A không suy biến thì AB và BA có cùng các giá trị riêng ,tức là: det(AB- λ I)=det(BA- λ I) (3.40) T d. Nếu A là ma trận đối xứng (A =A) thì các vector riêng ứng với những giá trị riêng khác nhau sẽ trực giao với nhau Trong Matlab ,sử dụng hàm eig(A) để tìm ma trận riêng và vector riêng. http://www.khvt.com
  9. Chương 3 : Điều khiển bền vững 3.1.3.7 Tính toán ma trận: Cho ma trận X = (xij) ∈ Cm × n là một ma trận thực (hoặc phức) và F(X) ∈ C là một vô hướng thực hoặc phức của X .Đạo hàm của F(X) đối với X được định nghĩa ∂ ⎡ ∂ ⎤ F(X ) = ⎢ F ( X )⎥ (3.41) ∂X ⎣⎢ ∂xij ⎦⎥ Cho A và B là những ma trận phức với không gian tương thích .Một số công thức đạo hàm : ∂ Trace ( AXB ) = AT B T (3.42) ∂X ∂ Trace ( X k ) = k(X k −1 T ) (3.43) ∂X ∂ Trace ( XBX T ) = 2 XB (B = BT ) (3.44) ∂X ∂ ( X T AX ) = AX + A T X (3.45) ∂X ∂ Trace ( AX T B ) = BA (3.46) ∂X 3.1.3.8 Chuẩn của ma trận: Người ta cần đến chuẩn của ma trận là nhằm phục vụ việc khảo sát tính giải tích của nó.Có nhiều chuẩn khác nhau cho một ma trận A=(aij) ,i=1,2,…,m;j=1,2,…,n. Những chuẩn thông thường được sử dụng: - Chuẩn 1 của ma trận A m A 1 = max ∑ aij (3.47) 1≤ j ≤ n i =1 - Chuẩn 2 của ma trận A A 2 = max λi ( A* A) (3.48) 1≤ i ≤ n - Chuẩn vô cùng của ma trận A Trang 9
  10. PGS.TS Nguyễn Thị Phương Hà Trang 10 n A ∞ = max ∑ a ij (3.49) 1≤i ≤ m j =1 - Chuẩn Euclide của ma trận A (chuẩn Frobenius) ∑∑ a 2 A F = ij = trace( AT A) (3.50) i j với A* là ma trận chuyển vị và lấy liên hiệp. λi ( A* A) là trị riêng của ma trận A* A là một số thực không âm. 3.1.4 Trị suy biến của ma trận – độ lợi chính(Principal gain) Trị suy biến của ma trận A(m x l) được ký hiệu là σ i ( A) được định nghĩa như sau: σ i ( A) = λi ( A* A) i = 1,2,...k (3.51) với k = min{m, l} . Nếu chúng ta biểu diễn ma trận A dưới dạng A(s) và đặt s = jω (0 ≤ ω < ∞) , thì trị suy biến của A( jω ) là một hàm của ω và được gọi là độ lợi chính của A(s). Ở đây chúng ta giả sử rằng σ i được sắp xếp theo thứ tự sao cho σ i ≥ σ i +1 . Như vậy, σ 1 là trị suy biến lớn nhất và σ k là trị suy biến nhỏ nhất. Ký hiệu σ là trị suy biến lớn nhất và σ là trị suy biến nhỏ nhất. Ta có: σ ( A) = max σ i ( A) = max λi ( A* A) = A 2 (3.52) Ax với A 2 = sup 2 . x 2 Độ lợi của hệ đa biến nằm giữa độ lợi chính lớn nhất và nhỏ nhất. Trong Matlab tìm trị suy biến của ma trận A dùng lệnh svd(A) Ví dụ: Cho ma trận A: http://www.khvt.com
  11. Chương 3 : Điều khiển bền vững ⎡9 4⎤ >> A = ⎢⎢6 8 ⎥⎥ ; ⎢⎣2 7 ⎥⎦ >> S =svd(A) S= [14.9359 5.1883] S: vector của các giá trị suy biến của ma trận A σ ( A) =14.9359 σ (A)=5.1883 3.1.5 Ổn định nội Ổn định nội là yêu cầu cơ bản đối với một hệ thống hồi tiếp thực. Ý nghĩa của ổn định nội là khi đầu vào hệ thống bằng không thì tất cả các trạng thái hệ thống đều phải về không từ mọi giá trị ban đầu. Mọi hệ thống tự động đều phải bảo đảm ổn định nội mới hoạt động được. + e1 w1 G + + K + w2 e2 Hình 3.2 : Sơ đồ hệ thống dùng để phân tích ổn định nội Định nghĩa : Hệ hồi tiếp hình 3.2 được gọi là ổn định nội nếu tất cả các hàm truyền đạt từ w1, w2 đến e1, e2 đều ổn định. Trang 11
  12. PGS.TS Nguyễn Thị Phương Hà Trang 12 Điều kiện ổn định nội chặt hơn điều kiện ổn định dựa trên hàm truyền vào- ra thông thường, vì nó tránh việc khử các cực và zero không ổn định giữa các khâu liên tiếp nhau. Khi thành lập hàm truyền vào-ra, có thể xảy ra hiện tượng khử cực và zero không ổn định của các khâu liên tiếp nhau. Như vậy, điều kiện ổn định nội bảo đảm các tín hiệu bên trong hệ thống đều hữu hạn khi tín hiệu vào là hữu hạn. Ví dụ, ta khảo sát điều kiện ổn định nội của hệ thống hình 3.2: e1 = w1 + Ke 2 = w1 + Kw2 + KGe1 ⇒ e1 = ( I − KG ) −1 w1 + ( I − KG ) −1 Kw2 e2 = w2 + Ge1 = w2 + Gw1 + GKe 2 ⇒ e2 = ( I − GK ) −1 Gw1 + ( I − GK ) −1 w2 Suy ra: ⎡ e1 ⎤ ⎡ ( I − KG ) −1 ( I − KG ) −1 K ⎤ ⎡ w1 ⎤ ⎢e ⎥ = ⎢ −1 −1 ⎥ ⎢ ⎥ ⎣ 2 ⎦ ⎣( I − GK ) G ( I − GK ) ⎦ ⎣ w2 ⎦ Điều kiện ổn định nội của hệ là các hàm truyền ( I − KG ) −1 , ( I − KG ) −1 K , ( I − GK ) −1 G , ( I − GK ) −1 đều ổn định. 3.1.6 Định lý độ lợi nhỏ (Small Gain Theorem) Cho hệ thống được biểu diễn như hình 3.3: Gọi λi là trị riêng của G r u y G - Hình 3.3 : Hệ thống hồi tiếp vòng kín Định lý độ lợi nhỏ được phát biểu như sau: Giả thiết rằng G(s) ổn định, ρ(G(jω)) là bán kính phổ của G(jω). Hệ thống vòng kín ổn định nếu ρ (G ( jω )) = max( λi ) < 1 , hoặc G ( jω ) < 1, ∀ω http://www.khvt.com
  13. Chương 3 : Điều khiển bền vững Đối với hệ SISO thì ρ (G ( jω )) = G ( jω ) < 1 (3.53) Định lý độ lợi nhỏ chỉ là điều kiện đủ để xét ổn định của hệ thống. Điểm mạnh của định lí này là nó không yêu cầu những thông tin chi tiết về hệ thống.Vì vậy nó không chỉ ứng dụng được cho hệ thống tuyến tính bất biến theo thời gian mà còn ứng dụng được cho hệ thống phi tuyến, thay đổi theo thời gian. 3.1.7 Ổn định bền vững 3.1.7.1 Định lý ổn định bền vững Đây là mô hình cơ bản dùng để phân tích tính ổn định bền vững của một hệ thống. Nếu hệ danh định ổn định thì M ổn định và Δ là sai số có thể làm cho hệ thống mất ổn định. Định lý sau thiết lập điều kiện của M để cho hệ thống vẫn ổn định dưới ảnh hưởng của Δ v w Δ M Hình 3.4 : Sơ đồ cấu trúc phân tích ổn định bền vững Định lý ổn định bền vững: Giả sử M và Δ ổn định, hệ thống vòng kín hình 3.4 sẽ ổn định khi và chỉ khi biểu đồ cực của đường cong Nyquist det(I-MΔ) không bao điểm gốc. Khi đó hệ thống vòng kín sẽ ổn định bền vững với mọi Δ (σ (Δ) ≤ 1) nếu và chỉ nếu khi một trong các điều kiện sau thỏa mãn: a. Det ( I − MΔ ( jω )) ≠ 0 ∀ω , ∀Δ (σ ≤ 1) (3.54) b. ρ ( MΔ( jω )) < 1 ∀ω , ∀Δ (σ ≤ 1) (3.55) c. M ∞ = σ ( M ( jω )) < 1 ∀ω (3.56) Trang 13
  14. PGS.TS Nguyễn Thị Phương Hà Trang 14 3.1.7.2 Điều kiện ổn định bền vững đối với sai số cộng: Với Δ A ( s ) = δ A ( s )Δ( s ), σ (Δ( jω )) ≤ 1 ∀ω , (3.57) v w Δ M δA + K G - Hình 3.5 : Sai số cộng Ta có: v( s ) = − K ( s )[δ A ( s ) w( s ) + G ( s )v( s )] (3.58) hay v( s ) = −[ I + K ( s )G ( s )]−1 K ( s )δ A ( s ) w( s ) (3.59) vậy K ( s )δ A ( s ) M ( s) = − (3.60) [ I + K ( s )G ( s )] Kết luận: Hệ thống vòng kín hình 3.5 ổn định bền vững khi và chỉ khi: K ( s )δ A ( s ) σ ( jω ) =||M(s)||∞=
  15. Chương 3 : Điều khiển bền vững 3.1.7.3 Điều kiện ổn định bền vững đối với sai số nhân ở đầu ra v w Δ M δ0 + K G - Hình 3.6 : Sai số nhân ở đầu ra Với Δ O ( s ) = δ O ( s )Δ( s ), σ (Δ( jω )) ≤ 1 ∀ω , (3.62) Ta có: v( s) = −G ( s) K ( s)[δ O ( s) w( s ) + v( s )] (3.63) hay v( s ) = −[ I + G ( s ) K ( s )]−1 G ( s ) K ( s )δ O ( s ) w( s ) (3.64) vậy G ( s ) K ( s )δ O ( s ) M =− (3.65) I + G ( s) K ( s) Kết luận: Hệ thống vòng kín hình 3.6 ổn định bền vững khi và chỉ khi: G ( s ) K ( s )δ O ( s )
  16. PGS.TS Nguyễn Thị Phương Hà Trang 16 3.2 Phương pháp LQG (Linear Quadratic Gaussian) 3.2.1 Đặt vấn đề Cho hệ thống x (t ) = Ax(t ) + Bu (t ) + γw(t ) y (t ) = Cx (t ) + v(t ) (3.67) z (t ) = Dx(t ) t∈R Ngõ ra y là ngõ ra hồi tiếp và đo được. Ngõ ra z là điều khiển được. Tín hiệu nhiễu w là nhiễu hệ thống và v là nhiễu đo . Tín hiệu v và w là những quá trình nhiễu trắng .Trạng thái ban đầu của x(0) được giả sử là một vector ngẫu nhiên . Nhiều sự giả sử khác nhau định nghĩa trạng thái x(t) t ∈ R và ngõ ra điều khiển được z(t),t ∈ R là những quá trình ngẫu nhiên .Biểu thức sai số toàn phương : z T (t )Qz (t ) + u T (t ) Ru (t ) t≥0 (3.68) là một quá trình ngẫu nhiên. Vấn đề của điều khiển hệ thống là giá trị mong đợi của tích phân : T ∫ E [z (t )Qz (t ) + u T (t ) Ru (t )]dt T (3.69) 0 là nhỏ. Đây là vấn đề điều khiển tuyến tính nhiễu loạn. Khoảng thời gian [0 T] là xác định nhưng thật sự chúng ta xem xét trường hợp T → ∞ . Tại bất kỳ thời gian t toàn bộ tín hiệu đo được ở quá khứ y(s) s ≤ t được giả sử có giá trị cho hồi tiếp. Hình (3.7) làm rõ trường hợp này : http://www.khvt.com
  17. Chương 3 : Điều khiển bền vững w z u SYSTEM + + y v CONTROLLER Hình 3.7 : Hồi tiếp LQG 3.2.2 Bộ quan sát Xem xét hệ thống quan sát : x (t ) = Ax(t ) + Bu (t ) y (t ) = Cx (t ) (3.70) t∈R Đây là hệ thống (3.67) nhưng không có nhiễu hệ thống w và nhiễu đo v. Trạng thái x của hệ thống (3.70) không thể sử dụng được trực tiếp bởi vì chỉ ngõ ra y là đo được. Xây dựng lại trạng thái với sự chính xác tùy ý bởi việc kết nối một bộ quan sát : xˆ (t ) = Axˆ (t ) + Bu (t ) + L[ y (t ) − Cxˆ (t )] t ∈ R (3.71) Tín hiệu xˆ (t ) là một ước lượng của trạng thái x(t).Nó thỏa mãn phương trình vi phân trạng thái của hệ thống (3.70) với thành phần thêm vào L [ y (t ) − Cxˆ (t )] .L là ma trận độ lợi quan sát cần được lựa chọn phù hợp. Sai số quan sát y(t) − Cxˆ (t ) là sự khác nhau giữa ngõ ra đo được thực tế y(t) và ngõ ra yˆ (t ) = Cxˆ (t ) .Thành phần thêm vào L [ y (t ) − Cxˆ (t )] cung cấp một sự điều chỉnh chủ động ngay khi sai số của sự quan sát là khác 0. Trang 17
  18. PGS.TS Nguyễn Thị Phương Hà Trang 18 u z SYSTEM y + L - yˆ SYSTEM MODEL xˆ Hình 3.8 : Cấu trúc của một bộ quan sát Hình (3.8) cho thấy cấu trúc của bộ quan sát .Định nghĩa : ~ x (t ) = x(t ) − xˆ (t ) (3.72) là sai số ước lượng trạng thái. Phương trình vi phân của ~ x nhận được sau khi trừ (3.70) cho (3.71) : ~ x (t ) = ( A − LC ) ~ x (t ) t ∈ R (3.73) Nếu hệ thống (3.70) được tìm thấy thì tồn tại ma trận độ lợi L mà sai số hệ thống (3.73) là ổn định. Nếu sai số hệ thống là ổn định thì ~ x (t ) → 0 khi t → ∞ cho bất kỳ sai số ~ x (0). Vì vậy →∞ xˆ (t ) ⎯t⎯ ⎯→ x(t ) (3.74) Trạng thái ước lượng hội tụ về trạng thái thực. Trong Matlab dùng hai lệnh acker và place để tính ma trận L của khâu quan sát trạng thái : L= acker(A’,C’,p) L= place(A’,C’,p) A’ : Chuyển vị của ma trận A C’ : Chuyển vị của ma trận C p : Khai báo các điểm cực mong muốn http://www.khvt.com
  19. Chương 3 : Điều khiển bền vững 3.2.3 Bộ lọc Kalman 3.2.3.1 Đặt vấn đề: Bộ lọc Kalman là một bộ quan sát được sử dụng cho các ứng dụng yêu cầu xây dựng lại hệ phương trình trạng thái khi tính đến ảnh hưởng của nhiễu đo được. Phương trình trạng thái của đối tượng : x =Ax+Bu+ γ w (3.75) y=Cx+v (3.76) với trạng thái x(t) ∈ R n ,ngõ vào điều khiển u(t) ∈ R m , và ngõ ra đo lường y(t) ∈ R p .Tín hiệu w(t) là nhiễu quá trình chưa biết trước tác động làm nhiễu hệ thống.Tín hiệu v(t) là một nhiễu đo không xác định được , làm suy giảm việc đo lường chẳng hạn như nhiễu cảm biến.Giá trị ban đầu x(0), nhiễu w(t) hoặc v(t) không biết được chính xác.Giả sử x(0), w(t) và v(t) đều trực giao qua lại với nhau. w(t) v(t) γ u x x y B ∫ C A Hệ thống ~ y L - xˆ xˆ yˆ B ∫ C A Bộ lọc Kalman Hình 3.9 : Bộ quan sát trạng thái của Kalman Trang 19
  20. PGS.TS Nguyễn Thị Phương Hà Trang 20 Gọi xˆ (t ) là ước lượng của x . Phương trình trạng thái của khâu lọc Kalman : xˆ = Axˆ + Bu + L( y − yˆ ) (3.77) yˆ = Cxˆ Mục tiêu của thiết kế bộ lọc Kalman : Tìm độ lợi ước lượng L để có sự ước lượng tối ưu trong sự hiện diện của nhiễu w(t) và v(t) Sai số ước lượng: ~x (t ) = x(t ) − xˆ (t ) (3.78) Độ lợi L sẽ được chọn sao cho giá trị trung bình của sai số ước lượng toàn phương là bé nhất . 3.2.3.2 Cơ sở toán học: Lý thuyết xác suất: Từ phương trình (3.75) được thêm vào bởi nhiễu quá trình, trạng thái x(t) bây giờ cũng là một quá trình ngẫu nhiên như là y(t). Để khảo sát những đặc tính thông thường của quá trình ngẫu nhiên cần nhắc lại một số khái niệm lý thuyết xác suất (Papoulis 1984). Mặc dù w(t) và v(t) là những đại lượng ngẫu nhiên không biết được, nhưng cần biết một vài đặc điểm để hổ trợ việc thiết kế các bộ điều khiển. Chẳng hạn như có thể biết được giá trị trung bình hoặc tổng năng lượng của chúng. Cho vector ngẫu nhiên z ∈ R n ,f z ( ξ ) là hàm mật độ xác suất (PDF) của z. Đại lượng PDF đặc trưng cho xác suất mà z lấy giá trị bên trong vùng vi phân d ξ đặt giữa ξ . Giá trị mong muốn của hàm g(z) của vector ngẫu nhiên được xác định như sau : ∞ E {g ( z )} = ∫ g (ξ ) ƒ( ξ )d ξ (3.79) −∞ Giá trị trung bình hay mong muốn của z được xác định như sau: ∞ E{z}= ∫ ξ ƒ z ( ξ )d ξ (3.80) −∞ được ký hiệu bằng z . Chú ý rằng z ∈ R n . http://www.khvt.com

CÓ THỂ BẠN MUỐN DOWNLOAD

 

Đồng bộ tài khoản