24小时热门版块排行榜    

查看: 2026  |  回复: 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
回复此楼

» 猜你喜欢

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

材料系博士在读/科技、艺术、推理小说爱好者
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zenghy789

铁虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
对导数直接拟合啊。具体的是,以差分代替微分,即(dr/dt)=(r2-r1)/(t2-t1)=Aexp(B/r);origin里面好像也可以直接画出导数的,好像。
2楼2015-05-09 14:36:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zenghy789

铁虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
Daisy915: 金币+20, 有帮助, 这个是实验数据,测量有浮动,没办法。 2015-05-09 14:52:55
还有,dr/dt = A exp (B/r)值要么正要么负,不可能有正有负,所以你的数据时有问题的
3楼2015-05-09 14:50:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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


怎么做微分方程的数据拟合-1
2.png

abcd
4楼2015-05-09 16:35:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^

【答案】应助回帖

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
abcd
5楼2015-05-09 16:36:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
1. 1stopt软件直接计算,需要高版本的1stopt软件,代码简单,效果好
2. MATLAB 编程:http://muchong.com/bbs/viewthread.php?tid=7603645
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
6楼2015-05-09 17:06:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Daisy915 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854电子信息求调剂 +3 α____ 2026-03-22 4/200 2026-03-24 16:22 by α____
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 086003食品工程求调剂 +4 淼淼111 2026-03-24 4/200 2026-03-24 11:53 by 544594351
[考研] 求调剂 +5 研研,接电话 2026-03-24 6/300 2026-03-24 11:39 by 星空星月
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 材料专硕找调剂 +3 哈哈哈吼吼吼哈 2026-03-23 3/150 2026-03-23 23:13 by peike
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +5 niko- 2026-03-22 5/250 2026-03-23 22:01 by fuyu_
[考研] 一志愿武理材料工程348求调剂 +6  ̄^ ̄゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +5 jiajunX 2026-03-22 5/250 2026-03-23 18:07 by YMU施老师
[考研] 工科0856求调剂 +5 沐析汀汀 2026-03-21 5/250 2026-03-23 17:56 by 海瑟薇-
[考研] 北科281学硕材料求调剂 +8 tcxiaoxx 2026-03-20 9/450 2026-03-23 12:16 by tcxiaoxx
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
信息提示
请填处理意见