A inferência de redes gênicas (GN) a partir de dados de expressão gênica temporal é um problema crucial e desafiador em Biologia Sistêmica. Os conjuntos de dados de expressão geralmente consistem em dezenas de amostras temporais e as redes consistem em milhares de genes, tornando inúmeros métodos de inferência inviáveis na prática. Para melhorar a escalabilidade dos métodos de inferência de GNs, esta tese propõe um arcabouço chamado GeNICE, baseado no modelo de redes gênicas probabilísticas. A principal novidade é a introdução de um procedimento de agrupamento de genes, com perfis de expressão relacionados, para fornecer uma solução aproximada com complexidade computacional reduzida. Os agrupamentos definidos são usados para reduzir a dimensionalidade permitindo uma busca exaustiva mais eficiente pelos melhores subconjuntos de genes preditores para cada gene alvo de acordo com funções critério multivariadas. GeNICE reduz consideravelmente o espaço de busca porque os candidatos a preditores ficam restritos a um gene representante por agrupamento. No final, uma análise multivariada é realizada para cada subconjunto preditor definido, visando recuperar subconjuntos mínimos para simplificar a rede gênica inferida. Em experimentos com conjuntos de dados sintéticos, GeNICE obteve uma redução substancial de tempo quando comparado a uma solução anterior sem a etapa de agrupamento, preservando a precisão da predição de expressão gênica mesmo quando o número de agrupamentos é pequeno (cerca de cinquenta) e o número de genes é grande (ordem de milhares). Para um conjunto de dados reais de microarrays de Plasmodium falciparum, a precisão da predição alcançada pelo GeNICE foi de aproximadamente 97% em média. As redes inferidas para os genes alvos da glicólise e do apicoplasto refletem propriedades topológicas de redes complexas do tipo \"mundo pequeno\" e \"livre de escala\", para os quais grande parte das conexões são estabelecidas entre os genes de um mesmo módulo e algumas poucas conexões fazem o papel de estabelecer uma ponte entre os módulos (redes mundo pequeno), e o grau de distribuição das conexões entre os genes segue uma lei de potência, na qual a maioria dos genes têm poucas conexões e poucos genes (hubs) apresentam um elevado número de conexões (redes livres de escala), como esperado. / Gene network (GN) inference from temporal gene expression data is a crucial and challenging problem in Systems Biology. Expression datasets usually consist of dozens of temporal samples, while networks consist of thousands of genes, thus rendering many inference methods unfeasible in practice. To improve the scalability of GN inference methods, this work proposes a framework called GeNICE, based on Probabilistic Gene Networks; the main novelty is the introduction of a clustering procedure to group genes with related expression profiles, to provide an approximate solution with reduced computational complexity. The defined clusters were used to perform an exhaustive search to retrieve the best predictor gene subsets for each target gene, according to multivariate criterion functions. GeNICE greatly reduces the search space because predictor candidates are restricted to one representative gene per cluster. Finally, a multivariate analysis is performed for each defined predictor subset to retrieve minimal subsets and to simplify the network. In experiments with in silico generated datasets, GeNICE achieved substantial computational time reduction when compared to an existing solution without the clustering step, while preserving the gene expression prediction accuracy even when the number of clusters is small (about fifty) relative to the number of genes (order of thousands). For a Plasmodium falciparum microarray dataset, the prediction accuracy achieved by GeNICE was roughly 97% on average. The inferred networks for the apicoplast and glycolytic target genes reflects the topological properties of \"small-world\"and \"scale-free\"complex network models in which a large part of the connections is established between genes of the same functional module (smallworld networks) and the degree distribution of the connections between genes tends to form a power law, in which most genes present few connections and few genes (hubs) present a large number of connections (scale-free networks), as expected.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-05-15T17:29:47Z No. of bitstreams: 1 bernardojacquesdelgadocosta.pdf: 7161081 bytes, checksum: 43a871813092712bd1059913da9e445f (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-05-17T16:01:38Z (GMT) No. of bitstreams: 1 bernardojacquesdelgadocosta.pdf: 7161081 bytes, checksum: 43a871813092712bd1059913da9e445f (MD5) / Made available in DSpace on 2017-05-17T16:01:38Z (GMT). No. of bitstreams: 1 bernardojacquesdelgadocosta.pdf: 7161081 bytes, checksum: 43a871813092712bd1059913da9e445f (MD5) Previous issue date: 2017-02-23 / Este trabalho propõe um método capaz de promover, em tempo real e de forma automática, a restauração de um sistema elétrico de distribuição radial frente a uma contingência em qualquer ponto, utilizando somente dados facilmente obtidos pelos equipamentos de proteção e controle disponíveis, tais como: valores de tensão, corrente, fator de potência e estado do equipamento. Devido à característica de operação autônoma do método, são considerados apenas os equipamentos telecontrolados instalados na rede. Muitos métodos têm sido propostos para o problema de restabelecimento de energia em sistemas de distribuição. Entretanto, em sua maioria, são dependentes de dados da rede que não podem ser facilmente obtidos, tais como os valores de impedância das linhas, demanda das cargas, entre outros dados. Além disso, normalmente são propostos planos de restabelecimento que deverão ser analisados pelo operador do Centro de Operações da Distribuição (COD) para que ele, posteriormente, execute as ações necessárias para o restabelecimento do sistema. O método proposto pode ser resumido como segue. Após a atuação do dispositivo de proteção frente a uma falta, o algoritmo entra em operação e identifica a área afetada, promovendo o isolamento dos trechos com defeito. Feito isso, é realizada uma análise das possibilidades de restauração do sistema, que se baseia na seleção de trechos pertencentes a alimentadores adjacentes à área afetada, que atendam algumas restrições. Após a seleção, é realizado um teste de todas as combinações possíveis de serem realizadas com esses trechos e a melhor solução é executada, podendo essa ser baseada no maior número de cargas restabelecidas ou na restauração do maior número de cargas prioritárias. / This work proposes a method capable of promoting, in real time and automatically, the restoration of a radial distribution system after a contingency at any point, using only data easily obtained by the available protection and control equipments, such as: voltage, current, power factor and equipment status. Due to the automatic operation feature of the method, only the equipments with remote control installed in the network are considered. Many methods have been proposed for the problem of distribution systems restoration. However, most rely on network data that cannot be easily obtained, such as line impedance, load demand, and other data. In addition, the proposed restoration plans should normally be reviewed by the operator of Distribution Operations Center (DOC), so that he can then take the necessary steps to restore the system. The proposed method can be summarized as follows. After the protection device acts against a fault, the algorithm starts and identifies the affected area, promoting the isolation of the faulty sections. Once this is done, an analysis is made of the possibilities of system restoration, which is based on the selection of sections belonging to adjacent feeders of the affected area, which meet some restrictions. After the selection, all possible combinations to be performed with these sections are tested and the best solution is performed, which can be based on the largest number of loads restored or the restoration of the largest number of priority loads.

There are some systems (or parts of systems) that are very critical and need especially good test suites to test them. For these critical systems, exhaustive testing may be a good way to test them. Thus, we have implemented two versions of bounded exhaustive search (BES) algorithms in a model-based testing tool called, Modbat. One of the BES versions (BESnL) visits each self-loop in the model only once. The other version (BESL) has no constraint or limitation on the number of time it visits each self-loop. We have then compared the two BES algorithms with each other and with an already implemented algorithm in Modbat called random search (RS). We have run the three mentioned algorithms (BESL, BESnL and RS) on five different models and compared their performance on these models in terms of time, coverage and finding faults. We have found that BESnL is faster than BESL and it can miss some faults that BESL can find. However, BESnL can find errors faster than BESL. BESL has sometimes better performance than BESnL in terms of branch coverage. In terms of other coverage criteria (like state coverage, transition coverage and instruction coverage), both BESL and BESnL has very similar performance. We have also found that running the RS algorithm is, in general, faster than both BES algorithms at generating tests (given the same total number of tests generated) if the model has a clear end state. RS may also be faster at finding faults than the BES algorithms. However, The BES algorithms and the RS algorithm have similar behaviours regarding coverage. Nevertheless, RS can sometimes reach higher coverage faster than the BES algorithms and with a smaller number of tests. / Det finns vissa system (eller delar av system) som är mycket essentiella och som behöver särskilt bra testsviter för att testa dem. För dessa essentiella system kan uttömmande tester vara ett bra sätt att testa dem. Således har vi implementerat två versioner av begränsad uttömmande sökning eller på engelska ”bounded exhuastive search” (BES) algoritmer i ett modellbaserat testverktyg kallat Modbat. En av BES-versionerna (BESnL) besöker varje självslinga i modellen bara en gång. Den andra versionen (BESL) har ingen begränsning av hur många gånger den besöker varje självslinga. Vi har sedan jämfört de två BES-algoritmerna med varandra och med en redan implementerad algoritm i Modbat som kallas slumpmässig sökning eller på engelska ”random search” (RS). Vi har kört de tre nämnda algoritmerna (BESL, BESnL och RS) på fem olika modeller och jämfört deras prestanda på dessa modeller när det gäller tid, täckning (coverage) och att hitta fel. Vi har funnit att BESnL är snabbare än BESL och det kan missa några fel som BESL kan hitta, men BESnL kan hitta fel snabbare än BESL. BESL har ibland bättre prestanda än BESnL när det gäller filialtäckning (branch-coverage). När det gäller andra täckningskriterier (som statlig täckning, övergångstäckning (tranintion-coverage) och instruktionstäckning) har både BESL och BESnL mycket liknande resultat. Vi har också funnit att körning av RS-algoritmen i allmänhet är snabbare än båda BES- algoritmerna vid generering av tester (givet samma totala antal genererade tester) om modellen har ett klart slutläge (end-state). RS kan också vara snabbare att hitta fel än BES-algoritmerna. BES-algoritmerna och RS-algoritmen har dock liknande beteenden när det gäller täckning. RS kan ibland nå högre täckning snabbare än BES-algoritmerna och med ett mindre antal tester.

O presente trabalho apresenta a identificação de parâmetros de motor de indução linear tubular para a extração de petróleo, denominado de MATÆOS -- Motor Assíncrono Tubular para Aplicação na Extração de Óleo Subterrâneo. O primeiro passo para a identificação de parâmetros do motor foi a implementação da bancada de testes com o seu sistema de controle e supervisório. O modelo adotado corresponde à situação de operação em regime estacionário e consiste no circuito contendo resistores e indutores tradicionalmente utilizado para representar motores de indução. A identificação dos parâmetros foi realizada por meio de um procedimento de otimização tomando por base o balanço de potência do motor. Para isso, em função dos parâmetros do modelo, definiu-se e minimizou-se uma medida de erro entre as potências medidas e as calculadas pelo modelo para um conjunto pré-fixado de escorregamentos. A solução desse problema de otimização foi obtida por meio de uma busca exaustiva, uma vez que o número de variáveis independentes do problema é pequeno e uma região contendo o ponto de mínimo é conhecida com razoável confiança. Os resultados obtidos indicam que os valores dos parâmetros nominais a quente apresentam algumas diferenças em relação àqueles produzidos pelo método de identificação. Essas diferenças se devem possivelmente a variações nas características térmicas do motor e a erros introduzidos pelos instrumentos de medida. O modelo identificado apresentou uma menor disponibilidade de força e um pior rendimento quando comparado com os valores de projeto do modelo nominal a quente. Apesar disto, verificou-se uma boa proximidade entre a curva característica do modelo e as obtidas experimentalmente. / This work presents the identification of parameters of the tubular linear induction motor for the extraction of oil, called MATÆOS - which is the Portuguese acronym for Asynchronous Tubular Motor for Application to Subterranean Oil Extraction. The first step to the identification of the motor parameters was the implementation of the supervisory and control system for the test bench. The model, which is the usual one containing resistors and inductors, represents the motor operating in steady-state. The parameter identification problem was written as an optimization problem based on the motor power balance. An error between the measured power and the one given by the model was defined and minimized with respect to the model parameters in correspondence to a given set of slipping frequencies. The numerical solution to this problem was obtained by an exhaustive search algorithm. This choice was made because the number of independent variables is small and a region containing the solution is known with a significant degree of confidence. The results obtained show that the nominal values of the model parameters differ slightly from those produced by the identification method. Possibly, such differences are due to the variations in the model caused by temperature changes and to the measurement errors. The identified model presented both a lower available force and a worst efficiency when compared to the nominal design values for the warm motor. Despite this the characteristic curve of the model was quite close to the nominal design one.

Road network of a region is of a paramount importance in the overall development. Management of road traffic is a key factor for the city authority and reducing the road traffic congestion is a significant challenge in this perspective. In this thesis, a Wireless Sensor Network (WSN) based road-traffic monitoring scheme with dynamic mathematical traffic model is presented that will not necessarily include all adjacent intersections of a block; rather the important major intersections of a city. The objective of this scheme is to reduce the congestion by re-routing the vehicles to better performing road-segments by informing the down-stream drivers through broadcasting the congestion information in a dedicated radio channel. The dynamic model can provide with the instantaneous status of the traffic of the road-network. The scheme is a WSN based multi-hop relay network with hierarchical architecture and composed of ordinary nodes, Cluster-Head nodes, Base Stations, Gateway nodes and Monitoring and Control Centers (MCC) etc. Through collecting the traffic information, MCC will check the congestion status and in defining the congestion, threshold factors have been used in this model. For the congested situation of a road-segment, a cost function has been defined as a performance indicator and estimated using the weight factors (importance) of these selected intersections. This thesis considered a traffic network with twelve major intersections of a city with four major directions. Traffic arrivals in these intersections are assumed to follow Poisson distribution. Model was simulated in Matlab with traffic generated through Poisson Random Number Generator and cost function was estimated for the congestion status of the road-segments over a simulation period of 1440 minutes starting from midnight. For optimization purpose we adopted two different approaches; in the first approach, performance of the scheme was evaluated for all threshold factor values iteratively one at a time, applying a threshold factor value to define threshold capacities of all the road segments; traffic was generated and relative cost has been estimated following the model specifications with the purpose of congestion avoidance. In the second approach, different values of threshold factor have been used for different road segments for determining the optimum set-up, and exhaustive search technique has been applied with a smaller configuration in order to keep computations reachable. Simulation results show the capacity of this scheme to improve the traffic performance by reducing the congestion level with low congestion costs.

O presente trabalho apresenta a identificação de parâmetros de motor de indução linear tubular para a extração de petróleo, denominado de MATÆOS -- Motor Assíncrono Tubular para Aplicação na Extração de Óleo Subterrâneo. O primeiro passo para a identificação de parâmetros do motor foi a implementação da bancada de testes com o seu sistema de controle e supervisório. O modelo adotado corresponde à situação de operação em regime estacionário e consiste no circuito contendo resistores e indutores tradicionalmente utilizado para representar motores de indução. A identificação dos parâmetros foi realizada por meio de um procedimento de otimização tomando por base o balanço de potência do motor. Para isso, em função dos parâmetros do modelo, definiu-se e minimizou-se uma medida de erro entre as potências medidas e as calculadas pelo modelo para um conjunto pré-fixado de escorregamentos. A solução desse problema de otimização foi obtida por meio de uma busca exaustiva, uma vez que o número de variáveis independentes do problema é pequeno e uma região contendo o ponto de mínimo é conhecida com razoável confiança. Os resultados obtidos indicam que os valores dos parâmetros nominais a quente apresentam algumas diferenças em relação àqueles produzidos pelo método de identificação. Essas diferenças se devem possivelmente a variações nas características térmicas do motor e a erros introduzidos pelos instrumentos de medida. O modelo identificado apresentou uma menor disponibilidade de força e um pior rendimento quando comparado com os valores de projeto do modelo nominal a quente. Apesar disto, verificou-se uma boa proximidade entre a curva característica do modelo e as obtidas experimentalmente. / This work presents the identification of parameters of the tubular linear induction motor for the extraction of oil, called MATÆOS - which is the Portuguese acronym for Asynchronous Tubular Motor for Application to Subterranean Oil Extraction. The first step to the identification of the motor parameters was the implementation of the supervisory and control system for the test bench. The model, which is the usual one containing resistors and inductors, represents the motor operating in steady-state. The parameter identification problem was written as an optimization problem based on the motor power balance. An error between the measured power and the one given by the model was defined and minimized with respect to the model parameters in correspondence to a given set of slipping frequencies. The numerical solution to this problem was obtained by an exhaustive search algorithm. This choice was made because the number of independent variables is small and a region containing the solution is known with a significant degree of confidence. The results obtained show that the nominal values of the model parameters differ slightly from those produced by the identification method. Possibly, such differences are due to the variations in the model caused by temperature changes and to the measurement errors. The identified model presented both a lower available force and a worst efficiency when compared to the nominal design values for the warm motor. Despite this the characteristic curve of the model was quite close to the nominal design one.

Road network of a region is of a paramount importance in the overall development. Management of road traffic is a key factor for the city authority and reducing the road traffic congestion is a significant challenge in this perspective. In this thesis, a Wireless Sensor Network (WSN) based road-traffic monitoring scheme with dynamic mathematical traffic model is presented that will not necessarily include all adjacent intersections of a block; rather the important major intersections of a city. The objective of this scheme is to reduce the congestion by re-routing the vehicles to better performing road-segments by informing the down-stream drivers through broadcasting the congestion information in a dedicated radio channel. The dynamic model can provide with the instantaneous status of the traffic of the road-network. The scheme is a WSN based multi-hop relay network with hierarchical architecture and composed of ordinary nodes, Cluster-Head nodes, Base Stations, Gateway nodes and Monitoring and Control Centers (MCC) etc. Through collecting the traffic information, MCC will check the congestion status and in defining the congestion, threshold factors have been used in this model. For the congested situation of a road-segment, a cost function has been defined as a performance indicator and estimated using the weight factors (importance) of these selected intersections. This thesis considered a traffic network with twelve major intersections of a city with four major directions. Traffic arrivals in these intersections are assumed to follow Poisson distribution. Model was simulated in Matlab with traffic generated through Poisson Random Number Generator and cost function was estimated for the congestion status of the road-segments over a simulation period of 1440 minutes starting from midnight. For optimization purpose we adopted two different approaches; in the first approach, performance of the scheme was evaluated for all threshold factor values iteratively one at a time, applying a threshold factor value to define threshold capacities of all the road segments; traffic was generated and relative cost has been estimated following the model specifications with the purpose of congestion avoidance. In the second approach, different values of threshold factor have been used for different road segments for determining the optimum set-up, and exhaustive search technique has been applied with a smaller configuration in order to keep computations reachable. Simulation results show the capacity of this scheme to improve the traffic performance by reducing the congestion level with low congestion costs.

This master´s thesis deals with the registration of ultrasound sequences using evolutionary algorithms. The theoretical part of the thesis describes the process of image registration and its optimalization using genetic and metaheuristic algorithms. The thesis also presents problems that may occur during the registration of ultrasonographic images and various approaches to their registration. In the practical part of the work, several optimization methods for the registration of a number of sequences were implemented and compared.

This diploma thesis deals with evolutionary algorithms used for travelling salesman problem (TSP). In the first section, there are theoretical foundations of a graph theory and computational complexity theory. Next section contains a description of chosen optimization algorithms. The aim of the diploma thesis is to implement an application that solve TSP using evolutionary algorithms.

No contexto da operação de sistemas de distribuição, um dos problemas com os quais os operadores lidam frequentemente é o de restabelecimento de energia. Este problema surge na ocorrência de uma falta permanente e pode ser tratado por meio de manobras em chaves presentes na rede primária. Uma vez que tais redes operam com topologia radial, a ocorrência de uma falta pode resultar no desligamento de consumidores saudáveis. Desta maneira, o problema consiste em definir, num curto intervalo de tempo, um número mínimo de chaves que devem ser operadas a fim de isolar a falta e restaurar o máximo de consumidores saudáveis desligados. Os esforços para a obtenção de ferramentas computacionais para fornecimento de soluções para o problema de restabelecimento têm sido intensificados nos últimos anos. Isto ocorre, em especial, devido aos enormes prejuízos causados pela falta de energia às companhias de eletricidade e a toda a sociedade. Neste sentido, o objetivo deste trabalho é a obtenção de um método para auxiliar o trabalho dos operadores através do fornecimento de planos adequados de restabelecimento em curtos intervalos de tempo. Os diferenciais deste método proposto são a sua capacidade de: lidar, em especial, com redes reais de grande porte com reduzido esforço computacional; considerar a existência de vários níveis de prioridade de atendimento entre os consumidores (note, por exemplo, que um hospital ou um centro de segurança pública devem ter maior prioridade de atendimento que um grande supermercado ou unidades residenciais) e priorizar o atendimento deles de acordo a sua prioridade; fornecer uma sequência por meio da qual as chaves possam ser operadas a fim de isolar os setores em falta e reconectar o maior número de consumidores saudáveis desligados executando-se o mínimo de manobras em chaves e priorizando os consumidores com maior prioridade; ser capaz de selecionar cargas menos prioritárias para permaneceram desligadas nas situações em que não é possível obter uma solução que restaure todas as cargas saudáveis fora de serviço; e, adicionalmente, priorizar a operação de chaves controladas remotamente, que, diferentemente das chaves controladas manualmente, podem ser operadas com menores custos e de maneira mais rápida. O método proposto consiste, de maneira sintética, na união de uma busca exaustiva aplicada localmente a um novo algoritmo evolutivo multi-objetivo em tabelas de subpopulação que faz uso de uma estrutura de dados eficiente denominada Representação Nó-Profundidade. Para avaliar a performance relativa do método proposto, simulações foram realizadas num sistema de distribuição de pequeno porte e os resultados foram comparados com os obtidos por um método de Programação Matemática. Na sequência, novos experimentos foram realizadas em diversos casos de falta na rede de distribuição da cidade de Londrina-PR e cidades adjacentes. As soluções fornecidas mostraram-se adequadas ao tratamento dos casos de falta, assim como as sequências de chaveamento associadas a elas, as quais foram capazes de priorizar o restabelecimento dos consumidores prioritários seguindo seus níveis de prioridade. Adicionalmente, estudos avaliaram a variação do tempo de processamento computacional do método proposto com a dimensão das redes de distribuições e também com o número de gerações realizadas pelo algoritmo evolutivo multi-objetivo proposto e os resultados mostraram-se satisfatórios às necessidades do problema Portanto, pode-se comprovar que o método proposto atingiu os objetivos especificados, em especial, o tratamento de aspectos práticos do problema. Além do próprio método proposto, algumas contribuições desta pesquisa são a proposição um novo algoritmo evolutivo multiobjetivo em tabelas de subpopulação e de um novo operador para manipulação de florestas de grafo armazenadas pela Representação Nó-Profundidade e voltado ao problema de restabelecimento. / In the context of distribution systems operation, service restoration is one of the problems with which operators constantly deal. It arises when a permanent fault occurs and is treated trough operations in switches at primary grid. Since distribution systems are usually radial, fault occurrence turns-off healthy customers. Thereby, the service restoration problem consists in defining, in a short processing time, the minimum amount of switches that must be operated for the isolation of the fault and reconnection of the maximum amount of healthy out-of-service customers. The efforts of developing computational tools for getting solution to this problems has increased in the last years. It is, in special, due to enormous losses caused to the utilities and to the whole society. In this sense, the main objective of this research is getting a method able to help the distribution system operator\'s work through providing service restoration plans quickly. The differentials of this research are its ability to: deal, in special, with large scale grids whit a reduced computational effort; consider costumers of several priority levels (note, for instance, a hospital has a higher supply priority in relation to a big supermarket) and prioritize the higher priority customers; provide a switching sequence able to isolate and reconnect the maximum amount of healthy out-of-service customer by the minimum amount of switching actions; select lower priority customers to keep out-of-service in order to reconnect higher priority customers when a it is not possible to restore all customers; and, additionally, prioritize switching operation in remotely controlled switches, whose operation is faster and cheapest than the operation of manually controlled switches. The proposed method mixes a local exhaustive search and a new multi-objective evolutionary algorithm in subpopulation tables that uses a data structure named Node-Depth Encoding. For evaluating the relative performance of proposed method, simulations were performed in small distribution systems and the performance was compared with the performance a Mathematical Programing method from literature. New experiments were performed a Mathematical Programing method from literature. New experiments were performed in several fault situations in the real and large-scale distribution system of Londrina-PR and adjacent cities. The solutions provided were appropriated to the treatment of such contingency situations. The same occurs with the switching sequences provided, which were able to prioritize the restoration of higher priority customers. Additional studies evaluated the variation of the running time with the size of grids and with the values adopted for the maximum number of generations of the evolutionary algorithm (which is an input parameter). The results expressed the running time of the proposed method is suitable to the problem needs. Therefore, it could be proved the proposed method achieved the specified objectives, in special, the treatment of practical aspects of the problem. Besides the proposed method, some contributions of this research are proposition of a new multi-objective evolutionary algorithm in subpopulation tables and a new reproduction operator to manipulate graph forests computationally represented by Node-Depth Encoding.