Second-order parameter estimation

This work provides a general framework for the design of second-order blind estimators without adopting anyapproximation about the observation statistics or the a prioridistribution of the parameters. The proposed solution is obtainedminimizing the estimator variance subject to some constraints onthe estimator bias. The resulting optimal estimator is found todepend on the observation fourth-order moments that can be calculatedanalytically from the known signal model. Unfortunately,in most cases, the performance of this estimator is severely limitedby the residual bias inherent to nonlinear estimation problems.To overcome this limitation, the second-order minimum varianceunbiased estimator is deduced from the general solution by assumingaccurate prior information on the vector of parameters.This small-error approximation is adopted to design iterativeestimators or trackers. It is shown that the associated varianceconstitutes the lower bound for the variance of any unbiasedestimator based on the sample covariance matrix.The paper formulation is then applied to track the angle-of-arrival(AoA) of multiple digitally-modulated sources by means ofa uniform linear array. The optimal second-order tracker is comparedwith the classical maximum likelihood (ML) blind methodsthat are shown to be quadratic in the observed data as well. Simulationshave confirmed that the discrete nature of the transmittedsymbols can be exploited to improve considerably the discriminationof near sources in medium-to-high SNR scenarios.

Trust region versus line search for computing the optical flow

We consider the numerical treatment of the optical flow problem by evaluating the performance of the trust region method versus the line search method. To the best of our knowledge, the trust region method is studied here for the first time for variational optical flow computation. Four different optical flow models are used to test the performance of the proposed algorithm combining linear and nonlinear data terms with quadratic and TV regularization. We show that trust region often performs better than line search; especially in the presence of non-linearity and non-convexity in the model.

Análisis Discriminante mediante SPSS.

El análisis discriminante es un método estadístico a través del cual se busca conocer qué variables, medidas en objetos o individuos, explican mejor la atribución de la diferencia de los grupos a los cuales pertenecen dichos objetos o individuos. Es una técnica que nos permite comprobar hasta qué punto las variables independientes consideradas en la investigación clasifican correctamente a los sujetos u objetos. Se muestran y explican los principales elementos que se relacionan con el procedimiento para llevar a cabo el análisis discriminante y su aplicación utilizando el paquete estadístico SPSS, versión 18, para el desarrollo del modelo estadístico, las condiciones para la aplicación del análisis, la estimación e interpretación de las funciones discriminantes, los métodos de clasificación y la validación de los resultados.

Modelo de comportamiento en sevicio para estructuras reticulares de hormigón. Parte I: Descripción del modelo y ejemplos de aplicación

A general method for instantaneous and time-dependent serviceability analysis of plane concrete frames is presented. The methodology is based in an extension of the classic matrix formulation for bars. The main aspects influencing the behaviour of the structural concrete are considered: cracking, creep, shrinkage or prestress losses. To simulate the effect of cracking a smeared model (developed in Part II) based on the modification of the tensile branch of the concrete stress-strain relationship is adopted. The general approach considered permits the application to different materials and constitutive laws. Sequential construction (sectional and structural), incorporation of reinforcement, consideration of the loads history; placing and removing shores, and restraining or releasing in boundary conditions are considered. Some examples are included to highlight the capabilities of the model

Modelo de comportamiento en sevicio para estructuras reticulares de hormigón. Parte II: Diagrama tensión- deformación para elementos fisurados a flexión

This paper presents a general methodology to adjust the coefficients defining a constitutive law for tensioned concrete to simulate, under serviceability conditions, the behaviour of cracked concrete members subjected to simple and combined bending (with axial compressive forces). Although experimental results could be used to obtain the mentioned objective, in this work the models proposed by the CEB are utilized, in order to obtain more general conclusions not subjected to specific experiments. A numerical parametrical analysis is performed to determine the more significant variables influencing the Tension-Stiffening effect and simplified closed-form solution is derived for rectangular sections. A systematic study has been to verify the ability of the method to simulate long-term effects. The methodology can be implemented as an initial process in a general model of analysis like the one presented in Part I

Markov chain montecarlo method applied to rounding zeros of compositional data: first approach

As stated in Aitchison (1986), a proper study of relative variation in a compositional data set should be based on logratios, and dealing with logratios excludes dealing with zeros. Nevertheless, it is clear that zero observations might be present in real data sets, either because the corresponding part is completely absent –essential zeros– or because it is below detection limit –rounded zeros. Because the second kind of zeros is usually understood as “a trace too small to measure”, it seems reasonable to replace them by a suitable small value, and this has been the traditional approach. As stated, e.g. by Tauber (1999) and by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000), the principal problem in compositional data analysis is related to rounded zeros. One should be careful to use a replacement strategy that does not seriously distort the general structure of the data. In particular, the covariance structure of the involved parts –and thus the metric properties– should be preserved, as otherwise further analysis on subpopulations could be misleading. Following this point of view, a non-parametric imputation method is introduced in Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000). This method is analyzed in depth by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2003) where it is shown that the theoretical drawbacks of the additive zero replacement method proposed in Aitchison (1986) can be overcome using a new multiplicative approach on the non-zero parts of a composition. The new approach has reasonable properties from a compositional point of view. In particular, it is “natural” in the sense that it recovers the “true” composition if replacement values are identical to the missing values, and it is coherent with the basic operations on the simplex. This coherence implies that the covariance structure of subcompositions with no zeros is preserved. As a generalization of the multiplicative replacement, in the same paper a substitution method for missing values on compositional data sets is introduced

Analysis of Pleistocene paleodrainage evolution in the Po Basin (Italy) by multivariate statistical techniques

In order to obtain a high-resolution Pleistocene stratigraphy, eleven continuously cored boreholes, 100 to 220m deep were drilled in the northern part of the Po Plain by Regione Lombardia in the last five years. Quantitative provenance analysis (QPA, Weltje and von Eynatten, 2004) of Pleistocene sands was carried out by using multivariate statistical analysis (principal component analysis, PCA, and similarity analysis) on an integrated data set, including high-resolution bulk petrography and heavy-mineral analyses on Pleistocene sands and of 250 major and minor modern rivers draining the southern flank of the Alps from West to East (Garzanti et al, 2004; 2006). Prior to the onset of major Alpine glaciations, metamorphic and quartzofeldspathic detritus from the Western and Central Alps was carried from the axial belt to the Po basin longitudinally parallel to the SouthAlpine belt by a trunk river (Vezzoli and Garzanti, 2008). This scenario rapidly changed during the marine isotope stage 22 (0.87 Ma), with the onset of the first major Pleistocene glaciation in the Alps (Muttoni et al, 2003). PCA and similarity analysis from core samples show that the longitudinal trunk river at this time was shifted southward by the rapid southward and westward progradation of transverse alluvial river systems fed from the Central and Southern Alps. Sediments were transported southward by braided river systems as well as glacial sediments transported by Alpine valley glaciers invaded the alluvial plain. Kew words: Detrital modes; Modern sands; Provenance; Principal Components Analysis; Similarity, Canberra Distance; palaeodrainage

Clustering compositional data trajectories

Our essay aims at studying suitable statistical methods for the clustering of compositional data in situations where observations are constituted by trajectories of compositional data, that is, by sequences of composition measurements along a domain. Observed trajectories are known as “functional data” and several methods have been proposed for their analysis. In particular, methods for clustering functional data, known as Functional Cluster Analysis (FCA), have been applied by practitioners and scientists in many fields. To our knowledge, FCA techniques have not been extended to cope with the problem of clustering compositional data trajectories. In order to extend FCA techniques to the analysis of compositional data, FCA clustering techniques have to be adapted by using a suitable compositional algebra. The present work centres on the following question: given a sample of compositional data trajectories, how can we formulate a segmentation procedure giving homogeneous classes? To address this problem we follow the steps described below. First of all we adapt the well-known spline smoothing techniques in order to cope with the smoothing of compositional data trajectories. In fact, an observed curve can be thought of as the sum of a smooth part plus some noise due to measurement errors. Spline smoothing techniques are used to isolate the smooth part of the trajectory: clustering algorithms are then applied to these smooth curves. The second step consists in building suitable metrics for measuring the dissimilarity between trajectories: we propose a metric that accounts for difference in both shape and level, and a metric accounting for differences in shape only. A simulation study is performed in order to evaluate the proposed methodologies, using both hierarchical and partitional clustering algorithm. The quality of the obtained results is assessed by means of several indices

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

A systematic and feasible method for computing nuclear contributions to electrical properties of polyatomic molecules

An analytic method to evaluate nuclear contributions to electrical properties of polyatomic molecules is presented. Such contributions control changes induced by an electric field on equilibrium geometry (nuclear relaxation contribution) and vibrational motion (vibrational contribution) of a molecular system. Expressions to compute the nuclear contributions have been derived from a power series expansion of the potential energy. These contributions to the electrical properties are given in terms of energy derivatives with respect to normal coordinates, electric field intensity or both. Only one calculation of such derivatives at the field-free equilibrium geometry is required. To show the useful efficiency of the analytical evaluation of electrical properties (the so-called AEEP method), results for calculations on water and pyridine at the SCF/TZ2P and the MP2/TZ2P levels of theory are reported. The results obtained are compared with previous theoretical calculations and with experimental values

Estimation of electronic coupling in π -stacked donor-bridge-acceptor systems: correction of the two-state model

Comparison of donor-acceptor electronic couplings calculated within two-state and three-state models suggests that the two-state treatment can provide unreliable estimates of Vda because of neglecting the multistate effects. We show that in most cases accurate values of the electronic coupling in a π stack, where donor and acceptor are separated by a bridging unit, can be obtained as Ṽ da = (E2 - E1) μ12 Rda + (2 E3 - E1 - E2) 2 μ13 μ23 Rda2, where E1, E2, and E3 are adiabatic energies of the ground, charge-transfer, and bridge states, respectively, μij is the transition dipole moments between the states i and j, and Rda is the distance between the planes of donor and acceptor. In this expression based on the generalized Mulliken-Hush approach, the first term corresponds to the coupling derived within a two-state model, whereas the second term is the superexchange correction accounting for the bridge effect. The formula is extended to bridges consisting of several subunits. The influence of the donor-acceptor energy mismatch on the excess charge distribution, adiabatic dipole and transition moments, and electronic couplings is examined. A diagnostic is developed to determine whether the two-state approach can be applied. Based on numerical results, we showed that the superexchange correction considerably improves estimates of the donor-acceptor coupling derived within a two-state approach. In most cases when the two-state scheme fails, the formula gives reliable results which are in good agreement (within 5%) with the data of the three-state generalized Mulliken-Hush model

Ecologia demogràfica de la perca americana (Micropterus salmoides) a l'Estany de Banyoles

L'Estany de Banyoles, sistema peculiar tant des del punt de vista de la seva formació geològica com de les seves característiques limnològiques, conté actualment una comunitat de peixos profundament modificada respecte de la comunitat original. La perca americana (Micropterus salmoides), introduïda a finals dels anys seixanta del segle XX, és avui una de les espècies dominants en aquesta comunitat, i ocupa sobretot l'hàbitat litoral de l'Estany. Es tracta d'una espècie molt ben estudiada a Nord Amèrica des de diverses disciplines de la biologia i des de fa diverses dècades, cosa que ha comportat que actualment es disposi d'un gran volum d'informació sobre ella. Amb tot, fora del seu continent d'origen ha rebut poca atenció, malgrat l'amplia expansió que ha experimentat arreu del món. En aquesta tesi doctoral s'han abordat, amb un enfocament descriptiu, aspectes fins ara desconeguts per a l'espècie a l'Estany de Banyoles, a la península ibèrica i fins i tot a Europa. Concretament, se n'ha analitzat la condició, el creixement i la demografia, així com les seves variacions temporals. Amb aquesta finalitat, s'ha dissenyat un mostreig composat de deu campanyes de pesca intensives més alguns petits mostrejos addicionals intercalats, mostreig que s'ha allargat des del juliol del 1997 i fins el novembre del 1999. La captura dels exemplars s'ha realitzat mitjançant una tècnica de pesca elèctrica amb una embarcació posada a punt expressament per a aquest estudi, la qual s'ha mostrat considerablement eficient malgrat les dificultats que ofereix el medi. S'ha realitzat un mostreig de marcatge-recaptura basat en la mutilació d'aletes i, en alguns casos, en el marcatge amb pintura acrílica. Només en la darrera campanya (novembre del 1999) s'ha sacrificat una part important de les captures a fi de retirar-ne els otòlits per a la determinació de l'edat. Pel que fa a l'anàlisi de les dades, s'ha aplicat un ampli ventall de mètodes i models per a cada un dels aspectes estudiats, a fi de contrastar-ne els resultats i validar-ne la seva fiabilitat. En el cas de la condició, s'han aplicat mètodes d'anàlisi de la covariància (ANCOVA) i altres mètodes anàlegs, així com, paral·lelament, regressions i anàlisis derivades a partir de la relació longitud-pes. En l'estudi del creixement, s'han realitzat ajustaments de diversos models mitjançant regressions sobre dades de mida a l'edat i sobre dades d'increments de mida observats per interval de temps. També s'han aplicat anàlisis de freqüències de longitud, i, finalment, s'han aplicat mètodes de retrocàlcul a partir dels increments anuals del radi observats en els otòlits. Finalment, en el cas de l'estudi de la demografia, s'han aplicat models de marcatge-recaptura per a l'estimació de la grandària poblacional i de la supervivència, i, a més, s'han ajustat diversos models continus de supervivència sobre aquestes estimacions prèvies. També s'han estimat les capturabilitats associades a la nova tècnica de captura. Per una altra banda, s'ha implementat i realitzat un mostreig sobre la població de pescadors esportius de l'Estany encarat a determinar, bàsicament, la pressió de pesca a què es veu sotmesa l'espècie. Els resultats mostren sobretot una alta estabilitat interanual en tots els aspectes estudiats, que s'explica per l'estabilitat ambiental que, al seu torn, és característica d'aquest ecosistema lacustre. Això reverteix en una longevitat màxima observada que iguala la màxima descrita a la literatura per a l'espècie. Alhora, també s'han descrit fortes oscil·lacions estacionals tant en la condició, com en el creixement, com també en la supervivència, les quals, però, presenten certes diferències en la seva temporalitat, cosa que indica una certa diferenciació en els factors que les regulen.

Desenvolupament de metodologia analítica per a la determinació de glifosat i adjuvants

El glifosat és un herbicida no selectiu àmpliament utilitzat en cultius agrícoles i zones no agrícoles que requereixen un control de la vegetació. Tot i la baixa toxicitat del glifosat en si mateix, la seva aplicació es realitza emprant formulacions comercials que contenen substàncies auxiliars que milloren l'eficàcia de l'herbicida, però modifiquen la toxicitat del producte final.Aquest treball s'ha centrat en la determinació de components de formulacions comercials d'herbicides i el seu comportament en el sòl. La primera part del treball s'ha dedicat al desenvolupament de mètodes d'anàlisi per a la determinació de glifosat i el seu principal metabòlit, AMPA, en matrius aquoses i sòls. A continuació s'ha desenvolupat un mètode d'anàlisi per les substàncies auxiliars addicionades freqüentment en les formulacions comercials, tensioactius no iònics. Finalment s'ha estudiat la mobilitat de l'herbicida en sòls de la zona mediterrània, incidint en els processos d'adsorció i degradació. També s'ha comparat el comportament de diversos herbicides comercials que contenen glifosat.

Modelització numèrica d’un sistema de descàrrega d’una columna de la catedral d Girona mitjançant el mètode dels elements finits (anàlisi no lineal)

Aquest projecte consisteix en aplicar el càlcul no lineal en la modelització volumètricanumèrica de l’estructura del sistema de descàrrega d’una columna del claustre de lacatedral de Girona mitjançant el mètode dels elements finits. A la Universitat de Gironas’ha fet diferents estudis del claustre de la catedral de Girona però sempre simulant uncomportament lineal de les característiques dels materials. El programa utilitzat és la versió docent del programa ANSYS disponible al Dept.d’EMCI i l’element emprat ha sigut el SOLID65. Aquest element permet introduircaracterístiques de no linealitat en els models i és adequat per a anàlisi no lineald’elements com la pedra de Girona

Anàlisi d'imatges estereoscòpiques utilitzant mètodes de baixa càrrega computacional

 Mitjançant imatges estereoscòpiques es poden detectar la posició respecte dela càmera dels objectes que apareixen en una escena. A partir de lesdiferències entre les imatges captades pels dos objectius es pot determinar laprofunditat dels objectes. Existeixen diversitat de tècniques de visió artificialque permeten calcular la localització dels objectes, habitualment amb l’objectiude reconstruir l’escena en 3D. Aquestes tècniques necessiten una gran càrregacomputacional, ja que utilitzen mètodes de comparació bidimensionals, i pertant, no es poden utilitzar per aplicacions en temps real.En aquest treball proposem un nou mètode d’anàlisi de les imatgesestereoscòpiques que ens permeti obtenir la profunditat dels objectes d’unaescena amb uns resultats acceptables. Aquest nou mètode es basa entransformar la informació bidimensional de la imatge en una informacióunidimensional per tal de poder fer la comparació de les imatges amb un baixcost computacional, i dels resultats de la comparació extreure’n la profunditatdels objectes dins l’escena. Això ha de permetre, per exemple, que aquestmètode es pugui implementar en un dispositiu autònom i li permeti realitzaroperacions de guiatge a través d’espais interiors i exteriors.