ZNN模型求解动态代数Riccati方程

最新推荐文章于 2023-06-27 17:16:28 发布

Bokman

最新推荐文章于 2023-06-27 17:16:28 发布

阅读量976

点赞数

分类专栏：数值算法文章标签：神经网络算法 matlab

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/bokman/article/details/121595933

版权

数值算法专栏收录该内容

16 篇文章 4 订阅

订阅专栏

动态代数Riccati方程的定义为：

$A^T(t)X(t)+X(t)A(t)-X(t)M(t)X(t)+H(t)=0$

其中矩阵 $A(t)\in\mathbb{R}^{n\times n}$ ， $X(t)\in\mathbb{R}^{n\times n}$ ， $M(t)\in\mathbb{R}^{n\times n}$ ， $H(t)\in\mathbb{R}^{n\times n}$ 。且 $M(t)=M^T(t),~H(t)=H^T(t)$ ， $M(t)$ 和 $H(t)$ 都是正定矩阵。

构造ZNN模型：

定义针对动态代数Riccati问题的误差方程：

$E(t)=A^T(t)X(t)+X(t)A(t)-X(t)M(t)X(t)+H(t)$

ZNN模型的演化公式：

$\dot e(t)=-\gamma e(t)$

结合以上两个等式可以得到求解动态代数Riccati方程的ZNN模型：

$A^T(t) \dot{X}(t)+ \dot{X}(t)A(t)- \dot{X}(t)M(t)X(t)-X(t)M(t) \dot{X}(t)\\ ~~~~~=X(t)\dot M(t)X(t)-\dot {A}^T(t)X(t)-X(t)\dot{A}(t)-\dot {H}(t)\\ ~~~~~-\gamma(A^T(t)X(t)+X(t)A(t)-X(t)M(t)X(t)+H(t))$

对上式进行化简，令：

$B(t)=I\otimes(A^T(t)-X(t)M(t))\in\mathbb{R}^{n^2\times n^2}$

$C(t)=(A(t)-M(t)X(t))^T\otimes I\in\mathbb{R}^{n^2\times n^2}$

$D(t)=I\otimes(X(t)\dot M(t)-\dot A^T(t))\in\mathbb{R}^{n^2\times n^2}$

$G(t)=\dot A^T(t)\otimes \in\mathbb{R}^{n^2\times n^2}$

$J(t)=A^T(t)\otimes I\in\mathbb{R}^{n^2\times n^2}$

$x(t)=\text{vec}X(t)\in\mathbb{R}^{n^2}$

$k(t)=\text{vec}(H(t))\in\mathbb{R}^{n^2}$

于是，以上复杂的等式即可转化为如下简单形式：

$(B(t)+C(t))\dot x(t)=(D(t)-G(t))x(t)-k(t)-\gamma((B(t)+J(t))x(t) + k(t))$

关注

0
点赞
踩
3

收藏

觉得还不错? 一键收藏
0
评论
ZNN模型求解动态代数Riccati方程

动态代数Riccati方程的定义为：其中矩阵，，，。且，和都是正定矩阵。构造ZNN模型：定义针对动态代数Riccati问题的误差方程：ZNN模型的演化公式：结合以上两个等式可以得到求解动态代数Riccati方程的ZNN模型：对上式进行化简，令：于是，以上复杂的等式即可转化为如下简单形式：...
复制链接

扫一扫

专栏目录

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。