计算机数学基础课件教学课件作者刘树利09课件.pptVIP

第九章数值计算初步 第九章　数值计算初步 后页 首页 前页 后页 首页 前页 基本要求、重点难点 9.1　数值计算中的误差 9.2　函数插值法 9.3　方程f(x)＝0的数值解法 9.4　数值积分 9.5　常微分方程的数值解法 9.6　演示与实验八 ? 基本要求 　　?掌握数值计算方法。 　　?了解函数插值法及数值积分。 　　?掌握方程f(x)＝0的数值解法及常微分方程的数值解法。 ? 重点难点 重点： ?无穷个数的问题——级数。 ?级数在各种领域的应用及幂级数。 难点： ?方程f(x)＝0的数值解法及常微分方程的数值解法。 ? 9.1　数值计算中的误差 　　1. 描述误差：在将实际问题转化为数学模型时，为了便于分析和计算，进 行合理简化，往往要忽略一些次要因素(如求物体下落速度时往往忽略空气阻力 等)，这样产生的误差为描述误差。　　2. 测量误差：数学模型中的已知参数，多数是通过测量得到的。而测量过 程受工具、方法、观测者的主观因素、不可预料的随机干扰等影响，得到的数 据与实际数据有一定的误差，称此类误差为测量误差。　　3. 舍入误差：计算机只能处理有限数位的小数运算，初始参数或中间结果 都必须进行四舍五入运算，这就必然产生舍入误差。　　4. 截断误差：当数学模型不能得到精确解时，通常用数值方法来求它的近 似解，其近似解与精确解之间的误差称为截断误差。 ? 9.2　函数插值法 ?9.2.1有关概念 定义9.1 　　设y＝f(x)在［a，b］上有定义，且已知在点a≤x0＜x1＜…＜xn≤b上的值y0，y1，y2，…，yn，若存在一简单函数P(x)使：P(xi)＝f(xi)成立，就称P(x)为f(x)的插值函数，点x0，x1，…，xn称为插值节点，包括插值节点的区间［a，b］称为插值区间。求插值函数P(x)的方法称为插值法。 　　从几何意义上：插值法就是求曲线y＝P(x)使其通过给定的n+1个点(xi，yi) (i＝0，1，…，n)，并用它近似已知曲线y＝f(x)，如图所示。 定理9.1 (魏尔斯特拉斯(Weierstrass)多项式逼近定理) 　　设y＝f(x)是闭区间［a，b］上的任意连续函数，ε＞0为任意正数，则存在多项式P(x)，使对一切x∈［a，b］有 ｜f(x)－P(x)｜≤ε。 ?9.2.2拉格朗日插值多项式 　　所谓插值法就是找一个多项式函数P(x)使在插值点xi (i＝0，1，2，…，n)上使P(xi )＝f(xi ) (i＝0，1，2，…，n)。 　　拉格朗日插值多项式: 定理9.2 　　满足P(xi )＝f(xi ) (i＝0，1，2，…，n)，且次数不高于n的多项式P(x) 是惟一的，即 　　　　　　　　　　　　　　　　　式为惟一满足此要求的多项式，且若f(x)在［a，b］上存在n+1阶导数，则存在ξ∈(a，b)，使 f(x)－P(x)＝f(n+1)(ξ)(n+1)!(x－x0)(x－x1)…(x－xn)。 ? 9.3　方程f(x)＝0的数值解法 ?9.3.1根的初始区间 　　例　考察f(x)＝x3－x－1＝0。　　注意到f(0)＝－1＜0，f(2)＝5＞0，知f(x)在区间(0，2)内至少有一个实根。 　　从上表可以看出区间［1.0，1.5］上，f(1.0)＜0，f(1.5)＞0，知区间(1.0，1.5)内必有一根，从而把f(x)在(0，2)上搜索根，很容易地缩小为(1.0，1.5)上搜索根。可见寻求方程根的初始区间很重要。　　用扫描法确定方程的根的初始区间，可给定一个扫描下界，如x＝a，然后按预先选定步长h，依次计算 f(a)，f(a+h)，f(a+2h)，…，f(a+nh)， 其中b＝a+nh是扫描区间的上界，同时检验每次函数值 f(a+kh)与f(a+(k－1)h) 是否异号，如果是异号，表示方程f(x)＝0在区间［a+(k－1)h，a+kh］内有根，即为一个初始有根区间。 ?9.3.2二分法 　　二分法：用二分有根区间求方程近似根的方法。　　二分法的算法步骤如下： ?9.3.3牛顿迭代法 　　牛顿(迭代)公式： xn+1＝xn－f(xn)/f ′(xn)。 　　牛顿迭代法的算法步骤如下： ? 9.4　数值积分 ?9.4.1矩形求积公式 ?9.4.2梯形求积公式 ?9.4.3抛物线法求积公式 ? 9.5　常微分方程的数值解法 ?9.5.1欧拉公式 　　在式y(xi+1)＝yi+h·f(xi+θh，y(xi+θh))中，取θ＝0，即取 f(xi+θh，y(xi+θh))≈f(xi，yi)， 于是有 y(xi+1)≈yi+h·f(xi，yi)， 记右端值为yi+1，则有 y(xi+1)≈yi+1＝y