I want to solve the Riccati equation symbolically. For example, I have 4 state variables in Q and 2 inputs in R lqr weights.
Q=diag([q11 q22 q33 q44]);
The Riccati equation is
A'*P+P*A-P*B*(inv( R ))*B'*P+Q=0
A:4x4 matrix , B:4x2 matrix
I want to solve that equation with respect to "P" matrix variable.
In other words, the answer should be like P=(...) and that should contain the variables q11 q22 q33 q44 and r11 r22. Is there anyone who can help? I appreciate all answers,