数值计算方法教案插值方法.docxVIP

复习： 数值计算方法的含义 误差及误差限 误差与有效数字 数值计算中应注意的问题 一．插值的含义 问题提出：  第二章 插值方法 已知函数 y ? f ?x?在 n+1 个点 x 0 , x , 1 , x 上的函数值 y n 0 , y , , y 1 n ，求任意一点 x? 的函数 值 f ?x??。 说明：函数 y ? f ?x?可能是未知的；也可能是已知的，但它比较复杂，很难计算其函数值 f ?x??。 解决方法： 构造一个简单函数 P ?x?来替代未知（或复杂）函数 y ? f ?x?，则用 P ?x??作为函数值f ?x??的近似值。 二、泰勒（Taylor）插值1.问题提出： 已知复杂函数 y ? f ?x? 在 x 0  点的函数值 f ?x 0  ?， 求 x 0  附近另一点 x 0  h 的函数值 f ?x 0 h?。 解决方法： 构造一个代数多项式函数 P n  ?x?，使得 P n  ?x?与 f ?x?在 x ? x 0  点充分逼近。 泰勒多项式为：  f ?n? ? f ?n? ?x ? ? 0 n! P ?x?? f ?x n 0 ?? f ??x 0 ??x ? x ?? 0 0 ?x ? x 2! 0 ?2 ? ?x ? x ?n 0 显然， P n ?x?与 f ?x?在 x ? x 0  点，具有相同的 i 阶导数值（i=0,1,…,n）。 几何意义为： P ?x?与 f ?x?都过点?x n 0 , f ?x 0 ??； P ?x?与 f ?x?在点?x n 0 P ?x?与 f ?x?在点?x n 0 , f ?x 0 , f ?x 0 ??处的切线重合； ??处具有相同的凹凸性； 其几何意义可以由下图描述，显然函数 f 3?x?能相对较好地在 x 点逼近 f ?x?。 0 f2(x) f3(x) f1(x) f(x) x 0 误差分析（泰勒余项定理）： f ?n?1? ?? ??n ?1?!P ?x? f ?n?1? ?? ? ?n ?1?! n 0 ?n?1 ，其中? 在 x 0 与 x 之间。 举例： x已知函数 f ?x?? x  ，求 f ?115?。 分析：本题理解为，已知“复杂”函数 f ?x?? 在 x =100 点的函数值为 f ?x x0 0 x ?? 10 ， 求 x 的附近一点 x 0 0 解： +15 的函数值 f ?x 0 ?15?。 构造 1 次泰勒多项式函数 P ?x?： P ?x?? f ?x ?? f ??x ??x ? x ?。 1 1 0 0 0 2其中 f ?x ?? f ?100?? 10 ， f ??x?? 1 x?1 ， f ??x ?? f ??100?? 1 ，则有： 2 0 2 0 20 P ?x?? 5 ? 0.05x 1 故有 f ?115?? P ?115?? 10.75 1 误差分析： P ?115?? f 1  ?115? ?  f ??? ?2!?115 f ??? ? 2! 函数 f ??x?在[100,115]区间绝对值的极大值为 f ??100?? 2.5?10?4 ， 则有： P ?115?? f ?115? ? 0.028125 ? 0.05 1 于是近似值 10.75 有三位有效数字。 几何意义：显然， P ?x?也过点（100,10），且 P ?x?就是函数 f ?x??  x在点（100,10） x 1 处的切线，如下图所示。 构造 2 次泰勒多项式函数 P 2 1 ?x?： P ?x?? f ?x 2 0 ?? f ??x 0 ??x ? x ?? 0 f ??x 0 2! ??x ? x 0 ?2 。 把 f ?100?? 10 ， f ??100?? 1 及 f ??100?? 2.5?10?4 代入，有 20 f ?115?? P 2 ?115?? 10.721875 。 分析误差 f ??? ?3!P ?115?? f ?115? ? ?115 ?10 f ??? ? 3! 2 函数 f ??x?在[100,115]区间绝对值的极大值为 f ??100?? 3.75?10?6 ，则有 P ?115?? f ?115? ? 0.0021093750.005 2 于是近似值 10.721875 有四位有效数字。运行文件 taylor.m： %已知函数 f(x)=x^(1/2)，求 f(115) %一次泰勒插值subplot(1,2,1); f=inline(x^(1/2)); p1=inline(5+0.05*x); fplot(f,[-50,300]); hold on fplot(p1,[-