Analysis of Multimodal Systems

The estimation of a multimodal linear system whose mode-to-mode transitions are described by a finite-state Markov chain is described. The problem has application in studying separation standards in an air traffic control environment. An optimal solution is formulated which is computationally infeasible. A suboptimal estimator is then derived which closely approximates the optimal estimator. An example is presented to illustrate the technique.     

基于CFD技术的鼻腔输药的流动特性研究

鼻腔输药作为一种新颖的输药方式,具有方便可靠的特点,不仅可以适用于局部而且适用于全身用药,已经受到广泛的重视。但由于人类的呼吸道复杂的立体结构使得鼻腔输药成为生物医药研究的一个长期的挑战。利用CFD(计算流体动力学,Computational Fluid Dynamics)技术对鼻腔输药的流动特性进行研究,对一个基于真实的人类鼻腔解剖的3D模型的输药的流动过程进行了分析,着重研究了不同的药物颗粒的流动和沉淀特性。利用CFD的商用软件FLUENT在层流情况下对不同的密度、直径和流动速率的药物颗粒在鼻腔中流动进行了数值模拟,并对计算结果进行了比较分析。CFD研究的结果为鼻腔输药这种新颖的药物输送方式的优化设计提供了依据。     

Probing the era of galaxy formation via TeV gamma ray absorption by the near infrared extragalactic background

We present models of the extragalactic background light (EBL) based on several scenarios of galaxy formation and evolution. We have treated galaxy formation with the Press-Schecter approximation for both cold dark matter (CDM) and cold+hot dark matter (CHDM) models, representing a moderate (z _f 3) and a late (z _f 1) era of galaxy formation respectively. Galaxy evolution has been treated by considering a variety of stellar types, different initial mass functions and star formation histories, and with an accounting of dust absorption and emission. We find that the dominant factor influencing the EBL is the epoch of galaxy formation. A recently proposed method for observing the EBL utilizing the absorption of 0.1 to 10 TeV gamma-rays from active galactic nuclei (AGN) is shown to be capable of discriminating between different galaxy formation epochs. The one AGN viewed in TeV light, Mrk 421, does show some evidence for a cutoff above 3 TeV; based on the EBL models presented here, we suggest that this is due to extinction in the source. The large absorption predicted at energies > 200 GeV for sources at z > 0.5 indicates that observations of TeV gamma-ray bursts (GRB) would constrain or eliminate models in which the GRB sources lie at cosmological distances.Now at University of Chicago, Dept. of Astronomy & Astrophysics.     

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.     

Specification and Performance of Passive Sonar Spectral Estimators

The application of existing estimation theory to the problem of specification and performance of passive sonar spectral estimators is considered. The classification function is addressed, so that the signal is assumed to be present, and so that the energy arrival angle is known. The spatial filter considered is a line array of M equally spaced omnidirectional hydrophones. Signal and ambient noise are both zero-mean, wide-sense, stationary Gaussian random processes that differ in their spatial correlation across the face of the array. The signal is a plane wave that can be made totally spacially corrected between array elements by inserting delays between sensors to invert the signal propagation delay. The noise correlation is a function of frequency, bandwidth, element separation, and the relative time delay between sensors. Under these assumptions, the Cramer-Rao lower bound is derived for the class of unbiased estimates of signal power in a narrow frequency band at the hydrophone in the presence of correlated ambient noise of known power. The bound is examined numerically, resulting in a threshold phenomenon with M that constitutes a new design consideration. In addition, there is a striking insensitivity to realistic values of ambient noise correlation, and there are ranges in signal-to-noise ratio for which one gains more by increasing M than by increasing the bandwidth-time product. Specific processors, including a new unbiased estimator when noise power is unknown, are developed.     

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 gentle process for the formation of Algols

This work is concerned with binary systems that we call ‘moderately close’. These are systems in which the primary (by which we mean the initially more massive star) fills its Roche lobe when it is on the giant branch with a deep convective envelope but before helium ignition (late case B). We find that if the mass ratio q(= M ₁/M ₂) < q _crit = 0.7 when the primary fills its Roche lobe positive feedback will lead to a rapid hydrodynamic phase of mass transfer which will probably lead to common envelope evolution and thence to either coalescence or possibly to a close binary in a planetary nebula. Although most Algols have probably filled their Roche lobes before evolving off the main-sequence we find that some could not have and are therefore ‘moderately close’. Since rapid overflow is unlikely to lead to an Algol-like system there must be some way of avoiding it. The most likely possibility is that the primary can lose sufficient mass to reduce q below q _crit before overflow begins. Ordinary mass loss rates are insufficient but evidence that enhanced mass loss does take place is provided by RS CVn systems that have inverted mass ratios but have not yet begun mass transfer. We postulate that the cause of enhanced mass loss lies in the heating of the corona by by magnetic fields maintained by an α — ω dynamo which is enhanced by tidal effects associated with corotation. In order to model the the effects of enhanced mass loss we ignore the details and adopt an empirical approach calibrating a simple formula with the RS CVn system Z Her. Using further empirical relations (deduced from detailed stellar models) that describe the evolution of red giants we have investigated the effect on a large number of systems of various initial mass ratios and periods. These are notable in that some systems can now enter a much gentler Algol-like overflow phase and others are prevented from transferring mass altogether. We have also investigated the effects of enhanced angular momentum loss induced by corotation of the wind in the strong magnetic fields and consider this in relation to observed period changes. We find that a typical ‘moderately close’ Algol-like system evolves through an RS CVn like system and then possibly a symbiotic state before becoming an Algol and then goes on through a red giant-white dwarf state which may become symbiotic before ending up as a double white dwarf system in either a close or wide orbit depending on how much mass is lost before the secondary fills its Roche lobe.     

Radar Identification of Naval Vessels

A scheme is presented for the identification of naval vessels via active multiple-frequency radar interrogation. A major virtue of the described method is the use of a response waveform synthesized using amplitude data only. A prediction correlation using natural resonances associated with substructures on the vessels is applied to synthetically generated matched-filter response waveforms. The identification scheme is tested using measured model data for 8 vessels on a simulated sea surface. A correct identification probability of roughly 77 percent is obtained.     

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.