(整理)测量平差教案

精品文档精品文档第一章绪论第一节观测误差一、观测值中为什么存在观测误差？观测条件对观测成果产生影响，不可避免产生观测误差。有观测就有误差的结论。二、观测误差的计算给出观测误差计算的纯量表达式和矩阵表达式。三、观测误差的分类及其处理1、分类给出误差分类的表达式，粗差、系统误差和偶然误差的定义。结合测角、测距和水准测量的全过程，让学生分析哪些因素引起的误差属于粗差，那些哪些因素引起的误差属于系统误差，那些哪些因素引起的误差属于偶然误差。2、处理总结粗差、系统误差和偶然误差的处理方法，让学生举例说明测量上哪些操作是为了消除系统误差影响的，那些计算改正为了消除系统误差影响的。四、测量平差的任务根据一系列含有观测误差的观测值求待定量的最佳估值。第二节测量平差学科的研究对象研究对象为含有观测误差的各类观测值。举例说明。第三节测量平差的简史和发展一、测量平差理论的发展1、经典平差理论的发展主要介绍高斯创立最小二乘原理和马尔可夫创立高斯-马尔可夫平差模型的历史背景和过程。2、近代平差理论的发展主要介绍二十世纪四十年代以后出现的近代平差理论，结合导线网平差和我国南极考察、建站，重点介绍方差分量估计和秩亏网平差的理论、方法及其用途。二、平差计算方法的发展1、手算阶段2、半自动平差阶段3、全自动平差阶段第四节测量平差的任务和内容一、任务讲授测量平差的基本理论和基本方法，为进一步学习和研究测量平差打下深入的基础。二、内容课本各章的内容。小结：本节介绍了观测条件的定义，观测条件与观测误差的关系，观测误差的定义、处理，以及测量平差的发展概况。第二章误差分布与精度指标第一节正态分布一、一维正态分布绘一维正态分布图，列出分布函数，讲解，强调两个分布参数的含义。二、n维正态分布讲解绘n维正态分布图，列出分布函数，讲解，强调两个分布参数的含义。第二节偶然误差的规律性一、偶然误差分布1、描述误差分布的三种方法（1）列表法（通过实例列表讲解）（2）绘图法（通过实例绘图讲解）（3）密度函数法（通过实例绘图讲解）二、偶然误差的分布特性（1）在一定的观测条件下，误差的绝对值不会超过一定的限值。（界限性）（2）绝对值较小的误差比绝对值较大的误差出现的概率要大。（小误差占优性）（3）绝对值相等的正负误差出现的概率相等。（对称性）三、两个重要概念（1）由偶然误差的界限性，可以依据观测条件来确定误差限值（2）由偶然误差的对称性知观测量的期望值就是其真值。小结：偶然误差有其统计规律，研究偶然误差的分布规律是为了更好的研究偶然误差的处理问题。第三节衡量精度的指标；第四节精度、准确度与精确度；第五节测量不确定度一、精密度指标(一)观测量的精密度指标1、观测条件与精密度配合误差分布曲线讲解精密度的定义和观测条件与精密度的关系。2、几种常用的精密度指标方差与标准差推导相应公式，给出其估值公式，讲解应用实例极限误差分析误差出现在某一范围内的概率的大小，给出极限误差定义公式相对误差给出相对精度的定义，用实例讲解其应用范围。平均误差与或然误差给出平均误差和或然误差的定义，讲解其在国际上应用的范围和地区，以及其与中误差的关系。(二)观测向量的精度指标1、n维随机向量的方差阵导出n维随机向量的方差阵表达形式，指出该阵是对称矩阵，并讲解矩阵中各元素的含义，同时给出当n维随机向量中各随机变量不相关时的矩阵形式。2、两随机向量的互协方差阵导出两个随机向量互协方差阵表达形式，并讲解矩阵中各元素的含义，同时给出当维随机向量不相关时的矩阵形式。二、准确度和精确度指标分别给出准确度和精确度的定义，及其数值指标，绘图讲解其几何意义。三、测量不确定度给出测量数据的不确定性、不确定度的概念，可测不确定度的计算方法，不可测不确定度的估计方法。小结：精度指标分为精密度指标、准确度指标和精确度指标三种，观测成果的质量应用精确度指标衡量，精密度指标中的方差、极限误差、相对误差几个指标应重点掌握。第三章协方差传播律及权第一节数学期望的传播律E(C)=C;E(CX)=CE(X);E(X+X+•••+X)二E(X)+E(X)+•••+E(X);12n12n当Xi相互独立时(i=1,2,・・・,)n,E(X,X，…,X)二E(X)E(X)…E(X)12n12n第二节协方差传播律协方差传播律是观测值(向量)与其函数(向量)之间精度传递的规律。一误差的传递1、线性函数误差的传递Y=fx+fx+...+fx+f1122nn0A二fA+fA+...+fA1x12x2nxn推导上述公式，讲解式中符号的含义2、非线性函数误差的传递Y=fClX2…TA二fA+fA+...+fA1x12x2nxn推导上述公式，讲解式中符号的含义3、函数向量误差的传递Y=FX+F0Y=F(X)y=fax讲解式中符号的含义，强调矩阵表达式与纯量表达式之间的相互表式二、协方差的传递1、基本公式函数向量Y=F(X)Z=K(X)其误差向量为y=faxZ=KAX则随机向量与其函数向量间的方差传递公式为

D=FDFtYXD=KDKZXD=FDKYZXD=KDFtZYX证明第一、第三式，并说明同理可证二、四式。=FDF=FDFtX=f2G2+f2G2+…+f2G11221讲解式中符号的含义，说明公式应用的条件，强调公式的重要性。3、分块向量函数向量的方差传递Z=t+r,1「X一t,1Yr,1「DDXXYDZ=t,tt,rDDYXXr,tr,r证明上式，对阵中元素加以说明，给出两向量不相关时该矩阵的形式。通过五个典型例题的讲解说明方差-协方差传播公式的应用方法和计算中需注意的问题。小结：协方差传播律是观测值（向量）与其函数（向量）之间精度传递的规律，用其解决观测值函数（向量）的精度评定问题。本节重点是利用协方差传播律解题的方法和步骤，以及只有一个观测值函数，且观测值之间不相关时的协方差传播公式的应用。第三节协方差传播律的应用1、水准测量的精度绘制具有N个测站的水准高差示意图，应用协方差传播公式导出高差中误差计算公式:&=JNGh站进一步导出S公里观测高差的中误差计算公式：&=品&hkm举例说明公式的应用。2、同精度独立观测值的算数平均值的精度由算术平均值公式，应用协方差传播公式导出其中误差计算公式x

x举例说明公式的应用。3、若干独立误差的联合影响A=A+A+…+Az12nO2=O2+O2+•…+O2Z12n即观测结果的方差，等于各独立误差所对应的方差之和。4、平面控制点的点位精度绘支导线略图，求未知点点位中误差，用两种方法求解。解法一：（1）、列函数式（2）线性化（3）应用协方差传播公式计算坐标方差（4）计算点位方差解法二：利用纵向方差和横向方差进行计算。小结：本节的重点内容为水准测量高差和同精度独立观测算数平均值的精度计算问题，应熟记计算公式，能熟练应用公式进行相关计算。第四节权与定权的常用方法一、权的定义权是衡量各观测值在平差结果中应起作用大小的数值。P=—oC=1,2,...,n）iO2iPi为观测值Li的权，O2是可以任意选定的比例常数。ii0观测值的权与观测值的方差成反比。二、单位权方差权的作用是衡量观测值的相对精度，称其为相对精度指标。确定一组权时，只能用同一个O令Oi=O0，则得：O2O2O2O=1O2i0上式说明O2是单位权（权为1）观测值的方差，简称为单位权方差。凡是方差等于O:的观测值，其权必等于1。权为1的观测值，称为单位权观测值。无论O2取何值，权之间的比例关系不变。0举例（例〔1〕、例〔2〕）讲解。三、测量中常用的定权方法1、水准测量的权（1）、用测站数定权（山地、起伏较大的丘陵）利用用测站数计算高差中误差的公式和权的定义式导出利用测站数定权的公式。解释式中符号的含义。（2）、用路线长度定权（平地）利用用路线长度计算高差中误差的公式和权的定义式导出利用路线长度定权的公式hS解释式中符号的含义。举例（例〔3〕、例〔4〕、例〔5〕）讲解。2、距离量测的权距离长度可通过钢尺丈量或测距仪测距得到。下面分别讨论两种情况下的定权方法1）钢尺量距的权解释式中符号的含义。2）测距仪测距的权P=^0s6SO=O+oS-10-6s标称标称解释式中符号的含义。3、等精度观测算术平均值的权利用等精度独立观测值算术平均值的方差计算公式和权的定义式导出利用观测次数定权的公式n牛C说明公式中符号的含义。小结：权是用来衡量观测成果的相对精度的，单位权方差可以根据计算方便任意选定，但观测值之间的比例关系不变。水准测量的权与测站数或路线长度成反比；钢尺量测的权与距离长度成反比，光电测距的权用定义式计算，其中测距方差由固定误差和比例误差两项组成；等精度算术平均值的权与观测次数成正比。应熟记定权公式，明确公式中各符号的含义，掌握利用公式解题的方法。第五节协因数和协因数传播律一、协因数

定义协因数Q=CC2=1/pl=1,2,...,nii厂ol权可表示为方差和标准差可表式为定义协因数Q=CC2=1/pl=1,2,...,nii厂ol权可表示为方差和标准差可表式为=C2•Q,0iiofQ0ll、协因数阵1、n维随机向量X的协因数阵定义互协因数：Q=C/6ljlj0利用方差协方差与协因数弧协因数的关系导出协因数阵C2Q11对称QQ1222Q]Q1n2nQnn上式矩阵中，Q=Q。当Qij=0(iMj)时，则Xi和Xj互不相关。ljjl2、分块向量的协因数阵_X—Z=Q=「QQXYYQ=「QQXYYX式中，Qx、Qy分别为X、Y向量的自协因数阵，而Qxy、Qyx分别为X向量关于Y向量的互协因数阵,Qxy与Qyx互为转置。当Qxy等于零时，表示X、Y互不相关。三、权阵「QQ...Q--1「PP..P1Q11Q12...Q1nP11P12..P1n21222n=21222nQQ...QPP..Pn1m2nnn1n2nnP=Q-1=观测值的权一般要通过对权阵求逆得到协因数阵，再利用权与协因数的倒数关系求权。当权阵为对角阵时，P=P。lll举例(例〔1〕例〔2〕讲解、分析四协因数传播律将协方差传播公式乘以+，并顾及Q=+D,Q=+D，即可得到观测向量X与其函数向量qY002YX㊁XY、Z之间的协因数传播公式。列出相应公式，以及只有一个函数，且观测值之间不相关时的协因数传播公式。举例（例〔3〕例〔4〕讲解、分析小结：权与协因数互为倒数关系，权阵与协因数阵互为逆阵关系，一般情况下给了观测值的权阵求观测值的权要先求权阵的逆阵得到其协因数阵，再利用权与协因数的关系求权；协因数传播律与协方差传播律公式相仿，只记住其中一套公式，再记住协因数阵与协方差阵的关系即可。第六节由真误差计算中误差及其实际应用一、利用不同精度的真误差计算单位权中误差的基本公式利用协因数传播律导出利用不同精度的真误差计算单位权中误差的基本公式oo-n~二、由真误差计算中误差的实际应用1、由三角形闭合差计算测角中误差利用协因数传播律导出由三角形闭合差计算测角中误差的公式■Ww]-矿说明公式的不严密性。2、利用双观测列之差求中误差（1）求单位权中误差利用协因数传播律导出利用双观测列之差求单位权中误差的公式不等精度观测Tpdd'r—I—'■・o7-2n-等精度观测0o:0o:[ddl说明公式中符号的含义。

2）求双观测列单次观测的中误差（3）求双观测列平均值的中误差利用协因数传播律导出相应公式不等精度观测等精度观测G_L等精度观测G_Lirr=G02PiQ=+b丄L.0\2i水准测量双观测平差应用例题小结：本节重点是利用双观测之差计算中误差的公式及其应用，该公式在测量中应用广泛，应重点掌握。第七节系统误差的传播一、观测值的系统误差与综合误差的方差1、观测值的系统误差——偏差导出偏差表达公式2、观测值的综合误差方差——可靠性如果系统误差部分是偶然中误差部分的三分之一或更小时，则可将系统误差的影响忽略不计。二、系统误差的传播导出传播公式三、系统误差与偶然误差的联合传播导出传播公式小结：了解系统误差的传播规律。第四章平差数学模型与最小二乘原理第一节测量平差概述一、测量控制网简介高程控制网（水准网或三角高程网）包括闭合水准网和符合水准网。绘出三组不同网形的水准网。网中元素：已知高程点，未知高程点和观测高差。2.平面控制网（1）三角网根据观测量的不同，三角网分为测角三角网、测边三角网和边角同测三角网。1）测角三角网包括独立三角网和符合三角网。绘出一组不同网形的三角网。网中元素：已知点，未知点和观测角度。2）测边三角网包括独立测边网和符合测边网。绘出一组不同网形的测边网。网中元素：已知点，未知点和观测边长。3）边角三角网包括独立边角网和符合边角网。绘出一组不同网形的边角网。网中元素：已知点，未知点，观测角度和边长。（2）导线网包括独立导线网和符合导线网。绘出一组不同网形的边角网。网中元素：已知点，未知点，观测角度和边长。还有三维网、GPS控制网、航测控制网、工程专用网等将在后续相应课程中介绍。二、必要起算数据确定几何（物理）图形的位置所必须具有的已知数据水准网（三角高程网）：一个已知点高程测站平差：一个已知方位测角网：一个已知点坐标，一个相邻已知方位，一个相邻已知边长或两个相邻点坐标。测边网和边角网：一个已知点坐标，一个相邻已知方位。各种控制网中少于等于必要起算数据的控制网成为独立网，多于必要起算数据的控制网成为非独立网或附合网。三、必要观测及其数目的确定确定几何、物理模型的形状、大小所必须进行的观测称为必要观测，其符号用符号t表示。高程网：t二p-q-丄测站平差：t二p-q-丄必要起算数据测角网：t=2p-q-4

测边网和边角网：t=2p-q-3P：总点数或总方向数（测站平差）；a：多余起算数据数必要起算数据之外的起算数据四、多余观测及其数目的确定必要观测之外的观测称为多余观测，其数目用符号r表示多余观测数=观测总数一必要观测数（r=n-1）五、必要观测和多余观测数目计算练习计算图3-1至图3-7的必要观测数和多余观测数。小结：本节介绍了测量控制网的类型，和各类控制网中应具备的必要起算元素，必要观测元素，应重点掌握必要观测元素数和多余观测元素数的计算。第二节函数模型FG)=0r,1~~FG)=0r,1~~L=FIXn,111,1丿2、间接（参数）平差法3、附有参数的条件平差法c,14、附有限制条件的间接（参数）平差法u,1u,1丿s,1用简单控制网图形举例说明。第三节函数模型的线性化设~~F=FIL,X

c,11n,1u,1丿用泰勒公式导出F的线性形式为F二FCX0)+AA+B~根据上述函数模型线性化过程，可将各种平差方法的函数模型线性化1、条件平差法AA-W二0式中A=-of,W=—F（L）L2、间接平差法

式中3、附有参数的条件平差法式中B式中3、附有参数的条件平差法式中B=-SFdXX0A=Bx-1,l=L—£o),AA+Bx-W二0W二—FC,Xo),4、附有限制条件的间接平差法A二B~—lc~—W-0x式中cs,ucs,udXX0，Wx第四节测量平差的数学模型1、各种平差方法的随机模型D=Q20n,nD=Q20n,nQn,n0n,n2、各种平差方法的数学模型各种平差方法函数模型的线性形式分别与平差的随机模型联立，即为相应平差方法的数学模型。小结：本次课所讲的各种平差方法的函数模型均能建立各观测值之间的函数关系式，正确建立这种关系式，是正确求得观测值最可靠结果的前提。第五节参数估计与最小二乘原理一、引例已知平面三角形三内角应满足或式中L+L+一、引例已知平面三角形三内角应满足或式中L+L+12L~—18003=0A+A+A—W=0123=1800—(L+L+L)111上方程中有三个未知数，是相容方程，只能在某一准则下求得式中未知数的估值。二、最小二乘准则:①=AtD—iA=min顾及方差阵与权阵的关系，并用△的估值V代替△又可得①二VtPV二min观测量真值向量的估值公式为：L=L+V式中L称为观测向量的“最或然值”向量或“观测值的平差值”向量；V称为改正数向量。三、最小二乘估计根据最小二乘准则进行的估计称为最小二乘估计，按此准则求得一组估值的过程，称为最小二乘平差，由此而得到的一组估值是满足方程的唯一解。如果方差阵D和权阵P是非对角阵，则表示观测值是相关的，按此准则进行的平差即称为相关观测平差。如果是对角阵，则表示观测值是彼此不相关的，此时称为独立观测平差。当观测值不相关，即P为对角阵时，则有①二VtPV二工PV2二PV2+PV2+...+PV2二minii1122nni=1当观测值不相关,并为等精度,即P=I时,则有:O=VtV=为v2=V2+V2+...+v2=mini12ni=1小结：最小二乘原理是测量平差的基本原理，按最小二乘准则求得的观测量及其函数的结果是最可靠的结果，后续所讲所有平差方法均按此准则求解。第五章条件平差第一节条件平差原理一、条件方程和改正数条件方程列出用观测值真值和真误差表示的条件平差函数模型导出用按最小二乘准则求得的观测值平差值和观测值改正数表示的条件平差的函数模型F'L）=0――条件方程r,1r,1AV-W=0――改正数条件方程r,nn,1r,1W=-F（L）――改正数条件方程常数项（闭合差）计算式举例（单三角形函数模型的建立）二、条件方程的纯量表达式和矩阵表达式r个条件方程的纯量表达式：

线性化后得改正数条件方程其中A=212n>av+avH线性化后得改正数条件方程其中A=212n>av+avHFav-w=01122nnabv+bvHFbv-w=01122nnbrv+rvHFrv-w=01122wawb=-F（L,L，…,L）12n=-F（L,L，…,L丿12nrr12naa•…awv12ni1bb…bwvi2n,W=2,V=2............•r,1n,1rr….rwv12nrn=-F（L,L，…,L）则改正数条件方程及其闭合差计算的矩阵表达式分别为AV-W=0W=-F（L）三、基础方程kT，称为联系数向量。组成函数r按求函数极值的拉格朗日乘数法，设其乘数为K=(kT，称为联系数向量。组成函数rabr,1①=VtPV-2Kt(av-W),对其求导整理得改正数V的计算公式V=P-1ATK=QATK——改正数方程当P为对角阵时，改正数方程的纯量形式为v=丄(ak+bkhfrk)i=1,2，…,nipiiaibir改正数条件方程与改正数方程联立，称为条件平差的基础方程。此时，方程的个数与未知数的个数相同，方程有唯一解。四、基础方程的解将改正数方程代入改正数条件方程，得AQAtK-W二0,令N=AQAt=AP-iAt,aa得NK-W二0――联系数法方程aa秩R（N）=R（AQAtLR（A）=r，即N是个r阶的满秩方阵，由此解出aaaaK=N-iWaa当P为对角阵时，法方程的纯量形式为raaak+raPab~k+pblk+•••+*br~k—W=0papbprb口…「1…•…nar~k+*bl~k+•••+rr~k—W=0papbprr-解出K，将其代入改正数方程，求出改正数V,在按L=L+V可求得平差值。五、条件平差步骤及示例用具有两个条件的符合水准网为例讲解。小结：本节应熟记条件方程，改正数条件方程，改正数条件方程闭合差计算式，法方程，改正数方程的表达形式，掌握用条件平差法平差的方法、步骤。第二节条件方程一、水准网（同§5.1中所述，略）二、测角网单三角形（同§5.1中所述，略）中心多边形以中心三边形为例，画出示意图，列出其条件方程和改正数条件方程的一般表达式。重点讲解极条件的列立方法和规律。举例（中心三边形实例）列条件方程和改正数条件方程。大地四边形画出示意图，列出其条件方程和改正数条件方程的表达式。重点讲解极条件的列立方法和规律。举例上图中，若以对角线交点为极列极条件，其极条件闭合差超限，说明角度观测存在问题，如何返工？先让让学生回答，然后教师讲解。三、测边网1.中心多边形画出测边中心三边形示意图。（1）列出以反算角表示的条件方程和改正数条件方程（2）建立反算角改正数与边改正数之间的关系（3）导出以边改正数表示的条件方程大地四边形画出测边大地四边形示意图。（1）列出以反算角表示的条件方程和改正数条件方程（2）建立反算角改正数与边改正数之间的关系（3）导出以边改正数表示的条件方程四、边角网如图，t=2p-q-3=8-3-3=2,r=n-1=8-2=6应列出6个条件方程条件分析：内角和条件2个正弦条件2个固定角条件1个规定边条件1个边角网条件方程列立例题讲解分析。小结：条件方程列立，首先应能正确确定应列的条件数目，保证方程之间不相关，其次应能分析条件类型，最后应掌握各类方程的列立规律，正确列出条件方程。第三节精度评定一、单位权方差估值计算VTPVO2=0rVTPV的计算：1、VtPV=PV2+PV2+…十PV2(权阵为对角阵时)1122nn2、VtPV=(QAtK)TP(QAtK)二KtAQPQAtK=KtNkaa3、VTPV=VTP(QATK)=VTATK=WTK二、协因数阵设r」ZT=LTWTKTVTLT列出各分块向量解的表达式及其微分式，利用协因数传播律导出各量的协因数阵和各量之间的互协因数阵的结果列于相应表中表中L与V、W、K的互协因数阵为零，说明L与V、W、K统计不相关证明：表中Q、Q的计算表达式。VVLL三、观测值平差值的精度评定D=£2QLL0LL四、平差值函数的精度评定1．平差值函数表达式及其协因数计算列出平差值函数表达式G=f(L,L，…,L)12n按泰勒公式展开，并按协因数传播律导出平差值函数协因数的计算公式Q=fTQf-(aqAn-iAQTaa珥(i=l,2,…,n)为偏导数值。

八八dL+ff1八八dL+•…+ff1八1%J1L=L血212L=L血丿nL=LdLnd(p=22．权函数式=f1dL1+f2dL2+…+fndLn—权函数式3．平差值函数的方差D=£2Q(p(p0(p(p小结:本节主要介绍了利用改正数计算单位权中误差的公式，各种平差量协因数和互协因数及方差协方差的计算，平差值函数式和权函数式的列立方法，平差值函数协因数和互协因数及方差协方差的计算方法，应重点掌握。第六章附有参数的条件平差一、概述设X=ZABD，又可列出1个极条件和一个固定边条件极条件为(以A点为极):固定边条件为(由AC一、概述设X=ZABD，又可列出1个极条件和一个固定边条件极条件为(以A点为极):固定边条件为(由AC边推算到AB边)：S=SABsinLsinG+L)^~~6〜8ACsinLsinX3sinLsinG+L),———6〜8=1S—ACSsinLsinXAB3由于选了一个参数，增加了一个条件，一般情况下，若选了U个参数，则条件方程的数目为c=r+u.从以上5个方程出发进行平差，就是附有参数的条件平差方法。二、基础方程观测量L和X的最佳估值L=L+V，X=X0+x，用奇表示的附有参数的条件平差函数模型为条件方程条件方程AV+Bx一W=0改正数条件方程c,nn,1c,uu,1c,1kbk卞，称为联系数向量。组成函数rW=_FC,Xkbk卞，称为联系数向量。组成函数r按求函数极值的拉格朗日乘数法，设其乘数为K二r,1①二VtPV-2Kt(Av+Bx-W),将①对V和x分别求一阶导数，并令其为零，导出改正数V的计算公式V二P-1AtK=QAtK—改正数方程附有参数的条件平差的基础方程为：AV+BX-W=0c,nn,1c,uu,1c,1V二P-1AtK二QAtKbtK=0方程的个数与未知数的个数相同，方程有唯一解。三、基础方程的解将改正数方程代入改正数条件方程，并令N=AQAt=Ap-1AT，则得法方程aaNK+BX-W二0aabtK=0法方程秩R(n)=R(AQAt)=R(A)=c，即N是个c阶的满秩方阵，顾及N=btN-1B由法方程aaaabbaa，可解出X二N-iBtN-iW,bbaaV二QAtN-1(W-Bx),aa四、精度评定(一)、单位权方差估值计算VTPVVTPVO2==0rc-uVTPV的计算：1、VtPV=PV2+PV2+…十PV2（权阵为对角阵时）1122nn2、VtPV二（QAtK）TPV二KtAV二KtW-KtBx=WtK3、VtPV=WtK=WtN-iW—WtN-1Bx=WtN-iW—GtN-iW）无aaaaaaaa（二）、协因数阵设r」Zt二*LtWtXtKtVtLt'列出各分块向量解的表达式及其微分式，利用协因数传播律导出各量的协因数阵和各量之间的互协因数阵的结果列于相应表中证明：表中Q农、Qvv的计算表达式。三）、观测值平差值的精度评定D=£2QLL0LLD=£2QXX0XX四）、平差值函数的精度评定对其全微分，得权函数式：式中d©=空dL+理dX=FTdL+FTdXdLdX按协因数传播律得6的协因数为:dL16L2帥6Ln帥dxnQ66=FT

FxT存XL

Q…］「FLXQ…FXX=FtQF+FtQF+FtQF+FtQFLXxxXLxXXx■/入入■/入入■/入入■/入入LXxxXLxXXxLL6的中误差为:0(p(pp0(p(p小结:掌握此种平差方法的应用范围，平差的方法步骤。第七章间接平差第一节间接平差原理、平差值方程与误差方程观测量L和X的最佳故值L=l+v,x=x0+x，用平差值和改正数表示间接平差的函数模型L=Fx—平差值方程（观测方程）n,1r,1V=BX-l—误差方程n,1n,tt,1n,1/=L-FG0L误差方程常数项（闭合差）计算式以测角单三角形为例，列出平差值方程和误差方程。二、方程的纯量表达式与矩阵表达式设有n个条件方程:TOC\o"1-5"\h\zL=FW,X，…,X八11(八12nL=FX,X,…,X2212nL=FX,X,…,Xnn12n线性化后得误差方程为v=aX+bXHFtx-111J1入21t1v=aX+bXFFtx-1221222t2>其中v=aX+bXFFtx-1其中nn1n2nt1=L-f\Xo,Xo,…,Xo11t12t1L-FXo,Xo,…,Xo22212t1L-F（Xo,Xo,,Xonnn12tB=a1a2bb2…t"1…t2,1=T112,V=v1v2ab…t1vnnnnn则误差方程的矩阵表达式为V二Bx-1误差方程常数项(闭合差)计算式的矩阵表达式为三、基础方程误差方程中未知数个数(n+1)大于方程个数n,方程有无穷多组解。根据最小二程原理可求得满足方程的唯一一组解。求VTPV的自由极值得基础方程JV=Bx-1|BTPV=0四、基础方程的解将基础方程第一式代入第二式，令N=BtPB,W=BtP1，得法方程bbNx-W=0法方程bb解上方程得：X=N-iWbb当P为对角阵时，法方程的纯量形式为[paa]x+tpablxFF[patk=[pa11[pabk+[pbbkFF[pbt]x=[pb1]12[patk+[pbtkHF[pttk=[ptl]12五、按间接平差法求平差值的计算步骤及示例用水准网例题讲解平差的方法步骤。小结：本节应熟记观测方程，误差方程，误差方程常数项计算式，法方程的表达形式，掌握用间接平差法平差的方法、步骤。第二节误差方程一、参数个数的确定与选取参数个数：等于必要观测数t；参数选取：水准网一般选择未知点高程为参数，也可选择观测高差为参数；平面控制网一般选择未知点坐标为

参数，也可选择观测角度等为参数。参数选择要求：足数；参数间线性无关。二、平差值方程及误差方程的列立1、观测高差平差值方程及误差方程的列立例1，以具有两个未知点的符合水准网为例讲解2、观测方向平差值方程及其误差方程的列立设X=0YXY大jjKK计算参数近似值Xo=Sx0Y0X0Y0jjkk平差值方程：入入入入Y-YL=—Z+ct=—Z+arctgk」jkj伙jX—Xkj—ztj—ztj+c0+8jkJKnL+V=—Z0—z+c0+8jkjkjjjkcJKnV=—zt+8—ljkjcJKjk其中—l=—Z0+c0—L——常数项计算式jkjjkjkY0—Y0c0=arctgkj—jkX0—X0kj〃QO=.QXjQO〃QO=.QXjQOQO入QOxH.yX+—jjQYjQXkQYX=xojX=xokX=XokX=X00p"Ao=~YtkVXSo2.JKft■八psinOo入=—XSo.JKp"Ao-ir^yjJKH八pCOSOLJKo.ySojJKo2kJKtt•tfpsinooo入p〃cosooS0JKo.ySokJK则观测方向的误差方程为：〃P〃Av=-z.+~^~

〃P〃Av=-z.+~^~

jkjSJKJKop"A「―Y卜JKo2JK0八xkp"Ao-1S02kjkJK〃v=jk八八”八-z+ax+by..k..k.ak匕「称」、k方向的方向系数，对于任一方向jm有：jk、jk〃八八『八v=-z+ax+by.m..m..m.坐标近似值的计算：-aXjkk-aXjmkjkyk-1jk+by-1jmkjm可用支导线法、前方交会法等方法计算。,nji,nji=1,2,Zo=oo—L,jiji/ji\工uo-L丿Zo=.n.误差方程列立规律：①符号；②系数；③特殊情况；④5"=5"TJKTKJ单位：坐标改正数为厘米时系数除100，...。3、观测角度平差值方程及其误差方程的列立平差值方程：误差方程：L=0L-OLijhjkJK〃=50Jh〃〃-5JK-50JK+O0-O0jhjkl=00-00-Lijhjki例2，以固定角内插一点得测角网为例讲解方程列立及求平差值的方法、步骤4、观测边长平差值方程及其误差方程的列立设：X=0YXjjk设：X=0YXjjkk平差值方程：sisi=S0+5isi其中6S5—i_入oS入qS入oSx+—壬y+—“x+—十j0丫jOXkoyX=X0jX=X0kX=X0kX=X05rsidXjA0A0A0A0误差方程：XJKS0jk-x—jYJKS0jky+jXJKS0jk-x+—ykS0kjkA0A0A0A0V—一XxJK一Y八y+X八X+Y八y-lJKJKJKSiS0jS0jS0kS0kijkjkjkjk=-cosT0x一sinT0y一cosT0x+sinT0y一ljkjjkjjkkjkki常数项：例3，小结：o一yo例3，小结：o一yo)(单位ixkj"kj以中心三边形内差一点的测边网为例讲解求未知点坐标的方法、l=S—S0=S—X0—X0iiiikj观测方程和误差方程的列立，首先应能正确确定应选参数数目，cm或dm)步骤。保证所选参数之间线形无关，其次应能掌握各类方程的列立规律，正确列出相应观测方程和误差方程。第三节精度评定、单位权方差估值计算VTPVO2=0rVTPV的计算：1.VtPV=PV2+PV2+…+PV2(权阵为对角阵时)1122nn2.VtPV=(BXc—l\PV=xtBtPV—ItPV=—ItP(BX—l)&及BtPV=0)=ItPI—ItPBx=ItPI—GtPlh3、在线性方程组解算表中计算二、协因数阵与互协因数阵设：(L(L1'E'f0、7KxN-1BtPN-1BtP=bbL—bbVbn-1BtP—EBN-1BtP—E八bbbbILJ.BN-1BtP,.BN-1BtP—E丿F(X0)按协因数传播导出各量的协因数阵和各量之间的互协因数阵的结果列于各量的协因数阵和各量之间的互协因数阵表中L与v和x与v的互协因数阵为零，说明L与v、x与v统计不相关证明表中Qxx，QVV，Qtt的计算表达式。三、参数的精度评定设所求量（如未知点高程或纵横坐标）为参数Xj,i=l,2,…，t,则&2=&2QX0XXii四、参数函数的精度计算设参数函数为：e设参数函数为：e=e（x）=e（X，1X，…x)2线性化得权函数式为：dgOX1dX+——kdgOX1dX+——k—1OX2dXO©+•…+——2OXt八X=x0=fdX+fdX+…+fdX=FTdX1122ttdXt由协因数传播律得：QJFtQXF五、各种平差量权函数式的列立1、高差平差值如图设未知点高程为参数，所求高差平差值的函数式为($=h=-X+Xijk其权函数式为d(9=dh=-dX+dXijk若j、k为已知点，其dX前的系数为零。2、方位平差值如图设未知点坐标为参数，所求方位平差值的函数式为a=

jkY-Yarctgja=

jkX-Xkj求全微分得其权函数式为宀P"AYo「

da—f—、jkdXjkfSob2.1oJjK式中加的单位为（"），3、角度平差值dX、p'AXo入p'AYo入p'AXo-f-、jkdY-f、jkdX+f-Vjk\Sob.10j\Sob.10k\sob.1jKjKjKdY的单位为分米，若j、k为已知点，其dX、"AXodYKdY前的系数为零如图设未知点坐标为参数，所求角度平差值的函数式为jX-X

kj入(p'AY0p'AY0|入da=|jX-X

kj入(p'AY0p'AY0|入da=|fj—f__、jk|dXjk丨\S0丄•10S0艾•10丿jjk""AX04dY•10Kjhp"AY0+^TS02•10jkp""AX0、jhS02•10jhp""AY0—f——\jhdX

S02•10jh\S0尢-10jkp'AXo+S0龙-10K£0龙-10h£0北•10jkjhjhdX、dY的单位为分米，若j、k为已知点，其dX、f/\4^dX—k4dYh式中叫的单位为（"）,4、边长平差值dY前的系数为零。如图设未知点坐标为参数，所求边长平差值的函数式为/\JI/\/\II/\/\IS=\：X—Xb+Y—Yk

jkkjkj求全微分得其权函数式为dSjkAX0AY0AX0AY0--pjkxdX-KjXdY+^-AdX+P^dY\S0jS0/S0kS0kjkjkdY的单位为分米,jkjh若j、k为已知点，其dX、dY前的系数为零。第八章附有限制条件的间接平差八入入Y—YY—YL=a—a二arctg4—arctg訣jhjkX—Xh求全微分得其权函数式为、概述如上图，选取i、k两点的坐标为未知数，可列出4个平差值方程。由于选定的未知数个数（u）多于必要观测数（t）,所以在所选定的未知数之间存在s=u-t个限制条件。—XKi—XKi—S=0KiikY-Y门arctg-a=0X-XjKi把上列两式线性化得cosaox+sina0y-cosa0x-sma0y-w=0jkKjkKjkKjkK1ax+by-ax-by-w=0jkKjkKjkKjkK(2w=S-‘0-X0)+V0—Y01ikkjkjY0-Y0w=a-arctgk—i—2jkX0-X0ki二、基础方程已知附有参数的条件平差法的函数模型s,1s,1其线性形式为A=Bx-1c~x-W=0x其中W=—①Qo)x由于n+svn+u，不能求得A和~的唯一解，只能按最小二乘原理求A和~的最佳故值v和x，从而求得观测量L和文的最佳故值L和攵，即为此，可用观测值平差值和参数平差值表示附有参数的条件平差的函数模型，即L=fQ>—平差值方程（观测方程）n，（1©乂厶0—限制条件方程s，1或用观测值改正数和参数改正数表示附有限制条件的间接平差法的函数模型，即V=Bx—l—误差方程n,1n,tt,1n,1CX-W=0一限制条件方程Xl二L-FG0L误差方程常数项（闭合差）计算式W=6^Xo）一限制条件方程常数项（闭合差）计算式X按求函数极值的拉格朗日乘数法’设其乘数为K叫J小，称为联系数向量。组成函数6=VtPV+2Kt（CX-W），sX将①对X求一阶导数，并令其为零，得36dV二=2VtP—+2KtC=2VtPB+2KtC=0，3x3xss转置得BtPV+CtK=0，s上式与误差方程和限制条件方程联立得附有参数的条件平差的基础方程：BtPV+CtK二0sV=BX-ln,1n,tt,1n,1CX-W二0x方程的个数与未知数的个数相同，方程有唯一解。三、基础方程的解将基础方程的第二式代入第一式与第三式联立，得NX+CtK-W二0，bbsCX-W二0—附有限制条件的间接平差法的法方程X将法方程第一式左乘CN-1与第二式相减，得bbCN-iCtK—Qn-iW-WL0bbsbbX令N二CN-iCtccbb则有NK—Qn-iW-WL0ccsbbX

式中N的秩R(N)-R(CN-1Ct)-R(C)-S,且Nt=Qn-1Ct)=CN-C，故N为sccccbbccbbbbcc阶满秩对称方阵。K=N-1Qn-iW—W)SccbbX将上式代入法方程第一式，可解得X=w-1—N-iCtN-1CN-1W+N-iCtN-iW,bbbbccbbbbccX代入误差方程可解出改正数V,从而可解出:四、精度评定一）、单位权方差估值计算VTPVVTPV02=0rn—u+sVTPV的计算:1.V1.VtPV=PV2+PV2+1122+PV2（权阵为对角阵时）nn2.VtPV=(Bxc-1》PV二XtBtPV-ItPV=—XtCtK—ItP(BX—l)G及BtPV=—CtK)SS=ITPl—XtCtK—ItPBXS=ItPI—WtK—WtXXS二）、协因数阵与互协因数阵3、在线性方程组解算表中计算令：Zt=VtWtKtXtVtLt)二）、协因数阵与互协因数阵S列出各分块向量解的表达式及其微分式,利用协因数传播律导出各量的协因数阵和各量之间的互协因数阵的结果列于相应表中,讲解。（三）、参数的精度评定设所求量（如未知点高程或纵横坐标）为参数Xi,i=l,2,…,t,则cy2=cy2QX0XXiii四）、参数函数的精度计算设参数函数为：（x）=o/\/X/X(X,X，…x)12t线性化得权函数式为：d(p=dt-dXi=fdX+fdX+•••+fdX1122tt由协因数传播律得：dX+-51dX2dXd©+…+^-2dXt=FTdXQrFTQXFD=02Qp0<Po=cy(Po小结：掌握此种平差方法的应用范围，平差的方法步骤。一、点位误差的概念及计算dXt第九章误差椭圆§9.1概述；§9.2点位误差1、点位真误差如图可得：A2=A2+A2,PxyA2=A2+A2,PSu无法求得（为什么?）2、点位方差及其计算由方差的定义式可得：Lx_E(x)LeLx-~)2LeL~—x)2]=Ey~y)2=eL~—y)2Eo2=E(xo2=E(y－y(A2)xA2)

y故有E(A2)=E(A2)+E(A2)=Q2+q2Pxyxy同理有：E(A2)=E(A2)+E(A2)=02+02PsuSu记02=E(A2)，则有：pP02=02+02=02+02――点位方差计算式Pxysu上式说明点位方差02的大小与坐标轴的方向无关,即与坐标系的选择无关P用点位方差衡量P点精度的缺陷：不能完善说明P点在任一各方向上的精度情况，不能确定P点在哪一个方向上的精度最好(最差)。二、二、由图可得下列关系式:Ap=TP=PP*+PP=cos申Ax+sin^Ay=cos=cospsinp_AxAy由协方差传播律得:0p2=02cos2p+02sin2p+02sin2ppxyxy02p=00Qp=02Qcos2p+Qsin2p+Qsin

0xyxy上式即为求任意方位角少方向上点位方差的计算公式。

三、位差的极值方向、极大值和极小值的确定由位差计算式可以看出，Qp随着少值的变化而改变，其具有最大值和最小值。函数有极值，其一阶导数等于零，设位差的极值方向为申，求导得出02Qxy—Q-Qxy将®0代入位差计算式得:o2=o2(Qcos2®+Qsin2®+Qsin2®)2tg®o—1+tg2®0TOC\o"1-5"\h\z®02tg®o—1+tg2®0=o2(Qcos2®+Qsin2®+Q0x0y0xy极值方向的判别方法:Q〉0,极大值在第I、III象限，极小值方向在第II、W象限；QV0,xyxy极大值在第II、W象限，极小值方向在第I、III象限位差极大值、极小值的计算：用®表示极大值方向、®=®+900表示极小值方向；用E、F分别表示位差的极大值和极小值。EFe则有E2=o2(Qcos2®+Qsin2®+Qsin2®)TOC\o"1-5"\h\z0xEyExyE>F2=o2(Qcos2®+Qsin2®+Qsin2®)0xFyFxyF把®代入位差计算式整理得0E2=1o2{q+Q)+K120xyIF2=1o2{Q+Q)-K}20xy‘其中K=^Q-Q)+4Q2xyxyo2与E2、F2有下面关系:Po2=E2+F2

P四、用E、F表示的任意方向屮上的位差

由图可知,任意方向在两个坐标系中的方位角有如下关系:把申=屮+9代入位差计算式整理得:EQ2=E2由图可知,任意方向在两个坐标系中的方位角有如下关系:把申=屮+9代入位差计算式整理得:EQ2=E2cos2屮+F2sin2屮例1如图，在固定三角形内插入一点P,经过平差后求得P点坐标的协因数阵为:(Q(3.81xxy[Q.Q丿yxy单位权方差估值为&2=1.96(2)。0+0.36、1+0.362.93丿试求(1)位差的极值方向9£和®f，⑵位差的极大值E与极小值F,(3)已算出PM的方位角aPM=75。29'，PM方向上的点位误差为多少,(4)P点的点位方差。例2如图,已知AT=345。18'AB为确定P点的位置,作如下观测0-8935'42〃土4〃.O,S=600.150m土10mm试确定P点位差的极大值及其方向。例3在例2中,平差后算得PA的方位角和边长^=248。15'30''訂；1827沁试求PA边的方位误差n受边长相对中误差。aAP按解算公式和相应方法解算，讲解。小结：点位方差°2的大小与坐标系的选择无关，位差可描述P点