区间约束平差模型的共轭梯度积极集算法

主要研究参数带有区间约束的平差算法,通过把平差问题转化成一个带有区间约束的二次规划问题,利用积极集对二次规划问题进行划分与重组,结合无约束共轭梯度优化算法,给出了带有区间约束的平差算法,并同时给出了参数解的精度评估。由于投影梯度法可以迅速改变积极约束集的构成,新的算法比经典的积极集法效率更高,可以降低模型的不适定性,保持参数先验信息中的统计、几何或物理意义,适合于求解大规模的带有区间约束的平差问题。     

一种附加边界约束和内约束的多中心ERP融合模型

不同技术、不同分析中心得到的地球自转参数（Earth rotation parameters,ERP）往往是不同的,为提供统一的ERP供用户使用,常需对ERP进行融合处理。提出了一种基于多分析中心ERP结果的附加边界约束和内约束融合模型,即先通过参数变换把各分析中心结果转换到相同时刻,考虑到相邻观测时段边界点处ERP应当一致这一特点,施加边界约束,然后对各分析中心的长期解施加转换参数内约束,最后得到多分析中心ERP的融合解。采用从2005—2011年共6 a的7个全球卫星导航系统（Global Navigation Satellite System,GNSS）分析中心的结果进行融合处理,并与IERS C04（International Earth Rotation and Reference Systems ServiceCombined 04）结果进行比较。结果表明,所提出的融合方法计算结果的精度有明显改善。     

带不等式约束的最小二乘方法比较

解决带不等式约束的平差问题已有较多方法,但不同方法之间的对比分析不足。本文讨论三种带不等式约束的最小二乘算法:简单迭代算法、虚拟误差方程法、迭代乘子法,三种算法都是基于有效约束理论将不等式约束转换为等式约束,通过迭代求解。采用同一组数据和蒙特卡罗仿真实验方法对三种算法计算结果进行对比分析,实验结果分析表明:迭代乘子法算法效率最高,推荐使用该算法。     

改进的约束总体最小二乘法

针对传统的约束最小二乘模型和总体最小二乘模型的局限性,该文提出了一种改进的约束总体最小二乘法。假设约束总体最小二乘问题中约束方程系数矩阵也存在误差,然后构造函数模型的广义拉格朗日函数,采用最小二乘法迭代求解非线性的法方程,最终获得了改进的约束总体最小二乘法的牛顿-高斯迭代公式和平差模型精度的无偏估计。该算法采用了更接近实际的平差模型,能够获得更加接近真值的估计参数,同时平差模型的精度更加接近模拟数据加入的噪声水平。实验结果表明,本文算法可有效解决对参数进行约束时的数据处理问题。     

一种三维激光点云中建筑物立面渐进分割方法

针对建筑物立面分割的问题,该文提出了一种三维激光点云中建筑物立面渐进分割方法。建立三维格网索引,通过分析建筑物立面在三维格网内的空间分布特征和二维平面格网内投影的线性分布特征,确定立面种子格网和投影线的拟合点,拟合投影线并基于种子格网约束生长完成每层格网中立面粗分割;使用RANSAC算法对粗分割后的立面点云进行面拟合,实现精细分割,并将各层格网中的立面进行合并,实现建筑物立面的完整分割。实验结果表明,该方法能有效地实现建筑物立面的精细化分割,有助于后续的建筑物精细三维模型重建。     

融合轨道轨迹约束的动态单历元整周模糊度解算方法对比分析

GPS动态单历元定位中,高精度定位的关键是整周模糊度的正确解算。针对这个问题,本文利用载波相位双差模型、载波相位联合伪距、宽巷载波联合宽巷伪距、融合轨道轨迹约束条件4种不同的方法解算分析模糊度的浮点解,并用LAMBDA法对模糊度浮点解进行固定,然后以LGO软件的解算结果作为真值,对比分析几种不同方法的实验结果。实验结果表明,解算模糊度的方法不同,模糊度浮点解的精度、分解速度和成功率均不同,宽巷组合的模糊度解算成功率整体高于载波相位联合伪距的解算成功率,而在宽巷组合的基础上增加轨迹作为约束条件后,其模糊度解算成功率进一步得到提高。融合轨道轨迹约束的模糊度分解方法得益于信息量的增加,在一定程度上能够加快模糊度分解的速度。     

等式约束病态模型的截断奇异值解及其统计性质

利用平差模型间合理的先验信息能够显著提高解的稳定性和精度。本文基于病态模型引入等式约束条件,并采用截断奇异值法重构了系数阵以削弱其病态性;建立了修正等式约束模型,导出了病态模型的约束截断奇异值解及其偏差、方差以及均方误差公式;分析了截掉奇异值所引起解的偏差引入量与方差下降量的关系,得到了确定截断参数的条件。数值算例和病态测边网算例分析结果表明,最小二乘解严重偏离真值,500次模拟实验的平均RMSE为6.693 5,正则化解和截断奇异值解精度较最小二乘解有所提高,平均RMSE分别为0.365 8和0.365 2;本文提出的约束截断奇异值解的精度最高,与约束正则化解精度相当,其平均RMSE仅为0.057 3。     

利用线性判别分析的LDA-SIFT影像匹配算法

针对SIFT算法特征向量维数较高,匹配实时性较差的问题,提出一种利用线性判别分析的LDA-SIFT算法。该算法首先利用SIFT算法提取特征点并生成特征向量矩阵;然后将特征向量矩阵转换为种子点向量矩阵并为数据设置标签;接着利用线性判别分析对种子点向量进行降维;最后在低维特征空间将种子点向量矩阵转换为特征向量矩阵并在欧式空间匹配,运用RANSAC算法剔除误匹配。采用多组图像验证该算法匹配的性能,实验结果表明,LDA-SIFT算法能够有效降低特征向量维度,缩短特征匹配时间,匹配精度与降维前相近。     

存在有限导流断层的复合油藏试井模型及井底压力动态分析

运用数学物理方法、渗流力学理论,考虑流体通过断层面和断层内部流体的流动,建立两区无限大线性复合油藏中存在有限导流断层的试井模型,绘制井底压力响应特征曲线并分析影响因素,并验证模型的正确性与实用性。结果表明:无因次断层导流能力越大,压力导数曲线下凹越多;界面表皮因数越大,压力导数曲线上翘越多,当界面表皮因数足够大时,表现为封闭断层的特征;无因次断层导压系数主要影响压力导数曲线下降幅度;流度比、厚度比、导压系数比主要影响压力导数曲线上升和下降幅度。该结果为存在类似边界条件的复合油藏试井解释提供指导。     

大别山及其周边现今垂直形变研究

基于1980年、2005年和2015年3期区域精密水准网观测资料,利用线性动态平差模型计算获取大别山地区水准网长期垂直运动速度场图像。研究发现,淮河平原地表下沉较为严重,大别山呈现弱隆升趋势运动,长江谷地边缘地区较大别山区呈现明显的隆升运动。跨郯庐断裂带水准剖面结果显示,垂直运动与地形呈负相关和弱相关。