Using the Cayley-Hamilton theorem with N-partitioned matrices

A definition is given for the characteristic equation of anN-partitioned matrix. It is then proved that this matrix satisfies its own characteristic equation. This can then be regarded as a version of the Cayley-Hamilton theorem, of use withN-dimensional systems.

Sparse block-Jacobi matrices with arbitrarily accurate Hausdorff dimension

We show that the Hausdorff dimension of the spectral measure of a class of deterministic, i.e. nonrandom, block-Jacobi matrices may be determined with any degree of precision, improving a result of Zlatos [Andrej Zlatos,. Sparse potentials with fractional Hausdorff dimension, J. Funct. Anal. 207 (2004) 216-252]. (C) 2010 Elsevier Inc. All rights reserved.

Um modelo de antecedentes para a cocriação de valores em serviços de geriatria na cidade de São Paulo: uma aplicação da modelagem de equações estruturais

Dissertação apresentada ao Programa de Pós-graduação em Administração da Universidade Municipal de são Caetano do Sul

Métodos numéricos para a integração das equações da cinética pontual de um reator nuclear

As equações da cinétiica pontual de um reator nuclear térmico são integradas numericamente, utilizando um método matricial de continuação analitica. Essas equações são essencialmente não-negativas e possuem um autovalor dominante vinculado à reatividade do sistema. Também, descrevem-se os métodos de Hansen e Porsching.

Análise e simulação de um sistema de equações químicas do tipo difusão-reação

Neste trabalho estudamos um sistema de equações diferenciais parabólicas que modelam um processo de difusão-reação em duas dimensões da mistura molecular e reação química irreverssível de um só passo entre duas espécies químicas A e B para formar um produto P. Apresentamos resultados analíticos e computacionais relacionados à existência e unicidade da solução, assim como estimativas do erro local e global utilizando elementos finitos. Para os resultados analíticos usamos a teoria de semigrupos e o principio do m´aximo, e a simulação numérica é feita usando diferenças finitas centrais e o esquema simplificado de Ruge-Kutta. As estimativas do erro local para o problema semi-discretizado são estabelecidas usando normas de Sobolev, e para estimar o erro global usamos shadowing finito a posteriori. Os resultados computacionais obtidos mostram que o comportamento da solução está dentro do esperado e concorda com resultados da referências. Assim mesmo as estimativas do erro local e global são obtidas para pequenos intervalos de tempo e assumindo suficiente regularidade sobre a velocidade do fluído no qual realiza-se o processo. Destacamos que a estimativa do erro global usando shadowing finito é obtida sob hipóteses a posteriori sobre o operador do problema e o forte controle da velocidade numa vizinhança suficientemente pequena.

Solução de equações intervalares

Este trabalho trata do tipo de dado intervalar e da importância da especificação de uma semântica para garantir a correção e a interpretação coerente de resultados gerados, tais como de soluções de equações envolvendo este tipo de dado. Para tanto, realiza um estudo comparativo das semânticas de envoltória intervalar de reais e de número-intervalo, procurando identificar a influência de cada uma sobre definições fundamentais, tais como as das operações aritméticas e a do tipo de solução encontrado. Uma vez caracterizadas as semânticas associadas ao tipo de dado intervalar, o trabalho apresenta resultados que permitem mapear algebricamente a operação de multiplicação de números-intervalo tanto na representação de extremo inferior e extremo superior como na representação por ponto médio e diâmetro. Com base nesses resultados apresenta os mapeamentos das expressões algébricas que definem as potências positivas inteiras tanto para a semântica de número-intervalo como para a de envoltória de reais. Conjugando os resultados obtidos com a semântica de número-intervalo, o trabalho apresenta procedimentos algorítmicos para a determinação de dois tipos de soluções de equações intervalares: solução própria, a obtida diretamente a partir da relação de igualdade estrutural algébrica entre intervalos, e envoltória intervalar de soluções reais, normalmente referenciada como a solução intervalar usual. Exemplos são apresentados para a validação dos procedimentos, bem como para a discussão do significado de cada tipo de solução sob o enfoque semântico.

Cálculo da complexidade exata de algoritmos do tipo divisão-e-conquista através das equações características

A equação de complexidade de um algoritmo pode ser expressa em termos de uma equação de recorrência. A partir destas equações obtém-se uma expressão assintótica para a complexidade, provada por indução. Neste trabalho, propõem-se um esquema de solução de equações de recorrência usando equações características que são resolvidas através de um "software" de computação simbólica, resultando em uma expressão algébrica exata para a complexidade. O objetivo é obter uma forma geral de calcular a complexidade de um algoritmo desenvolvido pelo método Divisão-e-Conquista.

Solução de sistemas de equações algébrico-diferenciais ordinárias de índice superior

A modelagem matemática de problemas importantes e significativos da engenharia, física e ciências sociais pode ser formulada por um conjunto misto de equações diferenciais e algébricas (EADs). Este conjunto misto de equações deve ser previamente caracterizado quanto a resolubilidade, índice diferencial e condições iniciais, para que seja possível utilizar um código computacional para resolvê-lo numericamente. Sabendo-se que o índice diferencial é o parâmetro mais importante para caracterizar um sistema de EADs, neste trabalho aplica-se a redução de índice através da teoria de grafos, proposta por Pantelides (1988). Este processo de redução de índice é realizado numericamente através do algoritmo DAGRAFO, que transforma um sistema de índice superior para um sistema reduzido de índice 0 ou 1. Após esta etapa é necessário fornecer um conjunto de condições inicias consistentes para iniciar o código numérico de integração, DASSLC. No presente trabalho discute-se três técnicas para a inicialização consistente e integração numérica de sistemas de EADs de índice superior. A primeira técnica trabalha exclusivamente com o sistema reduzido, a segunda com o sistema reduzido e as restrições adicionais que surgem após a redução do índice introduzindo variáveis de restrição, e a terceira técnica trabalha com o sistema reduzido e as derivadas das variáveis de restrição. Após vários testes, conclui-se que a primeira e terceira técnica podem gerar um conjunto solução mesmo quando recebem condições iniciais inconsistentes. Para a primeira técnica, esta característica decorre do fato que no sistema reduzido algumas restrições, muitas vezes com significado físico importante, podem ser perdidas quando as equações algébricas são diferenciadas. Trabalhando com o sistema reduzido e as derivadas das variáveis de restrição, o erro da inicialização é absorvido pelas variáveis de restrição, mascarando a precisão do código numérico. A segunda técnica adotada não tem como absorver os erros da inicialização pelas variáveis de restrição, desta forma, quando as restrições adicionais não são satisfeitas, não é gerada solução alguma. Entretanto, ao aplicar condições iniciais consistentes para todas as técnicas, conclui-se que o sistema reduzido com as derivadas das variáveis restrição é o método mais conveniente, pois apresenta melhor desempenho computacional, inclusive quando a matriz jacobiana do sistema apresenta problema de mau condicionamento, e garante que todas as restrições que compõem o sistema original estejam presentes no sistema reduzido.

Derivação e solução de equações modelo da dinâmica de gases rarefeitos

Neste trabalho, duas equações modedlo na área da dinâmica de gases rarefeitos, são derivadas a partir de algumas soluções exatas da equação linearizada de Boltzmann homogênea e não homogênea. Em adição, uma versão analítca do método de ordenadas discretas é usado para resolver problemas clássicos nesta área, descritos pelo "Modelo S". Resultados numéricos são apresentados para os problemas de fluxo de Couette, fluxo de Poiseuille, "Creep" Térmico, Deslizamento Térmico e problema de Kramers.