✨Điều khiển Gauss tuyến tính-bậc hai

Trong lý thuyết điều khiển tự động, bài toán điều khiển Gauss tuyến tính-bậc hai (LQG) là một trong những bài toán điều khiển tối ưu cơ bản nhất. Nó liên quan đến các hệ thống tuyến tính bị nhiễu bởi nhiễu phụ trắng Gaussian, có thông tin trạng thái không hoàn toàn (tức là không phải tất cả các biến trạng thái đo lường được và sẵn sàng phản hồi) và trải qua đối tượng điều khiển với chi phí bậc hai. Hơn nữa, lời giải là duy nhất và tạo thành một luật điều khiển phản hồi động học tuyến tính dễ dàng tính toán và thực hiện. Cuối cùng bộ điều khiển LQG cũng là nền tảng cho điều khiển tối ưu của các hệ thống phi tuyến nhiễu loạn.

Bộ điều khiển LQG chỉ đơn giản là sự kết hợp của một bộ lọc Kalman, tức là hàm ước lượng tuyến tính bậc hai (LQE), với một bộ điều chỉnh tuyến tính bậc hai (LQR). Nguyên tắc tách biệt đảm bảo rằng chúng có thể được thiết kế và tính toán độc lập. Điều khiển LQG áp dụng cho cả các hệ thống tuyến tính thời gian bất biến và các hệ thống tuyến tính thời gian biến đổi. Việc áp dụng vào các hệ thống tuyến tính biến đổi theo thời gian là nổi tiếng. Việc áp dụng đối với các hệ thống tuyến tính thời gian biến đổi cho phép thiết kế các bộ điều khiển phản hồi tuyến tính cho các hệ thống phi tuyến không ổn định.

Chính bộ điều khiển LQG là một hệ thống động học giống như hệ thống mà nó điều khiển. Cả hai hệ thống đều có cùng kích thước trạng thái. Do đó, việc thực thi bộ điều khiển LQG có thể có vấn đề nếu kích thước của trạng thái hệ thống là lớn. Bài toán LQG giảm bậc (bài toán LQG bậc cố định) sẽ khắc phục điều này bằng cách cố định một _tiên nghiệm _số lượng trạng thái của bộ điều khiển LQG. Bài toán này khó giải hơn vì nó không còn tách rời. Ngoài ra lời giải không còn là duy nhất. Mặc dù những thuật toán số học này là có sẵn  để giải các phương trình chiếu tối ưu liên quan tạo thành điều kiện cần và đủ cho một bộ điều khiển LQG giảm bậc tối ưu cục bộ. Sự ổn định mạnh mẽ của hệ thống vòng kín phải được kiểm tra riêng biệt sau khi bộ điều khiển LQG được thiết kế. Để thúc đẩy mạnh mẽ một số các thông số hệ thống có thể được giả thiết ngẫu nhiên thay vì xác định. Bài toán điều khiển liên quan khó hơn dẫn đến một bộ điều khiển tối ưu tương tự trong đó chỉ có các thông số điều khiển là khác nhau.

Mô tả toán học của bài toán và lời giải

Thời gian liên tục

Xem xét hệ thống động học tuyến tính thời gian liên tục : $\backslash dot\{\backslash mathbf\{x(t)\; =\; A(t)\; \backslash mathbf\{x\}(t)\; +\; B(t)\; \backslash mathbf\{u\}(t)\; +\; \backslash mathbf\{v\}(t),$ : $\backslash mathbf\{y\}(t)\; =\; C(t)\; \backslash mathbf\{x\}(t)\; +\; \backslash mathbf\{w\}(t),$ trong đó $\{\backslash mathbf\{x$ đại diện cho các biến trạng thái của hệ thống, $\{\backslash mathbf\{u$ là vector đầu vào điều khiển và $\{\backslash mathbf\{y$ là vector đầu ra được đo lường dùng để phản hồi. Cả nhiễu hệ thống Gauss phụ trắng $\backslash mathbf\{v\}(t)$ và nhiễu đo lường Gauss phụ trắng đều tác động tới hệ thống. Với hệ thống này, mục tiêu là phải tìm thấy lịch sử đầu vào điều khiển $\{\backslash mathbf\{u(t)$ tại mọi lúc $\{\backslash mathbf\{t$ có thể chỉ phụ thuộc vào các phép đo quá khứ  như vậy mà hàm chi phí sau đây được cực tiểu hóa: : $J\; =\; E\backslash left[\{\backslash mathbf\{x\}^\backslash mathrm\; T\}(T)F\{\backslash mathbf\{x(T)+\; \backslash int\_\{0\}^\{T\}\; \{\backslash mathbf\{x\}^\backslash mathrm\; T\}(t)Q(t)\{\backslash mathbf\{x(t)\; +\; \{\backslash mathbf\{u\}^\backslash mathrm\; T\}(t)R(t)\{\backslash mathbf\{u(t)\backslash ,dt\; \backslash right],$

: $F\; \backslash ge\; 0,\backslash quad\; Q(t)\; \backslash ge\; 0,\backslash quad\; R(t)\; >\; 0,$ trong đó $\{\backslash mathbf\{E$ là ký hiệu của giá trị mong muốn. Thời gian cuối cùng (chân trời) $\{\backslash mathbf\{T$ có thể là hữu hạn hoặc vô hạn. Nếu đường chân trời có xu hướng tiến đến vô cùng số hạng đầu tiên $\{\backslash mathbf\{x^\backslash mathrm\; T(T)F\{\backslash mathbf\{x(T)$ của hàm chi phí trở nên không đáng kể và không liên quan đến bài toán. Ngoài ra để giữ cho chi phí hữu hạn, hàm chi phí phải được thực hiện để có $\{\backslash mathbf\{J/T$.

Bộ điều khiển LQG mà giải quyết được bài toán điều khiển LQG được quy định bởi các phương trình sau đây: : $\backslash dot\{\backslash hat\{\backslash mathbf\{x\}(t)\; =\; A(t)\backslash hat\{\backslash mathbf\{x(t)\; +\; B(t)\{\backslash mathbf\{u(t)+K(t)\; \backslash left(\{\backslash mathbf\{y(t)-C(t)\backslash hat\{\backslash mathbf\{x(t)\; \backslash right),\; \backslash hat\{\backslash mathbf\{x(0)\; =\; E\; \backslash left[\; \{\backslash mathbf\{x(0)\; \backslash right],$

: $\{\backslash mathbf\{u(t)=\; -L(t)\; \backslash hat\{\backslash mathbf\{x(t).$ Ma trận $\{\backslash mathbf\{K(t)$ được gọi là độ lợi Kalman của bộ lọc Kalman liên quan được thể hiện bởi phương trình đầu tiên. Tại mỗi thời điểm $\{\backslash mathbf\{t$ bộ lọc này tạo ra ước tính $\backslash hat\{\backslash mathbf\{x(t)$ của trạng thái $\{\backslash mathbf\{x(t)$ sử dụng các phép đo và đầu vào trong quá khứ. Độ lợi Kalman $\{\backslash mathbf\{K(t)$ được tính toán từ các ma trận $\{\backslash mathbf\{A(t),\; C(t)$, hai ma trận cường độ $\backslash mathbf\{\}V(t),\; W(t)$ liên quan đến các nhiếu Gauss trắng $\backslash mathbf\{v\}(t)$ và $\backslash mathbf\{w\}(t)$ và cuối cùng $E\backslash left[\{\backslash mathbf\{x(0)\{\backslash mathbf\{x^\backslash mathrm\; T(0)\; \backslash right]$. Năm ma trận này xác định độ lợi Kalman thông qua ma trận phương trình vi phân Riccati liên quan sau đây: : $\backslash dot\{P\}(t)\; =\; A(t)P(t)+P(t)A^\backslash mathrm\; T(t)-P(t)C^\backslash mathrm\; T(t)\{\backslash mathbf\{W^\{-1\}(t)\; C(t)P(t)+V(t),$

: $P(0)=\; E\; \backslash left[\{\backslash mathbf\{x(0)\{\backslash mathbf\{x^\backslash mathrm\; T(0)\; \backslash right].$ Với các giải pháp $P(t),\; 0\; \backslash leq\; t\; \backslash leq\; T$ độ lợi Kalman bằng : $\{\backslash mathbf\{K(t)\; =\; P(t)C^\backslash mathrm\; T(t)W^\{-1\}(t).$ Ma trận $\{\backslash mathbf\{L(t)$ được gọi là ma trận độ lợi phản hồi. Ma trận này được xác định bởi các ma trận $\{\backslash mathbf\{A(t),\; B(t),\; Q(t),\; R(t)$và $\{\backslash mathbf\{F$ thông qua ma trận phương trình vi phân Riccati liên quan sau đây: : $\{\backslash mathbf\{S(T)\; =\; F.$ Với các lời giải $\{\backslash mathbf\{S(t),\; 0\; \backslash leq\; t\; \backslash leq\; T$ độ lợi phản hồi bằng : $\{\backslash mathbf\{L(t)\; =\; R^\{-1\}(t)B^\backslash mathrm\; T(t)S(t).$ Quan sát sự giống nhau của hai ma trận phương trình vi phân Riccati, phương trình đầu tiên chạy về phía trước theo thời gian, phương trình thứ hai chạy về phía ngược theo thời gian. sự tương đồng này được gọi là nhị nguyên. Ma trận phương trình vi phân Riccati đầu tiên giải quyết bài toán ước lượng tuyến tính bậc hai (LQE). Ma trận phương trình vi phân Riccati thứ hai giải quyết bài toán bộ điều chỉnh tuyến tính bậc hai (LQR). Những bài toán này là sóng đôi và chúng cùng nhau giải quyết bài toán tuyến tính bậc hai-Gaussian (LQG). Vì vậy, bài toán LQG chia tách thành bài toán LQE và LQR có thể được giải quyết một cách độc lập. Do đó, bài toán LQG được gọi là có thể tách rời.

Khi $\{\backslash mathbf\{A(t),\; B(t),\; C(t),\; Q(t),\; R(t)$ và các ma trận cường độ nhiễu$\backslash mathbf\{\}V(t)$, $\backslash mathbf\{\}W(t)$ không phụ thuộc vào $\{\backslash mathbf\{t$ và khi $\{\backslash mathbf\{T$ có xu hướng tiến đến vô cùng, bộ điều khiển LQG trở thành một hệ động lực biến đổi theo thời gian. Trong trường hợp đó cả hai ma trận phương trình vi phân Riccati có thể được thay thế bằng hai phương trình Riccati đại số liên quan.

Thời gian rời rạc

Do bài toán điều khiển LQG thời gian rời rạc tương tự với bài toán trong thời gian liên tục, mô tả dưới đây tập trung vào các phương trình toán học.

Các phương trình hệ thống tuyến tính thời gian rời rạc là : $\{\backslash mathbf\{x\_\{i+1\}\; =\; A\_i\backslash mathbf\{x\}\_i\; +\; B\_i\; \backslash mathbf\{u\}\_i\; +\; \backslash mathbf\{v\}\_i,$

: $\backslash mathbf\{y\}${i} = C{i} \mathbf{x}_i + \mathbf{w}i. Trong đó $\backslash mathbf\{\}i$ đại diện cho các chỉ số thời gian rời rạc và $\backslash mathbf\{v\}${i}, \mathbf{w}{i} đại diện cho nhiễu trắng Gauss thời gian rời rạc xử lý với các ma trận hiệp phương sai $\backslash mathbf\{\}V${i}, W_{i} tương ứng.

Hàm chi phí bậc hai được cực tiểu hóa là : $J\; =\; E\backslash left[\{\backslash mathbf\{x^\backslash mathrm\; T${N}F{\mathbf{x{N}+ \sum_{i=0}^{N-1}(\mathbf{x}_i^\mathrm T Q_i \mathbf{x}_i + \mathbf{u}_i^\mathrm T R_i \mathbf{u}_i)\right],

: $F\; \backslash ge\; 0,\; Q\_i\; \backslash ge\; 0,\; R\_i\; >\; 0.\; \backslash ,$ Bộ điều khiển LQG thời gian rời rạc là : $\backslash mathbf\{u\}\_i=-L\_i\backslash hat\{\backslash mathbf\{x\_i.\; \backslash ,$ Độ lợi Kalman bằng : $\{\backslash mathbf\{K\_i\; =\; P\_iC\; ^\backslash mathrm\; T\; \_i(C\_iP\_iC\; ^\backslash mathrm\; T\; \_i\; +\; W\_i)^\{-1\},$ trong đó $\{\backslash mathbf\{P\_i$ được xác định bởi ma trận phương trình vi phân Riccati sau đây mà chạy về phía trước theo thời gian:

Ma trận độ lợi phản hồi bằng : $\{\backslash mathbf\{L\_i\; =\; (B^\backslash mathrm\; T$iS{i+1}B_i + R_i)^{-1}B^\mathrm TiS{i+1}A_i trong đó $\{\backslash mathbf\{S\_i$ được xác định bởi ma trận phương trình vi phân Riccati sau đây theo hướng ngược thời gian: : $S\_i\; =\; A^\backslash mathrm\; T$i \left(S{i+1} - S_{i+1}B_i \left(B^\mathrm TiS{i+1}B_i+R_i \right)^{-1} B^\mathrm Ti S{i+1} \right) A_i+Q_i, \quad S_N=F. Nếu tất cả các ma trận trong việc xây dựng bài toán là thời gian bất biến và nếu đường chân trời $\{\backslash mathbf\{N$ có xu hướng tiến đến vô cùng, bộ điều khiển thời gian rời rạc LQG sẽ trở thành thời gian bất biến. Trong trường hợp đó, ma trận phương trình vi phân Riccati có thể được thay thế bằng phương trình đại số Riccati thời gian rời rạc liên quan của chúng. Những bộ ước lượng tuyến tính bậc hai bất biến theo thời gian và bộ điều chỉnh tuyến tính bậc hai thời gian bất biến được xác định trong thời gian rời rạc. Để giữ cho các chi phí hữu hạn thay vì $\{\backslash mathbf\{J$ ta phải xem xét $\{\backslash mathbf\{J/N$ trong trường hợp này.