Filters for linear sea-wave prediction

Deterministic sea-wave prediction (DSWP) models are appearing in the literature designed for quiescent interval prediction in marine applications dominated by large swell seas. The approach has focused upon spectral methods which are straightforward and intuitively attractive. However, such methods have the disadvantage that while the sea is aperiodic in nature, the standard discrete spectral processing techniques force an absolutely periodic structure onto the resulting sea surface prediction models. As it is the shape of the sea surface that is important in such applications, particularly near the end of the domain which is important, the standard windowing techniques used in signal processing work to reduce leakage artifacts cannot be employed. This has necessitated the use of end matching methods that can be both inconvenient and may reduce the fraction of the time for which legitimate predictions are available. As a result, an investigation has been undertaken of the use of finite impulse response prediction filters to provide the necessary dispersive phase shifting required in DSWP systems. The present work examines the theoretical basis for such filters and explores their properties together with their application to both long and short crested swell seas. It is shown that wide band forms of such filters are only convergent in the sense of distributions having both infinite duration impulse responses and asymptotically divergent first derivatives. However, appropriate band limitation can produce useful finite impulse responses allowing implementation via standard discrete convolution methods. It is demonstrated that despite the prediction filters having a non-causal impulse response such filters can be used in practice due to a combination of the asymmetric nature of the impulse response and the fundamental nature of the prediction process. The findings are confirmed against actual sea-wave data.     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

惯性／双星组合导航系统的滤波与仿真

介绍了惯性导航系统和双星导航系统各自的运行方式和优缺点；详细推导了组合导航系统的误差方程和量测方程；介绍了Kalman滤波的整个过程并最后进行了仿真计算。仿真结果显示，惯性导航系统和双星导航系统的组合弥补了各自的不足，具有较好的精度和可靠性。     

Parameter identification of inelastic structures under dynamic loads

The generalized model of differential hysteresis contains 13 control parameters with which it can curve‐fit practically any hysteretic trace. Three identification algorithms are developed to estimate the control parameters of hysteresis for different classes of inelastic structures. These algorithms are based upon the simplex, extended Kalman filter, and generalized reduced gradient methods. Novel techniques such as global search and internal constraints are incorporated to facilitate convergence and stability. Effectiveness of the devised algorithms is demonstrated through simulations of two inelastic systems with both pinching and degradation characteristics in their hysteretic traces. Owing to very modest computing requirements, these identification algorithms may become acceptable as a design tool for mapping the hysteretic traces of inelastic structures. Copyright © 2002 John Wiley & Sons, Ltd.     

井间地震串联滤波波场分离方法

本文根据井间地震波场响应以及多种波场分离方法的性能，采用基于Ｋ－Ｌ变换的特征向量滤波与中值滤波相结合的串联滤波方法分离井间地震观测到的上、下行反射波。数值计算结果表明：在波场分离过程中，由于同时综合考虑了地震波能量主分量分布以及反射波同相轴时差两个参数，使得分离效果明显、处理更为简便，具有滤波特征矢量的选取条件更为宽松、无需精确确定拉平轴的时间等特点。     

Simultaneous identification of time and space variant dynamic soil properties during the 1995 Hyogoken-Nanbu earthquake

Time and space variant soil properties at a liquefied site were simultaneously identified in the time domain by using borehole array strong motion records. During soil liquefaction at a site, soils usually show a wide variety of non-linear behavior along the depth as well as non-stationary behavior. Strong ground motion records were obtained at Port Island borehole array observatory, Kobe, during the 1995 Hyogoken-Nanbu earthquake. In this study, the instrumented soil was modeled by the equivalent linear MDOF system, and an extended Kalman filter with local iteration was employed for the identification of the soils. The identification process was successfully conducted, and the stress–strain relationships of the soils at the liquefied site were obtained from different depths all at once.     

锚固站对BDS自主导航精度的影响分析

介绍了自主导航的轨道确定及时间同步观测方程。以北斗仿真全星座为对象,通过采用仿真星间及卫星与地面锚固站间观测值,进行了60 d自主导航解算,分别探讨了锚固站数量及锚固站观测连续性对北斗卫星导航系统（BDS）3类卫星自主导航精度的影响。结果表明：锚固站数量及观测连续性对RERR及CERR无影响;加入1个锚固站即可显著改进URE结果精度,继续增加锚固站数量虽然可进一步提高URE精度但其改进效果有限;锚固站观测中断时间越长,其对应自主导航精度越低。因此,在BDS自主导航运行模式下应保持较高的锚固站观测频次以保证自主导航精度;另外,锚固站数量及观测连续性对北斗系统3类卫星自主导航精度的影响无显著差异。     

集合卡尔曼滤波同化中雷达位置的敏感性研究

针对对流尺度集合卡尔曼滤波（EnKF）雷达资料同化中雷达位置对同化的影响进行研究。为了考察强对流出现在雷达不同方位时集合卡尔曼滤波同化雷达资料的能力,以一个理想风暴为例,设计了8个均匀分布在模拟区域周围的模拟雷达进行试验。单雷达同化试验中,初期同化对雷达位置较敏感,而十几个循环后对雷达方位的敏感性降低。造成初期同化效果较差的雷达观测位于模拟区域正南和正北方向,这两部雷达与模拟区域中心的连线垂直于风暴移动方向（即环境气流的方向）。双雷达试验的结果表明,正东、正南、正西和正北方向的雷达组合观测会使同化初期误差较大,这说明并不是所有与风暴连线成90°的雷达组合都能在短时同化中得到合理的分析结果,还需要都处于模拟区域对角线上（即与环境气流成45°夹角）,同化效果才较好。短时同化后的确定性预报结果表明,较大分析误差也会导致较大预报误差。这些分析误差主要是由于同化初期不准确的集合平均场驱动出的不合理的背景误差协方差造成的。当背景场随着同化循环得到改进后,驱动出的合理的背景误差协方差使得不同位置雷达同化造成的差异逐步减小。基于上述结果,引入迭代集合均方根滤波（iEnSRF）算法,结果显示使用该算法后,雷达位置对同化效果的影响减小,同化不同位置的雷达资料均能有效降低分析和预报误差。      

实现遥感相机自主辨云的小波SCM算法

云噪声是光学卫星影像的常见问题,为了衡量云噪声对影像融合带来的影响,本文以高通滤波融合算法为例进行分析,指出云与地物的均值相差越大,云对影像融合的影响越大,并提出了一种针对含云影像的融合方法,即联合云检测与高通滤波的含云影像融合方法。该方法首先利用NIR/R-OTSU云检测算法实时进行云检测,判别出影像中的云覆盖区域;其次采用局部优化策略利用高通滤波融合方法分块对非云区域进行处理,得到融合影像。利用资源三号多光谱和正视全色影像进行融合实验,结果表明,本文算法比高通滤波融合方法、亮度色度饱和度（intensity hue saturation,IHS）变换融合方法、Pansharp融合方法更适用于含云影像的融合处理。     

基于DEnKF方法的考虑次网格变异性的MODIS雪盖同化

基于通用陆面模型(CoLM)和确定性集合卡尔曼滤波算法发展了一个考虑模型次网格变异性的MODIS雪盖同化方案,提高雪深模拟的估计精度。利用北疆阿勒泰地区5个气象站点2007年11月至2008年4月逐日雪深观测数据对同化结果进行了验证。结果表明,该同化方案不需要对MODIS雪盖观测数据进行扰动,能明显提高雪深模拟的精度。另外,雪深同化结果与地面观测雪深具有一致的时间变化趋势,能准确地反映积雪深度在各个不同时段的变化特性。