UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE CENTRO DE TECNOLOGIA PROGRAMA DE PÓS-GRADUAÇÃO EM ENGENHARIA ELÉTRICA E DE COMPUTAÇÃO PROPOSTA DE UMA NOVA ESTRUTURA MULTIFRACTAL COM BAIXA SENSIBILIDADE DA FREQUÊNCIA DE RESSONÂNCIA COM RELAÇÃO À PERIODICIDADE SERGIO ROBERTO DANTAS Orientador: Prof. Dr. Antônio Luiz Pereira de Siqueira Campos Dissertação de Mestrado apresentada ao Programa de Pós-graduação em Engenharia Elétrica e de Computação da UFRN (área de concentração: Telecomunicações) como parte dos requisitos para obtenção do título de mestre em ciências. Orientador: Prof. Dr. Antônio Luiz Pereira de Siqueira Campos Nº de ordem PPgEEC: M503 NATAL 2017 PROPOSTA DE UMA NOVA ESTRUTURA MULTIFRACTAL COM BAIXA SENSIBILIDADE DA FREQUÊNCIA DE RESSONÂNCIA COM RELAÇÃO À PERIODICIDADE Sergio Roberto Dantas Dissertação de mestrado aprovada em 06 de setembro de 2017 pela banca examinadora composta com os seguintes membros. Dedico este trabalho a minha mãe, Terezinha de Jesus Meneses, pois ela me amou e me ama desde meu nascimento. Tendo ela baixa renda, tive que morar com minha avó, mas o destino fez com que voltássemos a morar juntos, e lembro como hoje os momentos difíceis de nossas vidas. Essa memória muito recente em minha mente faz com que, vez por outra, eu lembre juntamente com ela dos momentos passados, o que nos traz alegria ao compararmos o passado e o presente. Assim, mãe querida, obrigado por tudo e que Jesus continue lhe dando saúde. Quero também dedicar este trabalho a minha segunda mãe, Ieda Maria Dantas, que de certa forma me educou e orientou em muitos momentos difíceis de minha vida. Vejo-a uma pessoa caridosa, sábia, compreensiva e muito mais. Assim, não poderia deixar de homenageá- la, sabendo que este trabalho é como se fosse dela. AGRADECIMENTOS A Deus primeiramente por sua infinita bondade e sabedoria. Ao meu orientador, professor Dr. Antônio Luiz Pereira de Siqueira Campos, por sua inestimável contribuição e orientações. Ao professor Dr. Érico Cadineli Braz pelas informações e dicas sobre os multifractais. Ao professor Dr. Alfredo Gomes Neto por ceder o laboratório, sem o qual não seria possível as medições e sua preciosa dica sobre como visualizar e configurar a malha no Ansoft Designer. Ao orientando José Jaime Guimarães Peixoto Neto pelas primeiras orientações de como manusear o Ansoft Designer. Ao colega de trabalho Aécio Vinicius Amorim Farias por entender e compreender a necessidade da flexibilização do horário de trabalho, sem o qual não seria possível pagar as disciplinas exigidas no mestrado. RESUMO Frequentemente, as superfícies seletivas em frequências têm recebido grande atenção por parte de pesquisadores em todo o mundo, devido a suas características de filtragem. Essas estruturas evoluíram em termos de fabricação. Assim, as geometrias utilizadas, inicialmente, eram fios em forma de grade e, posteriormente, diversos tipos de novas geometrias tais como: elementos retangulares, quadrados, circulares, anéis, espiras quadradas etc., foram utilizadas. Recentemente, elementos fractais e, mais recentemente, multifractais foram empregados. O emprego dessas geometrias é feito com o objetivo de se obter uma resposta em frequência específica. Parâmetros importantes, como largura de banda e frequência de ressonância são atendidos a partir do emprego de uma determinada geometria. Neste trabalho é proposta uma geometria multifractal com dupla similaridade, que tem como característica principal baixa sensibilidade à variação da periodicidade com relação à frequência de ressonância. Essa especificidade trará ganho no que diz respeito à elaboração do projeto, pois o controle da largura de banda é feito apenas com a alteração da periodicidade. Diversas simulações foram implementadas no Ansoft Designer e para efeito de comprovação, estruturas tradicionais foram simuladas, visando à demonstração do efeito da periodicidade na frequência de ressonância e, em seguida, dois protótipos foram construídos e feitas medições para validar experimentalmente as simulações. Palavras-chave: Superfícies seletivas em frequência, Fractal e Multifractal. ABSTRACT Often, frequency selective surfaces have received wide attention from researchers around the world because of their filtering characteristics. These structures evolved in terms of manufacturing. Thus, the geometries used, initially, were yarns in the form of a grid and later, several types of new geometries such as: rectangular, square, circular, rings, square, etc., were used. Recently, fractals and, more recently, multifractal elements have been employed. The use of these geometries is done in order to obtain a specific frequency response. Important parameters such as bandwidth and resonant frequency are met from the use of a given geometry. In this work a multifractal geometry with double similarity is proposed, which has as main characteristic low sensitivity to the variation of the periodicity with respect to the resonance frequency. This specificity will bring gain with regard to the elaboration of the project, because the control of the bandwidth is made only with the change of periodicity. Several simulations were implemented in Ansoft Designer and for the purpose of proving, traditional structures were simulated, aiming at demonstrating the effect of periodicity on the resonance frequency, and then two prototypes were constructed and measurements were taken to experimentally validate the simulations. Key words: Frequency selective surface, Fractal and Multifractal. LISTA DE FIGURAS Figura 1 – Elementos: (a) tipo aberturas e (b) tipo patch .......................................... 17 Figura 2 – Resposta em frequência: (a) FSS tipo aberturas e (b) FSS tipo patch..... 17 Figura 3 – Estrutura em cascata[20] ......................................................................... 17 Figura 4 – Grupo 1 .................................................................................................... 18 Figura 5 – Grupo 2 .................................................................................................... 19 Figura 6 – Grupo 3 .................................................................................................... 19 Figura 7 – Grupo 4 .................................................................................................... 19 Figura 8 – Tampa de forno como exemplo de aplicação de FSS[22] ........................ 21 Figura 9 – Papel bloqueador de rede sem fio [23] ..................................................... 22 Figura 10 – Vidro para janelas inteligentes [24] ........................................................ 23 Figura 11 – Sistema para medição de FSS ............................................................... 23 Figura 12 – Segmento de reta dividido em três partes. ............................................. 25 Figura 13 – Quadrado com quatro cópias de ½ ........................................................ 25 Figura 14 – Cubo de lado 1/2 .................................................................................... 26 Figura 15 – Conjunto de Cantor ................................................................................ 27 Figura 16 – Curva de Koch........................................................................................ 29 Figura 17 – Fractal tipo árvore .................................................................................. 30 Figura 18 – Fractal tipo Dürer .................................................................................... 31 Figura 19 – Fractal de Vicsek quadrados .................................................................. 32 Figura 20 – Fractal de Moran .................................................................................... 34 Figura 21 – Imagens digitalizadas do epitélio cervical [30] ....................................... 35 Figura 22 – Celular com antena fractal [30] .............................................................. 36 Figura 23 – Fractal usado em fibras óticas[30]. ......................................................... 37 Figura 24 – Multifractal de Cantor ............................................................................. 39 Figura 25 – Espectro multifractal de Cantor. ............................................................. 41 Figura 26 – Multifractal proposto ............................................................................... 42 Figura 27: Espectro multifractal DCB. ....................................................................... 43 Figura 28 – Medidas das geometrias:(a) espira quadrada, (b) dipolo cruzado e (c) patch quadrado ......................................................................................................... 45 Figura 29 – Resposta em frequência da transmissão para a FSS com espira quadrada e diferentes valores de p ........................................................................... 46 Figura 30 – Resposta em frequência da transmissão para a FSS com dipolo cruzado e diferentes valores de p ........................................................................................... 47 Figura 31 – Resposta em frequência da transmissão para a FSS com patch quadrado e diferentes valores de p ........................................................................... 48 Figura 32 – Geometria fractal com dupla razão de similaridade: (a) iniciador e (b) gerador ...................................................................................................................... 49 Figura 33 – Resposta em frequência da transmissão para a FSS com fractal proposto com razão de similaridade 𝒔𝟏 = 0,22 e 𝒔𝟐= 0,39 ........................................ 50 Figura 34 – Resposta em frequência da transmissão para a FSS com fractal proposto com razão de similaridade 𝒔𝟏 = 0,32 e 𝒔𝟐 = 0,34 ....................................... 51 Figura 35 – Resposta em frequência da transmissão para a FSS com fractal proposto com razão de similaridade 𝒔𝟏= 0,24 e 𝒔𝟐= 0,38. ........................................ 51 Figura 36 – Resposta em frequência da transmissão para a FSS com fractal proposto com razão de similaridade 𝒔𝟏= 0,246 e 𝒔𝟐= 0,377. .................................... 52 Figura 37 – Protótipos construídos ............................................................................ 53 Figura 38 – Setup de medição .................................................................................. 53 Figura 39 – Comparação entre resultados simulados e medidos para p = 14 mm. .. 54 Figura 40 – Comparação entre resultados simulados e medidos para p = 18 mm ... 54 Figura 41 – Resposta em frequência experimental da transmissão para as FSS construídas ................................................................................................................ 55 Figura 42 – Ilustração para dedução do fator de redução. ........................................ 67 LISTA DE TABELAS Tabela 1 – Conjunto de Cantor ................................................................................ 27 Tabela 2 – Curva de Koch ......................................................................................... 29 Tabela 3 – Fractal árvore .......................................................................................... 30 Tabela 4 – Fractal tipo Dürer ..................................................................................... 31 Tabela 5 – Fractal de Vicsek quadrado ..................................................................... 32 Tabela 6 – Fractal de Moran ..................................................................................... 34 Tabela 7 – Imagens digitalizadas .............................................................................. 36 LISTA DE ABREVIATURAS E SÍMBOLOS FDTD: Método das diferenças finitas no domínio do tempo. dB: Decibéis. 𝑝: Periodicidade. 𝑤: Largura do patch. 𝜆: Comprimento de onda no vácuo. 𝜙: Ângulo de incidência da onda eletromagnética. 𝜖: Permissividade elétrica. f: Frequência elétrica. µ: Permeabilidade magnética. Γ: Coeficiente de reflexão. 𝜏 : Coeficiente de transmissão. 𝐸𝑟 : Campo elétrico refletido. 𝐸𝑡 : Campo elétrico transmitido. 𝐸𝑖 : Campo elétrico incidente. d: Dimensão fractal ou dimensão local FSS: Superfície seletiva em frequência. fr: Frequência de ressonância BW: Largura de banda ρ1, ρ2 …: Similaridades 1,2… Wi-fi: Wireless Fidelity Bluetooth: Nome dado à tecnologia sem fio que usa ondas de radiofrequência. 3G: Terceira geração (padrão de tecnologia de telefonia móvel). Γ𝐴 Coeficiente de reflexão para FSS tipo abertura. Τ𝑃 Coeficiente de transmissão para FSS tipo patch. 𝑛 Nível do fractal. ( 𝑛 𝑝) Combinação binomial de 𝑛 elementos tomados 𝑝 a 𝑝. 𝑁(𝑟) Número de segmentos de comprimento r. 𝑟 Fator de redução. 𝑠 Dimensão Hausdorff. µ𝑘 Probabilidade na k- ésima iteração para um nível n. 𝑘 𝑘-ésima iteração, onde 𝑘 varia de 0 a n. 𝑠𝑘 Similaridades na k-ésima iteração #𝑠𝑘 Número de similaridades na k-ésima iteração. 𝑓𝑘 Função espetro multifractal na k-ésima iteração 𝑝1 Probabilidade no nível gerador. 𝑝2 Probabilidade no nível gerador. 𝑑𝑘 Dimensão local na k-ésima iteração EMI Interferência eletromagnética. RCS Seção transversal de radar. SUMÁRIO CAPÍTULO 1 INTRODUÇÃO .................................................................................... 14 CAPÍTULO 2 SUPERFÍCIES SELETIVAS EM FREQUÊNCIAS ............................... 16 2.1 INTRODUÇÃO .................................................................................................... 16 2.2 FORMATOS DOS ELEMENTOS PATCHES OU ABERTURAS ...................... 18 2.3 TÉCNICAS DE ANÁLISE ................................................................................. 20 2.4 APLICAÇÕES .................................................................................................. 21 2.5 SETUP DE MEDIÇÃO ..................................................................................... 23 CAPÍTULO 3 GEOMETRIAS FRACTAIS .................................................................. 24 3.1 INTRODUÇÃO ................................................................................................. 24 3.2 DIMENSÕES FRACTAIS ................................................................................. 24 3.3 FRACTAL CONJUNTO DE CANTOR.............................................................. 26 3.4 FRACTAL CURVA DE KOCH .......................................................................... 28 3.5 FRACTAL TIPO ÁRVORE ............................................................................... 30 3.6 FRACTAL TIPO DÜRER ................................................................................. 31 3.7 FRACTAL TIPO VICSEK ................................................................................. 32 3.8 FRACTAIS COM DUPLA SIMILARIDADE ....................................................... 33 3.9 APLICAÇÕES DE FRACTAIS ......................................................................... 35 CAPÍTULO 4 MULTIFRACTAIS ................................................................................ 38 4.1 INTRODUÇÃ AOS MULTIFRACTAIS. ............................................................. 38 4.2 ANÁLISES DO MULTIFRACTAL DE CANTOR .................................................. 39 4.3 ANÁLISE DO MULTIFRACTAL PROPOSTO ...................................................... 41 4.4 APLICAÇÕES DE MULTIFRACTAIS .................................................................. 44 CAPÍTULO 5 RESULTADOS .................................................................................... 45 5.1 RESPOSTAS EM FREQUÊNCIA DE GEOMETRIAS TRADICIONAIS ........... 45 5.2 GEOMETRIA PROPOSTA E SUA RESPOSTA EM FREQUÊNCIA ............... 48 CAPÍTULO 6 CONCLUSÕES ................................................................................... 56 APÊNDICE I .............................................................................................................. 61 APÊNDICE II ............................................................................................................. 66 14 CAPÍTULO 1 INTRODUÇÃO Superfícies seletivas em frequência (Frequency Selective Surface – FSS) são arranjos bidimensionais, de células unitárias, com patches metálicos depositados sobre um material dielétrico, ou o complemento desses patches, que seriam as aberturas. A FSS do tipo patch, geralmente, transmite parte das frequências e bloqueia outras frequências de ondas que se propagam no espaço livre. Ela age como um filtro rejeita-faixa [1] e são conhecidas, também, como filtros espaciais [2]. Os avanços em sistemas de comunicação modernos demandam novas FSS [3] e as aplicações das FSS incluem absorção de micro-ondas, redução da seção transversal de radar (RCS), redução de interferência electromagnética (EMI), aplicações em ondas milimétricas e submiliméttricas etc. [1], [2], [4]. Além disso, nos últimos anos, muitas pesquisas têm sido conduzidas em diferentes tipos de estruturas eletromagnéticas tradicionais baseadas em FSS. Essas estruturas podem ser utilizadas para processamento de sinais analógicos, em campos como estruturas electromagnéticas artificiais, para o controle e dispersão de ondas eletromagnéticas em tempo real. Para o projeto de FSS, existem vários tipos de geometrias de patches. Em [2], Munk dividiu geometrias em quatro grupos: N polos, loops, interior sólido, e combinações. Muitos autores têm investigado esses tipos de geometrias [5] - [10]. Recentemente, um novo grupo de geometrias foi aplicado para projetar FSS, os fractais [11] - [16]. Mais recentemente, estudos analisaram uma variação de fractais, os multifractals [17] - [19]. Todas essas investigações mostraram que a resposta em frequência das FSS tem uma forte dependência, em termos de frequência de ressonância e largura de banda, com o tipo de geometria, dimensões físicas da geometria e periodicidade da célula unitária. Na literatura, as pesquisas mostram que, para geometrias comuns, como quadrado sólido, espira quadrada, dipolo cruzado, entre outras, a alteração da periodicidade traz como consequência um deslocamento na frequência de ressonância, assim como na largura de banda. Isso implica que, mudanças de periodicidades, visando mudanças de largura de banda, implicam em mudança da frequência de ressonância. Caso se deseje alterar a largura de banda sem alterar a 15 frequência de ressonância, pode ser necessário reprojetar todas as dimensões físicas da FSS. Assim, neste trabalho propõe-se uma nova geometria multifractal, que permite a alteração da periodicidade sem alterar a frequência de ressonância em uma considerada faixa de periodicidade. Essa característica, não vista em geometrias tradicionais, até o presente momento, tem a vantagem de simplificar o projeto de filtros espaciais, pois com a alteração da periodicidade controla-se a largura de banda e sua frequência de ressonância permanece fixa. Assim, este trabalho está organizado da seguinte forma:  O Capítulo 2 descreve as definições de superfícies seletivas em frequência, parâmetros constitutivos de projetos, aplicações práticas, setup de medição e alguns métodos numéricos.  O Capítulo 3 conceitua os principais aspectos relacionados a geometrias fractais, processos de construção de fractais, aplicações e a noção de dimensão fractal.  O Capítulo 4 apresenta os conceitos sobre multifractais, as leis que regem essas estruturas, definições, aplicações e exemplo de cálculo de dimensão.  O Capítulo 5 apresenta os resultados simulados e medidos, para geometrias convencionais e para a geometria multifractal proposta.  Por fim, o Capítulo 6 apresenta as principais conclusões acerca do trabalho e propostas sobre a sua continuidade. 16 CAPÍTULO 2 SUPERFÍCIES SELETIVAS EM FREQUÊNCIAS 2.1 INTRODUÇÃO As FSS são estruturas periódicas feitas de patches ou elementos do tipo abertura [1] – [3]. Essas superfícies funcionam de forma semelhante a filtros de RF, com algumas diferenças no seu funcionamento. Enquanto os filtros de RF processam sinais elétricos com tensões ou correntes, as FSS processam ondas eletromagnéticas se propagando no espaço livre. Esses campos ao incidir nessas superfícies ressoam, transmitindo ou refletindo a onda eletromagnética. As FSS básicas funcionam como filtros passa-faixa ou rejeita-faixa, dependendo do tipo de elemento empregado na sua construção. A Figura 1 ilustra as estruturas físicas e a Figura 2 ilustra a resposta em frequência dessas FSS. A Figura 1(a) ilustra uma FSS com elementos do tipo abertura retangular (área branca), enquanto que a Figura 1(b) é constituída de elementos do tipo patch condutor retangular (área preta), sendo esses, normalmente, depositados sobre algum material dielétrico. Na Figura 2(a), tem-se a FSS do tipo abertura, que possui uma resposta em frequência do tipo passa-faixa, enquanto que na Figura 2(b) verifica-se que a FSS do tipo patch possui uma resposta em frequência do tipo rejeita-faixa. Além dessas configurações, tem-se diversas outras possibilidades para a construção de FSS, por exemplo, podem ser usados diferentes materiais dielétricos, o que pode alterar as características de transmissão e de reflexão. É possível usar superestratos dielétricos como em [1], ou seja, colocar o elemento condutivo no meio de dois dielétricos. Pode-se ainda empilhar, alternadamente, materiais dielétricos e grades condutivas sucessivamente. Enfim, as possibilidades de construção são inúmeras. A Figura 3 ilustra uma FSS cascateada fabricada com a alternância de dielétricos e grades condutivas, sendo que, nesse caso, foram usados dois dielétricos e duas grades condutivas. 17 Figura 1 – Elementos: (a) tipo aberturas e (b) tipo patch Fonte: imagem elaborada pelo autor, 2017. Figura 2 – Resposta em frequência: (a) FSS tipo aberturas e (b) FSS tipo patch Fonte: imagem elaborada pelo autor, 2017. Figura 3 – Estrutura em cascata[20] Fonte: imagem elaborada pelo autor, 2017. 18 Existe ainda uma propriedade importante com relação a FSS. Essa propriedade é conhecida como princípio de Babinet [4]. Estruturas complementares, ou seja, do tipo abertura e do tipo patch e sobre dielétricos de mesma característica, produzem respostas em frequência complementares. Não obstante, para que esse princípio ocorra é necessário que o condutor seja perfeito e infinitamente fino, tipicamente menor que 1/1000 do comprimento de onda. Outro limite é com relação ao dielétrico, ou seja, para uma espessura dielétrica maior que um quarto do comprimento de onda (λ/4), o comportamento é completamente diferente. Além disso, as características dos coeficientes de transmissão e reflexão também são completares na ressonância, ou seja, Γ𝐴 = Τ𝑃. Em que, Γ𝐴 é o coeficiente de reflexão da FSS do tipo abertura e Τ𝑃 é o coeficiente de transmissão da FSS do tipo patch. 2.2 FORMATOS DOS ELEMENTOS PATCHES OU ABERTURAS Quanto à apresentação das formas geométricas, a literatura apresentada por [1] – [3] classifica em quatros grupos, como segue: - Grupo 1 dos N-polos conectados pelo centro, como: dipolos, dipolos cruzados, cruz de Jerusalém, tripolos etc., como ilustrado na Figura 4. Figura 4 – Grupo 1 Fonte: imagem elaborada pelo autor, 2017. 19 - Grupo 2 das Espiras, como: Espiras quadradas, espiras quadradas duplas, espira quadrada com grade, anéis circulares concêntricos duplos etc., ilustrado na Figura 5. Figura 5 – Grupo 2 Fonte: imagem elaborada pelo autor, 2017. - Grupo 3 dos elementos de interiores sólidos, como: quadrado, retângulo, hexágono, círculo, e outros ilustrados na Figura 6. Figura 6 – Grupo 3 Fonte: imagem elaborada pelo autor, 2017. - Grupo 4 dos elementos combinados como os ilustrados na Figuras 7. Figura 7 – Grupo 4 Fonte: imagem elaborada pelo autor, 2017. 20 2.3 TÉCNICAS DE ANÁLISE Várias técnicas de análise têm sido usadas em análises de FSS. Há na literatura diversas técnicas de análise, efetuadas em anteparos periódicos. Fórmulas aproximadas e simples foram desenvolvidas por alguns autores, para se determinar as características de transmissão e de reflexão para uma FSS com patches condutores ou aberturas retangulares. Esse tipo de análise diminui o tempo computacional e produz resultados satisfatórios para determinadas estruturas [4], [5]. Entretanto, essas fórmulas têm suas limitações e, para determinadas aplicações, podem produzir resultados imprecisos. Um método simples e que produz resultados satisfatórios é o modelo do circuito equivalente. Nessa análise os vários segmentos de fita que formam o elemento patch em um arranjo periódico são modelados como componentes indutivos ou capacitivos em uma linha de transmissão. Da solução desse circuito, são encontradas as características de transmissão e reflexão da FSS. Essa técnica usa uma aproximação quase-estática para calcular os componentes do circuito e permite uma análise computacional muito rápida [6] – [8]. Outro método empregado é o da expansão modal [9], que permite uma análise capaz de fornecer detalhes das repostas da frequência e da polarização, junto com o entendimento físico da sua operação. Uma técnica bastante difundida, atualmente, é a técnica das diferenças finitas no domínio do tempo (FDTD). Essa técnica possibilita a análise de qualquer tipo de elemento, bem como a análise de perdas dielétricas e/ou magnéticas e a análise de estruturas não homogêneas [10]. A desvantagem dessa técnica é o grande esforço computacional despendido. Métodos híbridos foram empregados em 2001. Técnicas como interpolação racional de Krylov [11], [12], Método dos Momentos com o Método BI – RME [13] e Método dos Momentos em conjunto com elementos finitos [14], estão sendo muito usados na modelagem de elementos de forma complexa. Técnicas de inteligência artificial também têm sido usadas. Algoritmos genéticos [15], [16] podem ser empregados na análise e/ou síntese de superfícies seletivas de frequência. Outra técnica empregada é o método da linha de transmissão equivalente [17] em conjunto com o método de Galerkin [18]. Essa é uma análise de onda completa, 21 que produz resultados precisos, além de facilitar, relativamente, a manipulação matemática [19], [21]. 2.4 APLICAÇÕES A aplicação certamente mais conhecida de FSS é a grade perfurada que fica por trás do vidro e faz parte da tampa do forno de micro-ondas. A grade vazada em forma circular é uma FSS tipo abertura, através da qual a frequência da luz atravessa a tampa metálica, mas reflete as frequências de micro-ondas em todas as paredes, fazendo vibrar as moléculas da água e retirando água da comida ou das estruturas, preparando dessa forma as refeições. A Figura 8 mostra o forno e a tampa metálica vazada, ou a FSS propriamente dita. Figura 8 – Tampa de forno como exemplo de aplicação de FSS[22] Fonte: imagem elaborada pelo autor, 2017. Papéis de paredes construídos com pequenos cristais de prata são usados para bloquear sinais na faixa de frequência de redes sem fio e deixar passar outros sinais. Isso permite o confinamento do sinal da rede WLAN dentro do recinto. A Figura 9 ilustra essa estrutura proposta fabricada. Segundo o fabricante, se for pintado, o papel não perde sua propriedade. Portanto, é possível escolher o papel que quiser sem se preocupar com a estética das paredes [23]. O fabricante comercializa esses papéis para bloquear as frequências de 2,4GHz ou 5GHz. 22 Figura 9 – Papel bloqueador de rede sem fio [23] O MIT (Massachusetts Institute of Technology) publicou em 23 de outubro de 2015 a criação de um vidro desenvolvido por um grupo liderado por Delia Milliron, professora de engenharia química da universidade do Texas em Austin. Esse vidro pode bloquear seletivamente o componente de produção de calor, bem como a luz visível. Milliron diz que o desempenho agora é suficientemente bom para iniciar os planos de construir um protótipo de linha de produção com base nesses avanços recentes. A chave para o vidro da janela inteligente é um "quadro" de nanocristais feito de um material eletricamente condutivo incorporado em um material vítreo. Os nanocristais e o material vítreo têm propriedades ópticas distintas, que mudam quando os materiais são carregados ou descarregados eletronicamente. Os nanocristais podem bloquear a luz próximo ao infravermelho ou permitir que ela passe, enquanto o material vítreo pode fazer a transição entre um estado transparente que bloqueia a luz visível. A publicação informa que a produção do vidro iniciará em 2017 [24]. A Figura 10 ilustra esse material. 23 Figura 10 – Vidro para janelas inteligentes [24] 2.5 SETUP DE MEDIÇÃO Existem basicamente três tipos de setup na literatura. Cada uma com suas peculiaridades. Entretanto, neste trabalho foi usado um esquema semelhante ao apresentado na Figura 11. Conforme pode ser visto, as medições são feitas utilizando um analisador vetorial de redes, duas antenas cornetas de ganho padrões e um painel para minimizar as perdas por difração e reflexões. Figura 11 – Sistema para medição de FSS Fonte: imagem elaborada pelo autor, 2017. 24 CAPÍTULO 3 GEOMETRIAS FRACTAIS 3.1 INTRODUÇÃO No passado, matemáticos se preocupavam com conjuntos e funções com métodos clássicos. Conjuntos e funções não regulares eram considerados patológicos e não merecedores de estudos, [25]. Entretanto essa concepção mudou e percebeu- se que conjuntos fractais representam de maneira adequada diversos fenômenos naturais, como a deposição de sedimentos. Segundo o dicionário Aurélio, fractal é um conjunto geométrico ou objeto natural cujas partes têm a mesma estrutura (irregular e fragmentada) que o todo, mas em escalas diferentes [26], ou diz-se do objeto que é irregular e fragmentado. Bernoit Mandrebrot, matemático francês, criou em 1975 o termo fractal que significa quebrar, do adjetivo latino fractus, do verbo frangere [25], [27]. Outra definição de Mandelbrot é que um conjunto fractal tem dimensão de Hausdorff estritamente maior que a dimensão topológica [27]. Já o Falconer, assim como outros matemáticos, prefere evitar uma definição rígida e assinala cinco características de um fractal [25]: I. Tem uma estrutura fina, isto é, escalas menores mostram arbitrariamente detalhes do conjunto. II. É também irregular na descrição geométrica local e global. III. Frequentemente tem alguma forma irregular de autossimilaridade, talvez aproximada ou estatística. IV. Usualmente a dimensão fractal é maior que a dimensão topológica. V. Em muitos casos particulares é definido de forma simples, talvez recursivamente. 3.2 DIMENSÕES FRACTAIS A maneira mais simples de aprender a dimensionar objetos no espaço é usar os conceitos já formulados da geometria euclidiana, na qual um ponto tem dimensão zero, uma reta tem dimensão um, um quadrado tem dimensão dois e um volume tem dimensão três. Para dimensionar esses elementos inicia-se com uma linha de 25 tamanho unitário, conforme ilustrado na Figura 12. Em seguida divide-se em três segmentos e aplica-se: 𝑁(𝑟) = 1/𝑟𝑑 (1) onde 𝑁(𝑟) é a quantidade de segmentos de comprimento 𝑟, e 𝑟 é o comprimento de cada segmento ou fator de redução e o expoente de r é a dimensão. Para que 𝑁(𝑟) seja três, o expoente deve ser um e 𝑟, um terço. Figura 12 – Segmento de reta dividido em três partes. 0 1/3 2/3 1 Fonte: imagem elaborada pelo autor, 2017. Da mesma forma, para um quadrado de área unitária e lado unitário, conforme ilustrado, na Figura 13. Esse quadrado foi dividido lateralmente em 1/2 gerando quatro cópias idênticas. Figura 13 – Quadrado com quatro cópias de ½ ½ ½ Fonte: imagem elaborada pelo autor, 2017. Assim, tem-se quatro quadrados, menores de lado 1/2, que equivale a um expoente dois, conforme a expressão (1). Na Figura 14, tem-se um cubo de lado 1/2 e oito cópias com comprimento lateral de ½. Usando mais uma vez a expressão (1), o expoente deve ser três. De modo que, para um objeto de dimensão um, dois, três, quatro ou cinco, tem-se um expoente um, dois, três, quatro ou cinco. Esse expoente é normalmente identificado pela letra d 26 (dimensão do conjunto). A dimensão é obtida calculando o logaritmo em ambos os lados da expressão (1), resultando em: 𝑙𝑜𝑔𝑁(𝑟) = log⁡( 1 𝑟𝑑 ) 𝑑 = log⁡𝑁(𝑟))/ ⁡log ( 1 𝑟 ) (2) Figura 14 – Cubo de lado 1/2 ½ ½ Fonte: imagem elaborada pelo autor, 2017. 3.3 FRACTAL CONJUNTO DE CANTOR Georg Ferdinand Ludwig Philip Cantor nasceu em São Petersburgo, na Rússia, em 3 de março de 1845 e faleceu em Halle, Alemanha, em 6 de janeiro de 1918. A teoria dos conjuntos criada por Cantor é uma das mais notáveis inovações matemáticas dos últimos séculos. Apresentada em pleno século 19, foi combatida pelos contemporâneos do matemático, entre eles, seu maior inimigo, Leopol Kronecker [28]. A Figura 15 ilustra o conjunto de Cantor. O nível zero é composto de um segmento de reta de tamanho unitário. No nível um, esse segmento unitário é dividido em três partes iguais e retira-se a parte central. No nível dois, repete-se o processo para cada segmento obtendo quatro segmentos e assim sucessivamente. Assim, verifica-se que no enésimo termo o fator de redução é 3−𝑛 e o comprimento total do é ( 2 3 ) 𝑛 , conforme a Tabela 1. Dessa forma, quando n tende ao infinito, a linha torna-se praticamente um ponto, ou seja, para n tendendo ao infinito o comprimento do conjunto tende a zero. Isso pode ser verificado na terceira coluna da Tabela 1. Nessa, 27 a primeira coluna indica o nível fractal, na segunda, tem-se o fator de redução e na terceira, o comprimento total para cada nível. Figura 15 – Conjunto de Cantor Fonte: imagem elaborada pelo autor, 2017. Tabela 1 – Conjunto de Cantor Nível Fator de redução Comprimento total 0 1 1 1 1 3 2 3 2 ( 1 3 ) 2 ( 2 3 ) 2 3 ( 1 3 ) 3 ( 2 3 ) 3 𝑛 ( 1 3 ) 𝑛 ( 2 3 ) 𝑛 Fonte: tabela elaborada pelo autor, 2017. Já quanto a sua dimensão fractal, basta utilizar a equação (2) para calcular a dimensão do conjunto fractal, logo, para o conjunto de cantor da Figura 15, e utilizando-se do nível um, tem-se duas cópias e fator de redução de um terço. Portanto: 𝑑 = 𝑙𝑜𝑔2 𝑙𝑜𝑔3 𝑑 = 0,63 28 Então, verifica-se que o conjunto de Cantor tem uma dimensão fracionária, que é uma característica de um conjunto fractal. 3.4 FRACTAL CURVA DE KOCH O matemático Niels Fabian Helge von Koch nasceu no dia 25 de janeiro de 1870 em Estocolmo, Suécia. Koch estudou a autossimilaridade no final do século 19 e ficou conhecido em 1904, através de um artigo publicado sobre o processo de criação de curvas contínuas sem tangentes ou semitangentes em nenhum de seus pontos. Atualmente, essa curva contínua e sem nenhuma tangente é conhecida como curva de Koch. Essa curva, no nível zero, se inicia com um segmento de reta unitário. No primeiro nível, esse segmento de reta é dividido em três partes iguais, sendo que o segmento central é removido e adiciona-se dois segmentos formando um ângulo entre eles de 60º, assim tem-se quatro segmentos de comprimentos iguais conforme iteração um da Figura 16. É fácil perceber que o fator redução é 1/3 e comprimento total tem 4/3. No segundo nível, toma-se cada um dos quatros segmentos e repete- se o passo anterior. Logo o fator de redução é 1/9 e comprimento total é 16/9. Já no enésimo nível temos um comprimento de (4/3)n, conforme Tabela 2, onde 𝑛 é o nível fractal. Observando novamente, quando 𝑛 cresce, embora o comprimento de cada segmento tenda a zero, o comprimento total tende a infinito. 29 Figura 16 – Curva de Koch Iteração 0 Iteração 1 Iteração 2 Iteração 3 Fonte: imagem elaborada pelo autor, 2017. Tabela 2 – Curva de Koch Nível 0 1 2 3 𝑛 Número de segmentos 1 4 42 43 4𝑛 Comprimento de cada segmento 1 1 3 ( 1 3 ) 2 ( 1 3 ) 3 ( 1 3 ) 𝑛 Comprimento total da curva 1 4 3 ( 4 3 ) 2 ( 4 3 ) 3 ( 4 3 ) 𝑛 Fonte: Tabela elaborada pelo autor, 2017. Para o cálculo da dimensão da curva de Koch, na Figura 16, procede-se como no conjunto de Cantor, ou seja, toma-se no nível um a quantidade de seguimentos similares e o fator de redução, logo: 𝑑 = 𝑙𝑜𝑔4 𝑙𝑜𝑔3 𝑑 = 1,26 30 3.5 FRACTAL TIPO ÁRVORE Na Figura 17, mostram-se dois fractais com o mesmo fator de escala, mudando apenas o ângulo entre os dois segmentos. O fator de escala é 0,7. No nível zero tem- se um segmento unitário, assim como nos fractais anteriores. No nível um, foram acrescidos dois segmentos com ângulos de 90º entre si. No segundo nível tem-se mais quatro segmentos com fator de 0,49 e ângulos de 90º, no nível seguinte repete- se o procedimento indefinidamente. Já o fractal tipo árvore ao lado direito a sua construção é feita de forma semelhante, com a diferença que os segmentos adicionados têm um ângulo de 180º. Na Tabela 3 tem-se a sequência de formação. Figura 17 – Fractal tipo árvore Fonte: imagem elaborada pelo autor, 2017. Tabela 3 – Fractal árvore Nível 0 1 2 …… 𝑛 Número de segmentos 1 2 22 …… 2𝑛 Comprimento de cada segmento 1 (0,7) (0,7)2 ….. (0,7)𝑛 Fonte: tabela elaborada pelo autor, 2017. Novamente, para o fractal tipo árvore da Figura 17, têm-se duas cópias e fator de redução de 7/10, por conseguinte a dimensão fractal é: 𝑑 = 𝑙𝑜𝑔2 log⁡(10/7)⁡ 𝑑 = 1,94 Os fractais tipo árvore merecem uma observação, pois sua formação é por fronteira assim como a curva de Koch, mas um fato que chama a atenção é o caso para o fator de redução de 1 √2 . Nesse, tem-se a dimensão dois, o que caracteriza uma 31 dimensão inteira. Logo, poder-se-ia deduzir erroneamente que não seria um fractal, mas uma figura euclidiana. 3.6 FRACTAL TIPO DÜRER A Figura 18 ilustra três níveis do fractal tipo Dürer hexagonal. O nível zero tem lado unitário. O nível um tem-se o lado com fator de escala de 1/3 e seis cópias semelhantes. O nível dois tem-se o fator de 1/9 e doze cópias reduzidas. A Tabela 4 apresenta o processo de construção para esse fractal. De forma semelhante podem- se construir fractais tipo Dürer com pentágonos, heptágonos, octógonos etc. Figura 18 – Fractal tipo Dürer Fonte: imagem elaborada pelo autor, 2017. Tabela 4 – Fractal tipo Dürer Nível 0 1 2 …… 𝑛 Número de hexágonos 1 6 62⁡ …… 6𝑛 Comprimento do lado do hexágono 1 1 3 ( 1 3 ) 2 ….. ( 1 3 ) 𝑛 Fonte: tabela elaborada pelo autor, 2017. No caso do fractal de Dürer, apresentado na Figura 18, tem-se seis cópias e um fator de redução de um terço no nível um, assim: 𝑑 = 𝑙𝑜𝑔6 𝑙𝑜𝑔3 𝑑 = 1.63 32 3.7 FRACTAL TIPO VICSEK A Figura 19 mostra um fractal tipo Vicsek com elementos quadrados. O nível zero inicia-se com um quadrado de lado unitário. O nível um divide-se em nove quadrados retirando-se quatro quadrados com fator de escala 1/3. No nível dois, procede-se da mesma forma ficando com 25 quadrados e fator de escala 1/9. Assim, procede-se indefinidamente. Na Tabela 5 apresenta-se algumas características de construção. Figura 19 – Fractal de Vicsek quadrados Fonte: imagem elaborada pelo autor, 2017. Tabela 5 – Fractal de Vicsek quadrado Nível 0 1 2 …… 𝑛 Número de quadrados 1 5 52 …… 5𝑛 Comprimento do lado do quadrado 1 1 3 ( 1 3 ) 2 ….. ( 1 3 ) 𝑛 Fonte: tabela elaborada pelo autor, 2017. E para dimensionar o fractal de Vicsek da Figura 19, observa-se que o número de cópias são cinco e um fator de redução de um terço, logo: 𝑑 = 𝑙𝑜𝑔5 𝑙𝑜𝑔3 𝑑 = 1.46 33 3.8 FRACTAIS COM DUPLA SIMILARIDADE O multifractal com dupla similaridade apresentado será o fractal de Moran [29]. Esse fractal é construído a partir de um círculo de raio unitário e do ponto de vista do comprimento linear e de quantidades de cópias similares tem um comportamento diferente dos fractais anteriores. A Figura 20 ilustra três níveis desse fractal. No nível zero, tem-se um círculo de raio unitário. No nível, um têm-se quatro elementos e duas cópias, duas a duas com similaridades, 1 3 ⁡𝑒⁡ 1 2 . No nível, dois tem-se dezesseis elementos e três similaridades, 1 9 , 1 ⁡6 ⁡𝑒⁡ 1 4 . É simples perceber que o número de similaridades é 𝑛 + 1 e a quantidade de elementos é dada por 4𝑛. Para calcular as similaridades, consulte no apêndice I deste trabalho, cálculos de estruturas similares são feitos detalhadamente e servirão para o exemplo em questão. Outra informação é a quantidade de elementos similares que se repete k vezes e dado por: 2𝑛 ( 𝑛 𝑘 ), onde 𝑛 é o enésimo nível fractal e 𝑘 varia de 0 a 𝑛. Assim, para cada 𝑘 tem-se a quantidade de elementos similares. A Figura 20 mostra três níveis fractais do fractal de Moran e no capítulo sobre multifractais serão abordados ainda outros fatores, como as dimensões locais e o espectro multifractal, tópicos inerentes a esses tipos de estruturas. A Tabela 6 ilustra algumas características desse fractal. 34 Figura 20 – Fractal de Moran Fonte: imagem elaborada pelo autor, 2017. Tabela 6 – Fractal de Moran Nível 0 1 2 …… 𝑛 Número de círculos 1 4 42 …… 4𝑛 Número de similaridades 1 2 3 ….. 𝑛 + 1 Total de elementos similares para 𝑘(0⁡𝑎⁡𝑛) 1 2,2 4,8,4 ….. 2𝑛 ( 𝑛 𝑘 ) Fonte: tabela elaborada pelo autor, 2017. No entanto, para o fractal de Moran da Figura 20, o cálculo da dimensão fractal torna-se bastante simples quando se aplica a dimensão Hausdorff, que abrange os fractais com dupla ou mais autossimilares, conforme [29]: Na referência [25], o autor usa a fórmula ∑ 𝑐𝑖 = 1 𝑚 𝑖=1 que calcula a dimensão fractal, tanto para uma similaridade, como para duas, três ou mais. Nessa fórmula 𝑐 é a razão de similaridade. Dessa forma a dimensão do conjunto de Moran, fica: 2 ( 1 3 ) 𝑠 + 2( 1 2 ) 𝑠 = 1 𝑠 ≈ 1.60 Entretanto, deve-se ter o cuidado, pois o cálculo dessa dimensão caracteriza o fractal como um todo, porém para ter-se mais informação, deve-se usar distribuição de probabilidade, o que será abordado posteriormente no Capítulo 4. 35 3.9 APLICAÇÕES DE FRACTAIS Na medicina a dimensão fractal é usada como método de diagnóstico quantitativo de várias patologias. Um dos campos mais desenvolvidos é o diagnóstico do cancro. As evidências experimentais sugerem que se use a linha de investigação da detecção de núcleos atípicos para tumores [30]. A Figura 21 mostra algumas imagens digitalizadas e na tabela 7, tem a dimensão 0,97 representando uma possível patologia. Figura 21 – Imagens digitalizadas do epitélio cervical [30] 36 Tabela 7 – Imagens digitalizadas Tecido Dimensão fractal Núcleos de controle 0,97 Núcleos típicos 1,47 Fonte: tabela elaborada pelo autor, 2017. Outra aplicação de fractais é no projeto de antenas de microfita. A Motorola começou a usar fractais em antenas e anunciou uma eficiência de aproximadamente 25% com relação às tracionais que usam pedaços de fio condutor [30]. A Figura 22 ilustra uma antena usada em celular. Figura 22 – Celular com antena fractal [30] Por fim, o empacotamento apropriado de fibras ópticas produz guias de ondas com uma baixíssima distorção. Lee Coolk, da Galileo Electro-Optics Corp, mostrou, através do uso de pavimentações recursivas, que os melhores empacotamentos de fibras ópticas são aqueles que têm bordas fractais. Isso levou ao desenho de feixes de fibras ópticas fractais multi-multifibras, os quais exibem um melhor contraste de imagem [30]. O motivo da pavimentação (a forma dos azulejos) é construído recursivamente, conforme Figura 23. 37 Figura 23 – Fractal usado em fibras óticas[30]. 38 CAPÍTULO 4 MULTIFRACTAIS 4.1 INTRODUÇÃ AOS MULTIFRACTAIS. Os monofractais fractais, embora expliquem alguns fenômenos naturais, não servem para descrever alguns fenômenos, como: tremores de terra, turbulências em um fluído, redes neurais, sedimentos de rochas, variações de aplicações em bolsas de valores, distribuição de espaço no universo [27], [29] e outros são modelados com estruturas denominadas multifractais. Os multifractais, como o nome sugere, são estruturas compostas por vários fractais e possuem características semelhantes aos fractais, como: invariância de escalas e autossimilaridades. Os conjuntos multifractais têm tratamento matemático com distribuição de probabilidades, ou seja, enquanto os fractais são tratados como conjuntos de pontos, os multifractais são tratados como conjuntos de distribuição de probabilidade. A análise dos multifractais é dividida em análise fina e grossa. No caso da fina, a preocupação é nas distribuições locais ou na estrutura que surge de si mesmo. Na grossa, a preocupação é em relação às irregularidades e, para isso, é costume, ao invés de uma dimensão fractal haver um espectro multifractal. Ressalva-se ainda que, no livro escrito por Mandelbrot (1977-1983), ele descreveu objetos com formas irregulares. Seus exemplos incluem galáxias, comprimentos de linha costeira, dentre outros. Esses objetos têm características aparentemente bizarras, como linhas costeiras, que têm comprimentos infinitos e flocos de neve, que têm área infinita. Essa tendência ocorre quando o conjunto é irregular. Assim, uma maneira de descrever esses conjuntos é através de sua dimensão fractal. Para caracterizar um multifractal, considere uma massa µ espalhada sobre uma região, em tal caminho, onde a concentração de massa é irregular. Analisando esse conjunto em vários locais, percebe-se que essa concentração é regida por uma lei de potência µ(𝐵(𝑥, 𝑟)) = ⁡ 𝑟𝛼 e quando r tende a zero, determinam-se os vários índices, que são as dimensões locais e uma medida com tal riqueza é denominada de medida multifractal, ou simplesmente multifractal. Na fórmula anterior µ(𝐵(𝑥, 𝑟)), é a probabilidade de cobrir com uma estrutura com uma esfera de centro 𝑥 e raio 𝑟⁡, e α é a dimensão fractal. 39 4.2 ANÁLISES DO MULTIFRACTAL DE CANTOR Como um primeiro exemplo, será realizada a análise do multifractal de Cantor. Inicialmente a partir de um segmento unitário, retira-se um terço da parte central e aloca-se esse elemento em cima do elemento à direita a partir do nível um. No nível dois, os segmentos anteriores são novamente partidos e o segmento do meio são alocados no lado direito para cada segmento. Assim os próximos níveis serão formados da mesma forma. É fácil notar que a formação desse conjunto tem uma taxa de 1/3 para 2/3, sendo um terço para o segmento esquerdo e dois terços para o segmento direito. É possível também perceber que para cada nível temos 𝑛 + 1 probabilidades e 3−𝑛 similaridades formadas. Dessa forma, para o nível um, tem-se uma similaridade de um terço e a probabilidade desse elemento, também é um terço e o elemento à direita, dois terços. No nível dois, tem-se a razão de similaridade de um nono e as seguintes probabilidades da esquerda para direita: 1 9 , 2 9 , 𝑒 4 9 , conforme ilustrado na Figura 24. Figura 24 – Multifractal de Cantor 1/90 2/9 1/3 2/3 7/9 8/9 1 0 1/3 2/3 1 0 1 Fonte: imagem elaborada pelo autor, 2017. Continuando, percebe-se que esse conjunto tem um comportamento binomial. Para calcular essas probabilidades, pode-se usar a fórmula da distribuição binomial, como segue: µ𝑘 = (𝑝1) 𝑘(𝑝2) 𝑛−𝑘 (3) 40 onde,𝑛 é o nível do fractal e 𝑘 varia de 0 a 𝑛, 𝑝1e 𝑝2 são probabilidades no primeiro nível. Assim, para cada, tem-se uma probabilidade. É necessário também definir o número de similaridades que será usado para o cálculo da dimensão e pode ser dado por: #𝜇𝑘 = ( 𝑛 𝑘 ) (4) onde o termo dentro do parênteses é o número binomial dado por: 𝑛! 𝑘!(𝑛−𝑘)! . Dando continuidade à análise do multifractal, sua dimensão local pode ser calculada por: 𝑑𝑘 = log(µ𝑘) / log(𝑟 𝑛) (5) Onde 𝑟 é 1/3, 𝑛, o enésimo nível fractal e k varia de 0 a 𝑛 conforme já descrito. E por último a função de singularidade ou espectro multifractal é dado por: 𝑓𝑘 = log#(μk) −log⁡(rn) (6) onde todas incógnitas já foram definidas. Essa função está ilustrada na Figura 25 e percebe-se que quando, 𝑛 tende ao infinito a dimensão do conjunto é 0,63, ou seja, a dimensão do monofractal do conjunto de Cantor. Também se verifica que esse conjunto possui infinitas dimensões. Para um cálculo detalhado ver Apêndice I. 41 Figura 25 – Espectro multifractal de Cantor. Fonte: imagem elaborada pelo autor, 2017. 4.3 ANÁLISE DO MULTIFRACTAL PROPOSTO A Figura 26 ilustra o fractal com dupla similaridade proposto neste trabalho, cuja denominação dada foi multifracal DCB (Dantas Campos Braz). No nível inicial (ou zero), tem-se um círculo com raio unitário e um único elemento. No nível seguinte, tem-se duas similaridades, uma de um quinto e outra de dois quintos e um total de cinco elementos. No próximo nível, tem-se três similaridades: 1/25, 2/25, 4/25 e 25 elementos, de tal forma que o número de similaridades é dado por 𝑛 + 1,conforme já visto, em que 𝑛 é o nível fractal. Pode-se observar ainda, na Figura 26, que em termos de área tem-se as probabilidades 𝑝1 e 𝑝2 no nível um. Onde para estrutura proposta tem-se: 𝑝1 = 1/17 e 𝑝2 = ⁡4/17. 42 Para calcular as demais probabilidades, deve-se usar a fórmula da distribuição binomial, já apresentada anteriormente, porém as probabilidades eram de segmentos, sendo aqui, de áreas: µ𝑘 = (𝑝1) 𝑘⁡(𝑝2)𝑛−𝑘 (7) em que 𝑛 é o nível da estrutura fractal e 𝑘 varia de zero a 𝑛. Figura 26 – Multifractal proposto 2/5 1/ 5 Fonte: imagem elaborada pelo autor, 2017. Como esse multifractal inicialmente tem duas similaridades lineares em termos de raios, no enésimo nível, para k variando de zero a n, tem-se: 𝑠𝑘 = (𝑠1) 𝑘(𝑠2) 𝑛−𝑘 (8) Uma vez calculadas as similaridades em termos de raios, calcula-se as dimensões locais do conjunto por: 𝑑𝑘 = log(µ𝑘) / log(𝑠𝑘). (9) Em seguida deve-se calcular o número total de similaridades para o nível 𝑛, conforme a seguir: #𝑠𝑘 = 4 𝑘 ( 𝑛 𝑘 ) = 4𝑘( 𝑛! 𝑝!(𝑛−𝑘)! ) (10) onde 𝑛 é o nível fractal e 𝑘 varia de 0 a 𝑛 para cada nível. Uma vez tendo as várias quantidades de elementos pode-se calcular o espectro multifractal usando: 43 𝑓𝑘 = −𝑙𝑜𝑔(#𝑠𝑘)/𝑙𝑜𝑔𝑠𝑘. (11) A Figura 27 mostra o gráfico do espectro multifractal DCB com 50 pontos, ou seja, com n = 50, cujas similaridades lineares em termos de raio é 1/5 e 2/5. O gráfico sinaliza que o ponto máximo no eixo vertical tem valor de 1,567. O valor usando, Hausdorff [25] é 1,59 cuja diferença é aproximadamente dois centésimos. É possível mostrar que quanto maior a quantidade de pontos, esse valor tende a dimensão de 1,59 . Esse fato é importante, pois é uma maneira de comparar a dimensão Hausdorff com o método do espectro multifractal. Para um cálculo detalhado verificar Apêndice I. Figura 27: Espectro multifractal DCB. Fonte: imagem elaborada pelo autor, 2017. 44 4.4 APLICAÇÕES DE MULTIFRACTAIS Existem poucas aplicações de multifractais em problemas de Eletromagnetismo. Assim, em [31] os autores propõem uma técnica de sinais eletromiográficos aplicados na musculatura humana, onde surgem padrões multifractais e, de acordo com os autores, é possível classificar como musculatura fadigada e não fadigada usando um modelo multifractal. Uma junção da curva de Koch com o tapete de Sierpinski foi proposta em [32]. Os autores conseguiram uma estrutura cuja principal vantagem é projetar estruturas com resposta multibanda, tornando invariável em termos de controle de banda e frequência de ressonância. A propriedade de autossimilaridade de uma geometria de Cantor multifractal modificada foi explorada em [33] para projetar superfícies seletivas em frequência com resposta multibanda. A principal vantagem da estrutura proposta é projetar FSS multibanda, com relações de frequência múltiplas entre bandas adjacentes e estruturas de fácil construção. Uma nova geometria multifractal foi proposta em [34] e recebeu o nome de Ericampos. A geometria se mostrou de fácil construção e apresentou resposta dual- band, com estabilidade angular e insensibilidade de polarização. O artigo escrito por Li,D and J. Mao na revista Electromagnetic waves [1070- 4698] yr:2012 vol:130 pg:207 com o título Sierpinskized Koch-like sided multifractal dipole propõe uma antena com elementos multifractais e mostra suas vantagens em relação as antenas monofractais como: multibanda com taxas diferentes, menor tamanho e melhor diretividade. 45 CAPÍTULO 5 RESULTADOS 5.1 RESPOSTAS EM FREQUÊNCIA DE GEOMETRIAS TRADICIONAIS A estrutura proposta, ilustrada na Figura 32, foi simulada visando a análise do efeito da periodicidade da geometria na resposta em frequência da FSS. Para isso, foram analisadas as respostas em frequência de geometrias tradicionais, amplamente publicadas na literatura. Foram considerados três tipos de geometrias muito usados: A espira quadrada [35], o dipolo cruzado [36] e o patch quadrado [37]. Em todos os casos, a periodicidade, p, foi variada de 14 a 20 mm. Cada geometria pertence a um diferente grupo, dentre alguns dos classificados pelo Munk [38]. As geometrias estão ilustradas na Figura 28. Figura 28 – Medidas das geometrias:(a) espira quadrada, (b) dipolo cruzado e (c) patch quadrado Fonte: imagem elaborada pelo autor, 2017. A primeira geometria pertence ao grupo das espiras foi a espira quadrada. As medidas estão ilustradas na Figura 28(a). A espira tem uma largura, d, de 11 mm, com espessura de fita, w, de 1,6 mm. Na Figura 29 é ilustrada a resposta em frequência da potência transmitida. Pode-se observar que não apenas a largura de banda varia, mas a frequência de ressonância também. Isso significa que as dimensões precisam 46 ser reprojetadas, caso se deseje alterar a periodicidade e manter a frequência de ressonância. Para p = 14 mm a frequência de ressonância ocorre em 5,64 GHz e a largura de banda foi de 1,7 GHz. Para p = 20mm a frequência de ressonância foi de 6,1 GHz e a largura de banda foi de 900 MHz. Figura 29 – Resposta em frequência da transmissão para a FSS com espira quadrada e diferentes valores de p Fonte: imagem elaborada pelo autor, 2017. A segunda geometria pertence ao grupo dos n-polos conectados pelo centro foi o dipolo cruzado. As medidas são ilustradas na Fig. 28(b). O dipolo cruzado tem comprimento, d, de 11 mm, com espessura de fita, w, de 4 mm. Na Figura 30 ilustramos a resposta em frequência da potência transmitida. Pode-se observar que, como no caso da espira quadrada, não apenas a largura de banda varia, mas a frequência de ressonância também. Entretanto, nesse caso a variação da frequência de ressonância foi mais crítica, se comparada ao caso da espira quadrada. Para p = 14 mm a frequência de ressonância ocorre em 10,24 GHz e a largura de banda foi de 2,1 GHz. Para p = 20mm a frequência de ressonância foi de 9,1 GHz e a largura de banda foi de 800 MHz. 47 Figura 30 – Resposta em frequência da transmissão para a FSS com dipolo cruzado e diferentes valores de p Fonte: imagem elaborada pelo autor, 2017. A terceira geometria pertence ao grupo dos elementos de interior sólido foi o patch quadrado. As medidas são ilustradas na Fig. 28(c). O patch tem uma largura, d, de 11 mm. Na Figura 31 ilustramos a resposta em frequência da potência transmitida. Pode- se observar que não apenas a largura de banda varia, mas a frequência de ressonância também. Isso significa que as dimensões precisam ser reprojetadas, caso se deseje alterar a periodicidade e manter a frequência de ressonância. Para p = 14 mm, a frequência de ressonância ocorre em 15,3 GHz e a largura de banda foi de 7,8 GHz. Para p = 18mm, a frequência de ressonância foi de 11,1 GHz e a largura de banda foi de 2GHz. 48 Figura 31 – Resposta em frequência da transmissão para a FSS com patch quadrado e diferentes valores de p Fonte: imagem elaborada pelo autor, 2017. 5.2 GEOMETRIA PROPOSTA E SUA RESPOSTA EM FREQUÊNCIA Para tentar reduzir essa forte dependência da frequência de ressonância com a periodicidade, nós propomos uma geometria fractal, com dupla razão de similaridade. A Figura 32 ilustra a geometria fractal com dupla razão de similaridade proposta. Os autores chegaram a essa estrutura a partir da ideia de propor uma geometria fractal com dupla razão de similaridade, diferente das geometrias existentes na literatura; e baseando-se em conjuntos fractais oriundos de círculos como o tetracírculo, análise de monofractais tipo Dürer e o multifractal Ericampos. 49 Figura 32 – Geometria fractal com dupla razão de similaridade: (a) iniciador e (b) gerador Fonte: imagem elaborada pelo autor, 2017. A partir da geometria proposta, nós efetuamos uma análise paramétrica, para diferentes valores de raios, r1 e r2, que geram diferentes razões de similaridade, e diferentes periodicidades, com o objetivo de obter pouca ou nenhuma sensibilidade de fr com a periodicidade p. Na primeira estrutura, o elemento gerador tem raio de 4,5 mm. As similaridades usadas foram 𝑠1= 0,22 e 𝑠2 = 0,39. A Figura 33 ilustra a resposta em frequência do coeficiente de transmissão em função da frequência para quatro valores de p. Analisando o gráfico, percebe-se que a variação da periodicidade de 10 mm a 13 mm não altera a frequência de ressonância, mas a largura de banda de frequência diminui com o aumento da periodicidade. 50 Figura 33 – Resposta em frequência da transmissão para a FSS com fractal proposto com razão de similaridade 𝒔𝟏 = 0,22 e 𝒔𝟐= 0,39 Fonte: imagem elaborada pelo autor, 2017. Na segunda estrutura, o elemento gerador tem raio de 5,5 mm. As similaridades usadas foram 𝑠1= 0,32 e 𝑠2= 0,34. A Figura 34 ilustra a resposta em frequência do coeficiente de transmissão em função da frequência para quatro valores de p. Usamos quatro valores diferentes de p, daqueles usados na estrutura anterior. Mais uma vez, percebe-se que a variação da periodicidade de 12 mm a 15 mm não altera a frequência de ressonância, enquanto que a largura de banda apresenta o mesmo comportamento da estrutura anterior. Na terceira estrutura, o elemento gerador tem raio de 6,5 mm. As similaridades usadas foram 𝑠1= 0,24 e 𝑠2= 0,38. A Figura 35 ilustra a resposta em frequência do coeficiente de transmissão em função da frequência para quatro valores de p. Usamos quatro valores diferentes daqueles usados nas estruturas anteriores. Entretanto, percebe-se o mesmo comportamento das estruturas anteriores, tanto para fr quanto para BW. A variação da periodicidade de 17 mm a 20 mm. 51 Figura 34 – Resposta em frequência da transmissão para a FSS com fractal proposto com razão de similaridade 𝒔𝟏 = 0,32 e 𝒔𝟐 = 0,34 Fonte: imagem elaborada pelo autor, 2017. Figura 35 – Resposta em frequência da transmissão para a FSS com fractal proposto com razão de similaridade 𝒔𝟏= 0,24 e 𝒔𝟐= 0,38. Fonte: imagem elaborada pelo autor, 2017. Após a análise paramétrica, nós chegamos a um elemento gerador que tem raio de 6,5 mm. As razões de similaridade usadas foram 𝑠1= 0,246 e 𝑠2= 0,377. Simulamos a resposta em frequência da transmissão dessa estrutura, para periodicidade, p, variando de 14 a 20 mm. Na Figura 36 ilustramos a resposta da transmissão em função da frequência. Para os diferentes valores de periodicidade testados, a frequência de ressonância ocorre em, aproximadamente, 6,00 GHz. Com relação à 52 largura de banda, foi de 2,15 GHz, para p = 14mm, de 1,18 GHz para p = 16 mm, de 1,01 GHz para p = 18 mm, e de 825 MHz, para p = 20mm. Figura 36 – Resposta em frequência da transmissão para a FSS com fractal proposto com razão de similaridade 𝒔𝟏= 0,246 e 𝒔𝟐= 0,377. Fonte: imagem elaborada pelo autor, 2017. Para comprovar as simulações efetuadas, dois protótipos foram construídos e medições foram realizadas. A Figura 37 ilustra os protótipos construídos. O setup de medição pode ser visto na Figura 38. Ele é composto por um analisador vetorial de rede de duas portas Agilent E5071C, dois pares de antenas tipo corneta (SAS-571 de 0,7 - 18,0 GHz) e um anteparo para reduzir as influências do ambiente de medições. 53 Figura 37 – Protótipos construídos Fonte: imagem elaborada pelo autor, 2017. Figura 38 – Setup de medição Fonte: imagem elaborada pelo autor, 2017. O primeiro protótipo construído tem periodicidade de 14 mm nas direções x e y. Como consequência, 196 células foram impressas sobre um dielétrico de fenolite com (200 x 200) mm2 de área. A Figura 39 ilustra a comparação entre os resultados simulados e medidos. A boa concordância entre os dois resultados foram obtidos, ou seja, o simulado e o medido. 54 Figura 39 – Comparação entre resultados simulados e medidos para p = 14 mm. Fonte: imagem elaborada pelo autor, 2017. O segundo protótipo construído tem periodicidade de 18 mm nas direções x e y. Como consequência, 121 células foram impressas sobre o mesmo substrato do protótipo anterior. A Figura 40 ilustra a comparação entre os resultados. Novamente, a boa concordância entre os resultados foi obtido. Figura 40 – Comparação entre resultados simulados e medidos para p = 18 mm Fonte: imagem elaborada pelo autor, 2017. Por fim, geraram-se as curvas experimentais dos dois protótipos juntos e pode-se ver que não há variação na frequência de ressonância, mas a largura de banda varia de 2,46 GHz, com p = 14 mm, para 1,59 GHz com p = 18 mm. A Figura 41 ilustra esse comportamento: 55 Figura 41 – Resposta em frequência experimental da transmissão para as FSS construídas Fonte: imagem elaborada pelo autor, 2017. 56 CAPÍTULO 6 CONCLUSÕES Neste trabalho, foram apresentadas as definições de FSS para um embasamento de seu funcionamento, suas características, algumas técnicas de análise e formas geométricas mais comuns utilizadas nessas estruturas. Em seguida foram apresentados os fractais, alguns tipos de construções, análise de dimensão de conjuntos fractais e outros tópicos como fundamentos para as geometrias multifractais. Já no Capítulo 4, foram apresentados os multifractais, suas construções e cálculos do espetro de espectros multifractais que mostram o quanto essas estruturas são irregulares, e neste, foi proposta uma nova geometria multifractal para o projeto de FSS com pouca insensibilidade da frequência de ressonância, com relação à periodicidade da célula unitária. A nova geometria foi criada a partir de estudos das geometrias tradicionais fractais, como o fractal tetracírculo, fractal hexagonal de Dürer e o multifractal Ericampos, mas usando círculos com duas razões de similaridade. A FSS apresentou a vantagem de não precisar redimensionar suas medidas físicas para alterar a largura de banda, bastando alterar o afastamento (periodicidade) entre as células. A principal vantagem da estrutura proposta é a facilidade de ajustar a largura de banda, modificando apenas a periodicidade das células. A validação da estrutura proposta foi inicialmente verificada através de simulações em um software comercial e depois com medições apresentadas no capítulo de resultados. E conforme os gráficos apresentados foi obtida uma boa aproximação entre os resultados simulados e os medidos. As pequenas diferenças observadas são devidos aos processos de medidas, aos cortes imprecisos da estrutura e a resolução da ferramenta de corte na confecção do patch. Para continuidade deste trabalho sugere-se a análise dessa estrutura com relação à estabilidade angular e ao comportamento em outros níveis fractais. Sugere-se também pesquisar outras estruturas semelhantes para uma análise paramétrica detalhada. Por fim, sugere-se a construção da estrutura proposta como patch em antenas para averiguar se ocorrem características semelhantes às encontradas neste trabalho. 57 REFERÊNCIAS [1] T.K Wu, “Fundamentals of peridic strutures”, in Frequency Selective Surface and Grid Array, 1nd ed. vol. 1, E. U. A, John Wiley & Sons,INC., 1995. [2] A. L. P. S. Campos, “Formas de elementos”, superfícies seletivas em frequência análise e projeto, 1ed. Vol. 1,Brasil, editora IFRN, 2008) [3] E.C Braz, Tese de doutorado. Superfícies Seletivas em frequência com geometrias multifractais, Universidade federal do rio grande do norte, 2014. [4] S. W. Lee, G. Zarrillo, and C. L. Law, “Simple formulas for transmission through periodic metal grids or plates”, IEEE Transactions on Antennas and Propagation, vol. 30, no. 5, pp. 904-909, 1982. [5] G. Zarrilo and K. Aguiar, “Closed-Form low frequency solutions for electromagnetic waves through a frequency selective surface”, IEEE Transactions on Antennas and Propagation, vol. 35, no. 12, pp. 1406-1417, 1987. [6] R. J. Langley and E. A. Parker, “Equivalent circuit model for arrays of square loops”, Electronics Letters, vol. 18, no. 7, pp. 294-296, 1982. [7] R. J. Langley and E. A. Parker, “Double-square frequency selective surfaces and their equivalent circuit”, Electronic Letters, vol. 19, no. 17, pp. 675-677, 1983. [8] C. K. Lee and R. J. Langley, “Equivalent circuit models for frequency selective surfaces at oblique angles of incidence”, IEE Proceedings, vol. 132, no. 6, pp. 395- 399, 1985. [9] J. Huang, T. K. Wu, and S. W. Lee, “Tri-band FSS with circular ring elements”, IEEE Transactions on Antennas and Propagation, vol. 42, no. 2, pp. 166-175, 1994. [10] P. Harms, R. Mittra, and W. Ko, “Implementation of the periodic boundary condition in the finite-difference time-domain algorithm for FSS structures”, IEEE Transactions on Antennas and Propagation, vol. 42, no. 9, pp. 1317-1324, 1994. [11] D. S. Weile, E. Michielssen, and K. Gallivan, “Reduced order modeling of multiscreen frequency selective surfaces using Krylov–based rational interpolation”, IEEE Transactions on Antennas and Propagation, vol. 49, no. 5, pp. 801-813, 2001. 58 [12] D. S. Weile and E. Michielssen, “Analysis of frequency selective surfaces through the blazing onset using Krylov model order reduction and Woodbury singularity extraction”, IEEE Transactions on Antennas and Propagation, vol. 49, no. 10, pp. 1470-1478, 2001. [13] M. Bozzi, L. Perregrini, J. Weinzierl, and C. Winnewisser, “Efficient analysis of quase-optical filters by a hybrid MoM/BI-RME method”, IEEE Transactions on Antennas and Propagation, vol. 49, no. 7, pp. 1054-1064, 2001. [14] G. Pelosi, A. Cocchi, and S. Séller, “Electromagnetic scattering from infinite periodic structures with a localized impurity”, IEEE Transactions on Antennas and Propagation, vol. 49, no. 5, pp. 697-702, 2001. [15] G. Manara, A. Monorchio, and R. Mittra, “Frequency selective surfaces desing based on genetic algorithm”, Electronic Letters, vol. 35, no. 17, pp. 1400-1401, 1999. [16] S. Chakravarty, R. Mittra, and N. R. Williams, “On the application of the microgenetic algorithm to the design of broad band microwave absorbers comprising frequency selective surfaces embedded on multilayered dielectric media”, IEEE Transactions on Microwave Theory and Techniques, vol. 49, no. 6, pp. 1050-1059, 2001. [17] T. Itoh, “Spectral domain immitance approach for dispersion characteristics of generalized printed transmission lines”, IEEE Transactions on Antennas and Propagation, vol. 28, pp. 733-736, 1980. [18] T. Itoh e R. Mittra, “Spectral domain approach for calculating the dispersion characteristics of microstriplines”, IEEE Transactions on Microwave Theory and Techniques, vol. 21, pp. 496-499, 1973. [19] A. L. P. S. Campos, A. G. d’Assunção e M. A. B. Melo, “Frequency selective surfaces with anisotropic dielectric substrates”, International Journal of Infrared and Millimeter Waves, vol. 21, no. 3, pp. 461-475, 2000. 59 [20] T. L. Silva, “Análise de onda completa de superfície seletiva de frequência do tipo anteparo duplo”, dissertação de mestrado, Universidade Federal de Rio Grande do Norte, 2011 [21] A. L. P. S. Campos, M. A. B. Melo e A. G. d’Assunção, “Frequency selective surfaces with rectangular apertures on uniaxial anisotropic substrates”, Microwave and Optical Technology Letters, vol. 25, no. 2, pp. 126-129, 2000. [22] H. P. A, Ferreira, “Matrioska: Uma Proposta de geometria para fss multibandas”, Dissertação de mestrado, Instituto federal de educação, ciência e tecnologia da paraíba, 2014. [23] Eduardo Moreira. [online]. Avaliable:http://www.techtudo.com.br/noticias/noticia/2012/05/papel-de-parede-e- capaz-de-bloquear-sinal-da-rede-wi-fi.html. Accessed on: mar., 14, 2016. [24] Technology Review. [online]. Avaliable: http://www.technologyreview.com.br/printer_friendly_article.aspx?id=48498. Accessed on: mar., 16, 2017. [25] K, Falconer, ”Introduction”, Fractal Geometry Mathematical Foundations and Applications, 2a ed. vol. 1, John Willey and Sons, New York, 2003 [26] Dicionário do Aurélio. [online]. Avaliable: https://dicionariodoaurelio.com/fractal. Accessed on: mar., 12, 2017. [27] B. B. Madebroit, The fractal geometry of nature, 1983. Printed in the United States of America. [28] UOL Educação, Biografia de Gerg Cantor. [online]. Avaliable: https://educacao.uol.com.br/biografias/georg-cantor.htm. Accessed on: mar., 16, 2017. [29] D. Harte, Multifractals Theory and Applications, Chapman and Hall/CRC, May 2001. [30] Prisma: A luz da Física. Aplicações de fractais. [online]. Avaliable: http://cftc.cii.fc.ul.pt/PRISMA/capitulos/capitulo2/modulo4/topico9.php. Accessed on: mar., 18, 2017. 60 [31] Chakraborty, B. et al., “Multifractal Approach for Seafloor Characterization”, Geoscience and Remote Sensing Letters, vol. 11, no 1, pp. 54 – 58, 2014. [32] E. C. Braz e Campos, A. L. P. S. Campos, Dual/wide band multifractal frequency selective surface for applications in S- and X-band. Microwave and Optical Technology Letters, vol. 56, pp. 2217-2222, 2014. [33] E. C. Braz and Campos, A. L. P. S. Campos Multiband frequency selective surfaces with a modified multifractal cantor geometry. Journal of Microwaves, Optoelectronics and Electromagnetic Applications, vol. 13, pp. 111-121, 2014. [34] E. C. Braz; Campos, A. L. P. S., A New Multifractal Geometry for Design of Frequency Selective Surfaces with Dual Band Response. Applied Computational Electromagnetics Society Journal, vol. 30, pp. 76-81, 2015. [35] A. E. Yilmaz, M. Kuzuoglu. “Design of the Square Loop Frequency Selective Surfaces with Particle Swarm Optimization via the Equivalent Circuit Model”.Radioengineering, vol. 18, no. 2, pp. 95 – 102, 2009. [36] E. P. Santos et al., “Optimization of Frequency Selective Surfaces with Crossed Dipoles Using the Taguchi's Method”, Journal of Microwaves, Optoelectronics and Electromagnetic Applications, vol. 16, no. 2, pp. 371 – 384, 2017. [37] A. L. P. S. Campos, A. G. d’Assunção, andL. M. de Mendonça,“Scattering by FSS on Anisotropic Substrate forTE and TM Excitation”,IEEE Transactions on Microwave Theory and Techniques, vol. 50, no. 1, pp. 72 – 76, 2002. [38] B. A. Munk, “ Elements Types: A comparizon”, Frequency Selective Surfaces: Theory and Design, John Wiley & Sons, Nova Iorque, 2000. [39] R. A. Zanotto Estudo de Geometria Fractal Clássica, dissertação de mestrado, Universidade Federal de Goiás , 2015. 61 APÊNDICE I MULTIFRACTAL DE CANTOR CÁLCULO DA DIMENSÃO DETALHADA Usando o nível dois desse multifractal como exemplo, o primeiro passo é reconhecer que as probabilidades para o nível um é 1/3 e 2/3. Em seguida, utiliza-se 𝑛 = 2 e 𝑘 variando de 0 a 2. Assim, para o primeiro segmento à esquerda e os dois segmentos seguintes, aplica-se a equação (3). µ2 = ( 1 3 ) 2 ( 2 3 ) 2−2 = 1/9. µ1 = ( 1 3 ) 1 ( 2 3 ) 2−1 = 2 9 . µ0 = ( 1 3 ) 0 ∗ ( 2 3 ) 2−0 = 4 9 . O segundo passo é necessário calcular a quantidade de segmentos autossimilares usando análise combinatória conforme segue: Variando, novamente 𝑘 de 0 à 2 para os segmentos de probabilidades 1/9, 2/9 e 4/9, usa-se a equação (4). ( 𝑛 𝑘 ) = ( 2 2 ) = 1 ( 𝑛 𝑘 ) = ( 2 1 ) = 2 ( 𝑛 𝑘 ) = ( 2 0 ) = 1 Deve-se observar que os segmentos têm três probabilidades, 1, 2 e 1, ou seja, duas das probabilidades são iguais. 62 O passo posterior é calcular o comportamento local de cada barra ou segmento aplicando-se a equação (5) Para primeira barra, segunda e terceira, tem-se: 𝑑2 = log ( 1 9 ) log⁡( 1 9 ) 𝑑2 = 1 𝑑1 = log ( 2 9 ) log ( 1 9 ) 𝑑1 = 0.68 𝑑0 = log ( 4 9 ) log⁡( 1 9 ) 𝑑0 = − log3 ( 4 9 ) 2 𝑑0 = 0.36 Por último, calcula-se o espectro multifractal pela equação (6). Para a primeira barra, segunda e terceira, tem-se: 𝑓2 = log(1) −log⁡( 1 9 ) 𝑓2 = 0 63 𝑓1 = log(2) −log⁡( 1 9 ) 𝑓1 = log3 2 2 𝑓1 = 0.31 𝑓0 = log(1) −log⁡( 1 9 ) 𝑓0 = 0 MULTIFRACTAL DCB CÁLCULO DA DIMENSÃO DETALHADA. Para esse exemplo, usaremos também o nível dois desse multifractal, entretanto, conforme já mencionado, sendo as similaridades de raios diferentes, surge uma pequena diferença que será vista no exemplo que segue. Conforme a Figura 26, suas similaridades em raios são: 1/5 e 2/5. Suas respetivas áreas são: (1/5)2𝜋 e (2/5)2𝜋 e as probabilidades em áreas são, 1/17 e 4/17. Logo, as probabilidades para o nível dois, fazendo 𝑘 variando de 0 à 2 e aplicando a equação (7), tem-se: µ2 =⁡( 1 17 ) 2 ( 4 17 ) 2−2 µ2 = 1/289 µ1 = (1/17) 1⁡(4/17)2−1 64 µ1 = (1/17) ∗ (4/17) = 4/289 µ0 = (𝑝1) 0⁡(𝑝2)2−0 =⁡ (1/17)0⁡(4/17)2−0 µ0 = 16/289 No passo seguinte, calcula-se as similaridades em raios a partir da equação (8) repetida a seguir: 𝑠𝑘 = (𝑠1) 𝑘(𝑠2) 𝑛−𝑘 Para 𝑘 variando de 0 a 2, tem-se: 𝑠0 = ( 1 5 ) 0 ( 2 5 ) 2−0 𝑠0 = 4/25 𝑠1 = ( 1 5 ) 1 ( 2 5 ) 2−1 𝑠1 = 2/25 𝑠2 = ( 1 5 ) 2 ( 2 5 ) 2−2 𝑠2 = 1/25 No terceiro passo, acha-se o comportamento local para o nível 2, usando-se a equação (9) conforme segue: 𝑑𝑘 = log(µ𝑘) / log(𝑠𝑘). 65 𝑑0 = log(16/289) / log(4/25) = 1,57 𝑑1 = log(4/289) / log(2/25) =1,67 𝑑2 = log ( 1 289 ) / log ( 1 25 ) =1,76 Para o cálculo do número total de similaridades para o nível 2, do DBC aplica-se a equação (10). #𝑠0 = 4 0 ( 2 0 ) = 1 #𝑠1 = 4 1 ( 2 1 ) = 8 #𝑠2 = 4 2 ( 2 2 ) = 16 E o último passo, tem-se o cálculo do espectro multifractal, como segue a equação (11) abaixo transcrita: 𝑓𝑘 = −𝑙𝑜𝑔(#𝑠𝑘)/𝑙𝑜𝑔𝑠𝑘. Para 𝑘 variando também de 0 a 2, tem-se: 𝑓0 = − 𝑙𝑜𝑔1 log ( 4 25 ) = 0 𝑓1 = − 𝑙𝑜𝑔8 log ( 2 25 ) = 0,82 𝑓2 = − 𝑙𝑜𝑔1 log ( 1 25 ) = 0,86 66 APÊNDICE II Um ponto importante para construção de fractais tipo Dürer é a porcentagem de redução do lado iniciador para o gerador. Para determinar a porcentagem de redução dos lados, observe a Figura 42. Primeiramente, calcula-se o valor do ângulo interno do polígono em questão, por se tratar de polígonos regulares, usamos a fórmula da razão da soma dos ângulos internos pelo número de lados do polígono:𝑎𝑖 ⁡= 180(𝑛⁡−⁡2) 𝑛 onde n é o número de lados do polígono, assim, para o hexágono tem-se: 𝑎𝑖 = 180(6⁡−⁡2) 6 ⁡= ⁡120. Em seguida, usando o triângulo equilátero central formado em cada figura, estabelece-se os ângulos internos desse triângulo; Sabendo-se que a soma interna de um triângulo equilátero é 180. Dividindo 180 por 3, resulta em 60º. x = 2a + b. Para calcular a porcentagem de redução para o Dürer hexagonal, usando o triângulo equilátero formado nesse polígono, aplica-se a lei dos senos, como segue: 𝑠𝑖𝑛60𝑥⁡𝑎 = 𝑠𝑖𝑛60⁡𝑥𝑏 que resulta, 𝑎 = 𝑏. Aplicando em 𝑥 = 2𝑎 +b, fica 𝑥 = 2𝑎⁡ + 𝑎 = ⁡𝑎 = 𝑥/3. Portanto, conclui-se que o lado, em cada iteração, reduz a uma taxa de 1/3, ou 33,33%. De maneira análoga, para calcular a porcentagem de redução para os fractais Dürer com mais de quatro lados, deve-se calcular o ângulo interno do polígono pela fórmula já apresentada. Em seguida, percebe-se que as bases dos triângulos ficam sempre adjacentes aos lados do polígono formando 180. Para achar o ângulo da base, subtrai-se (180⁡ − ⁡𝑎𝑖). Para calcular o outro ângulo, subtrai-se (180 − 2(⁡180 − 𝑎𝑖⁡)). A seguir, usa-se a lei dos senos e por fim substitui-se a relação achada em x=2a+b para encontrar o fator de redução da primeira iteração. Usando um pouco de álgebra, chega-se: 𝑎 = 𝑥𝑠𝑖𝑛(180 − 𝑎𝑖) 2 + 𝑠𝑖𝑛𝛼 𝑎 = 𝑥𝑠𝑖𝑛(180 − 180(𝑛 − 2)/𝑛) 2 + 𝑠𝑖𝑛(180 − 2 (⁡180 − 180(𝑛−2) 𝑛 ))) Onde: 𝑎 é o fator de redução. 𝑥 é o lado do polígono iniciador. 67 (180 - ai) é o angulo adjacente. 𝛼 é o ângulo do vértice superior. b é a base do triângulo. 𝛼 é igual a (180-2( 180-ai )) como descrito anteriormente. Para uma outra descrição de como calcular o fator de redução, segue [39]. Figura 42 – Ilustração para dedução do fator de redução. xa b Fonte: imagem elaborada pelo autor, 2017.