Solve Riccati (DARE) with MATLAB

  • Context: MATLAB 
  • Thread starter Thread starter perplexabot
  • Start date Start date
  • Tags Tags
    Matlab
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
6 replies · 4K views
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:
[tex] \begin{equation}<br /> \begin{split}<br /> \Phi(p_\infty) &= \Phi_\infty = \Phi_{33}-\Phi_{32}(\Phi_\infty+\Phi_{R})^{-1}\Phi_{23}\\<br /> &=\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}[/tex]

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

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

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
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).