几种插值法简介.docVIP

举例来看：可以认为某水文要素T随时间t的变化是连续的，某一个测点的水文要素T可以看作时间的函数T=f(t)，这样在实际水文观测中，对测得的（n+1）个有序值进行插值计算来获取任意时间上的要素值。 ①平均值法：若求Ti和Ti+1之间任一点T，则直接取T为Ti和Ti+1的平均值。 插值公式为：T= EQ \F(Ti+Ti+1,2) ②拉格朗日（Lagrange）插值法：若求Ti和Ti+1之间任一点T，则可用Ti-1、T1、Ti+1三个点来求得，也可用Ti、Ti+1、Ti+2这三个点来求得。 前三点内插公式为：T= EQ EQ \F((t-ti)(t-ti+1),(ti-1-ti)(ti-1-ti+1)) Ti-1+ EQ \F((t-ti-1)(t-ti+1),(t-ti-1)(t-ti+1)) Ti+ EQ \F((t-ti)(t-ti-1),(ti+1-ti)(ti+1-ti-1)) Ti+1 后三点内插公式为：T= EQ EQ \F((t-ti+1)(t-ti+2),(ti-ti+1)(ti-ti+2)) Ti+ EQ \F((t-ti)(t-ti+2),(ti-ti)(ti-ti+2)) Ti+1+ EQ \F((t-ti)(t-ti+1),(ti+2-ti)(ti+2-ti+1)) Ti+2 为提高插值结果可靠性，可将前后3点内插值再进一步平均。 ③阿基玛（Akima）插值法：对函数T=f(t)的n+1个有序型值中任意两点Ti和Ti+1满足： f(ti)=Ti EQ \F(df,dt) |t-ti=ki f’(ti+1)=T’i EQ \F(df,dt) |t-ti+1=ki+1 式中ki，ki+1为曲线f(t)在这两点的斜率，而每点的斜率和周围4个点有关，插值公式为：T=P0+P1(t-ti)+P2(t-ti)2+P3(t-ti)3，来对Ti和Ti+1之间的一点T进行内差。 ④牛顿（Newton）插值法：若求Ti和Ti+1之间任一点T，插值公式为： T=f(x0)+(x-x0)f(x0,x1)+ (x-x0)(x-x1)f(x0,x1,x2)+…+(x-x0)(x-x1)…(x-xn-2)f(x0,x1,…,xn-1) 式中，f(x0,x1)，f(x0,x1,x2)，…f(x0,x1,…,xn-1)是函数f(x)的1到第n-1阶差商。 f(x0,x1)= EQ EQ \F(f(x0)-f(x1), x0-x1) f(x0,x1,x2)= EQ \F(f(x0,x1)- f(x1,x2),x0-x2) f(x0,x1,…,xn-1)= EQ \F(f(x0,x1,…,xn-2)- f(x0,x1,…,xn-1),x0-xn-1) 可以看出，每一阶的差商都可以从它的前一阶差商推出。通常按照水文数据的特点，选定牛顿插值的阶数3～4，然后计算各阶差商，按照插值公式计算插值点的值。 ⑤线性插值法：若求Ti和Ti+1之间任一点T，插值公式为： T=Ti+ EQ \F(Ti+1-Ti, ti+1-ti) (t-ti) 对于给定的插值点，首先查找插值区间，即查找与插值点相邻的两个实测点，然后根据上式进行插值计算，得出相应的插值结果[6,7]。 关于插值效果的分析： 判断各种插值方法插值效果的优劣，通过两种途径来判断：一是图形法，即利用插值曲线与实测曲线的拟合效果判断插值效果；二是最小二乘法，即将实际观测值Tr与插值计算值Tc之间的离差(Tr-Tc)或计算误差的平方和(∑(Tr-Tc)2)作为判断的根据，当∑(Tr-Tc)2最小时，说明反应曲线拟合效果最好，插值点与实测点的离散程度最低，∑(Tr-Tc)2越大，说明插值点与实测点的离散程度越大。