法方程。
令Ⅳｃｃ＝ｃⅣ品ｃＴ，由（７）、（８）两式可求得
●Ｘ●｜ｘ“－Ｘｎ
＃ｘ｛
Ｋ＝磁（ｃｚｖ；＂ｗ＋吼）
（９）
ｉ＝（以一聪矿磁ｃⅣ二）Ｗ一蚝－１ＬＴｚｖ∞－１吼（１０）
将未代入公式（３）可得ｙ，最后可求出观测量和参数
的平差值
￡＝Ｌ ４－Ｖ
（１１）
譬＝Ｘｏ ４－；
（１２）
３附有限制条件间接平差的虚拟观测算法
ｓｃｉｅｎｃｅ，２００６ｂ，１６（８）：８３３－８３７． １ １ Ｙａｎｇ Ｙ ａｎｄ Ｇａｏ Ｗ．Ａｎ ｏｐｔｉｍａｌ ａｄａｐｔｉｖｅ Ｋａｌｍａｎ ｆｉｌｔｅｒ
［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｃ，ｅｅｄｅｓｙ，２００６ａ，８０：１７７—１８３ ［１２］Ｙａｎｇ Ｙ ａｎｄ Ｇａｏ Ｗ．Ｉｎｆｌｕｅｎｃｅ ｃｏｍｐａｒｉｓｏｎ ｏｆ ａｄａｐｔｉｖｅ ｆａｃ—
［９］ Ｙａｎｇ Ｙ，Ｇａｏ Ｗ．Ｃｏｍｐａｒｉｓｉｏｎ ｏｆ ａｄａｐｔｉｖｅ ｆａｃｔｏｒ ｉｎ ｋａｌ— ｍａｎ ｆｉｌｔｅｒｓ ｏｎ ｎａｖｉｇａｔｉｏｎ ｒｅｓｕｌｔｓ［Ｊ］．Ｔｈｅ Ｊｏｕｒｎａｌ ｏｆ
Ｎａｖｉｇａｔｉｎ，２００５．５８：４７１－４７８．
［１０］Ｙａｎｇ Ｙ，Ｇａｏ Ｗ．Ａ ｎｅｗ ｌｅａｒｎｉｎｇ ｓｔａｔｉｓｔｉｃ ｆｏｒ ａｄａｐｔｉｖｅ ｆｉＩ－ ｔｅｒ ｂａｓｅｄ ｏｎ Ｐｒｅｄｉｃｔｅｄ ｒｅｓｉｄｕａｌｓ［Ｊ］．Ｐｒｏｇｒｅｓｓ ｉｎ ｎａｔｕｒａｌ
第３３卷第２期 ２００８年３月
测绘科学
Ｓｃｉｅｎｃｅ ｏｆ Ｓｕｒｖｅｙｉｎｇ ａｎｄ Ｍａｐｐｉｎｇ
ＶｏＬ ３３ Ｎｏ．２ Ｍａｒ．
附有限制条件间接平差的虚拟观测算法
赵海涛，郭广礼．查剑锋
（中国矿业大学环境与测绘学院，江苏徐州２２１００８）
【摘 要】提出了附有限制条件间接平差的虚拟观测算法，该算法计算简单，易于理解，讨论了虚拟观测值的权
的取值问题，最后用一算例验证了附有限制条件间接平差的虚拟观测算法的合理性和正确性。
【关键词】附有限制条件的问接平差；虚拟观测方程；虚拟观测值的权
【中图分类号】Ｐ２０７
【文献标识码】Ａ
【文章编号】１００９－２３０７（２００８）０２－００３３－０３
ＤＯＩ：１０．３７７１／ｊ．ｉｓｓｔｒ １００９－２３０７．２００８．０２．０１２
参考文献
［１］黄维彬．近代平差理论及其应用［Ｍ］．北京：解放军 出版社，１９９２．
［２］隋立芬，宋力杰．误差理论与测量平差基础［Ｍ］． 北京：解放军出版社。２００４．
［３］杨元喜．抗差估计理论及其应用［Ｍ］．北京：八一出 版社，１９９３．
［４］杨元喜．相关观测抗差最小二乘估计［Ｊ］．测绘通 报，１９９５，（３）：３６－３８．
在水准网（图１）中，Ａ和Ｂ是已知高程的水准点，矾 ＝５．０１６ｍ，以＝６．０１６ｍ，并设这些已知点高程无误差。图 １中Ｃ，Ｄ和Ｅ点是待定点。Ａ和Ｂ点高程、观测高差和相 应的水准路线长度见表１。试按附有限制条件的间接平差法 求各待定点的平差高程。
表１观测值与起始数据
路线号高銎？量，蒌船
ｌ
＋１．３５９
１．１
观测法与使用一般算法求得观测值的改正数若保留小数点
后两位，其结果相同；虚拟观测值的权取至１０００，１０００
时，使用虚拟观测法与使用一般算法求得观测值的改正数
若保留小数点后三位，其结果相同。由此可以看出，尽管
理论上虚拟观测值的权Ｐ’一∞，但是在使用虚拟观测法进
行计算时虚拟观测值的权取得足够大而且毋须取的非常大
此后即可采用间接平差的算法求解。
法方程为：
Ｎ，１８８ ｉ．一甲＝０
法方程系数Ⅳ么为满秩矩阵，解得
；．＝Ⅳ船。１矿
公式（１７）、（１８）中，胪船＝Ｂ胛册”，
（１６）ｘｌ
（１７） （１８）
万方数据
测绘科学
第３３卷
㈤黑㈦：活玎嬲∥朋。
将公式（１８）代入，公式（１６）可得ｖ，ｒ，即ｙ和ｙ，的
值，矿一定等于零。 观测量和参数的平差值：
两种算法求得的待定点的高程相同，巩＝６．３７４８ｍ，Ｈｏ ＝７．０２７９ｍ，％＝６．６１２１ｍ。
表２两种算法所求改正数对照表
改０正姗数Ｄ
一般 算法
设未知参数墨＝／／ｃ，丘＝‰，墨＝也，丘＝Ｊ｜ｌ，，墨 ＝丘，取各未知数的近似值为
霹＝６．３７５ｍ，雹＝７．０２５ｍ，霹＝６．６１３ｍ，冠＝１． ３５９ｍ，霹＝１．０１２ｍ。
ｔｏｒｓ ｏｎ ｎａｖｉｇａｔｉｏｎ ｒｅｓｕｌｔｓ［Ｊ］．Ｔｈｅ Ｊｏｕｒｎａｌ ｏｆ Ｎａｖｉｇａ· ｔｉｏｎ，２００５，５８（３）：４７１－４７８． ［１３］Ｙａｎｇ Ｙ，Ｈｅ Ｈ ａｎｄ Ｘｕ Ｇ Ａｄａｐｔｉｖｅｌｙ ｒｏｂｕｓｔ ｆｉｌｔｅｒｉｎｇ ｆｏｒ ｋｉｎｅｍａｔｉｃ ｇｅｏｄｅｔｉｃ ｐｏｓｉｔｏｎｉｎｇ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｄｅｓｙ， ２００１，７５（２／３）：１０９－１１６． ［１４］吴生武．自适应序贯抗差估计［Ｊ］．测绘通报，２００６， （１）：１４． ［１５］杨元喜．动态Ｋａｌｍａｎ滤波模型误差的影响［Ｊ］．测 绘科学，２００６，３１（１）． ［１６］王穗辉．顾及起算数据误差的附加基准平差［Ｊ］．大 地测量与地球动力学，２００５，２４（１）：７２—．７５． ［１７］陶本藻，许海威．大范围ＧＰＳ水准拟和模型误差的平 差补偿［Ｊ］．测绘通报，２００５，（７）：８－１０． ［１８］宋力杰．测量平差程序设计［Ｍ］．北京：解放军出版 社．１９９９．
Ｔｈｅ ｃｏｍｐｅｎｓａｔｉｏｎ ｗｉｔｈ ａｄｊｕｓｔｍｅｎｔ ａｎｄ ｉｎｆｌｕｅｎｃｅ ｏｆ ｔｈｅ ｏｕｔｌｉｅｒｓ ｉｎ ｔｈｅ ｐｒｉｏｒ ｐａｒａｍｅｔｅｒｓ Ａｂｓｔｒａｃｔ：Ｔｈｅ ｓｅｑｕｅｎｔｉａｌ ａｄｊｕｓｔｍｅｎｔ ｐｒｏｂａｂｌｙ ｍａｋｅｓ ｐａｒａｍｅｔｅｒｓ ａｎｄ ｔｈｅｉｒ ｐｏｓｔｅｒｉｏｒ ｅｏｖａｒｉａｎｃｅ ｕｎｒｅｌｉａｂｌｅ ｉｆ ｔｈｅｒｅ ａｒｅ ｏｕｆｌｉｅｒｓ ｉｎ ｔｈｅ ｐｒｉｏｒ ｐａｒａｍｅｔｅｒｓ．Ｉｎ ｔｈｉｓ ｐａｐｅｒ，ａｆｔｅｒ ａｎａｌｙｚｉｎｇ ｔｈｅ ｉｎｆｌｕｅｎｃｅ ｏｆ ｔｈｅ ｏｕｔｌｉｅｒｓ ｏｎ ｐａｒａｍｅｔｅｒｓ ｉｎ ｔｈｅ ｐｒｉｏｒ ｐａｒａｍｅｔｅｒｓ，ａｎｄ ｉｎ ｏｒｄｅｒ ｔｏ ｃｏｎｔｒｏｌ ｔｈｅｉｒ ｉｎｆｌｕｅｎｃｅ，ｔｈｅ ａｕｔｈｏｒ ｃｏｎｄｕｃｔ ａｎ ａｄｊｕｓｔｍｅｎｔ ｃｏｍｐｅｎｓａｔｉｏｎ ｍｅｔｈｏｄ ｉｎ ｄｅｔａｉＬ Ｔｈｅ ｅｘｔｒａ－ｐａｒａｍｅｔｅｒｓ ａｒｅ ａｄｄｅｄ ｉｎ ｔｈｅ ｆｕｎｅ－ ｔｉｏｎ ｍｏｄｅｌ ｉｆ ｏｕｔｌｉｅｒｓ ｉｎ ｔｈｅ ｐｒｉｏｒ ｐａｒａｍｅｔｅｒｓ ａ肥ｆｏｕｎｄ．Ａ ＧＰＳ ｅａｓｅ ｉｓ ｇｉｖｅｎ ｔｏ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｃｏｍｐｅｎｓａｔｉｏｎ ｍｅｔｈｏｄ Ｃａｎ ｃｏｎｔｒｏｌ ｔｈｅ ｉｎｆｌｕ·
本题ｎ＝７，ｔ＝３，ｒ＝４，ｕ＝５，列误差方程和限制条 件方程：
一１ １ Ｏ —ｌ ０ ｌ
． ｒｏ、
甄驯，
Ｏ ０ 一ｌ
巴： Ｏ 一１ ０Ｏ
萝 （？３）＝ｏ
以每公里观测高差为单位权观测，Ａ＝百１，观测值的
注：为对一般算法和虚拟观测算法所求结果进行精确比较，各改 正数取至小数点后四位。
由表２可以看出：
１）当虚拟观测值的权取至１０００００，１０００００时，虚拟
选取的虚拟观测值的权是否合适的依据：Ｖ’≠Ｏ，说明选取
的虚拟观测值的权偏小，应增大虚拟观测值的权，再求解；
矿＝０，说明虚拟观测值的权取值合适，所求未知参数的解
是正确的。为减少计算次数，虚拟观测值的权尽量取的大
一些。
４算例分析
取参考文献［１］例５－７的数据来验证附有限制条件的 间接平差虚拟观测算法的合理性和正确性。
求解结果就可以达到较高的精度。
（下转第３７页）
万方数据
第２期
曾安敏等参数先验信息异常的影响与平差补偿
３７
加参数的函数模型补偿方法，先对参数先验信息进行异常 判断。如某个参数先验信息有异常，则在误差方程中有针 对性地引入新的附加参数，从函数模型予以补偿，然后再 进行平差计算。采用附加参数的函数模型补偿法能很好地 控制异常信息对参数估值的影响，在控制异常先验信息的 同时，解出异常改正量，从而获得可靠的参数估值。
ｌ 引言
在一个平差问题中，多余观测数ｒ＝，Ｉ—ｔ，如果在平差 中选择的参数Ⅱ＞ｔ个，其中包含了ｔ个独立参数，则参数 间存在ｓ＝扯一‘个限制条件。平差时列出，１个观测方程和ｓ 个限制参数间关系的条件方程，以此为函数模型的平差方 法。就是附有限制条件的间接平差…。附有限制条件的间 接平差的一般算法，推导较为繁琐，公式较为复杂，本文 简述了计算简单，易于理解的附有限制条件的间接平差的 虚拟观测算法，并进行了算例分析。
￡＝Ｌ＋Ｖ
（１９）
譬＝ｒ＋；（２０、
附有限制条件间接平差的虚拟观测算法的关键是确定
虚拟观测值的权阵Ｐ’。Ｐ’内元素不相关，Ｐ’是一个对角阵。
因为每一个虚拟观测值的改正数等于零，方差等于零，所 以其权趋于无穷大，即．Ｐ７—００【４声Ｊ。在实际计算时，每一
个虚拟观测值的权取足够大的有限值即可求出所有未知参 数的正确解。Ｉ，’是否等于零（在一定精度内）可以作为判断
作者简介：赵海涛（１９８２．），男，山东 济宁人，硕士研究生，从事工程测量、 开采沉陷及其控制方面的研究。
Ｅ·ｍａｉｌ：ｚｈｔ２００２４１２＠１６３．ｏｏｍ
收稿日期：２００６．１２—１１
由公式（１３）、（１５）得总的误差方程
。。嚣。。。瓢。曼一。。量 （ｎ＋，）
（ｎ＋Ｊ）×Ⅱｎ ｘｌ
（ｎ＋１） ｘｌ
观测值改正数等于零（保留小数点后四位），使用虚拟观测
法与使用一般算法求得的观测值的改正数相同，说明使用
虚拟观测法进行附有限制条件的间接平差是合理的和正