| 查看: 2253 | 回复: 5 | ||||
Daisy915木虫 (职业作家)
|
[求助]
怎么做微分方程的数据拟合
|
|
实验中碰到一个微分函数无法积分成原函数，dr/dt = A exp (B/r)， 我有实验数据是关于r-t的，想要用这个函数进行拟合看看能不能匹配，并得到A和B的值。因为以前都是用Origin对公式直接拟合的，碰到微分函数就不知道该怎么办了。木虫里的大神请拔刀相助，谢谢啦！ 下面是我的实验数据 time (t) radius (r) 61.42 3.78881 61.57 3.86081 61.71 3.93148 61.85 3.79123 62 3.80332 62.14 3.76937 62.28 3.75596 62.42 3.75598 62.57 3.88922 62.71 3.80816 62.85 3.85962 63 3.79245 63.14 3.81537 63.28 3.6682 63.42 3.82258 63.57 3.72037 63.71 3.66821 63.85 3.74374 64 3.69438 64.14 3.42358 64.28 3.26156 64.42 3.3777 64.57 3.23189 64.71 3.23189 64.85 3.48075 65 3.36408 65.14 3.40612 65.28 3.32154 65.42 3.38448 65.57 3.42493 65.71 3.33669 65.85 3.45162 66 3.40746 66.14 3.35725 66.28 3.24888 66.42 3.29934 66.57 3.09405 66.71 3.02353 66.85 3.12944 67 3.11917 67.14 3.12357 67.28 3.03715 67.42 3.06573 67.57 3.022 67.71 2.98535 67.85 3.05373 68 3.09405 68.14 3.0068 68.28 2.72494 68.42 2.56355 68.57 2.47243 68.71 2.54738 68.85 2.38167 69 2.26916 69.14 2.08801 69.28 2.2732 69.42 1.83533 69.57 1.91132 69.71 1.84036 69.85 1.57131 70 1.14966 70.14 1.12135 70.28 0 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有127人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 常微分方程参数拟合问题，matlab和1stopt均可
已经有3人回复
- 请教一下各位大神怎么求解 Runge–Kutta method的微分方程的各个系数
已经有12人回复
- 跪求高手帮忙看看这组数据如何进行拟合
已经有7人回复
- 求助1stopt拟合动力学参数
已经有4人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- matlab拟合参数显示初始点是一个局部最小值
已经有18人回复
- 拟合方程并作图
已经有6人回复
- 怎样拟合数据求出x=f(A,B,C),y=g(A,B,C,D)中的A,B,C,D参数
已经有18人回复
- 两个自变量怎么拟合求参数值呢？
已经有7人回复
- 求助，matlab自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 请问这篇文章的拟合方程方程是什么？是怎么拟合的？用的什么软件？
已经有4人回复
- 通过实验数据拟合，求解公式中的参数
已经有14人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- 【公式拟合问题】如何拟合下面这样的一个公式，多谢大牛。
已经有6人回复
- 怎么由测试数据得到拟合方程
已经有11人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
2楼2015-05-09 14:36:18
3楼2015-05-09 14:50:00
ybkooo
至尊木虫 (著名写手)
^^
- 应助: 32 (小学生)
- 金币: 11056.5
- 散金: 178
- 红花: 15
- 帖子: 2312
- 在线: 545.4小时
- 虫号: 2710519
- 注册: 2013-10-09
- 性别: GG
- 专业: 控制理论与方法
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Daisy915: 金币+15, ★★★★★最佳答案 2015-05-09 19:10:36
Daisy915: 金币+20, ★★★★★最佳答案, 太感谢了，我在这边运行看看。。 2015-05-09 19:12:48
感谢参与，应助指数 +1
Daisy915: 金币+15, ★★★★★最佳答案 2015-05-09 19:10:36
Daisy915: 金币+20, ★★★★★最佳答案, 太感谢了，我在这边运行看看。。 2015-05-09 19:12:48
|
1stOpt代码 Parameters a,b; Variable t,r; ODEFunction r'=a*exp(b/r); Data; 61.42 3.78881 61.57 3.86081 61.71 3.93148 61.85 3.79123 62 3.80332 62.14 3.76937 62.28 3.75596 62.42 3.75598 62.57 3.88922 62.71 3.80816 62.85 3.85962 63 3.79245 63.14 3.81537 63.28 3.6682 63.42 3.82258 63.57 3.72037 63.71 3.66821 63.85 3.74374 64 3.69438 64.14 3.42358 64.28 3.26156 64.42 3.3777 64.57 3.23189 64.71 3.23189 64.85 3.48075 65 3.36408 65.14 3.40612 65.28 3.32154 65.42 3.38448 65.57 3.42493 65.71 3.33669 65.85 3.45162 66 3.40746 66.14 3.35725 66.28 3.24888 66.42 3.29934 66.57 3.09405 66.71 3.02353 66.85 3.12944 67 3.11917 67.14 3.12357 67.28 3.03715 67.42 3.06573 67.57 3.022 67.71 2.98535 67.85 3.05373 68 3.09405 68.14 3.0068 68.28 2.72494 68.42 2.56355 68.57 2.47243 68.71 2.54738 68.85 2.38167 69 2.26916 69.14 2.08801 69.28 2.2732 69.42 1.83533 69.57 1.91132 69.71 1.84036 69.85 1.57131 70 1.14966 70.14 1.12135 70.28 0  1.png  2.png |
4楼2015-05-09 16:35:41
ybkooo
至尊木虫 (著名写手)
^^
- 应助: 32 (小学生)
- 金币: 11056.5
- 散金: 178
- 红花: 15
- 帖子: 2312
- 在线: 545.4小时
- 虫号: 2710519
- 注册: 2013-10-09
- 性别: GG
- 专业: 控制理论与方法
【答案】应助回帖
|
Calculated r Target r 3.77210449133956 3.86081 3.75635769497891 3.93148 3.74045695991814 3.79123 3.72324464105732 3.80332 3.70701083370589 3.76937 3.69060893678837 3.75596 3.67403388789241 3.75598 3.65607672129263 3.88922 3.63912580050396 3.80816 3.62198454477008 3.85962 3.60340065950934 3.79245 3.58584506849759 3.81537 3.56807905033821 3.6682 3.55009532081249 3.82258 3.53057669247529 3.72037 3.5121169967516 3.66821 3.49341448921915 3.74374 3.47309619387611 3.69438 3.4538608490953 3.42358 3.43435260926003 3.26156 3.4145603736329 3.3777 3.39302578830312 3.23189 3.37260668132846 3.23189 3.35186453687286 3.48075 3.32926573620076 3.36408 3.3078067836049 3.40612 3.28597638963846 3.32154 3.26375624896833 3.38448 3.23949393639557 3.42493 3.21640161324879 3.33669 3.19285325329176 3.45162 3.16708773593726 3.40746 3.14251076721507 3.35725 3.11739141353791 3.24888 3.0916958058356 3.29934 3.06348278615381 3.09405 3.03647001492914 3.02353 3.00875266256845 3.12944 2.97821457584951 3.11917 2.94886547066328 3.12357 2.91863091563023 3.03715 2.88743558522358 3.06573 2.85284779077552 3.022 2.8193736152616 2.98535 2.78462968632804 3.05373 2.7458400296622 3.09405 2.70800502912818 3.0068 2.66839622632053 2.72494 2.62677769984569 2.56355 2.57962526958033 2.47243 2.53283767942265 2.54738 2.48288616570424 2.38167 2.42514910972137 2.26916 2.36642797214845 2.08801 2.30181572655192 2.2732 2.22953277025482 1.83533 2.14033846357106 1.91132 2.04077856429512 1.84036 1.9147536823837 1.57131 1.71574877516995 1.14966 1.06629364635781 1.12135 0.287010831972625 0 |
5楼2015-05-09 16:36:54
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
感谢参与，应助指数 +1
|
1. 1stopt软件直接计算，需要高版本的1stopt软件，代码简单，效果好 2. MATLAB 编程：http://muchong.com/bbs/viewthread.php?tid=7603645 |
6楼2015-05-09 17:06:38
投票:
10