高校(理工类)数学求积公式的误差教学(课堂讲义).pptVIP

§4.4 求积公式的误差 当f(x)是次数不超过n次的多项式时，f(x)的插值多项式pn(x)就是它本身，这时，如果不计舍入误差，则公式(1)将准确地成立，即有 特别地，令f(x)=1知： 求积公式的误差 如果数据f(xk)的舍入误差均不超过ε，则根据第1章第2节的讨论，积分值： 的舍入误差限为： 而在求积系数λk全为正数的情况下，上述误差限等于： 可见，只要数据f(xk)有足够多的有效数字，就可以控制积分值I的舍入误差充分小。由此得知，舍入误差对数值积分的影响不像数值微分那样显著。 求积公式的截断误差 分析求积公式的截断误差。 先考察梯形法则。假定f(x)的二阶导数在(a,b)上改变不大，即设f(x)近似地取某个定值C2。将f(x)在x=a处泰勒展开，有： 再对上式两端在(a,b)上求积分，得： （4.4.1） 截断误差 另一方面，注意到 代入梯形公式(2)，知 式(4.4.1)与(4.4.2)的前两项相同，相减即得 在长为h=(b–a)/n的每个子区间(xk-1,xk)上用梯形公式计算积分值，其误差按(4.4.3)式近似地取定值–C2h3/12，因此将它乘上n倍即得复化梯形公式(5)的余项： （4.4.2） （4.4.3） 辛卜生公式的误差 继续考察辛卜生公式。假定f(4)(x)在(a,b)上近似地取定值C4，将f(x)在(a,b)的中点c=(a+b)/2展开