(整理)EXCEL内置函数在测量平差中的应用..docxVIP

精品文档 精品文档 精品文档 精品文档 EXCEL内置函数在测量平差中的应用王海英 陈世梅 （南宁市勘测院 广西南宁 530021 ） 【摘 要】 测量平差是测量数据处理的重要环节，然而平差计算过程中，线性方程解算的问题却比较麻 烦。 本文利用 Excel 的内置函数简化了线性方程解算的问题， 并结合具体的实例， 讲述利用 Excel 内置函数 进行测量平差的方法。 【关键词】 矩阵 Excel 测量平差 0 引言 在现代测量中，对通过各种测量方法所采集得到的原始数据，往往需要根据误差理论的方法对存在的 各类误差进行平差处理从而取得最或然结果。而对于偶然误差的处理，利用最小二乘法准则进行平差计算 的过程， 通常都要对误差方程式或条件方程式进行整合处理求得法方程组， 然后解算法方程、 计算改正数、 精度评定等等， 一系列计算步骤的进行， 如果没有现成的专业软件， 而用手工的办法去完成平差计算工作， 将是一件很繁琐的事情。 其实，测量平差的过程简单的来说就是解算线性方程组的过程，而最令人烦恼的部分是线性方程的求 解。在处理一些小的工程项目平差计算的时候，如果手头没有像 MATLAB 等具有矩阵运算功能的应用软 件，利用 Microsoft Office Excel 的内置函数，同样的可以帮助计算者比较轻松的完成计算任务。 以下要阐述的就是利用 Excel 转置粘贴功能以及矩阵计算的函数 MMULT （矩阵乘） 、 MINVERSE （矩 阵求逆），实现测量平差之线性方程组解算的过程。 平差的数学模型 为了加快平差解算的作业效率，应该根据实际情况选择适当的数学模型。实际工作中，有两种数学模 型得到了较为广泛的应用，即间接平差模型和条件平差模型；它们的应用公式 [1] 分别如下： 间接平差 误差方程：V = B Sx+ J; J = B X0+ B0 L ; 法方程：BT P B Sx+ BT P J = 0 ; 改正数计算：Sx = Y BT P B) -1 BT P J ; PVV 计算：V-1 P V = JT P V + SxT BT P V ; 权函数式：S①=F x 平差值的权倒数计算： 1 / P①二FT (BT P B ) -1 F ; 公式中各种符号的含义：观测向量 L;相应的权阵 P;误差方程的系数矩阵 B;未知数的改正数 Sx 常数阵J;常数向量B0;未知数的近似值 X0;权函数的系数阵Fo 条件平差 条件方程式： A V + W= 0 ； W = A L + A0 ； 法方程： A P-1ATK + W = 0 ； 改正数计算： V = P-1 AT K ； PVV 计算： V-1 P V = — WT K ； 权函数式：S①=F ; S i 平差值的权倒数计算： 1 / P ①二 FTPI F — (A P-1 F) T (A P-1AT ) -1 A P-1 F ; 公式中各种符号的含义：观测向量 L;相应的权阵 P;条件方程的系数矩阵 A;条件改正数 V;条件 闭合差 W;法方程联系数矩阵 K;权函数的系数阵Fo 计算方法 不难看出，以上平差计算的数学模型，全都是以矩阵的形式给定的。如果掌握了矩阵计算的方法，所 有的计算将迎刃而解。 由于测量计算作业过程中采用不同的数学模型，以及图形条件的多样化，对于误差方程或者条件方程 的产生，在本文所述的解算方案当中，还不能自动完成；当然，可以利用 Excel 中的 VBA 功能编程处理， 但将可能给计算者增加了难度。在此仅就手工生成误差方程组或者条件方程组之后，如何使用 Excel 内置 函数进行矩阵计算的方案进行阐述。 数据的输入 将系数矩阵的数据填入电子表格的相应区域中。 矩阵的转置 （ 1 ）选中电子表格中系数矩阵的数据区域，点击右键，选击 ‘ 复制 ’ 菜单选项。 （ 2）在电子表格中点击系数矩阵数据区域以外的任意单元格，点击右键，选击 “选择性粘帖 ”菜单项。 （ 3）点击 “选择性粘帖 ”菜单项，出现 “ 选择性粘帖 ” 对话框；该对话框有三部分选项，第一部分 “粘帖 ”， 选中 “全部 ” 选项，第二部分 “ 运算 ”，选中 “无 ”选项，第三部分选中 “ 转置 ”选项，单击确定，即可完成系数矩 阵数据的转置操作。 矩阵乘运算 （ 1 ）选中电子表格中的空白区域，该区域的行数、列数等于矩阵相乘所得矩阵的行数、列数，点击工 具栏中的fx （函数）工具按钮。 （ 2）弹出插入函数的对话框，在该对话框的左边的函数分类中选择 “ 数学与三角函数 ” ，在右边的函数 名中选中 “ MMUL” T ，单击确定。 （ 3）弹出函数参数输入的对话框，该对话框提示输入两组参数，第一个参数 Array1 栏中输入第 4步 骤中转置矩阵区域的行列号， 在 Array2