我有一个约束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将保持恒定。
我想知道如何与时间t的值变化K3值的变化，对于本我想最初估计K3在t间隔t = 0：0.1：60 ..
任何帮助非常感谢
