(不推荐)连续时间代数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方程
这关心
函数还返回增益矩阵,
。
[x,l,g] =护理(a,b,q,r,s,e)
解决了更通用的Riccati方程
省略,R.
那S.
, 和E.
设置为默认值r = I.
那s = 0.
, 和E = I.
。以及解决方案X
那关心
返回增益矩阵
和一个矢量L.
闭环特征值,在哪里
L = EIG(A-B * G,E)
[x,l,g,报告] = care(a,b,q,......)
返回诊断报告
和:
当x无法存在时,此语法不会发出任何错误消息。
[x1,x2,d,l] =护理(a,b,q,...,'fing')
返回两个矩阵X1
那X2
和对角缩放矩阵D.
这样x = D *(x2 / x1)* d
。
矢量L包含闭环特征值。当相关的Hamiltonian矩阵有虚构轴上的特征值时,所有输出都是空的。
解决代数Riccati方程
给予
您可以解决Riccati方程式
经过
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
解决H-Infinity( )-like Riccati方程
重写它关心
格式为
您现在可以计算稳定解决方案 经过
B = [B1,B2] M1 =尺寸(B1,2)M2 =尺寸(B2,2)r = [-g ^ 2 *眼睛(m1)零(m1,m2);零(m2,m1)眼睛(m2)] x =护理(a,b,c'* c,r)
关心
实现所描述的算法[1]。当R是良好的时,它适用于Hamiltonian矩阵
;否则它使用扩展的汉密尔顿铅笔和QZ算法。
[1] Arnold,W.f.,III和A.J.洗衣,“代数Riccati方程的广义eIgenProbl算法和软件,”Proc。IEEE.®,72(1984),第1746-1754页