Correct biological timing in Arabidopsis requires multiple light-signaling pathways.

Circadian oscillators provide rhythmic temporal cues for a range of biological processes in plants and animals, enabling anticipation of the day/night cycle and enhancing fitness-associated traits. We have used engineering models to understand the control principles of a plant's response to seasonal variation. We show that the seasonal changes in the timing of circadian outputs require light regulation via feed-forward loops, combining rapid light-signaling pathways with entrained circadian oscillators. Linear time-invariant models of circadian rhythms were computed for 3,503 circadian-regulated genes and for the concentration of cytosolic-free calcium to quantify the magnitude and timing of regulation by circadian oscillators and light-signaling pathways. Bioinformatic and experimental analysis show that rapid light-induced regulation of circadian outputs is associated with seasonal rephasing of the output rhythm. We identify that external coincidence is required for rephasing of multiple output rhythms, and is therefore important in general phase control in addition to specific photoperiod-dependent processes such as flowering and hypocotyl elongation. Our findings uncover a fundamental design principle of circadian regulation, and identify the importance of rapid light-signaling pathways in temporal control.

Image deconvolution using a Gaussian scale mixtures model to approximate the wavelet sparseness constraint

This paper proposes to use an extended Gaussian Scale Mixtures (GSM) model instead of the conventional ℓ1 norm to approximate the sparseness constraint in the wavelet domain. We combine this new constraint with subband-dependent minimization to formulate an iterative algorithm on two shift-invariant wavelet transforms, the Shannon wavelet transform and dual-tree complex wavelet transform (DTCWT). This extented GSM model introduces spatially varying information into the deconvolution process and thus enables the algorithm to achieve better results with fewer iterations in our experiments. ©2009 IEEE.

Multi-channel signal separation

The separation of independent sources from mixed observed data is a fundamental and challenging problem. In many practical situations, observations may be modelled as linear mixtures of a number of source signals, i.e. a linear multi-input multi-output system. A typical example is speech recordings made in an acoustic environment in the presence of background noise and/or competing speakers. Other examples include EEG signals, passive sonar applications and cross-talk in data communications. In this paper, we propose iterative algorithms to solve the n × n linear time invariant system under two different constraints. Some existing solutions for 2 × 2 systems are reviewed and compared.

Vision-based global localization using a visual vocabular

This paper presents a novel coarse-to-fine global localization approach that is inspired by object recognition and text retrieval techniques. Harris-Laplace interest points characterized by SIFT descriptors are used as natural land-marks. These descriptors are indexed into two databases: an inverted index and a location database. The inverted index is built based on a visual vocabulary learned from the feature descriptors. In the location database, each location is directly represented by a set of scale invariant descriptors. The localization process consists of two stages: coarse localization and fine localization. Coarse localization from the inverted index is fast but not accurate enough; whereas localization from the location database using voting algorithm is relatively slow but more accurate. The combination of coarse and fine stages makes fast and reliable localization possible. In addition, if necessary, the localization result can be verified by epipolar geometry between the representative view in database and the view to be localized. Experimental results show that our approach is efficient and reliable. ©2005 IEEE.

Convergence of the auxiliary particle implementation of the PHD filter

Optimal Bayesian multi-target filtering is in general computationally impractical owing to the high dimensionality of the multi-target state. The Probability Hypothesis Density (PHD) filter propagates the first moment of the multi-target posterior distribution. While this reduces the dimensionality of the problem, the PHD filter still involves intractable integrals in many cases of interest. Several authors have proposed Sequential Monte Carlo (SMC) implementations of the PHD filter. However, these implementations are the equivalent of the Bootstrap Particle Filter, and the latter is well known to be inefficient. Drawing on ideas from the Auxiliary Particle Filter (APF), a SMC implementation of the PHD filter which employs auxiliary variables to enhance its efficiency was proposed by Whiteley et. al. Numerical examples were presented for two scenarios, including a challenging nonlinear observation model, to support the claim. This paper studies the theoretical properties of this auxiliary particle implementation. $\mathbb{L}_p$ error bounds are established from which almost sure convergence follows.

Auxiliary particle implementation of the probability hypothesis density filter

Optimal Bayesian multi-target filtering is, in general, computationally impractical owing to the high dimensionality of the multi-target state. The Probability Hypothesis Density (PHD) filter propagates the first moment of the multi-target posterior distribution. While this reduces the dimensionality of the problem, the PHD filter still involves intractable integrals in many cases of interest. Several authors have proposed Sequential Monte Carlo (SMC) implementations of the PHD filter. However, these implementations are the equivalent of the Bootstrap Particle Filter, and the latter is well known to be inefficient. Drawing on ideas from the Auxiliary Particle Filter (APF), we present a SMC implementation of the PHD filter which employs auxiliary variables to enhance its efficiency. Numerical examples are presented for two scenarios, including a challenging nonlinear observation model.

Quantifying the limits of HANPP and carbon emissions which prolong total species well-being

Anthropogenic climate and land-use change are leading to irreversible losses of global biodiversity, upon which ecosystem functioning depends. Since total species' well-being depends on ecosystem goods and services, man must determine how much net primary productivity (NPP) may be appropriated and carbon emitted so as to not adversely impact this and future generations. In 2005, man ought to have only appropriated 9.72 Pg C of NPP, representing a factor 2.50, or 59.93%, reduction in human-appropriated NPP in that year. Concurrently, the carbon cycle would have been balanced with a factor 1.26, or 20.84%, reduction from 7.60 Gt C/year to 5.70 Gt C/year, representing a return to the 1986 levels. This limit is in keeping with the category III stabilization scenario of the Intergovernmental Panel for Climate Change. Projecting population growth to 2030 and its associated basic food requirements, the maximum HANPP remains at 9.74 ± 0.02 Pg C/year. This time-invariant HANPP may only provide for the current global population of 6.51 billion equitably at the current average consumption of 1.49 t C per capita, calling into question the sustainability of developing countries striving for high-consuming country levels of 5.85 t C per capita and its impacts on equitable resource distribution. © Springer Science+Business Media B.V. 2009.

Development and validation of a new method for measuring friction between skin and nonwoven materials.

A new method for measuring the coefficient of friction between nonwoven materials and the curved surface of the volar forearm has been developed and validated. The method was used to measure the coefficient of static friction for three different nonwoven materials on the normal (dry) and over-hydrated volar forearms of five female volunteers (ages 18-44). The method proved simple to run and had good repeatability: the coefficient of variation (standard deviation expressed as a percentage of the mean) for triplets of repeat measurements was usually (80 per cent of the time) less than 10 per cent. Measurements involving the geometrically simpler configuration of pulling a weighted fabric sample horizontally across a quasi-planar area of volar forearm skin proved experimentally more difficult and had poorer repeatability. However, correlations between values of coefficient of static friction derived using the two methods were good (R = 0.81 for normal (dry) skin, and 0.91 for over-hydrated skin). Measurements of the coefficient of static friction for the three nonwovens for normal (dry) and for over-hydrated skin varied in the ranges of about 0.3-0.5 and 0.9-1.3, respectively. In agreement with Amontons' law, coefficients of friction were invariant with normal pressure over the entire experimental range (0.1-8.2 kPa).

Analytical approximations for the modal acoustic impedances of simply supported, rectangular plates.

Coupling of the in vacuo modes of a fluid-loaded, vibrating structure by the resulting acoustic field, while known to be negligible for sufficiently light fluids, is still only partially understood. A particularly useful structural geometry for the study of this problem is the simply supported, rectangular flat plate, since it exhibits all the relevant physical features while still admitting an analytical description of the modes. Here the influence of the fluid can be expressed in terms of a set of doubly infinite integrals over wave number: the modal acoustic impedances. Closed-form solutions for these impedances do not exist and, while their numerical evaluation is possible, it greatly increases the computational cost of solving the coupled system of modal equations. There is thus a need for accurate analytical approximations. In this work, such approximations are sought in the limit where the modal wavelength is small in comparison with the acoustic wavelength and the plate dimensions. It is shown that contour integration techniques can be used to derive analytical formulas for this regime and that these formulas agree closely with the results of numerical evaluations. Previous approximations [Davies, J. Sound Vib. 15(1), 107-126 (1971)] are assessed in the light of the new results and are shown to give a satisfactory description of real impedance components, but (in general) erroneous expressions for imaginary parts.

Optimal H2 order-one reduction by solving eigenproblems for polynomial equations

A method is given for solving an optimal H2 approximation problem for SISO linear time-invariant stable systems. The method, based on constructive algebra, guarantees that the global optimum is found; it does not involve any gradient-based search, and hence avoids the usual problems of local minima. We examine mostly the case when the model order is reduced by one, and when the original system has distinct poles. This case exhibits special structure which allows us to provide a complete solution. The problem is converted into linear algebra by exhibiting a finite-dimensional basis for a certain space, and can then be solved by eigenvalue calculations, following the methods developed by Stetter and Moeller. The use of Buchberger's algorithm is avoided by writing the first-order optimality conditions in a special form, from which a Groebner basis is immediately available. Compared with our previous work the method presented here has much smaller time and memory requirements, and can therefore be applied to systems of significantly higher McMillan degree. In addition, some hypotheses which were required in the previous work have been removed. Some examples are included.

Designing MPC controllers by reverse-engineering existing LTI controllers

This technical report presents a method for designing a constrained output-feedback model predictive controller (MPC) that behaves in the same way as an existing baseline stabilising linear time invariant output-feedback controller when constraints are inactive. The baseline controller is cast into an observer-compensator form and an inverse-optimal cost function is used as the basis of the MPC controller. The available degrees of design freedom are explored, and some guidelines provided for the selection of an appropriate observer-compensator realisation that will best allow exploitation of the constraint-handling and redundancy management capabilities of MPC. Consideration is given to output setpoint tracking, and the method is demonstrated with three different multivariable plants of varying complexity.

Asymptotic Recovery for Discrete-Time Systems

An asymptotic recovery design procedure is proposed for square, discrete-time, linear, time-invariant multivariable systems, which allows a state-feedback design to be approximately recovered by a dynamic output feedback scheme. Both the case of negligible processing time (compared to the sampling interval) and of significant processing time are discussed. In the former case, it is possible to obtain perfect. © 1985 IEEE.

H∞ sampled-data synthesis and related numerical issues

In this paper we present a new, compact derivation of state-space formulae for the so-called discretisation-based solution of the H∞ sampled-data control problem. Our approach is based on the established technique of continuous time-lifting, which is used to isometrically map the continuous-time, linear, periodically time-varying, sampled-data problem to a discretetime, linear, time-invariant problem. State-space formulae are derived for the equivalent, discrete-time problem by solving a set of two-point, boundary-value problems. The formulae accommodate a direct feed-through term from the disturbance inputs to the controlled outputs of the original plant and are simple, requiring the computation of only a single matrix exponential. It is also shown that the resultant formulae can be easily re-structured to give a numerically robust algorithm for computing the state-space matrices. © 1997 Elsevier Science Ltd. All rights reserved.