
图一 控制基标串测示意图图一中各符号的意义为:Pi:第i个导线点(控制基标)的点名。Si:第i-1点到第i 点的距离(观测值与理论值相差微小,以观测值表示)。α0、αn:附合导线两端的已知方位角。β’i:第i个控制点上的转折角理论值;βi:第i个转折角平差后的观测值。ui:第i点沿线路法线方向的归化改正数,含正负号。正值表示向观测角一侧改正。vi:转折角改正数,vi =β’i -βi;5.1 归化改正原理: 已知转折角改正数vi(i=0,1,2…n)和观测边si(i =1,2…n),在Vi剩=0(Vi剩 为转折角改正数残差;i=0,1,2…n)或S (vi剩^2)=min(最小)的条件下,求各归化改正数ui(i =0,1,2…n)。5.2数学模型的建立:

为求得U,下面讨论④式解的情况:考察④式的系数矩阵B的秩有:R(B)<=(n-1)<n+1,知B不可逆,方程④不是唯一解;当R(B|U)=R(B),相容方程④有解,且有无穷多解。当R(B|U)≠R(B),矛盾方程④无解。 针对以上两种情形,为求得④式的最优解,引入工程数学的“广义逆”(g逆)概念。设B的广义逆矩阵为B▔,最小范数g逆为Bm▔,最小二乘g逆为Bl▔;B的Moore-Penrose广义逆为B 。则(1)当R(B|U)=R(B)<n+1时,相容方程④的最小范数解不唯一,它的一个解为 U=(1/ρ)Bm▔V ⑤此时,U结果唯一,且满足 ||U||(U的范数)=最小。亦即横向归化改正值的平方和最小。 对于等边导线,设S1=S2=…=Sn=S,q=1/S 则④式中的B可写为:

可见,等边导线归化改正只有唯一解⑤’,它是⑤的特殊形式,同样满足||U||=最小。此时R(B)=R(B|U)=n-1。 此外等边导线具有两个重要规律:即满足两个公式(③和③’)。③’为等边导线所特有,也可用于检验等边导线观测值是否含有粗差。③是所有导线具有的规律,归化改正数的残差向量也符合这一规律,所以③除能检验观测值是否含有粗差外,还可以检验归化改正结果是否正确。(2)当R(B|U) ≠R(B)时,矛盾方程④无解,但可求最优近似解,即最小二乘解(不唯一),其一个解为 U=(1/ρ)Bl▔V ⑥⑥可使方程④残差向量的范数最小,即||ρBU||=最小。也就是附合导线的转折角改正数的残差向量的平方和为最小。综合(1)、(2)两种情况,即无论方程④有解或无解,均可得到最优解,并可统一写为下式: U=(1/ρ)B V ⑦ 在(1)的情形下,⑦是④的一个最小范数解;在(2)情形下,⑦是④唯一的最小二乘最小范数解。④、⑤、⑥三式或④、⑦二式即为控制基标归化改正的数学模型。 不论对于Bm-、Bl-还是B +,求解均有两种基本方法:满秩分解法、初等变换法。但精确解算十分繁琐,编程较难实现。实际编程是将V的残差向量(即转折角剩余改正数)的所有分量归化改正到满足限差要求为止。实践表明,该法可达到秒级精度要求。6.程序流程图程序流程图见图二
图二 程序流程图6.工程实例 本例点数包括两端4个固定点共16个,边数包括两条起算边共15条。用程序按5”限差计算,其中,折角改正数由转折角理论值减去观测值得到。计算成果见表一。表一:归化改正计算成果表
7.结束语 该控制基标归化改正程序,经过多年在北京、伊朗德黑兰地铁基标测量工程中使用,不仅可满足现行规范和施工要求,而且提高了工程效率和质量。归化改正程序必将在今后城市地下铁道、轻轨交通工程的铺轨基标工作中发挥其应有作用。同时,本归化改正程序对精密线形工程的高精度放样,也有一定的实用和参考价值。[注1]点属性:1表示固定点;0表示未知点(参与归化改正计算)。[注2]变换方法为:从第1行起,后两行加到前一行,依次类推到倒数第3行;再从倒数第2行起,每行加后边一行,依次类推至第1行。参考文献:[1] 秦长利。地铁铺轨基标测量方法探讨。铁路航测,1999(3)[2]北京城建勘察测绘院.《地下铁道、轻轨交通工程测量规范》(GB50308-1999).中国计划出版社, 2000年5月第一版[3] 陈永奇.《工程测量学》. 测绘出版社, 1995年5月第二版[4]於宗寿 鲁林成。《测量平差》。测绘出版社,1983年6月第二版作者简介:董伟东(1975-),男,北京城建勘测设计研究院工程师。Email:dwd2002@sohu.com
