Aspectos algébricos de sistemas dinâmicos

Este trabalho trata o problema genérico da obtenção analítica exata das variedades algébricas que definem domínios de estabilidade e multiestabilidade para sistemas dinâmicos dissipativos com equações de movimento definidas por funções racionais. Apresentamos um método genérico, válido para qualquer sistema dinâmico, que permite reduzir a análise de sistemas multidimensionais arbitrários à análise de um sistema unidimensional equivalente. Este método é aplicado ao mapa de Hénon, o exemplo paradigmático de sistema multidimensional, para estudar a estrutura aritmética imposta pela dinâmica das órbitas de períodos 4, 5, e 6, bem como seus domínios de estabilidade no espaço de parâmetros. Graças à obtençao de resultados analíticos exatos, podemos explorar pela primeira vez as peculariedades de cada um dos períodos mencionados. Algumas das novidades mais marcantes encontradas são as seguintes: Para período 4, encontramos um domínio de multiestabilidade caracterizado pela coexistência de duas órbitas definidas em corpos algébricos distintos. Observamos a existência de discontinuidades na dinâmica simbólica quando os parâmetros são mudados adiabáticamente ao longo de circulações fechadas no espaço de parâmetros e explicamos sua origem algébrica. Publicamos tais resultados em dois artigos: Physica A, 295, 285-290(2001) e Physical Review E, 65, 036231 (2002). Para período 5, obtivemos a variedade algébrica que define o "camarão" (shrimp) característico, obtemos uma expressão analítica para todas as órbitas de período 5, classificamos todas as singulariedades presentes no espaço de parâmetros e analisamos todas as mudanças que ocorrem ao circular-se em torno de tais singulariedades. Para período 6, da expressão analítica que fornece todas as órbitas, encontramos um resultado muito surpreendente, o mais notável desta dissertação: a possibilidade de coexistência de órbitas reais e complexas estáveis, para valores reais dos parâmetros físicos. Resultados preliminares parecem indicar serem tais órbitas complexas uma espécie de órbitas fantasmas, com semelhanças as órbitas encontradas por Gutzwiller para sistemas Hamiltonianos (não- dissipativos).

Solução LTSn para problemas de transferência radiativa com polarização em geometria plana

O método LTSN tem sido utilizado na resolução de uma classe abrangente de problemas de transporte de partículas neutras que são reduzidos a um sistema linear algébrico depois da aplicação da transformada de Laplace. Na maioria dos casos estudados os autovalores associados são reais e simétricos. Para o problema de criticalidade os autovalores associados são reais ou imaginários puros e simétricos, e para o o problema de multigrupo podem aparecer autovalores complexos. O objetivo deste trabalho consiste na generalização da formulação LTSN para problemas de transporte com autovalores complexos. Por esse motivo é focada a solução de um problema radiativo de transporte com polarização em uma placa plana. A solução apresentada fundamenta-se na aplicação da transformada de Laplace ao conjunto de equações SN dos problemas resultantes da decomposição da equação de transferência radiativa com polarização em série de Fourier, seguindo o procedimento de Chandrasekhar. Esse procedimento gera 2L + 2 sistemas lineares de ordem 4N dependentes do parâmetro complexo "s". Aqui, L é o grau de anisotropia e N a ordem de quadratura. A solução desse sistema simbólico é obtida através da aplicação da transformada inversa de Laplace depois da inversão da matriz simbólica pelo método da diagonalização. Para a obtenção das constantes de integração é assumido que os componentes do vetor de Stokes são reais e as matrizes dos autovalores e autovetores são separadas em suas partes real e imaginária. A solução LTSN para autovalores complexos é validada através da comparação da solução para uma placa com espessura unitária, grau de anisotropia L = 13, albedo de espalhamento simples $ = 0:99, coe ciente de re exão de Lambert ¸0 = 0:1 e N = 150, segundo dados da literatura consultada.

Proteção passiva de estruturas e sistemas elétricos contra incêndio aplicada a uma unidade de refino de petróleo

Nesta dissertação visa-se estudar e propor alternativas de solução para a proteção de estruturas e sistemas elétricos contra fogo numa unidade de craqueamento catalítico de uma refinaria de petróleo, por meio de proteção passiva. A proteção passiva tem por finalidade garantir a integridade das estruturas sujeitas a incêndio, durante um determinado período de tempo, para possibilitar, no caso da refinaria, a realização de procedimentos de parada da unidade de forma segura e controlar o incêndio a fim de diminuir a possibilidade de propagação do fogo para outras áreas. Com base em técnicas de análise de riscos fez-se a identificação de zonas potencialmente sujeitas a cenários de acidente envolvendo jato de fogo e/ou incêndio em poça. A delimitação das áreas onde haveria necessidade de proteção passiva foi realizada com base em modelos para jatos de fogo e incêndio em poça já estabelecidos na literatura. O dimensionamento da proteção passiva de estruturas e sistemas elétricos com o uso de diversos materiais usados comercialmente para este fim foi estimado com base em equações empíricas desenvolvidas por Jeanes, 1980, Stanzak, 1973 e PABCO, 1984, e, para alguns casos particulares foi feita uma verificação por solução numérica da equação da condução do calor em meio sólido.Assim, foram determinados quais os materiais mais adequados em cada caso de aplicação e qual a espessura em que deve ser aplicado para que a temperatura no elemento estrutural ou no sistema elétrico não atinja a sua determinada temperatura crítica em um período de tempo pré-determinado. Para os casos de elementos estruturais como colunas de sustentação da unidade de seção cilíndrica, o principal material para proteção passiva é a argamassa projetada e para perfil I, é o emprego de placas de gesso. Já para o caso de sistemas elétricos, podem ser utilizadas tanto tintas intumescentes quanto as mantas reforçadas com fibras minerais, esta escolha depende da geometria do sistema em que será empregado. Da comparação entre estes dois métodos pode-se concluir que o dimensionamento da proteção passiva fazendo o uso das correlações empíricas é menos conservativo que para o caso do uso da equação da difusão do calor resolvida por método numérico. Porém, os resultados diferem dentro de um limite considerado aceitável (em torno de 15%) levando-se em consideração os erros embutidos em cada método de cálculo. É importante mencionar que as correlações empíricas são de mais simples aplicação por possuir apenas operações matemáticas básicas. Usando as correlações empíricas para os perfis cilíndricos de aço (diâmetro de 0,1524 m e espessura de parede de 0,0254 m), a espessura de revestimento estimada com o uso das correlações empíricas necessária para garantir que a temperatura na interface entre os dois materiais não atinja 550°C em duas horas seria de 13,5 mm para argamassa projetada, 19,7 mm para vermiculita com silicato de sódio e 34,5 mm para recobrimento com concreto com proteção do tipo contorno. Fazendo o mesmo cálculo pelo método numérico proposto, os resultados foram de 15,53 mm para argamassa projetada, 22,06 mm para vermiculita com silicato de sódio e 38,98 mm para recobrimento com concreto com proteção do tipo contorno. Fazendo o mesmo cálculo pelo método numérico proposto, os resultados foram de 15,53 mm para argamassa projetada, 22,06 mm para vermiculita com silicato de sódio e 38,98 mm para recobrimento com concreto com proteção do tipo contorno. Cabe ressaltar que com a realização desta dissertação busca-se uma integração entre o mestrado acadêmico e o meio empresarial com o desenvolvimento de trabalhos de natureza acadêmica que tenham aplicação direta na prática. Espera-se assim permitir que Universidade dê retorno à sociedade que a mantém e propiciar que setores da sociedade possam usufruir da capacidade disponível na academia.

Sincronização idêntica em redes de osciladores caóticos

A presente dissertação versa sobre a sincronização idêntica em redes de osciladores caóticos. Uma perspectiva razoavelmente histórica sobre a literatura da área é apresentada . O conceito de caos é introduzido junto com outras idéias da dinâmica não-linear: sistemas dinâmicos, exemplos de sistemas, atratores, expoentes de Liapunov, etc. A integração numérica de equações diferenciais é largamente utilizada, principalmente, para o cálculo de expoentes e o desenho do diagrama de fases. A sincronização idêntica é definida, inicialmente, em redes que não passam de um par de osciladores. A variedade de sincronização (conjunto de pontos no espaço de fases no qual a solução do sistema é encontrada se há sincronização) é determinada. Diferentes variantes de acoplamentos lineares são enfocadas: acoplamento interno, externo, do tipo mestre-escravo e birecional, entre outras. Para detectar sincronização, usa-se o conceito de expoente de Liapunov transversal, uma extensão do conceito clássico de expoente de Liapunov que caracteriza a sincronização como a existência de um atrator na variedade de sincronização. A exposição é completada com exemplos e atinge relativo detalhe sobre o assunto, sem deixar de ser sintética com relação à ampla literatura existente. Um caso de sincronização em antifase que usa a mesma análise é incluído. A sincronização idêntica também é estudada em redes de osciladores idênticos com mais de dois osciladores. As possibilidades de sincronização completa e parcial são explanadas. As técnicas usadas para um par de osciladores são expandidas para cobrir este novo tipo de redes. A existência de variedades de sincronização invariantes é considerada como fator determinante para a sincronização. A sincronização parcial gera estruturas espaciais, analisadas sob a denominação de padrões. Algumas relações importantes entre as sincronizações são explicitadas, principalmente as degenerescências e a relação entre a sincronização parcial e a sincronização completa do respectivo estado sincronizado para alguns tipos de acoplamento. Ainda são objetos de interesse as redes formadas por grupos de osciladores idênticos que são diferentes dos osciladores dos outros grupos. A sincronização parcial na qual todos os grupos de osciladores têm seus elementos sincronizados é chamada de sincronização primária. A sincronização secundária é qualquer outro tipo de sincronização parcial. Ambas são exemplificadas e analisadas por meio dos expoentes transversais e novamente por meio da existência de invariantes de sincronização. Obtém-se, então, uma caracterização suficientemente ampla, completada por casos específicos.

Projetos de controladores para sistemas de potência utilizando LMI'S

A presente dissertação tem como objetivo estudar e aprimorar métodos de projetos de controladores para sistemas de potência, sendo que esse trabalho trata da estabilidade dinâmica de sistemas de potência e, portanto, do projeto de controladores amortecedores de oscilações eletromecânicas para esses sistemas. A escolha dos métodos aqui estudados foi orientada pelos requisitos que um estabilizador de sistemas de potência (ESP) deve ter, que são robustez, descentralização e coordenação. Sendo que alguns deles tiveram suas características aprimoradas para atender a esses requisitos. A abordagem dos métodos estudados foi restringida à análise no domínio tempo, pois a abordagem temporal facilita a modelagem das incertezas paramétricas, para atender ao requisito da robustez, e também permite a formulação do controle descentralizado de maneira simples. Além disso, a abordagem temporal permite a formulação do problema de projeto utilizando desigualdades matriciais lineares (LMI’s), as quais possuem como vantagem o fato do conjunto solução ser sempre convexo e a existência de algoritmos eficientes para o cálculo de sua solução. De fato, existem diversos pacotes computacionais desenvolvidos no mercado para o cálculo da solução de um problema de inequações matriciais lineares. Por esse motivo, os métodos de projeto para controladores de saída buscam sempre colocar o problema na forma de LMI’s, tendo em vista que ela garante a obtenção de solução, caso essa solução exista.

Identificação de sistemas contínuos utilizando o método dos momentos de Poisson

A maioria dos métodos de síntese e sintonia de controladores, bem como métodos de otimização e análise de processos necessitam de um modelo do processo em estudo. A identificação de processos é portanto uma área de grande importância para a engenharia em geral pois permite a obtenção de modelos empíricos dos processos com que nos deparamos de uma forma simples e rápida. Mesmo não utilizando leis da natureza, os modelos empíricos são úteis pois descrevem o comportamento específico de determinado processo. Com o rápido desenvolvimento dos computadores digitais e sua larga aplicação nos sistemas de controle em geral, a identificação de modelos discretos foi amplamente desenvolvida e empregada, entretanto, modelos discretos não são de fácil interpretação como os modelos contínuos pois a maioria dos sistema com que lidamos são de representação contínua. A identificação de modelos contínuos é portanto útil na medida que gera modelos de compreensão mais simples. A presente dissertação estuda a identificação de modelos lineares contínuos a partir de dados amostrados discretamente. O método estudado é o chamado método dos momentos de Poisson. Este método se baseia em uma transformação linear que quando aplicada a uma equação diferencial ordinária linear a transforma em uma equação algébrica evitando com isso a necessidade do cálculo das derivadas do sinais de entrada e saída Além da análise detalhada desse método, onde demonstramos o efeito de cada parâmetro do método de Poisson sobre o desempenho desse, foi realizado também um estudo dos problemas decorrentes da discretização de sinais contínuos, como por exemplo o efeito aliasing decorrente da utilização de tempos de amostragem muito grandes e de problemas numéricos da identificação de modelos discretos utilizando dados com tempos de amostragem muito pequenos de forma a destacar as vantagens da identificação contínua sobre a identificação discreta Também foi estudado um método para compensar a presença de offsets nos sinais de entrada e saída, método esse inédito quando se trata do método dos momentos de Poisson. Esse trabalho também comprova a equivalência entre o método dos momentos de Poisson e uma metodologia apresentada por Rolf Johansson em um artigo de 1994. Na parte final desse trabalho são apresentados métodos para a compensação de erros de modelagem devido à presença de ruído e distúrbios não medidos nos dados utilizados na identificação. Esses métodos permitem que o método dos momentos de Poisson concorra com os métodos de identificação discretos normalmente empregados como por exemplo ARMAX e Box-Jenkins.

Planejamento de perturbações para a identificação de modelos dinâmicos

A identificação de modelos é determinante no sucesso das modernas técnicas de controle avançado de processos. Um modelo para o sistema pode ser obtido através de modelagem rigorosa, baseada em equações governantes do sistema ou através da modelagem empírica e estimação de parâmetros. Embora mais rápida e fácil, a modelagem empírica necessita de alguns testes de identificação nos quais as variáveis manipuladas são variadas de modo que resultem em variações nas variáveis controladas. Os testes de identificação podem apresentar custos muito elevados tendo em vista que o sistema pode sair de seu ponto normal de operação, gerando produtos com folga de especificação. Este fato ocorre porque usualmente as perturbações aplicadas nas variáveis manipuladas nas indústrias de processos são independentes umas das outras, aumentando a duração do teste de identificação. Desta forma, neste trabalho foi desenvolvida uma nova metodologia de projeto de perturbações simultâneas para a identificação de modelos dinâmicos baseada na direcionalidade do sistema, com o objetivo de fornecer dados mais ricos para se capturar corretamente o comportamento multivariável do sistema e manter o processo no ponto de operação normal. As perturbações são projetadas conforme as características de um modelo simplificado do processo, ou pré-modelo. Este modelo inicial é obtido essencialmente de dados históricos de planta, selecionados através de uma sistemática análise de correlação desenvolvida neste trabalho A metodologia proposta é composta de duas partes: a primeira parte diz respeito à análise dos dados históricos de planta para obtenção de informações prelimirares as quais são utilizadas no planejamento de perturbações, tais como amplitude do ruído de medida, correlação entre as variáveis de processo, constante de tempo do sistema e matriz de ganhos. E a segunda parte consiste no cálculo da amplitude das perturbações baseado nos resultados da primeira etapa do planejamento. Para sistemas mal-condicionados verificou-se que as perturbações planejadas pela metodologia removem menos a planta de seu ponto de operação gerando resultados mais consistentes em relação às perturbações tradicionais. Já para sistemas bem-condicionados, os resultados são semelhantes. A metodologia foi aplicada em uma unidade piloto experimental e numa unidade de destilação da PETROBRAS, cujos resultados apontam pouca remoção dos sistemas do ponto de operação e modelos consistentes. A validação dos modelos também foi contemplada na dissertação, uma vez que foi proposto um novo critério de validação que considera a derivada dos dados de planta e a do modelo e não apenas os dados de planta e os dados da simulação das saídas do modelo.

Solução de uma classe de equações diferenciais parciais em espaços de Niko'skii através do método de Galerkin

Neste trabalho apresenta-se uma solu c~ao para um problema abstrato de Cauchy. Basicamente, d a-se uma formula c~ao abstrata para certos tipos de equa c~oes diferenciais parciais n~ao lineares de evolu c~ao em espa cos de Nikol'skii, tais espa cos possuem boas propriedades de regularidade e resultados de imers~ao compacta, num certo sentido s~ao intermedi arios entre os espa cos de Holder e os espa cos de Sobolev. Aplicando o m etodo de Galerkin, prova-se resultados de exist^encia global de solu c~oes fracas, como tamb em a exist^encia de solu c~oes fracas com a propriedade de reprodu c~ao. E impondo mais hip oteses sobre os operadores envolvidos demonstra-se unicidade de solu c~oes fracas.

Análise dos antecedentes à resistência a sistemas empresariais: uma abordagem explano-exploratória

O objetivo desse estudo explano-exploratório é determinar os antecedentes à resistência a sistemas empresariais. Com esse intuito, por meio de levantamento bibliográfico: (a) conceituou-se o comportamento de resistência, (b) identificaram-se os fatores que influenciam o comportamento de resistência, (c) definiu-se o objeto do comportamento de resistência, (d) elaborou-se um meta-modelo contendo os fatores que influenciam o comportamento de resistência e (e) desenvolveu-se o questionário estruturado como instrumento de coleta de dados. Durante o período de aplicação do instrumento de coleta de dados foram obtidos 169 questionários válidos, preenchidos por gestores de tecnologia da informação. A partir dos dados obtidos, a pesquisa valeu-se das técnicas de análise fatorial exploratória, análise fatorial confirmatória e modelagem de equações estruturais para confirmar o modelo proposto por teorias do campo de gerenciamento de sistemas de informação. No entanto, verificou-se que o meta-modelo de resistência a sistemas empresariais, apresentando: (a) as características idiossincráticas das pessoas ou grupos, (b) as características dos sistemas a serem implantados e (c) a interação das características pessoais e do sistema com o contexto organizacional, subdividida em interação sócio-técnica e interação poder e política, não se confirmou na realidade estudada. Dessa forma, após modificações no meta-modelo inicial, identificou-se um modelo onde os fatores “Sistemas” e “Inclinação Pessoal” explicaram aproximadamente 49% da variância do fenômeno de comportamento de resistência a sistemas empresariais.

Fatores antecedentes da satisfação do aluno e do uso de sistemas virtuais de aprendizagem

A educação a distância tem passado por grandes transformações, principalmente após o advento da internet e das tecnologias de informação e comunicação (TICs). Inúmeras perguntas sobre qualidade e resultados de aprendizagem em ambientes virtuais foram geradas com o crescimento da modalidade. Pesquisadores têm investigado métodos de avaliação dos benefícios promovidos pelo e-learning sob um número diversificado de perspectivas. O objetivo desta pesquisa é avaliar o impacto dos construtos qualidade do sistema, qualidade da informação e qualidade do serviço na satisfação do aluno e no uso de Sistemas Virtuais de Aprendizagem em ambientes de e-learning, utilizando como base teórica o modelo de Sucesso de e-learning, adaptado do modelo de Delone e McLean por Holsapple e Lee-Post. A metodologia de pesquisa tipo survey foi administrada por meio de um curso on-line ofertado a 291 estudantes de instituições públicas e privadas de todas as regiões do Brasil. Para o tratamento e análise dos dados, utilizaram-se técnicas de modelagem de equações estruturais e análise fatorial confirmatória. Os resultados demonstram que o uso do sistema é impactado pela variação dos construtos qualidade do sistema, qualidade da informação e qualidade dos serviços, já a satisfação do aluno é antecedida pela qualidade percebida da informação e do serviço. Muitos dos benefícios gerados pela educação a distância são causados pela satisfação do aluno e pela intensidade com que este utiliza o sistema de aprendizagem. Ao identificar os indicadores que antecedem estas variáveis, os gestores educacionais podem planejar seus investimentos visando atender às demandas mais importantes, além de utilizar a informação para lidar com um dos maiores problemas em EaD: a evasão.

Controle preditivo de sistemas mistos lógicos-dinâmicos

O comportamento dinâmico de um sistema é tradicionalmente descrito por um modelo, em geral associado a um conjunto de equações diferenciais ou de equações de diferenças, onde as variáveis representam grandezas físicas. Sistemas complexos, principalmente na indústria de processos, incorporam elementos com comportamento dinâmico lógico, tais como atuadores e sensores ON-OFF (estados aberto/fechado) ou proposições lógicas (estados verdadeiro/falso). Estes sistemas são denominados “Sistemas Híbridos”, “Sistemas Mistos Lógicos-dinâmicos” ou, simplesmente, “Sistemas Mistos”. Neste trabalho, são apresentadas técnicas que, associando variáveis lógicas a estes elementos, e mediante a incorporação de restrições sobre as variáveis, permitem obter um modelo matemático do sistema misto. Neste caso, técnicas clássicas de controle não permitem a incorporação destas novas variáves e restrições. Como opção de controle de sistemas mistos, é então proposta e estudada uma técnica de controle preditivo baseado em modelo. São apresentados tanto a formulação teórica do problema de controle, quanto exemplos e simulações bem como um estudo de caso de sua aplicação sobre um sistema de equalização de uma planta de tratamento de efluentes.

Integração semi-Lagrangeana num modelo oceano-atmosfera e ajuste geostrófico no domínio tempo

Este trabalho visa o uso da função de Green de valor inicial no ajuste geostrófico e do método Semi-Lagrangeano na integração de um modelo acoplado oceano-atmosfera descrito pelas equações de águas rasas. O ajuste geostrófico é considerado atravées de perturbações na pressão e do vento. No caso de sistemas sem rotação, é discutida a relação da equação hidrostática com ondas longas não-dispersivas. Com rotação, a conservação da vorticidade potencial permite escolher a elevação correspondente a um estado de equilíbrio geostrófico. O sistema de equações de águas rasas é desacoplado em equações de Klein-Gordon com valores iniciais e termos não-homogêneos acoplados. A resposta dinâmica formada pela resposta transiente e a resposta forçada é obtida para uma perturbação inicial da elevação. A ação do vento como forçante nas equações de momento 2D, através do transporte de Eckman, conduz a uma equação de águas rasas forçada. Uma decomposição da resposta forçada é realizada com uma resposta permanente, que satisfaz a equação de Helmholtz , e com o uso da base dinâmica gerada pela resposta impulso. Um modelo hidrodinâmico 3D introduzido por Casulli e governado por equações não-lineares de águas rasas é integrado na vertical para a obtenção de um modelo 2D. Com isto, as condições de contorno devido a tensão do vento e a fricção devido a topografia do fundo, transformam-se em forçantes do modelo. O modelo foi integrado com um método semi-implícito em diferenças finitas, utilizando-se o método Semi-Lagrangeano para a parte advectiva. Simulações simbólicas foram realizadas para o ajuste geostrófico devido a perturbações de duração infinita e finita para a elevação e para o efeito da tensão do vento. Foram realizadas simulações numéricas para variadas geometrias, em particular a Baia de Guanabara e a Lagoa do Patos.

Localização de faltas de alta impedância em sistemas de distribuição de energia : uma metodologia baseada em redes neurais artificiais

O presente trabalho apresenta uma nova metodologia de localização de faltas em sistemas de distribuição de energia. O esquema proposto é capaz de obter uma estimativa precisa da localização tanto de faltas sólidas e lineares quanto de faltas de alta impedância. Esta última classe de faltas representa um grande problema para as concessionárias distribuidoras de energia elétrica, uma vez que seus efeitos nem sempre são detectados pelos dispositivos de proteção utilizados. Os algoritmos de localização de faltas normalmente presentes em relés de proteção digitais são formulados para faltas sólidas ou com baixa resistência de falta. Sendo assim, sua aplicação para localização de faltas de alta impedância resulta em estimativas errôneas da distância de falta. A metodologia proposta visa superar esta deficiência dos algoritmos de localização tradicionais através da criação de um algoritmo baseado em redes neurais artificiais que poderá ser adicionado como uma rotina adicional de um relé de proteção digital. O esquema proposto utiliza dados oscilográficos pré e pós-falta que são processados de modo que sua localização possa ser estimada através de um conjunto de características extraídas dos sinais de tensão e corrente. Este conjunto de características é classificado pelas redes neurais artificiais de cuja saída resulta um valor relativo a distância de falta. Além da metodologia proposta, duas metodologias para localização de faltas foram implementadas, possibilitando a obtenção de resultados comparativos. Os dados de falta necessários foram obtidos através de centenas de simulações computacionais de um modelo de alimentador radial de distribuição. Os resultados obtidos demonstram a viabilidade do uso da metodologia proposta para localização de faltas em sistemas de distribuição de energia, especialmente faltas de alta impedância.

Redes de mapas acoplados em sistemas biológicos de reação-difusão

Em modelos em que a distribuição espacial da população não é con- siderada, isto é, quando se supõe que haja uma homogeneidade espacial, e se estuda a evolução temporal do sistema, há uma única variável independente: o tempo. Caso a população seja constituída de duas espécies, do tipo parasitóide-hospedeiro, e a variável independente tempo for considerada discreta, teremos um sistema de equações a diferenças, como por exemplo o modelo de Nicholson-Bailey cujas soluções são apresentadas neste trabalho. Populações espacialmente distribuídas, em um espaço de natureza discreta, juntamente com a dinâmica vital em tempo discreto, têm o seu comportamento estudado através de redes de mapas acoplados. Após estudar o modelo de Hassell (dinâmica vital de Nicholson-Bailey com movimentação por difusão) e o modelo planta-herbívoro com movimentação por taxia, deduzimos e simulamos um modelo incluindo movimentação por taxia, difusão e convecção. É também apresentado neste trabalho, um paralelo entre estes modelos de redes de mapas acoplados e aqueles com as equações diferenciais correspondentes.

Comparação de sistemas de programação

Este trabalho compara as soluções disponibilizadas pelos sistemas Derive 5.0, Maple 6 e Mathematica 4.0 para problemas que encontramos no ensino secundário e também nos primeiros anos da universidade. Procuramos destacar os aspectos distintos entre cada um dos programas ao mesmo tempo que fazemos referência aos pontos em que tudo se passa de forma semelhante. Esta dissertação aborda o cálculo numérico, o cálculo simbólico, a programação e os gráficos. Para cada um dos assuntos é estudada a forma como se podem resolver os problemas através dos três sistemas comparando-se estas soluções. Inicialmente, é feita uma abordagem que permite ao utilizador adquirir os conhecimentos básicos acerca dos diversos programas. Tratamos de seguida de algumas questões relacionadas com o cálculo numérico e com algumas funções nomeadamente da Teoria dos Números. Referimos listas e funções e são analisadas diversas formas de manipular listas e os seus elementos bem como algumas áreas da Análise Matemática das quais destacamos as equações, a derivação e a integração compreendendo cálculo numérico e cálculo simbólico. Examinamos um vasto conjunto de operações definidas sobre matrizes (representadas como listas de listas) e polinómios que abrangem as operações mais comuns de cada um dos campos. Analisamos também a programação recursiva, a programação imperativa, a programação funcional e a programação por regras de reescrita. A abordagem aqui adoptada foi a de fornecer ao utilizador as construções chave mais importantes que cada paradigma de programação utiliza bem como as informações básicas acerca do funcionamento de cada uma delas de modo a permitir a resolução dos problemas propostos. Por último os gráficos sobre os quais incidiu a nossa análise foram os de uma e de duas variáveis representados no referencial cartesiano, gráficos estes que são os mais utilizados quer ao nível do ensino superior quer ao nível do ensino secundário. A qualidade e a facilidade de obter rapidamente as representações dão outra dimensão ao estudo dos gráficos principalmente quando estamos a falar de gráficos a três dimensões. A ideia de animação gráfica é também aqui abordada sendo evidente os benefícios da utilização da mesma nos programas em que é possível efectuá-la. Concluímos que na programação o Mathematica destaca-se em relação aos demais o mesmo se passando no Maple no respeitante à representação gráfica. O Derive permite que durante o contacto inicial seja mais fácil trabalhar e aprender a linguagem própria.