Optimal data fitting on lie groups: A coset approach

This work considers the problem of fitting data on a Lie group by a coset of a compact subgroup. This problem can be seen as an extension of the problem of fitting affine subspaces in n to data which can be solved using principal component analysis. We show how the fitting problem can be reduced for biinvariant distances to a generalized mean calculation on an homogeneous space. For biinvariant Riemannian distances we provide an algorithm based on the Karcher mean gradient algorithm. We illustrate our approach by some examples on SO(n). © 2010 Springer -Verlag Berlin Heidelberg.

A Finite Element Numerical Algorithm for Modelling and Data Fitting in Complex Systems

Numerical modelling methodologies are important by their application to engineering and scientific problems, because there are processes where analytical mathematical expressions cannot be obtained to model them. When the only available information is a set of experimental values for the variables that determine the state of the system, the modelling problem is equivalent to determining the hyper-surface that best fits the data. This paper presents a methodology based on the Galerkin formulation of the finite elements method to obtain representations of relationships that are defined a priori, between a set of variables: y = z(x1, x2,...., xd). These representations are generated from the values of the variables in the experimental data. The approximation, piecewise, is an element of a Sobolev space and has derivatives defined in a general sense into this space. The using of this approach results in the need of inverting a linear system with a structure that allows a fast solver algorithm. The algorithm can be used in a variety of fields, being a multidisciplinary tool. The validity of the methodology is studied considering two real applications: a problem in hydrodynamics and a problem of engineering related to fluids, heat and transport in an energy generation plant. Also a test of the predictive capacity of the methodology is performed using a cross-validation method.

An Accelerated Chow and Liu Algorithm: Fitting Tree Distributions to High Dimensional Sparse Data

Chow and Liu introduced an algorithm for fitting a multivariate distribution with a tree (i.e. a density model that assumes that there are only pairwise dependencies between variables) and that the graph of these dependencies is a spanning tree. The original algorithm is quadratic in the dimesion of the domain, and linear in the number of data points that define the target distribution $P$. This paper shows that for sparse, discrete data, fitting a tree distribution can be done in time and memory that is jointly subquadratic in the number of variables and the size of the data set. The new algorithm, called the acCL algorithm, takes advantage of the sparsity of the data to accelerate the computation of pairwise marginals and the sorting of the resulting mutual informations, achieving speed ups of up to 2-3 orders of magnitude in the experiments.

Penalized Least Squares Fitting

* Supported by the Army Research Office under grant DAAD-19-02-10059.

Understanding mild acid pretreatment of sugarcane bagasse through particle scale modeling

Sugarcane bagasse is an abundant and sustainable resource, generated as a by-product of sugarcane milling. The cellulosic material within bagasse can be broken down into glucose molecules and fermented to produce ethanol, making it a promising feedstock for biofuel production. Mild acid pretreatment hydrolyses the hemicellulosic component of biomass, thus allowing enzymes greater access to the cellulosic substrate during saccharification. A particle-scale mathematical model describing the mild acid pretreatment of sugarcane bagasse has been developed, using a volume averaged framework. Discrete population-balance equations are used to characterise the polymer degradation kinetics, and diffusive effects account for mass transport within the cell wall of the bagasse. As the fibrous material hydrolyses over time, variations in the porosity of the cell wall and the downstream effects on the reaction kinetics are accounted for using conservation of volume arguments. Non-dimensionalization of the model equations reduces the number of parameters in the system to a set of four dimensionless ratios that compare the timescales of different reaction and diffusion events. Theoretical yield curves are compared to macroscopic experimental observations from the literature and inferences are made as to constraints on these “unknown” parameters. These results enable connections to be made between experimental data and the underlying thermodynamics of acid pretreatment. Consequently, the results suggest that data-fitting techniques used to obtain kinetic parameters should be carefully applied, with prudent consideration given to the chemical and physiological processes being modeled.

Control dynamics of severe acute respiratory syndrome transmission

Severe acute respiratory syndrome (SARS) is a serious disease with many puzzling features. We present a simple, dynamic model to assess the epidemic potential of SARS and the effectiveness of control measures. With this model, we analysed the SARS epidemic data in Beijing. The data fitting gives the basic case reproduction number of 2.16 leading to the outbreak, and the variation of the effective reproduction number reflecting the control effect. Noticeably, our study shows that the response time and the strength of control measures have significant effects on the scale of the outbreak and the lasting time of the epidemic.

One-dimensional models for organic magnetic materials

In the preparation of small organic paramagnets, these structures may conceptually be divided into spin-containing units (SCs) and ferromagnetic coupling units (FCs). The synthesis and direct observation of a series of hydrocarbon tetraradicals designed to test the ferromagnetic coupling ability of m-phenylene, 1,3-cyclobutane, 1,3- cyclopentane, and 2,4-adamantane (a chair 1,3-cyclohexane) using Berson TMMs and cyclobutanediyls as SCs are described. While 1,3-cyclobutane and m-phenylene are good ferromagnetic coupling units under these conditions, the ferromagnetic coupling ability of 1,3-cyclopentane is poor, and 1,3-cyclohexane is apparently an antiferromagnetic coupling unit. In addition, this is the first report of ferromagnetic coupling between the spins of localized biradical SCs.

The poor coupling of 1,3-cyclopentane has enabled a study of the variable temperature behavior of a 1,3-cyclopentane FC-based tetraradical in its triplet state. Through fitting the observed data to the usual Boltzman statistics, we have been able to determine the separation of the ground quintet and excited triplet states. From this data, we have inferred the singlet-triplet gap in 1,3-cyclopentanediyl to be 900 cal/mol, in remarkable agreement with theoretical predictions of this number.

The ability to simulate EPR spectra has been crucial to the assignments made here. A powder EPR simulation package is described that uses the Zeeman and dipolar terms to calculate powder EPR spectra for triplet and quintet states.

Methods for characterizing paramagnetic samples by SQUID magnetometry have been developed, including robust routines for data fitting and analysis. A precursor to a potentially magnetic polymer was prepared by ring-opening metathesis polymerization (ROMP), and doped samples of this polymer were studied by magnetometry. While the present results are not positive, calculations have suggested modifications in this structure which should lead to the desired behavior.

Source listings for all computer programs are given in the appendix.

溶胶-凝胶条波导Sagnac环的乙醇蒸气传感特性

采用有机/无机混合溶胶-凝胶法制作条形光波导,并将条波导接入光纤Sagnac 环中,测量了输出光功率随环境气氛中乙醇蒸气体积分数变化的特性,表明在实验研究的范围内,输出信号与乙醇蒸气体积分数呈正弦变化。根据Sagnac环结构输出特性的基本关系,反映了溶胶-凝胶条波导在乙醇蒸气气氛下产生了双折射效应。观察到双折射相移与乙醇体积分数的亚线性关系。对实验数据拟合,计算了偏振相移的线性项和二次项系数,得到所制备的条波导的双折射对乙醇体积分数的响应为Δn≈4.4×10-2。测量了信号变化的时间演变特性,典型的上升和

Mudanças de Estado e Multiplicadores Fiscais no Brasil entre 1999-2012

Este trabalho avalia o comportamento dos multiplicadores fiscais no Brasil entre 1999-2012. Para tanto, utiliza a metodologia desenvolvida por Sims, Waggoner e Zha (2008), que é um procedimento Bayesiano de estimação no qual os parâmetros do modelo mudam com alterações no estado da economia e os estados (regimes) seguem um processo de mudança de regime markoviano. Ou seja, foi estimado um modelo VAR Estrutural Bayesiano com mudança de regimes Markoviana (Markov Switching Structural Bayesian Vector Autoregression - MS-SBVAR). A base de dados é composta pelo consumo da administração pública, pela formação bruta de capital fixo da administração pública, pela carga tributária líquida e pelo Produto Interno Bruto (PIB), das três esferas do governo (federal, estadual, incluindo o Distrito Federal, e municipal). O software MATLAB/Dynare foi utilizado na estimação dos modelos e os resultados sugerem a ocorrência de 2 ou 3 regimes nos dois modelos que melhor se ajustaram aos dados. Os multiplicadores estimados apresentaram os sinais esperados e os diferentes tipos de multiplicadores fiscais calculados apresentaram valores maiores para a resposta do PIB a choques na formação bruta de capital fixo da administração pública que são eficazes, uma vez que possuem valores maiores do que um e impacto de longo prazo no PIB - quando comparado aos choques no consumo da administração pública, que possuem pouca persistência e são ineficazes (menores do que um), além de uma resposta negativa e persistente do PIB a choques na carga tributária líquida. Os resultados obtidos não indicam, ainda, multiplicadores fiscais maiores em regimes com maior variância nos resíduos do modelo.

Kinetic spectra of light-adaptation dark-adaptation and M-intermediate of BR-D96N

BR-D96N is a kind of genetically site-specific mutants of bacteriorhodopsin (BR) with obvious photochromic effect. Compared to the wild type BR, the lifetime of M state of BR-D96N is prolonged to several minutes so that the photochromic kinetics and the intermediates formation can be studied by the conventional spectra analysis. In the experiment, the absorption spectra of the sample at different time after light illumination are measured with spectrophotometer. By fitting and analyzing the variation of the spectra, we suppose that there are three main states in the, photochromic process, i.e., B state (light-adapted state), M state and D state (dark-adapted state). The absorption spectra of the B state, M state and D state are extracted from the experimental data based on this three-state model and the spectra at various time are fitted with the least-square method. So, the variations of population percentages of the M state, B state and D state are obtained and the M state and B state lifetimes are estimated. In another way, from the measurement of the absorption dynamics at 407 and 568 nm, the M state and B state lifetimes are also obtained by two exponential data fitting, which give coincident results with those of the spectra analysis. (C) 2003 Elsevier Science B.V. All rights reserved.

Multi-exponentially photoelectric response of bacteriorhodopsin

A thin oriented bacteriorhodopsin (bR) him is deposited on a stainless steel slide by use of the electrophoretic sedimentation method. A junction is made with electrolyte gels having a counterelectrode to construct a bR-based photoelectric detector;. The photoelectric response signal to a 10 ns laser pulse is measured. A theory on the photoelectric kinetics of bR is developed based on the concept of the charge displacement current and the bR photocycle rate equations. Comparison between the theoretical and experimental results proves that the bR photoelectric response to a short laser pulse is a multi-exponential process. The decay time constants and amplitudes of each, exponential component are obtained by data fitting.

Kinetics of picosecond laser pulse induced charge separation and proton transfer in bacteriorhodopsin

Bacteriorhodopsin (BR) films oriented by an electrophoretic method are deposited on a transparent conductive ITO glass. A counterelectrode of copper and gelose gel is used to compose a sandwich-type photodetector with the structure of ITO/BR film/gelose gel/Cu. A single 30-ps laser pulse and a mode-locked pulse train are respectively used to excite the BR photodetector. The ultrafast failing edge and the bipolar response signal are measured by the digital oscilloscope under seven different time ranges. Marquardt nonlinear least squares fitting is used to fit all the experimental data and a good fitting equation is found to describe the kinetic process of the photoelectric signal. Data fitting resolves six exponential components that can be assigned to a seven-step BR photocycle model: BR-->K-->KL-->L-->M-->N-->O-->BR. Comparing tests of the BR photodetector with a 100-ps Si PIN photodiode demonstrates that this type of BIR photocletector has at least 100-ps response time and can also serve as a fast photoelectric switch. (C) 2003 Society of Photo-Optical Instrumentation Engineers.

Photochromic kinetic spectra and intermediates of BR-D96N

BR-D96N is a kind of genetically site-specific mutant of bacteriorhodopsin (BR) with obvious photochromic effect. Compared to the wild type BR, the lifetime of M state of BR-D96N is prolonged to several minutes so that the photochromic kinetics and the intermediates formation can be studied by the conventional spectra analysis. In the experiment, the absorption spectra of the sample at different time after light illumination are measured with spectrophotometer. By fitting and analyzing the variation of the spectra, we suppose there-are three main states in the photochromic process, i.e. B state (light-adapted state), M state and D state (dark-adapted state). The absorption spectra of the B state, M state and D state are extracted from the experimental data based on this three-state model and the spectra at various time are fitted With the least square method. So, the variations of population percentages of the M state, B state and D state are obtained and the M state and B state lifetimes are estimated. In another,way, from measuring the absorption dynamics at 407 nm and 568 nm, the M state and B state lifetimes are also obtained by two exponential data fitting, which give coincident results with those of the spectra analysis.

Phase agile reflectarray cells based on liquid crystals

Experimental results at X-band are used to compare the electromagnetic scattering from a printed reflectarray cell, which is constructed on 500 mu m thick layers of three different nematic state liquid crystals. It is shown that a small voltage can be used to vary the permittivity of the tunable substrate and thereby control the phase of the reflected signals. Numerical results using Ansoft HFSS are compared with the measured phase, resonant frequencies and signal attenuation for two orientations of the liquid crystal molecules. Data fitting is employed to quantify the loss tangent and the permittivity values of the three anisotropic specimens. The performance trade-offs that are imposed by the use of commercially available materials are discussed, and the computer model is used to specify the electrical properties of a liquid crystal mixture, which can provide a signal loss of <1 dB and a dynamic phase range of 300 degrees from the patch elements at 10 GHz.