[自然科学]第二讲：插值法.ppt

在区间 上 （3） 为分段Hermite插值 §2.7 三次样条插值（简述） 1、样条的数学模型-样条函数 定义：如果函数 ，且在 上是三次多项式。 则称 是节点 上的 三次样条函数。 如果函数 满足 则称 是三次样条插值函数。 2、样条函数的求解 分析： 1、区间 上的三次多项式，共需待定系数 个。 2、已知条件有 个 个 个 个 共计 个 （1）补充条件 1、 边界条件 2、 特别 称自然边界条件 3、 周期性条件（略） （2）三转角方程 假定 在区间 上的Hermite插值函数为： 则 Hermite插值 ？ 所以 P36 （2.8.7） 所以 进而 整理得 P36 同理 第一步：在区间 上做Hermite插值函数 第二步：求 第三步：求 第四步：求 在 点的左极限得 由 得 化简得 其中 所谓的三转角 三转角方程： 的矩阵形式 第 一 类 边 界 条 件 如果已知 特别 则在区间 上 即 整理得 同理，在可得区间 上 所以 整理得 由 得 对于自然条件 有 因为 是三次多项式，所以 是线性函数。 （3）三弯矩方程 在区间 上令 积分两次得 事实上 因此 所以 类似用同样的方法在区间 做 经过两次积分，一次微分，并求左极限得 由 并整理得 三弯矩 其中 第 二 类 边 界 条 件 * * * * 例 求过点 所以 的二次插值多项式。 解一： 令 由 得 解二： 解三： 所以 已知等距节点 §2.4 Newton等距插值 1、差分 定义 简记为 简记为 简记为 向前差分 向后差分 中心差分 前差算子 后差算子 高阶向前差分 高阶向后差分 如 2、高阶差分 又如 3、前差与后差的关系 一般有 再定义 前移算子 不变算子 后移算子 则有 因此 4、差商与差分的关系 m阶向前差商与m阶向前差分的关系 m阶向后差商与m阶向后差分的关系 又 所以 5、差分的计算 6、等距节点的Newton插值 已知等距节点 得 令 由Newton插值公式 其中 参照 即 前插公式 同理可得后插公式（P27) §2.5 Hermite插值 已知数表 求多项式 满足 则 称为Hermite插值多项式 因为数表中有 个已知数，可确定一个 次 多项式。 当 较大时用待定系数法求 是困难的 令 为 次多项式 且满足 其中 且满足 令 为 次多项式 所以 为Hermite插值多项式。 Kronecker（克罗内克）符号 柏林科学院院士，巴黎科学院通讯院士，伦敦皇家学会外籍会员。 主张分析学应奠基于算术，而算术的基础是整数。克罗内克名言：“上帝创造了整数，其余都是人做的工作” 令 则 其中 又 则 由（1）（2） 得 所以 其中 则 所以 令 则 又 由 得 所以 例 已知 求三次多项式 满足 解 所以 验证： 练习：（ P43（19））求四次多项式 满足 令 解： 在点0，1，2上做Lagrange插值函数 则 因此 所以 考察函数 §2.6 分段低次插值 1、Runge现象 取点 做Lagrange插值 function y=lagrange(x0,y0,x); n=length(x0); m=length(x); for i=1:m s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(x(i)-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end end Lagrange函数 x0=-5:1:5; y0=1./(1+x0.^2); x=