第5章由离散点绘制光滑曲线的方法和程序实现.ppt

观测值改正数特点 6、求插值点坐标： P83,(5-47)式 观测值改正数特点 张力样条函数插值法程序 变量说明如下： N：已知离散点个数 X(N),Y(N):离散点坐标 XP(250),YP(250):先用于记录线性方程组的右端项， 后用于记录线性方程组的未知项。 S(250):累加弦长 A(250),B(250),C(250):三对角线矩阵的元素 SX1,SY1:首点导数 SXN,SYN:终点导数 SIGMA:张力系数 第五章测量误差基本知识 第五章 由离散点绘制光滑曲 线的方法和程序实现 第五章测量误差基本知识 线性迭代光滑法 分段三次多项式插值法 （五点光滑法） 二次多项式平均加权法 （正轴抛物线平均加权法） 张力样条函数插值法 ?迭代思路： ?迭代次数与插值点数的关系： §5-1 线性迭代法 A B C A’ B’ 4 3 2 1 1’ 2’ 4’ 3’ 5’ 6’ ?特点：曲线不通过原节点而向内收缩 ?特点： §5-2 分段三次多项式插值法 （五点光滑法） 每二个数据点之间建立一条三次多项式曲线 方程； 整条曲线具有连续的一阶导数； 各节点的一阶导数是以一点为中心，两边相 邻各二点来确定。 ?思路： 已知N个离散点，如果每个离散点上的导数是已知的，那么在任意两个相邻的离散点之间就有四个条件 采用以上四个条件，可以在两个离散点之间拟合一条三次多项式曲线。 ?由AKIMA的几何条件： D C B A 5 4 3 2 1 式中： 一、由于地图中的曲线是多值函数： 式中： 二.补点方法 ?开曲线： 2 3 4 5 n-4 n-3 n-2 n-1 n 1 1 5 n n-2 n-4 3 二.补点方法 ?闭曲线： 2 n-1 n-3 4 三.优点： 克服了线性迭代法曲线不通过数据点的缺点，所绘出的曲线严密地经过每一个数据点，整条曲线具有连续一阶导数的光滑性。 四.分段二次多项式插值法程序： 频率直方图 ?在2～3之间将产生两条 抛物线； ?通过前三个相邻点1、2、3 可以拟合一条抛物线，通过 后三个相邻点2、3、4又可 拟合一条抛物线； ?在重叠范围2～3点内， 用加权的办法获得平均 曲线作为最终的插值曲 线(虚线)。 x y §5-3 二次多项式平均加权法 1 2 3 4 §5-2 评定精度的标准 上式中，s是弦长定义为： 一.对多值函数的参数方程 (5-32) §5-2 评定精度的标准 对前三点（1，2，3）： §5-2 评定精度的标准 对后三点（2，3，4）： §5-2 评定精度的标准 仿单值函数的推导过程： 得到插值曲线： 二.相对误差 二.补点方法 ?开曲线： 1 2 3 4 n-3 n-2 n-1 n 二.相对误差 1 3 n n-1 n-2 2 二.补点方法 ?闭曲线： 三.容许误差 三.程序：二次多项式平均加权插值法 三.容许误差 四.缺点： 抛物线的顶点并不是三点之中间点（即第二点），中间点偏离正轴抛物线顶点越远，越不能客观地反映原有图形。 三.容许误差 五.弥补办法： 采用斜轴抛物线平均加权法； 可以做到曲线的明显转弯点位于节点。 算术平均值及改正数 §5-4 张力样条函数法 显著特征：具有一个张力系数σ 当σ→ 0时，张力样条函数就等同于三次样条 函数； 当σ→∞时，张力样条函数就退化成分段线性 函数，即节点与节点之间是折线 连接。 当选择合适的σ时，点与点之间的曲线尽量缩 短，好像在整条曲线的两端用一种作用力拉到合适的程度，既能消除可能出现的多余拐点，又保 持了曲线的光滑性。 观测值改正数特点 一. 定义： N个点（xi,yi）(i=1,2,…,n),σ,求f(x)使 从而，得到P78 (5-35)式：二阶非齐次常系数线性微分方程 观测值改正数特点 解为： 齐次方程 将此式代入（5－35），比较等式两边系数，求得a,b 非齐次方程的特解： 最后定出曲线f(x)的表达式 即P79 (5-36)式； 式中关键要求 观测值改正数特点 给出两个端点条件： 二. 节点关系与端点条件： 构成（5－44）线性方程组；解得 代入（5－36）得：张力样条插值函数 观测值改正数特点 规范化张力系数： 三. 张力系数的选择： 观测值改正数特点 参数方程： 四. 应用： 上式中，s是累加弦长定义为： 观测值改正数特点 1、判断曲线类型，计算首尾