控制网的优化设计 毕业论文 gis

第三章 GPS网的精度和可靠性 控制网的优化设计 摘要 优化设计是最优化理论和方法在设计中的应用，力求以最低的成本、最高的效率达到最优的目标。本文通过一系列的分析，对 控制网的优化方法进行分析，说明可行性。 为了解决控制网优化设计问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 ，本 论文 政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载 分两大部分，GPS网的优化设计和GPS网的精度和可靠性，在 GPS网形设计中，首先根据工程的特点和GPS网设计规范的要求，大致确定网的规模，用图论和树的有关算法推导出GPS网形中点、边、异步环之间的关系，然后给出一种生成网形的算法，自动生成初步网形，并用模拟法在顾及精度和可靠性准则下对初步网形进行优化设计，确定最终网形，并按最小路径方法生成观测 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 。 关键词： GPS控制网，优化设计，精度，可靠性 OPTIMIZING DESIGNING OF CONTROL NETWORK ABSTRACT The optimization design is a application of the most optimizative theory and method in the design. It is hopeful that the best target obtained through the expenditure of the fewest cost and the best efficient. The essay is a feasibility report of the optimization design of GPS control network’s methods by a series of analysis. This paper consists of two parts: Optimizing designing of GPS control network and the Precision and Reliability of GPS network. When designing a GPS control network ,its scale should be predicted as the project requested and the GPS surveying standard disciplined. According to the relationship among GPS points , edges and nonsynchronous loops, we can use an algorithm of Graphic Theory to produce a network when given the number of points and the maximum edges of each nonsynchronous loop, after being modified by using simulate optimizing method we can draw the ultimate network, then the observation plan can be gained by using the best way algorithm. KEYWORDS：gps control network, optimizing designing, precision, reliability 目录 TOC \o "1-3" \h \z \u 摘要 2 ABSTRACT 3 1绪论 5 1.1 控制网优化设计的发展历史 5 1.2 GPS控制网的发展历史 7 1.3 GPS网形的特点 8 1.4 GPS系统的应用前景 9 2 GPS网的精度和可靠性 11 2.1 GPS网平差的数学模型 11 2.2 误差的传递与转换 12 2.3 相对点位精度的合理评定 14 2.4 GPS控制网的可靠性 18 3 GPS控制网优化设计 22 3.1 GPS网优化设计的数学模型 22 3.2 一阶段优化设计 24 3.3 三阶段优化设计 27 结论 30 参考文献 31 致谢 32 前言 全球定位系统（Global Positioning System-GPS）是美国从20世纪70年代开始研制，历时20年，耗资200亿美元，于1994年全面建成，具有在海、陆、空、进行全方位实时三位导航与定位能力的新一代卫星导航与定位系统。最近十年我国测绘等部门的使用表明，GPS以全天候，高精度，自动化，高效益等显著特点，赢得广大测绘工作者的信赖，并成功的应用于大地测量、工程测量、航空摄影测量、运载工具导航和管制、地壳运动监测、工程变形监测，资源勘探、地壳运动学等多种学科，从而给测绘领域带来了一场深刻的技术变革。 随着全球定位系统的不断改进，硬、软件的不断完善，应用领域正在不断的开拓，目前已遍及国名经济各不呢不，并开始逐步深入人们的日程生活。 1绪论 1.1 控制网优化设计的发展历史 最优化问题是一个具有很久的研究历史和广阔的应用前景的领域，早在公元前500年毕达哥拉斯就发现了黄金长方形，即长方形的长与宽的最佳比率为1.61 8， 称为黄金分割比，在建筑或绘画中应用这个比例将使建筑和艺术最优美、协调。在微积分出现以前，己有许多人开始应用代数的或几何的方法来解决最优化问题;在微积分出现以后，利用求导法、变分法和拉格朗日乘数法，有效地解决了可导函数的极值问题，这类方法也称为古典最优化问题。近30年来，由于科学技术发展的需要，实践中许多最优化问题已无法用古典方法来解决，因此，许多新的最优化技术应运而生，为解决各种优化设计问题提供了有效的方法。目前最优化设计己普遍应用于国民经济的各个领域，如生产管理、运输调度、服务系统、信息系统等等。 在测绘学科方面，高斯所创立的最小二乘法至今仍是应用最广的数据处理方法。18年Helmert发表了关于“测量的合理性研究” 的论文，1882年Shreibe : 曾提出过著名的“线网最适当权分配”的方法，对基线扩大网的测角观测方案进行了优化设计，这可以看作是对测量控制网的优化设计所做的开创性工作。但是在相当长的一段时间内，这一问题未得到深入而系统的研究。直到六十年代，由于整个科学技术的发展和电子计算机的应用，测量控制网的优化设计问题引起了测绘学者的重视，尤其是近2 0年来，控制网的优化设一直是人们感兴趣的课题之一。值得提及的是Baarada对控制网的质量 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 问题进行了入的研究，提出了评价控制网质量的三项标准，即精度、可靠性(即抵抗观测粗差即其影响的能力)和经济性，并于1971年首先引入了准则矩阵的概念，提供了一种更为全面和密的精度标准。Grafarend对控制网优化设计问题作了广泛的理论研究，提出了现已为国际上所公认的四阶段优化分类概念，系统地应用了各种数学规划方法来解决控制网优化设计问题。在1972年Grafarend、Schaffin导出了具有均匀和各向同精度结构的准则矩阵，即著名的Taylor - Karman结构，简称TK结构。他们的工作为控制网的优化设计奠定了理论基础。 在GPS技术广泛应用于测量工作以前，控制网的优化设计理论与方法主要应用于传统的二角网、边角网、导线网以及变形监测网，其内容涉及:各种准则矩阵的构成，网形优化设计、观测权优化设计和附加观测值的优化设计，其目标函数主要有精度(包括整体精度与局部精度)、可靠性(内部可靠性和外部可靠性)、经济性、灵敏度等。采用的主要方法主要有解析法( Analytical Method)和机助模拟法(Computer Aided Simulative Method)，解析法设计的原理为在各阶段设计中根据固定参数及有关的约束条件，对待定参数求最优解，其函数模型可表示如下: 式中( )表示不等式 或不等式 或等式 求待定参数的方法主要是线性规划法，但由于实际的控制网约束条件复杂，不确定因素很多，很难计算出可行解，因此解析法主要用于小范围的精密工程控制网。 模拟法是指对于初步确定的网形与观测精度，模拟一组起始数据与观测值，输入计算机，然后根据平差原理计算出未知参数及其函数的精度，估算成本，或进一步计算观测值的可靠性、敏感性等信息，与预定的精度要求、成本约束、可靠性约束等相比较;根据计算机所提供的信息与设计者的经验，对控制网的基准、网形、观测精度等进行修正;然后重复计算，直至获得符合各项设计要求的较理想的设计方案。其工作流程可示如下: 显然，模拟法由于所得的最后结果会受到设计者经验的影响.只是一种相对的最优，而非严格最优方案，但是这种方法的数学模型不需建立。可以用于任何类型的设计优化问题，同时，它的最后结果一定是满足要求的、切实可行的近似最优方案。 1.2 GPS控制网的发展历史 近十年来，随着GPS技术的出现，由于其具有控制点间不需要相互通视、测量速度快、精度高、能全天侯作业等常规测量方法无法比拟的优点，在城市或工程控制网的测量中，它基本己取代了常规的测量方法，因此对GPS网的设计与优化进行研究具有较强的现实意义和广阔的应用前景。目前，人们对GPS网形设计主要形成了以外一些观点，一些人认为GPS网的图形矩阵A类似于水准网，而在水准网中的设计矩阵A与网的点位无关，也就是说GPS网的优化设计类似于相应的水准网的优化设计，也有文章认为GPS网的网形，即点位的分布， 则主要考虑点的观测条件及未来扩展的方便，观测权的优化设计对GPS网意义不大，也有一些人主张加测GPS基线以增加网的图形强度。但还未和常规控制网的优化设计一样，对GPS网的优化设计进行系统的研究，其主要原因在于和常规的测量方法相比较，GPS观测具有以下特点： 具有更为复杂的函数、随机模型，这些模型的先验信息在设计阶段难于获取和准确估计，即在实测前，各基线观测向量的误差(权逆阵)与模型误差一样属于非参数估计； 基线观测问题不受〔或几乎不受)通视条件限制，因而GPS网具有更为灵活，多样的布网形式，为优化设计提供了更为确实可行的条件。 这些特点决定了GPS网优化设计与实测方案难于完全吻合，即“最优”设计仅有理论上的意义而并非能在实测中完全准确实现，对GPS网的优化不仅对网形进行优化，而且还应对基线向量的观测时段、每时段的观测长度、地面点的几何精度因子(GDOP)等指标进行优选，因而导致了GPS网优化设计的复杂性。 1.3 GPS网形的特点 (1) GPS网的同步环与异步环 GPS同步观测是指使用 台GPS接收机，同时在相同的时段内连续跟踪接收相同的卫星组信号。通常称同步观测的时间段为时段或测段， 台GPS 接收机在一个时段内同步观测，可以由软件算得各点间的GPS基线向量(或简称GPS边)，这些GPS边可以构成一个完全图，图形中包括的GPS边的边数为: （1-3-1） 但其中仅有M-1 条是独立的GPS边，其余的都是非独立的GPS边，非独立的GPS边可以由独立的GPS 边推算得到。 当同步观测的GPS接收机数 时，独立的GPS边和非独立的GPS边将构成闭合的多边形环，称为多边形闭合环，同步闭合环的最少个数为: （1-3-2） 理论上，同步闭合环中各GPS边的坐标差之和（即闭合差）应为0，但因为有时各GPS接收机并不是严格同步和各种误差的影响，使得同步闭合环的闭合差并不等于零，只要同步闭合环的闭合差不超过限差，则认为此测段的GPS观测是合格的，但并不能说明GPS边的观测精度较高，也不能发现接收机的信号受到干扰而产生的粗差。 为了确定GPS观测效果的可靠性，有效地发现观测值的粗差，必须使GPS网的独立边构成一定的几何图形，这种几何图形可以由数条GPS独立边构成非同步多边形 （或称为非同步闭合环、异步环），如四边形、五边形…等，GPS网中的图形设计，就是根据所布设的GPS网的精度、可靠性及其它方面的要求，设计出由独立GPS边构成的多边形网（ 简称为GPS网）。 (2)GPS网中各要素之间的关系 GPS网形可以由图论来描述，从其组成要素上来看，它类似于GIS空间数据的构成要素，也是由点、线、面构成的，这里的点是指已选定位置的GPS点，线表示两点之间存在独立的GPS 边，面指由独立的GPS 边构成的异步环。GPS网形是靠一定数量的GPS 接收机在野外观测一定数量的时段而实现的，因此在设计网形、编制观测纲要时一定要建立这些量之间的关系。设 为GPS 点的个数， 为整个网中的独立边数， 为异步环的个数。在GPS测量规范中，有的要求网中异步环的边数小于某个给定的整数值，有的要求每个点上至少应设站两次，实际上这两个规定有其一致性。 （Ⅰ）当要求GPS网由m边形组成时，存在以下关系: （1-3-3） 其中： int( )表示取整数，mod ( )表示取余数，对于不同的 和m ，按上式计算出相应的 和 ，从上式可以看出，在根据测区范围和精度要求确定 和m以后，可以由上式求出需 要观测的独立的GPS边数 和 ，从而对布设GPS网的工作量有一个恰当的估计。 （Ⅱ）当要求每个点上至少设站 次时: 对于一个有 个点的控制网，如果在每个点上至少设站T次，则共应设站 次，如果用于观测的GPS接收机有M台，则需要观测的总时段数k为: （2-3-4） 此时，独立的GPS边和网中闭合环数为: （2-3-5） 1.4 GPS系统的应用前景 当初，设计GPS系统的主要目的是用于导航、手机情报等军事目的。但是，后来的应用开发表明，GPS系统不仅能够达到尚书目的，而且用GPS卫星发来的导航定位信号能够进行厘米级甚至毫米级精度的静态相对定位，米级至亚米级精度的动态定位，亚米级至厘米级精度的速度测量和毫微秒级精度的时间测量。应此GPS系统展现了极其广阔的应用前景。 其主要应用有： (1)在控制测量中的应用； (2)在精密工程测量及变形监测中的应用 (3)在航空摄影测量中的应用 (4)在线路勘测及隧道贯通测量中的应用 (5)在地形、地籍及房地产测量中的应用 (6)在海洋测绘中的应用 (7)在智能交通系统中的应用 (8)在地球动力学及地震研究中的应用 (9)在气象信息测量中的应用 (10)在航海航空导航中的应用 2 GPS网的精度和可靠性 2.1 GPS网平差的数学模型 大地测量控制网的优化设计是通过调整观测点的位置及观测方案来达到使目标函数最优的目的，为此，首先要建立起观测量、未知参数、观测量的权以及未知参数的方差、协方差矩阵之间的函数关系。 GPS接收机的初始观测量是瞬时载波相位差.而未知参数包括卫星坐标的改正量、地面点的坐标以及一些系统误差参数 (卫星钟差、接收机钟差、电离层、对流层的折射延迟以及整周模糊度等)，作为第一步，首先要解出基线向量，文献[2〕对GPS基线向量的解算已有详细叙述，本文不再述及。而以基线向量作为观测值的网平差的数学模型和随机模型正是本文所关心的。 “基线选择法”是GPS控制网平差最常用的方法，它通过选择独立基线来组成平差网形，分为两个过程，首先解算出同一时段的基线向量，然后从基线向量中选取独立的基线向量进行网平差。在网平差阶段，将地面点在的三维坐标当成待求的未知参数，从而得到以下的Gauss-Markov模型: (2-1-1) 这里A是由0，1组成的图形矩阵，X是待定点的近似坐标改正数向量，而△r则由下式给出 : （2-1-2） (2-1-3) 获得 (3-1-1)式基线向量方差的方式有两种，第一种方法是采用基线向量解算时计算出的方差和协方差阵作为验前协方差阵，第二种方法是在基线向量的经验公式。 的基础上采用验后最小范数二次无偏估计(MINQUE估计)得到，其中常数a、b分别为接收机的固定误差和比例误差系数，s为基线长度。 在GPS网的优化设计中，由于没有实际观测值，采用以下的经验公式作为基线向量各分量的近似精度: (2-1-4) 若某 GPS网有 n条基线，并不考虑同一基线各坐标分量之间的相关性，则所有基线分量的方差、协方差阵可表示为: （2-1-5） 其中diag 是一个由n个 的单位矩阵所组成的对角阵。 2.2 误差的传递与转换 GPS网精度估算是指在选定的基准条件下，通过误差传播律来近似计算待定点的方差和协方差、各点之间的相对中误差、或边长、方位角误差等。若以 表示GPS网中所有点的协方差阵，由于其它量 (如边长、方位角等)均为坐标未知数的函数，因此也可以用误差传播律来求坐标未知数函数的协方差( ),因此GPS网的精度评定主要是计算 精度估算与网平差的区别在于:精度估算中作为观侧值的基线向量是一个未知量，其先验的协方差未知，而在网平差中，基线向量的协方差阵己在基线解算中获得，其次，网平差中常数项可以准确地计算出，而在设计阶段，其值未知。在精度估算中，基线向量各分量的先验精度采用 ((3-1-4)式的计算方法获得，在网形设计阶段，我们最关心的是在独立坐标系下的精度，因此误差要经过以下一系列的计算和传播过程。 1 GPS网的三维平差 设任意两点i, j的GPS基线向量观测值为( ),它门是WGS-84坐标系中的空间直角坐标差，又设待定点在 WGS-84坐标系中的空间直角坐标为未知参数并记为: （2-2-1） 其中，( )为坐标近似值，( )为坐标平差值，( , , )为坐标改正数，容易得出GPS基线向量观测值( )与未知参数存在以下关系: (2-2-2) 式中，( )表示为( )的改正值，从而写成误差方程的形式: （2-2-3） 设某GPS网中共有m个点，观测了n条基线，则上式写成矩阵的形式为: （2-2-4） V: 3n 1阶矩阵，表示基线向量的改正数 A: 3n 3m阶矩阵，称为系数矩阵或图形矩阵，与网形有关 :3m 1阶矩阵，坐标未知数的改正数: L:常数项向量。 (2-2-4)式是不包含基准条件的误差方程，其系数矩阵是秩亏的，无法求解，对于三维控制网而言，本应包含三个位置基准，一个尺度基准，三个方位基准，在网平差时可以通过以下约束条件来获得网的基准: (1)固定约束 (即固定点的坐标值): (2)条件约束 (即通过使某些未知参数之间满足给定的条件); (3)权值约束 〔通过给某些观测量或未知量分配很大或很小的权)。 对于GPS网，作为观测值的GPS基线向量包含了尺度和方位信息，因此，在 WGS-84坐标系中平差时，可以取 GPS基线向量提供的尺度基准和方位基准，在此，我们取 b点的单点定位的三维坐标作为固定值，从而可以写出基准方程: （2-2-5） 写成矩阵形式为: （2-2-6） (2-2-4), (2-2-6)式联立，即为附加基准条件的GM模型，在最小二乘准则下求解得: (2-2-7) (2-2-8) 需要说明的是，在城市或大型工程控制网中，经常采用分级布网的方案，对于分级布设的GPS网，首级网一般采用只固定一个位置基准的无约束条件的三维平差，而次级网是在固定多个首级网点的情况下的有约束条件的三维平差，很明显，次级网的误差也包含首级网的误差，在次级网点的精度评定时，一定要考虑首级网的误差，具体的估算方法将在以后的章节中详细述及。 2 GPS 网空间直角坐标的协方差阵转换成大地坐标协方差阵 GPS网在空间直角坐标系下进行三维平差之后，需要在椭球下(WGS-84椭球、国家椭球或地方区域椭球) 将空间直角坐标( X Y Z ) 转换成大地坐标 ( B L H)，有以下熟知的公式: （2-2-9） 其中 为椭球的第一偏心率，N为卯酉面的曲率半径， 对 上式微分，得出大地坐标与空间直角坐标的微分关系: 3 GPS 网大地坐标协方差阵转换为高斯坐标协方差阵: GPS的大地坐标要变换成平面坐标，要通过高斯投影来得到，而高斯投影要先选择中央子午线，中央子午线可以按国家坐标系的 3度、6度带来选择，也可以选择任意的中央子午线，如在城市独立坐标系中，一般选择过城市中心区域的经线做为中央子午线。通过高斯投影的正算可以得到各网点的高斯坐标及其协方差阵。其协方差计算可以采用如下: 其中: M , N分别为子午面和卯酉面的曲率半径，经差 ( 为 中 央子午线经度),按误差传播公式，存在下式: 2.3 相对点位精度的合理评定 普通的( 绝对)点位精度是指控制点相对于一组起算数据(包括起算点的坐标、起算方位角和起算边长)的精度，而传统的相对点位精度评定方法是利用坐标差来确定两点之间的相对点位精度: (2-3-1) 这种相对点位精度的评定方法与绝对点位精度相比较，仅仅是改变了起算点的位置，而没有改变起算方位角和起算边长的位置，因而理论上是不严密的训，不能合理地评定两点之间的相对精度，而在实际工作中，有时需要根据不同位置的起算数据来分析控制网的精度。例如在工程控制网中，为确定离起算数据较远的控制点是否满足精度要求，有时并不要求这些点相对于原起算数据有很高的精度，而只要求它们相对于邻近某个点、某个方位和边长的相对精度，为了便于合理地评定相对点位精度，我们把平差基准分为坐标基准和方差基准。 1 坐标基准 在平差问题中，待估的未知参数往往不是观测值，例如在 GPS网平差中，观测值〔更准确地说应该是伪观测值)为基线向量，而取点的坐标为参数，如果没有足够的基准条件，这种坐标参数是无法确定的，这种起算数据称为平差问题的基准。为了便于和后面提出的方差基准区分，我们称此基准为坐标基准。GPS 测量的基准条件可统一表示成为下式: (2-3-2) 若在基准条件为 的条件下，则附基准条 件的GUASS-MARKOV模型表示为: (2-3-3) 在最小二乘准则下组成法方程:( 其中: K为附加的参数向量): (2-3-4) 上式中由于N秩亏，令 ， 对上面两个方程进行矩阵的初等变 换，可以证明N是满秩的，按文献的分块矩阵求逆公式，可得到(3-3-5) ，这与文献的附参数条件的间接平差计算结果是一致的。 = (2-3-5) 即： (2-3-6) (2-3-7) 在只固定一个基准点的情况下，上式可简化为： (2-3-8) (2-3-9) 但在城市控制网的首级网和工程控制网中，一般只取GPS基线的尺度基准，采用固定一个地面网的基准点和一个起始方位的无约束条件的三维平差，因此，在 WGS-84坐标系中对基线向量进行三维平差时应加入方位角信息。若以 为基 准点，以 的大地方位角为方位基准，由GPS三维坐标差求得的大地坐标方位角为: (2-3-10) 令 ， 对上式线性化，得 （2-3-11） 在加入方位角条件下，基准变为： （2-3-12） 为方便计算，可以将基准条件中的方位角条件作为伪观测值，取其权值为无穷大，在有m条基线的情况下，在附加第3 m + 1个条件方程(3-3-13) ，从而仍可以沿用(3-3-8)和(3-3-9)式: (2-3-13) 2 方差基准 在选定的坐标基准下，若同时默认起始数据的方差为零，即可得出重合于坐标基准的协方差基准，在此基准条件下得出的平差后网点坐标的协方差阵就是通常所说的相对于起算数据的绝对点位精度(绝对点位精度) 。如图，设4,7 点 为某隧道控制网的两个洞口点，设所选定的坐标基准条件 为固定1点即1-2方位，而我们所关心的并非洞口点的绝对点位精度，而是在固定7号点和7 -3 方位条件下(称为 )4号点的点位精度。在测量控制网中，我们关心的是平面坐标系F的点位精度，由于大地方位角和地面独立坐标系的起始方位角只差一个极小量的旋转角，因此在空间三维直角坐标系中固定大地方位角与在平面上固定坐标方位角是等价的，完全可以在二维的独立坐标系下对控制网进行协方差基准变换，若以 点为墓准点， 两点的方位角为起算方位角,(2-3-12)的三维基准条件与下面的二维基准条件是一致的: （2-3-14） 故在三维平差中，由三维基准条件与二维基准条件的一致性，将空间三维直角坐标下的坐标协方差阵转换成平面坐标协方差阵后.再沿用文献 列出的协方差基准的性质和转换公式，并由这些性质检验协方差基准的正确性，具体的公式在此不再推导。 3 在不同的方差基准下求相对点位梢度的实用公式 设在给定的坐标基准条件下，未知点的协方差阵经过三维平差以后，经过前已述及的一系列变换，得到了其高斯坐标的协因数阵 ，设新的方差基准条件 为固定 点及 之间的坐标方位角，取： （2-3-15） 式中 ， ， 。 条件 也可写为: （2-3-16） 如果求出所有点在新方差基准下的协因数阵，其计算量会相当庞大，也会占用较多的内存，而在实际工作中，我们只需要评定某两个点的相对点位精度，即某个点〔如K点)在新的方差基准下的点位精度，故可对式简化。 2.4 GPS控制网的可靠性 控制网的可靠性概念是为研究模型误差 (土要是粗差和系统误差)而提出来的，用来描述网本身发现某一模型误差的能力的指标称为网的内部可靠性，控制网抵抗模型误差的能力称为网的外部可靠性。 在工程实践中，人们往往倾向于注重控制网的精度，而忽略了控制网的可靠性，实际上可靠性与精度同样重要，而且_者不是等价的，高精度的控制网不一定有好的可靠性，而可靠性好的控制网也不一定具有较高的精度，因此在网的设计阶段一定要考虑到网的可靠性。对于GPS网而言，GPS网的原始观测值转换成为平差中的基线向量经过了若干次转换，也包含了较常规测量更多的系统误差源，在网形设计阶段充分考虑可靠性指标更有其必耍性。 和常规网相比，GPS基线向量网的数据的平差处理由两个阶段组成，一是对相位观测值进行基线向量解算，二是以基线向量为观测值进行网平差计算。第一阶段，可靠性研究侧重于如何通过残差分析剔除劣质相位观测值，包括消除小的周跳，确定好整周模糊度的偏差值等。第二阶段，主要研究以基线向量为观测值的可靠性问题。 1 网的可靠性矩阵 以GAUSS-MARKOV模型为例，在最小二乘准则下，未知参数的解为 （2-4-1） （2-4-2） （2-4-3） 若观测值中有一模型误差 ，则 对改正数的影响为 （2-4-4） 由式 ((3-4-3), (3-4-4)可以看出，任何一个模型误差一般对每一个改正数均有影响，其影响系数为矩阵R的元素，R的对角线上的元素 为 的模型误差 对 本身改正数 的影响系数，非对角线元素 表示 的模型误差 对 的改正数 的影响函数，因此，矩阵 R反映了模型误差对改正数的作用程度，而改正数是用来探测模型误差的主要依据，所以R矩阵包含了控制网的可靠性信息，称之为可靠性矩阵，对于常规的测量控制网和GPS网上式同样适用。 2 GPS网的内部可靠性 当某一观测值包含模型误差 式中H为系数矩阵，S为P维模型误差系数向量，则GM模型变为带模型误差的扩展模型: （2-4-5） 在最小二乘准则下，解出得： （2-4-6） 为了检查模型误差是否存在，提出以下的原假设和备选假设： 原假设: （2-4-7） 备选假设： （2-4-8） 在有多个粗差的情况下，能以一定的检验功效 ，通过显著水平为 的统计检验可发现的 (S为 的单位向量)方向可发现粗差的下界值为: （2-4-9） 式中 是 单 位 权中 误 差， 为 非中 心 化 参 数 。 针对GPS基线向量为观测值的G P S网，设基线向量的粗差向量的单位向量部分为: , , （2-4-10） 在网形设计阶段，我们讨论的是单基线的可靠性问题，其特点是各个单基线向量之间不相关，而一条基线的三个分量之间相关。因此 G P S基线向量观测值的权阵是一个分块对角阵。根据这一特点，对于GPS网，可以得出下式: （2-4-11） 式中， 和 分 别 是 权 阵和 阵的第 个 三 阶 主 子 块 。 由( 2 - 4 - 9 ) 式 可 求 得 某 一 基 中 式线向量的三个坐标分量方向上的内可靠性为: （2-4-12） 式中： （2-4-13） （2-4-15） （2-4-16） 上式中的 为第i 条基线向量三个坐标分量的多余分量主子块阵， 在无需 顾及方向的前提下，可按下式来求出单条基线向量可发现粗差的下界值: （2-4-17） 它是一种平均意义下的定义，即假设某条基线在空间各方向上可发现的粗差的下界值是相等的，因此单条基线向量可发现粗差的下界域也是一个球形域。若假设某基线向量的各分量是等权的，则 ( 2-4-12 )式可以变为: （2-4-18） 为了比较各种不同精度的观测值之间发现粗差能力的差别，将上式中的观测精度口 去掉，得到的无量纲的内部可靠性指标与文献的内部可靠性指标完全等价，这说明若将基线观测分量作为观测值，可以得到与常规网一样的内部可靠性指标。 3 GPS网的外部可靠性 外可靠性是指不可发现的最大模型误差对平差结果的影响，一般将这种模型误差对平差未知数的影响长度定义为外部可靠性指标。在 ( 2-4-5 )的数学模型下，外部可靠性可表述如下: （2-4-19） 式 中 : 类似于内部可靠性的讨论,可求得一个基线向量在三个坐标分量方向上的外部可靠性指标为: （2-4-20） 式中， 分别是该基线向量的权阵中的对角元素， 同样，定义一个基线向量平均意义下的外部可靠性指标为: 4 控制网可靠性的一些性质 矩阵 为可靠性矩阵，内部可靠性和外部可靠性都由其决定， 通过对它的分析，可以得到以下一些性质: (1)tr(R)=r二其中r为多余观测数，对于有m个点，n条基线的GPS网，则 ； (2)R的对角元素满足: ； (3)对于某个观测值 而言，其可靠性权阵 与自身的精度权阵 成反比， 而与其它观测值 的权阵 成正比； (4)增加新观测值，既能提高其它原有观测值的可靠性，又能提高网的总体可靠性； (5)衡量网的可靠性的准则，不仅取决于多余观测数值，也取决于各观测值的可靠性权阵分配是否均匀。 3 GPS控制网优化设计 3.1 GPS网优化设计的数学模型 (1) 常规网的优化设计模型 我们做任何一项工作时，总是希望在所有可行的方案中选择一种某种意义上最优的方案，这就是最优化设计问题，数学上的优化设计可以用以下的函数模型来表示: （3-1） 式中 为设计变量，可以由设计人员调整，其不同的取值 表示 不同的设计方案，是优化设计问题中最终要确定的变量。 为设计变量的函数，称为目标函数，其函数值的大小决定了设计方案的好坏，也称为优化设计的准则矩阵; 不等式约束条件， 为 等式约束条件，统称为约束条件。约束条件是为了确保设计方案能够满足使用上的要求而建立的，以使设计方案最终是确实可行的。满足约束条件的所有设计变量组成的集合称为可行域，可行域内的每一个点所代表的方案都是可行的。 数学最优模型应用到测量控制网中，用数学语言可以描述如下， （3-2） 其中 分别为精度、可靠性、经济指标， 分别为图形矩阵和权矩阵， 分别为未知数个数和观测值个数， 限制条件函数， 为常数。(3 - 2)即顾及精度、可靠性、经济性等综合指标的大地控制网优化设计的数学模型。按照 提出的，目前国际上广泛采用的分类方法，控制网的优化设计分为以下几个阶段: (1)零阶段设计问题(或称基准选择问题)，对一个已知图形结构和观测 计划 项目进度计划表范例计划下载计划下载计划下载课程教学计划下载 的自由网，为控制网点的坐标及其方差阵选择一个最优的坐标系。即在己知设计矩阵A和观测值权阵P 的条件下，确定网点的坐标向量X及其协因数阵 ， 使X的目标函数达到极值，因此，零阶段设计问题就是一个平差问题； (2) 一阶段设计问题(或称图形设计问题)，即在已知观测值的权阵P的条件下，确定设计矩阵A，使网中某些元素的精度达到预定值或最高精度，或者使坐标的协因数阵最佳逼近某个给定的矩阵 (准则矩阵)； (3)二阶段设计问题(或称观测值权的分配问题),即已知设计矩阵 A，确定观测值的权阵P，使网的精度(整体精度或局部精度)达到预定值或最高精度； (4)三阶段设计(或称为网的改造或加密方案的设计问题)，通过增加新点和新的观测值，以改善网的质量。在给定的改善质量前提下，使改造测量工作量最小。或者在改造费用一定的条件下，使改造方案的效果最佳。 以上各阶段的设计既相互区别，义是相互联系的，如三阶段设计问题可以看成是一阶段和二阶段设计问题的混合，各阶段设计问题的己 知量与设计量之间的变换关系可以用下表来表示: 类 别 已知量 设计变量 零阶段设计 一阶段设计 二阶段设计 三阶段设计 2 GPS网的优化设计 具体到GPS控制网,其数学模型与常规大地测量控制网一样采用(3- 2)式，其优化设计步骤也可以按上述的四个阶段逐步优化，但是，应该和常规网有所区别: GPS网平差广泛采用的方法是在 WGS-84空间直角坐标系中进行无约束条件的平差，然后将平差后的坐标及其协因数阵转换到高斯平面上，再在平面上建立坐标转换模型，根据GPS网和地面网的公共点的坐标及坐标约束条件。求解转换参数，然后将GPS点坐标及相关的精度信息，根据转换模型一次性转换到国家或地方坐标系下。若只固定一个点和一个已知方位，这称为无约束条件的GPS网平差，若固定点的个数在两个以上，则称为约束条件的GPS网平差。前者最大限度地保留了GPS技术高精度的特点，但地面公共点的坐标发生了变化，后者虽然可保持约束点的坐标不变，但约束平差后G P S网的尺度保持了常规网的尺度，其度可能是虚假的，特别是当约束点中存在较大误差或显著位移时。选择哪种平差法应根据建网的目的来考虑，当对城市控制网进行新建或大规模地扩建时，应选择无约束条件的网平差模型，当对原有地面网局部加密时。应采用约束转换。而对于工程控制网，大多采用无约束条件的网平差，以保持网的较高精度，如隧道网，它一般选择洞口点作为基准点，桥梁控制网，则选择其主轴线作为基准方向，轴线的某一端点作为基准点 ，对于变形监测网一般采用自由网平差或拟稳基准。故对于基准条件的选取受网的用途和目的的限制，其选择的自由度不大，故对于零阶段优化设计没有多大的实际意义，但要注意基准的选择引起GPS网点位误差变化的规律。 控制网优化设计的方法有许多种，目前广泛采用的方法主要是解析法(AnalyticalMethod )和机助模拟法(Computer Aided Simulative Method ) ，解析法是通过建立优化设计的数学模型，包括目标函数和约束条件.选择一种优算法,求出问 题的最优解;机助模拟法则是将电子计算机的计算功能和判别能力同设计者的知识和经验结合起来，通过一个凭经验拟订的初始观测方案，进行分析计算，求出各项质量指标，并对设计方案进行不断的修改，直到满足设计者的要求为止的一种方法。 对于GPS控制网，由于解析法优化设计模型比较复杂，难以解出可行解，目前编制出的控制网的优化设计软件一般都采用模拟法，需要指出的是模拟法的最终结果相对于解析法而言，在严格的数学意义上并非最优，而只是一种相对最优，下面将模拟一小型GPS网，在顾及精度和可靠性的准则下，对GPS网的优化设计进行探讨。 3.2 一阶段优化设计 1点位的移动及异步环的数盆对精度和可靠性的影响 一阶段设计即GPS网形设计，主要是指在控制点的个数。多余观测数大致确定的情况下，设计出最优的网形，以使整个控制网具有较好的精度和可靠性。 在常规测量中，点位的位置移动对精度和可靠性的影响较大，主要原因是点位的移动将改变传递三角形的角度，从而会改变点位误差椭圆的大小和方向，如三角网布网时就规定传递三角形的角度一般不得小于30度。而GPS观测类似于水准测量，只要不改变点间的连接方式，其误差方程就不会改变(即图形矩阵A不变)，点位的移动对权阵p会不会有影响呢?GPS基线向量的精度估算公式(3-2 -1)对基线长度微分得: (3-2-1) 在上式中，取s=20km ,ds=1km ,b=2ppm，计算出 ，在进行精 度估算时其变化量完全可以忽略，对于采用厂商提供的类似的标称精度公式: 来估算基线标准差的情况，同样可以证明由点位移动而引起基线标准差的变化量将会更小,所以，在不改变GPS边的连接关系的情况下，仅对点的平面位置进行小范围的调节是不会改变网的精度和可靠性的. 控制网网点的个数固定的情况下，异步环个数的改变是否会对网的精度和可靠性产生影响。规范中没有明确规定异步环中独立边长的数目，但在城市控制网的技术设计中一般说明异步环最多不能超过6条边，下面将以实例来分析:如图3-1 、图3-2为某模拟的GPS控制网，采用单频接收机观测(10m m+2ppm )， 平均边长约20公里，在相同的基准条件下(以228号点为固定点，以228-227的方向为起始方向)，而异步环基线向量数不同的情况下，得出的各点的误差如表4 -1 所示. 图3-1 图3-2 点名 网形一 网形二 Mx(cm) My(cm) M（cm） Mx(cm) My(cm) M（cm） 211 1.59 1.81 2.41 1.69 1.92 2.56 213 0.85 0.97 1.29 0.96 1.09 1.46 216 1.01 1. 15 1.53 1.28 1.46 1.94 218 1.84 2.09 2.78 1.86 2.12 2.82 220 1.93 2.19 2.92 1.96 2.23 2.97 226 1.41 1.60 2.14 1.83 2.07 2.76 227 0.73 0.83 1.11 0.79 0.90 1.19 249 1.56 1.77 2.36 1.92 2.17 2.90 244 1.41 1.60 2.14 2.06 2.34 3.12 245 1.01 1.15 1.53 1.41 1.60 2.1 3 247 1.00 1 .14 1.52 1.13 1.28 1.70 228 0.00 0.00 0.00 0.00 0.00 0.00 2 GPS网平差基本观测A的选取对精度的影响 作为网平差观测值的基线向量既包括独立的基线向量，也包括和独立基线向量相关的基线向量，原则上都可以参与平差计算，基本观测量的选取是否会对网点的精度产生影响?设整个网有n个观测值L，其中有 个 为函数独立的观测值， 个 与 函数相关的观测值，从而使平差模型变为具有相关观测值的平差问题，观测值L表示为: （3-3） 观测值的协因数阵为： (3-4) 误差方程 在最小二乘准则下: (3-5) (3-6) 单位权方差的无偏估计量为: (3-7) f为自由度，t 为必要观测数: （3-8） 将广义逆分解为特殊逆: (3-9) 如果只将独立观测只参与平差，则误差方程的解为: (3-10) 图3-3 对于高精度的 GPS网，基线解算一般采用高级软件，如 GAMIT软件等，一个同步图形的函数独立的基线向量，其观测协方差阵给出是完全的，严密的，故有 ， 而对于一般的商用GPS解算软件，它只能提供单一基线的协方差阵，对于如下同步基线向量，以下关系不一定成立: 故在选取独立基线时，若选取不同的独立基线，虽然图形矩阵A是等价的，但其观测值的协方差阵发生了变化，本文也提出了基线的优选问题。故提出选取使异步环闭合差最小的独立基线作为观测值平差，以此来提高平差后的点位精度，但选择独立观测值的算法很繁杂，对验后的精度改善也不是明显，还是沿用常规的按网形选择独立基线的做法。最值得注意的是在不增加总的观测值的情况下，而只在参加平差的独立基线向量中加入同步观测值，对精度是没有影响的，这一点由( 3-5 )，( 3-9 )，(3-10)式推导出。 3.3 三阶段优化设计 1增加观测值对精度、可靠性的影响 原有观测方案的误差方程: (3-12) 增加m个权为 的观测值 得到附加观测值的误差方程: （3-13） 将以上两式联立求解，即得到附加观测值后的网的精度矩阵 (3-14) 利用矩阵反演公式，将(4-13)式展开为: (3-15) 若用 表示网的整体精度，则增加观测值后网的整体精度变化为: (3-16) 同理也可以推导出增加观测值后的局部精度指标的变化量。 （3-17） 增加观测值后，(3-4-2)式的可靠性矩阵变化为: （3-18 在改造方案的优化设计中，自然是要优先加测那些对精度影响较大的观测值，以达到用最少数量的附加观测值来满足精度要求。由此，可以按以下的步骤来进行第三阶段的优化设计: (1)拟订附加的所有可能观测值及其观测值的精度，并给出改善精度的基本条件: (2)所有可能的附加观测值对精度、可靠性的影响值 或 (3)将影响最大的观测值加入网中，并计算新的精度矩阵; (4)检查精度指标是否满足要求，若满足要求则设计结束。 2 修改观测方案的实时算法 所谓实时计算方法就是利用前一次求得的结果来计算当前所需要的结果。增加一组观测值、 删除一组观测值和改变一组观测值的权是优化设计最常用的修改设计方案的手段，如果重新组成误差方程、重新计算未知数的改正数和坐标未知参数的协方差阵，必定会消耗大量的时间和占用更多的内存。在一般情况下，修改观测值的数目大大小于未知参数的数目，故可以采用实时算法。公式(4 -15)给出了增加观测值的算法，当减少观测值时，修改符号后该公式同样适用，该式可写成统一的形式: (3-19) 若是改变某些观测值的权，设 为改变权的观测值 所对应的 误差方程系数，则 的新权阵为: 若 将 在观测值的最后，可以得到观测值的新权阵: (3-20) 新权阵下的协因数阵为: (3-21) 利用矩阵反演公式得: (3-22) 在GPS网中，误差方程的系数矩阵的元素为分块单位对角阵，因此(3-19) 、 (3-21)可以大大得到简化，在设计阶段，不考虑各基线向量之间的相关性和每条基线的各分量之间的相关性，由(3-18),(3-20)式可以推导出附加基线向量只改变与基线相关联点的点位精度，而不会影响其它点位的精度。以图4-2为例。现拟增加基线244-245，45-228，16-226，228-211中的某两条，按实时算法估算增加基线对点位糟度的影响见表3-2. 表3-2 实时算法估算增加基线对点位糟度的影响 增加 基线 左端点 右端点 △Dα (cm2) △mx (cm) △my (cm) △mz (cm) △mx (cm) △my (cm) △mz (cm) 211-228 0.10 0.11 0.15 0.15 0.17 0.23 0.12 244-245 0.49 0.56 0.74 0.08 0.09 0.12 0.31 245-228 0.36 0.41 0.54 0.00 0.00 0.00 0.19 216-226 0.06 0.08 0.10 0.34 0.37 0.50 0.15 结论 本文针对城市、工程GPS控制网设计和控制网的精度和可靠性进行了探讨，在借鉴己有研究成果的基础之上，以GPS网的精度和可靠性为目标函数，对GPS网优化设计的各阶段进行了分析和例证;采用图论的方法生成GPS网形，并利用最短路径的方法生成观测方案; 尽管本论文已经完成，但是还还有不少关于GPS网方面的问题需要解决。通过这次的论文，我从中学到了很多知识同时也发现了自身存在的很多不足之处。此刻的我们即将迈入新的工作环境，通过这次的论文我们都从中学到很多宝贵的经验，我相信接下来在工作岗位中，我们应该会认真对待，努力补充自己。 参考文献 [1] 余学祥等.GPS 网平面坐标转换的约束平差.工程勘察1999 (3) [2] 张文卿,基于最优化区域性椭球的人城市GPS控制网数据处理及专用软件系 统.同济大学硕士论文.1999.2 [3] 刘人杰，施一民等.全球定位系统( GPS )的原理及数据处理.同济大学出版社.1996.3 [4] 周秋生,测量控制网优化设计.测绘出版社.1992.6 [5] 邹伟鸿,具有精度与可靠性准则的GPS网优化设计. GPS卫星定位的应用与数据处理.同济大学出版社.1994.9 [6] 于来法，论GPS 精密工程测量技术的优选.测绘通报.1994(3) [7] 李德仁,误差处理与可靠性理论.测绘出版社.1988 [8] 周江文，杨元喜等. 抗差估计论文集.测绘出版社.1992 [9] 陶本藻,关于高精度GPS 网平差基本观测量的选取问题.测绘通报.1999 [10] 吴素芹,GPS网可靠性强度分析.武汉测绘科技大学学报.1992(3) [11] 赵长胜,控制网设计与数据处理.北京科学教育出版社,2001 [12] 李征航、黄劲松，《GPS测量原理及应用》.武汉大学出版社2005 [13] Shanlong Kuang. A Strategy for GPS survey Planning:Choice of Optimum Baselines.Vol.54.No.4.1994.pp.187-201 [14] M. UNGUENDOLI.A Rational Approach to Use of a Large Number of GPS Receivers. Bulletin Geodesique. Vol.64.1990.pp.303-312 致谢 时光茬荐,岁月如梭,四年的大学生活一转眼就过去了，而在这四年中收获的一切将是我人生中最为宝贵的财富之一。 首先，我要感谢我的论文指导老师陈芳讲师。本文的工作自始至终都是在陈老师的关怀指导下完成的。从论文开题、资料收集、到论文的撰写与定稿，无不浸透着陈老师的关怀与智慧。陈老师豁达的人生态度、严谨的治学精神让我终身难忘。在今后的学习深造的问题上，陈老师都给予了我很好的指点，在此论文完成之际，谨向陈老师表示崇高的敬意和诚挚的谢意！ 感谢安徽理工大学地球与环境学院测绘与空间信息系的所有老师，感谢余学祥教授、杨世清教授、郑礼全副教授、王列平副教授，徐良骥感谢你们在平时教学工作中给予我的启发和指导！ 在论文的完成过程中，得到了一些同学在学习上的有益指导和生活上的热心帮助，在此表示感谢！ 感谢马翼、舒扬、姚锟、欧林峰等室友与我一起度过的快乐时光！ 感谢参与论文答辩会的各位老师，感谢你们在百忙之中给予的有益指导！ 感谢所有一直关心我、帮助我的亲人和朋友! 感谢论文引用参考文献中提及或未提及的著、译和作者! 谨以此文献给含辛茹苦养育我成人的父母和所有关心我的人们!