Dinâmica monetária eficiente sob encontros aleatórios: uma classe de métodos numéricos que exploram concavidade

A dificuldade em se caracterizar alocações ou equilíbrios não estacionários é uma das principais explicações para a utilização de conceitos e hipóteses que trivializam a dinâmica da economia. Tal dificuldade é especialmente crítica em Teoria Monetária, em que a dimensionalidade do problema é alta mesmo para modelos muito simples. Neste contexto, o presente trabalho relata a estratégia computacional de implementação do método recursivo proposto por Monteiro e Cavalcanti (2006), o qual permite calcular a sequência ótima (possivelmente não estacionária) de distribuições de moeda em uma extensão do modelo proposto por Kiyotaki e Wright (1989). Três aspectos deste cálculo são enfatizados: (i) a implementação computacional do problema do planejador envolve a escolha de variáveis contínuas e discretas que maximizem uma função não linear e satisfaçam restrições não lineares; (ii) a função objetivo deste problema não é côncava e as restrições não são convexas; e (iii) o conjunto de escolhas admissíveis não é conhecido a priori. O objetivo é documentar as dificuldades envolvidas, as soluções propostas e os métodos e recursos disponíveis para a implementação numérica da caracterização da dinâmica monetária eficiente sob a hipótese de encontros aleatórios.

Análise de métodos numéricos para precificação de opções

Analisa e compara, sob as óticas de acurácia e custo, os principais métodos numéricos (simulação de Monte carlo, métodos de diferenças finitas, e modelos baseados em lattice para precificação de opções financeiras, exóticas e reais. Aponta as situações onde cada método deve ser empregado, e apresenta os algoritmos necessários para implementação dos métodos numéricos na prática.

Uso do método das interfaces coesivas na análise da simulação do comportamento de fratura em materiais frágeis

Este trabalho apresenta uma abordagem numérica da fratura baseada no Método dos Elementos Finitos na qual são considerados materiais frágeis como objeto de estudo. As simulações realizadas são possíveis através de um algoritmo computacional que consiste na implementação de molas entre as faces dos elementos finitos. Esta composição constitui o método denominado Interfaces Coesivas. As molas fornecem o efeito de coesão e é representado por leis constitutivas do tipo bi-linear. Sua descrição é apresentada, bem como sua implementação na análise. Tanto o modo I como o modo II de ruptura são considerados. Com o objetivo de verificar o presente modelo com outros resultados e soluções, são analisados exemplos numéricos e experimentais extraídos da literatura. Tais exemplos tratam da nucleação, propagação e eventuais paradas de trincas em peças estruturais submetidas aos mais diversos carregamentos e condições de contorno. É ainda discutido como os parâmetros de fratura do material interferem no processo de ruptura. Uma tentativa é feita de determinar estes parâmetros, em especial no caso de concreto.

Identificação de sistemas contínuos utilizando o método dos momentos de Poisson

A maioria dos métodos de síntese e sintonia de controladores, bem como métodos de otimização e análise de processos necessitam de um modelo do processo em estudo. A identificação de processos é portanto uma área de grande importância para a engenharia em geral pois permite a obtenção de modelos empíricos dos processos com que nos deparamos de uma forma simples e rápida. Mesmo não utilizando leis da natureza, os modelos empíricos são úteis pois descrevem o comportamento específico de determinado processo. Com o rápido desenvolvimento dos computadores digitais e sua larga aplicação nos sistemas de controle em geral, a identificação de modelos discretos foi amplamente desenvolvida e empregada, entretanto, modelos discretos não são de fácil interpretação como os modelos contínuos pois a maioria dos sistema com que lidamos são de representação contínua. A identificação de modelos contínuos é portanto útil na medida que gera modelos de compreensão mais simples. A presente dissertação estuda a identificação de modelos lineares contínuos a partir de dados amostrados discretamente. O método estudado é o chamado método dos momentos de Poisson. Este método se baseia em uma transformação linear que quando aplicada a uma equação diferencial ordinária linear a transforma em uma equação algébrica evitando com isso a necessidade do cálculo das derivadas do sinais de entrada e saída Além da análise detalhada desse método, onde demonstramos o efeito de cada parâmetro do método de Poisson sobre o desempenho desse, foi realizado também um estudo dos problemas decorrentes da discretização de sinais contínuos, como por exemplo o efeito aliasing decorrente da utilização de tempos de amostragem muito grandes e de problemas numéricos da identificação de modelos discretos utilizando dados com tempos de amostragem muito pequenos de forma a destacar as vantagens da identificação contínua sobre a identificação discreta Também foi estudado um método para compensar a presença de offsets nos sinais de entrada e saída, método esse inédito quando se trata do método dos momentos de Poisson. Esse trabalho também comprova a equivalência entre o método dos momentos de Poisson e uma metodologia apresentada por Rolf Johansson em um artigo de 1994. Na parte final desse trabalho são apresentados métodos para a compensação de erros de modelagem devido à presença de ruído e distúrbios não medidos nos dados utilizados na identificação. Esses métodos permitem que o método dos momentos de Poisson concorra com os métodos de identificação discretos normalmente empregados como por exemplo ARMAX e Box-Jenkins.

Aerodinâmica de automóveis baseado no método dos contornos virtuais

Este trabalho desenvolve um método numérico para a solução de escoamentos bidimensionais em torno de geometrias automobilísticas utilizando o método de diferenças finitas. O código computacional resolve as equações de Navier-Stokes e de Euler para uma distribuição adequada dos pontos discretos na malha. O método de integração empregado baseia-se no esquema explícito de Runge-Kutta de 3 estágios para as equações da quantidade de movimento e no de sub-relaxações sucessivas para a pressão na base Gauss-Seidel. Utilizou-se a técnica dos contornos virtuais em coordenadas cartesianas para resolver o escoamento sobre uma geometria simplificada, com a superfície coincidente com a malha computacional, e uma geometria automobilística mais complexa (BMW). Para a certificação da técnica empregada, optou-se pela utilização da teoria do escoamento potencial e pela comparação com dados experimentais encontrados na literatura e outros coletados em túnel de vento em escala reduzida. Houve dificuldade nesta comparação devido à falta de artigos relativos às simulações numéricas de escoamentos sobre automóveis e na aplicação da técnica dos contornos virtuais em geometrias complexas. Os resultados foram satisfatórios, com boas perspectivas para trabalhos futuros, contribuindo assim para o desenvolvimento da área.

Análise de estruturas compostas por aduelas pré-moldadas de concreto com protenção externa, sob estado plano de tensão, através do método dos elementos finitos

Aduelas protendidas externamente utilizadas em pontes contrastam com as clássicas construções monolíticas, consistindo de “pequenos” segmentos pré-moldados protendidos por cordoalhas externas. Há muitas vantagens na utilização deste tipo de composição, como rapidez e versatilidade na construção, alto controle de qualidade e economia. No presente trabalho, discute-se uma formulação, baseada no método dos elementos finitos, para simular o comportamento de estruturas compostas por aduelas externamente protendidas. Esta formulação permite tanto a análise de situações de serviço, onde as seções permanecem completamente comprimidas, como também de situações últimas, onde se considera a abertura das juntas e a transferência de esforços através das mesmas. O modelo computacional utiliza dois procedimentos distintos, baseados em um modelo elasto-viscoplástico. No primeiro, determina-se a resposta da estrutura ao longo do tempo, considerando-se, nesta etapa, que o material possui comportamento viscoelástico. No segundo, busca-se a resposta da estrutura para um carregamento instantâneo, considerando-se o material com um comportamento elastoplástico. As juntas entre as aduelas são representadas através de elementos de contato unidimensionais. Os nós de cada elemento de contato coincidem com os nós do elemento de concreto que a ele foi associado.Um modelo de camadas superpostas é utilizado para representar o comportamento viscoelástico do concreto e da armadura protendida Para representar o concreto fissurado, utiliza-se um modelo de fissuras distribuídas, que leva em conta a contribuição do concreto entre fissuras. A armadura passiva é introduzida na análise através de um modelo incorporado e a armadura de protensão não aderente é disposta em elementos fictícios, paralelos aos elementos reais da estrutura. Para verificar a adequação do modelo computacional, compararam-se resultados numéricos com valores experimentais disponíveis na literatura.

Alguns aspectos numéricos e teóricos das equações de Navier-Stokes na modelagem do escoamento em torno de um vórtice

Neste trabalho desenvolvemos uma metodologia numérica para a solução do escoamento em torno de um vórtice. Como a análise completa deste tipo de fluxo não é uma tarefa fácil, simplificações quanto ao escoamento e ao método numérico são necessárias. Também investigamos o comportamento das soluções das equações governantes (Navier-Stokes) quando o tempo tende ao infinito. Nesse sentido, dividimos este trabalho em duas partes: uma numérica e outra analítica. Com o intuito de resolver numericamente o problema, adotamos o método de diferenças finitas baseado na formulação incompressível das equações governantes. O método numérico para integrar essas equações é baseado no esquema de Runge- Kutta com três estágios. Os resultados numéricos são obtidos para cinco planos bidimensionais de um vórtice com números de Reynolds variando entre 1000 e 10000. Na parte analítica estudamos taxas de decaimento das soluções das equações de Navier-Stokes quando os dados iniciais são conhecidos. Também estimamos as taxas de decaimento para algumas derivadas das soluções na norma L2 e comparamos com as taxas correspondentes da solução da equação do calor.

Solução analítica do problema bidimensional transiente de dispersão de poluentes atmosféricos pelo método GITT dupla

Neste trabalho é apresentada uma solução analítica de um problema bidimensional e transiente de dispersão de poluentes atmosféricos. O modelamento utilizado é conhecido na literatura como modelo Kzz para dispersão de poluentes atmosféricos e é representado por uma equação difusivo-advectiva com coeficientes de difusão e advecção variáveis. São utilizados três diferentes coeficientes de difusão nas simulações, bem como as componentes horizontal e vertical do vento são tomadas como variáveis. A solução analítica é gerada através da aplicação da técnica GITT (Generalized Integral Transform Technique) dupla com problema transformado resolvido por Transformada de Laplace e diagonalização de matrizes. Filtros matemáticos são usados para homogenizar as condições de contorno viabilizando o uso da técnica citada. Além disso, o tipo de filtro matemático utilizado permite a sensível diminuição do custo computacional. Resultados numéricos são obtidos e comparados com dados experimentais e outras soluções da literatura.

Evolução de camada limite planetária para dispersão de poluentes pelo método da GILTT

O objetivo deste trabalho é obter os parâmetros turbulentos para o crescimento da camada limite planetária (CLP), durante a realizaçãoo do experimento Olad (Overland along wind dispersion experiment), conduzido na transição da noite para o dia. Nesta hora a CLP exibe uma altura, geralmente, pequena, disponibilizando pouco volume para a dispersão dos poluentes. Assim, concentrações superficiais elevadas podem ocorrer, atacando materiais, plantas e a saúde da população. Logo, conhecer os parâmetros do crescimneto é de fundamental importância para o correto modelamento da dispersão atmosférica ao amanhecer. A validação dos parâmetros é realizada a partir da solução da equação da difusão-advecção bidimensional, pelo método da GILTT (Generalized Integral Laplace Transform Technique). São empregados coeficientes de difusão turbulenta (problema de fechamento) dependentes da estabilidade atmosférica. As concentrações superficiais tridimensionais são obtidas através do espalhamento lateral da pluma com distribuição gaussiana. Apresentam-se os resultados numéricos e estatísticos, comparando os resultados obtidos com os dados experimentais. O modelo proposto mostrou-se aceitável em relação aos dados do experimento.

“Método de Gauss”: inapropriado até no nome

Tal como ressaltado em de Faro e Guerra (2014), tem sido frequente em nossos tribunais, sentenças judiciais determinando que, relativamente ao caso de amortizações de dívidas com prestações constantes, a popular Tabela Price seja substituída por um sistema que, fundamentado em uma particular aplicação do regime de juros simples, vem sendo cognominado de “Método de Gauss” (cf. Antonick e Assunção, 2006 e Nogueira, 2013). E isso, frize-se, mantendo-se o valor numérico da taxa de juros especificada no contrato de financiamento (usualmente, habitacional). A par de ser totalmente inadequado, como discutido em de Faro (2014c), associar o nome do grande matemático alemão Johann Carl Friedrich Gauss (1777-1855) ao procedimento em questão, sucede que ao mesmo, como a qualquer outro que seja baseado no regime de juros simples, associam-se incontornáveis inconsistências. Como já anteriormente, amplamente evidenciado em de Faro (2013b e 2014a). Tomando a Tabela Price como base de comparação, o propósito do presente trabalho é o de aprofundar a análise das deficiências do que tem sido denominado como “Método de Gauss”. Em particular, dado que as sentenças judiciais costumam não alterar os valores numéricos das taxas contratuais de juros, substituindo tão somente o regime de juros compostos, que está implícito na Tabela Price, pela peculiar variante do regime de juros simples que está subjacente ao que se chama de “Método de Gauss”, buscar-se-á considerar a questão do ponto de vista do financiador.

Um método de linearização local com passo adaptativo para solução numérica de equações diferenciais estocásticas com ruído aditivo

Neste trabalho apresentamos um novo método numérico com passo adaptativo baseado na abordagem de linearização local, para a integração de equações diferenciais estocásticas com ruído aditivo. Propomos, também, um esquema computacional que permite a implementação eficiente deste método, adaptando adequadamente o algorítimo de Padé com a estratégia “scaling-squaring” para o cálculo das exponenciais de matrizes envolvidas. Antes de introduzirmos a construção deste método, apresentaremos de forma breve o que são equações diferenciais estocásticas, a matemática que as fundamenta, a sua relevância para a modelagem dos mais diversos fenômenos, e a importância da utilização de métodos numéricos para avaliar tais equações. Também é feito um breve estudo sobre estabilidade numérica. Com isto, pretendemos introduzir as bases necessárias para a construção do novo método/esquema. Ao final, vários experimentos numéricos são realizados para mostrar, de forma prática, a eficácia do método proposto, e compará-lo com outros métodos usualmente utilizados.

Método Fuzzy para auxílio ao diagnóstico de câncer de mama em ambiente inteligente de telediagnóstico colaborativo para apoio à tomada de decisão

Breast cancer, despite being one of the leading causes of death among women worldwide is a disease that can be cured if diagnosed early. One of the main techniques used in the detection of breast cancer is the Fine Needle Aspirate FNA (aspiration puncture by thin needle) which, depending on the clinical case, requires the analysis of several medical specialists for the diagnosis development. However, such diagnosis and second opinions have been hampered by geographical dispersion of physicians and/or the difficulty in reconciling time to undertake work together. Within this reality, this PhD thesis uses computational intelligence in medical decision-making support for remote diagnosis. For that purpose, it presents a fuzzy method to assist the diagnosis of breast cancer, able to process and sort data extracted from breast tissue obtained by FNA. This method is integrated into a virtual environment for collaborative remote diagnosis, whose model was developed providing for the incorporation of prerequisite Modules for Pre Diagnosis to support medical decision. On the fuzzy Method Development, the process of knowledge acquisition was carried out by extraction and analysis of numerical data in gold standard data base and by interviews and discussions with medical experts. The method has been tested and validated with real cases and, according to the sensitivity and specificity achieved (correct diagnosis of tumors, malignant and benign respectively), the results obtained were satisfactory, considering the opinions of doctors and the quality standards for diagnosis of breast cancer and comparing them with other studies involving breast cancer diagnosis by FNA.

Análise de superfícies seletivas de freqüência através do método dos potenciais vetorias de Hertz

This work presents a theoretical and numerical analysis of Frequency Selective Surfaces (FSS) with elements as rectangular patch, thin dipole and crossed dipole mounted on uniaxial anisotropic dielectric substrate layers for orientations of the optical axis along x, y and z directions. The analysis of these structures is accomplished by combination of the Hertz vector potentials method and the Galerkin's technique, in the Fourier transform-domain, using entire¬domain basis functions. This study consists in the use of one more technique for analysis of FSS on anisotropic dielectric substrate. And presents as the main contribution the introduction of one more project parameter to determinate the transmission and reflection characteristics of periodic structures, from the use of anisotropic dielectric with orientations of the crystal optical axis along x, y and z directions. To validate this analysis, the numerical results of this work are compared to those obtained by other authors, for FSS structures on anisotropic and isotropic dielectric substrates. Also are compared experimental results and the numerical correspondent ones for the FSS isotropic case. The technique proposed in this work is accurate and efficient. ln a second moment, curves are presented for the transmission and reflection characteristics of the FSS structures using conducting patch elements mounted on uniaxial anisotropic dielectric substrate layers with optical axis oriented along x, y and z directions. From analysis of these curves, the performance of the considered FSS structures as function of the optical axis orientation is described

Modelo 2,5D de predição de propagação para ambientes interiores utilizando o método do traçado de raios

A 2.5D ray-tracing propagation model is proposed to predict radio loss in indoor environment. Specifically, we opted for the Shooting and Bouncing Rays (SBR) method, together with the Geometrieal Theory of Diffrartion (GTD). Besides the line-of-sight propagation (LOS), we consider that the radio waves may experience reflection, refraction, and diffraction (NLOS). In the Shooting and Bouncing Rays (SBR) method, the transmitter antenna launches a bundle of rays that may or may not reach the receiver. Considering the transmitting antenna as a point, the rays will start to launch from this position and can reach the receiver either directly or after reflections, refractions, diffractions, or even after any combination of the previous effects. To model the environment, a database is built to record geometrical characteristics and information on the constituent materials of the scenario. The database works independently of the simulation program, allowing robustness and flexibility to model other seenarios. Each propagation mechanism is treated separately. In line-of-sight propagation, the main contribution to the received signal comes from the direct ray, while reflected, refracted, and diffracted signal dominate when the line-of-sight is blocked. For this case, the transmitted signal reaches the receiver through more than one path, resulting in a multipath fading. The transmitting channel of a mobile system is simulated by moving either the transmitter or the receiver around the environment. The validity of the method is verified through simulations and measurements. The computed path losses are compared with the measured values at 1.8 GHz ftequency. The results were obtained for the main corridor and room classes adjacent to it. A reasonable agreement is observed. The numerical predictions are also compared with published data at 900 MHz and 2.44 GHz frequencies showing good convergence

Aplicação do método LTT às estruturas retangulares e triangulares em multicamadas e empilhadas em substratos PBG para comunicações móveis

Recently the planar antennas have been studied due to their characteristics as well as the advantages that they offers when compared with another types of antennas. In the mobile communications area, the need for this kind of antennas have became each time bigger due to the intense increase of the mobile communications that needs of antennas which operate in multifrequency and wide bandwidth. The microstrip antennas presents narrow bandwidth due the loss in the dielectric generated by radiation. Another limitation is the radiation pattern degradation due the generation of surface waves in the substrate. In this work some used techniques to minimize the disadvantages (previously mentioned) of the use of microstrip antennas are presented, those are: substrates with PBG material - Photonic Bandgap, multilayer antennas and with stacked patches. The developed analysis in this work used the TTL - Transverse Transmission Line method in the domain of Fourier transform, that uses a component of propagation in the y direction (transverse to the direction real of propagation z), treating the general equations of electric and magnetic field as functions of y and y . This work has as objective the application of the TTL method to microstrip structures with single and multilayers of rectangular and triangular patches, to obtaining the resonance frequency and radiation pattern of each structure. This method is applied for the treatment of the fields in stacked structures. The Homogenization theory will be applied to obtaining the effective permittivity for s and p polarizations of the substrate composed of PBG material. Numerical results for the triangular and rectangular antennas with single layer, multilayers resonators with triangular and rectangular patches are presented (in photonic and isotropic substrates). Conclusions and suggestions for continuity of this work are presented