Drexel dragonThe Math ForumDonate to the Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.matlab

Topic: Solving Matrix Riccati
Replies: 3   Last Post: Apr 6, 2013 12:51 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Nasser Abbasi

Posts: 6,064
Registered: 2/7/05
Re: Solving Matrix Riccati
Posted: Apr 5, 2013 7:41 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

On 4/5/2013 6:20 PM, Jeremy wrote:
> I'm new to MATLAB and trying to use ODE45 to find the coefficients of the matrix
>Riccati equation for a controls problem.

Any reason why are you not using the specalized funtions for Riccati solver
in the control systems toolbox?



Continuous-time algebraic Riccati equation solution"

> I'm not sure what I'm doing wrong.
> I've defined my Riccati equation in a function 'regulator'
> %Define matrices
> function dS=regulator(t,S)
> A=[0 1;0 -2.972];
> B=[0;390];
> Q=[(1/(pi^2)) 0;0 30];
> R=1;
> dS=-1*(A'*S+S*A-S*B*(R^-1)*B'*S+Q);
> end
> Then I'm just going straight to ODE45, doing something like this:
> [t,S]=ode45('regulator',[0 5],[0 0;0 0])
> looking for a solution to matrix S.
> Can I solve like this or do I have to multiply out all the terms of S? Assuming S is
> S = [s11(t) s12(t) ; s12(t) s22(t)];
> Any help would be appreciated.

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum 1994-2015. All Rights Reserved.