Discrete-time observability and estimability analysis for bearings-only target motion analysis

Observability in the context of bearings-only tracking (BOT) is still the subject of important literature. Different from previous approaches, where continuous-time analysis was considered, our approach relies on discrete-time analysis. It is then shown that this allows us to use directly and efficiently the simple formalisms of linear algebra. Using the direct approach, observability analysis is essentially reduced to basic considerations about subspace dimensions. Even if this approach is conceptually quite direct, it becomes more and more complex as the source-encounter scenario complexity increases. For complex scenarios, the dual approach may present some advantages essentially due to the direct use of multilinear algebra. New results about BOT observability for maneuvering sources are thus obtained. Observability analysis is then extended to unknown instants of source velocity changes. Even if observability analysis provides thorough insights about the algebraic structure of the BOT problem, the optimization of the observer maneuvers is essentially a control problem. Basic algebraic considerations prove that a relevant cost functional for this control problem is the determinant of the Fisher information matrix (FIM). So, a large part of this work is devoted to the analysis of this cost functional. Using multilinear algebra, general approximations of this functional are given. In order to involve only directly estimable parameters, the source bearing-rates are examined. Using these approximations, a general framework for optimizing the observer trajectory is derived which allow us to approximate the optimal sequence of controls. It is worth stressing that our approach does not require the knowledge of the source trajectory parameters and is still valid for a maneuvering source.     

A New Mathematical Formulation for Strapdown Inertial Navigation

A differential equation is developed for the orientation vector relating the body frame to a chosen reference frame. The time derivative of this vector is the sum of the inertially measurable angular velocity vector and of the inertially nonmeasurable noncommutativity rate vector. It is precisely this noncommutativity rate vector that causes the computational problems when numerically integrating the direction cosine matrix. The orientation vector formulation allows the noncommutativity contribution to be isolated and, therefore, treated separately and advantageously. An orientation vector mechanization is presented for a strap down inertial system. Further, an example is given of the applica tion of this formulation to a typical rigid body rotation problem.     

固体火箭发动机寿命预估的一种考虑

利用弹性－粘弹性对应原理，分析长期贮存固体推进剂药柱由环境温度变化引起的热应力，计算药柱损伤，并对发动机寿命进行预估。理论分析和算法具有普遍性，可供工程设计参考。     

Medium PRF set selection using evolutionary algorithms

This paper presents a new and novel method of selecting multiple pulse repetition frequency (PRF) sets for use in medium PRF pulsed-Doppler radars. Evolutionary algorithms are used to minimise the blind areas in the range/Doppler space. The evolutionary algorithm allows optimal solutions to be generated quickly, far faster than with exhaustive searches, and is fully automatic, unlike existing techniques. The evolved solutions compare very favorably against the results of both an exhaustive search and existing published PRF set selection methods. This evolutionary approach to generation of PRF sets is a major advance in medium PRF radar design.     

A Correlation Polarimeter for Noise-Like Signals

Optimum estimation (tracking) of the polarization plane of a linearly polarized electromagnetic wave is determined when the signal is a narrow-band Gaussian random process with a polarization plane angle which is also a Gaussian random process. This model is Compared to previous work and is applicable to space communication. The estimator performs a correlation operation similar to an amplitude -comparison monopulse angle tracker, giving the name correlation polarimeter. Under large signal-to-noise ratio (SNR), the estimator is causal. Performance of the causal correlation polarimeter is evaluated for arbitrary SNR. Optimum precorrelation filtering is determined. With low SNR, the performance of this system is far better than that of previously developed systems. Practical implementation is discussed. A scheme is given to reduce the effect of linearly polarized noise.     

Magnetosphere Imaging Instrument (MIMI) on the Cassini Mission to Saturn/Titan

The magnetospheric imaging instrument (MIMI) is a neutral and charged particle detection system on the Cassini orbiter spacecraft designed to perform both global imaging and in-situ measurements to study the overall configuration and dynamics of Saturn’s magnetosphere and its interactions with the solar wind, Saturn’s atmosphere, Titan, and the icy satellites. The processes responsible for Saturn’s aurora will be investigated; a search will be performed for substorms at Saturn; and the origins of magnetospheric hot plasmas will be determined. Further, the Jovian magnetosphere and Io torus will be imaged during Jupiter flyby. The investigative approach is twofold. (1) Perform remote sensing of the magnetospheric energetic (E > 7 keV) ion plasmas by detecting and imaging charge-exchange neutrals, created when magnetospheric ions capture electrons from ambient neutral gas. Such escaping neutrals were detected by the Voyager l spacecraft outside Saturn’s magnetosphere and can be used like photons to form images of the emitting regions, as has been demonstrated at Earth. (2) Determine through in-situ measurements the 3-D particle distribution functions including ion composition and charge states (E > 3 keV/e). The combination of in-situ measurements with global images, together with analysis and interpretation techniques that include direct “forward modeling’’ and deconvolution by tomography, is expected to yield a global assessment of magnetospheric structure and dynamics, including (a) magnetospheric ring currents and hot plasma populations, (b) magnetic field distortions, (c) electric field configuration, (d) particle injection boundaries associated with magnetic storms and substorms, and (e) the connection of the magnetosphere to ionospheric altitudes. Titan and its torus will stand out in energetic neutral images throughout the Cassini orbit, and thus serve as a continuous remote probe of ion flux variations near 20R _S (e.g., magnetopause crossings and substorm plasma injections). The Titan exosphere and its cometary interaction with magnetospheric plasmas will be imaged in detail on each flyby. The three principal sensors of MIMI consists of an ion and neutral camera (INCA), a charge–energy–mass-spectrometer (CHEMS) essentially identical to our instrument flown on the ISTP/Geotail spacecraft, and the low energy magnetospheric measurements system (LEMMS), an advanced design of one of our sensors flown on the Galileo spacecraft. The INCA head is a large geometry factor (G ∼ 2.4 cm² sr) foil time-of-flight (TOF) camera that separately registers the incident direction of either energetic neutral atoms (ENA) or ion species (≥5^∘ full width half maximum) over the range 7 keV/nuc < E < 3 MeV/nuc. CHEMS uses electrostatic deflection, TOF, and energy measurement to determine ion energy, charge state, mass, and 3-D anisotropy in the range 3 ≤ E ≤ 220 keV/e with good (∼0.05 cm² sr) sensitivity. LEMMS is a two-ended telescope that measures ions in the range 0.03 ≤ E ≤ 18 MeV and electrons 0.015 ≤ E≤ 0.884 MeV in the forward direction (G ∼ 0.02 cm² sr), while high energy electrons (0.1–5 MeV) and ions (1.6–160 MeV) are measured from the back direction (G ∼ 0.4 cm² sr). The latter are relevant to inner magnetosphere studies of diffusion processes and satellite microsignatures as well as cosmic ray albedo neutron decay (CRAND). Our analyses of Voyager energetic neutral particle and Lyman-α measurements show that INCA will provide statistically significant global magnetospheric images from a distance of ∼60 R _S every 2–3 h (every ∼10 min from ∼20 R _S). Moreover, during Titan flybys, INCA will provide images of the interaction of the Titan exosphere with the Saturn magnetosphere every 1.5 min. Time resolution for charged particle measurements can be < 0.1 s, which is more than adequate for microsignature studies. Data obtained during Venus-2 flyby and Earth swingby in June and August 1999, respectively, and Jupiter flyby in December 2000 to January 2001 show that the instrument is performing well, has made important and heretofore unobtainable measurements in interplanetary space at Jupiter, and will likely obtain high-quality data throughout each orbit of the Cassini mission at Saturn. Sample data from each of the three sensors during the August 18 Earth swingby are shown, including the first ENA image of part of the ring current obtained by an instrument specifically designed for this purpose. Similarily, measurements in cis-Jovian space include the first detailed charge state determination of Iogenic ions and several ENA images of that planet’s magnetosphere.This revised version was published online in July 2005 with a corrected cover date.     

Deep Impact: Working Properties for the Target Nucleus – Comet 9P/Tempel 1

In 1998, Comet 9P/Tempel 1 was chosen as the target of the Deep Impact mission (A’Hearn, M. F., Belton, M. J. S., and Delamere, A., Space Sci. Rev., 2005) even though very little was known about its physical properties. Efforts were immediately begun to improve this situation by the Deep Impact Science Team leading to the founding of a worldwide observing campaign (Meech et al., Space Sci. Rev., 2005a). This campaign has already produced a great deal of information on the global properties of the comet’s nucleus (summarized in Table I) that is vital to the planning and the assessment of the chances of success at the impact and encounter. Since the mission was begun the successful encounters of the Deep Space 1 spacecraft at Comet 19P/Borrelly and the Stardust spacecraft at Comet 81P/Wild 2 have occurred yielding new information on the state of the nuclei of these two comets. This information, together with earlier results on the nucleus of comet 1P/Halley from the European Space Agency’s Giotto, the Soviet Vega mission, and various ground-based observational and theoretical studies, is used as a basis for conjectures on the morphological, geological, mechanical, and compositional properties of the surface and subsurface that Deep Impact may find at 9P/Tempel 1. We adopt the following working values (circa December 2004) for the nucleus parameters of prime importance to Deep Impact as follows: mean effective radius = 3.25± 0.2 km, shape – irregular triaxial ellipsoid with a/b = 3.2± 0.4 and overall dimensions of ∼14.4 × 4.4 × 4.4 km, principal axis rotation with period = 41.85± 0.1 hr, pole directions (RA, Dec, J2000) = 46± 10, 73± 10 deg (Pole 1) or 287± 14, 16.5± 10 deg (Pole 2) (the two poles are photometrically, but not geometrically, equivalent), Kron-Cousins (V-R) color = 0.56± 0.02, V-band geometric albedo = 0.04± 0.01, R-band geometric albedo = 0.05± 0.01, R-band H(1,1,0) = 14.441± 0.067, and mass ∼7×10¹³ kg assuming a bulk density of 500 kg m⁻³. As these are working values, {i.e.}, based on preliminary analyses, it is expected that adjustments to their values may be made before encounter as improved estimates become available through further analysis of the large database being made available by the Deep Impact observing campaign. Given the parameters listed above the impact will occur in an environment where the local gravity is estimated at 0.027–0.04 cm s⁻² and the escape velocity between 1.4 and 2 m s⁻¹. For both of the rotation poles found here, the Deep Impact spacecraft on approach to encounter will find the rotation axis close to the plane of the sky (aspect angles 82.2 and 69.7 deg. for pole 1 and 2, respectively). However, until the rotation period estimate is substantially improved, it will remain uncertain whether the impactor will collide with the broadside or the ends of the nucleus.     

Comparison of ≳40 keV Electron Events at High and Low Heliolatitudes: Ulysses/HI-SCALE and ACE/EPAM

We have performed a joint survey of anisotropic ≳40 keV electron events from August 1997 to September 2000 using the matched detectors on the Ulysses (ULS)/HI-SCALE and the ACE/EPAM instruments. A computer algorithm selected events with strong, statistically significant pitch-angle anisotropies. Electron pitch-angle distributions at ACE (∼1 AU) are often ‘beams’ that are strongly collimated along the local interplanetary magnetic field (IMF). These flare-associated impulsive injections can display rapid rise times (∼15 min) and slower decays, or more irregular intensity histories. At ULS, the electron intensities are lower and the time histories smoother, but strong anisotropies are still observable, indicating direct, nearly field-aligned propagation outward from the Sun. We focus on four event periods, selected from the survey, during times when the angle between the footpoints of the IMF lines intersecting ACE and ULS is small. These events span three full years and cover a wide range of distances and heliographic latitudes. We found one reasonably good association between impulsive electron events at ACE and ULS, and two events with small field-aligned gradients. This revised version was published online in August 2006 with corrections to the Cover Date.     

Orbits and Integrals in Self-Consistent Systems

We find the forms of the orbits in a self-consistent galactic model generated by a N-body simulation of the collapse of a protogalaxy. The model represents a stationary elliptical galaxy of type E5, which is approximately axisymmetric around its longest axis. The orbits are of three main types, box orbits (including box-like orbits), tube orbits and chaotic orbits. The box or box-like and tube orbits are represented by closed invariant curves on a Poincaré surface of section. The forms of the orbits and of the invariant curves can be explained by a third integral of motion I, that is given by the Giorgilli (1979) computer program. The nonresonant form of the third integral explains the box orbits, while a resonant form of this integral explains both the box orbits and the 1:1 tube orbits. The N-body model gives the distribution of velocities F, which is an exponential of the third integral.     

Spectral-domain covariance estimation with a priori knowledge

A knowledge-aided spectral-domain approach to estimating the interference covariance matrix used in space-time adaptive processing (STAP) is proposed. Prior knowledge of the range-Doppler clutter scene is used to identify geographic regions with homogeneous scattering statistics. Then, minimum-variance spectral estimation is used to arrive at a spectral-domain clutter estimate. Finally, space-time steering vectors are used to transform the spectral-domain estimate into a data-domain estimate of the clutter covariance matrix. The proposed technique is compared with ideal performance and to the fast maximum likelihood technique using simulated results. An investigation of the performance degradation that can occur due to various inaccurate knowledge assumptions is also presented