hxyonline
2012-04-16, 17:56
Ln[B/H(x)]=ln[(0.0048AEa)/G(a)*8.314]-1.0516[Ea/RT] （1）
G(a)函数形式未知，所以只能通过求取多个点斜率进行求解Ea
H（x）随着x(x=Ea/RT)的变化而变化，它们分别为式1进行迭代的定义函数，其定义式为 （2）
H(x)=[exp(-x)h(x)/x^2]/0.0048exp(-1.0516x) （2）
迭代过程分为3步：（1）给一初值Ea1（这个值由其他公式求算），4个T值（表示为T1,T2,T3,T4）。
（2）由Ea1和T可以求出H(x)1, H(x)2, H(x)3, H(x)4。
(3)将H(x)1, H(x)2, H(x)3, H(x)4的值代入式1或式2中，并以 对1/T进行线性回归，通过斜率可求的新的Ea2值。
（4）以Ea2代替Ea1重复步骤2,3直至Ea(i)-Ea(i-1)<0.01时，即可认为Ea值为真实值。
我是新手，不知道写的大家能不能看懂，请大家指正。以前没学过matlab，纯新手，数据分析，比较着急，先谢谢大家了。
G(a)函数形式未知，所以只能通过求取多个点斜率进行求解Ea
H（x）随着x(x=Ea/RT)的变化而变化，它们分别为式1进行迭代的定义函数，其定义式为 （2）
H(x)=[exp(-x)h(x)/x^2]/0.0048exp(-1.0516x) （2）
迭代过程分为3步：（1）给一初值Ea1（这个值由其他公式求算），4个T值（表示为T1,T2,T3,T4）。
（2）由Ea1和T可以求出H(x)1, H(x)2, H(x)3, H(x)4。
(3)将H(x)1, H(x)2, H(x)3, H(x)4的值代入式1或式2中，并以 对1/T进行线性回归，通过斜率可求的新的Ea2值。
（4）以Ea2代替Ea1重复步骤2,3直至Ea(i)-Ea(i-1)<0.01时，即可认为Ea值为真实值。
我是新手，不知道写的大家能不能看懂，请大家指正。以前没学过matlab，纯新手，数据分析，比较着急，先谢谢大家了。