[2018年最新整理]DEM表面建模.pptVIP

第六章 DEM表面建模及精度分析模型 6.1 数学基础 6.2 数字高程模型概念 6.3 DEM表面建模 6.4 DEM内插模型 6.5 基于正方形格网的DEM精度模型 6.6 基于正方形格网的高阶插值传递误差模型 6.7 基于TIN的DEM传递误差模型 6.8 地形描述误差模型 6.1 数学基础：一元函数插值与逼近（1） 一、引言： 函数y = f(x)以离散点的形式给出时，需构造一个函数y = P(x)近似代替y = f(x)，同时要求P（x）在点x0,x1,…,xn上满足 P（xi）=yi， i=0,1,…,n (7.1) 寻找满足上述条件的P（x）的问题称为插值问题。条件(7.1)称为插值条件；f(x)称为被插值函数；P（x）称为插值函数； x0,x1,…,xn称为插值结点。 f(x)- P（x）称为插值余项，记为R（x）。 6.1 一元函数插值与逼近（2） 二、n次代数插值多项式： 拉格朗日插值多项式： 拉格朗日基多项式： 当n=1时，用L1(x)近似代替f(x)称为线性插值，其线性公式为 即： 当n=2时，用L2(x)近似代替f(x)称为抛物插值，抛物插值公式为 即： 6.1 一元函数插值与逼近（3） 三、埃尔米特插值多项式： 埃尔米特插值是已知函数f(x)在n+1个互不相同的点上的相应的函数值及导数值，要求寻找次数不超过2n+1的代数多项式 利用插值条件，构造如下Hermite插值多项式： 特别当n=1时，得到三次埃尔米特插值多项式为： 6.1 一元函数插值与逼近（4） 假设有一函数： 龙格发现如果使用小于等于n阶多项式在-1与1之间按照 这样的等距点进行插值，那么在接近端点-1与1的地方插值结果就会出现震荡。在多项式阶数增高时插值误差甚至会趋向无穷大。 6.1 一元函数插值与逼近（5） 四、分段低次插值： 将插值区间分成若干个小区间，然后在每个小区间上进行了低次插值，这时在每个小区间上的插值取为次数较低的多项式，这样的插值方法称为分段低次插值。 1、分段线性插值： 将区间[a,b]划分为 a=x0x1…xn=b 在每一小区间[xj-1,xj]上作线性插值，得 分段线性插值的特点： ①用折线代替曲线，又叫折线插值。②在整个插值区间[a,b]上连续，但在[a,b]上一阶导数不连续。③其精度取决于取样点的密度。④取样点的分布可以是不规则的，可根据函数的变化调整。 6.1 一元函数插值与逼近（6） 2、分段三次多项式插值： 将区间[a,b]划分为 a=x0x1…xn=b 取四个点xi-1,xi,xi+1,xi+2,从而在区间[xj-1,xj+2]上得到一条三次多项式曲线，取其中间区间[xj,xj+1]的函数值。对于两点端点区间上的插值，可用线性插值代之，或将相邻间隔上的内插曲线延伸过去。 分段三次多项式插值的特点： ①在大区间[a,b]上连续，但不能保证插值函数的一阶导数在插值点的连续性，即曲线的光滑性。 ②较复杂，效率低。 ③若插值结点等距分布，则处理效率会高很多。 6.1 一元函数插值与逼近（7） 3、加权平均的抛物线拟合法： 基本思想：按数据点P0，P1，…，Pn的顺序，每过相邻三点作一条二次曲线，在每相邻两点之间的前后两曲线的重叠部分，用平均加权曲线作为最终的插值曲线。 设过P0，P1，P2的二次曲线方程为 再设过P1，P2，P3的二次曲线方程为 在两点之间的最终曲线上采用左右加权的曲线 6.1 一元函数插值与逼近（8） 3、加权平均的抛物线拟合法： 式（7.2）须满足条件： 权函数须满足条件： 权函数取为 在P0和P1之间、Pn-1和Pn之间没有中数可取，以原来的一次抛物线代之。 6.1 一元函数插值与逼近