24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1644  |  回复: 8

lyh315216

新虫 (初入文坛)

[求助] 用1stopt同时拟合两组数据,常微分方程,但是数据点不同 已有2人参与

用1stopt同时拟合两组数据,常微分方程,但是数据点不同,具体如下面
Parameter k[0,], KA[0,], KB[0,100];
Variable t, y;
ODEFunction
y'=-k*KA*KB^3*y*0.4^3/((1+KA*y)*(1+KB^0.5*0.4^0.5)^6);
第一组数据
240        1.15
360        1.062
480        0.97
600        0.886
720        0.823
960        0.746
1080 0.706
第二组数据
480        1.031
720        0.958
840        0.916
960        0.896
1080  0.866
1200  0.850
求助,怎么把这两组数据同时拟合啊??? @月只蓝
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有3个 )

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
月只蓝: 回帖置顶 2017-03-16 09:58:20
把“第一组数据”,”第二组数据”都换成“Data;”就可以直接计算了。

[ 发自手机版 http://muchong.com/3g ]
2楼2017-03-16 07:55:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
使用具体常微分拟合的版本,代码按dingd说的那样改就对了。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2017-03-16 09:58:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
5楼: Originally posted by lyh315216 at 2017-03-16 10:43:00
单独拟合都能拟合出来,但是两组数据一起拟合就无法拟合出来,请问是我的数据或者公式的问题吗?还是有什么需要设定的地方我没设定?
Title As 75℃;
Parameter k, KA, KB;
Variable t, y;
ODEFunction
y'=-k ...

复制你5楼的代码,不做任何修改,可以运行。
用1stopt同时拟合两组数据,常微分方程,但是数据点不同
20170316121734.png

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
8楼2017-03-16 12:18:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

lipenggg

专家顾问 (正式写手)

一团和气
3楼2017-03-16 09:49:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyh315216

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2017-03-16 07:55:04
把“第一组数据”,”第二组数据”都换成“Data;”就可以直接计算了。

单独拟合都能拟合出来,但是两组数据一起拟合就无法拟合出来,请问是我的数据或者公式的问题吗?还是有什么需要设定的地方我没设定?
Title As 75℃;
Parameter k[0,20], KA[0,], KB[0,];
Variable t, y;
ODEFunction
y'=-k*KA*KB^3*y*0.4^3/((1+KA*y)*(1+KB^0.5*0.4^0.5)^6);
Data;
240        1.15
360        1.062
480        0.97
600        0.886
720        0.823
960        0.746
1080        0.706
Data;
240        1.118
480        1.031
600        0.966
720        0.958
840        0.916
960        0.896
1080        0.866
1200        0.85
5楼2017-03-16 10:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
5楼: Originally posted by lyh315216 at 2017-03-16 10:43:00
单独拟合都能拟合出来,但是两组数据一起拟合就无法拟合出来,请问是我的数据或者公式的问题吗?还是有什么需要设定的地方我没设定?
Title As 75℃;
Parameter k, KA, KB;
Variable t, y;
ODEFunction
y'=-k ...

代码按dingd说的那样改完之后,我已经测试过,运行没问题。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
6楼2017-03-16 11:05:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyh315216

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by 月只蓝 at 2017-03-16 11:05:50
代码按dingd说的那样改完之后,我已经测试过,运行没问题。...

为啥我就拟合不出来呢,我的程序写的有问题吗请问?
7楼2017-03-16 11:50:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyh315216

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by 月只蓝 at 2017-03-16 12:18:59
复制你5楼的代码,不做任何修改,可以运行。

20170316121734.png
...

实验值和拟合值相差很大啊,那就是说是数据的问题了?
9楼2017-03-16 21:19:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lyh315216 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工学硕294分,求导师收留 +27 yzyzx 2026-04-12 30/1500 2026-04-14 05:42 by 西北望—风沙
[考研] 310求调剂 +15 666真好 2026-04-11 17/850 2026-04-13 22:27 by pies112
[考研] 279学硕食品专业求调剂院校 20+6 孤独的狼爱吃羊 2026-04-12 23/1150 2026-04-13 20:42 by biomen
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] 一志愿厦大生物学332求调剂 +11 池池池池池池 2026-04-08 11/550 2026-04-13 14:10 by 科研论
[考研] 一志愿华工085600 331分 +7 天下ww 2026-04-09 7/350 2026-04-13 09:01 by lhj2009
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen~
[考研] 药学专硕调剂 +8 ? 一路生?花? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[教师之家] 请问地理、遥感方面,可以做哪些横向项目啊,纵向完不成考核啊 +3 锦衣卫寒战 2026-04-07 5/250 2026-04-11 20:51 by 豫椒
[考研] 11408。358求调剂 +3 TMYzds 2026-04-07 3/150 2026-04-11 17:10 by 氮气气气
[考研] 269求调剂 +11 啊啊我我 2026-04-07 11/550 2026-04-11 16:45 by vgtyfty
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 298求调剂 +5 残荷新柳 2026-04-07 5/250 2026-04-11 11:02 by 紫曦紫棋
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 一志愿京区985,085401电子信息,本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 材料工程322 +18 哈哈哈吼吼吼哈 2026-04-07 19/950 2026-04-09 10:44 by cymywx
信息提示
请填处理意见