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.     

Matched Filter Criteria and Range Weighting for Weather Radar

The effective weighting function for weather radar is defined. This weighting function considers the effects of both the transmitted signal and the receiver filter. It is used to assign effective ranges to samples taken at prescribed times. For uniformly distributed targets it is shown that "signal"-to-noise ratio depends on the receiver filter, transmitted signal envelope, and receiver noise spectral density. Maximization of this signal-to-noise ratio when range resolution constraints are imposed is discussed, and a receiver design approach specifically adapted for Doppler weather radars is developed.     

Calibration of the MERIC full-polarimetric radar: theory and implementation

This paper describes the architecture and calibration design of the experimental ground based radar station MERIC. This full-polarimetric radar is conceived for the analysis and the recognition of non cooperative aircraft in flight.We carefully study how the full-polarimetric capability is obtained for a simultaneous transmit (simultaneous transmission of two linear FM with opposite slopes) radar system, using analogue deramping with a replica.The phase distortions of the signal propagating in the four polarimetric channels are carefully estimated. We define a phase calibrating method compatible with the outdoor measurements conditions and few constraints on reference targets.We show the phase accuracy obtained with the proposed calibration method on real measurements.     

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.     

An Optimum Phase Reference Detector for Fully Modulated Phase-Shift Keyed Signals

Many modern telemetry systems which use phase-shift keying (PSK) have receivers which derive a coherent reference from the fully modulated PSK signal itself and thus conserve the energy which otherwise would be allocated to a discrete reference signal. In this paper, an optimum receiver structure for estimating a phase reference from the PSK signal itself is derived and its realization discussed. It is shown that at low signal-to-noise ratios, the optimum detector can be realized with a Costas loop. Since a Costas loop and squaring loop exhibit identical performance, it follows that either of these simple devices gives optimum performance for low-input signal-to-noise ratios.     

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.     

A Possible Transfer Mechanism for the 11-Year Solar Cycle to the Lower Stratosphere

Observational evidence of the 11-year solar cycle (SC) modulation of stratosphere temperatures and winds from the ERA-40 dataset is reviewed, with emphasis on the Northern winter hemisphere. A frequency modulation of sudden warming events is noted, with warmings occurring earlier in solar minimum periods than in solar maximum periods. The observed interaction between the influence of the SC and the quasi biennial oscillation (QBO) on the frequency of sudden warmings is noted as a possible clue for understanding their mechanism of influence. A possible transfer route for the 11-year solar cycle from the equatorial stratopause region to the lowest part of the stratosphere is proposed, via an influence on sudden warming events and the associated induced meridional circulation. SC and QBO composites of zonal wind anomalies show anomalous wind distributions in the subtropical upper stratosphere in early winter. Mechanistic model experiments are reviewed that demonstrate a sensitivity of sudden warmings to small wind anomalies in this region. Various diagnostics from these experiments are shown, including EP fluxes and their divergence and also the synoptic evolution of the polar vortex, in order to understand the mechanism of the influence. Some recent GCM experiments to investigate the SC/QBO interaction are also described. They simulate reasonably well the observed SC/QBO interaction of sudden warming events and appear to support the hypothesis that tropical/subtropical upper stratospheric wind anomalies are an important influence on the timing of sudden warmings.     

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.