主要内容

关心

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

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

句法

[x,l,g] =护理(a,b,q)
[x,l,g] =护理(a,b,q,r,s,e)
[x,l,g,报告] = care(a,b,q,......)
[x1,x2,d,l] =护理(a,b,q,...,'fing')

描述

[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,报告] = care(a,b,q,......)返回诊断报告和:

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

  • -2当没有有限稳定的溶液时X

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

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

[x1,x2,d,l] =护理(a,b,q,...,'fing')返回两个矩阵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 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. 的) 对一定是稳定(即,所有不稳定模式都是可控的)。此外,相关的Hamiltonian矩阵或铅笔必须在假想轴上没有特征值。持有的充分条件是 问: 一种 的) 可检测的时间 S. = 0. R. > 0. , 或者

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

算法

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

兼容性考虑因素

展开全部

在R2019A开始不推荐

参考

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

也可以看看

在R2006A之前介绍