Numerically solving matrix Riccati ODE

  • Context: Graduate 
  • Thread starter Thread starter Avatrin
  • Start date Start date
  • Tags Tags
    Matrix Ode
Click For Summary

Discussion Overview

The discussion revolves around numerically solving a matrix Riccati ordinary differential equation (ODE) of the form $$\dot{X}(t) = FX(t) + X(t)F^T + B$$ with an initial condition. Participants explore various methods and resources for tackling this problem, including numerical techniques and existing literature.

Discussion Character

  • Exploratory, Technical explanation, Homework-related

Main Points Raised

  • One participant expresses uncertainty about how to proceed with solving the matrix Riccati ODE and seeks guidance.
  • Another participant suggests that the Wikipedia page on the algebraic Riccati equation provides a general overview and includes links to MATLAB and Python resources.
  • A participant mentions that while they have reviewed the Wikipedia page, they are unclear on how to apply the algebraic Riccati equation to the differential case.
  • One participant proposes using Euler's method as a potential approach but expresses a desire to understand how the continuous algebraic Riccati equation (CARE) might be applied to their problem.
  • Another participant recommends using MATLAB's Runge-Kutta solver for ODEs and systems of ODEs, while also noting that Euler's method, though simple, may lack stability.

Areas of Agreement / Disagreement

Participants do not reach a consensus on a specific method to solve the matrix Riccati ODE, and multiple approaches are discussed without agreement on the best solution.

Contextual Notes

Participants express varying levels of familiarity with the topic, and there are indications of missing connections between the algebraic and differential forms of the Riccati equation. The discussion reflects uncertainty regarding the stability of different numerical methods.

Who May Find This Useful

Readers interested in numerical methods for solving differential equations, particularly those involving matrix equations, may find this discussion relevant.

Avatrin
Messages
242
Reaction score
6
Hi

I need to solve an equation of the form $$\dot{X}(t) = FX(t) + X(t)F^T + B$$
All of these are matrices. I have an initial condition X(0)=X_0.

However, I have no idea how to proceed. How can I make any progress?
 
Physics news on Phys.org
scottdave said:
The Wikipedia page gives a general overview. https://en.wikipedia.org/wiki/Algebraic_Riccati_equation

But at the bottom there are some links to some MATLAB and Python resources which may help.
Well, I've already skimmed through that one. It is about the algebraic Riccati equation, and does mention that it can be applied to the differential Riccati equation. However, I just don't know how.

So, my best idea at the moment is just to use Euler's method. However, I wish I could find out how CARE can be used to solve my problem.
 
if you have matlab, just try the runge-kutta solver for solving ode's and systems of ode's. You can also use Euler's methods, as it is the simplest numerical method, but it is not always very stable.
 
  • Like
Likes   Reactions: scottdave and Avatrin

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 0 ·
Replies
0
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 24 ·
Replies
24
Views
4K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K