Sparse modeling using orthogonal forward regression with PRESS statistic and regularization

The paper introduces an efficient construction algorithm for obtaining sparse linear-in-the-weights regression models based on an approach of directly optimizing model generalization capability. This is achieved by utilizing the delete-1 cross validation concept and the associated leave-one-out test error also known as the predicted residual sums of squares (PRESS) statistic, without resorting to any other validation data set for model evaluation in the model construction process. Computational efficiency is ensured using an orthogonal forward regression, but the algorithm incrementally minimizes the PRESS statistic instead of the usual sum of the squared training errors. A local regularization method can naturally be incorporated into the model selection procedure to further enforce model sparsity. The proposed algorithm is fully automatic, and the user is not required to specify any criterion to terminate the model construction procedure. Comparisons with some of the existing state-of-art modeling methods are given, and several examples are included to demonstrate the ability of the proposed algorithm to effectively construct sparse models that generalize well.

A new frequency-uniform coercive boundary integral equation for acoustic scattering

A new boundary integral operator is introduced for the solution of the soundsoft acoustic scattering problem, i.e., for the exterior problem for the Helmholtz equation with Dirichlet boundary conditions. We prove that this integral operator is coercive in L2(Γ) (where Γ is the surface of the scatterer) for all Lipschitz star-shaped domains. Moreover, the coercivity is uniform in the wavenumber k = ω/c, where ω is the frequency and c is the speed of sound. The new boundary integral operator, which we call the “star-combined” potential operator, is a slight modification of the standard combined potential operator, and is shown to be as easy to implement as the standard one. Additionally, to the authors' knowledge, it is the only second-kind integral operator for which convergence of the Galerkin method in L2(Γ) is proved without smoothness assumptions on Γ except that it is Lipschitz. The coercivity of the star-combined operator implies frequency-explicit error bounds for the Galerkin method for any approximation space. In particular, these error estimates apply to several hybrid asymptoticnumerical methods developed recently that provide robust approximations in the high-frequency case. The proof of coercivity of the star-combined operator critically relies on an identity first introduced by Morawetz and Ludwig in 1968, supplemented further by more recent harmonic analysis techniques for Lipschitz domains.

Lot sizing and furnace scheduling in small foundries

A lot sizing and scheduling problem prevalent in small market-driven foundries is studied. There are two related decision levels: (I the furnace scheduling of metal alloy production, and (2) moulding machine planning which specifies the type and size of production lots. A mixed integer programming (MIP) formulation of the problem is proposed, but is impractical to solve in reasonable computing time for non-small instances. As a result, a faster relax-and-fix (RF) approach is developed that can also be used on a rolling horizon basis where only immediate-term schedules are implemented. As well as a MIP method to solve the basic RF approach, three variants of a local search method are also developed and tested using instances based on the literature. Finally, foundry-based tests with a real-order book resulted in a very substantial reduction of delivery delays and finished inventory, better use of capacity, and much faster schedule definition compared to the foundry`s own practice. (c) 2006 Elsevier Ltd. All rights reserved.

Hyperfine interaction studies with (181)Ta and (111)Cd probes in the compound Ti(2)Ag

By using the time-differential perturbed angular correlation technique, the electric field gradients (EFG) at (181)Hf/(181)Ta and (111)In/(111)Cd probe sites in the MoSi(2)-type compound Ti(2)Ag have been measured as a function of temperature in the range from 24 to 1073 K. Ab initio EFG calculations have been performed within the framework of density functional theory using the full-potential augmented plane wave + local orbitals method as implemented in the WIEN2k package. These calculations allowed assignments of the probe lattice sites. For Ta, a single well-defined EFG with very weak temperature dependence was established and attributed to the [4(e)4mm] Ti site. For (111)Cd probes, two of the three measured EFGs are well defined and correlated with substitutional lattice sites, i.e. both the [4(e)4mm] Ti site and the [2(a)4/mmm] Ag site.

On estimation and influence diagnostics for the Grubbs` model under heavy-tailed distributions

The Grubbs` measurement model is frequently used to compare several measuring devices. It is common to assume that the random terms have a normal distribution. However, such assumption makes the inference vulnerable to outlying observations, whereas scale mixtures of normal distributions have been an interesting alternative to produce robust estimates, keeping the elegancy and simplicity of the maximum likelihood theory. The aim of this paper is to develop an EM-type algorithm for the parameter estimation, and to use the local influence method to assess the robustness aspects of these parameter estimates under some usual perturbation schemes, In order to identify outliers and to criticize the model building we use the local influence procedure in a Study to compare the precision of several thermocouples. (C) 2008 Elsevier B.V. All rights reserved.

Análise de onda completa de superfície seletiva em freqüência do tipo anteparo duplo

This work presents a theoretical analysis and numerical and experimental results of the scattering characteristics of frequency selective surfaces, using elements of type patch perfectly conductor. The structures are composed of two frequency selective surfaces on isotropic dielectric substrates cascaded, separated by a layer of air. The analysis is performed using the method of equivalent transmission line in combination with the Galerkin method, to determine the transmission and reflection characteristics of the structures analyzed. Specifically, the analysis uses the impedance method, which models the structure by an equivalent circuit, and applies the theory of transmission lines to determine the dyadic Green's function for the cascade structure. This function relates the incident field and surface current densities. These fields are determined algebraically by means of potential incidents and the imposition of the continuity of the fields in the dielectric interfaces. The Galerkin method is applied to the numerical determination of the unknown weight coefficients and hence the unknown densities of surface currents, which are expanded in terms of known basis functions multiplied by these weight coefficients. From the determination of these functions, it becomes possible to obtain numerical scattered fields at the top and bottom of the structures and characteristics of transmission and reflection of these structures. At work, we present numerical and experimental results for the characteristics of transmission and reflection. Comparisons were made with other results presented in literature, and it was observed a good agreement in the cases presented suggestions continuity of the work are presented

Análise espectral de reflectarrays com substratos de duas camadas dielétricas anisotrópicas uniaxiais

Recently, an amazing development has been observed in telecommunication systems. Two good examples of this development are observed in mobile communication and aerospace systems. This impressive development is related to the increasing need for receiving and transmitting communication signals. Particularly, this development has required the study of new antennas and filters. This work presents a fullwave analysis of reflectarrays. The considered structures are composed by arrays of rectangular conducting patches printed on multilayer dieletric substrates, that are mounted on a ground plane. The analysis is developed in the spectral domain, using an equivalent transmission line method in combination with Galerkin method. Results for the reflection coefficient of these structures are presented and compared to those available in the literature. A good agreement was observed. Particularly, the developed analysis uses the transmission lines theory in combination with the incident potentials and the field continuity equations, at the structures interfaces, for obtaining the scattered field components expressions as function of the patch surface currents and of the incident field. Galerkin method is used to determine the unknown coefficients in the boundary value problem. Curves for the reflection coefficient of several reflectarray geometries are presented as function of frequency and of the structural parameters

Análise de antenas de microfita com patches circulares sobre substratos anisotrópicos usando o método dos potenciais de Hertz

This work consists on the theoretical and numerical analysis of some properties of circular microstrip patch antennas on isotropic and uniaxial anisotropic substrates. For this purpose, a full wave analysis is performed, using Hertz Vector Potentials method in the Hankel Transform domain. In the numerical analysis, the moment method is also used in order to determine some characteristics of the antenna, such as: resonant frequency and radiation pattern. The definition of Hertz potentials in the Hankel domain is used in association with Maxwell´s equations and the boundary conditions of the structures to obtain the Green´s functions, relating the components of the current density on the patch and the tangential electric field components. Then, the Galerkin method is used to generate a matrix equation whose nontrivial solution is the complex resonant frequency of the structure. In the analysis, a microstrip antenna with only one isotropic dielectric layer is initially considered. For this structure, the effect of using superconductor patches is also analyzed. An analysis of a circular microstrip antenna on an uniaxial anisotropic dielectric layer is performed, using the Hertz vector potentials oriented along the optical axis of the material, that is perpendicular to the microstrip ground plane. Afterwards, the circular microstrip antenna using two uniaxial anisotropic dielectric layers is investigated, considering the particular case in which the inferior layer is filled by air. In this study, numerical results for resonant frequency and radiation pattern for circular microstrip antennas on isotropic and uniaxial anisotropic substrates are presented and compared with measured and calculated results found in the literature

Análise espectral de reflectarrays com substrato de duas camadas dielétricas anisotrópicas uniaxiais

Recently, an amazing development has been observed in telecommunication systems. Two good examples of this development are observed in mobile communication and aerospace systems. This impressive development is related to the increasing need for receiving and transmitting communication signals. Particularly, this development has required the study of new antennas and filters. This work presents a fullwave analysis of reflectarrays. The considered structures are composed by arrays of rectangular conducting patches printed on multilayer dieletric substrates, that are mounted on a ground plane. The analysis is developed in the spectral domain, using an equivalent transmission line method in combination with Galerkin method. Results for the reflection coefficient of these structures are presented and compared to those available in the literature. A good agreement was observed. Particularly, the developed analysis uses the transmission lines theory in combination with the incident potentials and the field continuity equations, at the structures interfaces, for obtaining the scattered field components expressions as function of the patch surface currents and of the incident field. Galerkin method is used to determine the unknown coefficients in the boundary value problem. Curves for the reflection coefficient of several reflectarray geometries are presented as function of frequency and of the structural parameters

The application of interpolating MLS approximations to the analysis of MHD flows

The element-free Galerkin method (EFGM) is a very attractive technique for solutions of partial differential equations, since it makes use of nodal point configurations which do not require a mesh. Therefore, it differs from FEM-like approaches by avoiding the need of meshing, a very demanding task for complicated geometry problems. However, the imposition of boundary conditions is not straightforward, since the EFGM is based on moving-least-squares (MLS) approximations which are not necessarily interpolants. This feature requires, for instance, the introduction of modified functionals with additional unknown parameters such as Lagrange multipliers, a serious drawback which leads to poor conditionings of the matrix equations. In this paper, an interpolatory formulation for MLS approximants is presented: it allows the direct introduction of boundary conditions, reducing the processing time and improving the condition numbers. The formulation is applied to the study of two-dimensional magnetohydrodynamic flow problems, and the computed results confirm the accuracy and correctness of the proposed formulation. (C) 2002 Elsevier B.V. All rights reserved.

Lot sizing and furnace scheduling in small foundries

A lot sizing and scheduling problem prevalent in small market-driven foundries is studied. There are two related decision levels: (1) the furnace scheduling of metal alloy production, and (2) moulding machine planning which specifies the type and size of production lots. A mixed integer programming (MIP) formulation of the problem is proposed, but is impractical to solve in reasonable computing time for non-small instances. As a result, a faster relax-and-fix (RF) approach is developed that can also be used on a rolling horizon basis where only immediate-term schedules are implemented. As well as a MIP method to solve the basic RF approach, three variants of a local search method are also developed and tested using instances based on the literature. Finally, foundry-based tests with a real-order book resulted in a very substantial reduction of delivery delays and finished inventory, better use of capacity, and much faster schedule definition compared to the foundry's own practice. © 2006 Elsevier Ltd. All rights reserved.

Forced oscillations with continuum models of atomic force microscopy

The dynamics of the AFM-atomic force microscope follows a model based in a Timoshenko cantilever beam with a tip attached at the free end and acting with the surface of a sample. General boundary conditions arise when the tip is either in contact or non-contact with the surface. The governing equations are given in matrix conservative form subject to localized loads. The eigenanalysis is done with a fundamental matrix response of a damped second-order matrix differential equation. Forced responses are found by using a Galerkin approximation of the matrix impulse response. Simulations results with harmonic and pulse forcing show the filtering character and the effects of the tip-sample interaction at the end of the beam. © 2012 American Institute of Physics.

Método numérico-analítico generalizado para estimação do campo eletromagnético de linhas de transmissão de energia elétrica utilizando a teoria dos elementos finitos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

NDVI and vegetation phenology dynamics under the influence of sunshine duration on the Tibetan plateau

Climate change has significantly influenced vegetation dynamics on the Tibetan Plateau (TP). Past research mainly focused on vegetation responses to temperature variation and water stress, but the influence of sunshine duration on NDVI and vegetation phenology on the TP is not well understood. In this study, NDVI time series from 1982-2008 were used to retrieve spatiotemporal vegetation dynamics on the TP. Empirical orthogonal function (EOF) analysis was conducted to understand the spatiotemporal variations of NDVI. The Start of Season (SOS) was estimated from NDVI time series with a local threshold method. The first EOF, accounting for 35.1% of NDVI variations on the TP, indicates that NDVI variations are larger in areas with shorter sunshine duration. The needle-leaved forest and shrub in the southeastern TP are more sensitive to sunshine duration anomalies (p < 0.01) than broad-leaved forest, steppe, and meadow due to spatial and altitudinal distribution of sunshine duration and vegetation types. The decrease in sunshine duration for the growing season on the TP has resulted in a decreased NDVI trend in some areas of southeastern TP (p ranging from 0.32-0.05 with threshold ranging from 0.05 to 0.25) in spite of the overall NDVI increase. SOS dynamics in most parts of the TP were mainly related to temperature variability, with precipitation and sunshine duration playing a role in a few regions. This study enhances our understanding of vegetation responses to climatic change on the TP.

Solução das equações de Navier-Stokes e da equação de transporte por um método de elementos finitos estabilizado pela técnica de separação baseada na característica

Pós-graduação em Engenharia Mecânica - FEIS