Desenvolvimento de um controlador preditivo estocástico para processos da indústria química.

O sucesso de estratégias de controle preditivo baseado em modelo (MPC, na sigla em inglês) tanto em ambiente industrial quanto acadêmico tem sido marcante. No entanto, ainda há diversas questões em aberto na área, especialmente quando a hipótese simplificadora de modelo perfeito é abandonada. A consideração explícita de incertezas levou a importantes progressos na área de controle robusto, mas esta ainda apresenta alguns problemas: a alta demanda computacional e o excesso de conservadorismo são questões que podem ter prejudicado a aplicação de estratégias de controle robusto na prática. A abordagem de controle preditivo estocástico (SMPC, na sigla em inglês) busca a redução do conservadorismo através da incorporação de informação estatística dos ruídos. Como processos na indústria química sempre estão sujeito a distúrbios, seja devido a diferenças entre planta e modelo ou a distúrbios não medidos, está técnica surge como uma interessante alternativa para o futuro. O principal objetivo desta tese é o desenvolvimento de algoritmos de SMPC que levem em conta algumas das especificidades de tais processos, as quais não foram adequadamente tratadas na literatura até o presente. A contribuição mais importante é a inclusão de ação integral no controlador através de uma descrição do modelo em termos de velocidade. Além disso, restrições obrigatórias (hard) nas entradas associadas a limites físicos ou de segurança e restrições probabilísticas nos estados normalmente advindas de especificações de produtos também são consideradas na formulação. Duas abordagens foram seguidas neste trabalho, a primeira é mais direta enquanto a segunda fornece garantias de estabilidade em malha fechada, contudo aumenta o conservadorismo. Outro ponto interessante desenvolvido nesta tese é o controle por zonas de sistemas sujeitos a distúrbios. Essa forma de controle é comum na indústria devido à falta de graus de liberdade, sendo a abordagem proposta a primeira contribuição da literatura a unir controle por zonas e SMPC. Diversas simulações de todos os controladores propostos e comparações com modelos da literatura são exibidas para demonstrar o potencial de aplicação das técnicas desenvolvidas.

De rios e colinas: a cidade de São Paulo entre os séculos XVI a XVIII. Um estudo sobre a tradição urbanística de origem lusitana e suas transformações à época do iluminismo

Esta tese tem como objetivo analisar os processos históricos envolvidos na criação e difusão de uma cultura urbanística lusitana, a qual se formou alinhavando e reinterpretando inúmeras referências e que as transformou continuamente ao longo dos séculos e de sua expansão territorial. Tomamos o século XVI como referência da expansão lusitana e o caso de São Paulo de Piratininga tida pela bibliografia especializada como uma excepcionalidade na história urbanística e colonial portuguesa como exemplo de alinhamento com a tradição urbanística portuguesa. Essa inserção é visível através do cotejamento da história urbana de outras cidades portuguesas (como Coimbra, Sintra, Tomar, mas, sobretudo, Lisboa). De outro lado, contempla-se o século XVIII como um momento de profundas mudanças nessa tradição urbanística e a inserção de novas mentalidades, práticas, motivadas principalmente pelo surgimento de novos movimentos (como o iluminismo) e protagonistas ou, ainda, suas novas faces (como o Estado português, em processo de centralização, de aproximação com os interesses mercantis e de medidas modernizadoras). Procura-se, por fim, demonstrar a complexidade dos processos culturais que envolvem, ao longo dos séculos, a criação e transformação das cidades, bem como a intensa circulação de ideias, pessoas, práticas entre a Europa e as Américas.

Contexto e modularização em redes neurais recorrentes para aprendizagem de seqüências temporais

Este trabalho apresenta um sistema neural modular, que processa separadamente informações de contexto espacial e temporal, para a tarefa de reprodução de sequências temporais. Para o desenvolvimento do sistema neural foram considerados redes neurais recorrentes, modelos estocásticos, sistemas neurais modulares e processamento de informações de contexto. Em seguida, foram estudados três modelos com abordagens distintas para aprendizagem de seqüências temporais: uma rede neural parcialmente recorrente, um exemplo de sistema neural modular e um modelo estocástico utilizando a teoria de modelos markovianos escondidos. Com base nos estudos e modelos apresentados, esta pesquisa propõe um sistema formado por dois módulos sucessivos distintos. Uma rede de propagação direta (módulo estimador de contexto espacial) realiza o processamento de contexto espacial identificando a seqüência a ser reproduzida e fornecendo um protótipo do contexto para o segundo módulo. Este é formado por uma rede parcialmente recorrente (módulo de reprodução de sequências temporais) para aprender as informações de contexto temporal e reproduzir em suas saídas a seqüência identificada pelo módulo anterior. Para a finalidade mencionada, este mestrado utiliza a distribuição de Gibbs na saída do módulo para contexto espacial de forma que este forneça probabilidades de contexto espacial, indicando o grau de certeza do módulo e possibilitando a utilização de procedimentos especiais para os casos de dúvida. O sistema neural foi testado em conjuntos contendo trajetórias abertas, fechadas, e com diferentes situações de ambigüidade e complexidade. Duas situações distintas foram avaliadas: (a) capacidade do sistema em reproduzir trajetórias a partir de pontos iniciais treinados; e (b) capacidade de generalização do sistema reproduzindo trajetórias considerando pontos iniciais ou finais em situações não treinadas. A situação (b) é um problema de difícil ) solução em redes neurais devido à falta de contexto temporal, essencial na reprodução de seqüências. Foram realizados experimentos comparando o desempenho do sistema modular proposto com o de uma rede parcialmente recorrente operando sozinha e um sistema modular neural (TOTEM). Os resultados sugerem que o sistema proposto apresentou uma capacidade de generalização significamente melhor, sem que houvesse uma deterioração na capacidade de reproduzir seqüências treinadas. Esses resultados foram obtidos em sistema mais simples que o TOTEM.

Mineralogia do diamante da região do Alto Araguaia, Mato Grosso

Este trabalho apresenta as propriedades morfológicas, físicas e químicas do diamante da região do Alto Araguaia, situada na borda norte da Bacia do Paraná. Os cristais estudados provém dos garimpos dos rios Garças, Araguaia, Caiapó e seus tributários menores. As dimensões dos cristais embora variáveis concentram-se no intervalo 2-15 mm. Ao contrário do que se observa no Alto Paranaíba, são raros os achados de grande porte. Os fragmentos de clivagem são abundantes, indicando transporte prolongado. Cerca de metade dos indivíduos cristalinos são incolores, sendo os demais castanhos, amarelos, verdes e cinzas. Outras cores como rosa e violeta ocorrem excepcionalmente. A morfologia dos cristais é complexa e variada, caracterizando-se pelo acentuado predomínio do hábito rombododecaédrico, grande número de geminados de contato e ausência de formas cúbicas. As faces exibem grau variável de curvatura e diversos padrões de microestruturas resultantes da dissolução provocada por agentes oxidantes. O ataque parece ocorrer durante a colocação do kimberlito, quando o alívio de pressão eleva a temperatura favorecendo a corrosão. As microestruturas mais freqüentes são degraus, colinas e micro-discos em (110), depressões triangulares eqüiláteras em (111), e depressões quadráticas em (100). Alguns diamantes contém inclusões cristalinas, sendo as mais comuns olivina (forsterita), granada (piropo) e o próprio diamante. Esses minerais apresentam-se geralmente idiomorfos e circundados por anomalias ópticas (birrefringência anômala), evidenciando o caráter epigenético dessas inclusões. Anàlogamente ao que ocorre nos kimberlitos africanos e siberianos, essas inclusões são minerais característicos de altas pressões e temperaturas, sendo a olivina a variedade mais freqüente, seguida pelo diamante e granada. Essa é uma evidência de que o diamante do Alto Araguaia provém de matrizes ultrabásicas (kimberlíticas). As inclusões negras (carvões) são extremamente comuns e parecem constituir defeitos do retículo cristalino (clivagens internas, deslocamentos estruturais), resultando talvez de impactos sofridos pelos cristais durante as diversas fases de transporte. O comportamento ao infravermelho indica que a maior parte dos cristais (85%) contém impurezas de nitrogênio (tipo I), sob forma de placas submicroscópicas paralelas a (100) (tipo Ia), ou átomos dispersos no retículo cristalino (tipo Ib). Os espécimes desprovidos de nitrogênio (tipo II) são relativamente raros (6%), sendo os demais intermediários (9%) entre I e II. Além do nitrogênio, os espectros infravermelhos revelaram que alguns diamantes contém hidrogênio. A presença deste elemento juntamente com o carbono e nitrogênio, sugere uma derivação a partir de sedimentos ricos em matéria orgânica, os quais teriam sido incorporados ao magma kimberlítico por correntes de convecção subcrostais. Outras impurezas menores (elementos traços) são: Al, Ca, Si, Mg, Fe, Cu, Cr e Co. Estes elementos são os mesmos encontrados nos diamantes dos kimberlitos africanos. Observado sob luz ultravioleta, 36% dos cristais exibem cores de fluorescência, sendo as mais comuns o azul e verde, e mais raramente amarelo e castanho. Entre as variedades fluorescentes, 27% são também fosforescentes, sendo a cor mais comum o azul. Não foi observada nenhuma relação entre o comportamento luminescente e as demais propriedades estudadas. O peso específico varia entre 3,500 a 3,530. As principais causas dessa variação são as impurezas de nitrogênio, as inclusões minerais e alguns defeitos cristalinos. Os resultados desta pesquisa indicam claramente que as propriedades do diamante do Alto araguaia são semelhantes às dos diamantes africanos e siberianos, originários de matrizes kimberlíticas. Nenhum dado foi obtido que sugerisse uma origem \"sui generis\" para o material estudado. Em relação à localização das possíveis chaminés dispersoras dos cristais, duas alternativas se apresentam: 1) poderiam estar relacionadas ao magmatismo Neocretáceo, responsável por numerosos focos vulcânicos próximos à área; ou 2) ligadas a episódios mais antigos ocorridos no Pré-Cambriano. Em qualquer das alternativas, poderiam estar recobertas por sedimentos mais recentes ou alteradas superficialmente, sendo que, no caso de uma idade pré-cambriana, há ainda a possibilidade destas matrizes terem sido totalmente destruídas.

Aplicação de métodos geoestatísticos no estudo das distribuições espaciais de condutividade hidráulica em áreas do Sistema Aquífero Guarani (SAG) e do Sistema Aquífero Bauru (SAB) no Estado de São Paulo

A condutividade hidráulica (K) é um dos parâmetros controladores da magnitude da velocidade da água subterrânea, e consequentemente, é um dos mais importantes parâmetros que afetam o fluxo subterrâneo e o transporte de solutos, sendo de suma importância o conhecimento da distribuição de K. Esse trabalho visa estimar valores de condutividade hidráulica em duas áreas distintas, uma no Sistema Aquífero Guarani (SAG) e outra no Sistema Aquífero Bauru (SAB) por meio de três técnicas geoestatísticas: krigagem ordinária, cokrigagem e simulação condicional por bandas rotativas. Para aumentar a base de dados de valores de K, há um tratamento estatístico dos dados conhecidos. O método de interpolação matemática (krigagem ordinária) e o estocástico (simulação condicional por bandas rotativas) são aplicados para estimar os valores de K diretamente, enquanto que os métodos de krigagem ordinária combinada com regressão linear e cokrigagem permitem incorporar valores de capacidade específica (Q/s) como variável secundária. Adicionalmente, a cada método geoestatístico foi aplicada a técnica de desagrupamento por célula para comparar a sua capacidade de melhorar a performance dos métodos, o que pode ser avaliado por meio da validação cruzada. Os resultados dessas abordagens geoestatísticas indicam que os métodos de simulação condicional por bandas rotativas com a técnica de desagrupamento e de krigagem ordinária combinada com regressão linear sem a técnica de desagrupamento são os mais adequados para as áreas do SAG (rho=0.55) e do SAB (rho=0.44), respectivamente. O tratamento estatístico e a técnica de desagrupamento usados nesse trabalho revelaram-se úteis ferramentas auxiliares para os métodos geoestatísticos.

Uma formulação por média-variância multi-período para o erro de rastreamento em carteiras de investimento.

Neste trabalho, deriva-se uma política de escolha ótima baseada na análise de média-variância para o Erro de Rastreamento no cenário Multi-período - ERM -. Referindo-se ao ERM como a diferença entre o capital acumulado pela carteira escolhida e o acumulado pela carteira de um benchmark. Assim, foi aplicada a metodologia abordada por Li-Ng em [24] para a solução analítica, obtendo-se dessa maneira uma generalização do caso uniperíodo introduzido por Roll em [38]. Em seguida, selecionou-se um portfólio do mercado de ações brasileiro baseado no fator de orrelação, e adotou-se como benchmark o índice da bolsa de valores do estado de São Paulo IBOVESPA, além da taxa básica de juros SELIC como ativo de renda fixa. Dois casos foram abordados: carteira composta somente de ativos de risco, caso I, e carteira com um ativo sem risco indexado à SELIC - e ativos do caso I (caso II).

Estimação do índice de memória em processos estocásticos com memória longa: uma abordagem via ABC

Neste trabalho propomos o uso de um método Bayesiano para estimar o parâmetro de memória de um processo estocástico com memória longa quando sua função de verossimilhança é intratável ou não está disponível. Esta abordagem fornece uma aproximação para a distribuição a posteriori sobre a memória e outros parâmetros e é baseada numa aplicação simples do método conhecido como computação Bayesiana aproximada (ABC). Alguns estimadores populares para o parâmetro de memória serão revisados e comparados com esta abordagem. O emprego de nossa proposta viabiliza a solução de problemas complexos sob o ponto de vista Bayesiano e, embora aproximativa, possui um desempenho muito satisfatório quando comparada com métodos clássicos.

Modelagem estocástica da dispersão axial: aplicação em um reator tubular de polimerização.

Reatores tubulares de polimerização podem apresentar um perfil de velocidade bastante distorcido. Partindo desta observação, um modelo estocástico baseado no modelo de dispersão axial foi proposto para a representação matemática da fluidodinâmica de um reator tubular para produção de poliestireno. A equação diferencial foi obtida inserindo a aleatoriedade no parâmetro de dispersão, resultando na adição de um termo estocástico ao modelo capaz de simular as oscilações observadas experimentalmente. A equação diferencial estocástica foi discretizada e resolvida pelo método Euler-Maruyama de forma satisfatória. Uma função estimadora foi desenvolvida para a obtenção do parâmetro do termo estocástico e o parâmetro do termo determinístico foi calculado pelo método dos mínimos quadrados. Uma análise de convergência foi conduzida para determinar o número de elementos da discretização e o modelo foi validado através da comparação de trajetórias e de intervalos de confiança computacionais com dados experimentais. O resultado obtido foi satisfatório, o que auxilia na compreensão do comportamento fluidodinâmico complexo do reator estudado.

Hybrid convex combinations for IIR system identification.

The low complexity of IIR adaptive filters (AFs) is specially appealing to realtime applications but some drawbacks have been preventing their widespread use so far. For gradient based IIR AFs, adverse operational conditions cause convergence problems in system identification scenarios: underdamped and clustered poles, undermodelling or non-white input signals lead to error surfaces where the adaptation nearly stops on large plateaus or get stuck at sub-optimal local minima that can not be identified as such a priori. Furthermore, the non-stationarity in the input regressor brought by the filter recursivity and the approximations made by the update rules of the stochastic gradient algorithms constrain the learning step size to small values, causing slow convergence. In this work, we propose IIR performance enhancement strategies based on hybrid combinations of AFs that achieve higher convergence rates than ordinary IIR AFs while keeping the stability.

Modelagem do crescimento e produção aplicado ao manejo florestal na Amazônia brasileira

A floresta Amazônica possui um papel ambiental, social e econômico importante para a região, para o país e para o mundo. Dessa forma, técnicas de exploração que visam a diminuição dos impactos causados à floresta são essenciais. Com isso, o objetivo dessa tese é comparar a Exploração de Impacto Reduzido com a Exploração Convencional na Amazônia brasileira através de modelos empíricos de árvore individual de crescimento e produção. O experimento foi instalado na fazenda Agrossete, localizada em Paragominas - PA. Em 1993, três áreas dessa fazenda foram selecionadas para exploração. Na primeira área, 105 hectares foram explorados através da Exploração de Impacto Reduzido. Na segunda área, 75 hectares foram submetidos à Exploração Convencional. E, por fim, a terceira área foi mantida como área testemunha. A coleta de dados de diâmetro à altura do peito e a identificação das espécies dentro de uma parcela de 24,5 hectares, instalada aleatoriamente em cada área, foi realizada nos anos de 1993 (antes da colheita), 1994 (seis meses depois da colheita), 1995, 1996, 1998, 2000, 2003, 2006 e 2009. Dessa forma, as três áreas foram comparadas através do ajuste de um modelo de incremento diamétrico, considerando que efeito estocástico podia assumir outras quatro distribuições além da distribuição normal, de um modelo de probabilidade de mortalidade e de um modelo de probabilidade de recrutamento. O comportamento do incremento diamétrico indicou que as áreas que foram submetidas a exploração possuem o mesmo comportamento em quase todos os grupos de espécies, com exceção do grupo de espécies intermediárias. Os indivíduos que são submetidos a exploração possuem um maior crescimento em diâmetros quando comparados com área que não sofreu exploração. Além disso, assumir o efeito estocástico com distribuição Weibull melhorou o ajuste dos modelos. Em relação à probabilidade de mortalidade, novamente as áreas que sofreram exploração possuem comportamento semelhante quanto à mortalidade, mas diferente da área que não foi explorada, sendo que os indivíduos localizados nas áreas exploradas possuem uma maior probabilidade de morte em relação aos presentes na área não explorada. Os modelos de probabilidade de recrutamento indicaram diferença apenas entre as áreas exploradas e a área controle. Sendo que, as áreas exploradas apresentaram uma maior taxa de recrumento em comparação a área não explorada. Portanto, o comportamento individual das árvores após a exploração é o mesmo na Exploração Convencional e na Exploração de Impacto Reduzido.

Efeito da corrosão do chapeamento do fundo do casco sobre a confiabilidade estrutural de navios.

A engenharia é a ciência que transforma os conhecimentos das disciplinas básicas aplicadas a fatos reais. Nosso mundo está rodeado por essas realizações da engenharia, e é necessário que as pessoas se sintam confortáveis e seguras nas mesmas. Assim, a segurança se torna um fator importante que deve ser considerado em qualquer projeto. Na engenharia naval, um apropriado nível de segurança e, em consequência, um correto desenho estrutural é baseado, atualmente, em estudos determinísticos com o objetivo de obter estruturas capazes de suportar o pior cenário possível de solicitações durante um período de tempo determinado. A maior parte das solicitações na estrutura de um navio se deve à ação da natureza (ventos, ondas, correnteza e tempestades), ou, ainda, aos erros cometidos por humanos (explosões internas, explosões externas e colisões). Devido à aleatoriedade destes eventos, a confiabilidade estrutural de um navio deveria ser considerada como um problema estocástico sob condições ambientais bem caracterizadas. A metodologia probabilística, baseada em estatística e incertezas, oferece uma melhor perspectiva dos fenômenos reais que acontecem na estrutura dos navios. Esta pesquisa tem como objetivo apresentar resultados de confiabilidade estrutural em projetos e planejamento da manutenção para a chapa do fundo dos cascos dos navios, as quais são submetidas a esforços variáveis pela ação das ondas do mar e da corrosão. Foram estudados modelos estatísticos para a avaliação da estrutura da viga-navio e para o detalhe estrutural da chapa do fundo. Na avaliação da estrutura da viga-navio, o modelo desenvolvido consiste em determinar as probabilidades de ocorrência das solicitações na estrutura, considerando a deterioração por corrosão, com base numa investigação estatística da variação dos esforços em função das ondas e a deterioração em função de uma taxa de corrosão padrão recomendada pela DET NORSKE VERITAS (DNV). A abordagem para avaliação da confiabilidade dependente do tempo é desenvolvida com base nas curvas de resistências e solicitações (R-S) determinadas pela utilização do método de Monte Carlo. Uma variação estatística de longo prazo das adversidades é determinada pelo estudo estatístico de ondas em longo prazo e ajustada por uma distribuição com base numa vida de projeto conhecida. Constam no trabalho resultados da variação da confiabilidade ao longo do tempo de um navio petroleiro. O caso de estudo foi simplificado para facilitar a obtenção de dados, com o objetivo de corroborar a metodologia desenvolvida.