| 查看: 743 | 回复: 1 | ||
luckyman0125铁杆木虫 (职业作家)
|
[求助]
求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数
|
|
各位大神，请问如何用matlab/1stOpt 自定义微分方程来拟合实验数据来求方程中的参数。 本人实验数据如下 x = [0 600 1200 1800 2400 3600 5400 7200 14400 18000 43200 57600 86400 136800 151200] y =[0.026 0.022 0.02 0.019 0.018 0.0175 0.0165 0.0158 0.015 0.0145 0.013 0.0125 0.0118 0.0105 0.0103] 需要用 dy/dx=r^2/(0.1424*y) * (588*y - 2.16*10^(-3)*cos(t)/r)拟合来求出 r 和 t。 如果用公式 dy/dx=(1.58*10^(-7))/y * (588*y - 14.4*cos(t)), 拟合 t，又是多少 谢谢 [ Last edited by luckyman0125 on 2013-6-13 at 00:35 ] |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有3人回复
-
关于代码变化问题，想知道的进来
已经有16人回复
-
fileCode有新解读？
已经有18人回复
-
这样的filecode谁见过
已经有15人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- origin 非线性自定义拟合问题~~方程组编写及运行~~
已经有8人回复
- 求助，matlab自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- matlab高手帮忙做个自定义复杂函数拟合，含有复数，
已经有25人回复
- 用MATLAB或者ORIGIN怎么拟合复数数据
已经有10人回复
- matlab求解方程中的参数
已经有21人回复
- 请问matlab拟合二元方程的问题
已经有3人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- origin8.0拟合方程中参数怎么求
已经有5人回复
- 一阶微分方程求解，并优化参数值
已经有14人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- 学1stopt要多久？
已经有12人回复
- 高金求助matlab解微分方程组
已经有12人回复
- 求教origin拟合简单的渗流方程（又称逾渗方程），做复合材料的应该会遇到吧
已经有20人回复
- 怎么由测试数据得到拟合方程
已经有11人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 用1stopt进行数据拟合时，无法运行sharedmodel ,求助高手
已经有7人回复
- 求助Origin8.0拟合方程
已经有6人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
luckyman0125: 金币+20, ★★★很有帮助 2013-06-13 09:29:45
感谢参与，应助指数 +1
luckyman0125: 金币+20, ★★★很有帮助 2013-06-13 09:29:45
|
1stOpt求解： 1：dy/dx=r^2/(0.1424*y) * (588*y - 2.16*10^(-3)*cos(t)/r) 均方差(RMSE):0.00756626515276929 残差平方和(SSE):0.000801477157068152 相关系数(R): 0.792802636782712 相关系数之平方(R^2): 0.628536020889622 确定系数(DC): -3.76235514275764 F统计(F-Statistic): 2.51976168099286 参数 最佳估算 -------------------- ------------- r -1.38058887495168E-7 t -3.14159293861944 ====== 结果输出 ====== 文件: 数据文件-1 No x 目标 y 计算 y 1 600 0.022 0.0259516756931473 2 1200 0.02 0.0259032611459708 3 1800 0.019 0.0258547558513642 4 2400 0.018 0.0258061592974539 5 3600 0.0175 0.0257086903400129 6 5400 0.0165 0.0255617893810344 7 7200 0.0158 0.0254140384823027 8 14400 0.015 0.024814230435684 9 18000 0.0145 0.0245088170910798 10 43200 0.013 0.0222537567477591 11 57600 0.0125 0.0208558624037211 12 86400 0.0118 0.0177322184557397 13 136800 0.0105 0.0101712955998023 14 151200 0.0103 0.00656991795119879 2： dy/dx=(1.58*10^(-7))/y * (588*y - 14.4*cos(t)) 周期性函数会导致无穷多个t，效果不好！ |
2楼2013-06-13 09:02:58
20