Hans J. Han

Back

一、前言#

1、计算方法#

最小二乘法#

计算题【3-3】在甲苯中不同温度中测定偶氮二异丁睛的分解速率常数,数据如下求分解活化能。再求40℃和80℃下的半衰期,判断在这两个温度下聚合是否有效。

温度/℃5060.569.5
分解速率常数/s1s^{-1}2.64×10-61.16×10-53.78×10-5

kd=AeERTk_d=Ae^{-\frac{E}{RT}}

lnkd=lnAERT\ln k_d=\ln A-\frac{E}{RT} ,以lnkdlnk_d-1T\frac{1}{T} 作图,斜率为-ER\frac{E}{R} ,截距为 lnAlnA

一、最小二乘法的使用方法(简明版)#

适用场景:已知一组数据点 (x1,y1),(x2,y2),,(xn,yn)(x_1, y_1), (x_2, y_2), \dots, (x_n, y_n),需拟合直线方程 y=a+bxy = a + bx

步骤

  1. 确定变量

    • 自变量 xx(如示例中的 1/T1/T
    • 因变量 yy(如示例中的 lnkd\ln k_d
  2. 计算求和项

    • x\sum x:所有 xix_i 的和
    • y\sum y:所有 yiy_i 的和
    • xy\sum xy:所有 xiyix_i y_i 的和
    • x2\sum x^2:所有 xi2x_i^2 的和
  3. 代入公式求斜率和截距

    • 斜率 bbb=nxyxynx2(x)2b = \frac{n \sum xy - \sum x \sum y}{n \sum x^2 - (\sum x)^2}
    • 截距 aaa=ybxna = \frac{\sum y - b \sum x}{n}
  4. 写出方程

    y=a+bxy = a + bx

二、直观解释#

最小二乘法的核心思想是:通过最小化预测值与真实值的误差平方和,找到最佳拟合直线

  1. 定义误差平方和
    对于每个数据点,预测值为 y^i=a+bxi\hat{y}_i = a + bx_i,真实值为 yiy_i,则总误差为:

    S=i=1n(yiy^i)2=i=1n(yiabxi)2S = \sum_{i=1}^n (y_i - \hat{y}_i)^2 = \sum_{i=1}^n (y_i - a - bx_i)^2
  2. 最小化误差 SS

    • aabb 分别求偏导,令导数为零,找到极值点。
    • aa 求导Sa=2(yiabxi)=0\frac{\partial S}{\partial a} = -2 \sum (y_i - a - bx_i) = 0 化简得: yi=na+bxi\sum y_i = na + b \sum x_i
    • bb 求导Sb=2xi(yiabxi)=0\frac{\partial S}{\partial b} = -2 \sum x_i(y_i - a - bx_i) = 0 化简得: xiyi=axi+bxi2\sum x_i y_i = a \sum x_i + b \sum x_i^2
  3. 联立方程求解 aabb
    上述两式组成方程组:

    {na+bxi=yiaxi+bxi2=xiyi\begin{cases} na + b \sum x_i = \sum y_i \\ a \sum x_i + b \sum x_i^2 = \sum x_i y_i \end{cases}

    解此方程组即得到斜率和截距的公式。

简单叙述 数学方法 学习中
https://hansjack.com/blog/describe-mathematical-method
Author Hans J. Han
Published at April 5, 2025
Comment seems to stuck. Try to refresh?✨