Problemas inversos em engenharia financeira: regularização com critério de entropia

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.

Estudio de las ecuaciones de Boussinesq para la modelización de tsunamis

[es]Podemos encontrar las ecuaciones de Boussinesq en la descripción de playas, rios y lagos. Estas ecuaciones estudian la dinámica de las aguas poco profundas como las ecuaciones “ Korteweg-deVries (KdV)". Sin embargo, a pesar de ser más conocidas, las ecuaciones de KdV, no son capaces de modelar olas solitarias propagándose en distintas direcciones. Entre muchas otras aplicaciones de las ecuaciones de Boussinesq destaca la de modelar olas de tsunamis. Estos tipos de olas ya son perfectamente descritos por las ecuaciones de Navier Stokes, pero todavía no existen técnicas que permitan resolverlas en un dominio tridimensional. Para ello se usan las ecuaciones de Boussinesq, pensadas como una simplificación de las ecuaciones de Navier Stokes. Los años 1871 y 1872 fueron muy importantes para el desarrollo de las ecuaciones de Boussinesq. Fue en 1871 cuando Valentin Joseph Boussinesq recibió el premio de la “Academy of Sciences”, por su trabajo dedicado a las aguas poco profundas. Ahí fue donde Boussinesq introdujo por primera vez los efectos dispersivos en las ecuaciones de Saint-Venant. Por ello, se puede decir que las ecuaciones de Boussinesq son más completas físicamente que las ecuaciones de Saint-Venant. Las ecuaciones de Boussinesq contienen una estructura hiperbólica (al igual que las ecuaciones no lineales de aguas poco profundas) combinada con derivadas de orden elevado para modelar la dispersión de la ola. Las ecuaciones de Boussinesq pueden aparecer de muchas formas distintas. Dependiendo de como hayamos escogido la variable de la velocidad podemos obtener un modelo u otro. El caso más usual es escoger la variable velocidad en un nivel del agua arbitrario. La efectividad de la ecuación de Boussinesq seleccionada variará dependiendo de la dispersión. Una buena elección de la variable velocidad puede mejorar significativamente la modelización de la propagación de ondas largas. Formalmente, como veremos en el capítulo 1, podemos transformar términos de orden elevado en términos de menor orden usando las relaciones asintóticas. Esto nos proporciona una forma elegante de mejorar las relaciones de dispersi\'on. Las ecuaciones de Boussinesq más conocidas son las que resolveremos en el capítulo 2. En dicho capítulo veremos la ecuación cúbica de Boussinesq, que sirve para describir el movimiento de ondas largas en aguas poco profundas; las ecuaciones de Boussinesq acopladas, que describen el movimiento de dos fluidos distintos en aguas poco profundas (como puede ser el caso de un barco que desprende accidentalmente aceite, el aceite va creando una capa que flota encima de la superficie del agua); la ecuación de Boussinesq estándar, que describe un gran número de fenómenos de olas dispersivas no lineales como la propagaci\ón en ambas direcciones de olas largas en la superficie de aguas poco profundas. Pero en olas de longitud de onda corta presenta una inestabilidad y la ecuación es incorrecta para el problema de Cauchy, por ello Bogolubsky propuso la ecuación de Boussinesq mejorada. Esta ecuación es la última que estudiaremos en el capítulo 2 y es una ecuación físicamente estable, correcta para el problema de Cauchy y además como veremos en el capítulo 3, apropiada para las simulaciones numéricas. Como ya indicado, en el capi tulo 1 deduciremos las ecuaciones de Boussinesq a partir de las ecuaciones físicas del flujo potencial. El objetivo principal es deducir dos modelos de ecuaciones de Boussinesq acopladas y obtener su relación de dispersión. Para llegar a ello, se usa un método de la expansión asintótica de la velocidad potencial en términos de un pequeño parámetro. De esta manera conseguimos dos modelos distintos, cada uno asociado a uno de los dos modelo de disipación que hemos establecido. Por último dado que las ecuaciones siempre vienen dadas en variables dimensionales, volveremos a la notación dimensional para analizar la relación de dispersión de las ecuaciones de Boussinesq disipativas. En el capí tulo 2 pasaremos a su resolución analítica, buscando soluciones de tipo solitón. Introduciremos el método de la tangente hiperbólica, muy útil para encontrar soluciones exactas de ecuaciones no lineales. Usaremos este método para resolver la ecuación cúbica de Boussinesq, un sistema de ecuaciones acopladas de Boussinesq, la ecuación estandar de Boussinesq y la mejorada. Los sistemas que aparecen en la aplicación del método de la tangente hiperbólica estan resueltos usando el software Mathematica y uno de ellos irá incluido en el apéndice A. En el capíulo 3 se introduce un esquema en diferencias finitas, que sirve para convertir problemas de ecuaciones diferenciales en problemas algebraicos fácilmente resolubles numéricamente. Este método nos ayudaráa estudiar la estabilidad y a resolver la ecuación mejorada de Boussinesq numéricamente en dos ejemplos distintos. En el apéndice B incluiremos el programa para la resolución numérica del primer ejemplo con el Mathematica.

Problemas direto e inverso de processos de separação em leito móvel simulado mediante mecanismos cinéticos de adsorção

Diversas aplicações industriais relevantes envolvem os processos de adsorção, citando como exemplos a purificação de produtos, separação de substâncias, controle de poluição e umidade entre outros. O interesse crescente pelos processos de purificação de biomoléculas deve-se principalmente ao desenvolvimento da biotecnologia e à demanda das indústrias farmacêutica e química por produtos com alto grau de pureza. O leito móvel simulado (LMS) é um processo cromatográfico contínuo que tem sido aplicado para simular o movimento do leito de adsorvente, de forma contracorrente ao movimento do líquido, através da troca periódica das posições das correntes de entrada e saída, sendo operado de forma contínua, sem prejuízo da pureza das correntes de saída. Esta consiste no extrato, rico no componente mais fortemente adsorvido, e no rafinado, rico no componente mais fracamente adsorvido, sendo o processo particularmente adequado a separações binárias. O objetivo desta tese é estudar e avaliar diferentes abordagens utilizando métodos estocásticos de otimização para o problema inverso dos fenômenos envolvidos no processo de separação em LMS. Foram utilizados modelos discretos com diferentes abordagens de transferência de massa, com a vantagem da utilização de um grande número de pratos teóricos em uma coluna de comprimento moderado, neste processo a separação cresce à medida que os solutos fluem através do leito, isto é, ao maior número de vezes que as moléculas interagem entre a fase móvel e a fase estacionária alcançando assim o equilíbrio. A modelagem e a simulação verificadas nestas abordagens permitiram a avaliação e a identificação das principais características de uma unidade de separação do LMS. A aplicação em estudo refere-se à simulação de processos de separação do Baclofen e da Cetamina. Estes compostos foram escolhidos por estarem bem caracterizados na literatura, estando disponíveis em estudos de cinética e de equilíbrio de adsorção nos resultados experimentais. De posse de resultados experimentais avaliou-se o comportamento do problema direto e inverso de uma unidade de separação LMS visando comparar os resultados obtidos com os experimentais, sempre se baseando em critérios de eficiência de separação entre as fases móvel e estacionária. Os métodos estudados foram o GA (Genetic Algorithm) e o PCA (Particle Collision Algorithm) e também foi feita uma hibridização entre o GA e o PCA. Como resultado desta tese analisouse e comparou-se os métodos de otimização em diferentes aspectos relacionados com o mecanismo cinético de transferência de massa por adsorção e dessorção entre as fases sólidas do adsorvente.

Solução de problemas inversos de transferência radiativa em meios heterogêneos unidimensionais e uma e duas camadas utilizando o algoritmo dos vagalumes

Esta tese apresenta um estudo sobre modelagem computacional onde são aplicadas meta-heurísticas de otimização na solução de problemas inversos de transferência radiativa em meios unidimensionais com albedo dependente da variável óptica, e meios unidimensionais de duas camadas onde o problema inverso é tratado como um problema de otimização. O trabalho aplica uma meta-heurística baseada em comportamentos da natureza conhecida como algoritmo dos vagalumes. Inicialmente, foram feitos estudos comparativos de desempenho com dois outros algoritmos estocásticos clássicos. Os resultados encontrados indicaram que a escolha do algoritmo dos vagalumes era apropriada. Em seguida, foram propostas outras estratégias que foram inseridas no algoritmo dos vagalumes canônico. Foi proposto um caso onde se testou e investigou todas as potenciais estratégias. As que apresentaram os melhores resultados foram, então, testadas em mais dois casos distintos. Todos os três casos testados foram em um ambiente de uma camada, com albedo de espalhamento dependente da posição espacial. As estratégias que apresentaram os resultados mais competitivos foram testadas em um meio de duas camadas. Para este novo cenário foram propostos cinco novos casos de testes. Os resultados obtidos, pelas novas variantes do algoritmo dos vagalumes, foram criticamente analisados.

Aplicação de Inteligência Computacional para a Solução de Problemas Inversos de Transferência Radiativa em Meios Participantes Unidimensionais

Esta pesquisa consiste na solução do problema inverso de transferência radiativa para um meio participante (emissor, absorvedor e/ou espalhador) homogêneo unidimensional em uma camada, usando-se a combinação de rede neural artificial (RNA) com técnicas de otimização. A saída da RNA, devidamente treinada, apresenta os valores das propriedades radiativas [ω, τ0, ρ1 e ρ2] que são otimizadas através das seguintes técnicas: Particle Collision Algorithm (PCA), Algoritmos Genéticos (AG), Greedy Randomized Adaptive Search Procedure (GRASP) e Busca Tabu (BT). Os dados usados no treinamento da RNA são sintéticos, gerados através do problema direto sem a introdução de ruído. Os resultados obtidos unicamente pela RNA, apresentam um erro médio percentual menor que 1,64%, seria satisfatório, todavia para o tratamento usando-se as quatro técnicas de otimização citadas anteriormente, os resultados tornaram-se ainda melhores com erros percentuais menores que 0,04%, especialmente quando a otimização é feita por AG.

Estimativas a priori para problemas não lineares de condução de calor, com o uso da transformada de Kirchhoff.

Este trabalho apresenta uma estimativa a priori para o limite superior da distribuição de temperatura considerando um problema em regime permanente em um corpo com uma condutividade térmica dependente da temperatura. A discussão é realizada supondo que as condições de contorno são lineares (lei de Newton do resfriamento) e que a condutividade térmica é constante por partes (quando considerada como uma função da temperatura). Estas estimativas consistem em uma ferramenta poderosa que pode prescindir da necessidade de uma simulação numérica cara de um problema de transferência de calor não linear, sempre que for suficiente conhecer o valor mais alto de temperatura. Nestes casos, a metodologia proposta neste trabalho é mais eficaz do que as aproximações usuais que assumem tanto a condutividade térmica quanto as fontes de calor como constantes.

Simulação de problemas de transferência de calor em regime permanente com uma relação entre condutividade térmica e temperatura constante por partes.

Este trabalho estuda a transferência de calor por condução considerando a condutividade térmica como uma função constante por partes da temperatura. Esta relação, embora fisicamente mais realista que supor a condutividade térmica constante, permite obter uma forma explícita bem simples para a inversa da Transformada de Kirchhoff (empregada para tratar a não linearidade do problema). Como exemplo, apresenta-se uma solução exata para um problema com simetria esférica. Em seguida, propôe-se uma formulação variacional (com unicidade demonstrada) que introduz um funcional cuja minimização é equivalente à solução do problema na forma forte. Finalmente compara-se uma solução exata obtida pela inversa da Transformada de Kirchhoff com a solução obtida via formulação variacional.

A Região Noroeste Fluminense e o Projeto Rio Rural: Tendências e problemas da agricultura familiar

O objetivo deste trabalho é analisar o processo de regionalização do noroeste fluminense através das ações verticalidade sob controle do Estado. Verifica-se especificamente, o posicionamento da sociedade civil organizada e as relações de poder envolvidas diante de tais inferencias racionalizadoras do espaço, e de que forma esse posicionamento pode redimencionar o ordenamento espacial que se impõe através das hegemonias de poder. Entre as abordagens ao tema, o trabalho considera, primeiramente, o processo de regionalização do noroeste fluminense por meio de relações socioespaciais pretéritas a fundação da região enquanto unidade administrativa. Para além das atuais fronteiras da região em questão, o trabalho destaca o controle do Estado pelas frentes catequizadoras e colonizadoras, sobre o território dos Puri, Coroado e Coropó. A segunda intermediação do trabalho ao tema diz respeito à ação de atores locais, ligados principalmente a grandes pecuaristas, na luta por autonomia econômica em relação a Campos dos Goytacazes; que em suas relações de poder com o Estado do Rio de Janeiro desencadeiam a instituição político-administrativa do noroeste fluminense. Sucessivamente o trabalho grifa a atribuição de características regionais a área de estudo como: pobreza rural e representatividade da agricultura familiar, relacionadas a ações de intelectuais orgânicos locais para direcionar verbas do Governo Federal ao setor agropecuário. A terceira parte do trabalho incorpora ao tema a descentralização da gestão municipal/regional por parte de políticas de planejamento territorial como o Rio Rural. Visto que ações de verticalidade do Governo do Estado do Rio de Janeiro, através do Projeto Rio Rural, intervém no ordenamento sócio-econômico-espacial da região principalmente através do modo de produção da agricultura familiar, a problemática se estenderia ao posicionamento político dessa categoria, enquanto sociedade civil organizada; para consideração final das tendências estabelecidas pela agricultura familiar ao ordenamento espacial regional

Análise dos problemas direto e inverso em um modelo de evolução populacional através de transformações integrais e inferência bayesiana

Modelos de evolução populacional são há muito tempo assunto de grande relevância, principalmente quando a população de estudo é composta por vetores de doenças. Tal importância se deve ao fato de existirem milhares de doenças que são propagadas por espécies específicas e conhecer como tais populações se comportam é vital quando pretende-se criar políticas públicas para controlar a sua proliferação. Este trabalho descreve um problema de evolução populacional difusivo com armadilhas locais e tempo de reprodução atrasado, o problema direto descreve a densidade de uma população uma vez conhecidos os parâmetros do modelo onde sua solução é obtida por meio da técnica de transformada integral generalizada, uma técnica numérico-analítica. Porém a solução do problema direto, por si só, não permite a simulação computacional de uma população em uma aplicação prática, uma vez que os parâmetros do modelo variam de população para população e precisam, portanto, ter seus valores conhecidos. Com o objetivo de possibilitar esta caracterização, o presente trabalho propõe a formulação e solução do problema inverso, estimando os parâmetros do modelo a partir de dados da população utilizando para tal tarefa dois métodos Bayesianos.

Otimização global topográfica aplicada a problemas de otimização com restrições de igualdade e desigualdade

Os métodos de otimização que adotam condições de otimalidade de primeira e/ou segunda ordem são eficientes e normalmente esses métodos iterativos são desenvolvidos e analisados através da análise matemática do espaço euclidiano n-dimensional, o qual tem caráter local. Esses métodos levam a algoritmos iterativos que são usados para o cálculo de minimizadores globais de uma função não linear, principalmente não-convexas e multimodais, dependendo da posição dos pontos de partida. Método de Otimização Global Topográfico é um algoritmo de agrupamento, o qual é fundamentado nos conceitos elementares da teoria dos grafos, com a finalidade de gerar bons pontos de partida para os métodos de busca local, com base nos pontos distribuídos de modo uniforme no interior da região viável. Este trabalho tem como objetivo a aplicação do método de Otimização Global Topográfica junto com um método robusto e eficaz de direções viáveis por pontos-interiores a problemas de otimização que tem restrições de igualdade e/ou desigualdade lineares e/ou não lineares, que constituem conjuntos viáveis com interiores não vazios. Para cada um destes problemas, é representado também um hiper-retângulo compreendendo cada conjunto viável, onde os pontos amostrais são gerados.

Problemas da gestāo dos stocks de camarāo

The worldwide shrimp landings in 1988 were reported to be 2.484.000 tons an increase of about 460.000 tons compared with 1985. The majority of shrimp fishing areas located in the tropical and sub-tropical regions with a contribution of more than 2 million tone. The most important species are the shallow water Penaeid shrimps. This quantity of about 2.5 million tons represents approximately 3% of the world marine catch. In terms of value, it represents almost 30% of the world trade in fish products. Main management objectives include: long term resources conservation; to maximize physical catches; to maximize the total income from catches or foreign exchange; to maximize economic profits; to reduce the shrimp by catch or improve its utilization; other social and economic interests. These objectives are in part interconnected and some are in conflict. Resource conservation is a basic condition for all other management objectives. Management policy definition should be in accordance with national goals and based on available scientific knowledge of the resources and of the fishing industry. The definition and selection of management objectives is an important process, to which the scientists must contribute the best way they can, mainly in providing the necessary information and options for management.