一般概率测度的核估计

在测度弱收敛的意义下，研究了一般概率测度μ的核估计μn，得到了它的渐近无偏性、强相合性及在给定条件下的收敛速度。     

基于卫星观测数据的弹道导弹初始状态估计

利用卫星探测数据对大气层外飞行的弹道导弹弹道进行估计。在极大似然估计中,由于卫星探测导弹运动的可观测性较差,使初始点估计变得十分病态,传统的Gauss Newton迭代算法不再适用,而采用较为复杂Levenberg Marquardt算法,分别对单颗卫星和两颗卫星的情况进行了Monte Carlo仿真,得到了各自的归一化残差测量误差(代价函数),结果表明该算法有效。     

Study of ionospheric D region changes during solar flares using MF radar measurements

During solar flares, the X-ray radiation suddenly increases, resulting in an increase in the electron density of the atmospheric D region and a strong absorption of short-wave radio waves. Based on Langfang medium frequency (MF) radar, this paper analyzed the variation characteristics of D region in the lower ionosphere from 62 km to 82 km. The analysis focused on multiple C-level and M-level solar flare events before and after the large-scale flare event at 11:53 (UT) on September 6, 2017. The results show that it is difficult to detect the electron density over 70 km in Langfang during solar flares, but the electron density value can be obtained as low as 62 km, and the stronger the flare intensity, the lower the detectable electron density height. Besides, the equal electron density height, the received power of X and O waves will also be significantly reduced during the flares, and the reduction of equal electron density height has a weak linear relationship with flare intensity.     

Modeling and analysis of mover gaps in tubular moving-magnet linear oscillating motors

A tubular moving-magnet linear oscillating motor (TMMLOM) has merits of high efficiency and excellent dynamic capability. To enhance the thrust performance, quasi-Halbach permanent magnet (PM) arrays are arranged on its mover in the application of a linear electro-hydrostatic actuator in more electric aircraft. The arrays are assembled by several individual segments, which lead to gaps between them inevitably. To investigate the effects of the gaps on the radial magnetic flux density and the machine thrust in this paper, an analytical model is built considering both axial and radial gaps. The model is validated by finite element simulations and experimental results. Distributions of the magnetic flux are described in condition of different sizes of radial and axial gaps. Besides, the output force is also discussed in normal and end windings. Finally, the model has demonstrated that both kinds of gaps have a negative effect on the thrust, and the linear motor is more sensitive to radial ones.     

A new modeling scheme for powered parafoil unmanned aerial vehicle platforms: Theory and experiments

A novel framework is established for accurate modeling of Powered Parafoil Unmanned Aerial Vehicle (PPUAV). The model is developed in the following three steps: obtaining a linear dynamic model, simplifying the model structure, and estimating the model mismatch due to model variance and external disturbance factors. First, a six degree-of-freedom linear model, or the structured model, is obtained through dynamic establishment and linearization. Second, the data correlation analysis is adopted to determine the criterion for proper model complexity and to simplify the structured model. Next, an active model is established, combining the simplified model with the model mismatch estimator. An adapted Kalman filter is utilized for the real-time estimation of states and model mismatch. We finally derive a linear system model while taking into account of model variance and external disturbance. Actual flight tests verify the effectiveness of our active model in different flight scenarios.     

Orbit determination of BDS-3 satellite based on regional ground tracking station and inter-satellite link observations

The BeiDou global navigation satellite system (BDS-3) has established the Ka-band inter-satellite link (ISL) to realize a two-way ranging function between satellites, which provides a new observation technology for the orbit determination of BDS-3 satellites. Therefore, this study presents a BDS satellite orbit determination model based on ground tracking station (GTS) observations and ISL ranging observations firstly to analyze the impact of the ISL ranging observations on the orbit determination of BDS-3 satellites. Subsequently, considering the data fusion processing, the variance component estimation (VCE) algorithm is applied to the parameter estimation process of the satellite orbit determination. Finally, using the measured data from China’s regional GTS observations and BDS-3 ISL ranging observations, the effects of ISL ranging observations on the orbit determination accuracy of BDS-3 satellites are analyzed. Moreover, the impact of the VCE algorithm on the fusion data processing is evaluated from the aspects of orbit determination accuracy, Ka-band hardware delay parameter stability, and ISL ranging observation residuals. The results show that for China’s regional GTSs, the addition of BDS-3 ISL ranging observations can significantly improve the orbit determination accuracy of BDS-3 satellites. The observed orbit determination accuracy of satellite radial component is improved from 48 cm to 4.1 cm. In addition, when the initial weight ratio between GTS observations and ISL ranging observations is not appropriate, the various indicators which include orbit determination accuracy, ISL hardware delay, and ISL observation residuals were observed to have improved after the adjustment of the VCE algorithm. These results validate the effectiveness of the VCE algorithm for the fusion data processing of the GTS observations and ISL ranging observations.     

Thermosphere densities derived from Swarm GPS observations

After the detection of many anomalies in the Swarm accelerometer data, an alternative method has been developed to determine thermospheric densities for the three-satellite mission. Using a precise orbit determination approach, non-gravitational and aerodynamic-only accelerations are estimated from the high-quality Swarm GPS data. The GPS-derived non-gravitational accelerations serve as a baseline for the correction of the Swarm-C along-track accelerometer data. The aerodynamic accelerations are converted directly into thermospheric densities for all Swarm satellites, albeit at a much lower temporal resolution than the accelerometers would have been able to deliver. The resulting density and acceleration data sets are part of the European Space Agency Level 2 Swarm products.To improve the Swarm densities, two modifications have recently been added to our original processing scheme. They consist of a more refined handling of radiation pressure accelerations and the use of a high-fidelity satellite geometry and improved aerodynamic model. These modifications lead to a better agreement between estimated Swarm densities and NRLMSISE-00 model densities. The GPS-derived Swarm densities show variations due to solar and geomagnetic activity, as well as seasonal, latitudinal and diurnal variations. For low solar activity, however, the aerodynamic signal experienced by the Swarm satellites is very small, and therefore it is more difficult to accurately resolve latitudinal density variability using GPS data, especially for the higher-flying Swarm-B satellite. Therefore, mean orbit densities are also included in the Swarm density product.     

A comparative study of ionospheric IRIEup and ISP assimilative models during some intense and severe geomagnetic storms

Three-dimensional (3-D) electron density matrices, computed in the Mediterranean area by the IRI climatological model and IRIEup and ISP nowcasting models, during some intense and severe geomagnetic-ionospheric storms, were ingested by the ray tracing software tool IONORT, to synthesize quasi-vertical ionograms. IRIEup model was run in different operational modes: (1) assimilating validated autoscaled electron density profiles only from a limited area which, in our case, is the Mediterranean sector (IRIEup_re(V) mode); (2) assimilating electron density profiles from a larger region including several stations spread across Europe: (a) without taking care of validating the autoscaled data in the assimilation process (IRIEup(NV)); (b) validating carefully the autoscaled electron density profiles before their assimilation (IRIEup(V)).The comparative analysis was carried out comparing IRI, IRIEup_re(V), ISP, IRIEup(NV), and IRIEup(V) foF2 synthesized values, with corresponding foF2 measurements autoscaled by ARTIST, and then validated, at the truth sites of Roquetes (40.80°N, 0.50°E, Spain), San Vito (40.60°N, 17.80°E, Italy), Athens (38.00°N, 23.50°E, Greece), and Nicosia, (35.03°N, 33.16°E, Cyprus). The outcomes demonstrate that: (1) IRIEup_re(V), performs better than ISP in the western Mediterranean (around Roquetes); (2) ISP performs slightly better than IRIEup_re(V) in the central part of Mediterranean (around Athens and San Vito); (3) ISP performance is better than the IRIEup_re(V) one in the eastern Mediterranean (around Nicosia); (4) IRIEup(NV) performance is worse than the IRIEup(V) one; (5) in the central Mediterranean area, IRIEup(V) performance is better than the IRIEup_re(V) one, and it is practically the same for the western and eastern sectors.Concerning the overall performance, nowcasting models proved to be considerably more reliable than the climatological IRI model to represent the ionosphere behaviour during geomagnetic-ionospheric storm conditions; ISP and IRIEup(V) provided the best performance, but neither of them has clearly prevailed over the other one.     

Electron density inversed by plasma lines induced by suprathermal electron in the ionospheric modification experiment

Incoherent scatter radar (ISR) is the most powerful ground-based measurement facility to study the ionosphere. The plasma lines are not routinely detected by the incoherent scatter radar due to the low intensity, which falls below the measured spectral noise level of the incoherent scatter radar. The plasma lines are occasionally enhanced by suprathermal electrons through the Landau damping process and detectable to the incoherent scatter radar. In this study, by using the European Incoherent Scatter Association (EISCAT) UHF incoherent scatter radar, the experiment observation presents that the enhanced plasma lines were observed. These plasma lines were considered as manifest of the suprathermal electrons generated by the high-frequency heating wave during the ionospheric modification. The electron density profile is also obtained from the enhanced plasma lines. This study can be a promising technique for obtaining the accurate electron density during ionospheric modification experiment.     

Direction finding for two-dimensional incoherently distributed sources with Hadamard shift invariance in non-uniform orthogonal arrays

This paper proposes a novel algorithm for Two-Dimensional (2D) central Direction-of-Arrival (DOA) estimation of incoherently distributed sources. In particular, an orthogonal array structure consisting of two Non-uniform Linear Arrays (NLAs) is considered. Based on first-order Taylor series approximation, the Generalized Array Manifold (GAM) model can first be established to separate the central DOAs from the original array manifold. Then, the Hadamard rotational invariance relationships inside the GAMs of two NLAs are identified. With the aid of such relationships, the central elevation and azimuth DOAs can be estimated through a search-free polynomial rooting method. Additionally, a simple parameter pairing of the estimated 2D angular parameters is also accomplished via the Hadamard rotational invariance relationship inside the GAM of the whole array. A secondary but important result is a derivation of closed-form expressions of the Cramer-Rao lower bound. The simulation results show that the proposed algorithm can achieve a remarkably higher precision at less complexity increment compared with the existing low-complexity methods, which benefits from the larger array aperture of the NLAs. Moreover, it requires no priori information about the angular distributed function.