Bearing-only SLAM : a vision-based navigation system for autonomous robots

To navigate successfully in a previously unexplored environment, a mobile robot must be able to estimate the spatial relationships of the objects of interest accurately. A Simultaneous Localization and Mapping (SLAM) sys- tem employs its sensors to build incrementally a map of its surroundings and to localize itself in the map simultaneously. The aim of this research project is to develop a SLAM system suitable for self propelled household lawnmowers. The proposed bearing-only SLAM system requires only an omnidirec- tional camera and some inexpensive landmarks. The main advantage of an omnidirectional camera is the panoramic view of all the landmarks in the scene. Placing landmarks in a lawn field to define the working domain is much easier and more flexible than installing the perimeter wire required by existing autonomous lawnmowers. The common approach of existing bearing-only SLAM methods relies on a motion model for predicting the robot’s pose and a sensor model for updating the pose. In the motion model, the error on the estimates of object positions is cumulated due mainly to the wheel slippage. Quantifying accu- rately the uncertainty of object positions is a fundamental requirement. In bearing-only SLAM, the Probability Density Function (PDF) of landmark position should be uniform along the observed bearing. Existing methods that approximate the PDF with a Gaussian estimation do not satisfy this uniformity requirement. This thesis introduces both geometric and proba- bilistic methods to address the above problems. The main novel contribu- tions of this thesis are: 1. A bearing-only SLAM method not requiring odometry. The proposed method relies solely on the sensor model (landmark bearings only) without relying on the motion model (odometry). The uncertainty of the estimated landmark positions depends on the vision error only, instead of the combination of both odometry and vision errors. 2. The transformation of the spatial uncertainty of objects. This thesis introduces a novel method for translating the spatial un- certainty of objects estimated from a moving frame attached to the robot into the global frame attached to the static landmarks in the environment. 3. The characterization of an improved PDF for representing landmark position in bearing-only SLAM. The proposed PDF is expressed in polar coordinates, and the marginal probability on range is constrained to be uniform. Compared to the PDF estimated from a mixture of Gaussians, the PDF developed here has far fewer parameters and can be easily adopted in a probabilistic framework, such as a particle filtering system. The main advantages of our proposed bearing-only SLAM system are its lower production cost and flexibility of use. The proposed system can be adopted in other domestic robots as well, such as vacuum cleaners or robotic toys when terrain is essentially 2D.

Fundamental solution and discrete random walk model for time-space fractional diffusion equation

In this paper, we consider a time-space fractional diffusion equation of distributed order (TSFDEDO). The TSFDEDO is obtained from the standard advection-dispersion equation by replacing the first-order time derivative by the Caputo fractional derivative of order α∈(0,1], the first-order and second-order space derivatives by the Riesz fractional derivatives of orders β 1∈(0,1) and β 2∈(1,2], respectively. We derive the fundamental solution for the TSFDEDO with an initial condition (TSFDEDO-IC). The fundamental solution can be interpreted as a spatial probability density function evolving in time. We also investigate a discrete random walk model based on an explicit finite difference approximation for the TSFDEDO-IC.

The statistics of refractive error maps : managing wavefront aberration analysis without Zernike polynomials

The refractive error of a human eye varies across the pupil and therefore may be treated as a random variable. The probability distribution of this random variable provides a means for assessing the main refractive properties of the eye without the necessity of traditional functional representation of wavefront aberrations. To demonstrate this approach, the statistical properties of refractive error maps are investigated. Closed-form expressions are derived for the probability density function (PDF) and its statistical moments for the general case of rotationally-symmetric aberrations. A closed-form expression for a PDF for a general non-rotationally symmetric wavefront aberration is difficult to derive. However, for specific cases, such as astigmatism, a closed-form expression of the PDF can be obtained. Further, interpretation of the distribution of the refractive error map as well as its moments is provided for a range of wavefront aberrations measured in real eyes. These are evaluated using a kernel density and sample moments estimators. It is concluded that the refractive error domain allows non-functional analysis of wavefront aberrations based on simple statistics in the form of its sample moments. Clinicians may find this approach to wavefront analysis easier to interpret due to the clinical familiarity and intuitive appeal of refractive error maps.

Interleaved PC-OFDM to reduce peak-to-average power ratio

Parallel combinatory orthogonal frequency division multiplexing (PC-OFDM yields lower maximum peak-to-average power ratio (PAR), high bandwidth efficiency and lower bit error rate (BER) on Gaussian channels compared to OFDM systems. However, PC-OFDM does not improve the statistics of PAR significantly. In this chapter, the use of a set of fixed permutations to improve the statistics of the PAR of a PC-OFDM signal is presented. For this technique, interleavers are used to produce K-1 permuted sequences from the same information sequence. The sequence with the lowest PAR, among K sequences is chosen for the transmission. The PAR of a PC-OFDM signal can be further reduced by 3-4 dB by this technique. Mathematical expressions for the complementary cumulative density function (CCDF)of PAR of PC-OFDM signal and interleaved PC-OFDM signal are also presented.

Statistical modeling to determine resonant frequency information present in combustion chamber pressure signals

A statistical modeling method to accurately determine combustion chamber resonance is proposed and demonstrated. This method utilises Markov-chain Monte Carlo (MCMC) through the use of the Metropolis-Hastings (MH) algorithm to yield a probability density function for the combustion chamber frequency and find the best estimate of the resonant frequency, along with uncertainty. The accurate determination of combustion chamber resonance is then used to investigate various engine phenomena, with appropriate uncertainty, for a range of engine cycles. It is shown that, when operating on various ethanol/diesel fuel combinations, a 20% substitution yields the least amount of inter-cycle variability, in relation to combustion chamber resonance.

Computed success rates of various carrier phase integer estimation solutions and their comparison with statistical success rates

The success rate of carrier phase ambiguity resolution (AR) is the probability that the ambiguities are successfully fixed to their correct integer values. In existing works, an exact success rate formula for integer bootstrapping estimator has been used as a sharp lower bound for the integer least squares (ILS) success rate. Rigorous computation of success rate for the more general ILS solutions has been considered difficult, because of complexity of the ILS ambiguity pull-in region and computational load of the integration of the multivariate probability density function. Contributions of this work are twofold. First, the pull-in region mathematically expressed as the vertices of a polyhedron is represented by a multi-dimensional grid, at which the cumulative probability can be integrated with the multivariate normal cumulative density function (mvncdf) available in Matlab. The bivariate case is studied where the pull-region is usually defined as a hexagon and the probability is easily obtained using mvncdf at all the grid points within the convex polygon. Second, the paper compares the computed integer rounding and integer bootstrapping success rates, lower and upper bounds of the ILS success rates to the actual ILS AR success rates obtained from a 24 h GPS data set for a 21 km baseline. The results demonstrate that the upper bound probability of the ILS AR probability given in the existing literatures agrees with the actual ILS success rate well, although the success rate computed with integer bootstrapping method is a quite sharp approximation to the actual ILS success rate. The results also show that variations or uncertainty of the unit–weight variance estimates from epoch to epoch will affect the computed success rates from different methods significantly, thus deserving more attentions in order to obtain useful success probability predictions.

Single crystal Raman spectroscopy of brandholzite Mg[Sb(OH)6]2•6H2O and bottinoite Ni[Sb(OH)6]2•6H2O and the polycrystalline Raman spectrum of mopungite Na[Sb(OH)6]

The single crystal Raman spectra of minerals brandholzite and bottinoite, formula M[Sb(OH)6]2•6H2O, where M is Mg+2 and Ni+2 respectively, and the non-aligned Raman spectrum of mopungite, formula Na[Sb(OH)6], are presented for the first time. The mixed metal minerals comprise of alternating layers of [Sb(OH)6]-1 octahedra and mixed [M(H2O)6]+2 / [Sb(OH)6]-1 octahedra. Mopungite comprises hydrogen bonded layers of [Sb(OH)6]-1 octahedra linked within the layer by Na+ ions. The spectra of the three minerals were dominated by the Sb-O symmetric stretch of the [Sb(OH)6]-1 octahedron, which occurs at approximately 620 cm-1. The Raman spectrum of mopungite showed many similarities to spectra of the di-octahedral minerals informing the view that the Sb octahedra gave rise to most of the Raman bands observed, particularly below 1200 cm-1. Assignments have been proposed based on the spectral comparison between the minerals, prior literature and density field theory calculations of the vibrational spectra of the free [Sb(OH)6]-1 and [M(H2O)6]+2 octahedra by a model chemistry of B3LYP/6-31G(d) and lanl2dz for the Sb atom. The single crystal data spectra showed good mode separation, allowing the majority of the bands to be assigned a symmetry species of A or E.

Seeing the Wood for the Trees: A Critical Evaluation of Methods to Estimate the Parameters of Stochastic Differential Equations

Maximum-likelihood estimates of the parameters of stochastic differential equations are consistent and asymptotically efficient, but unfortunately difficult to obtain if a closed-form expression for the transitional probability density function of the process is not available. As a result, a large number of competing estimation procedures have been proposed. This article provides a critical evaluation of the various estimation techniques. Special attention is given to the ease of implementation and comparative performance of the procedures when estimating the parameters of the Cox–Ingersoll–Ross and Ornstein–Uhlenbeck equations respectively.

Stochastic modelling of T cell homeostasis for two competing clonotypes via the master equation

Stochastic models for competing clonotypes of T cells by multivariate, continuous-time, discrete state, Markov processes have been proposed in the literature by Stirk, Molina-París and van den Berg (2008). A stochastic modelling framework is important because of rare events associated with small populations of some critical cell types. Usually, computational methods for these problems employ a trajectory-based approach, based on Monte Carlo simulation. This is partly because the complementary, probability density function (PDF) approaches can be expensive but here we describe some efficient PDF approaches by directly solving the governing equations, known as the Master Equation. These computations are made very efficient through an approximation of the state space by the Finite State Projection and through the use of Krylov subspace methods when evolving the matrix exponential. These computational methods allow us to explore the evolution of the PDFs associated with these stochastic models, and bimodal distributions arise in some parameter regimes. Time-dependent propensities naturally arise in immunological processes due to, for example, age-dependent effects. Incorporating time-dependent propensities into the framework of the Master Equation significantly complicates the corresponding computational methods but here we describe an efficient approach via Magnus formulas. Although this contribution focuses on the example of competing clonotypes, the general principles are relevant to multivariate Markov processes and provide fundamental techniques for computational immunology.

N-Aryl stilbazolium dyes as sensitizers for solar cells

Eight new N-arylstilbazolium chromophores with electron donating –NR2 (R = Me or Ph) substituents have been synthesized via Knoevenagel condensations and isolated as their PF6− salts. These compounds have been characterized by using various techniques including 1H NMR and IR spectroscopies and electrospray mass spectrometry. UV–vis absorption spectra recorded in acetonitrile are dominated by intense, low energy π → π* intramolecular charge-transfer (ICT) bands, and replacing Me with Ph increases the ICT energies. Cyclic voltammetric studies show irreversible reduction processes, together with oxidation waves that are irreversible for R = Me, but reversible for R = Ph. Single crystal X-ray structures have been determined for three of the methyl ester-substituted stilbazolium salts and for the Cl− salts of their picolinium precursors. Time-dependent density functional theory calculations afford reasonable predictions of ICT energies, but greater rigour is necessary for –NPh2 derivatives. The four new acid-functionalized dyes give moderate sensitization efficiencies (ca. 0.2%) when using TiO2-based photoanodes, with relatively higher values for R = Ph vs Me, while larger efficiencies (up to 0.8%) are achieved with ZnO substrates.

Experimental and theoretical investigation of ligand effects on the synthesis of ZnO nanoparticles

ZnO nanoparticles with highly controllable particle sizes(less than 10 nm) were synthesized using organic capping ligands in Zn(Ac)2 ethanolic solution. The molecular structure of the ligands was found to have significant influence on the particle size. The multi-functional molecule tris(hydroxymethyl)-aminomethane (THMA) favoured smaller particle distributions compared with ligands possessing long hydrocarbon chains that are more frequently employed. The adsorption of capping ligands on ZnnOn crystal nuclei (where n = 4 or 18 molecular clusters of(0001) ZnO surfaces) was modelled by ab initio methods at the density functional theory (DFT) level. For the molecules examined, chemisorption proceeded via the formation of Zn...O, Zn...N, or Zn...S chemical bonds between the ligands and active Zn2+ sites on ZnO surfaces. The DFT results indicated that THMA binds more strongly to the ZnO surface than other ligands, suggesting that this molecule is very effective at stabilizing ZnO nanoparticle surfaces. This study, therefore, provides new insight into the correlation between the molecular structure of capping ligands and the morphology of metal oxide nanostructures formed in their presence.

Vapour phase assembly of a halogen bonded complex of an isoindoline nitroxide and 1,2-diiodotetrafluorobenzene

Vapour phase assembly has been used for the first time to prepare co-crystals in which the primary intermolecular interaction is halogen bonding. Co-crystals of the nitroxide 1,1,3,3-tetramethylisoindolin-2-yloxyl (TMIO) and 1,2-diiodotetrafluorobenzene (1,2-DITFB) are readily formed under standard sublimation conditions. Single crystal X-ray diffraction confirmed the structure of a 2:2 cyclic tetramer, (TMIO)2·(1,2-DITFB)2, which exhibits a new halogen bonding motif, with each nitroxide oxygen atom accepting two halogen bonds. Powder X-ray diffraction confirmed the homogeneity of the bulk sample. The crystalline complex was further characterized in the solid state using thermal analysis and vibrational spectroscopy (infrared and Raman). Density functional theory calculations were also used to evaluate the enthalpy of formation, electrostatic potential and unpaired electron density of the complex. These findings illustrate the preparation of co-crystals where solution state methodology is problematic and the potential of this approach for the formation of novel organic spin systems.

Utilising Reliability and Condition Monitoring Data for Asset Health Prognosis

The ability to forecast machinery health is vital to reducing maintenance costs, operation downtime and safety hazards. Recent advances in condition monitoring technologies have given rise to a number of prognostic models which attempt to forecast machinery health based on condition data such as vibration measurements. This paper demonstrates how the population characteristics and condition monitoring data (both complete and suspended) of historical items can be integrated for training an intelligent agent to predict asset health multiple steps ahead. The model consists of a feed-forward neural network whose training targets are asset survival probabilities estimated using a variation of the Kaplan–Meier estimator and a degradation-based failure probability density function estimator. The trained network is capable of estimating the future survival probabilities when a series of asset condition readings are inputted. The output survival probabilities collectively form an estimated survival curve. Pump data from a pulp and paper mill were used for model validation and comparison. The results indicate that the proposed model can predict more accurately as well as further ahead than similar models which neglect population characteristics and suspended data. This work presents a compelling concept for longer-range fault prognosis utilising available information more fully and accurately.