Maston ympärille asennetun antenniryhmän säteilykuvion laskeminen

Tämän insinöörityön tavoitteena oli tehdä MATLAB-sovellus, jolla voidaan laskea ja tulostaa yleisradiotoiminnassa käytettävien dipolipaneeliantennien säteilykuvioita. Toinen tavoite oli saattaa tämän työn valvojan Antti Koivumäen vuoden 1980 diplomityön laskenta ja grafiikka modernimmalle tasolle. Aluksi työssä tarkasteltiin yleisiä asioita VHF- ja UHF-taajuuksien lähetystavoista ja lähetinantennityypeistä. Alussa käytiin läpi myös antennijoukon ja -ryhmän käsite. Seuraavaksi tutustuttiin tarkemmin lähetinantennien rakenteeseen ja ominaisuuksiin. Antennien ominaisuuksista tarkasteltiin säteilykuviota, suuntaavuutta, vahvistusta, hyötysuhdetta ja polarisaatiota. Näistä säteilykuvio oli tarkimman tutkinnan kohteena. Säteilykuvion laskeminen esitettiin kaavojen ja havainnoillistavien kuvien avulla. Lopuksi perehdyttiin antennin säteilykuvion laskemiseen tietokoneella, etenkin MATLABilla. Osuudessa tarkasteltiin graafisen käyttöliittymän tekemistä, laskennan suorittamista ja horisontaalisen säteilykuvion esittämistä graafisesti. Työn tuloksena saatiin tehdyksi MATLAB-sovellus, joka laskee elementtiantennin horisontaalisen säteilykuvion annettujen parametrien avulla ja tulostaa sen.

Metamorphism and kinetics in the Torres del Paine contact aureole

Contact aureoles provide an excellent geologic environment to study the mechanisms of metamorphic reactions in a natural system. The Torres del Paine (TP) intrusion is one of the most spectacular natural laboratories because of its excellent outcrop conditions. It formed in a period from 12.59 to 12.43 Ma and consists of three large granite and four smaller mafic batches. The oldest granite is on top, the youngest at the bottom of the granitic complex, and the granites overly the mafic laccolith. The TP intruded at a depth of 2-3 km into regional metamorphic anchizone to greenschist facies pelites, sandstones, and conglomerates of the Cerro Toro and Punta Barrosa formations. It formed a thin contact aureole of 150-400 m width. This thesis focuses on the reaction kinetics of the mineral cordierite in the contact aureole using quantitative textural analysis methods. First cordierite was formed from chlorite break¬down (zone I, ca. 480 °C, 750 bar). The second cordierite forming reaction was the muscovite break-down, which is accompanied by a modal decrease in biotite and the appearance of k- feldspar (zone II, 540-550 °C, 750 bar). Crystal sizes of the roundish, poikiloblastic cordierites were determined from microscope thin section images by manually marking each crystal. Images were then automatically processed with Matlab. The correction for the intersection probability of each crystal radius yields the crystal size distribution in the rock. Samples from zone I below the laccolith have the largest crystals (0.09 mm). Cordierites from zone II are smaller, with a maximum crystal radius of 0.057 mm. Rocks from zone II have a larger number of small cordierite crystals than rocks from zone I. A combination of these quantitative analysis with numerical modeling of nucleation and growth, is used to infer nucleation and growth parameters which are responsible for the observed mineral textures. For this, the temperature-time paths of the samples need to be known. The thermal history is complex because the main body of the intrusion was formed by several intrusive batches. The emplacement mechanism and duration of each batch can influence the thermal structure in the aureole. A possible subdivision of batches in smaller increments, so called pulses, will focus heat at the side of the intrusion. Focusing all pulses on one side increases the contact aureole size on that side, but decreases it on the other side. It forms a strongly asymmetric contact aureole. Detailed modeling shows that the relative thicknesses of the TP contact aureole above and below the intrusion (150 and 400 m) are best explained by a rapid emplacement of at least the oldest granite batch. Nevertheless, temperatures are significantly too low in all models, compared to observed mineral assemblages in the hornfelses. Hence, an other important thermal mechanisms needs to take place in the host rock. Clastic minerals in the immature sediments outside the contact aureole are hydrated due to small amounts of expelled fluids during contact metamorphism. This leads to a temperature increase of up to 50 °C. The origin of fluids can be traced by stable isotopes. Whole rock stable isotope data (6D and δ180) and chlorine concentrations in biotite document that the TP intrusion induced only very small amounts of fluid flow. Oxygen whole rock data show δ180 values between 9.0 and 10.0 %o within the first 5 m of the contact. Values increase to 13.0 - 15.0 %o further away from the intrusion. Whole rock 6D values display a more complex zoning. First, host rock values (-90 to -70 %o) smoothly decrease towards the contact by ca. 20 %o, up to a distance of ca. 150 m. This is followed by an increase of ca. 20 %o within the innermost 150 m of the aureole (-97.0 to -78 %o at the contact). The initial decrease in 6D values is interpreted to be due to Rayleigh fractionation accompanying the dehydration reactions forming cordierite, while the final increase reflects infiltration of water-rich fluids from the intrusion. An over-estimate on the quantity and the corresponding thermal effect yields a temperature increase of less than 30 °C. This suggests that fluid flow might have contributed only for a small amount to the thermal evolution of the system. A combination of the numerical growth model with the thermal model, including the hydration reaction enthalpies but neglecting fluid flow and incremental growth, can be used to numerically reproduce the observed cordierite textures in the contact aureole. This yields kinetic parameters which indicate fast cordierite crystallization before the thermal peak in the inner aureole, and continued reaction after the thermal peak in the outermost aureole. Only small temperature dependencies of the kinetic parameters seem to be needed to explain the obtained crystal size data. - Les auréoles de contact offrent un cadre géologique privilégié pour l'étude des mécanismes de réactions métamorphiques associés à la mise en place de magmas dans la croûte terrestre. Par ses conditions d'affleurements excellentes, l'intrusion de Torres del Paine représente un site exceptionnel pour améliorer nos connaissances de ces processus. La formation de cette intrusion composée de trois injections granitiques principales et de quatre injections mafiques de volume inférieur couvre une période allant de 12.50 à 12.43 Ma. Le plus vieux granite forme la partie sommitale de l'intrusion alors que l'injection la plus jeune s'observe à la base du complexe granitique; les granites recouvrent la partie mafique du laccolite. L'intrusion du Torres del Paine s'est mise en place a 2-3 km de profondeur dans un encaissant métamorphique. Cet encaissant est caractérisé par un métamorphisme régional de faciès anchizonal à schiste vert et est composé de pélites, de grès, et des conglomérats des formations du Cerro Toro et Punta Barrosa. La mise en place des différentes injections granitiques a généré une auréole de contact de 150-400 m d'épaisseur autour de l'intrusion. Cette thèse se concentre sur la cinétique de réaction associée à la formation de la cordiérite dans les auréoles de contact en utilisant des méthodes quantitatives d'analyses de texture. On observe plusieurs générations de cordiérite dans l'auréole de contact. La première cordiérite est formée par la décomposition de la chlorite (zone I, environ 480 °C, 750 bar), alors qu'une seconde génération de cordiérite est associée à la décomposition de la muscovite, laquelle est accompagnée par une diminution modale de la teneur en biotite et l'apparition de feldspath potassique (zone II, 540-550 °C, 750 bar). Les tailles des cristaux de cordiérites arrondies et blastic ont été déterminées en utilisant des images digitalisées des lames minces et en marquant individuellement chaque cristal. Les images sont ensuite traitées automatiquement à l'aide du programme Matlab. La correction de la probabilité d'intersection en fonction du rayon des cristaux permet de déterminer la distribution de la taille des cristaux dans la roche. Les échantillons de la zone I, en dessous du lacolite, sont caractérisés par de relativement grands cristaux (0.09 mm). Les cristaux de cordiérite de la zone II sont plus petits, avec un rayon maximal de 0.057 mm. Les roches de la zone II présentent un plus grand nombre de petits cristaux de cordiérite que les roches de la zone I. Une combinaison de ces analyses quantitatives avec un modèle numérique de nucléation et croissance a été utilisée pour déduire les paramètres de nucléation et croissance contrôlant les différentes textures minérales observées. Pour développer le modèle de nucléation et de croissance, il est nécessaire de connaître le chemin température - temps des échantillons. L'histoire thermique est complexe parce que l'intrusion est produite par plusieurs injections successives. En effet, le mécanisme d'emplace¬ment et la durée de chaque injection peuvent influencer la structure thermique dans l'auréole. Une subdivision des injections en plus petits incréments, appelés puises, permet de concentrer la chaleur dans les bords de l'intrusion. Une mise en place préférentielle de ces puises sur un côté de l'intrusion modifie l'apport thermique et influence la taille de l'auréole de contact produite, auréole qui devient asymétrique. Dans le cas de la première injection de granite, une modélisation détaillée montre que l'épaisseur relative de l'auréole de contact de Torres del Paine au-dessus et en dessous de l'intrusion (150 et 400 m) est mieux expliquée par un emplacement rapide du granite. Néanmoins, les températures calculées dans l'auréole de con¬tact sont trop basses pour que les modèles thermiques soient cohérants par rapport à la taille de cette auréole. Ainsi, un autre mecanisme exothermique est nécessaire pour permettre à la roche encais¬sante de produire les assemblages observés. L'observation des roches encaissantes entourant les granites montre que les minéraux clastiques dans les sédiments immatures au-dehors de l'auréole sont hydratés suite à la petite quantité de fluide expulsée durant le métamorphisme de contact et/ou la mise en place des granites. Les réactions d'hydratation peuvent permettre une augmentation de la température jusqu'à 50 °C. Afin de déterminer l'origine des fluides, une étude isotopique de roches de l'auréole de contact a été entreprise. Les isotopes stables d'oxygène et d'hydrogène sur la roche totale ainsi que la concentration en chlore dans la biotite indiquent que la mise en place des granites du Torres del Paine n'induit qu'une circulation de fluide limitée. Les données d'oxygène sur roche totale montrent des valeurs δ180 entre 9.0 et 10.0%o au sein des cinq premiers mètres du contact. Les valeurs augmentent jusqu'à 13.0 - 15.0 plus on s'éloigne de l'intrusion. Les valeurs 5D sur roche totale montrent une zonation plus complexe. Les valeurs de la roche encaissante (-90 à -70%o) diminuent progressivement d'environ 20%o depuis l'extérieur de l'auréole jusqu'à une distance d'environ 150 m du granite. Cette diminution est suivie par une augmentation d'environ 20%o au sein des 150 mètres les plus proches du contact (-97.0 à -78%o au contact). La diminution initiale des valeurs de 6D est interprétée comme la conséquence du fractionnement de Rayleigh qui accompagne les réactions de déshydratation formant la cordiérite, alors que l'augmentation finale reflète l'infiltration de fluide riche en eau venant de l'intrusion. A partir de ces résultats, le volume du fluide issu du granite ainsi que son effet thermique a pu être estimé. Ces résultats montrent que l'augmentation de température associée à ces fluides est limitée à un maximum de 30 °C. La contribution de ces fluides dans le bilan thermique est donc faible. Ces différents résultats nous ont permis de créer un modèle thermique associé à la for¬mation de l'auréole de contact qui intègre la mise en place rapide du granite et les réactions d'hydratation lors du métamorphisme. L'intégration de ce modèle thermique dans le modèle numérique de croissance minérale nous permet de calculer les textures des cordiérites. Cepen¬dant, ce modèle est dépendant de la vitesse de croissance et de nucléation de ces cordiérites. Nous avons obtenu ces paramètres en comparant les textures prédites par le modèle et les textures observées dans les roches de l'auréole de contact du Torres del Paine. Les paramètres cinétiques extraits du modèle optimisé indiquent une cristallisation rapide de la cordiérite avant le pic thermique dans la partie interne de l'auréole, et une réaction continue après le pic thermique dans la partie la plus externe de l'auréole. Seules de petites dépendances de température des paramètres de cinétique semblent être nécessaires pour expliquer les don¬nées obtenues sur la distribution des tailles de cristaux. Ces résultats apportent un éclairage nouveau sur la cinétique qui contrôle les réactions métamorphiques.

Sistema informàtic d’ajut a la detecció i al seguiment del melanoma

Mitjançant les tècniques de visió per computador aquest projecte pretén desenvolupar una aplicació capaç de segmentar la pell, detectar nevus (pigues i altres taques) i poder comparar imatges de pacients amb risc de contreure melanoma preses en moments diferents. Aquest projecte pretén oferir diferents eines informàtiques als dermatòlegs per a propòsits relacionats amb la investigació. L’ objectiu principal d’ aquest projecte és desenvolupar un sistema informàtic que proporcioni als dermatòlegs agilitat a l’hora de gestionar les dades dels pacients amb les sevesimatges corresponents, ajudar-los en la realització de deteccions dels nevus d’aquestes imatges, i ajudar-los en la comparació d’exploracions (amb les deteccions realitzades)de diferents èpoques d’un mateix pacient

Reconeixement i seguiment d’objectes en vídeo a partir del seu moviment

L’objectiu d’aquest PFC és estudiar la branca de la detecció d’objectes en vídeos segons el seu moviment. Per fer-ho es crearà un algorisme que sigui capaç de tractar un vídeo, calculant el nombre d’objectes de l’escena i quina és la posició de cada un d’aquests. L’algorisme ha de ser capaç de trobar un conjunt de regions útils i a partir d’aquest, separar-lo en diferents grups, cada un representant un objecte en moviment. La finalitat d’aquest projecte és l’estudi de la detecció d’objectes en vídeo. Intentarem crear un algorisme que ens permeti dur a terme aquest estudi i treure’n conclusions. Pretenem fer un algorisme, o un conjunt d’algorismes, en Matlab que sigui capaç de donat qualsevol vídeo, pugui retornar un conjunt de imatges, o un vídeo, amb els diferents objectes de l’escena destacats. Es faran proves en diferents situacions, des de objectes sintètics amb un moviment clarament definit, fins a proves en seqüències reals extretes de diferents pel•lícules. Per últim es pretén comprovar l’eficiència d’aquest. Ja que el projecte s’emmarca en la línia de recerca de robòtica i visió per computador, la tasca principal serà la manipulació d’imatges. Per tant farem servir el Matlab, ja que les imatges no son res més que matrius i aquest programa permet el càlcul vectorial i matricial d’una manera senzilla i realment eficient

Estudi i implementació d’un mètode de reconstrucció 3D basat en SfM i registre de vistes 3D parcials

Aquest projecte es basarà en reconstruir una imatge 3D gran a partir d’una seqüència d’imatges 2D capturades per una càmera. Ens centrem en l’estudi de les bases matemàtiques de la visió per computador així com en diferents mètodes emprats en la reconstrucció 3D d’imatges. Per portar a terme aquest estudi s’utilitza la plataforma de desenvolupament MatLab ja que permet tractar operacions matemàtiques, imatges i matrius de gran tamany amb molta senzillesa, rapidesa i eficiència, per aquesta raó s’usa en moltes recerques sobre aquest tema. El projecte aprofundeix en el tema descrit anteriorment estudiant i implementant un mètode que consisteix en aplicar Structure From Motion (SFM) a pocs frames seguits obtinguts d’una seqüència d’imatges 2D per crear una reconstrucció 3D. Quan s’han creat dues reconstruccions 3D consecutives i fent servir un frame com a mínim en comú entre elles, s’aplica un mètode de registre d’estructures 3D, l’Iterative Closest Point (ICP), per crear una reconstrucció 3D més gran a través d’unir les diferents reconstruccions obtingudes a partir de SfM. El mètode consisteix en anar repetint aquestes operacions fins al final dels frames per poder aconseguir una reconstrucció 3D més gran que les petites imatges que s’aconsegueixen a través de SfM. A la Figura 1 es pot veure un esquema del procés que es segueix. Per avaluar el comportament del mètode, utilitzem un conjunt de seqüències sintètiques i un conjunt de seqüències reals obtingudes a partir d’una càmera. L’objectiu final d’aquest projecte és construir una nova toolbox de MatLab amb tots els mètodes per crear reconstruccions 3D grans per tal que sigui possible tractar amb facilitat aquest problema i seguir-lo desenvolupant en un futur

Three-dimensional patient-specific dosimetry in radioimmunotherapy with 90Y-ibritumomab-tiuxetan.

Aim of the present article was to perform three-dimensional (3D) single photon emission tomography-based dosimetry in radioimmunotherapy (RIT) with (90)Y-ibritumomab-tiuxetan. A custom MATLAB-based code was used to elaborate 3D images and to compare average 3D doses to lesions and to organs at risk (OARs) with those obtained with planar (2D) dosimetry. Our 3D dosimetry procedure was validated through preliminary phantom studies using a body phantom consisting of a lung insert and six spheres with various sizes. In phantom study, the accuracy of dose determination of our imaging protocol decreased when the object volume decreased below 5 mL, approximately. The poorest results were obtained for the 2.58 mL and 1.30 mL spheres where the dose error evaluated on corrected images with regard to the theoretical dose value was -12.97% and -18.69%, respectively. Our 3D dosimetry protocol was subsequently applied on four patients before RIT with (90)Y-ibritumomab-tiuxetan for a total of 5 lesions and 4 OARs (2 livers, 2 spleens). In patient study, without the implementation of volume recovery technique, tumor absorbed doses calculated with the voxel-based approach were systematically lower than those calculated with the planar protocol, with average underestimation of -39% (range from -13.1% to -62.7%). After volume recovery, dose differences reduce significantly, with average deviation of -14.2% (range from -38.7.4% to +3.4%, 1 overestimation, 4 underestimations). Organ dosimetry in one case overestimated, in the other underestimated the dose delivered to liver and spleen. However, both for 2D and 3D approach, absorbed doses to organs per unit administered activity are comparable with most recent literature findings.

Optimização de Células Fotovoltaicas

Perante os cenários do aumento da população mundial, da concentração de CO2, dos custos dos combustíveis, do consumo energético mundial e das alterações climáticas, surgiu a necessidade de encontrar fontes de energias alternativas. Neste contexto, a Energia Solar Fotovoltaica, fruto de investigações e investimentos realizados, teve um grande impacto na última década, registando um aumento significativo quer da produção de painéis fotovoltaicos ou de instalações de sistemas fotovoltaicos no Mundo. A Energia solar Fotovoltaica surge como uma energia alternativa limpa, inesgotável e que contribui para a diminuição do impacto ambiental, mas o elevado custo inicial é ainda um entrave à sua comercialização, sendo por isso importante conseguir uma optimização na produção dos painéis fotovoltaico bem como em instalações a fim de optimizar o seu rendimento. Um dos objectivos deste trabalho foi instalar e monitorizar um sistema fotovoltaico no telhado do laboratório de máquinas eléctricas do ISEL. Foi instalado um sistema com uma potência de 990 Watts. A monitorização dos módulos durante alguns períodos de 2011 e 2012 demonstraram um bom desempenho do sistema fotovoltaico instalado comparativamente aos valores estimados. Outro objectivo deste trabalho foi estudar a influência da temperatura no rendimento das células fotovoltaicas. A primeira fase deste estudo foi desenvolver um modelo matemático de uma célula fotovoltaica em Simulink/Matlab. As curvas obtidas da simulação numérica do modelo matemático permitiram observar e demonstrar a influência do aumento da temperatura das células fotovoltaicas na sua potência e rendimento. A segunda fase deste estudo tinha como objectivo comprovar experimentalmente o efeito da temperatura e analisar possíveis meios que permitissem refrigerar as células fotovoltaicas. Através de uma montagem experimental específica as células fotovoltaicas foram testadas num ambiente controlado. Os valores obtidos permitiram observar uma diminuição de cerca de 36% da temperatura das células utilizando refrigeração e consequente aumento do rendimento.

Otimização Não Linear de Mínimos Quadrados

O problema de otimização de mínimos quadrados e apresentado como uma classe importante de problemas de minimização sem restrições. A importância dessa classe de problemas deriva das bem conhecidas aplicações a estimação de parâmetros no contexto das analises de regressão e de resolução de sistemas de equações não lineares. Apresenta-se uma revisão dos métodos de otimização de mínimos quadrados lineares e de algumas técnicas conhecidas de linearização. Faz-se um estudo dos principais métodos de gradiente usados para problemas não lineares gerais: Métodos de Newton e suas modificações incluindo os métodos Quasi-Newton mais usados (DFP e BFGS). Introduzem-se depois métodos específicos de gradiente para problemas de mínimos quadrados: Gauss-Newton e Levenberg-Larquardt. Apresenta-se uma variedade de exemplos selecionados na literatura para testar os diferentes métodos usando rotinas MATLAB. Faz-se uma an alise comparativa dos algoritmos baseados nesses ensaios computacionais que exibem as vantagens e desvantagens dos diferentes métodos.

Desempenho Hidráulico de Exutores Submarinos

Nos centros urbanos próximos das orlas marítimas, tem-se optado pela descarga de águas residuais para o mar, aproveitando-se da elevada capacidade de depuração do oceano. Neste trabalho coloca-se o seguinte desafio: será possível um programa usado para a modelação de rede de distribuição de água ser capaz de modelar um exutor submarino? Para responder a este desafio, recorreu-se a um conjunto de artifícios, tendo por base conhecimentos da hidráulica. Foi construído um modelo unidimensional usando o EPANET cujas fronteiras são: a montante, a Câmara de Carga, a jusante os difusores. Recorreu-se aos componentes do EPANET, tais como: válvulas TCV (Throttle Control Valve) e reservatórios de nível fixo para simular a descarga e o nível de maré respetivamente. A modelação hidráulica de um exutor submarino tem implicações no impacto da descarga pelo que torna-se necessário o estudo da diluição inicial conseguida numa descarga. No presente trabalho desenvolveu-se um programa EPANET / MatLab com o objetivo de facilitar a aplicação desses artifícios e de estimar a diluição inicial. Aplicou-se o programa desenvolvido num caso de estudo e utilizou-se a versão beta do CORHYD, usado em obras de engenharia com alguma importância, para validar os resultados obtido através do programa EPANET /MatLab.

Análise do Desempenho de Amplificadores Ópticos Semiconductores Reflectivos em Redes Ópticas de Acesso da Próxima Geração

O constante crescimento de aplicações Internet têm vindo a contribuir para a necessidade do aumento de largura de banda na rede de acesso, o que fomenta a procura de novas soluções para redes de acesso de fibra óptica. Nesse contexto, soluções conhecidas como transmissores independentes do comprimento de onda (WIT wavelengthindependent transmitters) isto é, soluções para emissores que possam operar numa ampla gama de comprimentos de onda, poderão ser muito úteis. Sendo uma das mais promissoras o amplificador óptico semiconductor reflectivo (RSOA). Nesta dissertação foi implementado um modelo simplificado do RSOA que permite eficientemente determinar as suas características estáticas e dinâmicas. Foi também implementado um modelo de simulação, em Matlab, que permita estudar o desempenho de um sistema de comunicação por fibra óptica que utiliza RSOA e técnicas de modulação digital por desvio de fase com M estados (M-PSK M-ary Phase Shift Keying). Com o modelo de simulação foi avaliada qual a taxa de transmissão máxima atingível. A avaliação do desempenho foi feita através do cálculo da magnitude do vector erro (EVM Error Vector Magnitude), considerando os vários blocos implementados. Foi determinado o ponto óptimo de funcionamento do RSOA, sendo também avaliado do desempenho de sistema de transmissão, usando diferentes parâmetros do sinal M-PSK (M=4,16,64) com taxa de transmissão de 1Gbits/s e 10 Gbits/s considerando vários comprimentos de fibra óptica.

Solving heterogeneous-agent models by projection and perturbation

The paper proposes a numerical solution method for general equilibrium models with a continuum of heterogeneous agents, which combines elements of projection and of perturbation methods. The basic idea is to solve first for the stationary solutionof the model, without aggregate shocks but with fully specified idiosyncratic shocks. Afterwards one computes a first-order perturbation of the solution in the aggregate shocks. This approach allows to include a high-dimensional representation of the cross-sectional distribution in the state vector. The method is applied to a model of household saving with uninsurable income risk and liquidity constraints. The model includes not only productivity shocks, but also shocks to redistributive taxation, which cause substantial short-run variation in the cross-sectional distribution of wealth. If those shocks are operative, it is shown that a solution method based on very few statistics of the distribution is not suitable, while the proposed method can solve the model with high accuracy, at least for the case of small aggregate shocks. Techniques are discussed to reduce the dimension of the state space such that higher order perturbations are feasible.Matlab programs to solve the model can be downloaded.

Análise do Sistema de Tratamento de Efluentes Gasosos de uma Central de Valorização Energética de Resíduos Sólidos Urbanos

O sistema de tratamento dos gases formados durante a combustão de Resíduos Sólidos Urbanos, (RSU), é um dos componentes mais importantes de uma instalação de incineração. O seu controlo e optimização asseguram a redução dos gases nocivos formados para valores que não representam quaisquer perigos para a saúde pública e ambiental, o que leva a incineração a ser considerada uma das técnicas mais seguras e eficientes para o tratamento dos RSU. A combinação desta técnica com produção de energia eléctrica, por recuperação de calor da combustão dos resíduos, tem contribuído gradualmente para o desenvolvimento socioeconómico. Este trabalho surgiu da necessidade de optimização do sistema de controlo dos gases ácidos, (HCl, NOx, SO2 e HF), da Central de Tratamento de Resíduos Sólidos Urbanos, (CTRSU), da Empresa Valorsul e dos respectivos consumos dos reagentes utilizados no tratamento dos referidos gases, nomeadamente o leite de cal e a amónia. No entanto, procedeu-se apenas com uma análise detalhada do referido sistema de controlo, de modo a identificar as principais dificuldades e possíveis soluções teóricas para as mesmas. Os conhecimentos adquiridos poderão contribuir futuramente para a optimização do sistema de controlo em questão, de modo a garantir a dosagem óptima dos referidos reagentes, diminuindo assim o custo dos mesmos. O presente trabalho também tem como objectivos, a compreensão da tecnologia de valorização energética de resíduos no contexto do tratamento de RSU em Portugal e a análise do processo de formação dos gases durante a combustão dos RSU, principalmente a dos gases ácidos. Em Portugal o tratamento dos RSU é sobretudo por eliminação, deposição directa em aterro sanitário, tendo ocupado em 2012, uma percentagem de 54% do total dos RSU produzidos, no entanto a valorização tem vindo a desenvolver-se e actualmente cerca de 21% dos resíduos produzidos é valorizado energeticamente. Do mesmo modo, a Reciclagem e a valorização orgânica têm progredido nos últimos anos, embora lentamente. A formação dos gases ácidos é influenciada sobretudo pela composição física dos RSU, que contêm na sua matriz materiais constituídos por átomos de azoto, enxofre, cloro e fluor e por certas condições operatórias de incineração, nomeadamente temperatura elevadas, humidade e excesso de ar. A composição física dos resíduos é um parâmetro não controlável pelo operador do sistema de tratamento dos gases, pelo que as condições de incineração devem ser adequadamente controladas, no âmbito de prevenir ou reduzir a formação destes gases. A CTRSU da Valorsul apresenta, em geral, um sistema de tratamento dos gases eficiente, garantido assim que os valores limites de emissão dos gases se encontrem em conformidade com a lei. No entanto, para cumprir este requisito é necessário um consumo excessivo da amónia e leite de cal utilizados no tratamento destes gases. O elevado consumo destes reagentes deve-se sobretudo à dificuldade de optimização da dosagem dos mesmos, como VI consequência de um atraso verificado no sistema de controlo dos gases implementado na Central, “Sample and Hold PID”. Este atraso é de aproximadamente dois minutos e é causado pelo tempo em que os gases levam a percorrer o ponto onde é feita a injecção dos reagentes até ao ponto de extracção dos gases na chaminé, mais o tempo de atraso imposto pelo sistema de monitorização em contínuo das emissões. Através do software simulink do MATLAB realizou-se uma simples simulação do sistema de controlo de emissão dos gases NOx e caudal de amónia, em malha aberta com controlador proporcional, com o objectivo de reproduzir os resultados obtidos através de um controlo real em malha aberta realizada na CTRSU. Verificou-se que estes estão de acordo com o esperado, ou seja as saídas dos controladores do caudal de amónia e concentração do NOx, correspondem às variações efectuadas na válvula de regulação de amónia. Como possíveis soluções teóricas para a redução do atraso do sistema, sugeriu-se a introdução de um compensador de atraso baseado na técnica de Predição de Smith, no sistema de controlo já implementado, ou a implementação de um novo sistema de controlo baseado em sistemas de inteligências artificiais, nomeadamente as Redes Neuronais.

Calculation of radio electrical coverage in Medium-Wave Frequencies

The aim of this project is to accomplish an application software based on Matlab to calculate the radioelectrical coverage by surface wave of broadcast radiostations in the band of Medium Wave (WM) all around the world. Also, given the location of a transmitting and a receiving station, the software should be able to calculate the electric field that the receiver should receive at that specific site. In case of several transmitters, the program should search for the existence of Inter-Symbol Interference, and calculate the field strenght accordingly. The application should ask for the configuration parameters of the transmitter radiostation within a Graphical User Interface (GUI), and bring back the resulting coverage above a map of the area under study. For the development of this project, it has been used several conductivity databases of different countries, and a high-resolution elevation database (GLOBE). Also, to calculate the field strenght due to groundwave propagation, it has been used ITU GRWAVE program, which must be integrated into a Matlab interface to be used by the application developed.

Programació d'una aplicació gràfica per l'anàlisi de senyals usant LDB (Local Discriminant Bases) i MLDB (Modified LDB)

Tot seguit presentem un entorn per analitzar senyals de tot tipus amb LDB (Local Discriminant Bases) i MLDB (Modified Local Discriminant Bases). Aquest entorn utilitza funcions desenvolupades en el marc d’una tesi en fase de desenvolupament. Per entendre part d’aquestes funcions es requereix un nivell de coneixement avançat de processament de senyals. S’han extret dels treballs realitzats per Naoki Saito [3], que s’han agafat com a punt de partida per la realització de l’algorisme de la tesi doctoral no finalitzada de Jose Antonio Soria. Aquesta interfície desenvolupada accepta la incorporació de nous paquets i funcions. Hem deixat un menú preparat per integrar Sinus IV packet transform i Cosine IV packet transform, tot i que també podem incorporar-n’hi altres. L’aplicació consta de dues interfícies, un Assistent i una interfície principal. També hem creat una finestra per importar i exportar les variables desitjades a diferents entorns. Per fer aquesta aplicació s’han programat tots els elements de les finestres, en lloc d’utilitzar el GUIDE (Graphical User Interface Development Enviroment) de MATLAB, per tal que sigui compatible entre les diferents versions d’aquest programa. En total hem fet 73 funcions en la interfície principal (d’aquestes, 10 pertanyen a la finestra d’importar i exportar) i 23 en la de l’Assistent. En aquest treball només explicarem 6 funcions i les 3 de creació d’aquestes interfícies per no fer-lo excessivament extens. Les funcions que explicarem són les més importants, ja sigui perquè s’utilitzen sovint, perquè, segons la complexitat McCabe, són les més complicades o perquè són necessàries pel processament del senyal. Passem cada entrada de dades per part de l’usuari per funcions que ens detectaran errors en aquesta entrada, com eliminació de zeros o de caràcters que no siguin números, com comprovar que són enters o que estan dins dels límits màxims i mínims que li pertoquen.

Parametrització de contorns oberts

Un dels principals problemes quan es realitza un anàlisi de contorns és la gran quantitat de dades implicades en la descripció de la figura. Per resoldre aquesta problemàtica, s’aplica la parametrització que consisteix en obtenir d’un contorn unes dades representatives amb els mínims coeficients possibles, a partir dels quals es podrà reconstruir de nou sense pèrdues molt evidents d’informació. En figures de contorns tancats, la parametrització més estudiada és l’aplicació de la transformada discreta de Fourier (DFT). Aquesta s’aplica a la seqüència de valors que descriu el comportament de les coordenades x i y al llarg de tots els punts que formen el traç. A diferència, en els contorns oberts no es pot aplicar directament la DFT ja que per fer-ho es necessita que el valor de x i de y siguin iguals tan en el primer punt del contorn com en l’últim. Això és degut al fet que la DFT representa sense error senyals periòdics. Si els senyals no acaben en el mateix punt, representa que hi ha una discontinuïtat i apareixen oscil·lacions a la reconstrucció. L’objectiu d’aquest treball és parametritzar contorns oberts amb la mateixa eficiència que s’obté en la parametrització de contorns tancats. Per dur-ho a terme, s’ha dissenyat un programa que permet aplicar la DFT en contorns oberts mitjançant la modificació de les seqüencies de x i y. A més a més, també utilitzant el programari Matlab s’han desenvolupat altres aplicacions que han permès veure diferents aspectes sobre la parametrització i com es comporten els Descriptors El·líptics de Fourier (EFD). Els resultats obtinguts han demostrat que l’aplicació dissenyada permet la parametrització de contorns oberts amb compressions òptimes, fet que facilitarà l’anàlisi quantitatiu de formes en camps com l’ecologia, medicina, geografia, entre d’altres.