Optimal FIR filters for time-varying state-space models

An optimal FIR (finite impulse response) filter and smoother is introduced for the time-varying state-space model. The suggested filter has an FIR structure and utilizes finite observation. It is shown that the impulse response of the optimal FIR filter can be obtained by a simple Riccati-type matrix differential equation. Especially for time-invariant systems, this FIR filter reduces to previously known simple forms. For implementation, a recursive form of the optimal FIR filter and smoother is derived by using adjoint variables, and computational algorithms are suggested. It is also shown by sensitivity analysis that the proposed optimal FIR filter alleviates potential divergence characteristics of the standard Kalman filter     

Astrometric search for unseen stellar and sub-stellar companions to nearby stars and the possibility of their detection

Unseen companions to nearby stars are found astrometrically through perturbations in the proper motion from photographs taken with long-focus telescopes. The number of known unseen astrometric companions to nearby stars with photocentric orbits has grown by thirty percent in the last few years. Individual cases are discussed and optimum epochs given for resolution of the components. Orbital analysis of the photocentric positions on the photographic plates provides all information for accurate mass determination of the components except for m and angular separation, , of the two components which must come from another technique. There are potentially thirty low luminosity stars including some likely sub-stellar objects whose masses could be instantly found with the observations of these additional two parameters.A list of the stars known within five parsecs as of 1978 July is given and the status of unseen companions to these stars is discussed on the basis of long interval astrometric coverage.     

Slowly varying bias of magnetometer based on orbit geometry

A new method is proposed to correct the induced magnetic disturbance using a set of measured magnetic field data from the three-axis-magnetometer on orbit. The magnetic field correction is based on the orbit geometry and its algorithm excludes the geomagnetic field model. The considered satellite is supposed to operate in the Earth-point mode and the Sun-point mode in accordance with the mission requirements. Particularly, the magnetic field correction is performed via static estimation based on an average method and dynamic estimation with an estimation law. Besides, the usefulness of the proposed method is investigated throughout both the simulation and the real telemetry data of KOMPSAT-1.     

The carbonate-silicate cycle and CO2/climate feedbacks on tidally locked terrestrial planets

Atmospheric gaseous constituents play an important role in determining the surface temperatures and habitability of a planet. Using a global climate model and a parameterization of the carbonate-silicate cycle, we explored the effect of the location of the substellar point on the atmospheric CO(2) concentration and temperatures of a tidally locked terrestrial planet, using the present Earth continental distribution as an example. We found that the substellar point's location relative to the continents is an important factor in determining weathering and the equilibrium atmospheric CO(2) level. Placing the substellar point over the Atlantic Ocean results in an atmospheric CO(2) concentration of 7 ppmv and a global mean surface air temperature of 247 K, making ～30% of the planet's surface habitable, whereas placing it over the Pacific Ocean results in a CO(2) concentration of 60,311 ppmv and a global temperature of 282 K, making ～55% of the surface habitable.     

Flight performance analysis of GRACE K-band ranging instrument with simulation data

A dual one-way ranging (DOWR) system provides very high accuracy range measurements between two satellites. The GRACE satellite mission implements the DOWR, called KBR (K-band ranging), to measure very small inter-satellite range change in order to map the Earth gravity field. The flight performance of the KBR is analyzed by using a hybrid software simulator that incorporates actual satellite orbit data into a comprehensive KBR simulator, which was earlier used for computing the GRACE baseline accuracy. Three types of experiments were performed. First is the comparison of the flight data with the simulated data in spectral domain. Second is the comparison of double differenced noise level. Third is the comparison of the range-rate difference with GPS clock estimates. The analysis shows a good agreement with the simulation model except some excessive high frequency noise, e.g. 10⁻⁴ m/√Hz at 0.1 Hz. The range-rate difference shows 0.003 cyc/s discrepancy with the clock estimates. These analyses are helpful to refine the DOWR simulation model and can be benefit to future DOWR instrument development.     

基于非湍流脉动动能方程的高超声速转捩预测

通过求解基于稳定性理论的非湍流脉动动能输运方程,预测了来流Ma为5.91的裙锥绕流流动的层流-湍流转捩点,并采用SST(Shear-Stress Transport)湍流模型结合代数间歇因子模型对转捩流场进行了数值模拟.结果表明非湍流脉动动能输运方程模型可以很好地捕捉逆压梯度对转捩前期边界层内的不稳定波动频率的影响.等温壁面条件下计算得到的转捩位置与静风洞实验结果基本一致.壁面冷却促进第2模态不稳定波动的增长,使转捩提前发生.代数间歇因子模型模拟高马赫数流动的转捩区长度较短,温度峰值偏低.     

Algorithms for Power Converter Design Optimization

A practical optimization approach for power converters is established which allows conception of a design to meet all powercircuit performance requirements and concurrently to optimize a defined quantity such as weight or losses. In addition to facilitating a cost effective design, the computer-aided approach provides a means to readily assess a) the weight-efficiency tradeoff, b) impacts of converter requirements and component characteristics on a given design, and c) optimum power-system configurations. The following two popular algorithms for nonlinearly constrained optimization are utilized to design the power converter: 1) the sequential unconstrained minimization technique, SUMT Version 4; 2) the ALAG penalty function technique, ALAG5. These algorithms are compared, and suggestions are made for improving the efficiency of the optimization algorithms for power converter design.     

An LS-TLS-based sequential identification algorithm: Application toan accelerometer

Total least squares (TLS) parameter estimation is an alternative to least squares (LS) estimation when there are errors in both data matrix and observation vector. Especially, when some of the columns, not all, of the data matrix A are free of error, we call it a mixed LS-TLS problem. Accordingly, a sequential algorithm for solving a mixed LS-TLS problem is proposed here. The proposed algorithm employs an efficient algorithm to locate the minimum eigenpair, instead of singular value decomposition (SVD) which is computationally exacting. The proposed algorithm is applied to an accelerometer model to identify error parameters which are very important in inertial navigation systems (INS)     

Ulysses' race to the pole: Symposium Summary

This paper reproduces the Symposium Summary presentation at the 28th ESLAB Symposium on The High Latitude Heliosphere. The summary is not comprehensive but with considerable latitude covers the highlights of the conference. Featured are the observations of the high-speed wind from the growing southern polar coronal hole, new measurements of the interstellar gas and pickup ions, new measurements of cosmic ray modulation at high latitude, and our growing awareness of the three-dimensional heliosphere.     

Structure of reconnection layers in the magnetosphere

Magnetic reconnection can lead to the formation of observed boundary layers at the dayside magnetopause and in the nightside plasma sheet of the earth's magnetosphere. In this paper, the structure of these reconnection layers is studied by solving the one-dimensional Riemann problem for the evolution of a current sheet. Analytical method, resistive MHD simulations, and hybrid simulations are used. Based on the ideal MHD formulation, rotational discontinuities, slow shocks, slow expansion waves, and contact discontinuity are present in the dayside reconnection layer. Fast expansion waves are also present in the solution of the Riemann problem, but they quickly propagate out of the reconnection layer. Our study provides a coherent picture for the transition from the reconnection layer with two slow shocks in Petschek's model to the reconnection layer with a rotational discontinuity and a slow expansion wave in Levy et al's model. In the resistive MHD simulations, the rotational discontinuities are replaced by intermediate shocks or time-dependent intermediate shocks. In the hybrid simulations, the time-dependent intermediate shock quickly evolves to a steady rotational discontinuity, and the contact discontinuity does not exist. The magnetotail reconnection layer consists of two slow shocks. Hybrid simulations of slow shocks indicate that there exists a critical number,M _c, such that for slow shocks with an intermediate Mach numberM _IM _c, a large-amplitude rotational wavetrain is present in the downstream region. For slow shocks withM _I<M _c, the downstream wavetrain does not exist. Chaotic ion orbits in the downstream wave provide an efficient mechanism for ion heating and wave damping and explain the existence of the critical numberM _c in slow shocks.