一、前言#
1、计算方法#
最小二乘法#
计算题【3-3】在甲苯中不同温度中测定偶氮二异丁睛的分解速率常数,数据如下求分解活化能。再求40℃和80℃下的半衰期,判断在这两个温度下聚合是否有效。
温度/℃ | 50 | 60.5 | 69.5 |
---|
分解速率常数/s−1 | 2.64×10-6 | 1.16×10-5 | 3.78×10-5 |
kd=Ae−RTE
lnkd=lnA−RTE ,以lnkd-T1 作图,斜率为-RE ,截距为 lnA
一、最小二乘法的使用方法(简明版)#
适用场景:已知一组数据点 (x1,y1),(x2,y2),…,(xn,yn),需拟合直线方程 y=a+bx。
步骤:
-
确定变量:
- 自变量 x(如示例中的 1/T)
- 因变量 y(如示例中的 lnkd)
-
计算求和项:
- ∑x:所有 xi 的和
- ∑y:所有 yi 的和
- ∑xy:所有 xiyi 的和
- ∑x2:所有 xi2 的和
-
代入公式求斜率和截距:
- 斜率 b:
b=n∑x2−(∑x)2n∑xy−∑x∑y
- 截距 a:
a=n∑y−b∑x
-
写出方程:
y=a+bx
二、直观解释#
最小二乘法的核心思想是:通过最小化预测值与真实值的误差平方和,找到最佳拟合直线。
-
定义误差平方和:
对于每个数据点,预测值为 y^i=a+bxi,真实值为 yi,则总误差为:
S=i=1∑n(yi−y^i)2=i=1∑n(yi−a−bxi)2
-
最小化误差 S:
- 对 a 和 b 分别求偏导,令导数为零,找到极值点。
- 对 a 求导:
∂a∂S=−2∑(yi−a−bxi)=0
化简得:
∑yi=na+b∑xi
- 对 b 求导:
∂b∂S=−2∑xi(yi−a−bxi)=0
化简得:
∑xiyi=a∑xi+b∑xi2
-
联立方程求解 a 和 b:
上述两式组成方程组:
{na+b∑xi=∑yia∑xi+b∑xi2=∑xiyi
解此方程组即得到斜率和截距的公式。