导航战背景下的天文导航技术—天文导航技术的历史、现状及其发展趋势

综述近30年来我国内外天文导航技术及其相关理论和方法的进展，介绍世界各军事大国对天文导航的认识，并就天文导航在导航战背景下的显著优势进行论述。包括下列主要内容：（1）天文导航设备的历史和现状；（2）天文导航理论的新成果；（3）各军事大国对天文导航的态势；（4）天文导航的优势性；（5）天文导航自动化的主要关键技术，（6）天文导航自动化技术的发展趋势。     

CPLD技术及其在现代时统中的应用

对CPLD(复杂可编程逻辑器件)技术的基本特征和发展趋势作了简要介绍,揭示了该技术在现代数字系统中的重要地位及作用.利用CPLD对时统设备IRIG-B码产生器进行集成,其实验结果表明,集成了的B码产生器不但简单、可靠,而且便于调试,克服了以往硬件电路复杂的缺点.     

Time-spectra of chromospheric activity of old solar-type stars: detection of rotational signals from double wavelet analysis

We introduce a novel technique, called the double wavelet analysis (DWA), for the determination of stellar rotation periods from time serial data. This first paper aims narrowly at the discussion, introduction and application of the DWA technique to records of surface magnetism in solar-type (relatively old) lower main sequence stars that are obtained by the Mount Wilson Observatory (MWO) HK Project. The technique takes a series of careful steps that seek to optimize wavelet parameters and normalization schemes, ultimately allowing fine-tuned, arguably more accurate, estimates of rotation-modulated signals (with, e.g., periods of days to months) in records that contain longer periodicities such as stellar magnetic activity cycles (with, e.g., period of years). The apparent rotation periods estimated from the DWA technique are generally consistent with results from both “first-pass” (i.e., ordinary) global wavelet spectrum and earlier classical periodogram analyses. But there are surprises as well. For example, the rotation period of the ancient subdwarf Goombridge 1830 (HD 103095), previously identified as ≈31 days, suggests under the DWA technique a significantly slower period of 60 days. DWA spectra also generally reveal a shift in the cycle period toward high frequencies (hence shorter periods) compared to the first-pass wavelet spectrum. For solar-type stars analyzed here, the character of the DWA spectrum and slope of the first-pass global wavelet spectrum produce a classification scheme that allows a star's record to be placed into one of three categories.     

An Estimate for the Size of Cycle 23 Using Multivariate Relationships

A new technique based on multivariate analysis is described which allows for the prediction of the size or maximum amplitude of cycle 23. The technique uses the number of geomagnetic disturbances at selected times during the decline of cycle 22 (as precursors) to predict a maximum of about 152 (in terms of smoothed monthly mean sunspot numbers) for cycle 23. On the basis of this technique, hindcasts for cycles 17–22 are shown to agree with observed values within 5%.     

GPS掩星折射率剖面一维变分同化

近年来,GPS／LEO(全球定位系统／低地球轨道)卫星无线电掩星技术给出了地球大气探测的新途径．从LEO卫星观测到的掩星数据可以反演的地球大气的气压、水汽、温度等剖面;它们对气象和大气科学研究,是具有潜在价值的数据资源．掩星数据资料的同化技术可以有效地改进这些气象参数的剖面,从而改进目前的数值天气预报模式．在当前采用的一维变分同化反演技术中,可以用掩星观测资料的大气折射率或弯曲角剖面进行同化,来反演大气水汽和温度剖面以及海平面压强．作为独立自主开发的GPD／LEO掩星技术系统的一部分,以欧洲中尺度天气预报分析(ECMWF)资料为背景场,CHAMP 掩星观测得到的折射率剖面为观测值,采用Levenberg—Marquardt方法实行GPS掩星资料一维变分同化．在讨论中,用掩星观测点附近相应的探空气球资料来检验CHAMP掩星资料变分同化的结果．     

对称模式下的CHAMP弯曲角掩星数据同化

简单介绍了无线电掩星技术探测行星大气的发展史,列举了该技术中存在的若干问题。从 Eyre提出的统计学的最优估计反演方法,比较了用相位、弯曲角和折射率作为同化因子时出现的问题和各自的优缺点。对弯曲角同化因子,以欧洲中期天气预报中心(ECMWF)资料为背景场,运用一维变分技术,进行CHAMP掩星观测资料变分同化反演,从而获得水汽和温度剖面。将反演获得的气象剖面与非同化的剖面作比较,并且采用相应的探空气球资料作为验证,可以看出变分同化技术比传统的标准反演技术反演误差小。证实掩星数据资料的一维变分同化技术可以改进目前的数值天气预报模式。     

Elliptic Anomaly in Constructing Long-Term and Short-Term Dynamical Theories

The techniques of Brumberg and Brumberg (1999) based on the use of elliptic anomaly are specified in this paper in two aspects. The iteration technique (Broucke, 1969) to construct short-term semi-analytical theories of motion in rectangular coordinates in lines of Encke and Hill is reelaborated in terms of elliptic anomaly resulting in extending this technique for high-eccentricity orbits. In constructing long-term semi-analytical theories the key point is to integrate trigonometric functions of several angular arguments related to time by different differential expressions. This problem is reduced in the paper to linear algebraic recurrence relations admitting efficient solution by iterations.     

A technique for processing of planetary images with heterogeneous characteristics for estimating geodetic parameters of celestial bodies with the example of Ganymede

The new technique for generation of coordinate control point networks based on photogrammetric processing of heterogeneous planetary images (obtained at different time, scale, with different illumination or oblique view) is developed. The technique is verified with the example for processing the heterogeneous information obtained by remote sensing of Ganymede by the spacecraft Voyager-1, -2 and Galileo. Using this technique the first 3D control point network for Ganymede is formed: the error of the altitude coordinates obtained as a result of adjustment is less than 5 km. The new control point network makes it possible to obtain basic geodesic parameters of the body (axes size) and to estimate forced librations. On the basis of the control point network, digital terrain models (DTMs) with different resolutions are generated and used for mapping the surface of Ganymede with different levels of detail (Zubarev et al., 2015b).     

On 3D Reconstruction of Coronal Mass Ejections: I. Method Description and Application to SECCHI-COR Data

The data from SECCHI-COR1 and SECCHI-COR2 coronagraphs onboard the STEREO mission, which was launched in October 2006, provide us with the first-ever stereoscopic images of the Sun’s corona. These observations were found to be useful in inferring the three-dimensional structure of coronal mass ejections (CMEs) and their propagation direction in space. We apply four methods for reconstructing CMEs: i) Forward modeling technique; ii) Local correlation tracking (to identify the same feature in COR Ahead and COR Behind images) plus tie-point reconstruction technique; iii) Center of mass of the structures in a given epipolar plane plus tie-point reconstruction technique; iv) Polarization ratio technique. The four techniques are applied to three structured CMEs observed by COR1 and COR2 instruments, respectively, on 15 May 2007, 31 August 2007, and 25 March 2008. A comparison of the results obtained from the application of the four reconstruction algorithms is presented and discussed.     

Range-Doppler Mapping of Space-Based Targets Using the JRO 50 MHz Radar

The Jicamarca Radio Observatory (JRO) two-dimensional square array radar system operating at ~ 6-m wavelength was used to study the Moon and low Earth orbit satellites using the Range-Doppler inverse synthetic aperture radar technique also known as Delay-Doppler imaging. The radar data was collected on Oct 21, 2015. A circularly polarized coded pulse was transmitted from a quarter-array antenna segment during lunar transit over JRO. Dual-linear polarization receive systems were employed on two quarter-array segments and on two 1/64th array modules giving the longest possible baselines across the transit path. A Range-Doppler mapping technique that uses the rotational motion of the targets and an autofocusing motion and ionospheric delay compensation technique has been implemented to generate the two-dimensional maps of the point-target (Satellite) and range-spread target (Moon). A review of our technique and the maps obtained from these observations is presented herein. Range-Doppler maps of the Moon and satellites are instructive with regards to possible further improvement of the technique, especially regarding ionospheric compensation.     

Particle acceleration in reconnecting current sheets

We study motions of charged particles in reconnecting current sheets (CS) which have both transverse (perpendicular to the current sheet plane) and longitudinal (parallel to the electric current inside the sheet) components of the magnetic field. Such CS, called non-neutral, are formed in regions of magnetic field line reconnection in the solar atmosphere. We develop an analytical technique which allows us to reproduce previous results concerning the influence of transverse fields on particle motion and acceleration. This technique also allows us to evaluate the effect of the longitudinal field. The latter increases considerably the efficiency of particle acceleration in CS. The energizing of electrons during the main phase of solar flares can be interpreted as their acceleration in non-neutral CS.     

Correction of the ionospheric distortion on the MARSIS surface sounding echoes

Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) on the Mars Express (MEX) spacecraft has made numerous measurements of the Martian surface and subsurface. However, all of these measurements are distorted by the ionosphere and must be compensated before any analysis. We have developed a technique to compensate for the ionospheric distortions. This technique provides a powerful tool to derive the total electron content (TEC) and other higher-order terms of the limited expansion of the plasma dispersion function that are related to overall shape of the electron column profile. The derived parameters are fitted by using a Chapman model to derive ionospheric parameters like n₀, electron density primary peak (maximum for solar zenith angle (SZA) equal 0), and the neutral height scale H.

Our estimated ionospheric parameters are in good agreement with Mars Global Surveyor (MGS) radio-occultation data. However, since MARSIS does not have the observation geometry limitations of the radio occultation measurements, our derived parameters extend over a large range of SZA for each MEX orbit.

The first results from our technique have been discussed by Safaeinili et al. [2007, Estimation of the total electron content of the Martian ionosphere using radar sounder surface echoes. Geophys. Res. Lett. 34, L23204, doi:10.1029/2007GL032154].     

Status and Results from the HEGRA Air Shower Experiment

The HEGRA experiment investigates the relativistic (non thermal) universe by measuring air showers initiated in the atmosphere by photons and charged particles from the cosmos, using the imaging atmospheric Cherenkov light technique (> 500 GeV) and the showerfront sampling technique (> 20 TeV). New results concerning the search and study of gamma ray sources above 500 GeV, the measurement of the energy spectrum and chemical composition of cosmic rays above 1 PeV and the search of TeV counterparts of Gamma Ray Bursts are briefly described. This revised version was published online in July 2006 with corrections to the Cover Date.     

Automated Technique For Comparison Of Magnetic Field Inversion Lines With Filament Skeletons From The Solar Feature Catalogue

We present an automated technique for comparison of magnetic field inversion-line maps from SOHO/MDI magnetograms with solar filament data from the Solar Feature Catalogue created as part of the European Grid of Solar Observations project. The Euclidean distance transform and connected component labelling are used to identify nearest inversion lines to filament skeletons. Several filament inversion-line characteristics are defined and used to automate the decision whether a particular filament/inversion-line pair is suitable for quantitative comparison of orientation and separation. The technique is tested on 551 filaments from 14 Hα images at various dates, and the distributions of angles and distances between filament skeletons and line-of-sight (LOS) magnetic inversion lines are presented for six levels of magnetic field smoothing. The results showed the robustness of the developed technique which can be applied for a statistical analysis of magnetic field in the vicinity of filaments. The method accuracy is limited by the static filament detection which does not distinguish between filaments, fibrils, pre-condensations and filament barbs and this may increase the asymmetries in magnetic distributions and broadening in angular distributions that requires the incorporation of a feature tracking technique.     

Fast computation of 2D transfer functions from adaptive optics data

The use of atmospheric transfer functions is common in image reconstruction techniques such as speckle interferometry to calibrate the Fourier amplitudes of the reconstructed images. Thus, an accurate model is needed to ensure proper photometry in the reconstruction. The situation complicates when adaptive optics (AO) are used during data acquisition. I propose a novel technique to derive two‐dimensional transfer functions from data collected using AO simultaneously with the observations. The technique is capable to compute the relevant transfer functions within a short time for the prevailing atmospheric conditions and AO performance during data acquisition (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

DDS输出频谱特点及改善方法

利用非均匀采样模型得到了存在相位截断时直接数字频率合成（ＤＤＳ）的输出频谱,通过分析频谱特点,讨论了改善频谱分布的方法     

One-dimensional variational retrieval of refractivity profiles by GPS occultation

In recent years, the technique of Global Positioning System/Low Earth Orbit (GPS/LEO) radio occultation has opened up a new route for exploring the earth atmosphere. With this technique, the vertical profiles of pressure, temperature, water vapor of the earth atmosphere can be retrieved from GPS/LEO occultation data, thus making the latter a potentially valuable data resource for meteorological and atmospheric sciences. The technique may effectively improve the quality of the meteorological profiles, and so improve the current numerical weather forecasting model. In the one-dimensional variational (1DVAR) assimilation retrieval technique in current use, the profiles of atmospheric water vapor and temperature, as well as the sea level pressure are retrieved by assimilating the atmospheric refractivity or bending angle profile from GPS/LEO occultation data. In the system independently developed at Shanghai Observatory, the 1DVAR assimilation of CPS occultation data is accomplished by the Levenberg-Marquardt method with the analytical data from the European Center for Medium Range Weather Forecasting (ECMWF) as the background field, the refractivity profile obtained from CHAMP occultation data as the observational quantity. The result was checked with sounding balloon data near the occultation points.     

一种基于MUSIC算法的天地波识别方法

基于我国BPL长波脉冲信号的特征,利用MUSIC(多信号分类)算法对BPL天、地波延迟进行估计,实现天、地波识别。对传统谱估计IFFT(快速傅里叶逆变换)算法和现代谱估计MUSIC算法进行了仿真和比较,结果表明,这两种方法在较低信噪比条件下可有效分离天、地波,且识别误差都能控制在±5μs内,但MUSIC算法比IFFT算法具有更高的精度和分辨率。     

A comparative study of Newtonian,Kustaanheimo/Stiefel,and Sperling/Burdet optimal trajectories

An optimal trajectory problem is formulated in each of three sets of equations, and the resulting solutions are numerically compared. The three formulations are the classical Newtonian (N), the Kustaanheimo/Stiefel (K/S), and the Sperling/Burdet (S/B). The last two solutions are first regularized by the classical Sundman technique and the K/S solution is transformed before the optimization problem is posed. A novel technique is developed for generating initial control vectors for each solution. Numerically generated derivatives (central differences) are used by a type of gradient, Newton-Raphson iterator to converge the two-point boundary value problems. The results indicate that, although the K/S and S/B formulations are more difficult to express mathematically than the Newtonian formulation, the transformed solutions are significantly more numerically stable than the Newtonian solution when the perturbing acceleration is less than a minimum value (T/W _o=0.05 for the particular example problem treated).     

差分GPS在车辆导航系统中的应用

简单介绍了差分GPS和伪距差分技术在车辆导航系统中的应用原理,分析了车辆导航系统及其基准站与导航车辆的基本组成,并介绍了一种基于GPRS/ GSM通信网络的DGPS数据传输与处理过程。