24小时热门版块排行榜    

查看: 2004  |  回复: 6

swim8568

木虫 (小有名气)


[交流] 【求助】怎么用MATLAB求数值解啊(高手请帮忙啊)

求一个2×2的未知矩阵F。
方程为:F=A*(F+(b-F*C')*(inv(E(1)-C*F*C'))*(b-F*C')')*A'
除了F,其它都是已知量,请问,如何求矩阵F?(数值解,即近似解)

请懂的给个方法,给个程序最好了。实在搞不清楚怎么弄!

A = [1  , 0.001  ;   0   ,   1]
T = [-53 ; 2060]
C = [1 , 0]
E=[0.0519]

[ Last edited by swim8568 on 2010-11-12 at 16:32 ]
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

aaron1988

木虫 (正式写手)


swim8568(金币+2): 2010-11-11 09:05:41
数值方法的话, 有JACOBIAN MATRIX 解。
不然,能有ODE的方法来解。
我的建议是,先算出RHS(右手边)的矩阵等于什么先。。
现在没时间帮你解。。。不好意思
2楼2010-11-11 00:40:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swim8568

木虫 (小有名气)


引用回帖:
Originally posted by aaron1988 at 2010-11-11 00:40:09:
数值方法的话, 有JACOBIAN MATRIX 解。
不然,能有ODE的方法来解。
我的建议是,先算出RHS(右手边)的矩阵等于什么先。。
现在没时间帮你解。。。不好意思

要这么麻烦么?有没有其它简单点的方法?
我之前用fsolve,但解得很复杂,而且初值的变化对结果的影响太大
3楼2010-11-11 10:00:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aaron1988

木虫 (正式写手)


swim8568(金币+3): 2010-11-18 14:44:17
看起来麻烦,如果用JACOBIAN MATRIX应该不麻烦的。只要迭代进去计算就好。
如果可以的话,你应该给出A,b, C, 这些数值出来,那么大家可以看看具体怎么解..
4楼2010-11-11 13:09:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

swim8568

木虫 (小有名气)


A = [1  , 0.001  ;   0   ,   1]
b = [-53 ; 2060]
C = [1 , 0]
E(1)=0.0519
5楼2010-11-12 16:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
swim8568(金币+3): 2010-11-18 14:44:21
把未知矩阵拉成一个向量,看看能不能改写下方程。我觉得这么小的维数,弄成一个四元的方程组用solve求解应该也是可行的。fsolve的话还可以加些参数运行,比如自己指定jacobian,或者修改默认算法为线搜索的Gauss Newton,我觉得也值得尝试一下。说明文档里有介绍和例子,可以花点时间读读。
6楼2010-11-16 00:55:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

第一夫人

新虫 (初入文坛)


swim8568(金币+2): 2011-04-02 11:50:51
我来给楼主打气啊!~可惜我自大学开始就丢掉数学了~
7楼2011-03-22 21:46:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 swim8568 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见