非线性参数估计的自适应松弛正则化算法

非线性方程参数估计存在的弊端在于非线性观测方程存在不适定问题时，以线性化平差估计和高斯牛顿为代表的经典数值算法会产生较强的不稳定特征。因此，针对传统非线性最小二乘求解不稳定且可靠性低的特点，基于稳定泛函极小准则最优化思想，提出了一种自适应松弛正则化数值算法。该算法采用正则化参数几何递增计算方法和残差最小步长准则，实现了正则参数和迭代步长计算的完全自适应，提高了非线性迭代收敛效率。以病态仿真数据和水下实测数据为例，验证了该方法的数值收敛解优于线性平差估计解，收敛效率优于迭代Tikhonov正则化方法。     

GNSS接收机DCB短期时变特征分析和估计方法

提出一种使用非差非组合精密单点定位(PPP)估计和分析接收机DCB短时时变特征的方法。首先利用非差非组合PPP得到包含接收机DCB的重构电离层参数估值;然后通过IGS电离层GIMs格网模型内插剥离各历元站星斜向电离层距离延迟;最后通过最小二乘约束得到各历元接收机DCB解。由于格网本身精度(2～8 TECU)和插值精度限制,解算出来的接收机DCB并不能真实反映其短期时变特征。为此,提出利用站间单差或者历元间差分的方法还原其真实的变化态势。实验结果表明,所提出的方法能够正确估计接收机DCB,并能真实还原其短期时变特征,具有良好的适用性。     

基于ICSO的DGPS整周模糊度的求解方法

针对差分全球定位系统(DGPS)模糊度解算过程中效率低,搜索慢的问题,对鸡群优化算法(CSO)进行适应性改进,并将改进后的鸡群优化算法(ICSO)应用到整周模糊度的快速解算中,利用卡尔曼滤波求出双差模糊度的浮点解和协方差矩阵,采用Lenstra-Lenstra-Lovasz (LLL)降相关算法对模糊度的浮点解和方差协方差矩阵进行降相关处理,以降低模糊度各分量之间的相关性,在基线长度固定的情况下,利用ICSO搜索整周模糊度的最优解. 采用经典算例进行仿真,仿真结果表明,与已有文献相比在整周模糊度的解算过程中改进的鸡群优化算法能有效提高搜索速度和求解成功率.      

Helmert方差分量估计在GPS/GLONASS/BDS组合定位权比确定中的应用

因GNSS系统间观测噪声、轨道精度的差异，采用经验权比进行组合定位难以得到最优结果。基于此，在GPS/GLONASS/BDS组合定位中引入Helmert方差分量估计，对GPS/GLONASS/BDS组合单点定位和基线解算中各系统观测值进行合理定权。实验表明，采用该方法确定的伪距观测值最佳权比为5∶1∶1，相位观测值最佳权比为1∶1∶1，有效提高了GPS/GLONASS/BDS组合定位的精度和可靠性。     

单基站基于波达角度的刚体位置姿态最大似然估计

针对刚体定位同时估计目标的位置和姿态信息的问题，研究了在三维空间中利用单基站进行刚体定位（RBL）的框架.该框架采用单个基站对安装在刚性目标表面的小规模无线传感器信号的波达方向进行测量，并将其与传感器拓扑信息融合，提出了两种用于RBL的最大似然估计算子.利用改进的高斯-牛顿算法对旋转矩阵和平移向量的最大似然估计算子进行优化估计，解算出了目标物体的三维位置和姿态.仿真结果表明，文中提出的最大似然算子可以接近理论克拉美罗下限，并且在收敛成功率和运算成本方面具有较为出色的性能.      

基于SAR干涉点目标分析技术的城市地表形变监测

通过深入研究干涉点目标的相位模型,提出基于空间搜索的邻近点目标干涉相位差解缠方法,用以计算点目标的地形残差和线性形变,以及分离点目标大气延迟相位和非线性形变相位的时空域滤波方法,解决干涉点目标分析中的关键问题.最后,以苏州地区地表沉降监测为应用试验,利用形成的SAR干涉点目标形变信息提取技术,获取苏州市区1992-2002年间的地表沉降信息.研究结果与已有文献记录保持比较好的一致性,证明SAR干涉点目标技术完全可以发展成为应用于城市地表形变监测的实用化技术.

Abstract：

Interferometric point target analysis (IPTA) is one of the latest developments in radar interferometric processing. It is achieved by analysis of the interferometric phases of some individual point targets, which are discrete and present temporarily stable backscattering characteristics, in long temporal series of interferometric SAR images. This paper analyzes the interferometric phase model of point targets, and then addresses two key issues within IPTA process. Firstly, a spatial searching method is proposed to unwrap the interferometric phase difference between two neighboring point targets. The height residual error and linear deformation rate of each point target can then be calculated, when a global reference point with known height correction and deformation history is chosen. Secondly, a spatial-temporal filtering scheme is proposed to further separate the atmosphere phase and nonlinear deformation phase from the residual interferometric phase. Finally, an experiment of the developed IPTA methodology is conducted over Suzhou urban area. Totally 38 ERS-1/2 SAR scenes are analyzed, and the deformation information over 3 546 point targets in the time span of 1992-2002 are generated. The IPTA-derivecl deformation shows very good agreement with the published result, which demonstrates that the IPTA technique can be developed into an operational tool to map the ground subsidence over urban area.     

STUDY ON THE INTERACTION BETWEEN NDVI PROFILE AND THE GROWING STATUS OF CROPS

Daily and ten-day Normalized Difference Vegetation Index(NDVI) of crops were retrieved from meteorological statellite NOAA AVHRR images ,The temporal variations of the NDVI were analyzed during the whole growing season,and thus the principle of the interaction between NDIV profile and the growing status of crops was discussed,As a case in point,the relationship between integral NDVI and winter wheat yield of Henan Province in 1999 had been analyzed.By putting integral NDVI values of 60 sample counties into the winter wheat yield-integral NDVI coordination,scattering map was plotted. It demonstrated that integral NDVI had a close relation with winter wheat yield.These relation could be described with linear,cubic polynomial ,and exponential regression,and the cubic polynomial regression was the best way,In general ,NDVI reflects growing status of green vegetation ,so crop monitoring and crop yield estimation could be realized by using remote sensing technique on the basis of time serial NDVI data together with agriculture calendars.     

中国大陆科学钻探先导孔零偏VSP资料解释

用六级三分量检波器在中国大陆科学钻探先导孔中实施了零偏VSP测量.数据处理结果表明,中国大陆科学钻探孔区超高压变质岩石的地震波速度主要介于4500～7000m/s之间,显著高于一般的沉积岩地区,而且随深度变化不明显.声波测井速度系统地稍低于VSP层速度,可能是由于井壁处岩石的完整性受到破坏而造成的.地震波速度与岩石密度和岩性存在明确的对应关系,榴辉岩的密度和地震波速度均显著高于片麻岩类岩石;由榴辉岩退变生成的斜长角闪岩类岩石,其密度和地震波速度均呈现出较大的变化,主要与其退变质程度有关;超基性岩中的裂隙系统导致其密度和地震波速度大幅度下降.由于榴辉岩与其他岩石类型之间存在较大的波阻抗差异,因此用零偏VSP资料标定该区地震波的地质层位是有效的.关于地震波反射的原因,通过综合研究地震波(包括反射纵波、上行转换横波、井筒波)的特征、岩石速度和密度分布以及井径变化,认为主要是岩性分界面、韧性剪切带和断裂(带),但还有一些因素尚待进一步研究.     

抗差估计的初值选择

对L1 范估计的平差值和LS估计的平差值作为抗差估计的初值进行了比较 ,指出在 ρ函数是严凸函数时 ,二者无明显的区别。因为LS估计简单且易于计算 ,因此建议用LS估计的平差值作为抗差估计的初值。