Estudamos a otimização do método SOR clássico, para a resolução de um sistema linear Ax = b, com A não-singular, a partir dos resultados de Young [55, 57] e Varga [50, 51] para matrizes de blocos p-cíclicas consistentemente ordenadas. Num primeiro nível, a otimização refere-se à escolha do parâmetro de relaxação do SOR que produz a maior velocidade de convergência, e, num segundo nível, à escolha da p-ciclicidade que apresenta o melhor desempenho com os valores ótimos do parâmetro, e damos ênfase ao caso 2-cíclico. Além disso, descrevemos a otimização do parâmetro em três generalizações: a) num relaxamento das condições sobre o espectro da matriz de Jacobi associada a A; b) no método SOR para matrizes singulares; c) num novo método SOR, que substitui a decomposição A = D - L - U, onde D, L e U são a diagonal de A, a parte triangular inferior estrita de A e a parte triangular superior estrita de A, pela A = D - P - Q, onde P pertence a uma classe de matrizes constru ída a partir das matrizes-escada. Descrevemos também a aplicação do caso singular às cadeias de Markov, comentamos a computação paralela aplicada ao SOR, e apresentamos diversas simulações relativas à otimização desse método. / We study the optimization of the classic SOR method for solving a linear system Ax = b, where A is a nonsingular p-cyclic consistently ordered block matrix, based on the discoveries of Young [55, 57] and Varga [50, 51]. In a first levei, the optimization refers to the choice of the SOR relaxation parameter, which produces the greatest convergence speed and, in a second levei, to the p-cyclicity that presents the best performance with the optimal parameter values and emphasize the 2- cyclic case. Moreover we describe three SOR generalizations concerning optimization: a) by weakening the conditions on the spectrum of Jacobi matrix associated with A; b) by considering the SOR method for singular matrices; c) by approaching a new SOR, that replaces the splitting A = D - L - U, where O, L and U are the diagonal of A, the strict lower triangular part of A and the strict upper triangular part of A. respectively, by this one A = D - P - Q, where P is a stair matrix or a matrix even more general than a stair matrix. We also describe the application of the singular case to Markov chains, discuss parallel computing applied to SOR method, and present severa! simulations regarding the optimization of that method.

Estudamos a otimização do método SOR clássico, para a resolução de um sistema linear Ax = b, com A não-singular, a partir dos resultados de Young [55, 57] e Varga [50, 51] para matrizes de blocos p-cíclicas consistentemente ordenadas. Num primeiro nível, a otimização refere-se à escolha do parâmetro de relaxação do SOR que produz a maior velocidade de convergência, e, num segundo nível, à escolha da p-ciclicidade que apresenta o melhor desempenho com os valores ótimos do parâmetro, e damos ênfase ao caso 2-cíclico. Além disso, descrevemos a otimização do parâmetro em três generalizações: a) num relaxamento das condições sobre o espectro da matriz de Jacobi associada a A; b) no método SOR para matrizes singulares; c) num novo método SOR, que substitui a decomposição A = D - L - U, onde D, L e U são a diagonal de A, a parte triangular inferior estrita de A e a parte triangular superior estrita de A, pela A = D - P - Q, onde P pertence a uma classe de matrizes constru ída a partir das matrizes-escada. Descrevemos também a aplicação do caso singular às cadeias de Markov, comentamos a computação paralela aplicada ao SOR, e apresentamos diversas simulações relativas à otimização desse método. / We study the optimization of the classic SOR method for solving a linear system Ax = b, where A is a nonsingular p-cyclic consistently ordered block matrix, based on the discoveries of Young [55, 57] and Varga [50, 51]. In a first levei, the optimization refers to the choice of the SOR relaxation parameter, which produces the greatest convergence speed and, in a second levei, to the p-cyclicity that presents the best performance with the optimal parameter values and emphasize the 2- cyclic case. Moreover we describe three SOR generalizations concerning optimization: a) by weakening the conditions on the spectrum of Jacobi matrix associated with A; b) by considering the SOR method for singular matrices; c) by approaching a new SOR, that replaces the splitting A = D - L - U, where O, L and U are the diagonal of A, the strict lower triangular part of A and the strict upper triangular part of A. respectively, by this one A = D - P - Q, where P is a stair matrix or a matrix even more general than a stair matrix. We also describe the application of the singular case to Markov chains, discuss parallel computing applied to SOR method, and present severa! simulations regarding the optimization of that method.

Submitted by Cássia Santos (cassia.bcufg@gmail.com) on 2015-01-30T10:49:22Z No. of bitstreams: 2 license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Dissertação_Eduardo Koiti Hirota - 2014.pdf: 894506 bytes, checksum: 39a1f1c9a2e91954ecfdd1ef0513c5c0 (MD5) / Approved for entry into archive by Luciana Ferreira (lucgeral@gmail.com) on 2015-01-30T13:05:19Z (GMT) No. of bitstreams: 2 license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Dissertação_Eduardo Koiti Hirota - 2014.pdf: 894506 bytes, checksum: 39a1f1c9a2e91954ecfdd1ef0513c5c0 (MD5) / Made available in DSpace on 2015-01-30T13:05:19Z (GMT). No. of bitstreams: 2 license_rdf: 23148 bytes, checksum: 9da0b6dfac957114c6a7714714b86306 (MD5) Dissertação_Eduardo Koiti Hirota - 2014.pdf: 894506 bytes, checksum: 39a1f1c9a2e91954ecfdd1ef0513c5c0 (MD5) Previous issue date: 2014-10-09 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / Phenomenon that occur in the nature are essentially nonlinear and the dynamical systems theory aims to obtain a mathematical model that best represents the real physical systems, then nothing more coherent than the description or analysis of these natural phenomenon using models and techniques. In this dissertation, the technique of direct expansion for the development of two differential equations order to solve a nonlinear equation and the approximate determination of the roots of order algebraic equation higher or equal to two, was used. For this purpose, it was initially shown the development of a differential equation of motion subjected to a nonlinear damping, which is represented by the equation of Duffing – Van der Pol. Generally, it’s not easy to obtain an approximated analytical solution for this type equation, but this study was done with the purpouse of illustrating the technique used in the work, solving type solving a problem in which these techniques are routinely used to obtain a solution. Studied for application in basic education, it presents a way to obtain the approximate roots of equations of second and third degrees, using the technique of direct expansion for the sake of comparison. Since there are formulas for resolving this, It was proved that is possible to determine the roots of high-order equations by using the same technique. / Os fenômenos que ocorrem na natureza são essencialmente não lineares e a teoria de sistemas dinâmicos tem como objetivo obter um modelo matemático que represente melhor os sistemas físicos reais, então nada mais coerentes que a descrição ou análise desses fenômenos naturais usando modelos e técnicas não lineares. Nesta dissertação, foi utilizada a técnica da expansão direta para o desenvolvimento de equações diferenciais de ordem dois para resolução de uma equação não linear e na determinação aproximada de raízes de equações algébricas de ordem maior ou igual a dois. Com esse intuito, mostrou-se, inicialmente, o desenvolvimento de uma equação diferencial do movimento sujeito a um amortecimento não linear, que é representado pela equação de Duffing – Van der Pol. Geralmente, não é fácil obter uma solução analítica aproximada para esse tipo de equação, porém, este estudo é feito com a finalidade de ilustrar a técnica empregada no trabalho, resolvendo um tipo de problema no qual essas técnicas são corriqueiramente utilizadas para obter uma solução. Visando a aplicabilidade no ensino básico, apresenta-se uma forma de se obter as raízes aproximadas de equações do segundo e terceiro graus usando a técnica da expansão direta para efeito de comparação uma vez que existem fórmulas resolutivas para isso, provouse que é possível determinar as raízes de equações de ordem maior por meio da mesma técnica.