水准网测量平差matlab_基于MATLAB的测量平差计算

2013年第 7期 中州煤炭 总第 211期 基于 MATLAB的测量平差计算 高 思培 ，陈冠 宇 ，范新 华 ，王 耀鑫 ，王 强 昆 (1，河南理 工大学 测绘 与国土信息工程学院，河南 焦作 454000； 2．桂林理 工大学 测绘地理信 息学院，广西 桂林 541006；3．机械工业勘察设计研 究院，陕西 西安 710043) 摘要 ：运用 MATLAB软件对测量数据平差处理进 行 了尝试 ，给 出间接平差 中水准 网平差 的函数程 序和 实例 解答 ，并介 绍 了使 用 MATLAB优化工具箱处理 附线性不等式约束的最小二乘平差 问题的具体 步骤和 函数 程 序。MATLAB与其他算法语 言相 比，具有编程简单 、运 算速 度快的特点 ，特别是在 矩阵运 算方面，在最小 二乘 平差 中可 以发挥很好 的作用 ，能大大提高工程测 量中的数 据处 理效率。 关键 词 ：平 差；MATLAB优化箱 ；附线性不等式约束 中图分类号 ：TD172．3 文献标志码 ：B 文章 编号 ：1003—0506(2013)07—0053—03 1 MATLAB软件简介 MATLAB是 MathWorks公 司开 发 的科学 与工 程 计算软件 ，它以矩阵运算 为基础 ，把计算 、绘 图及动 态系统仿真等功能有机融合在一起 。MATLAB将高 性能的数值和符号计算功能、强大的绘图功能、程序 语言设计功能以及为数众多的应用工具箱集成在一 起 ，其核 心是 一个 基 于 矩 阵 运算 的快 速 解 释 处 理 程 序 ，它提供了一个开放式的集成环境 ，以交互式操作 接受用户输入的命令 ，然后输出结果 ，以满足用户需 求 。 2 测量平差原理和模型 2．1 间接 平差 2．1．1 原 理及模 型 在 一 个 平 差 问题 中 ，当所 选独 立 参 数 的个 数 等 于 必要 观测 数 时 ，可 以将 每 个 观 测 表 达成 这 个 参 数 的函数 ，组成观测方程 。设最或然值为 ￡，则给出间 接平 差 的线 性化模 型 为 ： L = B X d n ×1 × l f× 1 n × 1 平差时对 取近似值 ，令 =X。+ ，代人上 式 ，并 令 ，=L一(BX。+d)=L—L。，所 以可 得 以 下 方 程 。 误 差 方程 ： =曰 一Z 根据 最dxZ-乘 法原 理得 ： 收稿 日期：2013—05—13 作者 简介 ：高思培 (1987一 )，男 ，河 南永城人 ，河 南理工 大学在 读硕 士研究生 ，研究方 向为大 型建筑物变形监测及数据处理研究。 = ( 胎 ) B P1，P为 观测 值 的权 阵。 平差值 向量的协 因数阵 ： QLL=B(曰 PZ) 曰 平差值 向量的协 方差 阵 ： DLL=~oQLL= P 单位权 中误差 ： |v1Pv √— 根据平差模型，可 以明确 MATLAB矩 阵计算 中 的 目的函数 和约束 条 件 ，从 而根 据 条 件 对 测 量数 据 进行 相应 处理 。 2．1．2 间接 平差 实例解 算 为 了探讨 MATLAB在 平 差 计 算 中运 用 的 可 行 性 ，这 里 选 取 某 一 工 程 实 例 进 行 MATLAB平 差 解 算 ，并和传统平差方法进行对照(图 1)，比较两者 的 计算工作量和精度¨。]。 D 图 1 水 准路线 示意 C 在图 1所示 的水准 网中，已知水准点 的高程 Ha=237．483 m，为求 B、C、D三点的高程 ，进行了水 准测量 ，其结果 见表 1，试 按间接平差求定 B、