数值计算方法报告.pptVIP

代入(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, 正是线性拟合中法方程的解，即为使 最小的解。 上述极小问题也可表示为： 而 又是方程组 Aa = y 的最小二乘解。 定义 设 A=(aij)m×n, mn, 方程 ATAx = ATb 称为矛盾方程组  Ax = b 的法方程。 定理3.1 （1）若秩(A) = n，则法方程恒有解。 （2） x 是 的解，当且仅当x是法方程的解。 注：定理3.1告诉我们，求解拟合曲线的极小值问题与求解矛盾方程组的法方程等价。因此，求解拟合曲线的极小 问题可以转化为求解法方程：ATAx=ATb。 对于离散数据: (xi, yi), i =1,2,…,m, 用n次多项式来拟合曲线。设 求解极小问题 等价于求解矛盾方程组 Aa = y 的极小问题， 其中， 因此，求 a0, a1,…, an, 就是求解法方程： ATAa = ATy。 例3.3 用二次多项式函数拟合如下数据（见例3.2）。 解 m =7，经计算有： 设 p(x) = a0 + a1x + a2x2, 形成法方程： 约定 直接计算有： 于是得到法方程 解得 a0 = 0.66667, a1 = -1.39286,a2 = -0.13095, 所以 p(x)=0.66667-1.39286x-0.13095x2 ★ 非线性问题线性化 一般地，设 p(x) = a0 + a1x + a2x2 +…+ anxn. 作变换 将数据(xi,yi)由变换转化为(z(i),yi),再求解法方程：ATAa=ATy， 其中， 由法方程求出a0,a1,…,an, 再代回原变量x就求出了非线性曲线。 特别地，如果非线性函数为将给定数据(xi,yi)转换为(ui,vi),求出a,b,再代回原变量y, x，可求得原非线性拟合曲线。 例3.6 用给数据求经验公式：y = a ebx。 x 1 2 3 4 5 6 7 8 y 15.3 20.5 27.4 36.6 49.1 65.6 87.8 117.6 解 线性化。对经验公式取自然对数 ln y = ln a + bx （令u = ln y ，b0 = ln a ，则 u=b0+bx） 代入数据的矛盾方程组 由法方程 ATAB=ATy，B =（b0, b), 即 ∴ a = e2.4369 =11.4375. y =11.4375e0.2912x. 谢谢，请批评指正！ * 现代数值计算报告 2012年9月8日 主讲内容 误差及有效数字