数值4报告范本.pptVIP

} 例 给定区间[0,3]上3个点的函数值f(0)=0,f(1)=2,f(3)=4, 试求数 a,b,c,d,使函数 S(x)为给定点上的三次样条插值函数。 解 设 根据定义，由 得 d=0，故 则 由 得 由 得 由 得 由 得 求得 } 给定n+1个样点(xi , yi )(i=0,1,…,n) ，确定一个三次样条插值函数需要4n个独立条件。在定义中，已指定了4n-2个条件，即 所以，一般需补充指定2个边界条件。 } 4.4.2 三次样条插值函数的求法 三弯矩法：以各基点处的二阶导数作为待定参数。 在区间[xk , xk+1]上，对二阶导函数S″(x)利用拉格朗日线性插值公式得 由泰勒公式有 其中 用x=xk+1 代入上式可解得 令 } 得到 S(x) 在[xk , xk+1 ]上的表达式(记为Sk(x)) S(x) 在[xk-1 , xk ]上的表达式为 根据 经过一些推导整理可得到 其中 } 这是含有n+1个未知量m0 , m1 ,…,mn 的n-1个方程的方程组，要确定未知量还需补充两个边界条件。 常用的边界条件有三种: (1)第一种边界条件：给定端点处的一阶导数，即 (2)第二种边界条件：给定端点处的二阶导数m0 ,mn。特别取m0 = mn=0时称为自然边界条件，求得的S(x)称为自然样条函数。 (3)第三种边界条件：当y=f(x)是以b- a为周期的周期函数时，端点要满足 } 例如，当给定第二种边界条件时，方程组 化为 } 求三次样条插值函数的计算步骤归纳为: (1)计算 (2) 结合给定的边界条件得出下列方程组并求解： (3)将m0 , m1 ,…,mn代入 S(x) 的分段表达式： } 例 给定函数表，求自然三次样条插值函数，并求f(3)。 1 3 4 2 y 1 2 4 5 x 解 已知n=3, m0=m3=0. 求m1 , m2 的方程组形为 方程组化为 解得 * 第4章 插值法和曲线拟合 4.1 插值法的基本理论 4.2 拉格朗日插值多项式 4.3 牛顿均差插值多项式 4.4 三次样条插值 4.5 曲线拟合的最小二乘法 } 4.1.1 插值问题及代数多项式插值 函数y=f(x)给出一组函数值 x: x0 x1 x2 …… xn y: y0 y1 y2 …… yn 其中x0 ,x1,x2 ,…,xn是区间[a,b]上的互异点，要构造一个简单的函数φ(x) 作为f(x)的近似表达式，使满足 (插值原则、插值条件 ) 这类问题称为插值问题。 -----f(x)的插值函数, f(x) -----被插值函数 x0 ,x1,x2 ,…,xn -----插值基点 求插值函数的方法称为插值法。 若x∈[a,b]，需要计算f(x) 的近似值φ(x),则称x为插值点。 } } 当选择代数多项式作为插值函数时，称为代数多项式插值问题: 代数多项式插值问题: 设函数y=f(x)在[a,b]有定义, 且已知在n+1个点a≤x0x1……xn≤b上的函数值y0, y1,……,yn.,要求一个次数不高于n的多项式 使满足插值原则 称Pn(x)为f(x)的n次插值多项式。 这样的插值多项式是否存在、唯一？ } 定理4.1 在n+1个互异基点处满足插值原则且次数不超过n的多项式Pn(x)是存在并唯一的。 证 其系数行列式 因此方程组存在唯一的解 ，因此Pn(x)存在并唯一。 } 4.1.2 插值多项式的误差 截断误差： 也称为Pn(x)的余项。 定理4.2 设函数f(x)在包含基点x0 , x1 ,…, xn的区间[a,b]上连续，在(a,b)上具有n+1阶导数，Pn(x)为满足插值原则的n次插值多项式，则对任一点x∈[a,b]，总存在相应的点ξ∈(a,b)，使 其中 证 当x=xi (i=0,1,…,n) 时，由插值条件知Rn(xi)=0，故结论显然成立。 当x是[a,b]上任一个固定点，但不是插值基点时， } 作辅助函数 根据罗尔定理，在F(t)的两个零点之间至少有一点使导数为0，即 F(t)在区间[a,b]上至少有n+2个互异的零点x, x0 , x1 ,…, xn。 F(n+1)(ξ) = 0 从而结论成立。 依此类推，可知F(n+1)(t)在(a,b)内至少有一个零点ξ，则 F??(t) 在(a,b)内至少有n个互异零点 F?(t)在(a,b)内至少有n+1个互异的零点 } } 推论 当f(x)是次