[交流] 动力学参数拟合

19楼: Originally posted by feng198800 at 2012-06-05 22:42:55
谢谢了，可是我的动力学方程是每次实验的初始条件不同，也就是ode求解时的初值不同（温度T,流量V,进口湿基组成都不同）呀？...

function f = ObjFunc(k,P0,Pi)           % 目标函数
[m,n] = size(P0);
Pcal = zeros(m,n);
tspan =[0 0.0066];         % 即Wc/F0,kg.h/mol
for i = 1:m
[t PP] = ode45(@Euqations,tspan,P0(i,:),[],k);
Pcal(i,:) = PP(end,:);
end
f= Pcal-Pi;

% ------------------------------------------------------------------
function dPdt = Euqations(t, P, k)        % here t = Wc / F0
global keq1 keq2
m=k(6); n=k(7);
% 反应速率
r1 = k(1)*(P(6)*P(7)*P(7) - P(1)/keq1);
r2 = k(1)*P(1)*P(1);
r3 = k(3)*P(1)^(4-m)*P(2)^(m-2);
r4 = k(4)*P(6)^n*P(8)^(2*n+1);
r5 = k(5)*(P(6)*P(8) - P(5)*P(7)/keq2);

dPAdt = r1-2*r2-(4-m)*r3;
dPBdt = r2-(m-2)*r3;
dPCdt = r3;
dPDdt = r4;
dPEdt = r5;
dPFdt = r1+n*r4-r5;
dPGdt = -2*r1-(2*n+1)*r4+r5;
dPHdt = r2+r3+n*r4-r5;
dPdt = [dPAdt; dPBdt;dPCdt;dPDdt;dPEdt;dPFdt;dPGdt;dPHdt];

8楼: Originally posted by csgt0 at 2012-06-01 10:18:18
.*,.^不完全，行向量和列向量不分，请检查公式。