基于MATLAB的水准网平差程序设计.pdfVIP

70 号稿 来稿日期 2013-05-06 基于 MATLAB 的水准网平差程序设计 王鹏磊 254330952@ （西安科技大学,陕西省西安市雁塔中路 58 号，710054 ） 【摘要】MATLAB 是近年逐步发展起来的大型数学计算软件，在测绘数据处理中用途越来 越广泛，本文对水准网平差的算法进行分析，分别基于条件平差和间接平差两种方法，利用 MATLAB 实现了水准网的平差。 【关键词】MATLAB ；水准网；条件平差；间接平差 测绘科学是一门以大规模数据甚至是海量数据处理、分析与应用为基础的学科，其各 项工作都涉及大量的数值计算与分析，运用 MATLAB 来解算测绘数据处理问题具有其他程 序设计语言无可比拟的优越性。因此，本文鉴于 MATLAB 对测绘数据处理的优势，主要利 用其 GUI 功能，分别基于条件平差和间接平差两种模型实现了水准网的平差。 一、平差原理 1、条件平差 条件平差的基础方程为：AV +W 0 ；V P -1A T K 按最小二乘原理，组成法方程：AP-1A T K +W 0 解算法方程求出联系数 K ，将其代入改正数方程即可求出改正数，进而解算平差值并进 行精度评定。 2 、间接平差 间接平差的基础方程为：V Bx - l ；BT PV 0 解算此基础方程，消去 V ，得法方程BT PBx - BT Pl 0 ∧ 将上式求解 T -1 T ，进而可求改正数并进行相关精度评定。 x （B PB ）B Pl 二、水准网平差程序设计思路及关键程序 1、程序设计流程 程序设计具体流程示意图，如图 1 所示。 输入点号、高程和路线长度 选择平差方法 输入观测数、已知点数、未知点数 读取观测数据文件 建立传递系数阵 建立系数阵、闭合差向量 建立系数阵、常数项 平差计算 导出平差报表 结束 图 1 程序设计流程图 Figure 1 Flow chart of program design 2 、关键程序 1)因为条件平差中条件方程的多样性，要得出系数矩阵A 和闭合差向量 W 具有一定的 难度，引入了一个传递系数阵 T，来简化条件方程的建立。初始传递系数阵 T 为全 0 矩阵， 然后根据各条水准路线起点、终点的高程是否已知、近似高程是否已知四种情况，对矩阵 T 中相应的元素进行赋值或其他操作。在得到系数矩阵T 后，进而计算系数阵和常数项。 以下为条件平差中系数阵a 和常数项 w 的计算 a=zeros((gd-dd),gd); %gd、dd分别为总观测数和未知点数 w=zeros((gd-dd),1); r=0; while r(gd-dd) for k=1:gd i=qd(k); j=zd(k); if tip(k)~=1 r=r+1; a(r,1:gd)=t(i,1:gd)-t(j,1:gd); a(r,k)=1; w(r)=h0(i)+h1(k)-h0(j);