La simulation numérique de systèmes d’équations différentielles raides ordinaires ou algébriques est devenue partie intégrante dans le processus de conception des systèmes mécaniques à dynamiques complexes. L’objet de ce travail est de développer des méthodes numériques pour réduire les temps de calcul par le parallélisme en suivant deux axes : interne à l’intégrateur numérique, et au niveau de la décomposition de l’intervalle de temps. Nous montrons l’efficacité limitée au nombre d’étapes de la parallélisation à travers les méthodes de Runge-Kutta et DIMSIM. Nous développons alors une méthodologie pour appliquer le complément de Schur sur le système linéarisé intervenant dans les intégrateurs par l’introduction d’un masque de dépendance construit automatiquement lors de la mise en équations du modèle. Finalement, nous étendons le complément de Schur aux méthodes de type "Krylov Matrix Free". La décomposition en temps est d’abord vue par la résolution globale des pas de temps dont nous traitons la parallélisation du solveur non-linéaire (point fixe, Newton-Krylov et accélération de Steffensen). Nous introduisons les méthodes de tirs à deux niveaux, comme Parareal et Pita dont nous redéfinissons les finesses de grilles pour résoudre les problèmes raides pour lesquels leur efficacité parallèle est limitée. Les estimateurs de l’erreur globale, nous permettent de construire une extension parallèle de l’extrapolation de Richardson pour remplacer le premier niveau de calcul. Et nous proposons une parallélisation de la méthode de correction du résidu. / This PhD Thesis deals with the development of parallel numerical methods for solving Ordinary and Algebraic Differential Equations. ODE and DAE are commonly arising when modeling complex dynamical phenomena. We first show that the parallelization across the method is limited by the number of stages of the RK method or DIMSIM. We introduce the Schur complement into the linearised linear system of time integrators. An automatic framework is given to build a mask defining the relationships between the variables. Then the Schur complement is coupled with Jacobian Free Newton-Krylov methods. As time decomposition, global time steps resolutions can be solved by parallel nonlinear solvers (such as fixed point, Newton and Steffensen acceleration). Two steps time decomposition (Parareal, Pita,...) are developed with a new definition of their grids to solved stiff problems. Global error estimates, especially the Richardson extrapolation, are used to compute a good approximation for the second grid. Finally we propose a parallel deferred correction

Let F : Kn &rarr; Kp be a polynomial map, where K = R or C. Motivated by the characterization of the integral closure of ideals in the ring On by means of analytic inequalities proven by Lejeune-Teissier [46], we define the set Sp(F) of special polynomials with respect to F. The set Sp(F) can be considered as a counterpart, in the context of polynomial maps Kn &rarr; Kp, of the notion of integral closure of ideals in the ring of analytic function germs (~&lceil;+. In this work, we are mainly interested in the determination of the convex region S0(F) formed by the exponents of the special monomials with respect to F. Let us fix a convenient Newton polyhedron &lceil; + ~&sube; Rn. We obtain an approximation to S0</sub (F) when F is strongly adapted to ~&sube; +, which is a condition expressed in terms of the faces of ~&lceil;+ and the principal parts at infinity of F. The local version of this problem has been studied by Bivià-Ausina [4] and Saia [71]. Our result about the estimation of S0(F) allows us to give a lower estimate for the Lojasiewicz exponent at infinity of a given polynomial map with compact zero set. As a consequence of our study of ojasiewicz exponents at infinity we have also obtained a result about the uniformity of the ojasiewicz exponent in deformations of polynomial maps Kn &rarr; Kp. Consequently we derive a result about the invariance of the global index of real polynomial maps Rn &rarr; Rn. As particular cases of the condition of F being adapted to ~&lceil;+ there appears the class of Newton non-degenerate polynomial maps at infinity and pre-weighted homogeneous maps. The first class of maps constitute a natural extension for maps of the Newton non-degeneracy condition introduced by Kouchnirenko for polynomial functions. We characterize the Newton non-degeneracy at infinity condition of a given polynomial map F : Kn &rarr; Kp in terms of the set S0((F, 1)), where (F, 1) : Kn &rarr; Kp+1 is the polynomial map whose last component function equals 1. Motivated by analogous problems in local algebra we also derive some results concerning the multiplicity of F. / Seja F :Kn &rarr; Kp uma aplicação polinomial, onde K = C ou K = R. Motivados pela caracterização do fecho integral de ideais no anel On por meio de desigualdades analíticas provadas por Lejeune-Teissier [46], definimos o conjunto Sp(F) de polinomios especiais com respeito a F. O conjunto Sp(F) pode ser considerado como um homólogo, no contexto das aplicações polinomiais Kn &rarr; Kp, da noção de fecho integral de ideais no anel de germes de funções analíticas (Kn 0) &rarr; K. Neste trabalho, estamos interessados principalmente na determinação da região convexa S0 (F) formado pelos expoentes dos monômios especiais com respeito a F. Fixado um poliedro de Newton conveniente ~&lceil; + ~&sube; Rn, é obtida uma aproximação de S0(F), quando F é fortemente adaptada a &lceil; + o qual é uma condição expressada em termos das faces de ~&lceil; + e as partes principais no infinito de F. A versão local deste problema foi estudado por Bivià-Ausina [4] e Saia [71]. Nosso resultado sobre a estimativa de S0(F) nos permite dar uma estimativa inferior para o expoente Lojasiewicz no infinito de uma aplicação polinomial Kn &rarr; Kp, com conjunto F-1(0) compacto. Como uma consequência do estudo dos expoentes de Lojasiewicz no infinito também foi obtido um resultado sobre a uniformidade do expoente Lojasiewicz em deformações de aplicações polinomiais Kn &rarr; Kp e consequentemente, um resultado sobre a invariância do índice global de aplicações polinomiais reais Rn &rarr; Rn. Como casos particulares da condição de F ser adaptada a ~&lceil; + aparecem a classe de aplicações polinomiais Newton não degeneradas e as aplicações polinomiais pre-quase homogêneas. A primeira classe de aplicações constitui uma extensão natural da condição Newton não-degeneração introduzida por Kouchnirenko para funções polinomiais. Caracterizamos a condição Newton não-degeneração para uma determinada aplicação polinomial F : Kn &rarr; Kp em termos do conjunto S0((F, 1)), onde (F, 1) : Kn &rarr; Kp+1 é a aplicação polinomial cuja última função componente é igual a 1. Motivados por problemas análogos em álgebra local, também obtivemos alguns resultados sobre a multiplicidade de F.

Submitted by Cristiane Chim (cristiane.chim@ucpel.edu.br) on 2018-05-08T17:34:22Z No. of bitstreams: 1 Vagner Guidotti Furtado (1).pdf: 2942442 bytes, checksum: a43c18ecb28456284d4b6c622f11210d (MD5) / Made available in DSpace on 2018-05-08T17:34:22Z (GMT). No. of bitstreams: 1 Vagner Guidotti Furtado (1).pdf: 2942442 bytes, checksum: a43c18ecb28456284d4b6c622f11210d (MD5) Previous issue date: 2017-12-07 / The division operation in digital systems has its relevance because it is a necessary function in several applications, such as general purpose processors, digital signal processors and microcontrollers. The digital divider circuit is of great architectural complexity and may occupy a considerable area in the design of an integrated circuit, and as a consequence may have a great influence on the static and dynamic power dissipation of the circuit as a whole. In relation to the application of dividing circuits in circuits of the Digital Signal Processing (DSP) area, adaptive filters have a particular appeal, especially when using algorithms that perform a normalization in the input signals. In view of the above, this work focuses on the proposition of algorithms, techniques for reducing energy consumption and logical area, proposition and implementation of efficient dividing circuit architectures for use in adaptive filters. The Newton-Raphson and Goldschmidt iterative dividing circuits both operating at fixed-point were specifically addressed. The results of the synthesis of the implemented architectures of the divisors with the proposed algorithms and techniques showed considerable reduction of power and logical area of the circuits. In particular, the dividing circuits were applied in adaptive filter architectures based on the NLMS (Normalized least Mean Square) algorithm, seeking to add to these filters, characteristics of good convergence speed, combined with the improvement in energy efficiency. The adaptive filters implemented are used in the case study of harmonic cancellation on electrocardiogram signals / A operação de divisão em sistemas digitais tem sua relevância por se tratar de uma função necessária em diversas aplicações, tais como processadores de propósito geral, processadores digitais de sinais e microcontroladores. O circuito divisor digital é de grande complexidade arquitetural, podendo ocupar uma área considerável no projeto de um circuito integrado, e por consequência pode ter uma grande influência na dissipação de potência estática e dinâmica do circuito como um todo. Em relação à aplicação de circuitos divisores em circuitos da área DSP (Digital Signal Processing), os filtros adaptativos têm um particular apelo, principalmente quando são utilizados algoritmos que realizam uma normalização nos sinais de entrada. Diante do exposto, este trabalho foca na proposição de algoritmos, técnicas de redução de consumo de energia e área lógica, proposição e implementação de arquiteturas de circuitos divisores eficientes para utilização em filtros adaptativos. Foram abordados em específico os circuitos divisores iterativos Newton-Raphson e Goldschmidt ambos operando em ponto-fixo. Os resultados da síntese das arquiteturas implementadas dos divisores com os algoritmos e técnicas propostas mostraram considerável redução de potência e área lógica dos circuitos. Em particular, os circuitos divisores foram aplicados em arquiteturas de filtros adaptativos baseadas no algoritmo NLMS (Normalized least Mean Square), buscando agregar a esses filtros, características de boa velocidade de convergência, aliada à melhoria na eficiência energética. Os filtros adaptativos implementados são utilizados no estudo de caso de cancelamento de harmônicas em sinais de eletrocardiograma (ECG)

No description available.

Cette thèse est consacrée à la recherche des zéros d'un opérateur maximal monotone structuré, à l'aide de systèmes dynamiques dissipatifs continus et discrets. Les solutions sont obtenues comme limites des trajectoires lorsque le temps t tend vers l'infini. On s'intéressera principalement aux dynamiques obtenues par régularisation de type Levenberg-Marquardt de la méthode de Newton. On décrira aussi les approches basées sur des dynamiques voisines.Dans un cadre Hilbertien, on s'intéresse à la recherche des zéros de l'opérateur maximal monotone structuré M = A + B, où A est un opérateur maximal monotone général et B est un opérateur monotone Lipschitzien. Nous introduisons des dynamiques continues et discrètes de type Newton régularisé faisant intervenir d'une façon séparée les résolvantes de l'opérateur A (implicites), et des évaluations de B (explicites). A l'aide de la représentation de Minty de l'opérateur A comme une variété Lipschitzienne, nous reformulons ces dynamiques sous une forme relevant du théorème de Cauchy-Lipschitz. Nous nous intéressons au cas particulier où A est le sous différentiel d'une fonction convexe, semi-continue inférieurement, et propre, et B est le gradient d'une fonction convexe, différentiable. Nous étudions le comportement asymptotique des trajectoires. Lorsque le terme de régularisation ne tend pas trop vite vers zéro, et en s'appuyant sur une analyse asymptotique de type Lyapunov, nous montrons la convergence des trajectoires. Par ailleurs, nous montrons la dépendance Lipschitzienne des trajectoires par rapport au terme de régularisation.Puis nous élargissons notre étude en considérant différentes classes de systèmes dynamiques visant à résoudre les inclusions monotones gouvernées par un opérateur maximal monotone structuré M = $partialPhi$+ B, où $partialPhi$ désigne le sous différentiel d'une fonction convexe, semicontinue inférieurement, et propre, et B est un opérateur monotone cocoercif. En s'appuyant sur une analyse asymptotique de type Lyapunov, nous étudions le comportement asymptotique des trajectoires de ces systèmes. La discrétisation temporelle de ces dynamiques fournit desalgorithmes forward-backward (certains nouveaux ).Finalement, nous nous intéressons à l'étude du comportement asymptotique des trajectoires de systèmes dynamiques de type Newton régularisé, dans lesquels on introduit un terme supplémentaire de viscosité évanescente de type Tikhonov. On obtient ainsi la sélection asymptotique d'une solution de norme minimale. / This thesis is devoted to finding zeroes of structured maximal monotone operators, by using discrete and continuous dissipative dynamical systems. The solutions are obtained as the limits of trajectories when the time t tends towards infinity.We pay special attention to the dynamics that are obtained by Levenberg-Marquardt regularization of Newton's method. We also revisit the approaches based on some related dynamical systems.In a Hilbert framework, we are interested in finding zeroes of a structured maximal monotone operator M = A + B, where A is a general maximal monotone operator, and B is monotone and locally Lipschitz continuous. We introduce discrete and continuous dynamical systems which are linked to Newton's method. They involve separately B and the resolvents of A, and are designed to splitting methods. Based on the Minty representation of A as a Lipschitz manifold, we show that these dynamics can be formulated as differential systems, which are relevant to the Cauchy-Lipschitz theorem. We focus on the particular case where A is the subdifferential of a convex lower semicontinuous proper function, and B is the gradient of a convex, continuously differentiable function. We study the asymptotic behavior of trajectories. When the regularization parameter does not tend to zero too rapidly, and by using Lyapunov asymptotic analysis, we show the convergence of trajectories. Besides, we show the Lipschitz continuous dependence of the solution with respect to the regularization term.Then we extend our study by considering various classes of dynamical systems which aim at solving inclusions governed by structured monotone operators M = $partialPhi$+ B, where $partialPhi$ is the subdifferential of a convex lower semicontinuous function, and B is a monotone cocoercive operator. By a Lyapunov analysis, we show the convergence properties of the orbits of these systems. The time discretization of these dynamics gives various forward-backward splittingmethods (some new).Finally, we focus on the study of the asymptotic behavior of trajectories of the regularized Newton dynamics, in which we introduce an additional vanishing Tikhonov-like viscosity term.We thus obtain the asymptotic selection of the solution of minimal norm.

This thesis presents a study of impedance optimisation of active microstrip patch antennas to multiple frequency points. A single layered aperture coupled microstrip patch antenna has been optimised to match the source reflection coefficient of a transistor in designing an active antenna. The active aperture coupled microstrip patch antenna was optimised to satisfy Global Positioning System (GPS) frequency specifications. A rudimentary aperture coupled microstrip patch antenna consists of a rectangular antenna element etched on the top surface of two dielectric substrates. The substrates are separated by a ground plane and a microstrip feed is etched on the bottom surface. A rectangular aperture in the ground plane provides coupling between the feed and the antenna element. This type of antenna, which conveniently isolates any circuit at the feed from the antenna element, is suitable for integrated circuit design and is simple to fabricate. An active antenna design directly couples an antenna to an active device, therefore saving real estate and power. This thesis focuses on designing an aperture coupled patch antenna directly coupled to a low noise amplifier as part of the front end of a GPS receiver. In this work an in-house software package, dubbed ACP by its creator Dr Rod Waterhouse, for calculating aperture coupled microstrip patch antenna performance parameters was linked to HP-EEsof, a microwave computer aided design and simulation package by Hewlett-Packard. An ANSI C module in HP-EEsof was written to bind the two packages. This process affords the client the benefit of powerful analysis tools offered in HP-EEsof and the fast analysis of ACP for seamless system design. Moreover, the optimisation algorithms in HP-EEsof were employed to investigate which algorithms are best suited for optimising patch antennas. The active antenna design presented in this study evades an input matching network, which is accomplished by designing the antenna to represent the desired source termination of a transistor. It has been demonstrated that a dual-band microstrip patch antenna can be successfully designed to match the source reflection coefficient, avoiding the need to insert a matching network. Maximum power transfer in electrical circuits is accomplished by matching the impedance between entities, which is generally acheived with the use of a matching network. Passive matching networks employed in amplifier design generally consist of discrete components up to the low GHz frequency range or distributed elements at greater frequencies. The source termination for a low noise amplifier will greatly influence its noise, gain and linearity which is controlled by designing a suitable input matching network. Ten diverse search methods offered in HP-EEsof were used to optimise an active aperture coupled microstrip patch antenna. This study has shown that the algorithms based on the randomised search techniques and the Genetic algorithm provide the most robust performance. The optimisation results were used to design an active dual-band antenna.

Dans ce mémoire, on se propose de décrire certains résultats de l'auteur sur les propriétés $p$-adiques des fonctions $L$ associées à des caractères sur les corps finis, à la suite des travaux de Dwork, Robba, Adolphson et Sperber, Wan, entre autres. On parlera aussi de sommes de caractères (et de leurs fonctions $L$) définies sur certains anneaux locaux.

Esta dissertação aplica a regularização por entropia máxima no problema inverso de apreçamento de opções, sugerido pelo trabalho de Neri e Schneider em 2012. Eles observaram que a densidade de probabilidade que resolve este problema, no caso de dados provenientes de opções de compra e opções digitais, pode ser descrito como exponenciais nos diferentes intervalos da semireta positiva. Estes intervalos são limitados pelos preços de exercício. O critério de entropia máxima é uma ferramenta poderosa para regularizar este problema mal posto. A família de exponencial do conjunto solução, é calculado usando o algoritmo de Newton-Raphson, com limites específicos para as opções digitais. Estes limites são resultados do princípio de ausência de arbitragem. A metodologia foi usada em dados do índice de ação da Bolsa de Valores de São Paulo com seus preços de opções de compra em diferentes preços de exercício. A análise paramétrica da entropia em função do preços de opções digitais sínteticas (construídas a partir de limites respeitando a ausência de arbitragem) mostraram valores onde as digitais maximizaram a entropia. O exemplo de extração de dados do IBOVESPA de 24 de janeiro de 2013, mostrou um desvio do princípio de ausência de arbitragem para as opções de compra in the money. Este princípio é uma condição necessária para aplicar a regularização por entropia máxima a fim de obter a densidade e os preços. Nossos resultados mostraram que, uma vez preenchida a condição de convexidade na ausência de arbitragem, é possível ter uma forma de smile na curva de volatilidade, com preços calculados a partir da densidade exponencial do modelo. Isto coloca o modelo consistente com os dados do mercado. Do ponto de vista computacional, esta dissertação permitiu de implementar, um modelo de apreçamento que utiliza o princípio de entropia máxima. Três algoritmos clássicos foram usados: primeiramente a bisseção padrão, e depois uma combinação de metodo de bisseção com Newton-Raphson para achar a volatilidade implícita proveniente dos dados de mercado. Depois, o metodo de Newton-Raphson unidimensional para o cálculo dos coeficientes das densidades exponenciais: este é objetivo do estudo. Enfim, o algoritmo de Simpson foi usado para o calculo integral das distribuições cumulativas bem como os preços do modelo obtido através da esperança matemática. / This study aims at applying Maximum Entropy Regularization to the Inverse Problem of Option Pricing suggested by Neri and Schneider in 2012. They pointed out that the probability density that solves such problem in the case of calls and digital options could be written as piecewise exponentials on the positive real axis. The limits of these segments are the different strike prices. The entropy criteria is a powerful tool to regularize this ill-posed problem. The Exponential Family solution set is calculated using a Newton-Raphson algorithm, with specific bounds for the binary options. These bounds obey the no-arbitrage principle. We applied the method to data from the Brazilian stock index BOVESPA and its call prices for different strikes. The parametric entropy analysis for "synthetic" digital prices (constructed from the no-arbitrage bounds) showed values where the digital prices maximizes the entropy. The example of data extracted on the IBOVESPA of January 24th 2013, showed slippage from the no-arbitrage principle when the option was in the money: such principle is a necessary condition to apply the maximum entropy regularization to get the density and modeled prices. When the condition is fulfilled, our results showed that it is possible to have a smile-like volatility curve with prices calculated from the exponential density that fit well the market data. In a computational modelling perspective, this thesis enabled the implementation of a pricing method using the maximum entropy principle. Three well known algorithms were used in that extent. The bisection alone, then a combined bisection with Newton-Raphson to recover the implied volatility from market data. Thereafter, the one dimensional Newton-Raphson to calculate the coefficients of the exponential densities: purpose of the study. Finally the Simpson method was used to calculate integrals of the cumulative distributions and the modeled prices implied by the expectation.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A engenharia geotécnica é uma das grandes áreas da engenharia civil que estuda a interação entre as construções realizadas pelo homem ou de fenômenos naturais com o ambiente geológico, que na grande maioria das vezes trata-se de solos parcialmente saturados. Neste sentido, o desempenho de obras como estabilização, contenção de barragens, muros de contenção, fundações e estradas estão condicionados a uma correta predição do fluxo de água no interior dos solos. Porém, como a área das regiões a serem estudas com relação à predição do fluxo de água são comumente da ordem de quilômetros quadrados, as soluções dos modelos matemáticos exigem malhas computacionais de grandes proporções, ocasionando sérias limitações associadas aos requisitos de memória computacional e tempo de processamento. A fim de contornar estas limitações, métodos numéricos eficientes devem ser empregados na solução do problema em análise. Portanto, métodos iterativos para solução de sistemas não lineares e lineares esparsos de grande porte devem ser utilizados neste tipo de aplicação. Em suma, visto a relevância do tema, esta pesquisa aproximou uma solução para a equação diferencial parcial de Richards pelo método dos volumes finitos em duas dimensões, empregando o método de Picard e Newton com maior eficiência computacional. Para tanto, foram utilizadas técnicas iterativas de resolução de sistemas lineares baseados no espaço de Krylov com matrizes pré-condicionadoras com a biblioteca numérica Portable, Extensible Toolkit for Scientific Computation (PETSc). Os resultados indicam que quando se resolve a equação de Richards considerando-se o método de PICARD-KRYLOV, não importando o modelo de avaliação do solo, a melhor combinação para resolução dos sistemas lineares é o método dos gradientes biconjugados estabilizado mais o pré-condicionador SOR. Por outro lado, quando se utiliza as equações de van Genuchten deve ser optar pela combinação do método dos gradientes conjugados em conjunto com pré-condicionador SOR. Quando se adota o método de NEWTON-KRYLOV, o método gradientes biconjugados estabilizado é o mais eficiente na resolução do sistema linear do passo de Newton, com relação ao pré-condicionador deve-se dar preferência ao bloco Jacobi. Por fim, há evidências que apontam que o método PICARD-KRYLOV pode ser mais vantajoso que o método de NEWTON-KRYLOV, quando empregados na resolução da equação diferencial parcial de Richards. / Geotechnical Engineering is the area of Civil Engineering that studies the interaction between constructions carried out by man or natural phenomena with geological environment, which most of times is partially saturated soil. In this sense, work developing as stabilization, dam containing, retaining walls, foundations and highways are conditioned to a right prediction of water flow into the soil. However, considering the water flow, the studied region areas are commonly on the order of square kilometers, mathematical models solutions require computational meshes of large proportions, causing serious limitations linked to computational memory requirements and processing time. In order to overcome these limitations, efficient numerical methods must be used in the solution of the considered problem. Hence iterative methods for solving nonlinear and large sparse linear systems must be used in this type of application. In short, this study approached a solution to the Richard partial differential equation by the two dimensions finite volume method, bringing Picard and Newton method with greater efficiency. Linear system resolution iterative techniques based on Krylov space with pre-conditioners matrix were used. Portable Extensible Toolkit for Scientific Computation (PETSc) numerical library was a tool used during the task. The results indicate when a Richards equation is solved considering thr PICARD-KRYLOV method, no matter the soil evaluation model, the best combination for solving linear systems is the stabilized double gradient method and the SOR preconditioning. On the other hand, when the van Genuchten equations are used the gradients methods with the SOR preconditioning must be chosen. Adopting the NEWTON-KRYLOV method, the stabilized double gradient method is more efficient in soling Newton linear system, in relation to the preconditioning it must be giving preference to the Jacob block. Finally, there are strong indications that the PICARDKRYLOV method can be more effective than the NEWTON-KRYLOV one, when used for solving Richards partial differential equation.

Submitted by Maike Costa (maiksebas@gmail.com) on 2016-03-30T11:12:00Z No. of bitstreams: 1 arquivo total.pdf: 1989591 bytes, checksum: c1b3f2740144367fd7ef458d0603ba20 (MD5) / Made available in DSpace on 2016-03-30T11:12:00Z (GMT). No. of bitstreams: 1 arquivo total.pdf: 1989591 bytes, checksum: c1b3f2740144367fd7ef458d0603ba20 (MD5) Previous issue date: 2015-02-24 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / In several scienti c character problems, it is common to come across us with the need to obtain an approximate solution to nd roots of functions. At this point, this paper aims to conduct a study about some methods used to obtain an approximate solution of the functions of roots. The survey was made by means of a literature review, focusing on Numerical Methods Bisection, False Position, Fixed Point, Newton-Raphson and Secant. In order to illustrate the operation and application of these methods, numerical test problems taken from the literature were performed by implementing these. For each test performed were analyzed parameters that in uence each method and the convergence situation for the approximate solution of the analyzed problems. Although these methods do not always make available exact roots, they can be calculated with the precision that the problem needs. At this point, it is evident the importance of studying methods for nding such equations roots. Thus, the work is justi ed on the need to discuss the problems facing the nding roots of polynomial functions in the literature. In addition, this paper describes a comparison between the methods studied by applying mathematical problems. All this research material becomes adept and e ective for students and professionals from all areas that make use of them, or perhaps wish to extract it for enrichment of several sources of study. / Em diversos problemas de caráter cientí co, é comum depararmo-nos com a necessidade de obter uma solução aproximada para encontrar raízes de funções. Nesse ponto, este trabalho objetiva realizar um estudo acerca de alguns métodos utilizados para a obtenção de uma solução aproximada das raízes de funções. A pesquisa realizada deu-se por meio de uma revisão bibliográ ca, enfocando os Métodos Numéricos da Bisseção, Falsa Posição, Ponto Fixo, Newton-Raphson, Secante e Muller. Com o intuito de ilustrar o funcionamento e aplicação desses métodos, foram realizados testes numéricos de problemas extraídos da literatura por meio da implementação destes. Para cada teste realizado foram analisados os parâmetros que in uenciam cada método e a situação de convergência para a solução aproximada dos problemas analisados. Embora esses métodos, nem sempre, disponibilizem raízes exatas, estas poderão ser calculadas com a precisão que o problema necessite. Nesse ponto, ca evidente a importância de estudar métodos para encontrar tais raízes de equações. Diante disso, o trabalho se justi ca na necessidade de se discutir os problemas voltados a encontrar raízes de funções polinomiais, existentes na literatura. Além disso, o presente trabalho descreve um comparativo entre os métodos estudados mediante aplicação de problemas matemáticos. Todo esse material de pesquisa torna-se hábil e e caz para os estudantes e pro ssionais de todas as áreas que dele faça uso, ou, porventura, pretendam extraí-lo para enriquecimento de fontes diversas de estudo.