Resumo: Este projeto de pesquisa tem como objetivo a concepção, implementação, solução e teste de um modelo de Pré-Despacho de Geração (PDG) para o ambiente de mercados de energia, que supra as deficiências dos modelos de PDG adotados no Brasil. Assim, a abordagem proposta deve introduzir novos aspectos de modelagem, tais como: i) a introdução de aspectos associados aos mercados de energia internamente ao modelo de PDG; ii) a representação das inter-relações entre os mercados pool e bilateral em um único modelo de PDG; iii) a discretização do problema em base horária, possibilitando, de fato, a implementação de um mercado de curtíssimo prazo; iv) a avaliação da função de custo de oportunidade como base para a inserção de objetivos associados à otimização da produção de energia hidráulica no mercado pool / Abstract: This research aims at the conception, implementation, solution and testing of the proposed Short term Generation Scheduling Model (PDG), specific for the energy market environment. This model focuses on the improvement in the dispatch model used by the Brazilian energy sector. The proposed approach introduces brand new modeling aspects, such as: i) the introduction of modeling aspects associated with energy markets into the PDG model; ii) the representation of the interrelation between pool and bilateral markets within a single optimization problem; iii) the discretization of the problem is introduced in an hourly basis aiming at the implementation of an effective short time energy market; iv) the evaluation of the opportunity costs function as a basis for insertion of objectives associated with optimization of hydraulic energy production in pool market / Orientador: Leonardo Nepomuceno / Coorientador: Paulo Sérgio da Silva / Banca: Takaaki Ohishi / Banca: Edmea Cassia Baptista / Mestre

A detecção de dano estrutural é um assunto muito estudado em virtude do grau de responsabilidade envolvido na manutenção de certas estruturas por diversas décadas, tais como pontes, torres de transmissão de energia, ou mesmo em falhas mecânicas de componentes automotivos, navais, de aviação, etc. A identificação do dano em estruturas por meio de ensaios dinâmicos não destrutivos tem se tornado uma técnica atrativa, devido principalmente a redução de incertezas e queda nos preços dos equipamentos utilizados, tais como acelerômetros, placas de aquisição de sinais, computadores, etc. Simulated Annealing (Recozimento Simulado) é uma técnica baseada em algoritmos heurísticos utilizada para otimização (minimização ou maximização) de problemas de larga escala, sem necessidade de cálculos de gradientes, que contenham diversos ótimos locais, onde é necessário encontrar um ótimo global. No caso de detecção de danos, esta otimização pode ser usada para sintonizar o comportamento dinâmico de um modelo paramétrico com resultados experimentais modais, como por exemplo: freqüências naturais, amortecimento modal, modos, etc. Neste trabalho é avaliado o desempenho do emprego deste método na localização e quantificação de danos em estruturas simples. Alguns dos resultados foram comparados com uma outra técnica heurística de Algoritmos Genéticos. A avaliação das incertezas envolvidas na determinação dos parâmetros experimentais foi realizada de modo a se ter idéia da ordem de exatidão do sistema assim como limitações das previsões feitas por esta metodologia. Esta abordagem é de extrema importância, pois o sistema pode indicar falsos danos em virtude de falta de exatidão nas medições das freqüências naturais da estrutura estudada. Uma abordagem teórica nos casos da viga bi apoiada e engastada livre foi feita através da equação diferencial da viga. Esta abordagem demonstrou regiões onde certas freqüências não são alteradas independentemente da condição de dano. Esta é uma das maiores contribuições deste trabalho. O mesmo estudo foi feito de forma experimental, onde se verificou a confirmação dos lugares sugeridos pelos dois estudos anteriores, ou seja, determinadas freqüências não se alteram quando determinadas regiões da estrutura são danificadas. Concluiu-se que o método de simulated annealing tem eficiência equivalente aos algoritmos genéticos, pois em termos da eficiência em detectar corretamente danos, ele apresentou desempenho semelhante, o que foi comprovado pelas simulações numéricas e ensaios experimentais. As dificuldades de identificação de danos, como danos próximos as condições de contorno e danos menores de 10%, foram sentidos de forma parecida, sendo que para estes casos, o algoritmo não convergiu completamente. Os tempos de processamento em ambos os métodos de otimização foram similares, para todos os casos estudados. / The detection of structural damage has received an increasing interest on both academic research and industry application for several decades, including structures such as bridges, transmission lines towers, or even in automotive mechanical fault detection and naval and aviation components. The damage identification in structures by means of dynamic non-destructive tests became an attractive technique, since there was and increase in the measurement accuracy and cost decrease of devices used to perform this task, like accelerometers, data acquisition boards, computers, etc. Simulated Annealing is a technique that belongs to a class of the so called Heuristic Algorithms used in optimizations (maximizations or minimizations) of large scale problems. This technique does not use gradient evaluations and can be used with functions with several local optima, where it is desired to reach the global optimum. In damage detection, this technique can be used to tune the dynamic behavior of a parametric model with experimental modal results, such as: natural frequencies, modal damping, mode shapes, etc. In this work, the performance of this method in finding and quantificating damage in simple structures is evaluated. Some of the results were compared with another Heuristic technique: the Genetic Algorithm. The importance of uncertainties in the involved experimental parameters was investigated in order to estimate the system accuracy and limitations of the methodology. This approach is important, since the methodology could indicate false damages by virtue of lack of accuracy in the measurements. A theoretical approach for simple beams was accomplished by the differential equation of motion and a simple model for the crack. This study shows points where damage will not affect some frequencies, regardless of their extent. This is one of the main contributions of this work. The same study was carried experimentally, where there was a perfect agreement with the predicted sites for null interference suggested by the two previous studies, i.e., certain frequencies did not change when those sites where damaged. Finally it is concluded that the Simulated Annealing Technique has equivalent efficiency as the genetic algorithm, since the performance in detecting damage was similar, checked by the numerical simulations and experimental tests. The difficulties on damage identification lied in damage close to support conditions and damage about 10%, and for these cases, the algorithms did not converge completely. The processing times on both techniques used were similar, for all the studied cases.

Orientador: Jose Mario Martinez / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-07-23T01:06:27Z (GMT). No. of bitstreams: 1 Bielschowsky_RobertoHugo_D.pdf: 5956930 bytes, checksum: 49d1d2e75a72890d87aecbedd953965e (MD5) Previous issue date: 1997 / Resumo: Abordamos, nesta tese, o problema de obter pontos de mínimo local de funções diferenciáveis, definidas no lRn e sujeitas a restrições. Nosso ponto de partida reside numa aposta em algoritmos que têm muito em comum com algoritmos de pontos factíveis, tais como, por exemplo, o GRG e o Gradiente Projetado, porém relaxando de forma dinâmica as restrições de igualdade h(x) = 0. Ou seja, relaxaremos a condição h(x(k)) = 0, característica dos iterandos gerados em métodos de pontos factíveis, para uma na forma ||h(x(k))|| = O(||gp(x(k))||).gp(x) representa a projeção ortogonal do gradiente V¿(x), no espaço tangente às restrições N(h'(x)). No capítulo 1 situamos nossa abordagem. No capítulo 2 formulamos um algoritmo desenvolvendo-a para restrições de igualdade apenas, e que denominaremos de CDR (Controle Dinâmico das Restrições). Pensando em problemas de grande porte não estruturados formulamos uma versão adequada a tratar de forma inexata todos os subproblemas lineares envolvidos. Vale dizer, sem fatorações de matrizes. Ainda no segundo capítulo desenvolvemos uma teoria de convergência global para o método, e no terceiro uma teoria de convergência local. No quarto capítulo apresentamos os resultados de alguns testes preliminares com o algoritmo, realizados em colaboração com Francisco M. Gomes. No quinto capítulo e no apêndice tratamos de possíveis extensões de CDR, visando incluir também restrições de desigualdade. / Abstract: Not informed. / Doutorado / Doutor em Matemática Aplicada

Orientador: Jose Mario Martinez / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-07-26T03:15:42Z (GMT). No. of bitstreams: 1 Pilotta_ElvioAngel_D.pdf: 3939443 bytes, checksum: 1c19786cda4ad43e405851398d5d1a07 (MD5) Previous issue date: 2000 / Resumo: Nesta tese apresentamos métodos numéricos para problemas de minimização com restrições. O Capítulo 1 está baseado no artigo "Validation of an Augmented Lagrangian algorithm with a Gauss-Newton Hessian approximation using a set of Hard-Spheres Problems", de Krejié, Martínez, Mello e Pilotta. O Capítulo 2 está baseado no artigo "Inexact-Restoration algorithm for constrained optimization" , de Martínez e Pilotta, onde é considerado um novo método de tipo restauração inexata para um problema de minimização com restrições gerais. O Capítulo 3 estábaseado no artigo "Spectral Gradient method for linearly constrained optimization" , de Martínez, Pilotta e Raydan, onde é considerado um novo método para um problema de minimização com restrições lineares e canalizações usando gradiente espectral precondicionado e penalização exponencial. O Capítulo 4 está baseado no artigo "A limited-memory multipoint secant method for bound constrained optimization", de Burdakov, Martínez e Pilotta, onde é considerado um novo método para um problema de minimização com canalizações usando uma estratégia de restrições ativas e um método secante simétrico multipoint com memória limitada para resolver um subproblema quadrático em cada face / Abstract: We present numerical methods for constrained minimization problems. Chapter 1 is based on the paper "Validation of an Augmented Lagrangian algorithm with a Gauss-Newton Hessian approximation using a set of Hard-Spheres Problems", by Krejié, Martínez, Mello and Pilotta. Chapter 2 is based on the paper "InexactRestoration algorithm for constrained optimization" , by Martínez and Pilotta, where we introduce an inexact-restoration method for solving a general constrained minimization problem. Chapter 3 is based on the paper "Spectral Gradient method for linearly constrained optimization", by Martínez, Pilotta, and Raydan, where we introduce a new method for this problem which uses exponential penalization. Chapter 4 is based on the paper " A limited-memory multipoint secant method for bound constrained optimization", by Burdakov, Martínez and Pilotta, where we introduce a new method for bound constrained optimization that uses active set methods for solving a quadratic subproblem in each face / Doutorado / Doutor em Matemática Aplicada

Orientadores: Ivan de Araujo Simões Filho, Jesse Carvalho Costa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica / Made available in DSpace on 2018-07-26T15:19:33Z (GMT). No. of bitstreams: 1 Vanzeler_FranciscoJocleanAlves_M.pdf: 5979791 bytes, checksum: bdcf42ac30bd9b5de4f83039d6bc75a8 (MD5) Previous issue date: 1999 / Resumo: Este trabalho calcula os parâmetros elásticos e a densidade de um meio fraturado a partir da inversão dos coeficientes de reflexão de ondas planas qP entre dois meios anisotrópicos com pelo menos um plano de simetria paralelo a interface. O meio fraturado é modelado por um meio anisotrópico equivalente através da formulação de Hudson, neste caso um reservatório isotrópico com fraturas verticais é equivalente a um meio transversalmente isotrópico com eixo de simetria horizontal (Tlli). Sobre o reservatório, três tipos de rocha são considerados para a modelagem direta: 1) Folhelho lsotrópico; 2) Folhelho tranversalmente isotrópico com eixo de simetria vertical (TIV); 3) Folhelho com Fraturas verticais ou meio ortorrômbico. Os coeficientes de reflexão usados neste trabalho são calculados de forma exata. Conjuntos de dados sintéticos para múltiplos azimutes e ângulos de incidência, contaminados por ruído gaussiano, foram gerados para cada modelo considerado. Dois algoritmos de otimização global foram utilizados para inversão: Os algoritmos Genético (AG) e Evolução Diferencial (ED), que estimam os parâmetros do modelo através da minimização da norma 12do resíduo entre o coeficientes de reflexão "observado" e o calculado. Foi fornecido os valores dos parâmetros elásticos do meio superior exceto no caso onde a velocidade da onda cisalhante (Vs) do meio superior foi invertida. Cada parâmetro varia dentro dos limites observados durante a propagação de ondas sísmicas em rochas sedimentares. A inversão foi satisfatória nos três modelos estudados. Três azimutes (0°, 45° e 9Oj foram considerados na geometria de aquisição. Para o modelo 1,a incidência mínima para inversão foi de trinta e quatro graus (AG e ED). Para os modelos II e III, a incidência mínima para inversão foi de cinqüenta graus (AG). O ED inverteu o caso que inclui Vs do meio superior do modelo I com incidência mínima de quarenta e quatro graus. Esta abordagem mostra o potencial da análise AVO azimutal de ondas qP para caracterizar a anisotropia associada a meios fraturados / Abstract: In this work. we extract the elastic stifIDess and mass density from an multiazimutal qP-wave reflectioncoeficients at an interface separating two anisotropic media with monoclinic simmetry with at least one of its planes of simmetry parallel to the interface. This objective was reach by forward and inverse modeling. We calculate the qPwave reflection for three models (1. TI.Ill) of anisotropic equivalent medium: 1) isotropic mediumabove a TIH medium; 2) TIV medium above a TIH medium; and 3) orthorhombic medium above a TIH medium. The TIH medium is equivalent an isotropic fractured medium withequivalent elastic stiffness and mass density calculeted by the Hudson formulation. The reflection coeficients used was on its exact form and was generated for models I. n and m in multi-azimutaVincidenceangle s and contaminated by gaussian noise. In the inverse modeling we work with GA and with DE algoritms to calculate the inversion parameter (5 elastic stiffness and mass density for bottom media and Vs of upper isotropic media) by minimization of 12norm of difference between the "true" and sinthetic reflection coeficient. Weassume that we knwon the parameter of the upper media of the three models. except Vs for model one in especial case of inversion of upper media. The parameter to be determined by inverse modeling are parametrized in model space for values that is in acording with the value of the observed velocity of propagation of elastic waves in the earthcrust. and the resolution of measure, and constraints of elastic estability of the solid media. The GA and DE algoritms reached good inversion to the three models with at least three azimutal angles (0°,45° and 9Oj and incidence angles of 34° for model I, and 50° inverted only by GA for models 11and m; and the especial case take by DE tbat need at least 44° to invert the model I with the Vs of the upper media. From this resuhs we can see the potential to determine from qP-wave the elastic parameter of ftactured media and contribute to improve the use of azimutal AVO analises / Mestrado / Mestre em Engenharia de Petróleo

Orientador: Reginaldo Guirardello / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-07-28T22:45:04Z (GMT). No. of bitstreams: 1 Pissinatto_LucileneBeatriz_M.pdf: 2748990 bytes, checksum: 6d466001eb5fd7ba0fc5a7aec4d05338 (MD5) Previous issue date: 2001 / Resumo: O desenvolvimento adequado do layout de uma planta química é muito importante por motivos de economia, operação e segurança. A otimização global da planta deve incluir a otimização do posicionamento dos componentes e o comprimento e distribuição da tubulação, satisfazendo as restrições de segurança. O objetivo deste trabalho foi o desenvolvimento de modelos matemáticos mais eficientes, do ponto de vista computacional, para realizar a otimização do posicionamento dos componentes do layout de plantas químicas, através da formulação do problema em programação mista, inteira e linear (MILP), resolvendo-se com o auxílio do pacote comercial de otimização, GAMS/CPLEX. Um modelo mais eficiente implica em menor tempo de cálculo e resultados satisfatórios. Para testar os modelos foram propostos três casos de estudo, com diferentes números de componentes: uma unidade petroquímica (6 componentes); uma unidade de polimerização (7 componentes) e uma planta de recuperação de mercúrio a partir de efluentes eletrolíticos (9 componentes). O tempo computacional para resolver cada caso está intimamente ligado ao número de componentes, ou seja, quanto maior o número, tanto maior o tempo gasto para sua resolução. Os modelos testados variaram de caso para caso, isto é, com a adição de algumas restrições lógicas e/ou substituição de outras restrições relacionadas aos bocais e tubos, o tempo computacional diminuiu para a unidade petroquímica e para a planta de recuperação de mercúrio, no entanto, o modelo original (existente na literatura) apresentou melhores resultados para a planta de polimerização. Portanto, verificou-se que os modelos propostos são melhores em algumas situações, mas não em outras, dependendo do caso de estudo e das condições envolvidas (número de variáveis, número de iterações e nós inspecionados pelo algoritmo B&B, etc.) / Abstract: The adequate development of the layout of a chemical plant is very important because of economic, operational and safety reasons. The overall optimization of the plant must inc1ude the optimization of the placement of the components and the length and distribution of the pipes, satisfying safety constraints. The objective of this work was the development of more efficient mathematical models to accomplish the optimization of the placement of the components of the layout of chemical plants, through the formulation of the problem as mixed integer linear programming (MILP), solved with the aid of a commercial optimization package, GAMS/CPLEX. A more efficient model implies in smaller time of calculation and satisfactory results. To test the models, three case studies were proposed, with different numbers of components: a petrochemical unit (6 components); a suspension polymerization plant (7 components) and a recovery of mercury from the effluents of electrolytic plant (9 components). The computational time to solve each case is intimately linked to the number of components, in other words, as larger the number, so much larger the time spent for its resolution. The tested models showed mixed results, that is, with the addition of some logical constraints or substitution of other constraints related to the nozz1es and pipes, the computational time decreased for the petrochemical unit and for the plant of mercury recovery, but the original model (existent in the literature) presented better results for the suspension polymerization plant. Therefore, it was verified that the proposed models are better in some situations, but not in others, depending on the case study and of the involved conditions (number of variables, number of iterations and nodes inspected by the algorithm B&B, etc.) / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

Orientador: Jose Mario Martinez, Lucio T. Santos / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Matematica e Computação Cientifica / Made available in DSpace on 2018-07-24T10:01:35Z (GMT). No. of bitstreams: 1 Janesch_SilviaMartinideHolanda_D.pdf: 1787960 bytes, checksum: 222f3effcf1ad7501b5a884a0133db4b (MD5) Previous issue date: 1998 / Resumo: Apresentamos resultados gerais de penalização externa e exata. Estendemos o teorema clássico de penalização exata para o caso onde os subproblemas penalizados permanecem restritos. Introduzimos um algoritmo para resolver problemas de programação não linear baseado na função de penalização exata Li, onde penalizamos somente as restrições não lineares. Para resolver os subproblemas penalizados não suaves desenvolvemos um algoritmo de região de confiança. Ilustramos o método de penalização com região de confiança através de exemplos simples. Testes numéricos comparando o método de penalização com região de confiança com o algoritmo BOXQUACAN foram efetuados em 3 conjuntos de problemas. Abordamos o problema global de Lennard-Jones e propomos gerar bons pontos iniciais para este problema usando a solução de um subproblema restrito. / Abstract: We present the classical results for the exact and the exterior penalty problems. We extend the classic exact penalty function theorem for the case where the penalty subproblems remain constrained. We introduce an algorithm for solving nonlinear programming problems based on the L1 exact penalty function for which only the nonlinear constraints are penalized. For solving the nonsmooth penalty subproblems we develop a trust region algorithm. We illustrate the penalty method with trust region with simple examples. Numerical experiments comparing the penalty method with BOX-QUACAN algorithm were realized in three sets of problems. We attack Lennard Jones's global problem and we propose to generate good starting points for this problem using the solution of a constrained subproblem. / Doutorado / Doutor em Matemática Aplicada

Orientador: Jose Claudio Geromel / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de Computação / Made available in DSpace on 2018-07-25T07:21:34Z (GMT). No. of bitstreams: 1 Oliveira_MauricioCarvalhode_D.pdf: 9690579 bytes, checksum: cc4e519831f3c303a6a99c35e4695a0a (MD5) Previous issue date: 1999 / Resumo: Esta dissertação é dedicada ao estudo de métodos de análise de sistemas dinâmicos e de projeto de filtros e controladores para sistemas lineares, tendo como ferramenta básica as desigualdades matriciais lineares. Apresentamos uma pequena revisão de resultados conhecidos para análise de sistemas dinâmicos e desenvolvemos novas condições de estabilidade. Estes novos resultados nos permitem obter condições suficientes de estabilidade para diversos problemas, como a análise de estabilidade de sistemas lineares incertos com o auxílio de funções de Lyapunov dependentes dos parâmetros. Descrevemos parametrizações convexas para o projeto de filtros e controladores por realimentação do estado e por realimentação dinâmica da saída. Estas parametrizações são utilizadas na resolução de diversos problemas de filtragem e controle que ainda permaneciam sem solução como, por exemplo, o problema de controle positivo-real, o problema de filtragem descentralizada e robusta para sistemas incertos em donúnios poliedrais e o problema de controle H2 com seleção simultânea de atuadores e sensores. Outros problemas abordadosno texto são o controle robusto e descentralizadopor realimentação dinâmica da saída e o problema de controle com múltiplos objetivos. O texto é ilustrado com exemplos numéricos e simulações. ...Observação: O resumo, na íntegra, poderá ser visualizado no texto completo da tese digital / Abstract: This thesis is devoted to the analysis of dynamic systems and to the synthesis of filters and controllers for linear systems, all on the Linear Matrix Inequalities (LMI) framework. Wereview some well knownresults on dynamic systems analysisand developnew stability conditions. The new results enable us to draw sufficient stability conditions for several open problems, as the stability analysis of linear uncertain systems with parameter dependentLyapunov functions. We describe a set of convex parametrizations for the synthesis of filters and controllers under state and dynamic output feedback. These parametrizations provide solutions to several open problems as, for instance. the synthesis of positive-real control systems, the robust and decentralized filtering of uncertain polyhedral systems and the H2 control problem with simultaneous selection of sensors and actuators. Other problems addressed in the text are the robust and decentralized dynamic output feedback control problem and the multi-objective control problem. The text is illustrated by examples and simulations. ...Note: The omplete abstract is available with the full electronic digital thesis or dissertations / Doutorado / Doutor em Engenharia Elétrica

Orientador: Marko Antonio Rojas Medar / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-07-26T00:10:52Z (GMT). No. of bitstreams: 1 Santos_LucelinaBatistados_M.pdf: 6423420 bytes, checksum: b157a52e8273a28e48d361bc3e5ea356 (MD5) Previous issue date: 2000 / Resumo: Neste trabalho, um teorema de alternativa do tipo Gordan é utilizado no estudo de condições necessárias de otimalidade para o problema clássico de programação não linear (finito-dimensional). Condições suficientes são obtidas através de uma noção de convexidade generalizada (chamada invexidade). Além disso, sem hipóteses de convexidade (generalizada ou não) são obtidas condições suficientes de otimalidade via método de deformação. Resultados análogos são válidos para o problema de tempo contínuo (exceto o método de deformação). / Abstract: In this work, an alternative theorem (Gordon's type) is used to investigate necessary conditions of optimality for the classic nonlinear programming problem (finite dimensional). Sufficient conditions are obteined via a notion of generalized convexity (called invexity) . Moreover, without hipothesis of convexity (generalized or no) are obtained sufficient conditions via deformation method. Analogous results are valid for the continous time problem (except for the method of deformation). / Mestrado / Mestre em Matemática Aplicada

Orientador: Jose Mario Martinez Perez / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Científica / Made available in DSpace on 2018-07-14T15:46:51Z (GMT). No. of bitstreams: 1 Cetolo_LuciadeFatima_M.pdf: 1353850 bytes, checksum: ce075ee0ac0586d698c28b062b21a601 (MD5) Previous issue date: 1989 / Resumo: Este trabalho mostra os efeitos negativos do ajuste de funções quadráticas, quando as mesmas são utilizadas na minimização de funções em que suas primeiras derivadas não estão disponíveis, e a avaliação da função a ser minimizada é obtida experimentalmente, estando portanto, sujeita a erros de medição. No capítulo 1 fazemos uma descrição de um método para minimização de funções com as características acima, cuja proposta foi defendida em [6|, e consiste do ajuste de uma função quadrática a alguns pontos da função a ser minimizada e minimização da quadrática aproximada. No capítulo 2 apresentamos um contra-exemplo mostrando os efeitos negativos na abordagem do método proposto em [6]. No capítulo 3 formulamos o problema de maneira genérica e apresentamos as experiências realizadas com funções quadráticas como funções teste, na intenção de mostrar a confiabilidade do modelo em estudo. Com o intuito de realizar experiências com funções teste próximas da realidade, ou seja, funções com avaliação contendo certo erro de medição, realizamos no capítulo 4 experiências com funções quadráticas com perturbação. Finalmente, apresentamos no capítulo 5 as experiências numéricas realizadas com funções não quadráticas encontradas na literatura [8], visando aproximar o modelo ainda mais da realidade. / Abstract: Not informed. / Mestrado / Mestre em Matemática Aplicada