版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yupeifeng

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 960.3
散金: 34
红花: 1
帖子: 220
在线: 78.4小时
虫号: 1421266
注册: 2011-09-28
专业: 无机材料化学

[求助] 求高版本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

» 猜你喜欢

基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？已经有2人回复
B区学生调剂-兰州交通大学材料科学与工程学院已经有12人回复
物理学I论文润色/翻译怎么收费? 已经有66人回复
山西大同大学物理学专业还有调剂名额，欢迎调剂！已经有13人回复
桂林理工大学物理学专业招收调剂已经有18人回复
VASP 的一组 GPU / CPU 基准测试记录已经有0人回复
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名已经有0人回复
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树已经有0人回复
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）已经有4人回复

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

1stopt中使用sharedmodel拟合曲线出现问题，求指示！！已经有6人回复
1stOPT1.5版每次仿真结果不一样，求助！急！已经有3人回复
用1stopt拟合非线性方程结果与其他软件拟合结果差异大已经有5人回复
求高版本1stopt，解决一个问题，微分方程组参数拟合已经有3人回复
运算过程有复数出现，求高版本1stOpt 已经有6人回复
1stOpt程序的一个小问题已经有5人回复
大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？已经有8人回复
有1stOpt4.0版本的朋友帮忙拟合一下复数吧感激不尽呐已经有4人回复
1stopt4阶非线性常微分方程，帮小弟看一下，跪求。已经有7人回复
1STOPT可以拟合一下形式的方程组吗？已经有5人回复
1stOpt高版本求助！求高版本的帮忙走一下谢了！已经有16人回复
1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？已经有13人回复
跪求1stOpt2.5以上版本，或者请高人帮忙拟合包含复数的函数的参数已经有8人回复
学1stopt要多久？已经有12人回复
哪位大哥能帮我用1stopt3.0以上版本跑一个拟合程序已经有11人回复
用origin或1stopt拟合一条隐函数方程曲线已经有4人回复
用1stopt进行曲线拟合每次给出结果不一样已经有6人回复
用1stopt进行数据拟合时，无法运行sharedmodel ,求助高手已经有7人回复
求助关于1stopt拟合~··~~ 已经有4人回复
1stOPt拟合结果已经有4人回复

憧憬未来

1楼 2013-04-09 09:25:30

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68132.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

此外，第一个和四个方程，可以先积分为代数方程，再用低版本1stopt或者origin spss excel等拟合也可。

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

3楼2013-04-09 11:03:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 12 个回答

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68132.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★
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的回帖

yupeifeng

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 960.3
散金: 34
红花: 1
帖子: 220
在线: 78.4小时
虫号: 1421266
注册: 2011-09-28
专业: 无机材料化学

引用回帖:

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的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68132.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

★ ★
fegg7502: 金币+2, 鼓励交流 2013-04-11 10:35:08

引用回帖:

4楼: Originally posted by yupeifeng at 2013-04-09 21:31:23
可不可以让我看看你的程序怎么编写的啊，谢谢了...

第一个方程
%--------------------------start-----------------------------------------
function ode2
clear all;clc

format long

tspan=[0 1 2 4 12 24 48];
t1=length(tspan),
yexp=[6.45 6.35 6.2 5.96 5.85 5.84]';
t2=length(yexp)

y0=6.94;

k0=[2];
lb=[-1]*1e3;
ub=[1]*1e3;

yy=[y0 yexp'];

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp);
ci = nlparci(k,residual,jacobian);
fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n')
fprintf('\t反应速率常数 k = %.4f\n',k(1))
%fprintf('\t反应级数 n = %.4f\n',k(2))
fprintf('  The sum of the squares is: %.1e\n\n',resnorm)
ts=0:1:max(tspan);

[ts ys]=ode45(@KineticsEqs,ts,y0,[],k);
[ttt XXsim] = ode45(@KineticsEqs,tspan,y0,[],k);
y=XXsim(2:end);
xexp=yexp;
R2=1-sum((xexp-y).^2)./sum((xexp-mean(y)).^2);
fprintf('\n\t相关系数之平方R^2 = %.6f',R2);
plot(ts,ys,'b',tspan,yy,'or'),legend('计算值','实验值','Location','best'),

%---------------------------------------------------------
function f = ObjFunc(k,tspan,y0,yexp)          % 目标函数
[t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k);
ysim = Xsim(2:end);
f=ysim-yexp;
%----------------------------------------------------------

function dydt = KineticsEqs(t,y,k)
c0=6.9401;cp=5.8399;
dydt=-k(1)*(c0-y).^(1/3).*(y-cp);
%--------------------The end-------------------------------------------

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

5楼2013-04-09 21:47:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 12 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 求调剂 +4	苦命人。。。 2026-04-18	4/200	2026-04-19 02:01 by 烟雨流涯
[考研] 接受任何调剂 +6	也就是栗子 2026-04-17	7/350	2026-04-18 17:20 by 涵竹刘
[考博] 申博/考博 +3	啃面包的小书虫 2026-04-17	4/200	2026-04-17 23:54 by 阳阳阳^_^
[考研] 求调剂 +9	小聂爱学习 2026-04-16	11/550	2026-04-17 22:34 by chixmc
[考博] 求博导｜生物质基多孔碳/超级电容方向，已有相关成果，寻能源材料/碳材料方向老师 +3	猪猪人Zzz 2026-04-12	3/150	2026-04-17 19:10 by 阳阳阳^_^
[有机交流] 二苯甲酮酸类衍生物 50+3	小白爱主人 2026-04-17	6/300	2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6	tekuzu 2026-04-17	6/300	2026-04-17 13:48 by Espannnnnol
[论文投稿] 有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业 20+4	之护着 2026-04-16	5/250	2026-04-17 10:02 by bobvan
[考研] 291求调剂 +11	关忆北. 2026-04-14	11/550	2026-04-16 15:18 by jiahl2024
[考研] 297，工科调剂?河南农业大学本科 +14	河南农业大学-能 2026-04-14	14/700	2026-04-16 14:41 by dingyanbo1
[考研] 322求调剂 +8	123安康 2026-04-12	15/750	2026-04-16 11:07 by Espannnnnol
[考研] 327求调剂 +26	Xxjc1107. 2026-04-13	29/1450	2026-04-16 10:52 by Espannnnnol
[考研] 289 分105500药学专硕求调剂(找B区学校) +4	白云123456789 2026-04-13	4/200	2026-04-16 00:18 by 粉沁若尘
[考研] 药学求调剂 +11	RussHu 2026-04-12	13/650	2026-04-15 19:07 by zhuwenxu
[考研] 求调剂 +12	何气正 2026-04-13	13/650	2026-04-14 14:47 by zs92450
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6	王ymaa 2026-04-13	13/650	2026-04-14 11:33 by 王ymaa
[考研] 085600材料与化工329分求调剂 +24	叶zilin 2026-04-13	25/1250	2026-04-14 09:20 by 试管破裂
[考研] 2026硕士调剂_能动_河南农业大学 +4	河南农业大学-能 2026-04-12	4/200	2026-04-13 22:01 by bljnqdcc
[考研] 297工科，求调剂? +13	河南农业大学-能 2026-04-12	13/650	2026-04-13 14:12 by dingyanbo1
[考研] 一志愿浙大生物325分求调剂 +9	zysheng 2026-04-12	9/450	2026-04-12 22:31 by yuyin1233