计算方法插值方法-3-课堂教学版.pptVIP

可见，所谓分段插值，就是选取分段多项式作为插值函数。 * * S1(x)在每个子区间上表示为 在区间[a,b]上可表示为 Li(x)为基函数 分段线性插值基函数li(x)只在xi附近不为零，在其他地方均为零，这种性质称为局部非零性质。 * ab ≤(a+b)2/4 |ab|≤(|a|+|b|)2/4 * * 384=4！*42 * 是一种显示算法，其算法简单，收敛性能得到保证。只要节点间距充分小，分段插值法总能获得所要求的精度，而不会像高次插值那样出现龙格现象 如果修改某个数据，那么插值曲线仅仅在某个局部范围内受到影响，而代数插值却会影响到整个插值区间 分段线性插值没有导数，分段埃尔米特插值只有连续的一阶导数 这种插值要求给出各个节点上的导数值，所要提供的信息“太多” 前两条是优点，后两条是缺点。 改进这种插值以克服其缺点，就导致了所谓三次样条插值的提出。 分段低次插值的优点是具有收敛性与稳定性，缺点是光滑性较差，不能满足实际需要：例如高速飞机的机翼形线、船体放样形值线、精密机械加工等都要求有二阶光滑度，即二阶导数连续。 我们希望构造一种新的插值，既是分段低次插值，又能保证各段之间连接处的足够光滑性，既具有分段低次插值的一致性和计算的稳定性，又具有高次插值的整体的足够光滑性。这种双优兼备的插值，就是样条函数插值。 * 样条函数是美国数学家舒恩伯格(I.J.Schoenberg)于1946年首先提出的，60年代以来它的理论取得了重大发展和完善，并得到了广泛的应用，它是计算几何和计算机辅助设计的重要工具，今天在数值积分、数值微分、积分方程和微分方程等方面都有应用。本节主要介绍应用最广的三次样条插值函数，不涉及样条函数的其它深入理论和广泛应用问题。 * 几点说明： k次样条函数是k次多项式的推广。事实上，如果具有分划Δ的k次样条函数Sk(x)在每个节点xi(1≤i≤n-1)上具有连续的k阶导数，则它便退化为普通的k次多项式。 阶梯函数和折线函数分别是简单的零次样条函数和一次样条函数 特点：既是充分光滑的(Sk(x)直到k-1阶导数均连续)，又保留有一定的间断性(Sk(x)的k阶导数在节点处可能有间断)。光滑性保证了外形曲线的美观，间断性则使它能转折自如的被灵活运用。 今后样条函数常简称为样条 * * * 三次样条插值函数S3(x)是一个分段三次多项式，若用S3i(x)表示S3(x)在第i个子区间[xi,xi+1](i=0,1,…,n-1)上的表达式，则其形如 S3i(x)=ai+bix+cix2+dix3 其中ai、bi、ci、di为四个待定系数。由于子区间共有n个，所以要确定S3(x)需要确定4n个待定系数。 * 事实上，根据f(x)的周期性知y0=yn以及S3(xi)=yi，从而必有S3(x0+0)= S3(xn-0)。因此，只需要考虑后两个等式即可。 * 问题7理论上可以利用方程组(46)及边界条件求出4n个待定系数ai,bi,ci,di(i=0,1,…,n-1)，从而得到S3(x)在每个子区间[xi,xi+1]上的表达式S3i(x)。但是，由于这种做法的计算量相当大，在实际计算时很少用到。 假定了mi之后，就可以用mi构造每个子区间上的多项式，这个时候问题就转变为确定这n-1个未知的mi了。然后再加上内节点上二阶导数连续的n-1个条件即可确定mi，从而确定整个区间上的表达式。 * 三次样条插值是一种特殊的分段三次插值，选取节点上的导数值S3‘(x)=mi作为参数，参照分段三次插值公式，可列出样条函数的公式为： 可以看出不管参数mi怎么取值，这样构造出来的S3(x)在每个节点xi(1≤i≤n-1)上必定连续且有连续的一阶导数。现在的问题是，怎样选取参数mi的值使其二阶导数也连续。 * 为了保证二阶导数的连续性，要求 也就是上面这两个式子要相容 * * * 非零元素集中在三条对角线上而被称作是三对角型的。求解这类方程组的一种有效方法是追赶法。在后面会讲到具体的算法。 * * 我们看到，样条插值的逼近效果是令人满意的，它消除了前面所述的龙格现象。 * * 十二．样条函数 * 2．三次样条的定义 十二．样条函数 * 3．三次样条插值函数的边界条件 待定系数个数： 已知条件： 个方程 4n-2 补充条件：这两个条件通常在区间[a,b]的两个端点给出，称为边界条件 4n 十二．样条函数 * 3．三次样条插值函数的边界条件 3)当y=f(x)是周期为b-a的函数时，要求S3(x)及其导数也都是周期函数，相应的边界条件为 作为特例，当 ， 时称为自然边界条件。满足自然边界条件的样条函数称为自然样条函数。 1)给定一阶导数值 2)