24小时热门版块排行榜    

查看: 588  |  回复: 4

miaozaoyang

木虫 (正式写手)

么么哒


[求助] 如何求解这样的一个Riccati方程?交流一下

待解方程形式为:
S为待求解的,A、C都已知
关键是中间多了一项S,不太好搞定。

如何求解呢?
are()可求解的是:A'*X + X*A - X*B*X +C = 0 这样的方程;
还有没有其他的函数啊?
谢谢啦。在线等一下。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinhu1988

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
简单变形一下即可,将求解的方程变为

S(I+A)+(I+A)^TS-S-C=0

然后用are命令求解
君子务本,静水流深
2楼2012-12-14 02:59:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

miaozaoyang

木虫 (正式写手)

么么哒


引用回帖:
2楼: Originally posted by yinhu1988 at 2012-12-14 02:59:51
简单变形一下即可,将求解的方程变为

S(I+A)+(I+A)^TS-S-C=0

然后用are命令求解

这样还是不对啊。不过,谢谢你了。
后来我变形成这样的形式,用lyap()求解的。
S(0.5I+A)+(0.5I+A)^TS-C=0
还是要谢谢你了。
3楼2012-12-14 09:36:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yinhu1988

木虫 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by miaozaoyang at 2012-12-14 02:36:21
这样还是不对啊。不过,谢谢你了。
后来我变形成这样的形式,用lyap()求解的。
S(0.5I+A)+(0.5I+A)^TS-C=0
还是要谢谢你了。...

这是一个很简单的Riccati方程,变形后ARE是可以直接求解的。求解前先检查你的矩阵符合不符合求解标准型。只要你的矩阵满足要求,不会不可能求解的。
君子务本,静水流深
4楼2012-12-14 17:30:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

miaozaoyang

木虫 (正式写手)

么么哒


引用回帖:
4楼: Originally posted by yinhu1988 at 2012-12-14 17:30:27
这是一个很简单的Riccati方程,变形后ARE是可以直接求解的。求解前先检查你的矩阵符合不符合求解标准型。只要你的矩阵满足要求,不会不可能求解的。...

恩,可以的。现在已经解决的。谢谢啦
5楼2012-12-15 19:13:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 miaozaoyang 的主题更新
信息提示
请填处理意见