基于条件平差的卫星双向时间传递链路性能优化方法

卫星双向时间传递(Two-Way Satellite Time and Frequency Transfer, TWSTFT)是目前精度最高的时间传递方法之一,同时也是参与国际原子时计算的守时实验室之间比较原子时尺度的一种主要方法.提高TWSTFT链路的短期稳定度,降低周日效应对链路时间传递结果的影响,对优化TAI (International Atomic Time)的性能具有现实意义.提出了一种基于条件平差的TWSTFT链路性能优化方法,先依据TWSTFT链路测量噪声水平与谱分析结果建立TWSTFT链路性能优化网络(简称优化网络),再根据优化网络中各链路测量噪声分析结果设置权系数阵,建立条件平差模型.选取亚太地区的中国计量科学研究院(National Institute of Metrology, NIM)-中国科学院国家授时中心(National Time Service Center, NTSC)卫星双向时间传递链路作为待优化链路,以NTSC、NIM以及德国联邦物理技术研究所(Physikalisch-Technische Bundesanstalt,PTB)之间的TWSTFT链路组成优化网络,对优化网络的组网方法和条件平差模型进行实验验证.结果表明,平差后待优化链路短期稳定度得到了改善,同时其受周日效应的影响降低了约24.6%.使用该方法,能够有效提高待优化链路的时间传递性能.     

Non-stationary covariance function modelling in 2D least-squares collocation

Standard least-squares collocation (LSC) assumes 2D stationarity and 3D isotropy, and relies on a covariance function to account for spatial dependence in the observed data. However, the assumption that the spatial dependence is constant throughout the region of interest may sometimes be violated. Assuming a stationary covariance structure can result in over-smoothing of, e.g., the gravity field in mountains and under-smoothing in great plains. We introduce the kernel convolution method from spatial statistics for non-stationary covariance structures, and demonstrate its advantage for dealing with non-stationarity in geodetic data. We then compared stationary and non- stationary covariance functions in 2D LSC to the empirical example of gravity anomaly interpolation near the Darling Fault, Western Australia, where the field is anisotropic and non-stationary. The results with non-stationary covariance functions are better than standard LSC in terms of formal errors and cross-validation against data not used in the interpolation, demonstrating that the use of non-stationary covariance functions can improve upon standard (stationary) LSC.     

概括平差模型的解法探讨

分析了概括平差模型系数矩阵可能出现非列满秩情况的原因,研究了如何使系数矩阵为列满秩矩阵,避免了系数矩阵为非列满秩时的复杂解算过程。     

全极化SAR数据反演桥面高度

根据高分辨率SAR图像上建筑区的影像特征, 提出了基于灰度共生矩阵(gray-level cooccurrence Matrix, GLCM)纹理分析的建筑区提取方法, 该方法由初步定位和边界调整2个步骤组成, 均遵循特征计算、基于Bhattacharyya距离的特征选择和KNN分类流程, 所不同的是2个步骤中分别采用了逐块和逐点计算纹理特征的方式以兼顾纹理分析的效率和准确性。文中对不同SAR传感器获取的图像进行了实验。实验结果表明, 选用具有最大Bhattacharyya距离值的3或4个特征可以获得较好的初步定位结果, 建筑区的检测率超过80%, 虚警率低于10%;随着边界调整的进行, 检测到的建筑区边界逐渐接近于真实边界。实验结果验证了该算法的有效性。     

半参数估计与LS估计的有效性比较

建立半参数回归模型,将系统误差加入,推导出了模型正则化矩阵时参数平差的计算方法,并用直接法得到了参数、非参数的估计量,在MSE准则下与最小二乘估值(简称LS估值)进行比较,证明了半参数估值的有效性。     

网络分层用于最短路问题的算法研究

提出了一种基于Dijkstra方法的网络分层算法,实现了两点间节点数最少条件下最短通路的求取,并与传统Dijkstra算法进行了比较,得到了一些有益的结论。     

三维三分量CSAMT法有限元正演模拟研究初探

首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证.     

地震波的场方程矩阵和能量的正定二次型及其意义

场方程是能够表述半空间地震波场的整体特征.波动方程、速度方程和能量方程.通过分析可知每个场方程都具有各自的“场方程矩阵”.能量方程能够对所有场方程矩阵进行综合和贯通,给出了能量方程以“弹性矩阵”为核心的普适性表达形式.最后,运用矩阵的正定二次型理论阐述了“能量矩阵与弹性矩阵”之间一致的对称性和正定性.能量矩阵蕴含的动态力的平衡关系、速度的时间_空间分布和能量的传播及变化的物理意义,能够从能量矩阵的正定二次型特性表述出来.本文研究分析问题的方法完全适用于复杂介质模型,相关的认识和结论可以拓展到均匀黏弹性各向同性介质、均匀弹性各向异性介质、均匀黏弹性各向异性介质以及比奥饱和流体介质.     

盘锦湿地芦苇生态系统长期通量观测研究

针对2004年5月26日-2005年10月15日盘锦湿地芦苇生态系统碳通量、感热通量和潜热通量资料进行分析。结果表明:芦苇湿地具有较强的碳汇作用;2005年芦苇湿地固定二氧化碳为13.32 t/hm2,日平均感热通量和潜热通量分别为2 464 kJ/m2和3 880 kJ/m2。2004年和2005年6～9月的日累积值波文比平均值均为0.15。芦苇湿地碳通量、感热通量和潜热通量的日动态呈单峰单谷型变化,极值出现在中午前后,夜间线形平直。芦苇生长季白天通量绝对值远较夜间大,白天碳吸收,夜间碳排放。CO2浓度年平均日变化曲线亦为单谷单峰型,夜间浓度较高且逐渐升高,直到日出前达到峰值;日出后急剧下降,傍晚又开始逐渐增加。芦苇湿地感热通量昼正夜负,潜热通量与林地不同,全天为正。各通量季节变化明显,冬季CO2通量日变化不明显,趋近于零;感热通量总体向上输送,春季数值较大,生长季数值较小;潜热通量冬季最小,接近0,春季开始显著增加,生长季达到最大。     

Macroscopic criteria for Green type porous materials with spheroidal voids: application to double porous materials

Combined effects of matrix plastic compressibility and void shape are investigated for ductile porous materials. To this end, a spheroidal volume containing a confocal spheroidal (prolate or oblate) void subjected to uniform strain rate boundary conditions has been first studied. A Green type matrix is chosen as a prototype for investigating effects of plastic compressibility. This is carried out by using a kinematics limit analysis theory from which a closed‐form expression of the macroscopic criterion is established for the considered class of materials. These results are then extended to ductile porous materials made up of a green matrix containing randomly oriented spheroidal voids. In the framework of a two‐step homogenization procedure, the obtained results are implemented to describe the macroscopic behavior of double porous materials involving spherical voids at the microscale and randomly oriented and distributed spheroidal voids at the mesoscale. For validation purpose, the new derived criteria are assessed and validated by comparing their predictions to available upper bounds and numerical data from literature. Copyright © 2017 John Wiley & Sons, Ltd.