| 查看: 432 | 回复: 2 | ||
[求助]
弄动力学参数拟合参数，不知道怎么办了，哪个大神帮忙弄下代码
|
| 下面是我上传的问题和数据，看了几个月的实用化工计算机模拟没有进展，请高手帮忙一下， |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 问题.docx
- 附件 2 : 数据.xlsx
2015-03-02 08:46:05, 19.53 K
2015-03-02 08:46:11, 10.7 K
» 猜你喜欢
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有19人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这样的动力学方程能够有1stopt拟合出参数来么
已经有4人回复
- 使用MATLAB拟合动力学参数问题
已经有18人回复
- 1stopt数据拟合求助
已经有8人回复
- 跪求高手帮忙看看这组数据如何进行拟合
已经有7人回复
- 使用matlab最优化方法拟合获得多个动力学参数中的问题
已经有4人回复
- 求助1stopt拟合动力学参数
已经有4人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 拜求！matlab拟合动力学参数
已经有14人回复
- 关于两个参数的曲线拟合问题，向大神们求助
已经有3人回复
- 求助用matlab拟合动力学方程
已经有9人回复
- 使用lsqnonlin函数优化动力学参数，总是得不到合理的结果
已经有6人回复
- 求大家帮帮忙看看matlab拟合动力学方程的问题，已经走投无路了，求大家帮帮忙啊！
已经有19人回复
- 怎样用一组参数同时拟合两个曲线--matlab
已经有5人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 急求高手帮写个matlab程序拟合参数求动力学方程。
已经有4人回复
- 动力学参数拟合
已经有26人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
2楼2015-03-02 09:14:22
|
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  2.jpg |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Donglixue.mat
2015-03-02 16:23:20, 1 K
3楼2015-03-02 16:23:24
