ℓasso MPC: Smart regulation of over-actuated systems

In this paper, a novel MPC strategy is proposed, and referred to as asso MPC. The new paradigm features an 1-regularised least squares loss function, in which the control error variance competes with the sum of input channels magnitude (or slew rate) over the whole horizon length. This cost choice is motivated by the successful development of LASSO theory in signal processing and machine learning. In the latter fields, sum-of-norms regularisation have shown a strong capability to provide robust and sparse solutions for system identification and feature selection. In this paper, a discrete-time dual-mode asso MPC is formulated, and its stability is proven by application of standard MPC arguments. The controller is then tested for the problem of ship course keeping and roll reduction with rudder and fins, in a directional stochastic sea. Simulations show the asso MPC to inherit positive features from its corresponding regressor: extreme reduction of decision variables' magnitude, namely, actuators' magnitude (or variations), with a finite energy error, being particularly promising for over-actuated systems. © 2012 AACC American Automatic Control Council).

Improved bounds for subband-adaptive iterative shrinkage/thresholding algorithms.

This paper presents new methods for computing the step sizes of the subband-adaptive iterative shrinkage-thresholding algorithms proposed by Bayram & Selesnick and Vonesch & Unser. The method yields tighter wavelet-domain bounds of the system matrix, thus leading to improved convergence speeds. It is directly applicable to non-redundant wavelet bases, and we also adapt it for cases of redundant frames. It turns out that the simplest and most intuitive setting for the step sizes that ignores subband aliasing is often satisfactory in practice. We show that our methods can be used to advantage with reweighted least squares penalty functions as well as L1 penalties. We emphasize that the algorithms presented here are suitable for performing inverse filtering on very large datasets, including 3D data, since inversions are applied only to diagonal matrices and fast transforms are used to achieve all matrix-vector products.

Free energy surface reconstruction from umbrella samples using Gaussian process regression II: Multiple collective variables

We demonstrate how a prior assumption of smoothness can be used to enhance the reconstruction of free energy profiles from multiple umbrella sampling simulations using the Bayesian Gaussian process regression approach. The method we derive allows the concurrent use of histograms and free energy gradients and can easily be extended to include further data. In Part I we review the necessary theory and test the method for one collective variable. We demonstrate improved performance with respect to the weighted histogram analysis method and obtain meaningful error bars without any significant additional computation. In Part II we consider the case of multiple collective variables and compare to a reconstruction using least squares fitting of radial basis functions. We find substantial improvements in the regimes of spatially sparse data or short sampling trajectories. A software implementation is made available on www.libatoms.org.

Quantitative structure-property relationships for octanol-air partition coefficients of polychlorinated naphthalenes, chlorobenzenes and p,p '-DDT

The octanol-air partition coefficient (K-OA) is a key descriptor of chemicals partitioning between the atmosphere and environmental organic phases. Quantitative structure-property relationships (QSPR) are necessary to model and predict KOA from molecular structures. Based on 12 quantum chemical descriptors computed by the PM3 Hamiltonian, using partial least squares (PLS) analysis, a QSPR model for logarithms of K-OA to base 10 (log K-OA) for polychlorinated naphthalenes (PCNs), chlorobenzenes and p,p'-DDT was obtained. The cross-validated Q(cum)(2) value of the model is 0.973, indicating a good predictive ability of the model. The main factors governing log K-OA of the PCNs, chlorobenzenes, and p,p'-DDT are, in order of decreasing importance, molecular size and molecular ability of donating/accepting electrons to participate in intermolecular interactions. The intermolecular dispersive interactions play a leading role in governing log K-OA. The more chlorines in PCN and chlorobenzene molecules, the greater the log K-OA values. Increasing E-LUMO (the energy of the lowest unoccupied molecular orbital) of the molecules leads to decreasing log K-OA values, implying possible intermolecular interactions between the molecules under study and octanol molecules. (C) 2002 Elsevier Science Ltd. All rights reserved.

Quantitative relationships between molecular structures, environmental temperatures and octanol-air partition coefficients of PCDD/Fs

A concise quantitative model that incorporates information on both environmental temperature M and molecular structures, for logarithm of octanol-air partition coefficient (K-OA) to base 10 (logK(OA)) of PCDDs, was developed. Partial least squares (PLS) analysis together with 14 quantum chemical descriptors were used to develop the quantitative relationships between structures, environmental temperatures and properties (QRSETP) model. It has been validated that the obtained QRSETP model can be used to predict logK(OA) of other PCDDs. Molecular size, environmental temperature (T), q(+) (the most positive net atomic charge on hydrogen or chlorine atoms in PCDD molecules) and E-LUMO (the energy of the lowest unoccupied molecular orbital) are main factors governing logK(OA) of PCDD/Fs under study. The intermolecular dispersive interactions and thus the size of the molecules play a leading role in governing logK(OA). The more chlorines in PCDD molecules, the greater the logK(OA) values. Increasing E-LUMO values of the molecules leads to decreasing logK(OA) values, implying possible intermolecular interactions between the molecules under study and octanol molecules. Greater q(+) values results in greater intermolecular electrostatic repulsive interactions between PCDD and octanol molecules and smaller logK(OA) values. (C) 2002 Elsevier Science B.V. All rights reserved.

Quantitative structure-property relationships for octanol-air partition coefficients of polychlorinated biphenyls

Based on nine quantum chemical descriptors computed by PM3 Hamiltonian, using partial least squares analysis, a significant quantitative structure-property relationship for the logarithm of octanol-air partition coefficients (log K-OA) of polychlorinated biphenyls (PCBs) was obtained. The cross-validated Q(cum)(2) value of the model is 0.962, indicating a good predictive ability. The intermolecular dispersive interactions and thus the size of the PCB molecules play a key role in governing log K-OA. The greater the size of PCB molecules, the greater the log K-OA values. Increasing E-LUMO (the energy of the lowest unoccupied molecular orbital) values of the PCBs leads to decreasing log K-OA values, indicating possible interactions between PCB and octanol molecules. Increasing Q(Cl)(+) (the most positive net atomic charges on a chlorine atom) and Q(C)(-) (the largest negative net atomic charge on a carbon atom) values of PCBs results in decreasing log K-OA values, implying possible intermolecular electrostatic interactions between octanol and PCB molecules. (C) 2002 Elsevier Science Ltd. All rights reserved.

Quantitative structure-property relationships (QSPRs) on direct photolysis of PCDDs

By the use of partial least squares (PLS) method and 27 quantum chemical descriptors computed by PM3 Hamiltonian, a statistically significant QSPR were obtained for direct photolysis quantum yields (Y) of selected Polychlorinated dibenzo-p-dioxins (PCDDs). The QSPR can be used for prediction. The direct photolysis quantum yields of the PCDDs are dependent on the number of chlorine atoms bonded with the parent structures, the character of the carbon-oxygen bonds, and molecular polarity. Increasing bulkness and polarity of PCDDs lead to decrease of log Y values. Increasing the frontier molecular orbital energies (E-lumo and E-homo) and heat of formation (HOF) values leads to increase of log Y values. (C) 2001 Elsevier Science Ltd. All rights reserved.

Quantitative structure-property relationship study on reductive dehalogenation of selected halogenated aliphatic hydrocarbons in sediment slurries

In this study, by the use of partial least squares (PLS) method and 26 quantum chemical descriptors computed by PM3 Hamiltonian, a quantitative structure-property relationship (QSPR) model was developed for reductive dehalogenation rate constants of 13 halogenated aliphatic compounds in sediment slurry under anaerobic conditions. The model can be used to explain the dehalogenation mechanism. Halogenated aliphatic compounds with great energy of the lowest unoccupied molecular orbital (E-lumo), total energy (TE), electronic energy (EE), the smallest bond order of the carbon-halogen bonds (BO) and the most positive net atomic charges on an atom of the molecule (q(+)) values tend to be reductively dehalogenated slow, whereas halogenated aliphatic compounds with high values of molecular weight (Mw), average molecular polarizability (a) and core-core repulsion energy (CCR) values tend to be reductively dehalogenated fastest. (C) 2001 Published by Elsevier Science Ltd.

Quantitative structure-property relationship studies on n-octanol/water partitioning coefficients of PCDD/Fs

Based on some fundamental quantum chemical descriptors computed by PM3 Hamiltonian, by the use of partial least-squares (PLS) analysis, a significant quantitative structure-property relationship (QSPR) model for logK(ow) of polychlorinated dibenzo-p-dioxins and dibenzo-p-furans (PCDD/Fs) was obtained. The QSPR can be used for prediction. The intermolecular dispersive interactions and thus the bulkness of the PCDD/Fs are the main factors affecting the logK(ow). The more chlorines in the PCDD/F molecule, the greater the logK(ow) values. (C) 2001 Elsevier Science Ltd. All rights reserved.

QtUCP-A program for determining unit-cell parameters in electron diffraction experiments using double-tilt and rotation-tilt holders

A computer program, QtUCP, has been developed based on several well-established algorithms using GCC 4.0 and Qt (R) 4.0 (Open Source Edition) under Debian GNU/Linux 4.0r0. it can determine the unit-cell parameters from an electron diffraction tilt series obtained from both double-tilt and rotation-tilt holders. In this approach, two or more primitive cells of the reciprocal lattice are determined from experimental data, in the meantime, the measurement errors of the tilt angles are checked and minimized. Subsequently, the derived primitive cells are converted into the reduced form and then transformed into the reduced direct primitive cell. Finally all the patterns are indexed and the least-squares refinement is employed to obtain the optimized results of the lattice parameters. Finally, two examples are given to show the application of the program, one is based on the experiment, the other is from the simulation. (C) 2008 Elsevier B.V. All rights reserved.

Studies of tetragonal PbTiO3 subjected to uniaxial stress along the c-axis

Tetragonal PbTiO3 under uniaxial stress along the c-axis is investigated from first-principles. The structural parameters, polarization, and squares of the lowest optical phonon frequencies for E(1TO) and A(1)(1TO) modes at Gamma show abrupt changes near a stress sigma(c) of 1.04 GPa, which is related to the dramatic change of elastic constant c(33) resulting from the uniaxial stress applied along the c-axis. We also find that the uniaxial compressive stress could enhance the piezoelectric stress coefficients, whereas the uniaxial tensile stress could enhance the piezoelectric strain coefficients. It is also found that when the magnitude of uniaxial compressive stress sigma(33) is greater than 12 GPa, PbTiO3 is transformed to the paraelectric tetragonal phase.

Analysis of mode characteristics for equilateral-polygonal resonators with a center hole

Mode characteristics for equilateral triangles, squares, and hexagonal resonators with a center hole are numerically simulated by the finite-different time domain (FDTD) technique. The center hole does not break the symmetry behavior of the original resonators and can result in modification of the mode field patterns and mode Q factors. In an equilateral triangle resonator the center hole can suppress the symmetry state of degenerate states with the merit of single mode operation. In a square resonator, the Q factor can be enhanced for some modes with a suitable size of the hole. For a hexagonal resonator with a side length of 1 mu m and a refractive index of 3.2, the mode Q factors first gradually decrease with the increase of the hole diameter for modes at a wavelength of about 1500 nm, then the modes transform to that of a microdisk with a jump of the mode wavelength as the hole diameter approaches 0.7 mu m. Finally, the mode Q factors greatly enhance as the hole diameter reaches about 1 mu m. The results indicate that the center hole can greatly modify mode characteristics, especially that of the mode Q factor. (C) 2009 Optical Society of America

Calculation of light delay for coupled microrings by FDTD technique and Pade approximation

The Pade approximation with Baker's algorithm is compared with the least-squares Prony method and the generalized pencil-of-functions (GPOF) method for calculating mode frequencies and mode Q factors for coupled optical microdisks by FDTD technique. Comparisons of intensity spectra and the corresponding mode frequencies and Q factors show that the Pade approximation can yield more stable results than the Prony and the GPOF methods, especially the intensity spectrum. The results of the Prony method and the GPOF method are greatly influenced by the selected number of resonant modes, which need to be optimized during the data processing, in addition to the length of the time response signal. Furthermore, the Pade approximation is applied to calculate light delay for embedded microring resonators from complex transmission spectra obtained by the Pade approximation from a FDTD output. The Prony and the GPOF methods cannot be applied to calculate the transmission spectra, because the transmission signal obtained by the FDTD simulation cannot be expressed as a sum of damped complex exponentials. (C) 2009 Optical Society of America

Measurement of a reciprocal four-port transmission line structure using the 16-term error model

A new method to measure reciprocal four-port structures, using a 16-term error model, is presented. The measurement is based on 5 two-port calibration standards connected to two of the ports, while the network analyzer is connected to the two remaining ports. Least-squares-fit data reduction techniques are used to lower error sensitivity. The effect of connectors is deembedded using closed-form equations. (C) 2007 Wiley Periodicals, Inc.

A density-functional study of small titanium clusters

The atomic structures and electronic properties of small Ti-N (N=2-10) clusters have been studied by using the density-functional theory with a local spin density approximation. We find that the inner-shells (3s3p) of the titanium atom plays an important role in the formation of the small clusters. We have obtained the ground state of titanium clusters, Ti-7 is found to be a magic cluster, which is in good agreement with the experimental results. Starting with Ti-8 cluster some features of the electronic structure of the titanium bulk have been developed. The ionization potentials and magnetic moments for these small titanium clusters are also presented. (C) 2000 American Institute of Physics. [S0021- 9606(00)30544-X].