Análise modal baseada apenas na resposta: decomposição no domínio da frequência

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Microespectroscopia infravermelha por transformada de fourier: identificação e discriminação de cepas clínicas de Candida albicans e Candida glabrata

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

Desenvolvimento de método analítico para determinação do teor de etanol em óleo lubrificante usado proveniente de motores ciclo Otto

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

Análise e síntese de um processador digital wavelet

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

Influência de estruturas geológicas bidimensionais no campo geoeletromagnético na presença do eletrojato equatorial

A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE₁ e TE₂ e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θ_h entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θ_h é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θ_h. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θ_h se aproxima de 90°, os campos primários tornam-se variáveis para valores de θ_p diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE₂ e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10^-2 Hz. Para frequências da ordem de 10^-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE₂ e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10^-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10^-3 Hz) o acoplamento entre os modos TE₂ e TM não pode ser ignorado.

Aplicação de deconvolução homomórfica a dados sísmicos

Um registro sísmico é frequentemente representado como a convolução de um pulso-fonte com a resposta do meio ao impulso, relacionada ao caminho da propagação. O processo de separação destes dois componentes da convolução é denominado deconvolução. Existe uma variedade de aproximações para o desenvolvimento de uma deconvolução. Uma das mais comuns é o uso da filtragem linear inversa, ou seja, o processamento do sinal composto, através de um filtro linear, cuja resposta de frequência é a recíproca da transformada de Fourier de um dos componentes do sinal. Obviamente, a fim de usarmos a filtragem inversa, tais componentes devem ser conhecidas ou estimadas. Neste trabalho, tratamos da aplicação a sinais sísmicos, de uma técnica de deconvolução não linear, proposta por Oppenheim (1965), a qual utiliza a teoria de uma classe de sistemas não lineares, que satisfazem um princípio generalizado de superposição, denominados de sistemas homomórficos. Tais sistemas são particularmente úteis na separação de sinais que estão combinados através da operação de convolução. O algoritmo da deconvolução homomórfica transforma o processo de convolução em uma superposição aditiva de seus componentes, com o resultado de que partes simples podem ser separadas mais facilmente. Esta classe de técnicas de filtragem representa uma generalização dos problemas de filtragem linear. O presente método oferece a considerável vantagem de que não é necessário fazer qualquer suposição prévia sobre a natureza do pulso sísmico fonte, ou da resposta do meio ao impulso, não requerendo assim, as considerações usuais de que o pulso seja de fase-mínima e que a distribuição dos impulsos seja aleatória, embora a qualidade dos resultados obtidos pela análise homomórfica seja muito sensível à razão sinal/ruído, como demonstrado.

Filtragem de sinais de interferência em espectros dinâmicos solares através de suas componentes de Fourier

Through the data acquisition system of the instrument Brazilian Solar Spectroscope (BSS) at INPE, solar observations in the decimetric radio wave band (1000-2500 MHz) are regularly made. This data is showed as dynamic spectra using the software BSSView created for this purpose. The process of data acquisition can be influenced by various sources, dificulting the resulting dynamic spectrum analysis. The objective of this work is to create a computational routine that eliminates dynamic components of the spectrum attributed to interfering signals and integrate it into BSSView. It was done a preliminary study on the programming language Interactive Data Language (IDL), in which the BSSView was developed, and the Fourier transform, that is required for the application of the filter

Detecção de danos em estruturas guiadas usando ondas de alta frequência

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

Analysis Methodology for Vortex-Induced Motion (VIM) of a Monocolumn Platform Applying the Hilbert-Huang Transform Method

Vortex-induced motion (VIM) is a highly nonlinear dynamic phenomenon. Usual spectral analysis methods, using the Fourier transform, rely on the hypotheses of linear and stationary dynamics. A method to treat nonstationary signals that emerge from nonlinear systems is denoted Hilbert-Huang transform (HHT) method. The development of an analysis methodology to study the VIM of a monocolumn production, storage, and offloading system using HHT is presented. The purposes of the present methodology are to improve the statistics analysis of VIM. The results showed to be comparable to results obtained from a traditional analysis (mean of the 10% highest peaks) particularly for the motions in the transverse direction, although the difference between the results from the traditional analysis for the motions in the in-line direction showed a difference of around 25%. The results from the HHT analysis are more reliable than the traditional ones, owing to the larger number of points to calculate the statistics characteristics. These results may be used to design risers and mooring lines, as well as to obtain VIM parameters to calibrate numerical predictions. [DOI: 10.1115/1.4003493]

Fractal descriptors based on Fourier spectrum applied to texture analysis

This work proposes the development and study of a novel technique lot the generation of fractal descriptors used in texture analysis. The novel descriptors are obtained from a multiscale transform applied to the Fourier technique of fractal dimension calculus. The power spectrum of the Fourier transform of the image is plotted against the frequency in a log-log scale and a multiscale transform is applied to this curve. The obtained values are taken as the fractal descriptors of the image. The validation of the proposal is performed by the use of the descriptors for the classification of a dataset of texture images whose real classes are previously known. The classification precision is compared to other fractal descriptors known in the literature. The results confirm the efficiency of the proposed method. (C) 2012 Elsevier B.V. All rights reserved.

Tannin-phenolic resins: Synthesis, characterization, and application as matrix in biobased composites reinforced with sisal fibers

A tannin-phenolic resin (40 wt% of tannin, characterized by H-1 nuclear magnetic resonance (NMR) and C-13 NMR, Fourier transform infrared, thermogravimetry, differential scanning calorimetry) was used to prepare composites reinforced with sisal fibers (30-70 wt%). Inverse gas chromatography results showed that the sisal fibers and the tannin-phenolic thermoset have close values of the dispersive component and also have predominance of acid sites (acid character) at the surface, confirming the favoring of interaction between the sisal fibers and the tannin-phenolic matrix at the interface. The Izod impact strength increased up to 50 wt% of sisal fibers. This composite also showed high storage modulus, and the lower loss modulus, confirming its good fiber/matrix interface, also observed by SEM images. A composite with good properties was prepared from high content of raw material obtained from renewable sources (40 wt% of tannin substituted the phenol in the preparation of the matrix and 50 wt% of matrix was replaced by sisal fibers). (C) 2012 Elsevier Ltd. All rights reserved.

Hilbert transform

The Hilbert transform is an important tool in both pure and applied mathematics. It is largely used in the field of signal processing. Lately has been used in mathematical finance as the fast Hilbert transform method is an efficient and accurate algorithm for pricing discretely monitored barrier and Bermudan style options. The purpose of this report is to show the basic properties of the Hilbert transform and to check the domain of definition of this operator.

La trasformata veloce di Fourier (FFT): analisi e implementazione in C++

La trasformata di Fourier (FT) è uno strumento molto potente implementato, oggi, in un enorme numero di tecnologie. Il suo primo esempio di applicazione fu proprio il campionamento e la digitalizzazione di segnali analogici. Nel tempo l'utilizzo della FT è stato ampliato a più orizzonti in ambito digitale, basti pensare che il formato di compressione '.jpg' utilizza una FT bidimensionale, mentre uno degli ultimi esempi di applicazione si ha nell'imaging digitale in ambito medico (risonanza magnetica nucleare, tomografia assiale computerizzata TAC ecc...). Nonostante gli utilizzi della FT siano molto diversificati il suo basilare funzionamento non è mai cambiato: essa non fa altro che modificare il dominio di una funzione del tempo (un segnale) in un dominio delle frequenze, permettendo così lo studio della composizione in termini di frequenza, ampiezza e fase del segnale stesso. Parallelamente all'evoluzione in termini di applicazioni si è sviluppato uno studio volto a migliorare e ottimizzare la computazione della stessa, data l'esponenziale crescita del suo utilizzo. In questa trattazione si vuole analizzare uno degli algoritmi di ottimizzazione più celebri e utilizzati in tal senso: la trasformata veloce di Fourier (Fast Fourier Transformation o FFT). Si delineeranno quindi le caratteristiche salienti della FT, e verrà introdotto l'algoritmo di computazione tramite linguaggio C++ dedicando particolare attenzione ai limiti di applicazione di tale algoritmo e a come poter modificare opportunamente la nostra funzione in modo da ricondurci entro i limiti di validità.

DCT Type-III for Multicarrier Modulation

In this paper we propose the use of Discrete Cosine Transform Type-III (DCT3) for multicarrier modulation. There are two DCT3 (even and odd) and, for each of them, we derive the expressions for both prefix and suffix to be appended into each data symbol to be transmitted. Moreover, DCT3 are closely related to the corresponding inverse DCT Type-II even and odd. Furthermore, we give explicit expressions for the 1-tap per subcarrier equalizers that must be implemented at the receiver to perform the channel equalization in the frequency-domain. As a result, the proposed DCT3-based multicarrier modulator can be used as an alternative to DFT-based systems to perform Orthogonal Frequency-Division Multiplexing or Discrete Multitone Modulation

Accelerating FPGA-based evolution of wavelet transform filters by optimized task scheduling

Adaptive embedded systems are required in various applications. This work addresses these needs in the area of adaptive image compression in FPGA devices. A simplified version of an evolution strategy is utilized to optimize wavelet filters of a Discrete Wavelet Transform algorithm. We propose an adaptive image compression system in FPGA where optimized memory architecture, parallel processing and optimized task scheduling allow reducing the time of evolution. The proposed solution has been extensively evaluated in terms of the quality of compression as well as the processing time. The proposed architecture reduces the time of evolution by 44% compared to our previous reports while maintaining the quality of compression unchanged with respect to existing implementations. The system is able to find an optimized set of wavelet filters in less than 2 min whenever the input type of data changes.