|
|
1-初始量,2-溫度,3-壓力,4-甲醇摩爾分數(shù)x1,5-CO摩爾分數(shù)x2,6-O2摩爾分數(shù)X3,8-產(chǎn)物DMC摩爾數(shù)y1 9-產(chǎn)物CO2摩爾分數(shù)y2
上面的1、2、3、4、5、6、8、9在Donglixue.mat中
我的編程是:function yyh
clear all;clc
load Donglixue;
%輸入數(shù)據(jù)
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函數(shù)擬合
beta0=[10000 500 0 0 0 10000 10000 0.2 -0.3];
lb=[0 0 -2 -2 -2 0 0 -2 -2]; %參數(shù)的下限
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式子比較麻煩,所以看的式子比較亂,請大神們指導(dǎo)一下,
![弄動力學(xué)參數(shù)擬合參數(shù),不知道怎么辦了,哪個大神幫忙弄下代碼]()
1.jpg
![弄動力學(xué)參數(shù)擬合參數(shù),不知道怎么辦了,哪個大神幫忙弄下代碼-1]()
2.jpg
|
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋:小木蟲僅提供交流平臺,不對該內(nèi)容負責(zé)。
本內(nèi)容由用戶自主發(fā)布,如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題,其責(zé)任在于用戶本人,如對版權(quán)有異議,請聯(lián)系郵箱:xiaomuchong@tal.com
- 附件 1 : Donglixue.mat
2015-03-02 16:23:20, 1 K
|