三次样条插值算法详解.pptVIP

演示文稿三次样条插值算法详解目前一页\总数五十八页\编于二十三点

（优选）三次样条插值算法详解目前二页\总数五十八页\编于二十三点

3样条函数的定义定义4.1设区间[a,b]上给定一个节点划分a=x0x1……xn-1xn=b 如果存在正整数k使得[a,b]上的分段函数s(x)满足如下两条：(1)在[a,b]上有直到k-1阶连续导数。(2)在每个小区间[xi,xi+1]上是次数不大于k的多项式。则称分段函数s(x)是以(2.6)为节点集的k次样条函数。目前三页\总数五十八页\编于二十三点

4三次样条插值函数的定义并且关于这个节点集的三次样条函数s(x)满足插值条件：则称这个三次样条函数s(x)为三次样条插值函数。目前四页\总数五十八页\编于二十三点

5三次样条插值函数的边界条件插值条件：连续性条件：一阶导数连续条件：二阶导数连续条件：目前五页\总数五十八页\编于二十三点

6（1）因为s(x)在每个小区间上是一个次小于三次的多项式，故有四个未知系数；（2）因为s(x)有n分段，从而共有4n个未知系数！（3）但插值条件与样条条件仅给出4n-2个条件，无法定出4n个未知系数，还差2个条件！这2个条件我们用边界条件给出！目前六页\总数五十八页\编于二十三点

7通常我们对插值多项式在两端点的状态加以要求也就是所谓的边界条件:第一边界条件：由区间端点处的一阶导数给出即 目前七页\总数五十八页\编于二十三点

8第二边界条件：由区间端点处的二阶导数给出即 特殊情况为自然边界条件：由区间端点处的二阶导数恒为0给出即 目前八页\总数五十八页\编于二十三点

9这样三次样条插值问题就分成三类！其实不止这三类!第三类又称周期边界条件：由区间端点处的函数值或导数值满足周期条件给出 目前九页\总数五十八页\编于二十三点

10样条函数的例子容易验证：是满足如下数据的第一类边界样条插值问题解：x0123y0000y’10目前十页\总数五十八页\编于二十三点

11样条函数的例子目前十一页\总数五十八页\编于二十三点

12通常有三转角法、三弯矩法、B样条基函数法。三次样条插值函数的求法这三种方法的基本思想是类似的，都是通过待定某些参数来确定插值函数，但肯定不是待定4n个参数。而是利用已知条件将待定参数减小到最少。比如：待定一阶导数、待定二阶导数、采用基函数方法来确定插值函数。目前十二页\总数五十八页\编于二十三点

13三转角法:待定一阶数为了确定三次样条插值函数的表达式S(x)，我们采用待定系数法来求解，我们待定什么系数呢？考虑到带一阶导数的分段三次Hermite插值多项式目前十三页\总数五十八页\编于二十三点

14我们采用待定一阶导数的方法即设因为分段三次Hermite插值多项式已经至少是一阶连续可导了，为了让它成为三次样条函数只需确定节点处的一阶导数使这些节点处的二阶导数连续即可！目前十四页\总数五十八页\编于二十三点

15目前十五页\总数五十八页\编于二十三点

16由于在内部节点处二阶导数连续条件：整理化简后得：目前十六页\总数五十八页\编于二十三点

17称为三转角法基本方程组以上推导还没有考虑边界条件！针对不同类型的三次样条问题，就可以导出不同的方程组！目前十七页\总数五十八页\编于二十三点

18第一类三次样条插值问题方程组基本方程组化为n-1阶方程组由于已知：化为矩阵形式目前十八页\总数五十八页\编于二十三点

19这是一个严格对角占优的三对角方程组，用追赶法可以求解！目前十九页\总数五十八页\编于二十三点

20第二类三次样条插值问题的方程组由于已知：故得：目前二十页\总数五十八页\编于二十三点

21稍加整理得联合基本方程组得一个n+1阶三对角方程组，化成矩阵形式为：仍然是严格对角占优目前二十一页\总数五十八页\编于二十三点

22第三类样条插值问题的方程组立即可得下式：由于：目前二十二页\总数五十八页\编于二十三点

23其中：联合基本方程得一个广义三对角或周期三对角方程组：这个方程组的系数矩阵仍然是严格对角占优阵！目前二十三页\总数五十八页\编于二十三点

24求解这些不同类型的样条插值问题的方程组，我们可得所要待定的一阶导数：称为三次样条插值问题三转角公式！再代入S(x)的每一段表达式，就求得三次样条函数的表达式!目前二十四页\总数五十八页\编于二十三点

25例1.对于给定的节点及函数值解:这是自然边界条件下的样条问题。目前二十五页\总数五十八页\编于二十三点

26我们可以将上述计算列于表中：k0123xk1245yk1342mk????Mk0??0hk121*λk*2/31/3*μk*1/32/3*gk69/2-7/2-6目前二十六页\总数五十八页\编于二十三点

27由些得如下方程组：利用三转角公式：目前二十七页\总数五十八页\编于二