插值法综述《计算方法》学习报告.doc

插值法综述 一、插值法及其国内外研究进展 1.插值法简介 插值法是一种古老的数学方法，它来自生产实践，早在一千多年前，我国科学家在研究历法上就应用了线性插值与二次插值，但它的基本理论却是在微积分产生之后才逐渐完善的，其应用也日益增多，特别是在计算机广泛使用之后，由于航空、机械加工、自动控制等实际问题的需要，使插值法在实践和理论上都显得更为重要，并得到了空前的发展。 2.国内外研究进展 插值法在预测地基沉降的应用 插值法在不排水不可压缩条件下两相介质的两重网格算法的应用 拉格朗日插值法在地震动的模拟研究中的应用 插值法在结构抗震可靠性分析中的应用 插值法在应力集中应变分布规律实验分析中的应用 3.代表性文献 不等时距GM(1%2c1)模型预测地基沉降研究 秦亚琼 武汉理工大学学报(交通科学与工程版) 2008.2 不排水不可压缩条件下两相介质的两重网格算法 牛志伟 岩土力学 2008.3 基于拉格朗日插值法的地震动的模拟 白 可 山西建筑 2010.10 响应表面法用于结构抗震可靠性分析 张文元 世界地震工程 1997 小议应力集中应变分布规律的实验方法 查珑珑 淮海工学院学报 （自然科学版）2004.6 二、插值法的原理 【原理】 设有n+1个互不相同的节点（，） （i=0，1，2，...n）则存在唯一的多项式： 使得 证明：构造方程组 令： 方程组的矩阵形式如下： 由于所以方程组（4）有唯一解。 从而唯一存在。 三、常用插值法 3.1 Lagrange插值法 ，多项式 称为关于的次Lagrange插值多项式。 3.1.2 Lagrange插值多项式的构造 已知n+1个节点其中互不相同，不妨设要求形如： 的插值多项式。 若n次多项式在n+1个节点上满足条件： 就称这n+1个n次多项式为节点的n次插值基函数。 3.1.3 Lagrange插值法的程序设计 f[x_]:=Exp[x] A=Table[{x,f[x]},{x,0,0.8,0.2}]//N g1=ListPlot[Table[A],Prolog-AbsolutePointSize[18]]; Interpolation[A,InterpolationOrder-3] g2=Plot[%[x],{x,0,0.8}] Show[g1,g2] N[%%%[0.12],20] N[%%%%[0.72],20] N[f[0.12],20] N[f[0.72],20] 3.1.4 Lagrange插值法典型例题及其解法 已知，构造二次拉格朗日插值多项式。 （1）计算； （2）估计误差并与实际误差相比较。 解 （1）以插值点(27,3), (64,4), (125,5)代入插值公式，得 = (2) 由误差公式有 记在[27，125]上是单调递减函数。 实际误差：。 3.1.5 Lagrange插值法误差估计 3.2 Newton插值法 称为Newton插值多项式。 3.1.2 Newton插值多项式的构造 由插值条件 当时,. 当时,,推得 . 当时， 则引入记号： 依次递推可得的一般表达式： 3.1.3 Newton插值法的程序设计 {x[0],x[1],x[2],x[3],x[4]}={10,11,12,13,14}; y[k_]:=Log[x[k]] Table[y[k],{k,0,4}]//N; MatrixForm[%] f[i_,j_]:=(y[j]-y[i])/(x[j]-x[i]) Table[f[i,i+1],{i,0,3}]//N; MatrixForm[%] f[i_,j_,k_]:=(f[j,k]-f[i,j])/(x[k]-x[i]) Table[f[i,i+1,i+2],{i,0,2}]//N; MatrixForm[%] f[i_,j_,k_,l_]:=(f[j,k,l]-f[i,j,k])/(x[l]-x[i]) Table[f[i,i+1,i+2,i+3],{i,0,1}]//N; MatrixForm[%] f[i_,j_,k_,l_,m_]:=(f[j,k,l,m]-f[i,j,k,l])/(x[m]-x[i]) Table[f[i,i+1,i+2,i+3,i+4],{i,0,0}]//N; MatrixForm[%] A={{y[0],y[1],y[2],y[3],y[4]},{0,f[0,1],f[1,2],f[2,3],f[3,4]}, {0,0,f[0,1,2],f[1,2,3],f[2,3,4]},{