应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求高版本1stopt拟合，
51/1返回列表
查看: 1907  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yupeifeng

金虫 (小有名气)

[求助] 求高版本1stopt拟合，

给出拟合后的图像及参数、相关系数值，以便确定选择哪种模型。程序已经编好，只需运行后，把拟合后的图像，参数值及相关系数R复制过来就行了，谢谢！！
一、T=40℃
第一个程序：
Constant c0=6.9401,cp=5.8399;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(1/3)*(c-cp)
Data;
t,c
0        6.94
1        6.45
2        6.35
4        6.2
12        5.96
24        5.85
48        5.84

第二个程序
Constant c0=6.9401,cp=5.8399;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(2/3)*(c-cp)^p
Data;
t,c
0        6.94
1        6.45
2        6.35
4        6.2
12        5.96
24        5.85
48        5.84

第三个程序
Constant c0=6.9401,cp=5.8399;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(4/3)*(c-cp)^p
Data;
t,c
0        6.94
1        6.45
2        6.35
4        6.2
12        5.96
24        5.85
48        5.84

第四个程序
Constant c0=6.9401,cp=5.8399;
Variable t,c;
ODEFunction c'=-k*(c-cp)^n
Data;
t,c
0        6.94
1        6.45
2        6.35
4        6.2
12        5.96
24        5.85
48        5.84

二、T=30℃
第一个程序：
Constant c0=6.9401,cp=6.0499;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(1/3)*(c-cp)
Data;
t,c
0        6.94
2        6.39
4        6.3
8        6.2
12        6.1
48        6.05

第二个程序
Constant c0=6.9401,cp=6.0499;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(2/3)*(c-cp)^p
Data;
t,c
0        6.94
2        6.39
4        6.3
8        6.2
12        6.1
48        6.05

第三个程序
Constant c0=6.9401,cp=6.0499;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(4/3)*(c-cp)^p
Data;
t,c
0        6.94
2        6.39
4        6.3
8        6.2
12        6.1
48        6.05

第四个程序
Constant c0=6.9401,cp=6.0499;
Variable t,c;
ODEFunction c'=-k*(c-cp)^n
Data;
t,c
0        6.94
2        6.39
4        6.3
8        6.2
12        6.1
48        6.05

三、T=50℃
第一个程序：
Constant c0=6.7901,cp=5.7899;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(1/3)*(c-cp)
Data;
t,c
0        6.79
1        6.39
2        6.25
4        6.1
12        5.9
22        5.8
48        5.79
第二个程序：
Constant c0=6.7901,cp=5.7899;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(2/3)*(c-cp)^p
Data;
t,c
0        6.79
1        6.39
2        6.25
4        6.1
12        5.9
22        5.8
48        5.79
第三个程序：
Constant c0=6.7901,cp=5.7899;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(4/3)*(c-cp)^p
Data;
t,c
0        6.79
1        6.39
2        6.25
4        6.1
12        5.9
22        5.8
48        5.79
第四个程序：
Constant c0=6.7901,cp=5.7899;
Variable t,c;
ODEFunction c'=-k*(c-cp)^n
Data;
t,c
0        6.79
1        6.39
2        6.25
4        6.1
12        5.9
22        5.8
48        5.79[ Last edited by fegg7502 on 2013-4-11 at 10:34 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 1stopt.doc
    • 2013-04-09 09:21:26, 35 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
憧憬未来
1楼 2013-04-09 09:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hd2015

铜虫 (初入文坛)
★
fegg7502: 金币+1, 鼓励交流 2013-04-11 10:35:34
文件：数据文件-1
No        t        目标c        计算c
1        1        6.39        6.42163541587917
2        2        6.25        6.24143714580419
3        4        6.1        6.06833016513733
4        12        5.9        5.88880471181152
5        22        5.8        5.84059209475139
6        48        5.79        5.81017603235905
赞一下 回复此楼
9楼2013-04-10 21:43:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★
csgt0: 金币+2, 应助指数+1, 谢谢 2013-04-09 11:07:07
手头没有高版本的1stopt，用MATLAB作了一下，因为一般自编的MATLAB程序拟合结果被初值影响很大，第三情况方程的刚性比较大，结果仅供参考。T=40℃的情况，分别见附图1,2,3,4。目前看来，第四个方程比较好。

附图1.jpg



附图2.jpg



附图3.jpg



附图4.jpg
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2013-04-09 10:58:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
此外，第一个和四个方程，可以先积分为代数方程，再用低版本1stopt或者origin spss excel等拟合也可。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2013-04-09 11:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
3楼: Originally posted by 月只蓝 at 2013-04-09 11:03:06
此外，第一个和四个方程，可以先积分为代数方程，再用低版本1stopt或者origin spss excel等拟合也可。

可不可以让我看看你的程序怎么编写的啊，谢谢了
赞一下 回复此楼
憧憬未来
4楼2013-04-09 21:31:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿上海交大生物与医药专硕324分，求调剂 +3 jiajunX 2026-03-22 3/150 2026-03-22 19:32 by brblmd
[考研] 生物学071000 329分求调剂 +5 我爱生物生物爱� 2026-03-17 5/250 2026-03-22 16:42 by tcx007
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +3 niko- 2026-03-22 3/150 2026-03-22 15:15 by 杨杨杨紫
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 材料与化工（0856）304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 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
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
关闭
请填处理意见
关闭