1

我有一个带有约束 K3 的化学动力学模型（2-组织隔室模型）（其中 K3 是速率常数）

我已经对等离子体函数和化学动力学模型进行了建模，以绘制输出特性

我想从下面的代码中估计速率常数 k3

function c_t = output_function_constrainedK3(t, a1, a2, a3,b1,b2,b3,td, tmax,k1,k2,k3)

DV_free= k1/(k2+k3);
K3 = k3*((k1/k2)/DV_free);
K_1   = (k1*k2)/(k2+K3);
K_2   = (k1*K3)/(k2+K3);

c_t = zeros(size(t));

ind = (t > td) & (t < tmax);

c_t(ind)= conv(((t(ind) - td) ./ (tmax - td) * (a1 + a2 + a3)),(K_1*exp(-(k2+K3)*t(ind)+K_2)),'same');

ind = (t >= tmax);

c_t(ind)= conv((a1 * exp(-b1 * (t(ind) - tmax))+ a2 * exp(-b2 * (t(ind) - tmax))) + a3 * exp(-b3 * (t(ind) - tmax)),(K_1*exp(-(k2+K3)*t(ind)+K_2)),'same');

plot(t,c_t);
figure

%plot(t,c_tnp);
axis([0 50 -2000 80000]);
xlabel time ;
ylabel concentration ;

end

下面列出了所有参数的初始估计值

t=0:0.1:60;
td =0.3;
tmax=0.8;
a1=2501;
a2=18500;
a3=65000;
b1=0.5;
b2=0.7;
b3=0.3;
k1=0.014;
k2=0.051;
k3=0.07;

请建议我一种从上述非线性动力学模型代码中估计 K3 参数的方法

在上述函数中，a1、a2、a3、b1、b2、b3、td、tmax、k1、k2 的参数值将保持不变。

我想知道 K3 值如何随着时间 t 值的变化而变化，为此我想最初以 t 间隔 t=0:0.1:60 估计 K3..

非常感谢任何帮助

4

0 回答 0