误差理论与测量平差基础(全套课件299P).ppt

1、南京林业大学土木工程学院误差理论与测量平差基础误差理论与测量平差基础测绘工程专业核心课程Errors Theory and Foundation of Surveying Adjustment主讲教师：隋铭明主讲教师：隋铭明南京林业大学土木工程学院课程结构课程结构绪论绪论误差分布与精度指标误差分布与精度指标协方差传播律及权协方差传播律及权最小二乘原理最小二乘原理条件平差条件平差附有参数的条件平差附有参数的条件平差间接平差间接平差附有限制条件的间接平差附有限制条件的间接平差概括平差函数模型概括平差函数模型误差椭圆误差椭圆2南京林业大学土木工程学院Ch1 绪论绪论v课程基本情况 v教材 误差理论与

2、测量平差基础 误差理论与测量平差基础习题集 武汉大学出版社总学时总学时讲课讲课上机上机习题习题725012103南京林业大学土木工程学院Ch1 绪论绪论v怎样学好测量平差怎样学好测量平差预习、复习加习题练习预习、复习加习题练习独立思考并推导公式独立思考并推导公式平差思想和解题思路平差思想和解题思路高数高数 线代线代 概率概率习题练习公式推导数学基础习题练习公式推导平差思想平差思想数学基础4南京林业大学土木工程学院Ch1 绪论绪论v为什么要学测量平差？为什么要学测量平差？ 1. 测量过程中可能会出现 照错目标 读错数 如何避免错误或及时发现错误？ 解决方法：增加多余观测。 2. 有多余观测，如何

3、消除不符，求出最优值？ 5南京林业大学土木工程学院Ch1 绪论绪论测量平差的任务和意义测量平差的任务和意义q 任务 1）消除不符值，寻求未知参数的最佳估值； 2）评定结果的精度。q 意义 所有观测数据只有通过平差才能使用，即测量平差是测绘科学和技术的基础和灵魂。6南京林业大学土木工程学院Ch1 绪论绪论v测量平差的作用和地位测量平差的作用和地位1）解决测量工作中的实际问题，对测量数据进行处理，求出最佳估值。2）是测绘学科的基础理论，是对仪器操作和基本测量方法的主要补充。3）其核心知识是后续专业课程的重要基础，如大地测量、GPS测量原理、变形监测等。4）是测绘工程专业研究生入学考试课程，是硕士和

4、博士阶段的重要课程。7南京林业大学土木工程学院Ch1 绪论绪论v课程结构参见目录8章节章节主要内容主要内容Ch1绪论绪论Ch2- Ch3平差基础知识平差基础知识Ch4 平差基本原则平差基本原则Ch5- Ch8Ch5- Ch8四种经典平差方法四种经典平差方法Ch9平差方法总结平差方法总结Ch10Ch10点位精度讨论点位精度讨论Ch11统计假设检验统计假设检验Ch12Ch12近代平差简介近代平差简介南京林业大学土木工程学院Ch1 绪论绪论v基本概念 误差对未知量进行测量的过程称为观测，测量所得的结果称为观测值。观测值与其真实值（真值）之间的差异称为测量误差或观测误差，通常称真误差，简称误差。 测量

5、平差测量平差是测量数据调整的意思。其定义是，依据某种最优化准则，由一系列带有观测误差的测量数据，求定未知量的最佳估值及精度的理论和方法。9南京林业大学土木工程学院1.1 观测误差观测误差 一、误差来源测量仪器：仪器精密度；仪器轴线关系引起。观测者：操作水平，工作态度，使用习惯。外界环境：温度，湿度，风力，大气折光等。10南京林业大学土木工程学院1.1 观测误差观测误差二、误差分类 偶然误差在相同误差在大小和符号上表现出偶然性 系统误差误差在大小和符号上表现出系统性，或按一定规律变化 粗差即错误11南京林业大学土木工程学院1.1 观测误差观测误差误差名称误差名称误差特点误差特点消除或削弱的办法消

6、除或削弱的办法举例举例偶然误差Random error单个误差没有规律性，整体具有统计规律，服从或近似服从正态分布采用测量平差的方法照准误差对中误差估读误差系统误差Systematic error误差在大小和符号上表现出系统性，或按一定规律变化，或为常数采用适当的观测方法校正仪器计算加改正尺长误差i角误差粗差Gross error即大的偏差或错误重复观测严格检核发现舍弃或重测大数读错输入错误照错目标12南京林业大学土木工程学院1.2 测量平差的研究对象测量平差的研究对象研究对象：带有误差的观测值经典测量平差：只含有偶然误差的观测值近代测量平差：观测值除了含有偶然误差，还含有系统误差或粗差，或两

7、种兼有。平差问题的解决思路：13分析观测值选择平差准则确定平差模型解算模型精度评定南京林业大学土木工程学院1.3 测量平差简史及发展测量平差简史及发展v1794年，C.F. Gauss从概率统计角度提出了最小二乘法v1806年，A.M. Legendre从代数角度提出了最小二乘法v1809年，Gauss在天体运动的理论一文中发表，称为Gauss- Legendre方法v1912年，A.A. Markov，对最小二乘原理进行了证明，形成数学模型（函数模型+随机模型）v近代发展v现在的国内相关专家1412020)(, 0lim)(PQAXLEnEAXLnPLAPAAXTT1)(南京林业大学土木工程

8、学院1.4 本课程的任务和内容本课程的任务和内容v本书主要为经典测量平差内容，即只讨论带有偶然误差的观测值。（1）偶然误差理论。偶然误差特性，传播；精度指标及估计；权。（2）测量平差的函数模型和随机模型，最小二乘原理。（3）测量平差的基础方法。条件平差，附有未知参数的条件平差，间接平差，附有限制条件的间接平差。平差计算模型及精度评定公式，各种平差方法的概括及联系。（4）测量平差中的统计假设检验方法。15南京林业大学土木工程学院16南京林业大学土木工程学院Ch2 误差分布与精度指标误差分布与精度指标偶然误差的规律性偶然误差的规律性1正态分布正态分布2精度及其衡量精度指标精度及其衡量精度指标3本章

9、总结及习题本章总结及习题417南京林业大学土木工程学院2.1偶然误差的规律性偶然误差的规律性基本假设：基本假设：系统误差已消除，粗差不存在，即观测系统误差已消除，粗差不存在，即观测误差仅为随机误差。误差仅为随机误差。偶然误差：偶然误差：单个误差在误差大小及符号上没有明显单个误差在误差大小及符号上没有明显的规律，表现出随机性，称为偶然误差。但对大的规律，表现出随机性，称为偶然误差。但对大量误差进行统计具有明显的规律。量误差进行统计具有明显的规律。寻找偶然误差之规律性的方法寻找偶然误差之规律性的方法(统计分析统计分析)： 1、统计表、统计表 2、直方图、直方图 3、误差分布、误差分布iiiLL 1

10、8南京林业大学土木工程学院统计表统计表误差区间+ +个数K频率K/n(K/n)/d个数K频率K/n(K/n)/d0.000.20450.1260.630460.1280.6400.200.40400.1120.560410.1150.5750.400.60330.0920.460330.0920.4600.600.80230.0640.320210.0590.2950.801.00170.0470.235160.0450.2251.001.20130.0360.180130.0360.1801.201.4060.0170.08550.0140.0701.401.6040.0110.05520.

11、0060.0301.60000000和1811810.5050.5051771770.4950.495o例1：在相同的条件下独立观测了358个三角形的全部内角，三角形内角和应为180度，但由于误差的影响往往不等于180度，计算各内角和的真误差，并按误差区间的间隔0.2秒进行统计。19南京林业大学土木工程学院(K/n)/d00.40.60.8-0.8-0.6-0.4闭合差概率密度函数曲线面积= (K/n)/d* d= K/n所有面积之和=k1/n+k2/n+.=1直方图直方图202.1偶然误差的规律性偶然误差的规律性南京林业大学土木工程学院偶然误差的特性偶然误差的特性由统计分析可以看出，偶然误差

12、具有下列特性：1、有界性：在一定的观测条件下，偶然误差的绝对值有一定的限值，即超过一定限值的偶然误差出现的概率为零2、聚中性：绝对值较小的偶然误差比绝对值较大的偶然误差出现的概率大；3、对称性：绝对值相等的正负偶然误差出现的概率相同；4、抵偿性：偶然误差的理论平均值为零，即01lim1niiin212.1偶然误差的规律性偶然误差的规律性南京林业大学土木工程学院l例2：在相同的条件下独立观测了421个三角形的全部内角，每个三角形内角之和应等于180度，但由于误差的影响往往不等于180度，计算各内角和的真误差，并按误差区间的间隔0.2秒进行统计。误差区间+个数K频率K/n(K/n)/d个数K频率K

13、/n(K/n)/d0.000.20400.0950.475460.0880.4400.200.40340.0810.405410.0850.4250.400.60310.0740.370330.0690.3450.600.80250.0590.295210.0640.3200.801.00200.0480.240160.0430.2151.001.20160.0380.190130.0400.200.2.402.6010.0020.01020.0050.00252.60000000和2102100.4990.4992112110.5010.50122南京林业大学土木工程学院 频数/d00.40

14、.6 0.8-0.8-0.6-0.4闭合差0.630 频数/d00.40.6 0.8-0.8-0.6-0.4闭合差0.475 频数/d00.40.6 0.8-0.8-0.6-0.4闭合差 00.40.6 0.8-0.8-0.6-0.4闭合差提示：观测值定了其分布也就确定了，因此一组观测值对应相同的分布。不同的观测序列，分布不同。但其极限分布均是正态分布。22221)(ef图1图223南京林业大学土木工程学院n 频数/d00.40.6 0.8-0.8-0.6-0.4闭合差 00.40.6 0.8-0.8-0.6-0.4闭合差22221)(ef 当偶然误差的个数 时，偶然误差出现的频率就趋于稳定。

15、此时，若把偶然误差区间的间隔无限缩小，则直方图将分别变为如图所示的两条光滑的曲线。242.2正态分布正态分布南京林业大学土木工程学院 由概率论知，该曲线是正态分布正态分布的概率分布曲线。高斯在研究误差理论时最先使用了这一分布，所以正态分布又称为高斯分布高斯分布。测量上通常将正态分布作为偶然误差的理论分布。或者说偶然误差服从正态分布。其密度函数为：式中： 和 为参数。,)(21exp21)(22f2.2正态分布正态分布25南京林业大学土木工程学院由密度函数知，偶然误差 为正态随机变量。所以又称偶然误差为随机误差。 下面来看参数 和 是什么。 对正态随机变量 求数学期望： ,)(21exp21)(

16、22fddfE22)(21exp21)()(2.2正态分布正态分布26南京林业大学土木工程学院作变量代换，令 得因tdttdtttdtttE22221exp221exp221exp)(21)(221exp,021exp22dttdttt2.2正态分布正态分布27南京林业大学土木工程学院2.2正态分布正态分布所以再求 的方差 。同样作变量代换，可得：22)(E)(DddfED2222)(21exp)(21)()()(2222)(D28南京林业大学土木工程学院 由以上推导知，参数 和 分别是随机误差 的数学期望和方差。它们确定了正态分布曲线的形状。 由 知，随机误差 的数学期望等于零。 由正态分布

17、知，正态分布曲线具有两个拐点，这两个拐点在横轴上的坐标为 方差的几何意义是：方差是正态分布曲线的拐点横坐标。01lim)(1niinnE拐292.2正态分布正态分布南京林业大学土木工程学院2.3精度及其衡量精度指标精度及其衡量精度指标 观测值的质量取决于观测误差（偶然误差、系统误差、粗差）的大小。 1、精度：、精度：指误差分布的密集或离散程度，可利用方差指误差分布的密集或离散程度，可利用方差协方差阵描述。协方差阵描述。 2、准确度：、准确度：描述系统误差和粗差，可用观测值的真值描述系统误差和粗差，可用观测值的真值与观测值的数学期望之差来描述，即：与观测值的数学期望之差来描述，即： 3、精确度：

18、、精确度：是精度和准确度的合成，描述偶然误差、是精度和准确度的合成，描述偶然误差、系统误差和粗差的集成，精确度可用观测值的均方误差系统误差和粗差的集成，精确度可用观测值的均方误差来描述，即：来描述，即： 当当 ，即观测值中不存在系统误差和粗差时，即观测值中不存在系统误差和粗差时，亦即观测值中只存在偶然误差时，均方误差就等于方差，亦即观测值中只存在偶然误差时，均方误差就等于方差，此时精确度就是精度。此时精确度就是精度。)(LEL 222)()()(LLELLELMSELLLE)(30南京林业大学土木工程学院精度、准确度和精确度的形象描述精度、准确度和精确度的形象描述2.3精度及其衡量精度指标精度

19、及其衡量精度指标31精度准确度精确度南京林业大学土木工程学院 4、衡量精度的指标、衡量精度的指标 精度虽然可以通过直方图或分布曲线的形状来描精度虽然可以通过直方图或分布曲线的形状来描述，但在实际工作中很麻烦，且不能用一个数字来衡述，但在实际工作中很麻烦，且不能用一个数字来衡量其高低。为此，人们希望通过一个数字来偶然误差量其高低。为此，人们希望通过一个数字来偶然误差的离散程度。能反映偶然误差的离散程度的数字称为的离散程度。能反映偶然误差的离散程度的数字称为衡量精度的指标。这样的数字很多，比如：衡量精度的指标。这样的数字很多，比如： 4.1、方差和中误差、方差和中误差 设在相同的观测条件下得到一设

20、在相同的观测条件下得到一组独立观测误差组独立观测误差 ，则其方差定义为：，则其方差定义为： ndfEDniin12222lim)()()(i2.3精度及其衡量精度指标精度及其衡量精度指标32南京林业大学土木工程学院2.3精度及其衡量精度指标精度及其衡量精度指标 方差的算术平方根定义为中误差，即 在实际工作中，n总是有限的，由有限个观测值的真误差只能求得方差和中误差的估值： 和 nniin12limnnii12nnii12233南京林业大学土木工程学院 4.2、平均误差 设在相同的观测条件下得到一组独立观测误差 ，则其平均误差由 之绝对的数学期望定义，即： 因为 所以 2.3精度及其衡量精度指标

21、精度及其衡量精度指标iindfEniin1lim)()(dfdf)(2)(0547979. 0245253. 1234南京林业大学土木工程学院 由上式知，不同的 ，对应着不同的 ，于是就对应着不同的误差分布曲线。所以平均误差 也可作为衡量精度的指标。 在实际工作中，既可通过以上等量关系来计算平均误差的估值： 也可由下式计算之： 2.3精度及其衡量精度指标精度及其衡量精度指标nnii1nnii12545435南京林业大学土木工程学院 4.3、或然误差 当观测误差出现在 之间的概率等于二分之一时，称 为或然误差（如图），即令 ，则有 由概率积分表可查得，当概率为二分之一时，积分限为0.6745，于

22、是可得中误差与或然误差的理论关系： 2.3精度及其衡量精度指标精度及其衡量精度指标),(21)(df1/41/41/20t212exp212)(2/0dttdf234826. 1,326745. 036南京林业大学土木工程学院2.3精度及其衡量精度指标精度及其衡量精度指标 中误差、平均误差和或然误差都可以作为衡量精度的中误差、平均误差和或然误差都可以作为衡量精度的指标，但由于指标，但由于当当n不大时，中误差比平均误差更能反映大误差的影响不大时，中误差比平均误差更能反映大误差的影响中误差具有明确的几何意义（分布曲线的拐点坐标）中误差具有明确的几何意义（分布曲线的拐点坐标）平均误差和或然误差都与中

23、误差存在理论关系平均误差和或然误差都与中误差存在理论关系 所以，世界上各国都采用中误差作为衡量精度的指所以，世界上各国都采用中误差作为衡量精度的指标，我国也统一采用标，我国也统一采用中误差中误差作为衡量精度的指标。作为衡量精度的指标。37南京林业大学土木工程学院2.3精度及其衡量精度指标精度及其衡量精度指标4.4、极限误差 由中误差的定义知，中误差是一组同精度观测误差的平方的平均值的平方根的极限。既然是平均值，就会有的观测误差的绝对值比中误差大，有的观测误差的绝对值比中误差小。那么，绝对值比中误差小的观测误差出现的概率是多少？绝对值比中误差大的观测误差出现的概率又是多少呢？由下图，通过积分38

24、南京林业大学土木工程学院2.3精度及其衡量精度指标精度及其衡量精度指标 可得观测误差 出现在给定区间 内的概率为：)(kkdkkPkk222exp21)(作变量代换，得k=1，2，3时的概率分别为：%7 .99)33(%5 .95)22(%3 .68)(PPP 上式表明：绝对值大于中误差的观测误差出现的概率为31.7%；绝对值大于二倍中误差的观测误差出现的概率为4.5%；绝对值大于三倍中误差的观测误差出现的概率仅为0.3%。即观测误差的绝对值一般不会大于三倍中误差。因此，实际工作中通常以三倍中误差作为观测误差的极限，并称为极限误差，用 表示。限39南京林业大学土木工程学院4.5、相对误差 观测

25、值的中误差与观测值本身之比，称为相对误差，常用表示 。N12.3精度及其衡量精度指标精度及其衡量精度指标40南京林业大学土木工程学院本章小结本章小结几个名词几个名词误差误差观测误差观测误差偶然误差偶然误差 随机误差随机误差精度精度精确度精确度系统误差系统误差准确度准确度粗差粗差真误差真误差方差中误差方差中误差平均误差平均误差或然误差或然误差极限误差极限误差相对误差相对误差绝对误差绝对误差衡量精度的指标衡量精度的指标1、几个基本概念及相互关系41南京林业大学土木工程学院2、一个事实、一个事实 不论观测条件如何，观测误差总是不可避免的。不论观测条件如何，观测误差总是不可避免的。3、基本假设、基本假

26、设 在本课程中，我们假设观测误差为偶然误差，即不在本课程中，我们假设观测误差为偶然误差，即不含系统误差和粗差。换句话说，我们假设观测误差为服含系统误差和粗差。换句话说，我们假设观测误差为服从正态分布的随机误差。从正态分布的随机误差。4、统计规律、统计规律 在一定的观测条件下，偶然误差的绝对值有一定的在一定的观测条件下，偶然误差的绝对值有一定的限值，即超过一定限值的偶然误差出现的概率为零；限值，即超过一定限值的偶然误差出现的概率为零；绝对值较小的偶然误差比绝对值较大的偶然误差出现的绝对值较小的偶然误差比绝对值较大的偶然误差出现的概率大；概率大；绝对值相等的正负偶然误差出现的概率相同；绝对值相等的

27、正负偶然误差出现的概率相同；偶然误差的理论平均值为零。偶然误差的理论平均值为零。本章小结本章小结42南京林业大学土木工程学院43南京林业大学土木工程学院Ch3协方差传播律及权协方差传播律及权44南京林业大学土木工程学院3.1观测向量及其方差观测向量及其方差协方差阵协方差阵 作为衡量精度的指标，中误差可衡量一组观测值的精度。在实际工作中，我们得到的观测值往往是由多个观测值所构成的观测向量。 因此，需要引入观测向量矩阵和方差协方差矩阵。一、协方差一、协方差 对于变量X、Y，其协方差为：)()()()(XEXYEYEYEYXEXEYXXY45YXXYnYXXY南京林业大学土木工程学院3.1观测向量及

28、其方差观测向量及其方差协方差阵协方差阵二、协方差阵二、协方差阵 设有n维观测向量为 则其方差协方差阵定义为：特点：特点：对称；正定；互不相关时为对角矩阵，对角线元素相等时，为等精度观测。46TnnlllL)(2112222122121211)()(nnnnnlllllllllllllllTnnLLLELLELED南京林业大学土木工程学院3.1观测向量及其方差观测向量及其方差协方差阵协方差阵三、互协方差阵三、互协方差阵 设有两组观测向量为， n维的X，r维的Y。则，它们的互协方差阵为：47rnnnrrYXYXYXYXYXYXYXYXYXXYD212221212111YXZ思考：思考：若求DZZ？

29、南京林业大学土木工程学院3.2协方差传播律协方差传播律1、协方差传播律的作用、协方差传播律的作用 计算观测向量函数的方差协方差矩阵，从而评定观测向量函数的精度。2、预备公式、预备公式当随机变量 两两独立时，有 )()()(, )()(,)(YEXEYXEXCECXECCE)()()()(2121nnXEXEXEXXXEnXXX,21)()()()(2121nnXEXEXEXXXE48南京林业大学土木工程学院3.2协方差传播律协方差传播律3、观测向量线性函数的方差、观测向量线性函数的方差 设观测向量X及其期望和方差为：观测向量线性函数为 式中： 为常数。TnTnXEXEXEXEXXXX)()()

30、()(,)(21210kKXZ021,kkkkKnTXXnnnnnTXXDXEXXEXED2212221211221)()(49南京林业大学土木工程学院3.2协方差传播律协方差传播律Z的期望为Z的方差为即展开成纯量形式：00)()()(kXKEkKXEZETXXTTTTZZKKDKEXEXXEXKEkXKEkKXkXKEkKXEZEZZEZED)()()()()()()(0000TXXZZKKDD50nnnnnnnnZZkkkkkkkkkD, 111112212222222121222南京林业大学土木工程学院v例题1v例题2v例题351南京林业大学土木工程学院3.2协方差传播律协方差传播律4、

31、多个观测向量线性函数的协方差阵、多个观测向量线性函数的协方差阵 若观测向量的多个线性函数为则令0221120222212121012121111tntntttnnnnkXkXkXkZkXkXkXkZkXkXkXkZ0201010212222111211211,tttnttnnntttkkkKkkkkkkkkkKZZZZ52南京林业大学土木工程学院3.2协方差传播律协方差传播律于是，观测向量的多个线性函数可写为 。故有 式中： 为对称方阵。 若还有观测向量的另外r个线性函数 其矩阵形式为：0KKXZTXXZZKKDD0221120222212121012121111rnrnrrrnnnnfXfX

32、fXfYfXfXfXfYfXfXfXfY0FFXYTZZZZDD53南京林业大学土木工程学院3.2协方差传播律协方差传播律则有：而同理：TYYTXXrrYYDFFDDTYZTXXrtZYDFKDDTXXTTTTtrYZKFDKEXEXXEXFEkXKEkKXFXFEFFXEZEZYEYED)()()()()()()(000054南京林业大学土木工程学院3.2协方差传播律协方差传播律5、多个观测向量非线性函数的协方差阵、多个观测向量非线性函数的协方差阵 基本思想：a、用全微分代替全增量，得到函数误差表达式（线性近似）；b、应用协方差传播律。 设观测向量的t个非线性函数为： 对上式求全微分，得nt

33、tnnXXXfZXXXfZXXXfZ,212122211155南京林业大学土木工程学院令 则由误差传播定律得： nnttttnnnndXXfdXXfdXXfdZdXXfdXXfdXXfdZdXXfdXXfdXXfdZ22112222112212211111ntttnnntXfXfXfXfXfXfXfXfXfKdXdXdXdXdZdZdZdZ2122212121112121,TXXZZKKDD56南京林业大学土木工程学院3.2协方差传播律协方差传播律v由以上推导知，求非线性函数的方差协方差矩阵比求线性函数的方差协方差矩阵只多一个求全微分的步骤。v例题6v例题757南京林业大学土木工程学院6、应用

34、协方差传播律时应注意的问题、应用协方差传播律时应注意的问题（1）根据测量实际，正确地列出函数式；（2）全微分所列函数式，并用观测值计算偏导数值；（3）计算时注意各项的单位要统一；（4）将微分关系写成矩阵形式；（5）直接应用协方差传播律，得出所求问题的方差协方差矩阵。3.2协方差传播律协方差传播律58南京林业大学土木工程学院3.2协方差传播律协方差传播律v协方差传播律的应用1、水准测量的精度2、算术平均值的精度3、若干独立误差的联合影响4、平面控制点的点位精度点位方差：59222222222SSuSyxp南京林业大学土木工程学院 权的概念权的概念权是表征精度的相对指标，指观测值所占的比重，精度越

35、高，比重越大。 权的定义权的定义权与方差成反比权的意义，不在于其数值的大小，重要的是它们之间的比例关系。 示例示例1220iip3.3权及定权的常用方法权及定权的常用方法602222122022202120211:1:1:nnnppp南京林业大学土木工程学院3.3权及定权的常用方法权及定权的常用方法 权的特点权的特点（1）选定一个 ，即有一组对应的权；（2） 不同，权不同，但权之间的比值不变；（3）同一个问题中只能选一个 ，不能选多个，否则就破坏了权之间的比例关系。（4）只要事先给定观测条件，就可确定权的数值。 单位权中误差的概念单位权中误差的概念 权为1的观测值所对应的中误差，称为单位权中误

36、差，即 。 610000南京林业大学土木工程学院3.3权及定权的常用方法权及定权的常用方法 定权的常用方法定权的常用方法1、水准测量的权1）按测站数确定2）按路线长度确定2、同精度观测值之算术平均值的权62iiNCp nnNNNppp1:1:1:2121nnSSSppp1:1:1:2121CNpiiN：每段路线的测站数S：各水准路线的长度N：观测值的观测次数南京林业大学土木工程学院1、协因数与协因数阵、协因数与协因数阵 协因数即为权倒数。3.4协因数和协因数传播律协因数和协因数传播律632021iiiipQXXXXDQ201XXXXQD20nnnnnnXXQQQQQQQQQQ212222111

37、211特点：I 对称，对角元素为权倒数 II 正定 III 各观测量互不相关时，为对角矩阵。当为等精度观测，为单位阵。南京林业大学土木工程学院3.4协因数和协因数传播律协因数和协因数传播律即有：协因数阵也称为权逆阵。64nnLLLLPPPDQ1.00.0.10.01.000.0.012120220222021201,LLLLLLLLQPQPnLLPPPP.00.0.0.021南京林业大学土木工程学院3.4协因数和协因数传播律协因数和协因数传播律2、协因数传播律、协因数传播律称为协因数传播律，或权逆阵传播律。与协方差传播律合称为广义传播律。3、权倒数传播律、权倒数传播律65TXXZZKKQQ0K

38、KXZ0FFXYTYYTXXrrYYQFFQQTYZTXXrtZYQFKQQnnZnPLfPLfPLfPLLLfZ1)(1)(1)(1),(222212121南京林业大学土木工程学院3.4协因数和协因数传播律协因数和协因数传播律全微分例题10 ：算术平均值之权等于观测值之权的n倍。 例题11 ：带权平均值的权等于各观测值权之和。例题12： 66nnZpkpkpkpKdLdZ11112222121nppXniiXpp121KXZFXYTXXYZKFQQ21南京林业大学土木工程学院4、单位权中误差的计算、单位权中误差的计算 用不同精度的真误差计算单位权中误差的公式如下： 实际应用1）由三角形闭合差

39、求测角中误差67npniiin120limnpniii120nniin12lim3nnii312菲列罗公式南京林业大学土木工程学院本章小结：本章小结：1、方差、方差协方差矩阵的定义协方差矩阵的定义2、协方差传播律、协方差传播律（线性和非线性）3、应用协方差传播律所应注意的问题、应用协方差传播律所应注意的问题4、权与定权的常用方法、权与定权的常用方法5、协因数和协因数传播律、协因数和协因数传播律TXXnnnnnXXDD2212221211221TXXZZKKDDCh3协方差传播律及权协方差传播律及权68南京林业大学土木工程学院 测试题测试题3-1 已知单位权方差为 、观测值 的权矩阵为试求：1、

40、 的方差2、 的方差3、 与 的协方差420200042026P523211LLLF22F2F1F21FF12310,10,5LLL2222112132233111111105510510FLL LL LLL LL69南京林业大学土木工程学院测试题测试题3-2 某地块由一梯形和一个半圆形组成，如图所示。已知观测值a=12m、b=8m、c=10m的方差协方差矩阵为：试求该地块的面积S的方差 。（注：取 ）2S32401103110000114LLDm70南京林业大学土木工程学院71南京林业大学土木工程学院721平差函数模型平差函数模型2平差数学模型平差数学模型3参数估计与参数估计与最小二乘最小二

41、乘条件平差间接平差附参数条件平差附条件间接平差随机模型函数模型函数模型线性化参数估计最优估计的性质最小二乘原理Ch4最小二乘原理最小二乘原理南京林业大学土木工程学院734.1平差函数模型平差函数模型函数模型函数模型 描述观测量与未知量间的数学函数关系的模型 目的：最优估计函数模型的未知量 函数模型如何构造？必要观测、多余观测必要观测、多余观测1 1）确定平面三角形的形状确定平面三角形的形状观测三个内角的任意两个即可,称其必要元素个数为2，必要元素有 种选择23C南京林业大学土木工程学院74必要观测、多余观测必要观测、多余观测2）确定平面三角形的形状与大小6个元素中必须有选择地观测三个内角与三条

42、边的三个元素，因此，其必要元素个数为3。任意2个角度+1个边、2个边+1个角度、三个边。4.1平差函数模型平差函数模型s1s3s2南京林业大学土木工程学院75必要观测、多余观测必要观测、多余观测3）确定如图四点的相对高度关系必须有选择地观测6个高差中的3个，其必要元素个数为3。h1、h5、h6或h1、h2、h3或h1、h2、h4等ADCBh1h6h5h2h4h34.1平差函数模型平差函数模型必要观测必要观测: 能够唯一确定一个几何模型所必要的观测 一般用t表示。特点: 给定几何模型，必要观测及类型即定，与观测无关。 必要观测之间不能有任何函数关系，即相互独立。南京林业大学土木工程学院76多余观

43、测: 观测值的个数n与必要观测个数t之差 一般用r表示，r=n-t。观测值: 为了确定几何模型中各元素的大小进行的实际 观测，称为观测值，观测值的个数一般用n表示。nt,，可以确定模型，还可以发现粗差。4.1平差函数模型平差函数模型南京林业大学土木工程学院必要观测可以唯一确定模型，其相互独立。可见若有多余观测必然可用这t个元素表示，即形成r个条件。123tnrtn180ADCBh1h6h5h2h4h3336tnrtn0621hhh0432hhh0546hhh实际上：1801804.1平差函数模型平差函数模型南京林业大学土木工程学院78函数模型函数模型1 1、条件平差、条件平差 以条件方程为平差

44、函数模型的平差方法2 2、间接平差、间接平差 以观测方程为平差函数模型11010rrnnrALA1110rrnnrWA111nttnndXBL111nttnnlxB或或4.1平差函数模型平差函数模型南京林业大学土木工程学院794.1平差函数模型平差函数模型3 3、附有参数的条件平差、附有参数的条件平差 以含有参数的条件方程为平差函数模型4 4、附有条件的间接平差、附有条件的间接平差 以观测方程和约束参数的条件方程为平差函数模型110110ccuucnncAXBLAtusWxClxBVsxsusnuunn,01111111110ccuucnncWXBA或南京林业大学土木工程学院条件方程的综合形式

45、为：条件方程的综合形式为：),(1 ,1 ,1 ,uncXLFF 为了线性化，取X的近似值：0X取 的初值： LLxXX0LL将F按台劳级数在X0，L处展开，并略去二次以及以上项：4.2平差平差数学模型数学模型非线性函数线性化非线性函数线性化南京林业大学土木工程学院xXFLFXLFxXLFFXLXL00,0),(),(0,212221212111,XLnnnnnncLFLFLFLFLFLFnLFLFLFLFA0,212221212111,XLunnnuuucXFXFXFXFXFXFXFXFXFXFBBxAXLFxXLFF),(),(0模型形式与线性函数类似。南京林业大学土木工程学院82用平差值

46、代替真值11010rrnnrALA1110rrnnrWVA11110ccuucnncWxBVA111nttnnlxBVtusWxClxBVsxsusnuunn,01111112020PQD4.2平差数学模型平差数学模型南京林业大学土木工程学院834.2平差数学模型平差数学模型函数模型随机模型测量平差数学模型平差数学模型是平差函数模型和随机模型的综合体。表达模型并用于求未知量最佳估值用于评定精度方差 协方差南京林业大学土木工程学院4.3参数估计与最小二乘参数估计与最小二乘v不论何种平差方法，平差最终目的都是对参数和观测量作出某种估计，并评定其精度，统称为对平差模型的参数进行估计。一、参数估计测量

47、平差的参数估计，是要在众多的解中，找出一个最为合理的解，作为最终估计。最终估计值应具有最优的统计性质。84南京林业大学土木工程学院4.3参数估计与最小二乘参数估计与最小二乘二、最优估计的性质1、无偏性 为参数 的估计量，若有 则称 是 的无偏估计量2、一致性 若估计量同时满足 则称 是 的严格一致估计量)(E1)(limP0)(lim)(2EEn85南京林业大学土木工程学院4.3参数估计与最小二乘参数估计与最小二乘3、有效性 若 是 的无偏估计，具有无偏性的估计量并不唯一。如果对于两个无偏估计量 和 ，有 则称 比 有效。 若 此时为最有效估计量。12)()(21DD12min)(D86南京林

48、业大学土木工程学院4.3参数估计与最小二乘参数估计与最小二乘三、最小二乘原理例：作匀速运动的质点在时刻 的位置是 ，函数如下： 在不同时刻 测定质点位置，得一组观测值 由运动方程可得： 或用图解表示如图： yyYXBVn.,21nyyy.,21iiiyvoyiiyy87南京林业大学土木工程学院4.3参数估计与最小二乘参数估计与最小二乘 从图中看到，由于存在观测误差，由观测数据绘出的点不成直线。 采用什么准则对参数 和 进行估计，从而使估计直线最佳地拟合于观测点？ 一般应用的是最小二乘原理，使各观测点到该曲线的偏差的平方和达到最小，即： 或 或 满足上式的估计称为最小二乘估计，此方法称为最小二乘

49、法。min12niiVminVVTminPVVT88南京林业大学土木工程学院89南京林业大学土木工程学院Ch5 条件平差条件平差条件平差原理条件平差原理条件方程的列立条件方程的列立条件平差精度评定条件平差精度评定公式汇编及示例公式汇编及示例90南京林业大学土木工程学院v在测量中，为了能及时发现错误和提高测量成果的精度，常做多余观测。v如果一个几何模型中有r个多余观测，就产生了r个条件方程，以条件方程为函数模型的平差方法，就是条件平差。5.1条件平差原理条件平差原理91南京林业大学土木工程学院5.1条件平差原理条件平差原理v函数模型 或v随机模型v未知数个是n，由于nr，所以条件方程是不定方程，

50、如何求解？v平差准则 11010rrnnrALA1110rrnnrWVA0ALAW12020PQDminTPVV如何求如何求V值？值？92南京林业大学土木工程学院5.1条件平差原理条件平差原理一、基础方程及其解一、基础方程及其解按求函数极值的拉格朗日乘数法，构造新函数求偏导，并令其等于0，以获得极值 导出 改正数方程1、基础方程 )(2WAVKPVVTT022AKPVdVdTTKQAVTKQAVWAVT0（1） 两式称为条件平差的基础方程（2）n+r 个方程，n+r个未知数93南京林业大学土木工程学院5.1条件平差原理条件平差原理2、法方程及其组成 令3、求解并计算观测量最佳估值01WKAAP