如何求recatti方程的稳态解

12次观看(过去30天)
大家好!
我想求如下代数方程的解:
-A_Star ' * x * q + x * G *……x是这个方程的解
其中A_Star, A Q和G的定义如下:
N = 50;
a1 = 1;a2 = 1;
τ= 1;
= 0 (N + 1, N + 1);
=诊断接头(- N /τ*的(N + 1, - 1)) +诊断接头(N /τ*的(N, 1), 1);
(1) = a1;
(1, N + 1) = a2;
B = 0 (N + 1);
B (1) = 1;
Q = 0 (N + 1, N + 1); Q (1, - 1) = 1;
R = 1;
W = 0 (N + 1, N + 1);
W =诊断接头(τ)/ N (N + 1) *的);
W (1, - 1) = 1;
G = B * R * B ';
A_Star =发票(W) *的* W;
提前感谢你的帮助,

接受的答案

Torsten
Torsten 2022年4月16日
如果是里卡提方程,用"icare"或"idare"。
7评论

登录发表评论。

更多答案(1)

山姆翟
山姆翟 2022年4月16日
你可以用隐式代数里卡提方程求解器求出x的解:
[X, K, L] = icare(A_Star, [], Q, [], [], [], G)
对于较老版本的MATLAB(以前 R2019a ),然后用这个:
[X, L, G] = care(A_Star, B, Q)
2的评论

登录发表评论。

s manbetx 845


释放

R2019a

社区寻宝

在MATLAB中央找到宝藏,并发现社区可以如何帮助你!

开始狩猎!