稳健最优不变二次无偏估计

导出了稳健最优不变二次无偏估计、稳健最小范数二次无偏估计、稳健Ｈｅｌｍｅｒｔ估计，并说明了最优不变二次无偏估计、最小范数二次无偏估计以及Ｈｅｌｍｅｒｔ估计等均是稳健最优不变二次无偏估计的特例。     

A distribution-free alternative to least-squares regression and its application to Rb/Sr isochron calculations

A distribution-free estimator of the slope of a regression line is introduced. This estimator is designated S_m and is given by the median of the set of n(n – 1)/2 slope estimators, which may be calculated by inserting pairs of points (X _i, Y_i)and (X _j, Y_j)into the slope formula S _i = (Y _i – Y_j)/(X _i – X_j),1 i < j n Once S _m is determined, outliers may be detected by calculating the residuals given by R_i = Y_i – S_mX_i where 1 i n, and chosing the median R_m. Outliers are defined as points for which |R_i – R_m| > k (median {|R _i – R_m|}). If no outliers are found, the Y-intercept is given by R_m. Confidence limits on R_m and S_m can be found from the sets of R_i and S_i, respectively. The distribution-free estimators are compared with the least-squares estimators now in use by utilizing published data. Differences between the least-squares and distribution-free estimates are discussed, as are the drawbacks of the distribution-free techniques.     

神经网络辅助的GPS/INS组合导航自适应滤波算法

首先利用预报残差构造的最优自适应因子设计GPS/INS组合导航自适应滤波器。并针对BP神经网络存在的训练速度慢、容易陷入局部极小等问题,给出网络的改进算法。利用神经网络对自适应滤波器状态方程的预报值进行在线修正,给出神经网络辅助的GPS/INS组合导航自适应滤波算法。最后,利用实测数据进行验证。结果表明,改进的神经网络算法明显提高网络收敛速度;两种自适应滤波算法相对标准组合导航算法都能够可靠地反映载体运动轨迹;神经网络辅助的GPS/INS组合导航自适应滤波算法相对GPS/INS组合导航自适应滤波算法在精度和可靠性方面又有明显提高。     

有偏迭代估计在短基线集形变模型中的应用

进行SBAS技术时,选择的干涉对一般会出现一景影像和其他几景影像干涉的现象,导致构建的线性形变模型病态。在这种情况下,最小二乘法求解SBAS线性形变模型已不太适合。针对这一问题,本文结合Liu估计的有偏迭代法在求解病态线性模型时能够克服病态性对结果的影响,提出将此方法应用到求解SBAS线性形变模型中。以覆盖济宁地区的13景ENVISAT ASAR数据展开了SBAS技术在城市地面沉降监测中的应用实验,分别利用LS法、基于Liu估计的有偏迭代法解算形变模型。结果表明,基于Liu估计的有偏迭代法获得研究区的年平均沉降速率标准差比LS法获得的低1.4,且其求解的结果相较于LS法获得的结果更准确、稳定。     

瞬时属性估计器在探地雷达数据处理中的应用

为了提高浅层埋藏目标识别的瞬时物理属性的精度,根据极平坦频率特性有限脉冲相应滤波器的思想,首次提出用它作为微分器和延时器组成解析信号瞬时属性估计器（FIFM）,对探地雷达信号的瞬时属性进行提取。结果表明,瞬时属性剖面的分辨率得到了明显的提高;采用瞬时属性估计器进行探地雷达信号瞬时参数的研究是一种非常有效的方法,可用于公路、桥梁、隧道等结构物的无损检测数据分析。      

Stochastic inversion of pneumatic cross-hole tests and barometric pressure fluctuations in heterogeneous unsaturated formations

The distributions of permeability and porosity are key factors that control airflow and gas phase transport in unsaturated formations. To understand the behavior of flow and transport in such formations, characterization procedure is a typical approach that has been widely applied to laboratories and fields. As is recognized by most investigations, this approach relies on accurate measurements, and more importantly, an adequate tool to interpret those measurements from experiments. This study presents a pneumatic inverse model that is capable to estimate the distributions of permeability (k) and porosity () with high resolution in heterogeneous unsaturated formations. Based on the concept of sequential successive linear estimator (SSLE), the developed model accounts for compressibility and density of air and estimates the geologic parameters using air pressure measurements from sequential cross-hole pneumatic pumping or injection tests. Four synthetic examples, including a one-dimensional well-posed, a horizontally two-dimensional, and two three-dimensional problems, are used to evaluate the developed model in estimating the distributions of permeability and porosity in unsaturated formations. Results of the numerical experiments are promising. The developed pneumatic inverse model can reconstruct the property (i.e., permeability and porosity) fields if the well-defined conditions are met. With a relatively small number of available measurements, the proposed model can accurately capture the patterns and the magnitudes of estimated properties for unsaturated formations. Results of two complex three-dimensional examples show that the proposed model can map the fracture connectivity using a small number of subsurface pressure measurements and estimate k and in shallow soil layers using spatial variations of barometric pressure.     

Spatial interpolation of marine environment data using P-MSN

ABSTRACT

When a marine study area is large, the environmental variables often present spatially stratified non-homogeneity, violating the spatial second-order stationary assumption. The stratified non-homogeneous surface can be divided into several stationary strata with different means or variances, but still with close relationships between neighboring strata. To give the best linear-unbiased estimator for those environmental variables, an interpolated version of the mean of the surface with stratified non-homogeneity (MSN) method called point mean of the surface with stratified non-homogeneity (P-MSN) was derived. P-MSN distinguishes the spatial mean and variogram in different strata and borrows information from neighboring strata to improve the interpolation precision near the strata boundary. This paper also introduces the implementation of this method, and its performance is demonstrated in two case studies, one using ocean color remote sensing data, and the other using marine environment monitoring data. The predictions of P-MSN were compared with ordinary kriging, stratified kriging, kriging with an external drift, and empirical Bayesian kriging, the most frequently used methods that can handle some extent of spatial non-homogeneity. The results illustrated that for spatially stratified non-homogeneous environmental variables, P-MSN outperforms other methods by simultaneously improving interpolation precision and avoiding artificially abrupt changes along the strata boundaries.     

基于ARMA模型的BDS卫星钟差异常值探测及其短期预报

星载原子钟在运行过程中会受到恶劣空间环境与设备老化等因素的影响,使得卫星钟差数据中经常存在异常值,其中AO(additive outlier)类异常值是钟差序列中常见的一类异常值.结合最大期望算法与自回归滑动平均(autoregressive moving average, ARMA)模型,提出一种AO类异常值探测算法...     

Sequential kriging and cokriging: Two powerful geostatistical approaches

A sequential linear estimator is developed in this study to progressively incorporate new or different spatial data sets into the estimation. It begins with a classical linear estimator (i.e., kriging or cokriging) to estimate means conditioned to a given observed data set. When an additional data set becomes available, the sequential estimator improves the previous estimate by using linearly weighted sums of differences between the new data set and previous estimates at sample locations. Like the classical linear estimator, the weights used in the sequential linear estimator are derived from a system of equations that contains covariances and cross-covariances between sample locations and the location where the estimate is to be made. However, the covariances and cross-covariances are conditioned upon the previous data sets. The sequential estimator is shown to produce the best, unbiased linear estimate, and to provide the same estimates and variances as classic simple kriging or cokriging with the simultaneous use of the entire data set. However, by using data sets sequentially, this new algorithm alleviates numerical difficulties associated with the classical kriging or cokriging techniques when a large amount of data are used. It also provides a new way to incorporate additional information into a previous estimation.