Distributed Estimation Over an Adaptive Incremental Network Based on the Affine Projection Algorithm

We study the problem of distributed estimation based on the affine projection algorithm (APA), which is developed from Newton`s method for minimizing a cost function. The proposed solution is formulated to ameliorate the limited convergence properties of least-mean-square (LMS) type distributed adaptive filters with colored inputs. The analysis of transient and steady-state performances at each individual node within the network is developed by using a weighted spatial-temporal energy conservation relation and confirmed by computer simulations. The simulation results also verify that the proposed algorithm provides not only a faster convergence rate but also an improved steady-state performance as compared to an LMS-based scheme. In addition, the new approach attains an acceptable misadjustment performance with lower computational and memory cost, provided the number of regressor vectors and filter length parameters are appropriately chosen, as compared to a distributed recursive-least-squares (RLS) based method.

The Cramer-Rao bound for initial conditions estimation of chaotic orbits

We derive the Cramer-Rao Lower Bound (CRLB) for the estimation of initial conditions of noise-embedded orbits produced by general one-dimensional maps. We relate this bound`s asymptotic behavior to the attractor`s Lyapunov number and show numerical examples. These results pave the way for more suitable choices for the chaotic signal generator in some chaotic digital communication systems. (c) 2006 Published by Elsevier Ltd.

DETECTABILITY AND ERROR ESTIMATION IN ORBITAL FITS OF RESONANT EXTRASOLAR PLANETS

We estimate the conditions for detectability of two planets in a 2/1 mean-motion resonance from radial velocity data, as a function of their masses, number of observations and the signal-to-noise ratio. Even for a data set of the order of 100 observations and standard deviations of the order of a few meters per second, we find that Jovian-size resonant planets are difficult to detect if the masses of the planets differ by a factor larger than similar to 4. This is consistent with the present population of real exosystems in the 2/1 commensurability, most of which have resonant pairs with similar minimum masses, and could indicate that many other resonant systems exist, but are currently beyond the detectability limit. Furthermore, we analyze the error distribution in masses and orbital elements of orbital fits from synthetic data sets for resonant planets in the 2/1 commensurability. For various mass ratios and number of data points we find that the eccentricity of the outer planet is systematically overestimated, although the inner planet`s eccentricity suffers a much smaller effect. If the initial conditions correspond to small-amplitude oscillations around stable apsidal corotation resonances, the amplitudes estimated from the orbital fits are biased toward larger amplitudes, in accordance to results found in real resonant extrasolar systems.

Advanced functionality for radio analysis in the Offline software framework of the Pierre Auger Observatory

The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs ""radio-hybrid"" measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluorescence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request. (C) 2011 Elsevier B.V. All rights reserved.

O uso de sondagens dipolo-dipolo em estudos hidrogeológicos e de depósitos de resíduos

Ensaios de sondagem e caminhamento elétrico têm sido realizados globalmente na avaliação de uma grande variedade de problemas em exploração mineral e hidrogeologia, bem como na caracterização de locais com risco de contaminação. A despeito da utilidade de ensaios geoelétricos em estudos hidrogeológicos e na caracterização de depósitos de resíduos desativados, um problema na utilização de sondagens elétricas ocorre em conseqüência da falta de espaço para aquisição de dados nessas áreas. Em áreas urbanas, isto constitui uma limitação importante devido à existência de construções, o que muitas vezes interrompe o desenvolvimento de ensaios de campo com arranjos de comprimentos adequados. Embora investigações rasas tendam a ser eficazes em caracterização geoambiental, a estimativa de parâmetros-chave pode depender da investigação de porções mais profundas. Entretanto, uma profundidade de investigação ligeiramente maior é obtida pelo arranjo dipolo-dipolo em comparação aos arranjos Schlumberger e Wenner de mesmo comprimento. Essa característica do arranjo dipolo-dipolo pode ser observada em vários estudos, mas provavelmente devido à tradição e à relação sinal-ruído, a grande maioria das sondagens elétricas é feita com os arranjos Schlumberger e Wenner.O arranjo dipolo-dipolo é mais utilizado em ensaios 2D. Até onde se sabe, essa maior capacidade de investigação em profundidade do arranjo dipolo-dipolo 1D não foi explorada em estudos geoambientais e nem em estudos hidrogeológicos em áreas urbanas. Este trabalho apresenta resultados de ensaios de campo que ressaltam a utilidade de sondagens dipolo-dipolo na investigação de tais áreas.

PENCALC: a program for penetrance estimation in autosomal dominant diseases

We present a computer program developed for estimating penetrance rates in autosomal dominant diseases by means of family kinship and phenotype information contained within the pedigrees. The program also determines the exact 95% credibility interval for the penetrance estimate. Both executable (PenCalc for Windows) and web versions (PenCalcWeb) of the software are available. The web version enables further calculations, such as heterozygosity probabilities and assessment of offspring risks for all individuals in the pedigrees. Both programs can be accessed and down-loaded freely at the home-page address http://www.ib.usp.br/~otto/software.htm.

Imagem por ressonância magnética na investigação da cabeça de cães

A imagem por ressonância magnética (IRM) é o método de diagnóstico por imagem não invasivo mais sensível para avaliar as partes moles, particularmente o encéfalo, porém trata-se de uma técnica onerosa. O método fundamenta-se no fenômeno da ressonância magnética nuclear que ocorre quando núcleos atômicos com propriedades magnéticas presentes no corpo são submetidos a um campo magnético intenso, sendo posteriormente excitados por energia de radiofrequência e gerando, por sua vez, um sinal de onda de radiofrequência capaz de ser captado por uma antena receptora, passando por um processo matemático, chamado Transformada de Fourier, para posterior formação da imagem. Esse estudo objetivou realizar 10 exames completos da cabeça em cadáveres de cães normais à IRM e confeccionar um Atlas com as estruturas identificadas. As imagens foram adquiridas em um aparelho de ressonância magnética Gyroscan S15/HP Philips com campo magnético de 1,5Tesla. Os cadáveres foram posicionados com a cabeça no interior de uma bobina de cabeça humana e foram submetidos a cortes iniciais sagitais a partir de onde se planejou os cortes transversais e dorsais nas sequências de pulso spin-eco T1, T2 e DP. Em T1 utilizou-se TR=400ms e TE=30ms, T2 utilizou-se TR=2000ms e TE=80ms e na DP utilizou-se TR=2000ms e TE=30ms. A espessura do corte foi de 4mm, o número de médias foi igual a 2, a matriz foi de 256x256, o fator foi igual a 1,0 e o campo de visão foi de 14cm. A duração do exame completo da cabeça foi de 74,5minutos. As imagens obtidas com as sequências utilizadas e com a bobina de cabeça humana foram de boa qualidade. Em T1 a gordura tornou-se hiperintensa e o líquido hipointenso. Em T2 a gordura ficou menos hiperintensa e o líquido hiperintenso. A cortical óssea e o ar foram hipointensos em todas as sequências utilizadas devido a baixa densidade de prótons. A sequência DP mostrou o melhor contraste entre a substância branca e cinzenta quando comparada a T2 e a T1. T2 evidenciou o líquido cefalorraquidiano tornando possível a distinção dos sulcos e giros cerebrais. Através do exame de IRM foi possível, pelo contraste, identificar as estruturas ósseas componentes da arquitetura da região, músculos, grandes vasos venosos e arteriais e estruturas do sistema nervoso central, além de elementos do sistema digestório, respiratório e estruturas dos olhos entre outras. Nesse estudo as IRM adquiridas nas sequências T1, DP e T2 foram complementares para o estudo dos aspectos anatômicos da cabeça de cães demonstrando-os com riqueza de detalhes. O tempo requerido para o exame completo da cabeça é compátivel para uso em animais vivos desde que devidamente anestesiados e controlados. Os resultados obtidos por esse trabalho abrem caminho em nosso meio, para o estudo de animais vivos e para o início da investigação de doenças, principalmente as de origem neurológica, visto ser esta técnica excelente para a visibilização do encéfalo.

Imagem por ressonância magnética: princípios básicos

A ressonância magnética é a propriedade física exibida por núcleos de determinados elementos que, quando submetidos a um campo magnético forte e excitados por ondas de rádio em determinada freqüência (Freqüência de Larmor), emitem rádio sinal, o qual pode ser captado por uma antena e transformado em imagem. A imagem por ressonância magnética (IRM) é o método de diagnóstico por imagem não-invasivo mais sensível para avaliar partes moles, particularmente o encéfalo, porém trata-se de uma técnica onerosa. Ela apresenta grande potencial diagnóstico, poucos efeitos deletérios e muitos benefícios a serem obtidos com o seu uso. Além disso, a IRM fornece informações anatômicas acuradas, imagens em qualquer plano do corpo, bom contraste e resolução espacial e por si só pode sugerir um diagnóstico. Porém, não permite um diagnóstico histológico específico e deve ser interpretada em contexto com outros achados clínicos e patológicos. Esta revisão teve como objetivos mostrar as bases físicas da ressonância magnética e propiciar mais conhecimento aos veterinários.

Load Ratio Estimation Through Striation Height and Spacing Analysis of an Aerospace Al Alloy 7475-T7351

It is well known that striation spacing may be related to the crack growth rate, da/dN, through Paris equation, as well as the maximum and minimum loads under service loading conditions. These loads define the load ratio, R, and are considered impossible to be evaluated from the inter-spacing striations analysis. In this way, this study discusses the methodology proposed by Furukawa to evaluate the maximum and minimum loads based on the experimental fact that the relative height of a striation, H, and the striation spacing, s, are strongly influenced by the load ratio, R. Fatigue tests in C(T) specimens were conducted on SAE 7475-T7351 Al alloy plates at room temperature and the results showed a straightforward correlation between the parameters H, s, and R. Measurements of striation height, H, were performed using scanning electron microscopy and field emission gun (FEG) after sectioning the specimen at a large inclined angle to amplify the height of the striations. The results showed that for increasing R the values of H/s tend to increase. Striation height, striation spacing, and load ratio correlations were obtained, which allows one to estimate service loadings from fatigue fracture surface survey.

Estimation and prediction of parameters and breeding values in soybean using REML/BLUP and Least Squares

The aim of this study was to compare REML/BLUP and Least Square procedures in the prediction and estimation of genetic parameters and breeding values in soybean progenies. F(2:3) and F(4:5) progenies were evaluated in the 2005/06 growing season and the F(2:4) and F(4:6) generations derived thereof were evaluated in 2006/07. These progenies were originated from two semi-early, experimental lines that differ in grain yield. The experiments were conducted in a lattice design and plots consisted of a 2 m row, spaced 0.5 m apart. The trait grain yield per plot was evaluated. It was observed that early selection is more efficient for the discrimination of the best lines from the F(4) generation onwards. No practical differences were observed between the least square and REML/BLUP procedures in the case of the models and simplifications for REML/BLUP used here.

Continental passive microwave-based rainfall estimation algorithm: Application to the Amazon Basin

This paper presents a new statistical algorithm to estimate rainfall over the Amazon Basin region using the Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI). The algorithm relies on empirical relationships derived for different raining-type systems between coincident measurements of surface rainfall rate and 85-GHz polarization-corrected brightness temperature as observed by the precipitation radar (PR) and TMI on board the TRMM satellite. The scheme includes rain/no-rain area delineation (screening) and system-type classification routines for rain retrieval. The algorithm is validated against independent measurements of the TRMM-PR and S-band dual-polarization Doppler radar (S-Pol) surface rainfall data for two different periods. Moreover, the performance of this rainfall estimation technique is evaluated against well-known methods, namely, the TRMM-2A12 [ the Goddard profiling algorithm (GPROF)], the Goddard scattering algorithm (GSCAT), and the National Environmental Satellite, Data, and Information Service (NESDIS) algorithms. The proposed algorithm shows a normalized bias of approximately 23% for both PR and S-Pol ground truth datasets and a mean error of 0.244 mm h(-1) ( PR) and -0.157 mm h(-1)(S-Pol). For rain volume estimates using PR as reference, a correlation coefficient of 0.939 and a normalized bias of 0.039 were found. With respect to rainfall distributions and rain area comparisons, the results showed that the formulation proposed is efficient and compatible with the physics and dynamics of the observed systems over the area of interest. The performance of the other algorithms showed that GSCAT presented low normalized bias for rain areas and rain volume [0.346 ( PR) and 0.361 (S-Pol)], and GPROF showed rainfall distribution similar to that of the PR and S-Pol but with a bimodal distribution. Last, the five algorithms were evaluated during the TRMM-Large-Scale Biosphere-Atmosphere Experiment in Amazonia (LBA) 1999 field campaign to verify the precipitation characteristics observed during the easterly and westerly Amazon wind flow regimes. The proposed algorithm presented a cumulative rainfall distribution similar to the observations during the easterly regime, but it underestimated for the westerly period for rainfall rates above 5 mm h(-1). NESDIS(1) overestimated for both wind regimes but presented the best westerly representation. NESDIS(2), GSCAT, and GPROF underestimated in both regimes, but GPROF was closer to the observations during the easterly flow.

Populations III.1 and III.2 gamma-ray bursts: constraints on the event rate for future radio and X-ray surveys

Aims. We calculate the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive first-generation (Population III; Pop III) stars. The Pop III GRBs could be super-energetic with the isotropic energy up to E(iso) greater than or similar to 10(55-57) erg, providing a unique probe of the high-redshift Universe. Methods. We consider both the so-called Pop III.1 stars (primordial) and Pop III.2 stars (primordial but affected by radiation from other stars). We employ a semi-analytical approach that considers inhomogeneous hydrogen reionization and chemical evolution of the intergalactic medium. Results. We show that Pop III.2 GRBs occur more than 100 times more frequently than Pop III.1 GRBs, and thus should be suitable targets for future GRB missions. Interestingly, our optimistic model predicts an event rate that is already constrained by the current radio transient searches. We expect similar to 10-10(4) radio afterglows above similar to 0.3 mJy on the sky with similar to 1 year variability and mostly without GRBs (orphans), which are detectable by ALMA, EVLA, LOFAR, and SKA, while we expect to observe maximum of N < 20 GRBs per year integrated over at z > 6 for Pop III.2 and N < 0.08 per year integrated over at z > 10 for Pop III.1 with EXIST, and N < 0.2 for Pop III.2 GRBs per year integrated over at z > 6 with Swift.

Theory of small-signal ac response of a dielectric liquid containing two groups of ions

The analysis of Macdonald for electrolytes is generalized to the case in which two groups of ions are present. We assume that the electrolyte can be considered as a dispersion of ions in a dielectric liquid, and that the ionic recombination can be neglected. We present the differential equations governing the ionic redistribution when the liquid is subjected to an external electric field, describing the simultaneous diffusion of the two groups of ions in the presence of their own space charge fields. We investigate the influence of the ions on the impedance spectroscopy of an electrolytic cell. In the analysis, we assume that each group of ions have equal mobility, the electrodes perfectly block and that the adsorption phenomena can be neglected. In this framework, it is shown that the real part of the electrical impedance of the cell has a frequency dependence presenting two plateaux, related to a type of ambipolar and free diffusion coefficients. The importance of the considered problem on the ionic characterization performed by means of the impedance spectroscopy technique was discussed. (c) 2008 American Institute of Physics.

Radio-frequency spectroscopy of (6)Li p-wave molecules: Towards photoemission spectroscopy of a p-wave superfluid

We study rf spectroscopy of a lithium gas with the goal to explore the possibilities for photoemission spectroscopy of a strongly interacting p-wave Fermi gas. Radio-frequency spectra of quasibound p-wave molecules and of free atoms in the vicinity of the p-wave Feshbach resonance located at 159.15G are presented. The spectra are free of detrimental final-state effects. The observed relative magnetic-field shifts of the molecular and atomic resonances confirm earlier measurements realized with direct rf association. Furthermore, evidence of molecule production by adiabatically ramping the magnetic field is observed. Finally, we propose the use of a one-dimensional optical lattice to study anisotropic superfluid gaps as most direct proof of p-wave superfluidity.

Estimation of the thickness and the optical parameters of several stacked thin films using optimization

The reverse engineering problem addressed in the present research consists of estimating the thicknesses and the optical constants of two thin films deposited on a transparent substrate using only transmittance data through the whole stack. No functional dispersion relation assumptions are made on the complex refractive index. Instead, minimal physical constraints are employed, as in previous works of some of the authors where only one film was considered in the retrieval algorithm. To our knowledge this is the first report on the retrieval of the optical constants and the thickness of multiple film structures using only transmittance data that does not make use of dispersion relations. The same methodology may be used if the available data correspond to normal reflectance. The software used in this work is freely available through the PUMA Project web page (http://www.ime.usp.br/similar to egbirgin/puma/). (C) 2008 Optical Society of America