Optimización de aperturas mediante técnicas de algoritmos evolutivos

Los sistemas de imagen por ultrasonidos son hoy una herramienta indispensable en aplicaciones de diagnóstico en medicina y son cada vez más utilizados en aplicaciones industriales en el área de ensayos no destructivos. El array es el elemento primario de estos sistemas y su diseño determina las características de los haces que se pueden construir (forma y tamaño del lóbulo principal, de los lóbulos secundarios y de rejilla, etc.), condicionando la calidad de las imágenes que pueden conseguirse. En arrays regulares la distancia máxima entre elementos se establece en media longitud de onda para evitar la formación de artefactos. Al mismo tiempo, la resolución en la imagen de los objetos presentes en la escena aumenta con el tamaño total de la apertura, por lo que una pequeña mejora en la calidad de la imagen se traduce en un aumento significativo del número de elementos del transductor. Esto tiene, entre otras, las siguientes consecuencias: Problemas de fabricación de los arrays por la gran densidad de conexiones (téngase en cuenta que en aplicaciones típicas de imagen médica, el valor de la longitud de onda es de décimas de milímetro) Baja relación señal/ruido y, en consecuencia, bajo rango dinámico de las señales por el reducido tamaño de los elementos. Complejidad de los equipos que deben manejar un elevado número de canales independientes. Por ejemplo, se necesitarían 10.000 elementos separados λ 2 para una apertura cuadrada de 50 λ. Una forma sencilla para resolver estos problemas existen alternativas que reducen el número de elementos activos de un array pleno, sacrificando hasta cierto punto la calidad de imagen, la energía emitida, el rango dinámico, el contraste, etc. Nosotros planteamos una estrategia diferente, y es desarrollar una metodología de optimización capaz de hallar de forma sistemática configuraciones de arrays de ultrasonido adaptados a aplicaciones específicas. Para realizar dicha labor proponemos el uso de los algoritmos evolutivos para buscar y seleccionar en el espacio de configuraciones de arrays aquellas que mejor se adaptan a los requisitos fijados por cada aplicación. En la memoria se trata el problema de la codificación de las configuraciones de arrays para que puedan ser utilizados como individuos de la población sobre la que van a actuar los algoritmos evolutivos. También se aborda la definición de funciones de idoneidad que permitan realizar comparaciones entre dichas configuraciones de acuerdo con los requisitos y restricciones de cada problema de diseño. Finalmente, se propone emplear el algoritmo multiobjetivo NSGA II como herramienta primaria de optimización y, a continuación, utilizar algoritmos mono-objetivo tipo Simulated Annealing para seleccionar y retinar las soluciones proporcionadas por el NSGA II. Muchas de las funciones de idoneidad que definen las características deseadas del array a diseñar se calculan partir de uno o más patrones de radiación generados por cada solución candidata. La obtención de estos patrones con los métodos habituales de simulación de campo acústico en banda ancha requiere tiempos de cálculo muy grandes que pueden hacer inviable el proceso de optimización con algoritmos evolutivos en la práctica. Como solución, se propone un método de cálculo en banda estrecha que reduce en, al menos, un orden de magnitud el tiempo de cálculo necesario Finalmente se presentan una serie de ejemplos, con arrays lineales y bidimensionales, para validar la metodología de diseño propuesta comparando experimentalmente las características reales de los diseños construidos con las predicciones del método de optimización. ABSTRACT Currently, the ultrasound imaging system is one of the powerful tools in medical diagnostic and non-destructive testing for industrial applications. Ultrasonic arrays design determines the beam characteristics (main and secondary lobes, beam pattern, etc...) which assist to enhance the image resolution. The maximum distance between the elements of the array should be the half of the wavelength to avoid the formation of grating lobes. At the same time, the image resolution of the target in the region of interest increases with the aperture size. Consequently, the larger number of elements in arrays assures the better image quality but this improvement contains the following drawbacks: Difficulties in the arrays manufacturing due to the large connection density. Low noise to signal ratio. Complexity of the ultrasonic system to handle large number of channels. The easiest way to resolve these issues is to reduce the number of active elements in full arrays, but on the other hand the image quality, dynamic range, contrast, etc, are compromised by this solutions In this thesis, an optimization methodology able to find ultrasound array configurations adapted for specific applications is presented. The evolutionary algorithms are used to obtain the ideal arrays among the existing configurations. This work addressed problems such as: the codification of ultrasound arrays to be interpreted as individuals in the evolutionary algorithm population and the fitness function and constraints, which will assess the behaviour of individuals. Therefore, it is proposed to use the multi-objective algorithm NSGA-II as a primary optimization tool, and then use the mono-objective Simulated Annealing algorithm to select and refine the solutions provided by the NSGA I I . The acoustic field is calculated many times for each individual and in every generation for every fitness functions. An acoustic narrow band field simulator, where the number of operations is reduced, this ensures a quick calculation of the acoustic field to reduce the expensive computing time required by these functions we have employed. Finally a set of examples are presented in order to validate our proposed design methodology, using linear and bidimensional arrays where the actual characteristics of the design are compared with the predictions of the optimization methodology.

How the owl resolves auditory coding ambiguity

The barn owl (Tyto alba) uses interaural time difference (ITD) cues to localize sounds in the horizontal plane. Low-order binaural auditory neurons with sharp frequency tuning act as narrow-band coincidence detectors; such neurons respond equally well to sounds with a particular ITD and its phase equivalents and are said to be phase ambiguous. Higher-order neurons with broad frequency tuning are unambiguously selective for single ITDs in response to broad-band sounds and show little or no response to phase equivalents. Selectivity for single ITDs is thought to arise from the convergence of parallel, narrow-band frequency channels that originate in the cochlea. ITD tuning to variable bandwidth stimuli was measured in higher-order neurons of the owl’s inferior colliculus to examine the rules that govern the relationship between frequency channel convergence and the resolution of phase ambiguity. Ambiguity decreased as stimulus bandwidth increased, reaching a minimum at 2–3 kHz. Two independent mechanisms appear to contribute to the elimination of ambiguity: one suppressive and one facilitative. The integration of information carried by parallel, distributed processing channels is a common theme of sensory processing that spans both modality and species boundaries. The principles underlying the resolution of phase ambiguity and frequency channel convergence in the owl may have implications for other sensory systems, such as electrolocation in electric fish and the computation of binocular disparity in the avian and mammalian visual systems.

Cellular mechanisms for resolving phase ambiguity in the owl's inferior colliculus

Both mammals and birds use the interaural time difference (ITD) for localization of sound in the horizontal plane. They may localize either real or phantom sound sources, when the signal consists of a narrow frequency band. This ambiguity does not occur with broadband signals. A plot of impulse rates or amplitude of excitatory postsynaptic potentials against ITDs (ITD curve) consists of peaks and troughs. In the external nucleus (ICX) of the owl's inferior colliculus, ITD curves show multiple peaks when the signal is narrow-band, such as tones. Of these peaks, one occurs at ITDi, which is independent of frequency, and others at ITDi ± T, where T is the tonal period. The ITD curve of the same neuron shows a large peak (main peak) at ITDi and no or small peaks (side peaks) at ITDi ± T, when the signal is broadband. ITD curves for postsynaptic potentials indicate that ICX neurons integrate the results of binaural cross-correlation in different frequency bands. However, the difference between the main and side peaks is small. ICX neurons further enhance this difference in the process of converting membrane potentials to impulse rates. Inhibition also appears to augment the difference between the main and side peaks.

A fast marching level set method for monotonically advancing fronts.

A fast marching level set method is presented for monotonically advancing fronts, which leads to an extremely fast scheme for solving the Eikonal equation. Level set methods are numerical techniques for computing the position of propagating fronts. They rely on an initial value partial differential equation for a propagating level set function and use techniques borrowed from hyperbolic conservation laws. Topological changes, corner and cusp development, and accurate determination of geometric properties such as curvature and normal direction are naturally obtained in this setting. This paper describes a particular case of such methods for interfaces whose speed depends only on local position. The technique works by coupling work on entropy conditions for interface motion, the theory of viscosity solutions for Hamilton-Jacobi equations, and fast adaptive narrow band level set methods. The technique is applicable to a variety of problems, including shape-from-shading problems, lithographic development calculations in microchip manufacturing, and arrival time problems in control theory.

Emission-line galaxy surveys with GTC: The Hα Universe

Emission line galaxies are the most easily detected and studied objects in the high redshift Universe. They are being used to trace the evolution of critical observables of the Universe such as Star Formation Rate densities, starburst properties and abundances. Most of the research is being done using [OII]3727 and UV lines, but Hα is still one the best tracers for Star Formation Rate and physical properties of current star-forming galaxies. As a complementary contribution to studies of galaxy evolution, our team has focused into a long-term project to study the population of Hα-selected star-forming galaxies of the Universe at different redshifts. In 1995 we first determined the local Hα luminosity function, and from it the Star Formation Rate density (SFRd) of the local Universe. We then, using narrow-band imaging in the optical, extended this measurement to z ≃ 0.24 and z ≃ 0.4. Working in the near-infrared, GTC will be a very powerful tool to study the evolution of the Hα emission-line galaxy populations at different redshifts. We will both quantify the SFRd evolution and characterize the star-forming galaxy populations by directly comparing the same observables at all redshifts up to z ≃ 2.5.

Estudo e aplicação de guias de ondas integrados ao substrato em frequências de micro-ondas.

Esse trabalho de pesquisa apresenta um estudo detalhado sobre guias de ondas integrados ao substrato (SIW) operando em frequências de micro-ondas com base na teoria de guias de ondas retangulares (RWG). O estudo sobre guias SIW associa equações apresentadas na literatura e utiliza simulações eletromagnéticas para desenvolver um procedimento de projeto bem definido. É considerada a integração entre guias SIW e linhas de transmissão de microfita, projetando-se transições entre essas duas estruturas com o propósito de prover casamento de impedância e de viabilizar a caracterização em frequências de micro-ondas. São apresentadas considerações sobre processos de fabricação de circuitos SIW em substratos constituídos por laminados de alta frequência. Uma vez estabelecidos os procedimentos de fabricação e os critérios de projeto, a tecnologia SIW é aplicada ao projeto de três guias de ondas SIW nas bandas S e X, que foram fabricados empregando laminados de alta-frequência. Foram projetados dois filtros SIW passa-faixa empregando cavidades ressonantes e postes metálicos indutivos. Os dois filtros operam na frequência central de 10,61 GHz, sendo que um deles tem banda de passagem de 7,5%e é de 3ª ordem e o outro filtro tem banda de passagem de 15%, sendo de 5ª ordem. Foram realizadas comparações entre o desempenho simulado e experimental das estruturas SIW projetadas. Os resultados de simulações eletromagnéticas e experimentais demonstraram boa concordância. Os projetos em tecnologia SIW apresentados neste trabalho de pesquisa possuem perdas de retorno melhores que 10 dB na banda de operação e perdas por inserção de 1,0 dB a 1,5 dB. É apresentada a análise da sensibilidade do desempenho dos guias de ondas e filtros SIW projetados a desvios dimensionais típicos do processo de fabricação por microfresagem mecânica. Com os resultados experimentais e de simulação foi possível validar os procedimentos de projeto e de fabricação de circuitos SIW operando em frequências de micro-ondas.

Role of potassium orbitals in the metallic behavior of K3picene

Detailed electronic structure calculations of picene clusters doped by potassium modeling the crystalline K3picene structure show that while two electrons are completely transferred from potassium atoms to the lowest-energy unoccupied molecular orbital of pristine picene, the third one remains closely attached to both material components. Multiconfigurational analysis is necessary to show that many structures of almost degenerate total energies compete to define the cluster ground state. Our results prove that the 4s orbital of potassium should be included in any interaction model describing the material. We propose a quarter-filled two-orbital model as the most simple model capable of describing the electronic structure of K-intercalated picene. Precise solutions obtained by a development of the Lanczos method show low-energy electronic excitations involving orbitals located at different positions. Consequently, metallic transport is possible in spite of the clear dominance of interaction over hopping.

The Hi content of compact groups of galaxies

The Hi content of Hickson Compact Groups in the southern hemisphere is measured using data from the Hi Parkes All-Sky Survey (HIPASS), and dedicated observations using the narrow band filter on the Multibeam instrument on the Parkes telescope. The expected Hi mass of these groups was estimated using the luminosity, diameter, and morphological types of the member galaxies, calibrated from published data. Taking careful account of non-detection limits, the results show that the compact group population that has been detected by these observations has an Hi content similar to that of galaxies in the reference field sample. The upper limits for the undetected groups lie within the normal range; improvement of these limits will require a large increase in sensitivity.

A worldwide correlation for exponential bed particle size variation in subaerial aqueous flows

The particle size of the bed sediments in or on many natural streams, alluvial fans, laboratory flumes, irrigation canals and mine waste deltas varies exponentially with distance along the stream. A plot of the available worldwide exponential bed particle size diminution coefficient data against stream length is presented which shows that all the data lie within a single narrow band extending over virtually the whole range of stream lengths and bed sediment particle sizes found on Earth. This correlation applies to both natural and artificial flows with both sand and gravel beds, irrespective of either the solids concentration or whether normal or reverse sorting occurs. This strongly suggests that there are common mechanisms underlying the exponential diminution of bed particles in subaerial aqueous flows of all kinds. Thus existing models of sorting and abrasion applicable to some such flows may be applicable to others. A comparison of exponential laboratory abrasion and field diminution coefficients suggests that abrasion is unlikely to be significant in gravel and sand bed streams shorter than about 10 km to 100 km, and about 500 km, respectively. Copyright (C) 1999 John Wiley & Sons, Ltd.

Induced visual illusions and gamma oscillations in human primary visual cortex

Using magnetoencephalography, we studied the spatiotemporal properties of cortical responses in terms of event-related synchronization and event-related desynchronization to a range of stripe patterns in subjects with no neurological disorders. These stripes are known for their tendency to induce a range of abnormal sensations, such as illusions, nausea, dizziness, headache and attacks of pattern-sensitive epilepsy. The optimal stimulus must have specific physical properties, and maximum abnormalities occur at specific spatial frequency and contrast. Despite individual differences in the severity of discomfort experienced, psychophysical studies have shown that most observers experience some degree of visual anomaly on viewing such patterns. In a separate experiment, subjects reported the incidence of illusions and discomfort to each pattern. We found maximal cortical power in the gamma range (30-60 Hz) confined to the region of the primary visual cortex in response to patterns of 2-4 cycles per degree, peaking at 3 cycles per degree. This coincides with the peak of mean illusions and discomfort, also maximal for patterns of 2-4 cycles per degree. We show that gamma band activity in V1 is a narrow band function of spatial frequency. We hypothesize that the intrinsic properties of gamma oscillations may underlie visual discomfort and play a role in the onset of seizures.

A study of field emission from semiconductors

In an attempt to clarify the behaviour of semi-conductor field emitters the properties of a narrow band gap material were investigated. A retarding potential analyser was built and tested using a tungsten emitter. The energy distribution of electrons emitted from single crystals of lead telluride (band gap 0.3 eV) and gallium phosphide (band gap 2.26 eV) were measured. The halfwidths of the distributions are discussed with respect to the relevant parameters for the materials. Methods of tip preparation had to be developed. The halfwidth of the energy distribution of electrons field emitted from carbon fibres was measured to be 0.21 ± 0.01 eV. A mechanism explaining the long lifetime of the emitters in poor vacuua is proposed.

Estimation of functional connectivity from electromagnetic signals and the amount of empirical data required

An increasing number of neuroimaging studies are concerned with the identification of interactions or statistical dependencies between brain areas. Dependencies between the activities of different brain regions can be quantified with functional connectivity measures such as the cross-correlation coefficient. An important factor limiting the accuracy of such measures is the amount of empirical data available. For event-related protocols, the amount of data also affects the temporal resolution of the analysis. We use analytical expressions to calculate the amount of empirical data needed to establish whether a certain level of dependency is significant when the time series are autocorrelated, as is the case for biological signals. These analytical results are then contrasted with estimates from simulations based on real data recorded with magnetoencephalography during a resting-state paradigm and during the presentation of visual stimuli. Results indicate that, for broadband signals, 50-100 s of data is required to detect a true underlying cross-correlations coefficient of 0.05. This corresponds to a resolution of a few hundred milliseconds for typical event-related recordings. The required time window increases for narrow band signals as frequency decreases. For instance, approximately 3 times as much data is necessary for signals in the alpha band. Important implications can be derived for the design and interpretation of experiments to characterize weak interactions, which are potentially important for brain processing.

Gasification of biomass in a downdraft reactor

The objective of this study was to design, construct, commission and operate a laboratory scale gasifier system that could be used to investigate the parameters that influence the gasification process. The gasifier is of the open-core variety and is fabricated from 7.5 cm bore quartz glass tubing. Gas cleaning is by a centrifugal contacting scrubber, with the product gas being flared. The system employs an on-line dedicated gas analysis system, monitoring the levels of H2, CO, CO2 and CH4 in the product gas. The gas composition data, as well as the gas flowrate, temperatures throughout the system and pressure data is recorded using a BBC microcomputer based data-logging system. Ten runs have been performed using the system of which six were predominantly commissioning runs. The main emphasis in the commissioning runs was placed on the gas clean-up, the product gas cleaning and the reactor bed temperature measurement. The reaction was observed to occur in a narrow band, of about 3 to 5 particle diameters thick. Initially the fuel was pyrolysed, with the volatiles produced being combusted and providing the energy to drive the process, and then the char product was gasified by reaction with the pyrolysis gases. Normally, the gasifier is operated with reaction zone supported on a bed of char, although it has been operated for short periods without a char bed. At steady state the depth of char remains constant, but by adjusting the air inlet rate it has been shown that the depth of char can be increased or decreased. It has been shown that increasing the depth of the char bed effects some improvement in the product gas quality.

Pooling and suppression in human spatial vision

A distinct feature of several recent models of contrast masking is that detecting mechanisms are divisively inhibited by a broadly tuned ‘gain pool’ of narrow-band spatial pattern mechanisms. The contrast gain control provided by this ‘cross-channel’ architecture achieves contrast normalisation of early pattern mechanisms, which is important for keeping them within the non-saturating part of their biological operating characteristic. These models superseded earlier ‘within-channel’ models, which had supposed that masking arose from direct stimulation of the detecting mechanism by the mask. To reveal the extent of masking, I measured the levels produced with large ranges of pattern spatial relationships that have not been explored before. Substantial interactions between channels tuned to different orientations and spatial frequencies were found. Differences in the masking levels produced with single and multiple component mask patterns provided insights into the summation rules within the gain pool. A widely used cross-channel masking model was tested on these data and was found to perform poorly. The model was developed and a version in which linear summation was allowed between all components within the gain pool but with the exception of the self-suppressing route typically provided the best account of the data. Subsequently, an adaptation paradigm was used to probe the processes underlying pooled responses in masking. This delivered less insight into the pooling than the other studies and areas were identified that require investigation for a new unifying model of masking and adaptation. In further experiments, levels of cross-channel masking were found to be greatly influenced by the spatio-temporal tuning of the channels involved. Old masking experiments and ideas relying on within-channel models were re-elevated in terms of contemporary cross-channel models (e.g. estimations of channel bandwidths from orientation masking functions) and this led to different conclusions than those originally arrived at. The investigation of effects with spatio-temporally superimposed patterns is focussed upon throughout this work, though it is shown how these enquiries might be extended to investigate effects across spatial and temporal position.

Call admission control algorithms in OFDM-based wireless multiservice networks

Orthogonal frequency division multiplexing (OFDM) is becoming a fundamental technology in future generation wireless communications. Call admission control is an effective mechanism to guarantee resilient, efficient, and quality-of-service (QoS) services in wireless mobile networks. In this paper, we present several call admission control algorithms for OFDM-based wireless multiservice networks. Call connection requests are differentiated into narrow-band calls and wide-band calls. For either class of calls, the traffic process is characterized as batch arrival since each call may request multiple subcarriers to satisfy its QoS requirement. The batch size is a random variable following a probability mass function (PMF) with realistically maximum value. In addition, the service times for wide-band and narrow-band calls are different. Following this, we perform a tele-traffic queueing analysis for OFDM-based wireless multiservice networks. The formulae for the significant performance metrics call blocking probability and bandwidth utilization are developed. Numerical investigations are presented to demonstrate the interaction between key parameters and performance metrics. The performance tradeoff among different call admission control algorithms is discussed. Moreover, the analytical model has been validated by simulation. The methodology as well as the result provides an efficient tool for planning next-generation OFDM-based broadband wireless access systems.