| 查看: 2529 | 回复: 11 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
求高版本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
» 猜你喜欢
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有251人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
-
3月 中国石油和化学工业联合会科学技术奖 马上开始，申获奖要点
已经有15人回复
-
中国石油和化学工业联合会科学技术奖的评审标准是什么？
已经有5人回复
-
经费如何使用，花不完怎么办？
已经有5人回复
-
兼职
已经有12人回复
-
海外Gustau Catalan课题组情况了解
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
4楼2013-04-09 21:31:23
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68184.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.5小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2013-04-09 10:58:55
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68184.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.5小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2013-04-09 11:03:06
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68184.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.5小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
★ ★
fegg7502: 金币+2, 鼓励交流 2013-04-11 10:35:08
fegg7502: 金币+2, 鼓励交流 2013-04-11 10:35:08
|
第一个方程 %--------------------------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------------------------------------------- |
5楼2013-04-09 21:47:34
