24小时热门版块排行榜    

查看: 351  |  回复: 2

清辉12365

新虫 (初入文坛)

[求助] 弄动力学参数拟合参数,不知道怎么办了,哪个大神帮忙弄下代码

下面是我上传的问题和数据,看了几个月的实用化工计算机模拟没有进展,请高手帮忙一下,
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 问题.docx
  • 2015-03-02 08:46:05, 19.53 K
  • 附件 2 : 数据.xlsx
  • 2015-03-02 08:46:11, 10.7 K

» 猜你喜欢

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

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

清辉12365

新虫 (初入文坛)

帮顶,,
2楼2015-03-02 09:14:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清辉12365

新虫 (初入文坛)

1-初始量,2-温度,3-压力,4-甲醇摩尔分数x1,5-CO摩尔分数x2,6-O2摩尔分数X3,8-产物DMC摩尔数y1   9-产物CO2摩尔分数y2
上面的1、2、3、4、5、6、8、9在Donglixue.mat中
我的编程是:function   yyh
clear all;clc
load Donglixue;
%输入数据
NO=Donglixue(:,1);  %输入初始量
ym0=Donglixue(:,4);   %输入甲醇摩尔比
yco0=Donglixue(:,5);   %co的摩尔比
yo2=Donglixue(:,6);    %o2的摩尔比
ydmc=Donglixue(:,8);    %实验的DMC的摩尔比
yco2=Donglixue(:,9);    %实验的CO2的摩尔比
t=Donglixue(:,2);    %输入温度
T=t+273.15;           %温度换成K
P=Donglixue(:,2);     %输入压力
%lsqnonlin函数拟合
beta0=[10000 500 0 0 0 10000 10000 0.2 -0.3];
lb=[0 0 -2 -2 -2 0 0 -2 -2]; %参数的下限
Y0=[0 0];
tspan=[0,0.375];     %积分限
[beta,resonrm,residual,exitflag,output,lambda,jacobian]=...
    lsqnonlin(@fun1,beta0,lb,[],[],tspan,Y0,NO,ym0,yco0,yo2,ydmc,yco2,T,P)
function f=fun1(beta,tspan,Y0,NO,ym0,yco0,yo2,ydmc,yco2,T,P)
   [t,Y]=ode45(@fun2,tspan,Y0,[],beta);
   f=[Y(:,1)-ydmc;Y(:,2)-yco2];
    function ff=fun2(t,Y,NO,ym0,yco0,yo2,ydmc,yco2,T,P)
   ff(1)=((1+1.5*Y(:,1)+0.5*Y(:,2)).*(1+1.5*Y(:,1))*...
                                       (ym0+Y(:,1).*(1.5*ym0-2)+0.5*ym0.*Y(:,2)).^beta(3)...
       .*P.^(beta(3)+beta(4)+beta(5)).*(yco0+Y(:,1).*(1.5*yco0-1)+Y(:,2).*(0.5*yco0-1)).^beta(4).*...
                                       (yo2+0.5.*Y(:,1)*(3*yco0-1)+0.5*Y(:,2)*(yo2-1)).^beta(5)...
                                       .*beta(1)*exp(-beta(2)./T))
                                       +0.5*Y(:,1).*beta(6)*exp(-beta(7)./T).*P.^(beta(8)+beta(9)).*...
                                         (yco0+Y(:,1).*(1.5*yco0-1)+Y(:,2).*(0.5*yco0-1)).^beta(8).*...
                                         (yo2+0.5.*Y(:,1)*(3*yco0-1)+0.5*Y(:,2)*(yo2-1)).^beta(9)))./NO
ff(2)=(1+1.5*Y(:,1)+0.5*Y(:,2)).* ((0.5*Y(:,2)+1).*beta(6)*exp(-beta(7)./T).*P.^(beta(8)+beta(9)).*...
                                                   (yco0+Y(:,1).*(1.5*yco0-1)+Y(:,2).*(0.5*yco0-1)).^beta(8).*...
                                                   (yo2+0.5.*Y(:,1)*(3*yco0-1)+0.5*Y(:,2)*(yo2-1)).^beta(9)+...
                                                   1.5*Y(:,2)*...
                                                   (ym0+Y(:,1).*(1.5*ym0-2)+0.5*ym0.*Y(:,2)).^beta(3).*...
                                                   (yco0+Y(:,1).*(1.5*yco0-1)+Y(:,2).*(0.5*yco0-1)).^beta(4).*...
                                                   (yo2+0.5.*Y(:,1)*(3*yco0-1)+0.5*Y(:,2)*(yo2-1)).^beta(5).*P.^(beta(3)+beta(4)+beta(5).*...
                                                   beta(1)*exp(-beta(2)./T)./NO

因为r1、r2式子比较麻烦,所以看的式子比较乱,请大神们指导一下,
弄动力学参数拟合参数,不知道怎么办了,哪个大神帮忙弄下代码
1.jpg


弄动力学参数拟合参数,不知道怎么办了,哪个大神帮忙弄下代码-1
2.jpg

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Donglixue.mat
  • 2015-03-02 16:23:20, 1 K
3楼2015-03-02 16:23:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 清辉12365 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考研] 316求调剂 +5 梁茜雯 2026-03-19 5/250 2026-03-20 21:26 by 你好你好666
[考研] 279分求调剂 一志愿211 +10 chaojifeixia 2026-03-19 11/550 2026-03-20 20:32 by JourneyLucky
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 8/400 2026-03-20 15:58 by babero
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 307求调剂 +9 冷笙123 2026-03-17 9/450 2026-03-19 22:44 by 学员8dgXkO
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 0817调剂 +3 没有答案_ 2026-03-14 3/150 2026-03-19 09:51 by Xu de nuo
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
信息提示
请填处理意见