如何求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;
提前感谢你的帮助,
0评论
更多答案(1)
山姆翟
2022年4月16日
你可以用隐式代数里卡提方程求解器求出x的解:
[X, K, L] = icare(A_Star, [], Q, [], [], [], G)
对于较老版本的MATLAB(以前
R2019a
),然后用这个:
[X, L, G] = care(A_Star, B, Q)