24小时热门版块排行榜    

查看: 1937  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

tomangel

新虫 (初入文坛)

[求助] matlab微分方程组参数拟合,以周为单位求解,汇总后以年为单位进行数值比较 已有1人参与

我现在要做一个微分方程组的参数和初值拟合,求解的时间单位是周,但我目前有的真实数据是以年为单位的,现在的思路就是对微分方程组求解结果按每52周进行汇总,再与相应的年实际值相减,以残差平方和最小为最优,用的是matlab中的lsqnonlin,fmincon也尝试了一下,但都是运行1个小时都出不来结果,也没提示错误。如果我直接用一组可能的参数代入微分方程组求解的话只需不到1分钟就能求出数值解。不知道问题出在哪儿,我是新手,请高手帮忙分析下原因呢。另外,lsqnonlin和fmincon能否用于这种情况,不行的话有没有其它思路可以做?先谢谢各位啦!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
dbb627: 金币+1, 谢谢应助! 2012-04-12 09:08:50
tomangel: 金币+1, 有帮助, 确实对初值依赖很严重,改后就能拟合出结果了 2012-04-13 23:28:22
微分方程拟合一般是比较耗时的,Matlab的lsqnonlin和fmincon函数对初值依赖严重,猜的不好更会严重降低计算效率,可以试试1stOpt,本版有一些其做微分方程拟合的帖子,搜索参考下。
2楼2012-04-12 08:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tomangel 的主题更新
信息提示
请填处理意见