Extended Kalman filters using explicit and derivative-free local linearizations

We propose three variants of the extended Kalman filter (EKF) especially suited for parameter estimations in mechanical oscillators under Gaussian white noises. These filters are based on three versions of explicit and derivative-free local linearizations (DLL) of the non-linear drift terms in the governing stochastic differential equations (SDE-s). Besides a basic linearization of the non-linear drift functions via one-term replacements, linearizations using replacements through explicit Euler and Newmark expansions are also attempted in order to ensure higher closeness of true solutions with the linearized ones. Thus, unlike the conventional EKF, the proposed filters do not need computing derivatives (tangent matrices) at any stage. The measurements are synthetically generated by corrupting with noise the numerical solutions of the SDE-s through implicit versions of these linearizations. In order to demonstrate the effectiveness and accuracy of the proposed methods vis-à-vis the conventional EKF, numerical illustrations are provided for a few single degree-of-freedom (DOF) oscillators and a three-DOF shear frame with constant parameters.

A Joint Source-Channel Coding Scheme for Transmission of Correlated Discrete Sources over a Gaussian Multiple Access Channel

We consider the problem of transmission of correlated discrete alphabet sources over a Gaussian Multiple Access Channel (GMAC). A distributed bit-to-Gaussian mapping is proposed which yields jointly Gaussian codewords. This can guarantee lossless transmission or lossy transmission with given distortions, if possible. The technique can be extended to the system with side information at the encoders and decoder.

Two-Stage Extended Kalman Filters with Derivative-Free Local Linearizations

This paper proposes a derivative-free two-stage extended Kalman filter (2-EKF) especially suited for state and parameter identification of mechanical oscillators under Gaussian white noise. Two sources of modeling uncertainties are considered: (1) errors in linearization, and (2) an inadequate system model. The state vector is presently composed of the original dynamical/parameter states plus the so-called bias states accounting for the unmodeled dynamics. An extended Kalman estimation concept is applied within a framework predicated on explicit and derivative-free local linearizations (DLL) of nonlinear drift terms in the governing stochastic differential equations (SDEs). The original and bias states are estimated by two separate filters; the bias filter improves the estimates of the original states. Measurements are artificially generated by corrupting the numerical solutions of the SDEs with noise through an implicit form of a higher-order linearization. Numerical illustrations are provided for a few single- and multidegree-of-freedom nonlinear oscillators, demonstrating the remarkable promise that 2-EKF holds over its more conventional EKF-based counterparts. DOI: 10.1061/(ASCE)EM.1943-7889.0000255. (C) 2011 American Society of Civil Engineers.

On the Gaussian Many-to-One X Channel

In this paper, the Gaussian many-to-one X channel (XC), which is a special case of general multiuser XC, is studied. In the Gaussian many-to-one XC, communication links exist between all transmitters and one of the receivers, along with a communication link between each transmitter and its corresponding receiver. As per the XC assumption, transmission of messages is allowed on all the links of the channel. This communication model is different from the corresponding manyto- one interference channel (IC). Transmission strategies, which involve using Gaussian codebooks and treating interference from a subset of transmitters as noise, are formulated for the above channel. Sum-rate is used as the criterion of optimality for evaluating the strategies. Initially, a 3 x 3 many-to-one XC is considered and three transmission strategies are analyzed. The first two strategies are shown to achieve sum-rate capacity under certain channel conditions. For the third strategy, a sum-rate outer bound is derived and the gap between the outer bound and the achieved rate is characterized. These results are later extended to the K x K case. Next, a region in which the many-to-one XC can be operated as a many-to-one IC without the loss of sum-rate is identified. Furthermore, in the above region, it is shown that using Gaussian codebooks and treating interference as noise achieve a rate point that is within K/2 -1 bits from the sum-rate capacity. Subsequently, some implications of the above results to the Gaussian many-to-one IC are discussed. Transmission strategies for the many-to-one IC are formulated, and channel conditions under which the strategies achieve sum-rate capacity are obtained. A region where the sum-rate capacity can be characterized to within K/2 -1 bits is also identified. Finally, the regions where the derived channel conditions are satisfied for each strategy are illustrated for a 3 x 3 many-to-one XC and the corresponding many-to-one IC.

Twist phase in Gaussian-beam optics

The recently discovered twist phase is studied in the context of the full ten-parameter family of partially coherent general anisotropic Gaussian Schell-model beams. It is shown that the nonnegativity requirement on the cross-spectral density of the beam demands that the strength of the twist phase be bounded from above by the inverse of the transverse coherence area of the beam. The twist phase as a two-point function is shown to have the structure of the generalized Huygens kernel or Green's function of a first-order system. The ray-transfer matrix of this system is exhibited. Wolf-type coherent-mode decomposition of the twist phase is carried out. Imposition of the twist phase on an otherwise untwisted beam is shown to result in a linear transformation in the ray phase space of the Wigner distribution. Though this transformation preserves the four-dimensional phase-space volume, it is not symplectic and hence it can, when impressed on a Wigner distribution, push it out of the convex set of all bona fide Wigner distributions unless the original Wigner distribution was sufficiently deep into the interior of the set.

Gaussian Schell-model beams and general shape invariance

We present two six-parameter families of anisotropic Gaussian Schell-model beams that propagate in a shape-invariant manner, with the intensity distribution continuously twisting about the beam axis. The two families differ in the sense or helicity of this beam twist. The propagation characteristics of these shape-invariant beams are studied, and the restrictions on the beam parameters that arise from the optical uncertainty principle are brought out. Shape invariance is traced to a fundamental dynamical symmetry that underlies these beams. This symmetry is the product of spatial rotation and fractional Fourier transformation.

Three-dimensional 3d-4f heterometallic polymers containing both azide and carboxylate as co-ligands

The hydrothermal reaction of Ln(NO3)(3), Ni(NO3)(2), NaN3, and isonicotinic acid (L) yielded two novel 3-D coordination frameworks (1 and 2) of general formula [Ni(2)Ln(L)(5)(N-3)(2)(H2O)(3)] center dot 2H(2)O (Ln = Pr(III) for 1 and Nd(III) for 2), containing Ni-Pr or Ni-Nd hybrid extended three-dimensional networks containing both azido and carboxylate as co-ligands. Both the compounds are found to be isostructural and crystallize in monoclinic system having P2(1)/n space group. Here the lanthanide ions are found to be nonacoordinated. Both bidentate and monodentate modes of binding of the carboxylate with the lanthanides have been observed in the above complexes. Variable temperature magnetic studies of the above two complexes have been investigated in the temperature range 2-300 K which showed dominant antiferromagnetic interaction in both the cases and these experimental results are analyzed with the theoretical models. (c) 2008 Elsevier B.V. All rights reserved.

Anisotropic gaussian beams

Anisotropic gaussian beams are obtained as exact solutions to the parabolic wave equation. These beams have a quadratic phase front whose principal radii of curvature are non-degenerate everywhere. It is shown that, for the lowest order beams, there exists a plane normal to the beam axis where the intensity distribution is rotationally symmetric about the beam axis. A possible application of these beams as normal modes of laser cavities with astigmatic mirrors is noted.

Exact evaluation of Gaussian path integrals involving non-local potentials

An exact expression for the calculation of gaussian path integrals involving non-local potentials is given. Its utility is demonstrated by using it to evaluate a path integral arising in the study of an electron gas in a random potential.

Invited paper Application of the least mean fourth (LMF) adaptive algorithm to signals associated with gaussian noise

This paper considers the applicability of the least mean fourth (LM F) power gradient adaptation criteria with 'advantage' for signals associated with gaussian noise, the associated noise power estimate not being known. The proposed method, as an adaptive spectral estimator, is found to provide superior performance than the least mean square (LMS) adaptation for the same (or even lower) speed of convergence for signals having sufficiently high signal-to-gaussian noise ratio. The results include comparison of the performance of the LMS-tapped delay line, LMF-tapped delay line, LMS-lattice and LMF-lattice algorithms, with the Burg's block data method as reference. The signals, like sinusoids with noise and stochastic signals like EEG, are considered in this study.

Validation-Based Sparse Gaussian Process Classifier Design

Gaussian processes (GPs) are promising Bayesian methods for classification and regression problems. Design of a GP classifier and making predictions using it is, however, computationally demanding, especially when the training set size is large. Sparse GP classifiers are known to overcome this limitation. In this letter, we propose and study a validation-based method for sparse GP classifier design. The proposed method uses a negative log predictive (NLP) loss measure, which is easy to compute for GP models. We use this measure for both basis vector selection and hyperparameter adaptation. The experimental results on several real-world benchmark data sets show better orcomparable generalization performance over existing methods.

Analysis of randomly time varying systems by gaussian closure technique

The Gaussian probability closure technique is applied to study the random response of multidegree of freedom stochastically time varying systems under non-Gaussian excitations. Under the assumption that the response, the coefficient and the excitation processes are jointly Gaussian, deterministic equations are derived for the first two response moments. It is further shown that this technique leads to the best Gaussian estimate in a minimum mean square error sense. An example problem is solved which demonstrates the capability of this technique for handling non-linearity, stochastic system parameters and amplitude limited responses in a unified manner. Numerical results obtained through the Gaussian closure technique compare well with the exact solutions.

Gaussian-Maxwell beams

Gaussian-beam-type solutions to the Maxwell equations are constructed by using results from relativistic front analysis, and the propagation characteristics of these beams are analyzed. The rays of geometrical optics are shown to be the trajectories of energy flow, as given by the Poynting vector. The longitudinal components of the field vectors in the direction of the beam axis, though small, are shown to be essential for a consistent description.

Anisotropic Gaussian Schell-model beams: Passage through optical systems and associated invariants

Anisotropic Gaussian Schell-model (AGSM) fields and their transformation by first-order optical systems (FOS’s) forming Sp(4,R) are studied using the generalized pencils of rays. The fact that Sp(4,R), rather than the larger group SL(4,R), is the relevant group is emphasized. A convenient geometrical picture wherein AGSM fields and FOS’s are represented, respectively, by antisymmetric second-rank tensors and de Sitter transformations in a (3+2)-dimensional space is developed. These fields are shown to separate into two qualitatively different families of orbits and the invariants over each orbit, two in number, are worked out. We also develop another geometrical picture in a (2+1)-dimensional Minkowski space suitable for the description of the action of axially symmetric FOS’s on AGSM fields, and the invariants, now seven in number, are derived. Interesting limiting cases forming coherent and quasihomogeneous fields are analyzed.