- 1、本文档共26页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
一、三次样条的产生和背景 第二章 插值法 §5 三次样条插值 二、三次样条函数的定义 三、三转角方程 四、三弯矩方程 预备知识 已知:4个条件 xi xk xk+1 yi = f(xi) yk yk+1 求:一个次数不超过3的多项式H3(x) Hermite插值: 结论: 其中 实际中有许多计算问题对插值函数的光滑性有较高的要求,例如飞机机翼外形、发动机进、排气口都要求有连续的二阶导数。 一、 三次样条的产生和背景 1.问题的产生 显然我们前面介绍的方法已不能解决这个问题。 2.样条的概念(Spline) 样条是工程设计中使用的一种绘图工具,它是富有弹性的细木条或细金属条。绘图员利用它把一些已知的点连接成一条光滑曲线称为样条曲线,样条曲线在连接点处有连续的曲率(即连续的二阶导数),它实际上是分段三次曲线拼接而成,在连接点上要求二阶导数连续。 二、 三次样条函数的定义 若函数S(x)∈C2[a,b],且在每个小区间[xj,xj+1]上是三次多项式,其中 a =x0<x1 <… <xn=b 是给定节点,则称S(x)是节点x0,x1, … ,xn上的三次样条函数。 1.三次样条的定义 a.S(x)∈C2[a,b] b.S(x)在[xj,xj+1]上是三次多项式 即: 三次样条函数 2.三次样条插值函数的定义 + S(xi) = yi 3.求解三次样条插值函数的已知条件数和未知条件数 未知参数个数 4n 已知条件个数 插值条件: n+1 S(x)∈C2[a,b] :3(n-1) 共 计: 4n-2 缺少条件,通常在插值区间的端点给出,称为边界条件。 4.常用的三种边界条件 1°已知两端的一阶导数值,即: 2°已知两端的二阶导数值,即: 3°当f(x)是以xn-x0为周期的周期函数时,则要求S(x)也是周期函数,即 周期样条 三、 求解方法之一:三转角方程 设在[a,b]上给出插值条件: 1.条件 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn 求三次样条插值函数 S(x) 思路: (1)首先要补条件:每个区间上构造三次多项式需要四个条件,但现在最多有三个,故要补充条件,形成四个; (2)补什么条件:或函数值,或一阶导数值,或二阶导数值。这里选一阶导数较合适; (3)如何补?若随意给,则只能保证构造出的插值函数的函数值和一阶导数值连续,但不一定能保证二阶导数值连续,故只能选那组使二阶导数连续的一阶导数值。 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn x1处: 得到与m0,m1,m2有关的等式 x2处: 得到与m1,m2,m3有关的等式 共n-1个等式 设在[a,b]上给出插值条件: 1.条件 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn 求三次样条插值函数 S(x) 设法求出 求解过程具体如下: 2.求解 mj 的思路 由内部节点上的二阶导数连续求出 考虑S(x)在[xj , xj+1]上的表达式 hj=xj+1-xj 对S(x)求二阶导数得: 于是 同理可得S(x)在区间[xj-1 , xj]上的二阶导数: 于是 由条件 可得 进一步简化为 写成矩阵形式为 四、求解方法之二:三弯矩方程 设在[a,b]上给出插值条件: 1.条件 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn 求三次样条插值函数 S(x) 思路: (1)首先要补条件:每个区间上构造三次多项式需要四个条件,但现在最多有三个,故要补充条件,形成四个; (2)补什么条件:或函数值,或一阶导数值,或二阶导数值。这里选二阶导数较合适; (3)如何补?若随意给,则只能保证所构造出的函数的函数值和相邻两段在公共点的二阶导数的极限值连续,但不能保证一阶导数连续,故只能选那组使一阶导数连续的二阶导数值。 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn x1处: 得到与M0,M1,M2有关的等式 x2处: 得到与M1,M2,M3有关的等式 共n-1个等式 求解方法之二:三弯矩方程 设在[a,b]上给出插值条件: 1.条件 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn 求三次样条插值函数 S(x) 2.求解S(x)的思路及求解 1)首先确定S(x)与二阶导数值的关系 2)求出中间节点上的一阶导数值 1)首先确定S(x)与二阶导数值的关系 由于S(x)在区间[xj ,xj+1]上是三次多项式, 故S〃(x)在[xj ,xj+1]上是线性函数,可表示为 对S〃(x)积分两次并利用S(xj)=yj 及S(xj+1)=yj+1
文档评论(0)