24小时热门版块排行榜    

查看: 1214  |  回复: 2

wangathlon

新虫 (初入文坛)

[交流] 关于matlab的fsolve函数 已有1人参与

手里有个模型:
  y=c*[(x-a)^2-b^2]/(x-a), 其中a,b,c可以通过实验得出。我取了两组实验值,然后通过fsolve解x,y的数值解,每次exitflag都是返回-2,而且结果y对初始值相当敏感。之后直接用手解两个方程,发现消去y以后,是一个一元三次的方程,应该是有确切的解。不知道为何会产生这样的问题,另外,我原先是想尝试一下lsqcurvefit的,如果有三个解的话,是同时搜索到三个解(如果三个都是实数解)还是程序只搜到第一个解就返回呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bitinging

金虫 (小有名气)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiegangmai: 金币+2, 谢谢参与! 2012-05-23 22:18:36
fsolve求解的方法是用最小二乘法,给定的初始点、方程形式、系数差距对求解都有很大影响。建议可以多试试几个初值然后用fsolve精细控制来控制fsolve求解,或者既然可以化成一元,那可以考虑fzero求解,fzero采用的求解方法比fsolve多很多。
none
2楼2012-05-23 09:15:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bitinging

金虫 (小有名气)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
xiegangmai: 金币+1, 谢谢参与! 2012-05-23 22:18:46
另外,matlab求根算法基本都是只能找到初值附近的解,不会求出多个解的。你可以用三元画图大概看看函数图形,看可能有几个解和解的大致区域。
另外解这类方程据说1st opt比matlab强,对初值要求很小,可以试试。
none
3楼2012-05-23 09:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangathlon 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 15/750 2026-09-01 12:23 by icm639
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 有没有仍没收到信息的 +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
[基金申请] 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
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +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
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
信息提示
请填处理意见