24小时热门版块排行榜    

查看: 698  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+3 winsaint 2026-08-27 7/350 2026-08-28 07:06 by hejicker
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 05:13 by winsaint
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见