计算方法_第三章节.pptVIP

以 对式（2－7）两边作内积，得 将 的表达式，代入式（7）， 得 若记 则得递推公式 （10） 另外，误差得递推关系也不难得到。由式（6）可得 例3 利用正交多项式对例1中的数据作二次拟合。 解 按递推公式（8）和（9）计算权函数 为1的正交多项式 按公式（4）计算广义多项式系数 代入式（5）中，得到所求的二次函数 两边取对数，得 则得 令 两边取自然对数，得 令 则得 §3 非线性最小二乘拟合 （1） （2） （3） 两边取对数，得 则得 令 令 则得 （4） （5） 令 则得 （6） 令 则得 例4 给定实验数据 x 1.00 1.25 1.50 1.75 2.00 y 5.10 5.79 6.53 7.45 8.46 试求形如 的拟合函数。 解 对拟合函数的两边取自然对数，即 令 则上式 成为关于A，b 的线性函数 根据数据(x , y) 算出对应的(z , w) , 得下表 z 1.00 1.25 1.50 1.75 2.00 w 1.6292 1.7561 1.8764 2.0082 2.1353 建立法方程 解得 因此，所求的拟合函数为 §4 线性矛盾方程组的最小二乘法 方程个数大于未知量个数的方程组称为矛盾方程组，一般形式为 即 Ax＝b A是 n ×m阶的列满秩矩阵, x是 m维的列向量, b是 n维的列向量, 剩余向量 （11） （12） 因为 及A的m个列向量线性无关，易知 是m x m阶对称正定矩阵，而且上式右端最后一项是 正定二次型，同时其它两项与x无关。因此，欲使式 （12）成立，必须有 亦即 该式称为方程组Ax＝b 的法方程。因此，求解n阶矛盾 方程组的问题转化求解m阶线性方程组的问题。 如果考虑各个方程可靠性不同的因素，应对各方程 引进权系数 设W是由权系 数形成的对角矩阵，欲使 由 必有 其中W＝diag 这是考虑权系数方程组Ax＝b的法方程 例5 已知观测数据(1，-5)，(2，0)，(4，5)，(5，6)，试用最小二乘法求形如 上的经验公式。 解：记 按题意，得矛盾方程组， 写成矩阵形式，为 其中 其法方程为 即 解得 于是所求拟合曲线为 最小二乘法方曲线拟和是实验数据处理的常用方法。但当正规方程阶数较高时，往往出现病态。因此必须谨慎对待和加以巧妙处理。有效方法之一是引入正交多项式以改善其病态性。 最小二乘法方法评注 第三章习题 1 用最小二乘法求一个形如 的经验公式，使它与下列数据相拟合。 （1） 用普通多项式求解（2）用解矛盾方程组求解。 x 1 2 3 4 y 2 5 12 20 2、给定表中数据 试求形如 的拟合函数。 x 1 2 4 y 0.25 0.2 0.125 实例讲解 某种合成纤维的强度与其拉伸倍数有直接关系，下表是实际测定的24个纤维样品的强度与相应拉伸倍数的记录。 提示：将拉伸倍数作为x, 强度作为y,在座标纸上标出各点，可以发现什么? 第三章 数据拟合与最小二乘法 数据表格 已知的离散数据yi=f(xi) (i=0,1,2, …,n)往往是通过观测而得到的，经常带有观测误差。 曲线拟合：希望找到—条曲线，它既能反映结定数据的总体分布形式，又不致于出现局部较大的波动。这种逼近方式．只要所构造的逼近函数?(x)与被逼近函数f(x)在区间[a,b]上的偏差满足其种要求即可。 设给定数据点 (xi,yi), (i=0,1,2, …,n),记 并称ei为偏差。 曲线拟合的最小二乘法：以使得偏差的平方和 最小为标难 1.1 直线拟合 假设所给的数据点(xi，yi), (i=0,1,2, …,n)的分布 大致呈直线，故可选择线性函数作拟合曲线。 §1 线性最小二乘拟合 【问题1】 对于给定的数据点(xi，yi), (i=0,1,2, …,n) ， 求作一次式y＝a+bx，使总误差为最小。 【解】由微积分的知识可知，这一问题的求解，可归结为求二元函数E(a,b)的极值，即 这是关于a,b的线性力程组，称为法方程。 1.2 多项式拟合 有时所给数据点的分布并不一定近似地呈一条直线， 这时若仍用直线似合显然是不合适的。对于这种情 况，可以考虑用多项式拟合。 【问题2】 对于给定的数据点(xi，yi), (i=0,1,2, …,n) ， 求作多项式，