大地测量中不适定问题的正则化解法研究
大地测量中病态性问题研究

大地测量中病态性问题研究高立恒(山东理工大学建工学院测绘工程山东淄博 255000) 摘要:目前,随着大地测量技术的迅速发展,在大地测量界对病态问题处理的精度和准确度的要求越来越高,人们所采用的数据处理手段也应随之发展。
如何充分考虑大地测量实际,探讨有偏估计和观测方程直接解算的新方法,以满足大地测量发展的需要,是大地测量学者面临的一个新任务。
关键词:大地测量;病态性问题;测量误差;数据处理引言在大地测量中,大部分实际问题都归结为解线性方程组问题。
当方程组系数矩阵呈病态时,测量平差系统所得参数估值性质明显变差。
病态性问题在测量工程中大量存在,而且有时是不可避免的,大地测量中是很常见的,如在GPS快速定位、某些控制网平差、大地测量反演、航天飞行器的精密轨道解算以及重力场的向下延拓等方面存在病态问题。
探索病态问题处理的新方法,是现代测量误差理论及数据处理研究中一项重要课题国内为许多的专家学者已经提出不少解决病态问题的理论及方法。
其中著名的有岭估计、广义岭估计、主成分估计法、遗传算法、加权迭代改善算法、奇异值分解法(SVD)等。
1 病态性产生的原因病态性在数据处理中危害性是很大的。
人们必须对病态性产生的原因和机理进行准确的分析和掌握,这将有助于对病态性进行正确的诊断,它是研究如何削弱和克服病态影响的基础。
病态性的产生与参数选取、观测以及计算方法的选择密切相关。
1.1参数选择的原因。
人们在建立—个平差模型时,在平差前,由于不可能对所有参数的规律全部了解掌握,往往造成过度参数化,使得参数之间存在着某些程度的复共线性,从而导致模型病态。
例如,在测边网中,要用多台测距仪测量边长,为了补偿不同仪器的系统误差,在数据处理中,经常要附加参数,这样常常会造成过度参数化,引起模型病态。
1.2 观测的原因。
观测的原因是指统计分析中的采样不足或测量当中观测量不够。
在统计分析中,由于某些条件限制,子样采样为局部结果。
此时设计矩阵表现出严重的复共线性,其实质为采样不足而引起的模型病态。
一种适用任意距离的大地主题解算方法的研究和试验

下,计算不收敛,无法解算出正确解。若两大地点同在赤
道,即B1=B2=0,由反算①和④式可得:“l=“2=“。=0,
由此反算⑤式无法计算,从而计算不收敛。此时由图1并根
据正算⑨式计算可知:o-,=0,o-m=曲/2。因此只需在反算
⑤式计算时首先判断u。是否为0,若为0则按o-m=曲/2计
算,否则还是根据反算⑤式计算%。
明显的缺陷u J,主要表现在: 1)应用范围受到一定的限制。尽管目前很多算法在理
论上适用于任意大地距离,但实际上受数值精度的影响, 它们几乎都只适用于一定范围而不能解算任意大地距离的 主题。
2)奇异问题。当两大地点位于同一子午圈或同在赤道 时,此算法不能计算,必须采用别的方法解算。
这些问题的存在,使得大地主题计算问题变得很复杂。 大地上的坐标点有远有近,如果所使用的解算公式不能计 算任意距离的正反算,则对于某些计算结果会产生很大误 差,或者选择不同的算法公式而使问题(特别是对于大批量 的计算)变得复杂化。下面介绍一种嵌套系数法,此算法的 基础是贝塞尔的两个微分公式,但算法不同,我们也可以 把它看作是贝塞尔的一种改进算法。
43 00 55.8771 118 10 02.999 32 218 11 26.795 88
43 00 55.878 4 11810 03.000 0 218 11 26.797
55 00 0.D00 0
9 59 59.999 9 216 45 07.391 55 00 00.000 79 9 59 59.997 20 216 45 07.389 67 55 00 00.00683 9 59 59.964 71 216 45 07.382
血=厶/(尼1b)+dka
(4)
公式(4)就是球面大圆弧如与大地线△之间的换算公
现代平差理论与测量平差

现代测量与现代平差技术摘要:本文首先简述了现代测量平差中的各种理论与经典测量平差之间的关系,指出现代测量平差与数据处理理论仍然是以高斯-马尔柯夫模型为核心,通过该模型在不同层面上的扩充、发展形成了若干新理论、新方法,并以图描述了经典测量与现代测量数据处理中各种理论之间的关系。
然后分别阐述了现代测量数据处理中粗差理论、系统误差的处理、病态问题的处理、非线性问题的处理、不等式约束的平差等的发展,最后综述了其他数据处理理论的一些发展情况。
最后讲了整体平差法是一个严格而又有效的平差方法,其应用与现代计算机技术密切相关。
具体介绍了整体平差法的基本原理,并以实测GPS控制网的布设为例,探讨了它在现代测量控制网建立中的具体应用及其技术优势。
关键词:经典测量平差;现代测量平差;高斯-马尔柯夫误差模型;误差模型扩展整体平差分级平差GPS控制网Abstract: This paper described the relationship between the theories in modern surveying adjustment and the traditional surveying adjustment. It pointed out that the theories of modern surveying adjustment and the data processing should be still based on Gauss-Markov error model. Through enlargement and development in different aspects of the model, new theories and methods are worked out. A figure showing such relationship is given.Meanwhile, the theories on blunder detection, systematic error processing, ill-pose problem, nonlinear model,inequality constraints are elaborated. At the last the progresses of other theories on data processing are summarized.Key words: traditional surveying adjustment; modern surveying adjustment; Gauss-Markov error model;extension of error model1、现代测量平差与数据处理理论发展概述经典的测量平差与数据处理是以高斯-马尔柯夫模型为核心[1]:L=AX+Δ(1a)E(Δ) = 0,D(Δ) =σ20Q=σ20P-1(1b)Rnk(A) =n,R(Q) =R(P) =n(1c)这里L为观测向量,Δ为误差向量,X为未知参数向量,A为X的系数矩阵,E(·)为数学期望,σ2为单位权方差,P为观测权矩阵,Q为协因素矩阵,n为观测个数。
大地电磁测深阶段式自适应正则化反演

2020年8月第55卷 第4期 *江西省南昌市经开区广兰大道418号东华理工大学地球物理与测控技术学院,330013。
Email:dwjhtj@hotmail.com本文于2019年12月27日收到,最终修改稿于2020年2月14日收到。
本项研究受国家自然科学基金项目“基于线性相关约束的地球物理联合反演研究”(41604104)和“圆锥型场源矿井瞬变电磁法理论基础研究”(41974086)、东华理工大学教育部核技术应用工程研究中心开放基金项目“宽范围岩石物性约束技术的改进及其在相山地区MT和重力正则化联合反演中的应用”(HJSJYB2016-7)及江西省自然科学基金项目“城市多参数可控电磁勘探理论研究”(20192BAB202012)联合资助。
·非地震·文章编号:1000-7210(2020)04-0906-09大地电磁测深阶段式自适应正则化反演郭一豪①② 陈 晓*①② 杨海燕①② 张志勇①② 曾志文①② 周 勇①②(①东华理工大学核技术应用教育部工程研究中心,江西南昌330013;②东华理工大学地球物理与测控技术学院,江西南昌330013)郭一豪,陈晓,杨海燕,张志勇,曾志文,周勇.大地电磁测深阶段式自适应正则化反演.石油地球物理勘探,2020,55(4):906-914.摘要 如何合理确定正则化因子一直是地球物理反演的研究热点和难点。
从提高反演稳定性的角度对正则化因子的选择进行考量,是确定正则化因子取值的新思路。
此外,以往的正则化反演研究对非线性优化算法的随机性考虑不足。
基于此,以Zhdanov提出的自适应算法为框架,提出了一种新的自适应正则化算法,即阶段式自适应算法,按照“阶段”自适应地调整正则化因子。
将此方法分别应用于大地电磁测深(MT)的共轭梯度和差分进化算法(DE)反演。
模型试验表明该算法可以提高反演的稳定性,在一定程度上降低衰减因子的影响,而且该算法具有同时适用于线性和非线性优化算法的特点。
大地电磁反演中改进的自适应正则化因子选取_向阳

A n I m r o v e d A d a t i v e R e u l a r i z e d P a r a m e t e r p p g S e l e c t i o n i n M a n e t o t e l l u r i c I n v e r s i o n g
] 4 5 - , 该方法能给出估计值和真 偏风险估计预测方法 [
: r a d i e n t A b s t r a c t B a s e d o n t h e c o n u a t e m e t h o d t o s o l v e t h e g j g ,a r o b l e m o f l a e r e d m a n e t o t e l l u r i c( MT)m o d e l s i n v e r s e d p y g c o m a r a t i v e s t u d i s m a d e b e t w e e n t h e i n v e r s i o n r e s u l t s a n d p y a r a m e t e r s e v e r a l o t h e r r e u l a r i z a t i o n s e l e c t i o n m e t h o d s u n d e r p g i n i t i a l m o d e l s i n o r d e r t o a n a l z e t h e i r r e s e c t i v e d i f f e r e n t y p b e h a v i o r a n d c o n d i t i o n s .T h e r e s u l t s s h o w t h a t t h e a d a t i v e p r e u l a r i z a t i o n e t a l o r i t h m c a n a s i m i l a r i n v e r s i o n r e s u l t i n a g g g , f a s t e r a n d e a s i e r w a .M o r e o v e r w i t h t h e u r o s e t o r e d u c e t h e y p p l i m i t d u e t o d e e n d i n o n t h e i n i t i a l m o d e l a n d i m r o v e t h e p g p , s t a b i l i t o f i n v e r s i o n a n i m r o v e d a d a t i v e r e u l a r i z a t i o n p p g y u t m e t h o d i s f o r w a r d b s e l e c t i n a l a r e r a t i o o f t h e m i s f i t p y g g f u n c t i o n t o t h e s t a b i l i z i n f u n c t i o n a s t h e i n i t i a l v a l u e o f g , a r a m e t e r a r e u l a r i z a t i o n n d a n a u t o m a t i c d e c a c o e f f i c i e n t p g y i s r o o s e d t o r e d u c e t h e a r a m e t e r w i t h i t e r a t i o n . s c h e m e p p p
GPS快速定位中相位模糊度动态解算的一种正则化方法

GPS快速定位中相位模糊度动态解算的一种正则化方法段荣;赵修斌;庞春雷;李媛;伍劭实【期刊名称】《中国惯性技术学报》【年(卷),期】2015(023)005【摘要】针对GPS快速定位中少数历元组成的法方程存在严重病态性的问题,研究了GPS单频整周模糊度快速解算的新方法:即首先采用改进SVD分解获得系数矩阵的精确奇异值,避免了较小奇异值抖动的影响;然后基于改进SVD分解结合法矩阵病态性特点,设计了一种改进Tikhonov正则化方法,合理构造了正则化矩阵,有效抑制法矩阵的病态性.算例表明:与LS估计-LAMBDA方法和Tikhonov正则化-LAMBDA方法相比,新方法能够有效降低法矩阵的条件数约3个数量级,仅解算4个历元数据,浮点解偏离真值的方差从41.89减小到1.04,可以快速获得准确、稳定的模糊度浮点解.模糊度固定性能结果进一步表明,新方法显著地提高模糊度搜索效率和成功率,解算成功率提高100%.【总页数】7页(P624-630)【作者】段荣;赵修斌;庞春雷;李媛;伍劭实【作者单位】空军工程大学信息与导航学院,西安710077;空军工程大学信息与导航学院,西安710077;空军工程大学信息与导航学院,西安710077;94153部队,陕西咸阳712200;空军工程大学信息与导航学院,西安710077【正文语种】中文【中图分类】V557【相关文献】1.GPS动态定位中整周模糊度的快速解算 [J], 刘少鹏;卢艳娥;周亚飞2.单频GPS动态定位中整周模糊度的一种快速解算方法 [J], 刘宁;熊永良;冯威;徐韶光3.GPS动态定位中相位模糊度的解算方法 [J], Sauer.,K;王爱朝4.GPS接收机RTK定位中整周模糊度的快速解算方法 [J], 冯仲科;韩熙春5.差分GPS载波相位整周模糊度快速解算方法 [J], 朱志宇;刘维亭;张冰因版权原因,仅展示原文概要,查看原文内容请购买。
不适定逆成像问题的多准则正则化求解方法

第39卷第1期20009-2月复旦学报(自然科学版)Journa]ofFudanUniversity(NaturalScience)V01.39No.1Feb.2000文章编号:042771r】|Iz()OO)Ot0026Oq不适定逆成像问题的多准则正则化求解方法戴伟辉,高汝熹(管理学院)搐要:不适定性(illposedness)是图像重建(x射线CT的投影重建、心电ECG厦脑电EEG信号的逆成像重建)逆问题求解中非常普遍的情况,由于对观测数据及正问题数字离散化过程中所产生的误差撮为敏感,它在很^程度上将影响到重建凰像的质量、稳定性和求解精度.通过对由投影重建图像逆问题求解过程中的不适定性丹析,给出r求解不适定逆成像问题的多准则正则化方法及其理论基础,实验结果表明诫方法此传统的单准则正则化求解方法具有更高的成像精度与稳定性.该方法已成功地应用于CT图像重建、NMR围像重建、心电逆同题求鹪以及大气湍流图像重建等领域.关键词:图像重建;不适定逆问题;多准则正则化中图分类号:TP39l41文献标识码:A图像重建过程中.由于逆算子的不适定性.成像质量和精度对观测数据及正问题数字离散化过程中所产生的误差极为敏感.解决这一问题的通常途径是采用Tikhonov和Miller形式的正则化求解方法“~.近年来,人们在图像重建领域内研究的新成果表明,一-:采取多个准则函数同时优化的正则化求解方法和某些更加精细而有效地描述待重建对象属性的先验知识时.可以获得精度大大高于常规方法下的重建图像,为求解不适定逆成像问题提供了有敬的途径.本文通过对由投影重建图像逆问题求解过程中的不适定性分析,给出了一组受ghost函数影响下的重建图像示例,阐述了求解不适定逆成像问题的多准则正则化方法及其理论基础.实验结果表明,该方法比传统的单准则正则化求解方法具有更高的成像精度与稳定性.并已成功地应用于CT图像重建、NMR图像重建、心电逆问题求解以及大气湍流图像重建等领域。
有界椭球不确定性平差算法在病态问题中的应用

有界椭球不确定性平差算法在病态问题中的应用夏玉国; 宋迎春【期刊名称】《《大地测量与地球动力学》》【年(卷),期】2019(039)009【总页数】5页(P956-959,965)【关键词】不确定性; 椭球约束; 病态问题; 岭估计【作者】夏玉国; 宋迎春【作者单位】中南大学地球科学与信息物理学院长沙市麓山南路932号 410083; 中南大学有色金属成矿预测与地质环境监测教育部重点实验室长沙市麓山南路932号 410083【正文语种】中文【中图分类】P207病态问题广泛存在于大地测量的数据处理中,如GPS 快速定位、大地反演中常存在不适定问题,其原因是观测信息不充分或者对参数物理信息和先验信息不了解[1-2]。
为解决病态问题,学者们提出一系列有偏估计方法,如岭估计法、Tikhonov正则化法、截断奇异值法,这些方法主要是利用数学手段缩小病态观测方程的系数阵(或法矩阵)的条件数,但很少进行病态原因的分析和从机制上采取有效措施,没有重视算法在实际问题中的物理意义是否合理[3-4]。
解决病态问题的有效途径是充分利用参数的附加信息或先验信息对部分参数进行约束,从而保证参数解的唯一性和稳定性,如等式约束、不等式约束、区间约束以及椭球约束等[5-8]。
等式约束和不等式约束在大地测量实际问题中是一种简单的约束信息,相关文献较多,如谢建等[6]建立附等式约束的病态模型,推导出参数估值的计算公式;冯光财等[7]给出不等式约束和广义岭参数之间的转换方法。
有学者提出利用椭球集合描述不确定先验信息,通过集员估计方法进行参数估计,如Fogel等[8]提出的椭球集员估计算法。
利用椭球约束描述复杂的先验信息,在测绘数据处理中是一种新的尝试,由于带有椭球约束的线性模型与岭估计有较强的相关性[7,9],因此用集员估计方法处理病态问题比较合适,此方法已经在其他领域得到广泛应用,但在测量平差数据处理中还需要进一步研究[10]。
本文将利用椭球约束描述观测向量和参数的有界不确定先验信息,把它们融入到测量平差模型中,建立平差准则对不适定性进行抑制,研究新的有界椭球不确定性平差方法在病态问题中的应用。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
大地测量中不适定问题的正则化解法研究
摘要:为了解决大地测量中的不适定问题,人们提出了正则化解法,并期望通
过对正则解法的不断研究从而彻底解决大地测量中的不适定问题。
论文对大地测
量中不适定问题的正则化解法研究进行详细论述,给相关人士提供参考。
关键词:大地测量;?不适定问题;?正则化解法;?系统误差;
大地测量是一项对地球的相关数据进行测量的活动。
大地测量活动的开展不
但可以有效提升地形测图以及工程测量的精准度,同时还可以促进国家空间科学
以及国防建设的发展。
此外,随着大地测量的不断深入,人们可以对地壳运动以
及地震等地质活动进行预测,从而降低地震等自然灾害对于人类的危害。
然而在
大地测量中,时常会遇到一些不适定问题。
例如,测量中所存在的控制网平差、GPS无法快速定位等。
这些大地测量中的不适定问题虽然表现形式不同,但却有
着一些相同点。
首先,这些不适定问题一般解均不唯一。
再者,这些不适定问题
有时还会出现无解的状况。
此外,这些不适定问题常常还会出现解不稳定的现象。
这些不适定问题的出现严重影响了大地测量的进行与发展,因此,为了解决大地
测量中的不适定问题,对其解决方法进行了深入的研究,并将其逐步演变为正则
化解法。
通过正则化解法,可以有效地解决大地测量中的不适定问题,并针对病
态性的算法进行改进,从而促进大地测量的快速发展。
1 推导了大地测量不适定问题解的统一表达
为对大地测量中不适定问题开展正则化解法研究,最初研究推导了大地测量
中不适定问题解的同意表达。
旨在分析大地测量中不适定问题常用的一些数学模型,研究表明在该阶段常见的数学模型主要有拟合推估模型、自由网平差模型、
病态模型和半参数模型等。
经计算显示,这些数学模型的解可以用某个数学关系
式统一表达,而令研究者所震惊的是这些数学模型都能够在TIKHONOV正则化原
理下推导出。
实际推导过程中,为保证计算结果的准确度,研究者要把握好这些
数学模型之间的共性问题,尽可能地分析出他们的个性,求解时既要考虑数学模
型的基本计算理论,又要寻求合适的优化求解方案,以此来深化研究。
2 克服病态性的改进算法研究
在克服病态性的改进算法研究中,从以下3步展开论述:首先,针对一些难
以确定的岭参数,系统会主动选择研究确定的岭参数L曲线。
为使L曲线的效果
能够更加清晰地展现出来,该算法研究采用对比法,将L曲线法同传统的岭迹法
相比较,以此来得出全新的结论。
其次,研究还提出了克服病态性的两步解法,
需重点研究了两步解法的计算原理和相关数据性质以及相应的计算适应条件等。
同常规的克服病态性改进算法研究方案相比,该方案更为优异。
最后,研究提出
了一种新的奇异值修正方案,该方案的核心是将奇异值分为2个部分进行分别修
正处理。
实践证明这种方案是很有研究效果的,同其他克服病态性的改进算法相
比该方案的结算结果更为精准。
3 单频GPS快速定位中减弱病态性的新方法研究
本次研究,主要论述了单频GPS快速定位中减弱病态性的新方法,能够在较
短的时间内实现快速GPS定位。
为此,首先分析了关于GPS快速定位的矩阵的结
构特性。
在正则化原理的前提下,有针对性地提出了以下2种正则化矩阵的构造
方法。
利用这2种新的方案,可以在很大程度上减弱传统法矩阵的病态性,利用
较短的时间就可以得出较为准确的结论。
为此,对这2种新型的减弱矩阵病态性
方案进行了列示:减弱方程病态性的MINEI方案、减弱法方程病态性的MINEII方案。
4 单频GPS快速定位中ARCE方法的改进
ARCE方案的提出是为了实现单频GPS的快速定位,该计算方法主要是以LS
估计、零空间的思想为理论基础而提出的。
实际运算中主要适用于单频接收机,
是一种快速结算整周模糊度的全新方案。
在过去ARCE方法的性能还不够完善,
进而在单频GPS快速定位中使用该方法只能适用于一些观测时间段位几分钟的情况,这将会在很大程度上阻碍GPS的定位工作。
为此,对该方案提出了具体的改
进措施,针对传统ARCE方法下,法矩阵病态性状况严重,引起检测数据结果不
可靠的现象,利用TIKHONOV正则化原理,设计了一种正则化矩阵的构造方法,
大大削弱了法矩阵的病态性,使得出的数据更为准确,大大缩小了模糊度的搜索
范围,之后利用ARCE方法可以结算出整周模糊度的原理固定整周模糊度,进而
成功解算出精准的模糊度。
5 半参数模型中正则化矩阵的选取方法
在半参数模型的计算中,为保证求出半参数模型的解,需要选择合适的正则
化矩阵,在本文的介绍中主要将半参数模型中的信号分为随机量和非随机量2种。
这一分类是为了更好地选择相对应的正则化矩阵进行相关的数据探讨工作。
根据
二者的不同计算方案也有所不同,最终利用这2个方案可以更好地发挥正则化矩
阵方法的研究效果。
6 高精度GPS基线处理中系统误差的分离
新的正则化矩阵的选取方法的提出是为了实现高精度GPS基线处理中系统误
差的分离工作。
在一定条件下,该方案采用的参数变换和选取系统误差方法有了
全新的改变,能够从真正意义上消除或削弱系统误差的影响,真正做到改善GPS
基线解的精密程度,同时做好分析系统误差的分析工作,加速高精度GPS基线处
理中系统误差分析工作的进行。
7 大地测量数据共享环境及其关键技术
在大地测量技术的应用过程中,需要保证信息数据的快速共享,为数据信息
共享环境提供必要的保障措施。
所谓的技术环境主要针对的是影响某种工作开展
的相关技术类型,比如网络技术、数据共享技术以及网络信息技术等。
7.1 开放性集成化技术使用
基于开放性和集成化技术的使用过程中,需要建立起一个稳定的数据共享环境,同时通过开放性与集成化技术的融合运用,有效提高了大地测量技术的工作
效率以及信息收集效率,可以为用户提供更加精确和科学的测量定位参数,以此
来保证后续测量工作的高质量开展。
在我国大地测量工作领域当中应用非常广泛。
在该项工作的开展过程中,需要建立起一个质量更高的大地数据共享环境,并且
对其中所涉及的关键性技术和问题展开分析和研究。
7.2 大地测量数据的建模技术
随着大地测量技术的不断升级和优化,需要在大地测量数据的发展基础之上,有效总结出测量信息的使用方式,并且通过系统软件工作为基础,为大地测量技
术提供必要的平台服务和技术支撑,实现大地测量技术发展的和谐与统一。
在实
际测量工作中,需要基于最初的工作设计、施工方法以及管理工作等技术运用,
本质上是通过多个不同的技术相互之间交织配合工作,所形成的一个综合性过程
网络,要针对各种综合信息数据来开展建模工作,
8 结语
综上所述,大地测量不但可以有效提高地形测图的准确性,同时还可以促进
国防建设以及地震预报等行业的发展。
然而,在进行大地测量的过程中时常会遇
到不适定问题,从而影响测量的准确度以及测量的进行。
通过正则化解法,不但
可以解决大地测量中的不适定问题,同时还可以克服病态性算法并对其进行改进。
此外,通过正则化解法,还可以减弱单频GPS快速定位中的病态性,并使单频GPS快速定位中的ARCE方法进一步改善。
而且,通过正则化解法,可以处理高
精度GPS基线中的系统误差分离问题。
因此,加强对于大地测量中不适定问题正
则化解法的研究是极为必要的。
相信随着对大地测量中不适定问题正则化解法研
究的深入,我国的地质测量技术、国防建设以及地质活动预测技术均会达到一个
新的高度。
参考文献
[1]冯宝宾.不适定问题的正则化解算方法设计及应用[D].成都:成都理工大
学,2010.
[2]王振杰.大地测量中不适定问题的正则化解法研究[J].测绘学报,2005(2):185.
[3]王振杰.大地测量中不适定问题的正则化解法研究[D].武汉:中国科学院研究
生院(测量与地球物理研究所),2003.。