function dXdt = mRiccati(t,X,A,B,R,Q)
X = reshape(X, size(A)); %Convert from "n^2"-by-1 to "n"-by-"n"
dXdt = A.'*X + X.*A - X*B.*inv(R).*B'.*X + Q; %Determine derivative
dXdt = dXdt(:); %Convert from "n"-by-"n" to "n^2"-by-1

Waleed new 2018-3-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549643

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549643

在 MATLAB Online 中打开

my programm is :

% diff equat
R=0.1*eye(3,3);
Q=eye(3,3);
B=1/2.*[0,0.3604,-0.4397;-0.3604,0,0.2226;0.4397,-0.2226,0]+0.8224*eye(3,3);
X0 = [1; 1; 1; 1;1;1];
[T X] = ode45(@(t,X)mRiccati(t,X,A,B,R,Q), [0 10], X0)

Waleed new 2018-3-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549644

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549644

it returns this message error:

Torsten 2018-3-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549652

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549652

编辑：Torsten 2018-3-26

I don't see that you define A somewhere.

And X0 must by 9x1, not 6x1.

Waleed new 2018-3-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549669

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549669

在 MATLAB Online 中打开

yes it works :

% diff equat
A=[0,0,0;0,0,0;0,0,0];
R=0.1*eye(3,3);
Q=eye(3,3);
B=1/2.*[0,0.3604,-0.4397;-0.3604,0,0.2226;0.4397,-0.2226,0]+0.8224*eye(3,3);
X0 = [1; 1; 1; 1;1;1;1;1;1];
[T X] = ode45(@(t,X)mRiccati(t,X,A,B,R,Q), [0 10], X0)
the result :

Waleed new 2018-3-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549672

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549672

thank you, how do you know that x0 is 9x1 ?

Torsten 2018-3-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549695

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549695

A is 3x3, and X has the same dimension as A.

ODE45 expects X to be in vector form, thus of length 3x3=9.

Waleed new 2018-3-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549783

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549783

thank you Torsten it's very helpfull

Torsten 2018-3-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549790

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549790

编辑：Torsten 2018-3-26

在 MATLAB Online 中打开

There are quite a few errors in the line where you calculate dXdt. This is the correct formula:

dXdt = A.'*X + X*A - X*B*inv(R)*B.'*X + Q; %Determine derivative

Best wishes

Torsten.

Waleed new 2018-3-26

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549823

⋮

链接

此评论的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/390593-how-can-i-solve-matrix-riccati-differential-equation#comment_549823

i thought x is a matrix so i changed the product () with this (.) i got this figure

</matlabcentral/answers/uploaded_files/110184/2018-03-26_162737.png> but now with the correct formula :

</matlabcentral/answers/uploaded_files/110185/2018-03-26_162259.png> thank you Torsten ,I'm trying to simulate equations that i found it in a research article and all what i do is executing programs and compare the result with figures of the article 8).

请先登录，再进行评论。

请先登录，再回答此问题。