[理学]计算物理课件 Chapter3.ppt

第 3 章 曲线拟合的最小二乘法 3.1 拟合曲线 3.2 线性拟合和二次拟合函数 3.3 解矛盾方程组 3.4 程序示例 3.1 拟合曲线 已知一组数据(xi, yi), y = f(xi)，i = 1,2,…, m。 构造插值函数φ(x) 来逼近 f(x), 则有 φ(xi) = f(xi) = yi, i = 1,2,…, m。 (0) 或记 Q =(φ(x1) , φ(x2) ,…,φ(xm) ), Y = (y1, y2,…,ym), 则有 Q = Y. 拟合函数：如果数据不能同时满足某个特定函数，而要求所求的逼近函数“最优地”靠近数据点，即向量Q与Y 的误差或距离最小。按 Q 与Y 的误差最小原则作为最优标准所构造出的函数，我们称为拟合函数。 在多个景点之间修一条直线公路。 定义Q与Y 之间的距离： 其中，R 称为均方误差。 最小二乘法：按均方误差达到极小构造拟合曲线的方法。 确定拟合曲线的方法：（1）选择曲线类型；（2）若曲线类型难以确定，画散点图x；（3）用多种曲线类型拟合，选择最小二乘法意义下误差最  小的拟合曲线。 如在多个景点之间修一条主干道。 已知景点(xi, yi), i=1,2,…, m. 设 φ(x) = a +b x, 求 a，b使其均方误差达到最小。记δi = φ(xi)–yi， 3.2 线性拟合和二次拟合函数 ★ 线性拟合 给定数据(xi, yi), i=1,2,…, m. 求拟合曲线 p(x) = a +b x, 确定 a，b，使其均方误差 达到最小。由多元函数极值的必要条件，Q(a, b)达到极小，a, b满足 (3.3)成为拟合曲线的法方程(正规方程组)。 解之得 a，b。 代入 p(x) = a +b x, 即得所求的拟合曲线。 例3.1 给出21组数据，用线性函数拟合鱼的种类和鱼的数 量的关系，m = 21。 解 设 p(x) = a + b x, 经计算： 法方程： ★ 二次拟合函数 给定数据(xi, yi), i=1,2,…, m. 求二次拟合曲线 p(x) = a0 +a1 x+a2x2, 确定 a0 ，a1 ，a2，使其均方误差达到最小： 由极值的必要条件得 约定 ，得到法方程如下： 例3.2 用二次多项式函数拟合如下数据: (例3.3) 解 m =7。经计算有： 代入(3.5)得到法方程： 解得 a0= 0.66667, a1=-1.39286, a2=-0.13095, 所以 p(x)=0.66667-1.39286x-0.13095x2. 拟合曲线的均方误差： 3.3 解矛盾方程组 设 A=(aij)m×n (mn), x=(x1, x2,…, xn), b=(b1, b2,…, bn), 考虑 Ax=b，即 由代数知识知，若 秩(A,b) = 秩(A),则方程组有解； 若 秩(A,b) ≠ 秩(A),则方程组无解,此时称为矛盾方程组。 矛盾方程组的解，即最小二乘解，是指在均方误差极小意 义下的解，即 min‖Ax-b‖22。 我们已经知道，用直线 p(x) = a0+a1x 拟合给定数据:(xi, yi), i = 1,2,…,m, 把数据代入直线方程得： 考察方程组：ATAa = ATy，即 于是我们得到法方程：(3.3) 上述表明， ATAa = ATy 的解a = (a0, a1)T, 正是线性拟合中法方程的解，即为使 最小的解。 上述极小问题也可表示为： 而