24小时热门版块排行榜    

查看: 2241  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 14/700 2026-08-17 14:38 by Ldrop2023
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6GojgJvkDudM 2026-08-16 6/300 2026-08-17 13:58 by EYK67XJMj64E
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 13:49 by EYK67XJMj64E
[基金申请] filecode=后面第一个是大写字母 +7 wangze12014 2026-08-14 8/400 2026-08-17 12:45 by a441914426
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
请填处理意见