How can I solve the Riccati equation with a variable?

2 次查看(过去 30 天)
I have the fixed value of A,B,Q and a variable R. How can I determine the solution of the Riccati equation P as a function of r?
A=[0 1 0;0 0 1;0 0 0];
B=[0;0;0];
Q=[1 0 0;0 0 0;0 0 0];
R=r;
eqn= A'*P+P*A+Q-P*B*R^(-1)*B'*P==0

回答(1 个)

Jaynik
Jaynik 2024-4-4
For determining the solution of the Riccati equation, you can use the "care" or "icare" function from the "Control System Toolbox". These functions are used to compute unique solution for the given equation.
As per the documentation, starting in R2019a, it is recommended to use the "icare" function to solve the equation as it has improved accuracy through better scaling. Here is a sample code you can refer:
A = [0 1 0; 0 0 1; 0 0 0];
B = [0; 0; 0];
Q = [1 0 0; 0 0 0; 0 0 0];
R = r;
% Solving the Riccati equation using care
[P1, ~, ~] = care(A, B, Q, R);
% Solving the Riccati equation using icare
[P2,~,~,info] = icare(A, B, Q, R);
You can refer the following documentation to know more about each functions:

类别

Help CenterFile Exchange 中查找有关 Matrix Computations 的更多信息

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by