基于均值漂移的粒子滤波算法设计及其在导航数据处理中的应用

针对标准粒子滤波算法中存在的计算量大和粒子的权值退化的缺陷,将均值漂移算法和PF算法进行融合,设计基于均值漂移搜索算法的粒子滤波新算法。该算法仍遵从粒子滤波算法的计算框架,基本原理是利用MS算法对粒子的聚类作用,将均值漂移思想融合到粒子滤波算法的重要性采样过程中,对粒子集进行确定性搜索,使每个粒子收敛于局部最优值,这样粒子的状态表示更接近真实的状态分布,因此只需较少的粒子数便可达到未嵌入MS的使用大量粒子数的粒子滤波状态估计的性能,从而在缓解粒子的权值退化的同时提高粒子滤波算法的实时性。大量的数值试验和对GPS/DR组合导航数据处理的结果验证了该方法的有效性。     

Cassini SAR, radiometry, scatterometry and altimetry observations of Titan’s dune fields

Large expanses of linear dunes cover Titan’s equatorial regions. As the Cassini mission continues, more dune fields are becoming unveiled and examined by the microwave radar in all its modes of operation (SAR, radiometry, scatterometry, altimetry) and with an increasing variety of observational geometries. In this paper, we report on Cassini’s radar instrument observations of the dune fields mapped through May 2009 and present our key findings in terms of Titan’s geology and climate. We estimate that dune fields cover ∼12.5% of Titan’s surface, which corresponds to an area of ∼10 million km², roughly the area of the United States. If dune sand-sized particles are mainly composed of solid organics as suggested by VIMS observations (Cassini Visual and Infrared Mapping Spectrometer) and atmospheric modeling and supported by radiometry data, dune fields are the largest known organic reservoir on Titan. Dune regions are, with the exception of the polar lakes and seas, the least reflective and most emissive features on this moon. Interestingly, we also find a latitudinal dependence in the dune field microwave properties: up to a latitude of ∼11°, dune fields tend to become less emissive and brighter as one moves northward. Above ∼11° this trend is reversed. The microwave signatures of the dune regions are thought to be primarily controlled by the interdune proportion (relative to that of the dune), roughness and degree of sand cover. In agreement with radiometry and scatterometry observations, SAR images suggest that the fraction of interdunes increases northward up to a latitude of ∼14°. In general, scattering from the subsurface (volume scattering and surface scattering from buried interfaces) makes interdunal regions brighter than the dunes. The observed latitudinal trend may therefore also be partially caused by a gradual thinning of the interdunal sand cover or surrounding sand sheets to the north, thus allowing wave penetration in the underlying substrate. Altimetry measurements over dunes have highlighted a region located in the Fensal dune field (∼5° latitude) where the icy bedrock of Titan is likely exposed within smooth interdune areas. The hemispherical assymetry of dune field properties may point to a general reduction in the availability of sediments and/or an increase in the ground humidity toward the north, which could be related to Titan’s asymmetric seasonal polar insolation. Alternatively, it may indicate that either the wind pattern or the topography is less favorable for dune formation in Titan’s northern tropics.     

On the Collision Nature of Two Coronal Mass Ejections: A Review

Observational and numerical studies have shown that the kinematic characteristics of two or more coronal mass ejections (CMEs) may change significantly after a CME collision. The collision of CMEs can have a different nature, i.e. inelastic, elastic, and superelastic processes, depending on their initial kinematic characteristics. In this article, we first review the existing definitions of collision types including Newton’s classical definition, the energy definition, Poisson’s definition, and Stronge’s definition, of which the first two were used in the studies of CME–CME collisions. Then, we review the recent research progresses on the nature of CME–CME collisions with the focus on which CME kinematic properties affect the collision nature. It is shown that observational analysis and numerical simulations can both yield an inelastic, perfectly inelastic, merging-like collision, or a high possibility of a superelastic collision. Meanwhile, previous studies based on a 3D collision picture suggested that a low approaching speed of two CMEs is favorable for a superelastic nature. Since CMEs are an expanding magnetized plasma structure, the CME collision process is quite complex, and we discuss this complexity. Moreover, the models used in both observational and numerical studies contain many limitations. All of the previous studies on collisions have not shown the separation of two colliding CMEs after a collision. Therefore the collision between CMEs cannot be considered as an ideal process in the context of a classical Newtonian definition. In addition, many factors are not considered in either observational analysis or numerical studies, e.g. CME-driven shocks and magnetic reconnections. Owing to the complexity of the CME collision process, a more detailed and in-depth observational analysis and simulation work are needed to fully understand the CME collision process.     

泥河湾盆地虎头梁剖面河湖相沉积物磁组构特征及其环境意义

泥河湾盆地在第四纪环境变化和古人类演化研究领域占有重要地位。虎头梁剖面记录了中更新世以来泥河湾盆地逐渐消亡的过程。本文对虎头梁天然露头剖面进行了详细的沉积学、岩石磁学和磁组构研究。剖面沉积物中的磁性矿物以假单畴（PSD）的磁铁矿为主,还含有一定量的赤铁矿,磁组构类型为正常沉积磁组构。根据虎头梁剖面沉积相的划分和磁组构参数特征,将泥河湾盆地中-晚更新世晚期的演化过程划分为3个阶段：下部为第Ⅰ阶段（60~53m）,中部为第Ⅱ阶段（53~12m）,上部为第Ⅲ阶段（12~0m）,包含3个完整湖侵旋回,每个旋回由下部灰黄色粗粒三角洲相到中部黄灰色滨湖相再到顶部灰绿色细粒湖相组成,相应的磁组构参数值由大到小的变化规律,中部第Ⅱ阶段推断由气候因素主导,上部第Ⅲ阶段可能主要由构造驱动,长期稳定的优势古水流方向以NW向为主,而东坡遗址时期的古人类可能生活在湖侵阶段的相对短期湖退的滨岸地带。

     

MULTI-DIMENSIONAL NUMERICAL STUDIES OF CONVECTION-RELATED STELLAR PROBLEMS

1　IntroductionManybottlenecksinthetheoreticalstudyofstellarevolutionarisefromthelackofunder standingofsomehighlycomplicatedprocessesthatcannotbetreatedwithanalyticaltech niques .Alargeportionoftheseprocessesisrelatedtohydrodynamics,particularlyconvec tion .Theapplicationofmultidimensionalnumericalcomputationstotheinvestigationofthesedifficultproblemshasmadegreatprogressinrecentyears .Here ,weprovideasurveyoftherecentnumericalstudiesofconvectionrelatedproblems.Theobjectiveistoassistthosewhoar…     

Comparative Analysis of Seismic Response Characteristics of Pile-Soil-Stnicture Interaction System

The study on the earthquake-resistant performance of a pile-soil-structure interaction system is a relatively complicated and primarily important issue in civil engineering practice. In this paper, a computational model and computation procedures for pile-supported structures, which can duly consider the pile-soil interaction effect, are established by the finite element method. Numerical implementation is made in the time domain. A simplified approximation for the seismic response analysis of pile-soil-structure systems is briefly presented. Then a comparative study is performed for an engineering example with numerical results computed respectively by the finite element method and the simplified method. Through comparative analysis, it is shown that the results obtained by the simplified method well agree with those achieved by the finite element method. The numerical results and findings will offer instructive guidelines for earthquake-resistant analysis and design of pile-supported structures.     

Numerical simulation of the diurnal cycle of rainfall in SW Amazon basin during the 1999 rainy season: the role of convective trigger function

In continental areas, the maximum rainfall simulated with the Brazilian developments on the Regional Atmospheric Modeling System (BRAMS) occurs around 4?h earlier than the one observed with rain gauges. This work presents the successful implementation of a new convective trigger function (CTF) in the convective parameterization scheme used in BRAMS that corrects this misfit between model and observations. The importance of the CTF formulation on the diurnal cycle of rainfall over the Amazon Basin is reflected by the following numbers: Over Rondonia (SW Amazonia), the original version of BRAMS simulates the maximum rainfall at 1400 UTC (1000 LST), with the new CTF maximum shifting to 1800?UTC (1400?LST), while the S-band radar rainfall maximum is at 1900?UTC (1500?LST). This is attributed to two factors: (1) the new CTF is now coupled to the sensible and latent heat fluxes at surface; (2) during the early morning, the convective available potential energy is reduced.     

西苕溪流域不同土地类型下氮元素输移过程

以西苕溪流域为研究对象,选择最有代表性的5种土地类型,模拟天然大暴雨,通过3次重复实验研究不同形态氮素随暴雨径流及径流沉积物的迁移过程,估算氮素在流域内不同土地利用/土地覆被条件下的损失率。研究结果表明,在相同的降雨条件下,氮素的流失速率和流失量随土地利用/土地覆被类型的不同表现出明显差异,地表径流水相总氮的流失量桑林最大,水田最小。水相不同形态的氮素流失量亦有所不同,悬浮颗粒态氮占地表径流水相总氮的70 %~90 %,水相溶解态氮的流失量以松林为最高,竹林、桑林和水田接近而且较低,不同类型的水相溶解态氮也随土地利用类型的不同表现出各自的特征。各土地类型单位面积、表层10 cm土壤氮素流失高达4.66~9.40 g·m-2,其中随径流沉积物相迁移的氮素占绝大部分(90 % 以上)。估算出的各土地利用类型总氮流失速率,地表径流水相为2.68~14.48 mg·m-2·min-1,径流沉积物相高达100.01~172.67 mg·m-2·min-1。     

中国城市地理学近期进展

中国城市地理学的兴盛时期始于70年代中期,至80年代进入历史以来最旺盛的发展时期。城市地理学主要研究城市体系、城市化、城市地域结构与具体的大小城镇。其重要的实践领域则主要通过城市规划和区域规划、国土规划。特别是70年代中期借助于我国城市规划工作的复兴,为城市地理学的发展提供了广阔的实践基础。     

Prompt optical emission and synchrotron self-absorption constraints on emission site of GRBs

We constrain the distance of the gamma-ray burst (GRB) prompt emission site from the explosion centre R , by determining the location of the electron's self-absorption frequency in the GRB prompt optical-to-X/γ-ray spectral energy distribution, assuming that the optical and the γ-ray emissions are among the same synchrotron radiation continuum of a group of hot electrons. All possible spectral regimes are considered in our analysis. The method has only two assumed parameters, namely the bulk Lorentz factor of the emitting source Γ and the magnetic field strength B in the emission region (with a weak dependence). We identify a small sample of four bursts that satisfy the following three criteria: (1) they all have simultaneous optical and γ-ray detections in multiple observational time intervals, (2) they all show temporal correlations between the optical and γ-ray light curves and (3) the optical emission is consistent with belonging to the same spectral component as the γ-ray emission. For all the time intervals of these four bursts, it is inferred that   R ≥ 10¹⁴  (Γ/300)^3/4 ( B /10⁵ G)^1/4  cm. For a small fraction of the sample, the constraint can be pinned down to   R ≈ 10¹⁴–10¹⁵ cm  for  Γ∼ 300  . For a second sample of bursts with prompt optical non-detections, only upper limits on R can be obtained. We find no inconsistency between the R -constraints for this non-detection sample and those for the detection sample.