计算方法课件第三章.ppt

* * * * * * * * * * */47 §3.2 什么是最小二乘法 常见做法： ?使 最小 /* minimax problem */ 偏差最大绝对值最小 使 最小 偏差绝对值之和 使 最小 /* Least-Squares method */ 偏差平方和最小 太复杂? 不可导，求解困难? */47 §3.2 什么是最小二乘法 三. 最小二乘原则 1. 最小二乘原则 使偏差平方和最小（上页中方法3）的原则称为最小二乘原则； 2. 最小二乘法 按最小二乘原则选择拟合曲线y=φ(x)的方法 P71 */47 §3.2 什么是最小二乘法 四. 最小二乘法 P72 */47 §3.2 什么是最小二乘法 */47 §3.3 最小二乘解的求法 本节内容 一. 最小二乘解的求法 二. 最小二乘拟合多项式的存在唯一性 三. 一般最小二乘拟合 返回章节目录 */47 §3.3 最小二乘解的求法 一. 最小二乘解的求法 */47 §3.3 最小二乘解的求法 P73 3.3.2 */47 §3.3 最小二乘解的求法 */47 §3.3 最小二乘解的求法 P74 m次多项式拟合 */47 §3.3 最小二乘解的求法 （1）直线拟合（一次多项式拟合） 若 ， a0, a1满足法方程组 即a0, a1是法方程组的解。 */47 §3.3 最小二乘解的求法 例1 已知一组试验数据 试用直线拟合这组数据. (计算过程保留3位小数)。 解 设直线y＝a0+a1x，那么a0, a1满足的法方程组公式为 */47 故法方程组为 解得a0=1.229 a1=1.483 所求直线方程为 y=1.229+1.483x §3.3 最小二乘解的求法 计算列表如下： */47 §3.3 最小二乘解的求法 （2）二次多项式拟合 若 满足法方程组 即a0, a1 , a2是法方程组的解。 */47 §3.3 最小二乘解的求法 例2 已知一组试验数据 试用二次多项式拟合这组数据. (计算过程保留2位小数) 解 设 满足的法方程组 */47 §3.3 最小二乘解的求法 计算列表如下： */47 §3.3 最小二乘解的求法 故法方程组为 解得a0=5.05 a1=-4.04 a2=1.01 所求二次多项式为 y=5.05－4.04x+1.01x2 */47 §3.3 最小二乘解的求法 二. 最小二乘拟合多项式的存在唯一性 定理：设点x0,x1,…xm互异，则法方程组 存在唯一的一组解 证明：用反证法（略） */47 §3.3 最小二乘解的求法 三. 一般最小二乘拟合 */47 §3.3 最小二乘解的求法 两种非线性最小二乘问题的求解途径 （1）化为线性最小二乘问题 部分可化为线性拟合问题的常见函数类型见下页表 （2）马奎特（Marquardt）方法 是在计算机上求解非线性最小二乘拟合问题的最为常用和有效的方法之一。（略） */47 §3.3 最小二乘解的求法 */47 §3.4 加权最小二乘法 本节内容 一. 加权最小二乘法 二. 例 返回章节目录 */47 §3.4 加权最小二乘法 一. 加权最小二乘法 重度：即权重或者密度，统称为权系数。它的大小反映了数据（xi,yi）地位的强弱。 定义加权平方误差为： */47 §3.4 加权最小二乘法 */47 §3.4 加权最小二乘法 二. 例 例3 已知一组实验数据（xi , yi）及权 Wi 如表所示。若 x 与 y 之间有线性关系y＝a + bx，试用最小二乘法确定系数 a 和 b。 P84 */47 §3.4 加权最小二乘法 解 因拟合曲线为一次多项式曲线（直线）φ1 (x) = a + bx，故相应的法方程组形式如式。将表中各已知数据代入即得法方程组 解之得 */47 小结 §3.1 引言