多变量模型参考自适应控制及在飞行控制中的应用

讨论了多变量模型参考自适应控制及其在飞行控制中的应用问题.针对相对阶任意的多输入多输出线性时不变系统,提出了一种模型参考自适应控制方案.基于SDU分解,仅要求高频增益矩阵的顺序主子式的符号已知,从而放宽了对高频增益矩阵的限制.通过构造合适的Lyapunov函数,证明了闭环系统的全局稳定性和跟踪误差的渐进收敛性.将所提方案应用到了飞行控制系统,仿真结果验证了所提方案的有效性.     

基于空域划分的超视距空战态势威胁评估

编队超视距空战(BVR,Beyond Visual Range)已成为现代空战的主要模式.在空战优势区域与劣势区域判断的基础上对整个空域进行划分,并给出4种特定空域态势.从空中态势和编队作战能力两方面对空战态势进行分析.使用主成分分析法选取输入变量分析编队作战能力,降低评估过程中收集数据的复杂度.应用遗传神经网络对影响BVR各因素进行效能评估,将遗传算法(GA,Genetic Algorithms)与多层前馈(BP,Back Propagation)网络结合,利用GA的全局搜索优化BP网络的结构参数,有效克服BP算法的局部收敛等问题.结果表明:该模型能在综合分析空战各指标后给出红蓝双发的态势评估指标,该模型可有效减少评估中的人为因素,使评估结果更为客观可信.     

Comparison of ionosphere characteristic parameters obtained by ionosonde with IRI-2007 model over Southeast Asia

In this work, the foF2 and hmF2 parameters at the conjugate points near the magnetic equator of Southeast Asia are studied and compared with the International Reference Ionosphere (IRI) model. Three ionosondes are installed nearly along the magnetic meridian of 100°E; one at the magnetic equator, namely Chumphon (10.72°N, 99.37°E, dip angle 3.0°N), and the other two at the magnetic conjugate points, namely Chiang Mai (18.76°N, 98.93°E, dip angle 12.7°N) and Kototabang (0.2°S, 100.30°E, dip angle 10.1°S). The monthly hourly medians of the foF2 and hmF2 parameters are calculated and compared with the predictions obtained from the IRI-2007 model from January 2004 to February 2007. Our results show that: the variations of foF2 and hmF2 predicted by the IRI-2007 model generally show the similar feature to the observed data. Both parameters generally show better agreement with the IRI predictions during daytime than during nighttime. For foF2, most of the results show that the IRI model overestimates the observed foF2 at the magnetic equator (Chumphon), underestimates at the northern crest (Chiang Mai) and is close to the measured ones at the southern crest of the EIA (Kototabang). For hmF2, the predicted hmF2 values are close to the hmF2(M3000F2_OBS) during daytime. During nighttime, the IRI model gives the underestimation at the magnetic equator and the overestimation at both EIA crests. The results are important for the future improvements of the IRI model for foF2 and hmF2 over Southeast Asia region.     

Digital signal processing and numerical analysis for radar in geophysical applications

Numerical solutions for signal processing are described in this work as a contribution to study of echo detection methods for ionospheric sounder design. The ionospheric sounder is a high frequency radar for geophysical applications. The main detection approach has been done by implementing the spread-spectrum techniques using coding methods to improve the radar’s range resolution by transmitting low power. Digital signal processing has been performed and the numerical methods were checked. An algorithm was proposed and its computational complexity was calculated.     

Systematic examination of the geomagnetic storm sudden commencement using multi resolution analysis

Geomagnetic storm sudden commencements (SSC) contain a wealth of information which is useful in many applications. It is important to point out that the SSCs used in this study are sampled at the rate of one sample per second in order get use of such high resolution data. In this paper, two studies are made on the geomagnetic SSC and an SSC onset automatic detection algorithm is introduced. The first study is about finding the relationship between the SSC rise time and its amplitude. Where it is found that there is a positive correlation between the amplitude and the amplitude gradient which is the amplitude divided by rise time. The second study is checking the spectrum of the SSC, starting from its onset until the end of the SSC rise time. This check had proved that the SSC contains both low and high frequency regions. This led us to introduce a new term, namely the SSC variation rate (VR). This VR is defined as the maximum rate of change of the field in the higher-frequency region of the SSC. These two studies were the guide to build an SSC automatic detector of one sample per second data using multi resolution analysis (MRA) of the discrete wavelet transform (DWT). The data set contains 134 SSCs with different VRs that were collected from the Circum-pan Pacific Magnetometer Network (CPMN). It is found that the standard deviation of the detection error is 41 s and that the average error is 9 s. From the calculated error distribution function, it is found that the detection error is within the range of −1.5 to 3 min. The detection process, as will be shown in the article, takes 70 s for one station and 3 min if the decision is related to the detection(s) of other stations. These results demonstrate the superiority of the proposed algorithm over other algorithms in which the detection error ranges between −8 and 5 min and the detection process takes 2–10 min. In addition to being faster and more accurate than the other algorithms, the proposed algorithm is the first algorithm that automatically detects the SSC onset times from high-resolution data unlike previous studies that focused on determining the SSC times automatically using one-minute resolution data.     

Modelling and observing Jovian electron propagation times in the inner heliosphere

The propagation of Jovian electrons in interplanetary space was modelled by solving the relevant transport equation numerically through the use of stochastic differential equations. This approach allows us to calculate, for the first time, the propagation time of Jovian electrons from the Jovian magnetosphere to Earth. Using observed quiet-time increases of electron intensities at Earth, we also derive values for this quantity. Comparing the modelled and observed propagation times we can gauge the magnitude of the transport parameters sufficiently to place a limit on the 6 MeV Jovian electron flux reaching Earth. We also investigate how the modelled propagation time, and corresponding Jovian electron flux, varies with the well-known ∼13 month periodicity in the magnetic connectivity of Earth and Jupiter. The results show that the Jovian electron intensity varies by a factor of ∼10 during this cycle of magnetic connectivity.     

Testing the three-dimensional IRI-SIRMUP-P mapping of the ionosphere for disturbed periods

This paper describes the three-dimensional (3-D) electron density mapping of the ionosphere given as output by the assimilative IRI-SIRMUP-P (ISP) model for three different geomagnetic storms. Results of the 3-D model are shown by comparing the electron density profiles given by the model with the ones measured at two testing ionospheric stations: Roquetes (40.8°N, 0.5°E), Spain, and San Vito (40.6°N, 17.8°E), Italy. The reference ionospheric stations from which the autoscaled foF2 and M(3000)F2 data as well as the real-time vertical electron density profiles are assimilated by the ISP model are those of El Arenosillo (37.1°N, 353.3°E), Spain, Rome (41.8°N, 12.5°E), and Gibilmanna (37.9°N, 14.0°E), Italy. Overall, the representation of the ionosphere made by the ISP model is better than the climatological representation made by only the IRI-URSI and the IRI-CCIR models. However, there are few cases for which the assimilation of the autoscaled data from the reference stations causes either a strong underestimation or a strong overestimation of the real conditions of the ionosphere, which is in these cases better represented by only the IRI-URSI model. This ISP misrepresentation is mainly due to the fact that the reference ionospheric stations covering the region mapped by the model turn out to be few, especially for disturbed periods when the ionosphere is very variable both in time and in space and hence a larger number of stations would be required. The inclusion of new additional reference ionospheric stations could surely smooth out this concern.     

Estimation of ballistic coefficients of low altitude debris objects from historical two line elements

This paper presents a new method for estimating ballistic coefficients (BCs) of low perigee debris objects from their historical two line elements (TLEs). The method uses the drag perturbation equation of the semi-major axis of the orbit. For an object with perigee altitude below 700 km, the variation in the mean semi-major axis derived from the TLE is mainly caused by the atmospheric drag effect, and therefore is used as the source in the estimation of the ballistic coefficient. The method is tested using the GRACE satellites, and a number of debris objects with external ballistic coefficient values, and agreements of about 10% are achieved.     

Seasonally steady planetary disturbances in the troposphere and stratosphere as seen in 30 years of NCEP reanalysis data

Zonal velocity and temperature daily global reanalysis data of 30 years are used to search seasonally steady planetary disturbances in the middle troposphere (400 hPa) and middle stratosphere (10 hPa). Significant wavenumber 1, 2 and 3 modes are found. Constant phase lines of zonal velocity 1 modes exhibit significant inclination angles with respect to the meridians. The winter hemisphere generally shows a more significant presence of structures. The Northern Hemisphere (NH) exhibits all over the year a larger amount of structures and more intense amplitudes than the Southern Hemisphere (SH). Middle latitudes exhibit the most significant cases and low latitudes the least significant ones. Longitudinally oriented land–sea transitions at ±$\pm$65° and −35° latitudes appear to play a significant role for the presence of steady planetary modes. The stratosphere exhibits a much simpler picture than the troposphere. Large scale structures with respectively NE–SW (NH) and NW–SE (SH) tilts in the observed temperature and zonal velocity constant phase lines recall the quasi-stationary Rossby wave trains that favor the poleward transport of angular momentum.