04_多元模型回归和分析.pptVIP

* * 使用MATLAB进行参数估计示例 例1：等温积分反应器的参数估计 在一等温积分反应器的动力学实验中，发生如下反应 ： 已知：当t=0时，CA0=1，CB0=0，CC0=0。实验数据如下表： ti 0 5 10 15 20 50 80 100 200 CAi,exp 1.000 0.840 0.679 0.503 0.425 0.137 0.030 0.006 0.000 CBi,exp 0.000 0.155 0.327 0.426 0.477 0.571 0.502 0.385 0.145 模型： 因为CA+CB+CC=1，所以描述反应器出口反应物浓度变化的模型方程只需要两个： 初始条件为：t=0， CA0=1，CB0=0，CC0=0。 * 使用MATLAB的dsolve求积分 s=dsolve(Dy1=-k1*y1, Dy2=k1*y1-k2*y2, y1(0)=1, y2(0)=0) s = y1: [1x1 sym] y2: [1x1 sym] s.y1 s.y2 simplify(s.y2) ans = exp(-k1*t) ans = -(-k1+k2)/(k1-k2)^2*k1*exp(-k2*t)-1/(k1-k2)*k1*exp(-k1*t) ans = k1*(exp(-