主要内容

关心

(不推荐)连续时间代数Riccati方程解决方案

关心不推荐。用我在乎反而。有关更多信息,请参阅兼容性的考虑

句法

[x,l,g] =护理(a,b,q)
[x,l,g] =护理(a,b,q,r,s,e)
[X, L, G,报告]=保健(A, B, Q,…)
(X1, X2, D, L) =保健(A, B, Q,…,“因素”)

描述

[x,l,g] =护理(a,b,q)计算唯一的解决方案X连续时间代数Riccati方程

一种 T. X + X 一种 - X B. B. T. X + 问: = 0.

关心函数还返回增益矩阵, G = R. - 1 B. T. X E.

[x,l,g] =护理(a,b,q,r,s,e)解决了更通用的Riccati方程

一种 T. X E. + E. T. X 一种 - E. T. X B. + S. R. - 1 B. T. X E. + S. T. + 问: = 0.

省略,R.S.,E.是否设置为默认值r = I.s = 0.,E = I.。以及解决方案X关心返回增益矩阵 G = R. - 1 B. T. X E. + S. T. 和一个矢量L.的闭环特征值,其中

L = EIG(A-B * G,E)

[X, L, G,报告]=保健(A, B, Q,…)返回一个诊断报告和:

  • -1当关联的时候Hamiltonian铅笔在虚构的轴上有特征值(失败)

  • -2当没有有限稳定解时X

  • 相对残差的Frobenius规范X存在并且是有限的。

当x无法存在时,此语法不会发出任何错误消息。

(X1, X2, D, L) =保健(A, B, Q,…,“因素”)返回两个矩阵X1X2和对角缩放矩阵D.这样x = D *(x2 / x1)* d

矢量L包含闭环特征值。当关联的Hamiltonian矩阵有虚构轴上的特征值时,所有输出都是空的。

例子

例1

解决代数Riccati方程

给予

一种 = [ - 3. 2 1 1 ] B. = [ 0. 1 ] C = [ 1 - 1 ] R. = 3.

你可以解出Riccati方程

一种 T. X + X 一种 - X B. R. - 1 B. T. X + C T. C = 0.

经过

a = [-3 2; 1 1] b = [0;1] c = [1 -1] r = 3 [x,l,g] =护理(a,b,c'* c,r)

这就产生了解决方案

x x = 0.5895 1.8216 1.8216 8.8188

您可以通过比较特征值来验证此解决方案是否确实稳定一种A-B * g

[Eig(a)eig(a-b * g)] ans = -3.4495 -3.5026 1.4495 -1.4370

最后,请注意变量L.包含闭环特征值eig (a - b * g)

L = -3.5026 -1.4370

例2.

解决H-Infinity( H )-like Riccati方程

来解决 H 例如Riccati方程式

一种 T. X + X 一种 + X γ - 2 B. 1 B. 1 T. - B. 2 B. 2 T. X + C T. C = 0.

重写它关心格式为

一种 T. X + X 一种 - X [ B. 1 B. 2 ] } B. [ - γ 2 一世 0. 0. 一世 ] } R. - 1 [ B. 1 T. B. 2 T. ] X + C T. C = 0.

您现在可以计算稳定解决方案 X 经过

B = [B1,B2] M1 =尺寸(B1,2)M2 =尺寸(B2,2)r = [-g ^ 2 *眼睛(m1)零(m1,m2);零(m2,m1)眼睛(m2)] x =护理(a,b,c'* c,r)

限制

一种 B. 对一定是稳定(即所有不稳定模式都是可控的)。此外,相关的哈密顿矩阵或铅笔在虚轴上必须没有特征值。它成立的充分条件是 问: 一种 可检测的时间 S. = 0. R. > 0. , 或者

[ 问: S. S. T. R. ] > 0.

算法

关心实现算法中描述的算法[1]。当R是良好的时,它适用于Hamiltonian矩阵 E. = 一世 ;否则它使用扩展的汉密尔顿铅笔和QZ算法。

兼容性的考虑

展开全部

不推荐在R2019A开始

参考文献

[1] Arnold,W.f.,III和A.J.替换,“代数Riccati方程的广义eIgenProble算法和软件,”Proc。IEEE.®,72(1984),第746-1754页

也可以看看

之前介绍过的R2006a