基于增强型GPS的自适应UKF实时星间相对定位方法

高精度的星间实时相对定位是卫星编队飞行的一项关键技术。本文以双星编队为例,提出一种基于GPS双频P码、双频载波相位以及星间距离观测信息的自适应UKF、实时相对定位方法。仿真结果表明:相比传统的EKF方法,该算法能有效地提高滤波稳定性及相对定位精度,而且,星间测距信息的引入能大大减少整周模糊度判定所需的历元数。     

基于ICCP算法的重力辅助惯性导航

迭代最近等值线算法(ICCP)是一种重要的匹配导航算法,文中首先介绍ICCP算法的基本原理,随后在0.2′×0.2′重力异常数据库的基础上,利用ICCP算法进行仿真计算得到最佳匹配位置。最后为了验证匹配位置是否可用于修正惯导误差,提出将匹配位置误差作为观测量,用卡尔曼滤波对惯导系统误差进行最优估计。由最后的仿真结果可以看出,ICCP算法可有效抑制惯导纬度误差的增长,且最大纬度误差不超过2,′以匹配位置误差作为观测量可以用来估计惯导方位误差角。     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

Managing sustainable farmed landscape through 'alternative' food networks: a case study from Italy

This paper focuses on a case study of an 'alternative' food network based in the Abruzzo National Park, Italy, to explore how ideas of sustainable farmland management can be expressed through broader understandings of developing networks of care concerned with local economies and societies, high-quality specialist food products, particular 'traditional' farming practices and livestock breeds, as well as the ecology of a farmed landscape. The scheme allows customers, internationally as well as in Italy, to 'adopt' a milking sheep on a large mountain farm. In return, adopters are sent food products from the farm. The adoption scheme is inter-twined with an agri-tourism project which provides accommodation, runs a restaurant and engages in educational activities. The scheme is the result of the individual initiative of its founder, and is associated with a strongly expressed ethical position concerning the value of sustaining valued local rural landscapes and lifestyles, and the importance of 'reconnecting' urban dwellers with rural areas, farming and 'quality' food production. Yet the localness of the scheme is sustained through wider national and international networks: volunteer and paid workers are drawn from several European countries, funding has been acquired from the EU LEADER programme, and internet and transport technologies are essential in connecting with and supplying an international customer base. The broader economy of care instanced in this case study draws attention to a need to develop strategies for sustainable farmland management constructed around wider programmes of social, economic and cultural, as well as environmental, concern.     

基于卡尔曼滤波的时间序列影像背景估计

文章将卡尔曼滤波理论应用于时间序列影像的背景估计,引入衰减记忆因子解决滤波的发散问题,最终实现了对背景信号的估计。     

一种新的综合Kalman滤波及其在星载GPS 低轨卫星定轨中的应用

针对将Kalman滤波方法应用到星载GPS定轨时，由于动态噪声和观测噪声确定不准而造成滤波的发散、污染观测值造成Kalman滤波估值的扭曲及计算舍入误差可能带来协方差阵的不正定性等缺陷，提出了一种新的综合Kalman滤波方法。该方法用拟准检定法准确地探测和修正量测方程中存在的粗差；用UD分解算法克服了数值的不稳定性，改进了计算精度；用Sage自适应滤波器克服滤波器的发散。算例结果表明，这种综合卡尔曼滤波方法具有数值稳定性好、较强的自适应性和较好地削弱粗差影响等优点。     

MRI分割的几何活动轮廓模型的半隐式AOS迭代方法

几何活动轮廓模型以其计算稳定等优点，被广泛用于图像分割。给出一种基于AOS(additive operator splitting)算法的快速活动轮廓模型，将AOS方法与几何活动轮廓模型相结合构造一种半隐式迭代模型，它对迭代步长没有限制，从而可以选择较大的迭代步长，提高了模型的时间效率，并使用区域信息构造系数矩阵，使模型具有全局性。对MR(magnetic resonance)图像分割的实验表明，该模型可以得到较好的分割结果，同时时间效率也有大幅提高，方便了实时应用。     

伴随模式同化系统在修正模式地形中的应用

任何模式都只是实际大气的一种近似模拟，存在误差。鉴于传统的四维资料伴随模式同化系统都是假设模式完全正确仅对初始场进行修正。将伴随模式同化系统应用于修正模式地形误差，通过对不同初始地形的修正试验表明，MM5伴随模式同化系统能很好地对地形误差进行修正，能够反演出与初始气象要素场分辨率相匹配、与模式更协调的地形场，为伴随模式同化系统的广泛应用提供了一种新的思路。     

变分同化方法在Lorenz系统中的简单应用研究

利用Lorenz模式作变分同化数值试验，通过对一个简单系统的讨论，介绍四维变分同化方法。对初值敏感性和观测点的个数及观测值作了对比试验，发现随着模式对初值敏感性的增加，同化效果会越来越差；观测点越少，观测值误差越大，这些都会影响同化效果，甚至导致同化失败。