12-MATLAB在测量数据处理中的应用.pptVIP

一个条件平差的例子 function [kk,vv,dert0]=tjpc(aa,pp,ww) pni=inv(pp); naa=aa*pni*aa; nni=inv(naa); kk=nni*ww; vv=pni*aa*kk; vtpv=kk*ww; [rr,nn]=size(aa); dert0=sqrt(vtpv/rr); * 第七章 MATLAB在测量数据处理中的应用 7.1测量平差中的应用：条件平差、间接平差、附有未知数的条件平差、附有未知数条件的间接平差。(要求熟悉模型) 7.2高程网平差设计。(理解模型、用会) 7.3导线网(或GPS三维平差) 。 (理解模型、用会程序计算) 7.4其它应用：大地测量计算、工程测量、图像处理等。(根据自己情况自由发展) 7.2高程网平差设计。 一、主要任务： 水准网概算及平差(掌握) 三角高程网平差：高差计算、(选择未知点高程和折光差C（K）为未知数） 动态水准网平差. 二、建立模型 1、一段高差为一个结构元素(一条记录)，建立观测值结构矩阵 2、一个点的高程为一个结构元素(一条记录)，建立高程结构矩阵 3、定权： 4、起算点(已知点)统计，用t0表示 5、未知点高程近似值推算(核心部分) 推算原理： ①未知点高程值均为0(或-9999) ②当一段高差的一点已经有高程大于0，另一个高程为0(或-9999)，则高程要推算出，其它情况不需要推算。 ③每计算出一个高程未知点个数加1 。把计算出一个高程放到高程结构矩阵中。 ④直到所有未知点高程推算结束 。 6、平差计算：列误差方程式或直接列出法方程(课堂上不讲，有意者课外辅导) 高差改正数=-起点高程改正数+终点高程改正数-常数项(观测值-观测值计算值) 7、精度评定 高程平差值精度； 高差精度。(见下面) 8、平差报告 输入值 输出值和精度 精度分析 高差精度。 高差平差值 高差两端点的协因数 水准网平差综合例子 掌握： 1、水准网平差的结构和步骤； 2、数据准备和使用程序解决实际问题； 3、要读懂程序很难，希望同学们不要在本课程结束前花费太多的时间。如果有兴趣的话毕业设计再仔细研究。 %结构体survey存放高差观测值及相关信息 qdm＝起点名 zdm＝终点名 gcha＝高差观测值 line=水准路线长度 qdgc=起点高程 zdgc=终点高程 qdh=起点号 zdh=终点号 quan=权 gcgj=高差估计值 surveyw=高差平差值中误差 %结构体height存放已知点高程 dm=点名 gc=高程 surveyw=高程中误差 %程序结果分别放进文件height.clk,surveying.clk中 7.3边角网导线网平差算例 本例选自于《误差理论与测量平差基础习题集》70页7.10.65 总信息：观测值个数n=8+2=10,t=2*2=4 设角观测值中误差为单位权中误差，所有角度的权为1，边长的权为 此例中度分秒表示的角度要转化为弧度计算 已有函数deg2rad(dd.mmss) 在列立角度和边长误差方程式要用到方位角和边长的近似值（计算值），所以有自定义函数 [rf,jl]=pmzbfs(x1,y1,x2,y2) 此例暗示坐标改正数的单位应选厘米，使误差方程单位一致。边长改正数的单位应选厘米 误差方程式的列立：引入符号，设J到k点的方向系数为： 其中 边长误差方程 给定未知点近似坐标情况下的导线网(边角网)平差程序设计的步骤： 1、定义或从文本文件中读取已知数据信息； 2、定义或从文本文件中读取角度观测和边长值信息； 3、定权：角度的权和边长的权 4、列角度误差方程 5、列边长误差方程 6、平差计算 7、精度评定 8、平差报告 一个点的坐标协因数是 特征方程是 特征值是 时一三象限 边角网(导线网)平差例子 掌握： 1、边角网平差程序的结构和步骤； 2、数据准备和使用程序解决实际问题； 3、本程序简单，读懂程序不是很难。 4、本程序要改进的地方：用结构矩阵存放观测值、已知数据 本程序要改进的地方： 用结构矩阵存放观测值、已知数据； 未知点坐标的推算最好让计算机完成； 增加输出每条边的相对误差椭园参数、边长中误差、方位角中误差； 每个观测值的多余观测分量和中误差。 实验十测量平差程序设计 从《误差理论与测量平差基础习题集》中选一个例题：如水准网平差或边角网，用例题中的数据(权、条件方程系数和常数项 、误差方程和常数项)，选择条件平差、间接平差、附有未知数的条件平差、附有未知数条件的间接平差的一种进行平差计算，验证书上例题是否正确。建议同学们选择有挑战性的例题，用较为复杂的数学模型。 *