24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1723  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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

» 猜你喜欢

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

憧憬未来
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见