m教材《数值分析》陈晓江主编武汉理工大学出版社Y参考书目《数值分析》李庆扬、王能超、易大义编清华大学出版社施普林格出版社《科学和工程计算基础》施妙根、顾丽珍编清华大学出版社上页下页返园
上页 下页 返回 教材 《数值分析》 陈晓江主编 武汉理工大学出版社 参考书目 《数值分析》 李庆扬、王能超、易大义编 清华大学出版社 施普林格出版社 《科学和工程计算基础》 施妙根、顾丽珍编 清华大学出版社
第二章插值法第一节问题的提出第二节拉格朗日插值第三节牛顿插值第四节埃尔米特插值第五节分段低次插值第六节三次样条插值上页下页返园
上页 下页 返回 第二章 插值法 第一节 问题的提出 第二节 拉格朗日插值 第三节 牛顿插值 第四节 埃尔米特插值 第五节 分段低次插值 第六节 三次样条插值
86三次样条插值样条(spline),是早期飞机、造船工业中绘图员用来画光滑曲线的细木条或细金属丝,绘图时,为了将一些已知的点连成光滑曲线,绘图员用压铁把样条固定在相邻的若干点上,样条具有弹性,形成通过这些点的光滑曲线,沿着它就可画出所需的曲线.数学上仿此得出的函数便称为样条函数,由于分段线性插值和分段埃尔米特插值的光滑性不够(例如船体放样、飞机的外形曲线设计常需二阶可导),况且节点处的导数值难以获得。样条插值是用分段低次多项式去逼近被逼近函数,并且能满足对光滑性的要求,又无需给出每个节点处的导数值.它除了要求给出各个节点处的函数值之外,只需提供两个边界节点处的导数上页下页信息返园
上页 下页 返回 §6 三次样条插值 由于分段线性插值和分段埃尔米特插值的光滑性不够(例 如船体放样、飞机的外形曲线设计常需二阶可导),况且节点 处的导数值难以获得
一、样条函数设有节点a=x<x,<<x,=b,分段函数s(x)满足:(1)于每个区间x,xk+il上是一个次数不超过m的实系数代数多项式。(2)s(x)于[a,b上具有m-1阶连续导数。则称s(x)为[a,b上的m次样条函数。对m=3,便得到三次样条函数(1)S(x)在每个小区间x,Xk+il上都是三次多项式(2) S"(x)在区间a,b]上连续,即S(x)C[a,b]则称S(x)为区间a,bl上的三次样条函数此时曲线S(x)处处均有连续曲率(阶导数连续)上页一工程必须的。下页返圆
上页 下页 返回 设有节点a x0 x1 xn b,分段函数s(x)满足: 一、样条函数 2 ( ) [ , ] , ( ) [ , ] 2 ( )S x 在区间a b 上连续 即S x C a b (1)S(x)在每个小区间[xk , xk1 ]上都是三次多项式 则称 S(x)为区间[a,b]上的三次样条函数 工程必须的。 此时曲线 处处均有连续曲率(二阶导数连续) S(x) 系数代数多项式。 (1)于每个区间[xk , xk1 ]上是一个次数不超过m的 实 (2)s(x)于[a,b]上具有m 1阶连续导数。 则称s(x)为[a,b]上的m次样条函数。 对m 3,便得到三次样条函数:
如果函数f(x)在节点xo,xi,,x,处的函数值为f(x,)=yj,j=0,l,...,n而三次样条函数S(x)满足(1)S(x,)=yj,j=0,l,..",n则称S(x)为f(x)在[a,b上的三次样条插值函数注:三次样条与分段埃尔米特插值的根本区别在于S(x)自身光滑,不需要知道f的导数值(除了在2个端点可能需要);而埃尔米特插值依赖于f在所有插值点的导数值。H(x)f(x)S(x)上页下页返园
上页 下页 返回 注:三次样条与分段埃尔米特插值的根本区别在于S(x)自 身光滑,不需要知道f 的导数值(除了在2个端点可能需 要);而埃尔米特插值依赖于f 在所有插值点的导数值。 f(x) H(x) S(x) 如果函数f (x)在节点x0 , x1 , , xn 处的函数值为 f (xj ) yj , j 0,1, ,n 而三次样条函数S(x)满足 S(xj ) yj , j 0,1, ,n 则称S(x)为f (x)在[a,b]上的三次样条插值函数 -(1)