现代科学与工程计算-孟大志-第一篇1-4.pptVIP

* ? 问 题 ? 分段低次插值虽然具有简单、收敛性、整体连续性及数值计算的稳定性等优点，但在节点处常有“尖点”出现，光滑性较差。特别是需要给出节点处的导数值，这在多数问题中是不实际的。如何在没有节点导数数据时也能达到上述目的？为此引入样条插值函数。 1.4 三次样条插值 图1.4.1 早期机翼下轮廓的放样 如图1.4.1所示，在早期的板材曲线切割时，常把富有弹性的细长木条（样条）固定在样点上，其它地方让其自由弯曲，然后画出长条的曲线称为样条曲线，由此启发设计整体连续光滑的样条插值函数。 1.4.1 三次样条插值函数 定义1.4.1 给定函数f(x)在 [a, b]上一组节点： 及节点上的函数值 ， 函数S(x)是满足下列条件的函数 的三次样条插值 (1) S(xi) = f(xi) = yi， (2) S(x)在子区间 上是三次多项式，记为 ； (3) S(x)在插值节点处连续，即 (4) 在插值节点处连续，即 ； ； (5) ， 。 要保证S(x)的存在唯一性，必须附加两个边界条件。例如，满足下列三种边界条件中的任意一个： (1) 固支边界条件： (2) 自然边界条件： (3) 周期边界条件： 上述三种边界条件都有它们的实际意义，从力学角度看，附加 边界条件相当于在细梁两端加上约束。工程中常用自然边界条 件求样条插值函数，这类插值函数称为自然样条函数，利用插 值条件和连续线性条件列出线性方程组并求解，是一种构造样 条的基本方法。 1. 三次样条插值函数的构造 (1) 构造思想: 通过构造含待定参数的分段三次Hermite插值多项式来构造三次样条插值函数。构造Hermite插值多项式需要知道被逼近函数f(x)的导数，而导数通常是不知道的。三次样条插值函数的构造则不需要知道f(x)的导数值，直接将其作为待定参数，利用各节点在连接处的光滑性与连续性条件，建立关系式来确定待定参数，从而构造插值多项式。 (2) 构造方法：设样条函数在插值节点处的一阶导数为 ，由 于三次样条插值本质上是分段三次插值多项式，故可通过分段Hermite插值来构造三次样条插值多项式 ，在每个子区间 上有 （1.4.1） 其中 是Hermite插值的基函数，这 里的形式与 Hermite插值的区别在于 中 是未知的，因 此必须设法求出 。 (3) 追赶法解 关系式: 令 ，对 求二阶导数得 在点 的值为 同理可以得到 在 上的二阶导数在点 的值为 由于 ，则 ，即 上式称为三次样条插值的 关系式。 合并同类项可得方程组 （1.4.2） 其中 ， ， 这是关于n +1个未知量 的 个线性方程组。如果 加上边界条件则可解，例如自然边界条件 （1.4.3） 其中， 。 合并(1.4.2)式和(1.4.3)式并用矩阵形式表示方 程组为 （1.4.4） 追赶法求解方程组(1.4.4)，由方程组第一式可得 ， ，其中， 。再将 代入方程组第二式整 理得 ，其中， 。 依次类推，可得递推关系 其中 令 ，则由递推公式 1.4.2三次样条插值的计算步骤 (1) 根据已知条件顺次计算方程组系数 ， 求出三次样条插值的 关系式； (2) 由给定边界条件，确定 ； (3) 用追赶法求解方程组(1.4.4)，先顺次计算出 ， 再使用逆序回代求出 ； (4) 将求出的 代入(1.4.1)式，求出 上的三次样条插值 函数 。 对给定的点 ，须首先确定 所在的子区间 ，然后按三次样条插值的计算步骤计算 ，进而求出 。 1.4.3 收敛性与误差估计 1．收敛性 对三次样条插值函数来说，当插值节点逐渐加密时，插值函数收敛于 ，且插值函数的导数也收敛于 。这种性质 要比多项式插值优越得多。 2. 误差估计 设被插函数 ，插值区间的划分为： 则当最大步长 时，对一切 ， 三次样条插值函数 具有误差估计 其中， 。 例4 给定插值条件 0 0 0 0 3 2 1 0 边界条件为： ，试求满足上述条件 的三次样条插值 函数的表达式。 解： ， ， 代入m 关系式(1.4.2) 得 将 代入上面方程，化简得 解得 ，利用公式 求得三次样条函数如下： 1.4.4 汽车刹车距离求解 使用三次样条插值来预测车辆的停止距离，若给定自然边界条件进行三次样条插值，则有 ， 由公式 可计算出 代入m 关系式(1.4.2)可得 由追赶法可求得 代入公式 求得三次样条函数如表1.4.1所示，插值图像如图1.4.2 所示。 表1.4.1 图1.4.2 刹车距离问题的三次样条插