Solve Riccati (DARE) with MATLAB

  • Context: MATLAB 
  • Thread starter Thread starter perplexabot
  • Start date Start date
  • Tags Tags
    Matlab
Click For Summary
SUMMARY

The forum discussion centers on the use of MATLAB's dare() function for solving the discrete algebraic Riccati equation (DARE). The user defines variables such as E = I, B = I, R = \Phi_R, Q = \Phi_{33}, S = \Phi_{32}, and A = 2\Phi_R^{-1}S^T. They seek confirmation on the correctness of their DARE formulation and express confusion regarding the output of their code. The discussion highlights the importance of validating results against expected outcomes.

PREREQUISITES
  • Understanding of discrete algebraic Riccati equations (DARE)
  • Familiarity with MATLAB programming and syntax
  • Knowledge of matrix operations and inverses
  • Basic control theory concepts related to state-space representation
NEXT STEPS
  • Review MATLAB's dare() function documentation for detailed usage
  • Learn about the implications of matrix dimensions in DARE formulations
  • Explore MATLAB tutorials on care() and dare() for continuous and discrete cases
  • Investigate common pitfalls in DARE implementations and their solutions
USEFUL FOR

Control engineers, MATLAB users, and researchers working on discrete-time control systems who need to solve Riccati equations effectively.

perplexabot
Gold Member
Messages
328
Reaction score
5
Hey all, I was wondering if you can help me define the arguments for the dare() function of matlab. Here is the DARE form I have:
<br /> \begin{equation}<br /> \begin{split}<br /> \Phi(p_\infty) &amp;= \Phi_\infty = \Phi_{33}-\Phi_{32}(\Phi_\infty+\Phi_{R})^{-1}\Phi_{23}\\<br /> &amp;=\Phi_{33}+\Phi_{32}\Phi_{R}^{-1}\Phi_\infty\Phi_{R}^{-1}\Phi_{23}-\Phi_{32}\Phi_{R}^{-1}\Phi_\infty(\Phi_\infty+\Phi_{R})^{-1}\Phi_\infty\Phi_{R}^{-1}\Phi_{23}-\Phi_{32}\Phi_{R}^{-1}\Phi_{23}<br /> \end{split}<br /> \end{equation}<br />

Here is how MATLAB's dare() function works.

Here is how I defined my variables: E = I, B = I, R = \Phi_R, Q = \Phi_{33}, S=\Phi_{32}, , A = 2\Phi_R^{-1}S^T.

Is this correct?
Thank you : )

PS: Please let me know if this thread belongs elsewhere, as I would like to maximize feedback from the PF community.
 
Physics news on Phys.org
perplexabot said:
That's the same link I provided in my original post : (

Sorry I didn't see as it was hidden in the "Here" text.
 
jedishrfu said:
Sorry I didn't see as it was hidden in the "Here" text.
No problem! I still appreciate the fact that you are trying to help. Thank you.
 
Have you tried to run your code? If so, what was wrong with the answer?

You could try taking the answer and testing to see if it is the solution to gain confidence in what Matlab returned.
 
jedishrfu said:
Have you tried to run your code? If so, what was wrong with the answer?

You could try taking the answer and testing to see if it is the solution to gain confidence in what Matlab returned.
I have run my code. I got an answer but I am confused with the result as it is not as I expected. I am trying to figure out why the answer is what it is. I was wondering whether my DARE formulation was correct (hence, why I have this post).
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
1K
  • · Replies 2 ·
Replies
2
Views
835
  • · Replies 4 ·
Replies
4
Views
7K
Replies
1
Views
1K
  • · Replies 0 ·
Replies
0
Views
3K
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 7 ·
Replies
7
Views
6K
  • · Replies 5 ·
Replies
5
Views
5K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 9 ·
Replies
9
Views
2K