第一章 数据处理(插值法).ppt

* * 1.1 插值法 1.1.4 差分与等距节点插值公式 实际应用时经常遇到等距节点的情形，这时插值公式可进一步简化，计算也简单很多。为了得到等距节点的插值公式，先介绍差分的概念。 （1）差分及其性质 一阶向前差分 一阶向后差分 一阶中心差分 * * 利用一阶差分可定义二阶差分为 （1）差分及其性质 一般地可定义阶差分为 1.1 插值法 1.1.4 差分与等距节点插值公式 * * （2）差分与差商的关系 1.1 插值法 1.1.4 差分与等距节点插值公式 因此， 同理，差商与相向后差分关系为 * * 1.1 插值法 1.1.4 差分与等距节点插值公式 （3）等距节点插值公式 引入 牛顿前插公式 * * 1.1 插值法 1.1.4 差分与等距节点插值公式 （3）等距节点插值公式 同理，引入（p为负值） ，可得 当引入 ，可得牛顿前插公式 牛顿后插公式 * * （3）等距节点插值公式 例 利用下列数据，求液氮在 下的蒸气压。 解：由所给数据做差分表 * * （3）等距节点插值公式 前插 后插 前插外推 后插外推可得到相同的结果。 * * 1.1 插值法 1.1.5 分段插值法 （1）分段线性插值 求一个分段函数P(x), 使其满足: （1）P(xi)=yi (i=0,1, ..., n) （2） 在每个子区间[xi，xi+1] 上是线性函数. 称满足上述条件的函数P(x) 为分段线性插值函数. * * （2）分段抛物线插值 (2) 在每个子区间[xi-1, xi+1] 上，L(x)是次数不超过2的 多项式. 称满足上述条件的函数L(x)为分段抛物线插值函数. （1）L(xi)=yi (i=0,1, ..., n); 对 求一个分段函数L(x), 使其满足: 即将区间[a, b]分为小区间[xi-1, xi+1] (i=1,2, …,n) * * 1.1 插值法 1.1.6 三次样条插值函数 （1）定义 给定区间[a,b]的一个划分 a=x0x1…xn=b， yi=f (xi) (i=0,1,…,n),如果函数S(x)满足： （1）S(xi )=yi (i=0,1,…,n); （2） 在每个小区间[xi, xi+1] (i=0,1,...,n-1)上是次数不超过3的多项式; (3) 在每个内节点xi (i=1,2,...,n-1)上具有二阶连续导数, 则称 S(x) 为关于上述划分的一个三次样条插值函数。 * * ① ③ 在任一区间 上是次数不高于3次的多项式。 ② 连续性 （1）定义 S(x)满足下列条件： 1.1 插值法 1.1.6 三次样条插值函数 * * 1.1 插值法 1.1.6 三次样条插值函数 （2）求解函数表达式 A 待定系数法 由三次样条函数的定义可得： 另外2个方程由边界条件给出。 * * （2）求解函数表达式 A 待定系数法 例 解：由已知条件知，由两个区间就有两个待定方程。 首先利用插值条件 * * （2）求解函数表达式 A 待定系数法 即可求出 然后连续条件 边界条件 * * （2）求解函数表达式 B 三弯矩法 思 路： 是一个次数不超过三次的多项式，则 为一个线性函数，对 进行积分，即可求出 。 积分可得： 合并可得： * * （2）求解函数表达式 B 三弯矩法 * * （2）求解函数表达式 B 三弯矩法 三对角矩阵 * * （2）求解函数表达式 B 三弯矩法 例： 解： * * （2）求解函数表达式 B 三弯矩法 同理可得 例：给定数据表如下： （2）求解函数表达式 B 三弯矩法 （2）求解函数表达式 B 三弯矩法 （2）求解函数表达式 B 三弯矩法 * * 第一章 数据处理 1.2 数值微分 1.2.1 用差商近似微商 在微积分中函数的导数是通过极限定义的。当函数以表格给出时，就不可能用定义求它的导数，但工程上又需要求列表函数在节点和非节点处的导数值，这就是数值微分所要解决的问题。 在微积分中，微商是极限的概念 * * 1.2 数值微分 1.2.1 用差商近似微商 第一章 数据处理 取其达到极限以前的形式，即可得到微商的差商近似式。 一阶向前差商 一阶向后差商 一阶中心差商 * * 1.2 数值微分 1.2.2 用插值函数计算微商 第一章 数据处理 * * 1.2 数值微分 1.2.2 用插值函数计算微商 第一章 数据处理 t