24小时热门版块排行榜    

查看: 2029  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

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

至尊木虫 (著名写手)

^^

【答案】应助回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 085602 289分求调剂 +5 WWW西西弗斯 2026-03-24 5/250 2026-03-24 18:51 by jhhcooi
[考研] 招08考数学 +7 laoshidan 2026-03-20 16/800 2026-03-24 17:52 by 乌拉儿山脉
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 327求调剂 +5 prayer13 2026-03-23 5/250 2026-03-23 22:11 by 星空星月
[考研] 材料专硕英一数二306 +8 z1z2z3879 2026-03-18 8/400 2026-03-23 20:49 by baobaoye
[考研] 求调剂材料学硕080500,总分289分 5+3 @taotao 2026-03-19 21/1050 2026-03-23 10:17 by 冠c哥
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
信息提示
请填处理意见