Extended solution to multiple maneuvering target tracking

An improved algorithm for tracking multiple maneuvering targets is presented. This approach is implemented with an approximate adaptive filter consisting of the one-step conditional maximum-likelihood technique together with the extended Kalman filter and an adaptive maneuvering compensator. In order to avoid the extra computational burden of considering events with negligible probability, a validation matrix is defined in the tracking structure. With this approach, data-association and target maneuvering problems can be solved simultaneously. Detailed Monte Carlo simulations of the algorithm for many tracking situations are described. Computer simulation results indicate that this approach successfully tracks multiple maneuvering targets over a wide range of conditions     

The effects of space flight on some rat liver enzymes regulating carbohydrate and lipid metabolism

We have examined, in the livers of rats carried aboard the Cosmos 936 biosatellite, the activities of about 30 enzymes concerned with carbohydrate and lipid metabolism. In addition to the enzyme studies, the levels of glycogen and of the individual fatty acids in hepatic lipids were determined. Livers from flight and ground control rats at recovery (R₀) and 25 days after recovery (R₂₅) were used for these analyses.

For all parameters measured, the most meaningful comparisons are those made between flight stationary (FS) and flight centrifuged (FC) animals at R₀. When these two groups of flight rats were compared at R₀, statistically significant decreases in the activity levels of glycogen phosphorylase, -glycerol phosphate acyl transferase, diglyceride acyl transferase, aconitase and 6-phosphogluconate dehydrogenase and an increase in the palmitoyl CoA desaturase were noted in the weightless group (FS). The significance of these findings was strengthened by the fact that all enzyme activities showing alterations at R₀ returned to normal 25 days postflight. When liver glycogen and total fatty acids of the two sets of flight animals were determined, significant differences that could be attributed to reduced gravity were observed. The weightless group (FS) at R₀ contained, on the average, more than twice the amount of glycogen than did the centrifuged controls (FC) and a remarkable shift in the ratio of palmitate to palmitoleate was noted. These metabolic alterations, both in enzyme levels and in hepatic constituents, appear to be characteristic of the weightless condition. Our data seem to justify the conclusion that centrifugation during flight is equivalent to terrestrial gravity.     

Geolocation of a known altitude object from TDOA and FDOA measurements

Most satellite systems for locating an object on Earth use only time difference of arrival (TDOA) measurements. When there are relative motions between an emitter and receivers, frequency difference of arrival (FDOA) measurements can be used as well. Often, the altitude of an object is known (it is zero, for example) or can be measured with an altimeter. Two sets of geolocation solutions are proposed which exploit the altitude constraint to improve the localization accuracy. One is for TDOAs alone and the other for the combination of TDOA and FDOA measurements. The additional complexity by imposing the constraint is a one-dimensional Newton's search and the rooting of a polynomial. The covariance matrices of the new estimators are derived under a small measurement noise assumption and shown to attain the constrained Cramer-Rao lower bound (CRLB). When there is a bias error in the assumed altitude, using the altitude constraint will introduce a bias to the solution. Since applying the constraint decreases the variance, there is a tradeoff between variance and bias in the mean square error (MSE). The maximum allowable altitude error such that the constraint solution will remain superior to the unconstraint is given. Simulation results are included to corroborate the theoretical development.     

一维非定常流的气体格子模型

本文研究气体格子模型的一维流动。阐述一般的简单气体格子族。利用2-速度系数模型研究了激波的形成、反射以及相互干扰。计入了不同的初始条件和边界条件。一维气体格子模型考虑了可压缩性和二维、三维模型的基本特性。     

The Solar Origin of Corotating Interaction Regions and Their Formation in the Inner Heliosphere

Corotating Interaction Regions (CIRs) form as a consequence of the compression of the solar wind at the interface between fast speed streams and slow streams. Dynamic interaction of solar wind streams is a general feature of the heliospheric medium; when the sources of the solar wind streams are relatively stable, the interaction regions form a pattern which corotates with the Sun. The regions of origin of the high speed solar wind streams have been clearly identified as the coronal holes with their open magnetic field structures. The origin of the slow speed solar wind is less clear; slow streams may well originate from a range of coronal configurations adjacent to, or above magnetically closed structures. This article addresses the coronal origin of the stable pattern of solar wind streams which leads to the formation of CIRs. In particular, coronal models based on photospheric measurements are reviewed; we also examine the observations of kinematic and compositional solar wind features at 1 AU, their appearance in the stream interfaces (SIs) of CIRs, and their relationship to the structure of the solar surface and the inner corona; finally we summarise the Helios observations in the inner heliosphere of CIRs and their precursors to give a link between the optical observations on their solar origin and the in-situ plasma observations at 1 AU after their formation. The most important question that remains to be answered concerning the solar origin of CIRs is related to the origin and morphology of the slow solar wind. This revised version was published online in August 2006 with corrections to the Cover Date.     

Corotating Interaction Regions at High Latitudes

Ulysses observed a stable strong CIR from early 1992 through 1994 during its first journey into the southern hemisphere. After the rapid latitude scan in early 1995, Ulysses observed a weaker CIR from early 1996 to mid-1997 in the northern hemisphere as it traveled back to the ecliptic at the orbit of Jupiter. These two CIRs are the observational basis of the investigation into the latitudinal structure of CIRs. The first CIR was caused by an extension of the northern coronal hole into the southern hemisphere during declining solar activity, whereas the second CIR near solar minimum activity was caused by small warps in the streamer belt. The latitudinal structure is described through the presentation of three 26-day periods during the southern CIR. The first at ∼24°S shows the full plasma interaction region including fast and slow wind streams, the compressed shocked flows with embedded stream interface and heliospheric current sheet (HCS), and the forward and reverse shocks with associated accelerated ions and electrons. The second at 40°S exhibits only the reverse shock, accelerated particles, and the 26-day modulation of cosmic rays. The third at 60°S shows only the accelerated particles and modulated cosmic rays. The possible mechanisms for the access of the accelerated particles and the CIR-modulated cosmic rays to high latitudes above the plasma interaction region are presented. They include direct magnetic field connection across latitude due to stochastic field line weaving or to systematic weaving caused by solar differential rotation combined with non-radial expansion of the fast wind. Another possible mechanism is particle diffusion across the average magnetic field, which includes stochastic field line weaving. A constraint on connection to a distant portion of the CIR is energy loss in the solar wind, which is substantial for the relatively slow-moving accelerated ions. Finally, the weaker northern CIR is compared with the southern CIR. It is weak because the inclination of the streamer belt and HCS decreased as Ulysses traveled to lower latitudes so that the spacecraft remained at about the maximum latitudinal extent of the HCS. This revised version was published online in August 2006 with corrections to the Cover Date.     

Comparison of various periodograms for sinusoid detection andfrequency estimation

With the advent of the fast Fourier transform (FFT) algorithm, the periodogram and its variants such as the Bartlett's procedure and Welch method, have become very popular for spectral analysis. However, there has not been a thorough comparison of the detection and estimation performances of these methods. Different forms of the periodogram are studied here for single real tone detection and frequency estimation in the presence of white Gaussian noise. The threshold effect in frequency estimation, that is, when the estimation errors become several orders of magnitude greater than the Cramer-Rao lower bound (CRLB), is also investigated. It is shown that the standard periodogram gives the optimum detection performance for a pure tone while the Welch method is the best detector when there is phase instability in the sinusoid. As expected, since the conventional periodogram is a maximum likelihood estimator of frequency, it generally provides the minimum mean square frequency estimation errors     

Dynamically adaptable m-best 2-D assignment algorithm andmultilevel parallelization

In recent years, there has been considerable interest within the tracking community in an approach to data association based on the m-best two-dimensional (2D) assignment algorithm. Much of the interest has been spurred by its ability to provide various efficient data association solutions, including joint probabilistic data association (JPDA) and multiple hypothesis tracking (MHT). The focus of this work is to describe several recent improvements to the m-best 2D assignment algorithm. One improvement is to utilize a nonintrusive 2D assignment algorithm switching mechanism, based on a problem sparsity threshold. Dynamic switching between two different 2D assignment algorithms, highly suited for sparse and dense problems, respectively, enables more efficient solutions to the numerous 2D assignment problems generated in the m-best 2D assignment framework. Another improvement is to utilize a multilevel parallelization enabling many independent and highly parallelizable tasks to be executed concurrently, including 1) solving the multiple 2D assignment problems via a parallelization of the m-best partitioning task, and 2) calculating the numerous gating tests, state estimates, covariance calculations, and likelihood function evaluations (used as cost coefficients in the 2D assignment problem) via a parallelization of the data association interface task. Using both simulated data and an air traffic surveillance (ATS) problem based on data from two Federal Aviation Administration (FAA) air traffic control radars, we demonstrate that efficient solutions to the data association problem are obtainable using our improvements in the m-best 2D assignment algorithm     

Evaluation of gravity-dependent membrane potential shift in Paramecium

It is still debated whether or not gravity can stimulate unicellular organisms. This question may be settled by revealing changes in the membrane potential in a manner depending on the gravitational forces imposed on the cell. We estimated the gravity-dependent membrane potential shift to be about 1 mV G⁻¹ for Paramecium showing gravikinesis at 1–5 G, on the basis of measurements of gravity-induced changes in active propulsion and those of propulsive velocity in solutions, in which the membrane potential has been measured electrophysiologically. The shift in membrane potential to this extent may occur from mechanoreceptive changes in K⁺ or Ca²⁺ conductance by about 1% and might be at the limit of electrophysiological measurement using membrane potential-sensitive dyes. Our measurements of propulsive velocity vs membrane potential also suggested that the reported propulsive force of Paramecium measured in a solution of graded densities with the aid of a video centrifuge microscope at 350 G was 11 times as large as that for −29 mV, i.e., the resting membrane potential at [K⁺]_o = 1 mM and [Ca²⁺]_o = 1 mM, and, by extrapolation, that Paramecium was hyperpolarized to −60 mV by gravity stimulation of 100- G equivalent, the value corrected by considering the reduction of density difference between the interior and exterior of the cell in the graded density solution. The estimated shift of the membrane potential from −29 mV to −60 mV by 100- G equivalent stimulation, i.e., 0.3 mV G⁻¹, could reach the magnitude entirely feasible to be measured more directly.     

Modelling solar irradiance variations: Comparison with observations,including line-ratio variations

Solar irradiance variations show a strong temporal and spectral dependence. The progression of the Sun through its activity cycle as well as solar rotation are mirrored in the irradiance variations. The spectral dependence is such that the variations are several magnitudes larger in the EUV than in the visible or infrared.We present a simple 3-component model that is based on the assumption that changes in the solar flux are exclusively due to changes in spot and facular coverage. We compare our model to observations of the spectral solar irradiance variations.Despite its simplicity, we find that the agreement between our model and the observations is surprisingly good. We also explore the reliability and the limitations of our approach by comparing observations of the solar facular contrast and of the changes in spectral line depths with our calculations.