卡尔曼滤波和自相关分析方法短期预报电离层f0F2的比较

电离层参数f₀F₂的预报是电离层研究的一个重要方向.本文选取2011年北京、长春、青岛和苏州四个常规观测站的f₀F₂数据,分别采用卡尔曼滤波、自相关分析法和国际参考电离层模型（IRI）对f₀F₂实施短期预报（提前1h预报）,并通过与实际观测数据的对比,对三种方法预报f₀F₂的性能进行了比较.研究结果表明,在磁平静时期,采用卡尔曼滤波方法进行f₀F₂预报的均方误差为0.532MHz,相对误差8.11%,比国际电离层参考模型的均方误差和相对误差分别降低1.47MHz和14.58%;采用自相关分析方法进行f₀F₂预报的均方根误差为0.967MHz,相对误差11.46%,比国际电离层参考模型的均方误差和相对误差分别降低1.035MHz和11.23%.比较结果说明二者对f₀F₂短期预报的精度相对于国际电离层参考模型均有大幅提升.对磁暴期间三种方法的预测性能做了进一步比较,试验结果表明卡尔曼滤波短期预报性能总体上优于自相关分析法,这为f₀F₂短期预报的方法选择提供了一定指导.      

Contribution analysis of QZSS to single-frequency PPP of GPS/BDS/GLONASS/Galileo

The Quasi-Zenith Satellite System (QZSS) established by the Japan Aerospace Exploration Agency mainly serves the Asia-Pacific region and its surrounding areas. Currently, four in-orbit satellites provide services. Most users of GNSS in the mass market use single-frequency (SF) receivers owing to the low cost. Therefore, it is meaningful to analyze and evaluate the contribution of the QZSS to SF precise point positioning (PPP) of GPS/BDS/GLONASS/Galileo systems with the emergence of GNSS and QZSS. This study compares the performances of three SF PPP models, namely the GRoup and PHase Ionospheric Correction (GRAPHIC) model, GRAPHIC with code observation model, and an ionosphere-constrained model, and evaluated the contribution of the QZSS to the SF PPP of GPS/BDS/GLONASS/Galileo systems. Moreover, the influence of code bias on the SF PPP of the BDS system is also analyzed. A two-week dataset (DOY 013–026, 2019) from 10 stations of the MGEX network is selected for validation, and the results show that: (1) For cut-off elevation angles of 15, 20, and 25°, the convergence times for the static SF PPP of GLONASS + QZSS are reduced by 4.3, 30.8, and 12.7%, respectively, and the positioning accuracy is similar compared with that of the GLONASS system. Compared with the BDS single system, the convergence times for the static SF PPP of BDS + QZSS under 15 and 25° are reduced by 37.6 and 39.2%, the horizontal positioning accuracies are improved by 18.6 and 14.1%, and the vertical components are improved by 13.9 and 21.4%, respectively. At cut-off elevation angles of 15, 20, and 25°, the positioning accuracy and precision of GPS/BDS/GLONASS/Galileo + QZSS is similar to that of GPS/BDS/GLONASS/Galileo. And the convergence times are reduced by 7.4 and 4.3% at cut-off elevation angles of 20 and 25°, respectively. In imitating dynamic PPP, the QZSS significantly improves the positioning accuracy of BDS and GLONASS. However, QZSS has little effect on the GPS-only, Galileo-only and GPS/BDS/GLONASS/Galileo. (2) The code bias of BDS IGSO and MEO cannot be ignored in SF PPP. In static SF PPP, taking the frequency band of B1I whose multipath combination is the largest among the frequency bands as an example, the vertical component has a systematic bias of approximately 0.4–1.0 m. After correcting the code bias, the positioning error in the vertical component is lower than 0.2 m, and the positioning accuracy in the horizontal component are improved accordingly. (3) The SF PPP model with ionosphere constraints has a better convergence speed, while the positioning accuracy of the three models is nearly equal. Therefore the GRAPHIC model can be used to get good positioning accuracy in the absence of external ionosphere products, but its convergence speed is slower.     

135.6nm夜气辉峰值电子密度反演算法及误差

根据夜间135.6nm大气辉光光强与F₂层峰值电子密度N_mF₂平方成正比的物理机制,在前期夜间135.6nm气辉辐射激发模型研究的基础上建立了峰值电子密度的反演算法,把全球经纬度分成若干格点,每个格点的电离层及中性成分信息分别由IRI2000和MSISE90提供,将电离层及中性成分廓线输入夜气辉辐射激发模型,计算每个格点135.6nm气辉的辐射强度,然后将各个格点的135.6nm气辉辐射强度与电离层廓线输入的N_mF₂平方拟合得到气辉强度与N_mF₂的转换因子.利用此方法可获得不同地方时、季节和太阳活动周期的转换因子组成查算表,进而根据实际探测的135.6nm气辉辐射强度反演相应时空的N_mF₂.最后对该算法的反演误差进行了综合分析,为该算法适用的时空特性提供重要理论支撑.      

一种双捷联微小型惯性组合导航系统设计研究

介绍了一种使组合导航系统具有系统级冗余的基于MEMS-IMU的双SINS／GPS／Doppler组合导航系统，讨论了系统的体系结构及分层联合卡尔曼滤波器的设计和计算方法，给出了基于滤波器测量残差检验的硬故障FDIR算法和基于状态估计的软故障FDIR算法。计算机仿真结果表明，该系统具有结构灵活、计算方便、可靠性高的优点。     

空间遥感器热光学分析的研究与应用

论述空间遥感器热光学分析与应用现状及其发展趋势;分析当前遥感器热光学分析的意义,认为这是提高长寿命遥感器可靠性的一种重要手段;针对某长寿命遥感相机,提出进行热光学分析研究的总体方案、具体实施途径、技术关键和解决方法。     

信阳陨石热历史与原始氩特征

信阳陨石于1977年12月1日18时57分, 陨落在我国河南省信阳县肖王公社。陨落后, 破裂为X-Ⅰ(48kg)及X-Ⅱ(27.5kg)两大块和几小块碎屑。陶克捷等对该陨石做过化学组成、矿物组分及结构特征研究, 指出属H6型。本文试图根据40Ar/39Ar快中子活化资料及原始捕获氩同位素特征, 探讨信阳陨石的演化特点。      

一种并行式CPL/LHP冷凝器的流动/传热分析

针对航天器的可展开式辐射器 ,对比分析了辐射器内冷凝管路的两种布局方式 ,并利用SINDA/FLUINT对微重力条件下并行式冷凝器内的两相流动 /传热进行了计算分析 ;研究了冷凝器后端节流管对冷凝管内汽 /液界面位置的影响 ,为冷凝器的设计和提高辐射器性能提供了相应的依据     

Interconnection of high-latitude and low-latitude boundary layers when IMF BY is dominant

The antiparallel merging model places the location of the reconnection region for a dominant interplanetary magnetic field (IMF) B_Y at high latitudes at the dayside magnetopause and predicts that the low-latitude boundary layer (LLBL) is located on open field lines of the magnetospheric flanks. Interball-1 data obtained in the wide local time range near the low-latitude magnetopause makes it possible to analyze the LLBL plasma population and to find a link between possible reconnection at high latitudes and LLBL occurrence. We found that no boundary layer was observed in the regions which have no topological connection with the merging site. All cases of LLBL observations are located downstream from a specific boundary. This boundary coincides with the first magnetospheric field line touching the reconnection region and can be located in a wide local time region depending on the instant IMF direction. Even the LLBL on closed field lines shows the tendency to be concentrated in the vicinity of this boundary. Thus we show that all types of observed LLBLs are linked to reconnection sites predicted by the antiparallel merging model.     

Estimating inter-system biases for tightly combined Galileo/BDS/GPS RTK

The overlapping carrier frequencies L1/E1, L5/E5a and B2/E5b from GPS/Galileo/BDS allow inter-system double-differencing of observations, which shows a clear advantage over differencing of the observations of each constellation independently. However, the inter-system biases destroy the integer nature of the inter-system double-differencing ambiguities. Two methods of direct rounding and parameter estimation are used to determine the ISB value. By analyzing data collected from Curtin University from 2015 to 2018, the phase and code inter-system bias (ISB) are related to the receiver type, firmware version and the selected overlapping frequency. Upgrade of receiver firmware version results in changes of ISB values. For example, the upgrade of Javad firmware in Dec, 15, 2017 causes the difference of 0.5 cycles ISB between BDS GEO and non-GEO satellites. By comparing the three dynamic models which include white noise process, random walk process, and random constant in the parameter estimation method, the ISB determined by the random constant model is consistent with the value obtained by the direct rounding method. After the calibration of ISBs, the performances of tightly combined positioning are assessed. The success rate of ambiguity resolution and accuracy of positioning for the tight combination (TC) are significantly improved in comparison with that for the loose combination (LC) over short baselines. For L5/E5a, on which only few satellites can be observed, the maximum increase in success rates of ambiguity resolution can reach 31.7%, i.e., from 54.9% of LC to larger than 86.6% of TC, and the positioning accuracies can even be increased by 0.13 m, i.e., from 0.208 of LC to 0.074 m of TC in East direction for the mix-receiver TRIMBLE NETR9-SEPT POLARX4 in 2018.     

Aura/MLS与TIMED/SABER观测全球重力波特性

大气重力波是临近空间环境主要大气波动之一,对全球环流具有重要影响。卫星上搭载的临边探测器能够探测临近空间大气温度,可用于临近空间大气重力波研究。利用2012－2014年Aura的微波临边探测器（MLS）和TIMED的红外临边探测器（SABER）的探测数据,对20~50 km高度的大气重力波扰动分布特征开展了分析研究,两种观测重力波活动基本一致,重力波随季节、纬度及高度的变化显著。冬季半球高纬度重力波扰动较强,赤道和夏季半球近赤道地区上空也存在明显重力波活动区域,夏季半球高纬度重力波扰动最弱。重力波扰动强度随高度增加。TIMED/SABER重力波扰动强度数值比 Aura/MLS略强。      

The model of turbulent plasma sheet during IMF Bz > 0

Theory of the plasma sheet with medium-scale developed turbulence gives the possibility to explain the main processes of plasma sheet bifurcation and theta-aurora formation during IMF B_z > 0. The model suggests that during IMF B_z > 0 small bulge structure in the plasma sheet center is formed. The polarization of the bulge due to dawnward electron motion and duskward ion motion decreases the large-scale electric field in the bulge region. The decrease of the large-scale field in the conditions of constant coefficient of diffusion leads to the bulge growth. The results of plasma sheet bifurcation and theta-aurora formation modelling are presented.     

Tight integration of BDS-3/BDS-2/GPS/Galileo observations considering the new overlapping DISBs and its application in obstructed environments

Tight integration can enhance the model strength and positioning performance by considering the characteristic of differential inter-system bias (DISB), especially in obstructed environments. However, limited work emphasizes the comprehensive analysis of five-frequency DISBs between BDS-3 and other systems considering the receiver type, receiver configuration, and antenna type. In addition, the overlapping DISBs between BDS-3 and BDS-2 are also in great demand for further investigation since they are often regarded as one system. In this study, one DISB-float model is introduced to estimate the DISBs, and one DISB-fixed model and one DISB-free model are formulated to enhance the model strength of tight integration. Four dedicated datasets were collected to estimate the DISBs, which are also comprehensively analyzed considering the receiver type, receiver configuration, and antenna type. The results show that the DISBs between BDS-3 and other systems are rather stable over a certain period and are related to the receiver type and receiver configuration, whereas are not related to the antenna type. More interestingly, the B1I code DISB between BDS-3 and BDS-2 exhibits significant magnitude with a mean value of ?1.44 m for the baseline composed of two different receivers. In this case, the B1I code DISB must be considered and the tight integration between BDS-3 and BDS-2 considering its calibration can improve the positioning performance. Besides, the tight integration of the DISB-fixed model can significantly improve the positioning accuracy between multiple GNSS. Compared to the loose integration, the improvement of 60.6 %, 56.6 %, and 61.2 % can be obtained in the E, N, and U directions, when only two satellites are available for each system. In real obstructed environments, the tight integration of the DISB-free model can also improve the positioning performance in terms of positioning availability and accuracy, as well as the ambiguity resolution performance.     

Near real-time PPP-based monitoring of the ionosphere using dual-frequency GPS/BDS/Galileo data

Ionosphere delay is very important to GNSS observations, since it is one of the main error sources which have to be mitigated even eliminated in order to determine reliable and precise positions. The ionosphere is a dispersive medium to radio signal, so the value of the group delay or phase advance of GNSS radio signal depends on the signal frequency. Ground-based GNSS stations have been used for ionosphere monitoring and modeling for a long time. In this paper we will introduce a novel approach suitable for single-receiver operation based on the precise point positioning (PPP) technique. One of the main characteristic is that only carrier-phase observations are used to avoid particular effects of pseudorange observations. The technique consists of introducing ionosphere ambiguity parameters obtained from PPP filter into the geometry-free combination of observations to estimate ionospheric delays. Observational data from stations that are capable of tracking the GPS/BDS/GALILEO from the International GNSS Service (IGS) Multi-GNSS Experiments (MGEX) network are processed. For the purpose of performance validation, ionospheric delays series derived from the novel approach are compared with the global ionospheric map (GIM) from Ionospheric Associate Analysis Centers (IAACs). The results are encouraging and offer potential solutions to the near real-time ionosphere monitoring.     

Optical/radio/pulsars integrated navigation for Mars orbiter

In this paper, we address the issue of the integrated navigation algorithm with different combination of measurements for Mars orbiter. First, system dynamic model and navigation measurement models using optical measurement information, radio measurement information and X-ray pulsars measurement information are respectively established. Second, optical/radio/pulsars integrated navigation algorithm is proposed, and observability analysis of the integrated navigation system is also conducted. Third, adaptive extended Kalman filter is adopted to fuse measurement information and suppress measurement and process noise to optimally estimate the state of Mars orbiter. Monte Carlo simulation results show that optical/radio/pulsars integrated navigation can effectively improve the navigation accuracy and satisfy the navigation requirements of Mars orbiter.     

A different method to update monthly median hmF2 values

The height, h_mF2, and the electron density, N_mF2, of the F2 peak are key model parameters to characterize the actual state of the ionosphere. These parameters, or alternatively the propagation factor, M3000F2, and the critical frequency, f_oF2, of the F2 peak, which are related to h_mF2 and N_mF2, are used to anchor the electron density vertical profile computed with different models such as the International Reference Ionosphere ( Bilitza, 2002), as well as for radio propagation forecast purposes. Long time series of these parameters only exist in an inhomogeneous distribution of points over the surface of Earth, where dedicated instruments (typically ionosondes) have been working for many years. A commonly used procedure for representing median values of the aforementioned parameters all over the globe is the one recommended by the ITU-R ( ITU-R, 1997). This procedure, known as the Jones and Gallet mapping technique, was based on ionosondes measurements gathered from 1954 to 1958 by a global network of around 150 ionospheric stations (  and ). Even though several decades have passed since the development of that innovative work, only few efforts have been dedicated to establish a new mapping technique for computing h_mF2 and N_mF2 median values at global scale or to improve the old method using the increased observational database. Therefore, in this work three different procedures to describe the daily and global behavior of the height of the F2 peak are presented. All of them represent a different and simplified method to estimate h_mF2 and are based on different mathematical expressions. The advantages and disadvantages of these three techniques are analyzed, leading to the conclusion that the recommended procedure to represent h_mF2 is best characterized by a Spherical Harmonics expansion of degree and order equal to 15, since the differences between the h_mF2 values obtained with the Jones and Gallet technique and those obtained using the abovementioned procedure are of only 1%.     

基于IMU/FADS/无线电的火星上升器组合导航方案

以火星采样返回任务中火星表面上升为背景,研究了基于惯性测量单元（Inertial Measurement Unit, IMU）、嵌入式大气数据传感系统（Flush Air Data Sensing System, FADS）和无线电信标的组合导航方法。首先,在传统的IMU导航框架中加入由无线电测量获得的相对距离、速度信息,以及由FADS获取的动压、温度数据,建立了基于IMU、无线电和FADS的导航观测模型;然后,基于无迹卡尔曼滤波（Unscented Kalman Filter, UKF）技术对测量信息进行了融合,并压制了过程噪声和测量噪声,从而对上升器的状态进行了联合估计;最后,在数值仿真中,将UKF与自适应无迹卡尔曼滤波（Adaptive Unscented Kalman Filter, AUKF）技术进行了对比,在比较不同滤波器性能的同时,验证了组合导航方法的有效性。     

Radio/FADS/IMU integrated navigation for Mars entry

Supposing future orbiting and landing collaborative exploration mission as the potential project background, this paper addresses the issue of Mars entry integrated navigation using radio beacon, flush air data sensing system (FADS), and inertial measurement unit (IMU). The range and Doppler information sensed from an orbiting radio beacon, the dynamic pressure and heating data sensed from flush air data sensing system, and acceleration and attitude angular rate outputs from an inertial measurement unit are integrated in an unscented Kalman filter to perform state estimation and suppress the system and measurement noise. Computer simulations show that the proposed integrated navigation scheme can enhance the navigation accuracy, which enables precise entry guidance for the given Mars orbiting and landing collaborative exploration mission.     

A survey of radio jets with Chandra and HST

Relativistic jets are a common property of radio-loud Active Galactic Nuclei (AGN). Understanding jet physical properties is an essential precursor to understanding the mechanisms of energy transport, and ultimately, how energy is extracted from the central black hole. In this paper, I highlight recent developments from Chandra and HST observations of kpc-scale jets in AGN, with particular emphasis on our survey of 17 radio jets in a sample of FRII radio galaxies. These observations show that (1) X-ray and optical emission is common from kpc-scale jets, (2) a large fraction of the bolometric luminosity is emitted at X-rays, and (3) in most sources, a candidate emission process for the X-rays is inverse Compton scattering of the Cosmic Microwave Background off the relativistic electrons in the jet. If the latter scenario holds, the implication is that jets are still relativistic on kpc scales.     

氢氧发动机推力室化学反应流场计算

采用弱耦合点隐式方法的MacCormack格式对氢氧火箭发动机推力室化学反应粘性流场进行数值模拟.数值模拟时采用6种组分、8个反应有限速率的化学反应模型,湍流模型采用Baldwin-Lomax代数湍流模型.数值模拟得到了流场参数在燃烧室和喷管中的分布.结果分析表明,得到的数值模拟结果与理论分析一致,说明结果是可靠的.本文的工作为氢氧火箭发动机喷管设计提供了依据,并为进一步开展火箭发动机推力室有化学反应的两相流动的数值模拟打下了基础.