秩亏自由网平差的研究
刘 阳
(江苏师范大学,城建学部,江苏 徐州 )
摘要:秩亏自由网是因为控制网中没有足够的起始数据, 即缺乏基准的平差问题,因此按间接平差进行平
差时, 其误差方程的系数阵 B 不能满足列满秩的要求, 相应的法方程系数阵Nbb?BTPB 是秩亏阵.
为了求定未知参数的唯一确定解, 除了遵循最小二乘准则外, 还需增加新的基准约束条件 , 从而得到未知参数的唯一确定解.
本文主要利用MATLAB从传统的测量平差的观点出发, 来计算例题,分析,和论述亏秩自由网平差
之解的性质,讨论了附加矩阵S的形式了确定的方式,讨论了秩亏自由网平差之解与传统自由网平差之解的关系, 给出了详细的解答过程,并且比较了俩种方法的各自的优缺点,给出总结。 关键词:秩亏自由网;平差;间接平差
Research Rank Defect Free Network
Adjustment
Liuyang
(School of Urban construction and design, Jiangsu Normal University, 221116)
Abstract:Rank Defect Free Network control network because of not enough initial data,That lack of adjustment problems benchmark.Therefore, when carried out by indirect adjustment adjustment, the coefficient matrix B error equation does not meet the requirements of full rank.Corresponding normal equation coefficient matrix is rank deficient matrix.In order to find a unique set of unknown parameters to determine the solution, in addition to following the least squares criterion, the need to add a new benchmark constraints, resulting in a unique solution to determine the unknown parameters.
The main advantage of MATLAB article from the traditional viewpoint of Surveying Adjustment,Analysis of the nature of the calculation examples, and discusses the loss of rank free net adjustment of the solution,Additional discussion of the form of the S matrix determined, discusses the relationship between solutions of rank defect free network adjustment of the solution with the traditional free network adjustment, the process gives a detailed answer, and compare the two methods of their advantages and disadvantages.Gives summary.
Key words: Rank-defect free net adjustment; adjustment; condition comparison
引言
在现代测量数据处理过程中,秩亏自由网平差在近几十年得到了广泛应用,是重要的数据处理方法之一,特别是在变形监测、最优化设计中,秩亏自由网平差都展现出其优势。
1 秩亏自由网平差
1.1 秩亏自由网平差的提出
在经典间接平差中,必须有足够的起算数据.当控制网中仅含必要的起算数据,通常称为自由网.用经典方法平差这种网,俗称经典自由网平差.当控制网除必要的起算数据,还有多余的起算数据的网称为附合网,在间接平差时,不论是自由网还是附合网,当所选的参数不存在函数关系时,误差方程系数矩阵B总是列满秩的,即R(B)=t(t为必要观测).由此得到的法方程系数阵的秩
R(Nbb)?R(BTPB)?R(B)?t 法方程具有唯一解
.
在图(1)水准网中,假定P3的高程已知为H3,待定点P1、
??X0?x??X0?x?1,X?2.各段路P2的高程平差值为X1122线长度为S,高差为等权观测,误差方程
V?3?13?22?1??l 图(1) Bx3?1?v1??10??l1??x????11??1???l? v的显式为??2????x?2??2??????v3????0?1???l3??法方程及其显式为
??BTl BTBx?1??w1??2?1??x ????w? ??x??12???2??2?在误差方程系数阵B中,存在一个二阶行列式不等于零,如
10?1,故B的?11秩R(B)=2,即B为列满秩阵.由此法方程系数的秩R(N)=R(B)=2,所以法方程
有唯一解为
??(BTB)?1BTl x这就是经典自由网平差情况.
上述间接平差函数模型还可以用下面方式组成:先设P3点的平差值
0??X0??x,?,?3?0作为参数的条件方程,然后另X3将xX?X333参与列误差方程,3于是其函数模型为
V
T1?3?3?1??l Bx3?33?13?1??0 Cx3?1式中CT??001?,其显式为
?1??l1??v1??10?1??x????110??x?2???l2? v ?2??????????3??v3????0?11????x???l3???1??x?3?0) ?2??0 (即x ?001??x????3??x?可见俩种模型等价,平差结果相同.
在这种情况下,误差方程的行列式等于零,即
10?1 ?110?0
0?11其中有二阶行列式不等于零,故R(B)=2,数2为网中必要观测数,B为秩亏阵,其列亏数d=3-2=1,表示缺少一个起始高程,因此给定条件式,转化成附有限制条件的间接平差问题,可求其唯一解.
秩亏自由网的法方程系数阵N奇异,即N?0,故N的凯利逆N?1不存在,法方程有无穷解.如何合理的求解这类平差问题,就是本文要讨论的秩亏自由网平差问题.
产生秩亏的原因是控制网中没有起算数据,所以d就是网中必要起算数据的个数。
1.水准控制网:水准观测值是两个控制点之间的高差。为了确定水准网中各水准点的高程,就必须至少有一个已知高程点作为全网起算数据。所以,水准网的基准数据的个数 d?1。
2平面测角三角网:要推导各待定点的坐标,就必须有一个起算点坐标,还需要一个起算方位和一条起算边,或者有两个起算点坐标。所以,测角三角网的基准数据个数d=4。
2.测边三角网、边角同测三角网和导线网:必要起算数据有一个点的两个坐标和一个方位角,即这三种控制网的基准数据个数d=3。
3.三维控制网:三维控制网需要一个起算点(X,Y,Z)、三个已知定向角(?X,?Y,?Z)和一条空间已知边长布S,即需要7个起算数据;如果三维控制网中的观测值包括空间边长,则必要起算数据的个数为6,即三维控制网的基准 ,或d=6(含空间边长观测值)。 d=7(不含空间边长观测值)
4.GPS控制网:为确定各待定点在空间坐标系中的三维坐标(X,Y,Z),GPS控制网的必要起算数据只为一个点的三维坐标,所以GPS控制网得基准数为d=3。
1.2 秩亏自由网附加条件平差原理
附加条件法的基本思想:由于网中没有起算数据,平差时多选了d个未知参数,因此在u个参数之间必定满足d个附加条件式,即在原平差函数模型中需要加入d个未知参数间的限制条件方程,从而可以按附有条件的间接平差法求解。
秩亏自由网平差的函数模型为 L?BX?Δ,n,1n,uu,1n,1R(B)?t?u (1-1)
B的列亏数d?u?t,随机模型为
22?1 D??0 Q??0P (1-2)
按最小二乘原理VTPV?min,P为非奇异,法方程具有无穷多组解。
一般,为了获得未知数参数的惟一解,给定基准条件为 STd?uu?uu?1??0 (1-3)Px

