数值分析课件第六章数值插值方法.pptVIP

(2)分段二次Lagrange插值的公式为第66页,共97页，星期六，2024年，5月第67页,共97页，星期六，2024年，5月样条：是指飞机或轮船等的制造过程中为描绘出光滑的外形曲线(放样)所用的工具。样条本质上是一段一段的三次多项式拼合而成的曲线，在拼接处，不仅函数是连续的，且一阶和二阶导数也是连续的。1946年，Schoenberg将样条引入数学，即所谓的样条函数。因分段线性插值导数不连续，为克服在节点处不光滑的缺陷，故引入样条插值概念。三、三次样条插值第68页,共97页，星期六，2024年，5月一、三次样条插值函数的定义定义：给定区间[a,b]上的一个划分： a=x0x1…xn=b，已知函数f(x)在点xj上的函数值为 f(xj)=yj,(j=0,1,2,···,n)如果存在分段函数第69页,共97页，星期六，2024年，5月(1)S(x)在每一个子区间[xj-1,xj](j=0,1,2,···,n)上是一个三次多项式；(2)S(x)在每一个内接点xj(j=1,2,···,n-1)上具有直到二阶的连续导数；则称S(x)为节点x0,x1,…,xn上的三次样条函数。若S(x)在节点x0,x1,…,xn上还满足插值条件：(3)S(xj)=yj(j=0,1,2,···,n)则称S(x)为三次样条插值函数。（即全部通过样点的二阶连续可微的分段三次多项式函数）满足下述条件：第70页,共97页，星期六，2024年，5月三次样条插值多项式的确定：由(1)知，S(x)在每一个小区间[xj-1,xj]上是一三次多项式，若记为Sj(x)，则可设要确定函数S(x)的表达式，须确定4n个未知系数{aj，bj，cj，dj}(j=1,2,…,n)。由(2)知，S(x)，S`(x)，S``(x)在内节点x1,x2,…,xn-1上连续，则j=1,2,…,n-1第71页,共97页，星期六，2024年，5月可得3n-3个方程，又由条件(3)j=0,1,…,n得n+1个方程，共可得4n-2个方程。要确定4n个未知数，还差两个方程。通常在端点x0=a,xn=b处各附加一个条件，称边界条件，常见有两种：(1)二阶边界条件：(2)一阶边界条件：共4n个方程，可唯一地确定4n个未知数。第72页,共97页，星期六，2024年，5月例已知f(x)：f(-1)=1,f(0)=0,f(1)=1，求f(x)在[-1,1]上的三次自然样条插值函数，f(x)在-1,1处的二阶导数为0。解设由插值条件和函数连续条件得：第73页,共97页，星期六，2024年，5月由一阶及二阶导数连续得：由自然边界条件得：联立上面8个方程，求解得故第74页,共97页，星期六，2024年，5月用二阶导数值构造三次样条插值函数（又称三弯矩算法）记hj=xj-xj－1，(j=1,2,…,n)设S``(xj)=Mj，(j=0,1,2,…,n)因为S(x)是三次多项式，故S``(x)是线性函数。按线性插值公式可得积分两次得二、三次样条插值函数的建立第75页,共97页，星期六，2024年，5月其中c1，c2为积分常数。将S(xj-1)=yj-1，S(xj)=yj代入上式，可确定c1，c2。故第76页,共97页，星期六，2024年，5月上式称M表达式，只需确定Mi，即可确定三次样条插值函数。将M表达式两端对x求导，得令x=xj，得左导数令x=xj-1，得右导数第77页,共97页，星期六，2024年，5月故因S(x)一阶导数连续，即整理得第78页,共97页，星期六，2024年，5月记则得Mj的n-1个方程：(j=1,2,…,n-1)对固定边界条件：第79页,共97页，星期六，2024年，5月得或记为其中第80页,共97页，星期六，2024年，5月则求Mi的方程的矩阵形式为上述方程组是关于Mj(j=0,1,…,n)的三对角方程组。该方程组是严格对角占优的三对角方程组，有唯一解，可用追赶法求解。第81页,共97页，星期六，2024年，5月上述方程组是关于Mj(j=1,…,n-1)的三对角方程组。Mj在力学上解释为细梁在xj截面处的弯矩，称为S(x)的矩，故称三弯矩方程组。该方程组是严格对角占优的三对角方程组，有唯一解，可用追赶法求解。对于二阶边界条件：第82页,共97页，星期六，2024年，5月6.5埃尔米特插值（Hermite）Newton插值和Lagrange插值虽然构造比较简单，但都存在插值