平差数学模型与最小二乘原理电子教案
- 格式:doc
- 大小:709.00 KB
- 文档页数:54
2.1 测量平差概述1、函数模型是怎样定义的?试举例说明函数模型的作用。
2、以下图2-1为例,说明必要元素、多余观测及其观测量的概念及其三者之间的关系。
图2-13、平差值和改正数是怎么定义的? 2.2 测量平差的数学模型1、试按条件平差法列出下列图形的函数模型已知点:A 、B 已知点:A 、B 观测值:31~ββ、1S 、2S 观测值:51~h h5h2、试按附有参数的平差法列出下列图形的函数模型已知点:A 、B 已知点:A 、B 观测值:41~h h 观测值:61~L L平差参数:C 点的高程 平差参数:角度ABC ∠、DBC ∠3、试按间接平差法列出下列图形的函数模型已知点:A 已知点:A 、B 、C 观测值:61~h h 观测值:31~S S平差参数:B 、C 、D 点的高程 平差参数:P 点的坐标4、试按附有限制条件的间接平差法列出下列图形的函数模型 已知值:矩形的对角边S 已知点:A 观测值:41~L L 观测值:41~h h平差参数:321~~~L L L 、、 平差参数:43~~h h H B 、、5、在下图所示的水准网中,A 为已知点,B 、C 、D 、E 为待定点,观测了9条路线的高差91~h h ,列出下列四种情况下的函数模型,并指出方程的个数。
(1) 条件平差法的函数模型;(2) 选取B 、C 、D 三点的高程平差值为参数; (3) 选取51~h h 的高差平差值为参数; (4) 选取85~h h 的平差值为参数。
6、试用表格的形式总结四种基本平差方法函数的异同。
7、四种基本平差方法的随机模型是什么?有什么作用?8、同精度观测了下图中的5个角度i L ,A 、B 为已知点,C 点为待定点,CD 边的方位角CD 为已知方位角,试列出条件平差的函数模型。
9、在下图所示的直角三角形中,我们观测了三角形的三个边长1L 、2L 、3L ,选取边长1~L 、2~L 为平差参数,试列出间接平差的观测方程。
授课题目:第二章 平差数学模型与最小二乘原理教学方法:理论讲授 教学手段:多媒体课件教学;以电子课件为主,投影及板书相结合为辅,使学生能够充分利用课堂有效的时间了解尽可能多的相关知识。
本章教学时数:4学时内容提要:主要介绍必要观测、多余观测、不符值、独立参数概念;测量平差的函数模型及两种平差的基本方程:条件方程和误差方程式;其它函数模型:附有参数的条件平差、附有限制条件的间接平差,以及平差的随机模型的概念及形态;平差基本方程的线性化,最小二乘原理。
教学要求:理解必要观测、多余观测、不符值、独立参数概念,掌握条件方程和误差方程式含义和最小二乘原理,会进行平差基本方程--条件方程和误差方程式的线性化。
本章重点:重点掌握测量平差数学模型的类型、建立方法,平差随机模型的意义和形态,以及最小二乘原理在测量平差中的应用。
教学难点:教学难点是对平差函数与随机模型含义与建立方法的理解。
本章教学总的思路:地理空间几何图形内部存在着严格的数学关系,测绘获得的是地理空间几何图形的基本元素,如角度(或方向值)、边长、高差的最佳估值,必须满足地理空间几何图形的基本数学关系,这是建立测量平差基本方程--条件方程和误差方程式的基础,在讲清楚这一点的基础上讲解基础方程的建立,进而推开讲解附有参数的条件方程、附有限制条件误差方程模型,并说明平差的随机模型的概念。
为解算的需要必须线性化条件方程式和误差方程式,其基本方法是利用泰勒级数展开基本方程并取其至一次项,从而完成线性化;在解释天然的平差模型为什么没有唯一解的原因基础上,讲解最小二乘原理,并举例验证,以此突破本课程难点内容的教学。
最后对教学重点内容作概括性总结,使学生加深理解与认知的程度。
§1测量平差概述本节教学时数:0.5学时本节重点:(1)测量元素-—角度(方向)、长度、高差、几何图的数学关系(2)观测值个数、必要观测数、多余观测数及其作用;(3)观测值、改正数、最优改正数、最优估值,平差的概念本节教学思路:以日常生活中最常见到的简单几何图三角形为例,说明测量观测值、平差值、几何图数学关系,平差模型与平差的概念,为下一节的讲讲解作好知识铺垫。
2平差数学模型与最小二乘原理2.1 参数估计及其最优性质几何模型:包括水准网和平面控制网(包括测角网、测边网、边角网)。
每种几何模型都包含有不同的几何元素,如水准网中包括点的高程、点间的高差,平面网中包含角度、边长、边的坐标方位角以及点的二维或三维坐标等元素。
这些元素都被称为几何量。
在诸多几何量中,有的可以直接测量,有的是间接求出。
几何模型不同,它所需要知道的元素的个数与类型也不同,目标是确定几何模型的唯一性。
1.如图2-1的三角形ABC中,为了确定它的形状,只需要知道三个内角中的任意两个内角的大小就可以了。
它们都是同一类型的元素。
2.要确定该三角形的大小和形状,就必须知道三个不同的元素,即任意的一边两角、任意的两边一角或者是三边。
它们中间都至少包含一条边长该情况包含角度和边长两类元素。
3.要确定该三角形的大小、形状和它在一个特定坐标系中的位置和方向,则必须知道图中15个元素(6个坐标元素,3个内角元素,3个边长元素,3个方位角元素)中的6个不同的元素,这6个元素可以构成更多的组合,至少要包含一个点的坐标和一条边坐标方位角,它们的改变只相当于整个网在坐标系中发生了平移和旋转,并不影响该三角形的内部形状和大小。
如果A、B两点都是已知点,为确定三角形的大小、形状、位置和方向,则只需要任意两个元素就行了,如两角、两边或一边一角等。
我们把能够唯一地确定一个几何模型所必要的元素,称为必要观测元素。
必要观测个数用t表示。
例如,确定三角形的形状,必要观测元素个数t=2;确定三角形的大小和形状,必要观测元素个数t=3;确定三角形的大小、形状、位置和方向,必要观测元素个数t=6。
对于后两种情况,不仅要考虑必要观测元素的个数,还要考虑到元素的类型,否则就无法唯一地确定模型。
必要起算数据个数用d表示,水准网为1,测角网为4,测边网和边角网为3。
观测值个数用n个表示。
当n <t 时,显然无法确定模型的解;当n =t 时,则可唯一地确定该模型,但对观测结果中含有的粗差和错误都将无法发现;当n >t 时,能及时发现测量中的粗差和错误,提高观测成果的精度和可靠性。
令多余观测个数t n r-=,在统计学中称r 为自由度。
一个几何模型能通过t 个必要而独立的量唯一的确定下来,当模型中有r 个多余观测量,一定存在着r 个这样的函数关系式。
例如在上述2中,如果观测了角度1L 、2L 、3L ,即n =3,t =2,则r =1,它们的真值之间存在如下关系式0180~~~321=-++︒L L L有r 个多余观测,就会有r 个这样的关系式(条件方程)。
由于观测不可避免地含有误差,所以0180321≠-++︒L L L 为了消除矛盾,通常用另一组被称为“观测值估值”(又叫平差值、最或是值、最或然值)Lˆ来代替观测值L ,即 ii i V L L +=ˆ i V 称为观测值的改正数,未知数个数>方程式个数,无数多组,所以问题的关键点是:在无数多组解中求得唯一的一组最优改正数。
测量平差的任务就是对参数(未知数)及其方差(协方差)进行估计,即对平差数学模型的参数进行估计(点估计和区间估计)。
由于多余观测而产生的平差数学模型,都不可能直接解方程而求得唯一解,测量平差中的参数估计,就是要在无数多组解中,找到一组最优的解作为平差参数的最终估计,为此,必须对平差数学模型附加某种约束条件,实现满足最优性质的参数唯一解,其中最广泛采用的平差准则是最小二乘准则。
最优估计量主要有以下3个性质。
1.一致性满足0])ˆ[(lim 2=-∞→X X E n 的估计量Xˆ为参数X 的一致性估计量。
2.无偏性 满足X XE =)ˆ( 则称Xˆ为X 的无偏估计量。
同时满足X XE =)ˆ( 0])ˆ[(lim 2=-∞→X X E n则称Xˆ为X 的严格一致性估计量。
3.有效性具有无偏性的估计量并不唯一,但毫无疑问方差最小的估计量是最优的。
设有估计量Xˆ1和Xˆ2,如果D(Xˆ1)<D(Xˆ2),则称Xˆ1比Xˆ2有效。
其中D(Xˆ)=min,为X的最有效估计量,称为最优无偏估计量。
数理统计理论证明,具有无偏性、最优性的估计量一定是一致性估计量,所以测量平差中参数的最佳估计值要求是最优无偏估计量。
2.2 最小二乘原理在测量工作及其它科学工程领域,应用最早也最广泛的就是所谓的“最小二乘准则”min =∆∆=ΦP T测量工作中习惯上用符号V 代替∆ˆ min ==ΦPV V T当P 为非对角阵,表示观测值相关,按min =PV VT进行的平差称为相关观测平差。
当P为对角阵,表示观测值不相关,此时最小二乘准则可表示为纯量形式,即min 2222211=+++==Φnn Tv p v p v p PV V特别地,当观测值不相关且等精度时,权阵P 为单位阵,此时最小二乘准则可表示为min 22221=+++==ΦnTv v v PV V其实,估计的准则有许多种,最小二乘准则是其中的一种,还有一种常用的估计叫做最大似然估计,其概率分布密度函数为⎭⎬⎫⎩⎨⎧∆∆-=∆=Φ-121221exp )2(1)(D Df T n π所谓极大似然估计,就是要在概率分布密度函数达到极大的条件下来对真误差∆进行估计。
显然,当)(1∆∆-D T达到极小时,概率分布密度函数可取得极大值,仍用V 表示对∆的估计结果,即要求:min 1==Φ-V D V T相当于min ==ΦPV V T显然,当观测向量服从正态分布时,极大似然估计与最小二乘估计的结果是一致的。
例[2-1] 设对某量X ~进行了n 次同精度独立观测,得观测值1,n L ,试按最小二乘准则求该量的估计值。
解:设该量的估计值为xˆ,误差方程式为 n n L xv L x v L xv -=-=-=ˆ................ˆˆ2211 写成矩阵形式1121211ˆˆ111⨯⨯⨯-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=n n n n n L xB L L L x v v v V按最小二乘准则,顾及E P =,得min ==V V PV V TT将上式对x ˆ取一阶导数,并令其为零,得211122ˆ2ˆ1==⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡===∑ni T T T T v V B V x d dVV x d V dV 将i i L xv -=ˆ代入上式得 0ˆ)ˆ(111=-=-=∑∑∑nininiL x n L xv解得n L L n x n i ][1ˆ1==∑2.3 测量平差的数学模型用数学关系式来描述对象的某种特征或内在的联系,这种数学关系式就称为数学模型。
在研究任何平差方法时,平差数学模型由函数模型和随机模型组成。
1.函数模型函数模型是描述观测量与待求量之间的数学函数关系的模型。
不同的函数模型有与之对应的平差方法。
函数模型有理论模型与实用模型之分,理论模型中的待估量用x X L ~,,~,~∆表示,实用模型中的待估量用x v X Lˆ,,ˆ,ˆ表示。
函数模型有线性与非线性之分,测量平差通常是基于线性函数模型,当函数模型为非线性时总是要将其线性化。
下面简述各种经典平差方法的线性函数模型。
1)条件平差法在图中,观测了三个内角,n =3,t =2,则r =n -t =1,存在一个函数关系式,也称为条件方程,可以表示为180~~~)~(321=-++=︒L L L L F令[]111A 31=⨯, []TL L L L 32113~~~~=⨯, ][-180=A 0︒ 则上式为0~0=+A L A一般而言,如果有n 个观测值1⨯n L ,必要观测个数为t ,则应列出r =n -t 个条件方程,即0)~(1=⨯L F r如果条件方程为线性形式,则可以直接写为~101=+⨯⨯⨯r n n r A L A 将∆+=L L ~代入,并令)(0A AL W +-=则0=-∆W A式即为条件平差的函数模型。
以此模型为基础的平差计算称为条件平差法。
2)间接平差法(参数平差法、坐标平差法)一个几何模型可以由t 个独立的必要观测量唯一的确定下来,因此,平差时若把这t 个量都选作参数,即t u =(这是独立参数的上限),那么通过这t 个独立参数就能唯一地确定该几何模型。
选择几何模型中t 个独立量为平差参数,将每一个观测量表达成所选参数的函数,共列出n t r u r =+=+个这种函数关系式,以此作为平差的函数模型的平差方法称为间接平差。
如图三角形ABC 中,观测了三个内角1L 、2L 、3L ,3=n ,2=t ,1=-=t n r ,平差时选A ∠、B ∠为平差参数21~~X X 、,即[]T X X X 21~~~=,2=u ,共需列出3=+u r 个函数关系式,列立方法是将每一个观测量表达成所选参数的函数,由图知:⎪⎪⎭⎪⎪⎬⎫+--===︒180~~~~~~~2132211X X L X L X L 方程的个数恰好等于观测值的个数。
令[]T X X X 21~~~=,[]T L L L L 321~~~~=,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=111001B ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=︒18000d 则可写为 13122313~~⨯⨯⨯⨯+=d X B L 一般而言,如果某一平差问题中,观测值个数为n ,必要观测个数为t ,多余观测个数为t n r -=,再增选u 个独立参数1~⨯u X ,t u =,则总共应列出n u r c =+=个函数关系式,其一般形式为)~(~1X F L n =⨯如果这种表达式为线性的,一般为111~~⨯⨯⨯⨯+=n t t n n d X B L 将∆+=L L ~和x X X o ~~+=代入上式,并令 )(d BX L l o+-=则可写为111~⨯⨯⨯⨯-=∆n t t n n l x B 以上就是间接平差的函数模型。
3)附有参数的条件平差法在平差问题中,设观测值个数为n ,必要观测个数为t ,则可以列出r =n -t 个条件方程,现又增设了u 个独立量作为未知参数,且t u <<0,每增加一个参数应增加一个条件方程,因此,共需列出u r +个条件方程,以含有参数的条件方程为平差函数模型的平差方法,称为附有参数的条件平差法。
如图2-2的三角形ABC 中,观测了三个内角1L 、2L 、3L ,3=n ,2=t ,1=-=t n r ,平差时选A ∠为平差参数X ~,即1=u ,此时条件方程个数应为2=+u r 个,它们可以写成:0180~~~321=-++︒L L L0~~1=-X L令⎥⎦⎤⎢⎣⎡=001111A ,⎥⎦⎤⎢⎣⎡-=10B ,⎥⎦⎤⎢⎣⎡-=︒01800A 则上式可写成0~~12011121332=++⨯⨯⨯⨯⨯A X B L A一般而言,在某一平差问题中,观测值个数为n ,必要观测个数为t ,多余观测个数为t n r -=,再增选u 个独立参数,t u <<0,则总共应列出u r c +=个条件方程,其一般形式为0)~,~(1=⨯X L F c 如果条件方程是线性的,其形式为0~~1011=++⨯⨯⨯⨯⨯c u u c n n c A X B L A 将∆+=L L ~和x X X o ~~+=代入上式,并令 )(0A BX AL W o++-=则得11110~⨯⨯⨯⨯⨯⨯=-+∆c c u u c n n c W x B A 为附有参数的条件平差的函数模型。