Resumos
O fenômeno de caos é investigado em um experimento relacionado com circuitos RLC forçados, com o objetivo de estudar conceitos básicos e possíveis aplicações práticas. Neste circuito, a substituição do capacitor por um diodo, cuja capacitância varia não-linearmente com a sua voltagem, pode levar a oscilações caóticas se a tensão aplicada ao circuito oscilar com a freqüência nescessária. Programas de apoio são utilizados de maneira que os estudantes possam simular as trajetórias de caos.
caos; bifurcação; mapa logístico
The chaos phenomenon is investigated in an experiment related to RLC circuits with the aim of studying the basic concepts and possible practical applications. In this circuit, the substitution of the capacitor by the diode provides a variable capacitance leading to chaotic oscillations if the frequency of the applied tension is adequate. Support softwares are used in such a way that the students can simulate the trajectories leading to chaos.
chaos; bifurcation; logistic map
ARTIGOS GERAIS
Comportamento caótico em um circuito RLC não-linear
Chaotic behavior in a non-linear RLC circuit
N. Carlin1; E.M. Szanto; W.A. Seale; F.O. Jorge; F.A. Souza; I.H. Bechtold; L.R. Gasques
Instituto de Física, Universidade de São Paulo, São Paulo, SP, Brasil
RESUMO
O fenômeno de caos é investigado em um experimento relacionado com circuitos RLC forçados, com o objetivo de estudar conceitos básicos e possíveis aplicações práticas. Neste circuito, a substituição do capacitor por um diodo, cuja capacitância varia não-linearmente com a sua voltagem, pode levar a oscilações caóticas se a tensão aplicada ao circuito oscilar com a freqüência nescessária. Programas de apoio são utilizados de maneira que os estudantes possam simular as trajetórias de caos.
Palavras-chave: caos, bifurcação, mapa logístico.
ABSTRACT
The chaos phenomenon is investigated in an experiment related to RLC circuits with the aim of studying the basic concepts and possible practical applications. In this circuit, the substitution of the capacitor by the diode provides a variable capacitance leading to chaotic oscillations if the frequency of the applied tension is adequate. Support softwares are used in such a way that the students can simulate the trajectories leading to chaos.
Keywords: chaos, bifurcation, logistic map.
1. Introdução
Historicamente, o estudo de caos teve início em Matemática e Física, sendo que mais tarde foi estendido para outras áreas como Engenharia, Informática, Ciências Sociais, Biologia e Medicina [1]. Nos últimos anos, tem aumentado o interesse em aplicações comerciais e industriais. Como exemplo pode-se mencionar o uso do caos em processamento de informações, análises financeiras e predição e controle de atividade cardíaca [2, 3].
O presente experimento é parte de uma série de atividades laboratoriais relacionadas a eletricidade, magnetismo e óptica para estudantes de segundo e terceiro anos do curso de bacharelado no Instituto de Física e Instituto de Astronomia, Geofísica e Ciências Atmosféricas, ambos situados na Universidade de São Paulo.
O desenvolvimento e interpretação de experimentos é a base inicial do curso, juntamente com um estudo formal dos fenômenos em questão e uma introdução de aplicações com abordagem mais quantitativa. A infraestrutura adequada foi preparada, com computadores acoplados a interfaces para aquisição de dados, câmeras CCD para gravação de imagens e o softwares para realização de cálculos, simulações e gráficos em grupos de dois a três estudantes.
Este experimento é realizado em três aulas, cada uma com quatro horas de duração. A primeira aula é dedicada para um estudo teórico básico de ressonância em circuitos RLC e na segunda e terceira aulas são introduzidos os conceitos básicos de caos, sendo que neste artigo iremos nos concentrar na parte relacionada com caos. Ao final do experimento espera-se que os estudantes tenham adquirido algumas habilidades para reconhecer comportamentos caóticos, bem como alguns conhecimentos sobre sistemas não-lineares e suas equações diferenciais.
No circuito RLC o capacitor é substituído por um diodo, o que gera uma variação não-linear da capacitância em função da voltagem aplicada, o que pode resultar no aparecimento do caos. Os estudantes fazem observações experimentais da duplicação de período, que é uma das possíveis rotas para o caos, e coletam dados para construção do diagrama de bifurcações, sendo que deste diagrama é possível calcular o número de Feigenbaum [4]. O retrato de fase deste circuito é observado e comparado com o de um circuito RLC normal e simulações computacionais são realizadas para fixar a compreensão.
2. Base teórica
Os conceitos básicos e nomenclaturas necessárias para o entendimento deste experimento por parte dos estudantes são expostos de uma maneira qualitativa. Um sistema dinâmico está associado a uma lei de evolução, de forma que a partir de um conjunto de variáveis que descrevem a configuração num dado instante pode-se determinar valores dessas variáveis em instantes posteriores, e isso pode ser feito através de um conjunto de equações diferenciais. Um exemplo simples seria considerar as leis de Newton, as quais permitem a determinação de velocidades e posições de um conjunto de partículas a partir dos seus valores iniciais. Se essas equações diferenciais forem não-lineares, o sistema que elas descrevem é dito não-linear e pode, eventualmente, ter evoluções caóticas.
Neste trabalho iremos estudar as equações diferenciais que descrevem o circuito utilizado, que consiste de um gerador de função, um resistor R, um indutor L e um diodo D conectados em série. Devido à presença do diodo, o circuito se comporta como um oscilador forçado com uma resposta não-linear. A relação entre corrente e voltagem em um diodo foi estudada em um experimento anterior dos estudantes. Devido às propriedades da junção p-n, o diodo apresenta uma capacitância não-linear (C), descrita por [5]:
onde VD é a voltagem no diodo, e é a carga do elétron, T é a temperatura em Kelvin e k = 1,38 x 10-23 J/K é a constante de Boltzmann. Portanto, o diodo pode ser modelado como uma capacitância dependente da voltagem em paralelo com um diodo ideal (ver lado direito da Fig. 3).
Considerando a carga q e a corrente i como sendo as variáveis dinâmicas, as equações que descrevem um circuito RLC com uma voltagem aplicada V0senwt são:
Para o circuito deste experimento, com a troca do capacitor por um diodo, as equações dinâmicas podem ser escritas como:
onde f(q) e g(q) levam em conta o comportamento não-linear do diodo.
Em geral, a evolução de um sistema físico depende de um ou mais parâmetros, denominados de parâmetros de controle. O comportamento dinâmico de um sistema pode ser bastante diferente se os parâmetros de controle forem alterados. Por exemplo, no caso de um circuito RLC, tanto a freqüência como a voltagem aplicada podem ser considerados parâmetros de controle.
Quando um parâmetro de controle µ é variado, podem aparecer novos padrões de comportamento ou seqüências de novos estados para o sistema. Neste caso, diz-se que ocorreram bifurcações e µj é o valor do parâmetro de controle para o qual ocorre cada bifurcação.
O tipo mais comum de bifurcação é a chamada bifurcação de duplicação de período. Nesse caso, uma solução estável do sistema perde a estabilidade com a variação de um parâmetro de controle e aparece ao mesmo tempo uma nova solução estável com o dobro do período da solução anterior. Se o parâmetro de controle for variado continuamente, esse processo pode se repetir indefinidamente levando a soluções com período infinito, o que geralmente leva ao caos. Chamando de µj o valor de µ na j-ésima bifurcação, o período da solução imediatamente após a bifurcacao será 2j. Os estados assintóticos do sistema são geralmente denominados de atratores e dizemos que o atrator sofre a duplicação de período. Geralmente, conforme varia-se o parâmetro de controle, regiões com soluções caóticas podem ser observadas entre regiões com atratores estáveis, que usualmente são chamadas de janelas de caos. Foi verificado que os valores de µ j, para os quais as bifurcações ocorrem, obedecem à seguinte lei de universalidade:
A constante d quantifica quanto decresce a diferença entre os valores do parâmetro de controle para duas bifurcações sucessivas e é conhecida como número de Fiegenbaum (d = 4,6692016091029...), sendo universal para sistemas que apresentam duplicação de período [4]. Experimentalmente, quanto maior o valor de j, mais próximo d estará do valor acima.
3. Experimento numérico
Em muitos casos, é possível descrever a evolução de um sistema físico por um conjunto de equações discretas ou mapas. O exemplo mais conhecido é o do mapa logístico [6], usado por exemplo na análise do controle de populações, utilizando a equação abaixo:
onde xn-1 e xn representam os estados inicial e posterior, e n = 1,2,3,... é o número de iterações. O comportamento assintótico do sistema é determinado pelo parâmetro de controle r. Para estudar a evolução de uma trajetória com a Eq. (5), é necessário conhecer a condição inicial x0 e o parâmetro r. Pode-se utilizar uma planilha de cálculo para gerar uma tabela de dados com os valores assintóticos de x em função de r, os quais podem ser posteriormente transformados em um gráfico. Esse gráfico é conhecido como diagrama de bifurcações.
Na Fig. 1(a) é mostrado o valor assintótico da coordenada x do sistema considerando-se a condição inicial x0 = 0,5 na Eq. (5). Para valores r = 3,25 e n > 100, o sistema oscila entre dois valores distintos de x (atrator de período 2) e para r = 3,52 observam-se quatro valores assintóticos de x, representando um atrator de período quatro (ver Figs. 1(b) e 1(c)). O diagrama de bifurcações montado com os valores calculados para 1,00 < r < 4,00 está apresentado na Fig. 1(d), onde no detalhe de 3,50 < r < 4,00 é possível observar as janelas de caos.
Neste mesmo contexto os estudantes foram estimulados a introduzir pequenas variações na condição inicial x0 para um valor fixo de r. Na Fig. 2 apresentamos simulações para r = 3,99 (onde a função já apresenta um regime caótico) e valores de x0 = 0,50000 e x0 = 0,50001. Com isso, os estudantes puderam perceber que a evolução do sistema é bastante sensível às condições iniciais, pois observa-se que com esta variação de x0 o sistema muda de comportamento a partir de 24 iterações. A importância prática disso é que o controle de um sistema caótico requer uma precisão grande no controle das condições iniciais do sistema, caso contrário a partir de um baixo número de iterações torna-se impossível prever o estado do sistema.
4. Procedimento experimental
O experimento é iniciado com a montagem do circuito esquematizado na Fig. 3, cujos componentes são: R = 10 W, L = 1 mH e um diodo 1N4007. Na Fig. 4(a) é apresentada uma fotografia do circuito montado com as ligações do osciloscópio e o gerador de tensão senoidal com freqüência variável. Nota-se que um dos canais do osciloscópio mede a tensão aplicada ao circuito, enquanto o outro canal mede a queda de tensão apenas no diodo. A capacitância de um diodo típico em função da voltagem aplicada é mostrada na Fig. 4(b). Neste experimento, utilizamos a freqüência como parâmetro de controle, com uma tensão aplicada constante de amplitude V0 = 2 V.
Inicialmente, os estudantes utilizaram os dois canais do osciloscópio para acompanhar a tensão aplicada no circuito e a tensão no diodo, variando-se a freqüência entre 40 e 200 kHz. Na Fig. 5(a) é possível observar o sinal dos dois canais do osciloscópio para diferentes freqüências, onde fica clara a evolução do sistema, apresentando duplicação do período e regiões de caos na tensão do diodo (o sinal de entrada do circuito é dado pela curva superior e o do diodo pela curva inferior nas imagens da Fig. 5(a)). O diagrama de bifurcações pode ser construído acompanhando-se a variação nas amplitudes dos picos presentes na resposta do diodo em função da freqüência. Na Fig. 5(b) apresentamos o diagrama de bifurcações obtido, através do qual observa-se as janelas de caos. Com as freqüências correspondentes às bifurcações, o número de Feigenbaum pode ser estimado para comparação com o valor conhecido, onde determinou-se que d = (µ4 - µ2)/(µ8 - µ4) = 4,7±0,2. A incerteza neste valor foi determinada através de propagação de erros e deve-se basicamente à incerteza na medida da freqüência de duplicação do período que foi de ±0,1 kHz. Entretanto, é preciso ressaltar que o número de Feigenbaum (d = 4,6692016091029...) é o limite obtido para essa razão após se observar um grande número de bifurcações.
As trajetórias para o caos também podem ser visualizadas no modo x-y do osciloscópio (x ® tensão de entrada e y ® tensão no diodo), como apresentado na Fig. 6, o que permite a observação da duplicação do período de maneira mais nítida, onde o número de períodos é dado pelo número de picos. Nesta figura, estão identificadas as freqüências aplicadas e observa-se que a primeira janela de caos ocorre para f = 57,0 kHz.
No estudo dos retratos de fase, a observação experimental das variáveis de interesse q e i é acompanhada no modo x-y do osciloscópio pela tensão no diodo (que é proporcional a q) e no resistor (que é proporcional a i). Como forma de comparação, os estudantes inicialmente observaram o retrato de fase de um circuito RLC linear, pela troca do diodo por um capacitor de 1 µF. Pelo conhecimento anterior de fasores, os estudantes facilmente entenderam que o padrão correspondente devia ser uma elipse cuja forma varia com a variação da freqüência aplicada, conforme Fig. 7(a). Então, o capacitor foi novamente substituído pelo diodo e um novo padrão era observado em função da freqüência, como mostrado na Fig. 7(b). Para complementar este estudo, os estudantes utilizaram o programa Maple [7] para simular os retratos de fase para os casos linear e não-linear considerados acima, utilizando as Eqs. (2) e (3), respectivamente. Os parâmetros das equações foram variados para melhor compreensão de sua influência no comportamento das rotas que levam para o caos. Na Fig. 8 são apresentados algumas simulações que coincidem com os resultados experimentais apresentados na Fig. 7.
5. Comentários finais
O objetivo do experimento apresentado é o estudo de circuitos RLC lineares e não-lineares. A não-linearidade introduzida por um diodo, gerando uma capacitância não-linear, possibilita a observação de caos. O experimento é realizado em três períodos de quatro horas cada, iniciando com uma descrição mais formal do fenômeno de ressonância em um circuito RLC linear, seguido por um estudo mais qualitativo da duplicação de períodos, bifurcação, caos e os retratos de fase correspondentes.
A parte experimental é acompanhada por simulações computacionais que auxiliam a compreenção por parte dos estudantes, onde vale destacar a utilização da equação logística (Eq. (5)) que simula perfeitamente os parâmetros importantes na rota para o caos. A partir do diagrama de bifurcações obtido experimentalmente foi calculado o número de Feigenbaum d = 4,7±0,2, que considerando-se o experimento em questão está bem próximo do valor esperado d = 4,6692016091029... Com este estudo, é possível os estudantes observarem e entenderem na prática o fenômeno de caos a partir de um circuito RLC.
A filosofia descrita é aplicada em experimentos para estudantes de segundo e terceiro anos do Instituto de Física e Instituto de Astronomia, Geofísica e Ciências Atmosféricas da Universidade de São Paulo. Os experimentos são relacionados à eletricidade, magnetismo e óptica. O contato dos estudantes com equipamentos e técnicas experimentais modernas e o aprendizado de variáveis e fenômenos das atividades do cotidiano aumenta significativamente seu conhecimento, motivação e criatividade.
Agradecimentos
Os autores agradecem ao Prof. Dr. Iberê Luiz Caldas pela leitura cuidadosa do manuscrito, pelas várias sugestões e discussões.
Recebido em 14/9/2004; Aceito em 16/12/2004
Referências bibliográficas
- [1] Nelson F.-Ferrara, Caos: Uma Introduçăo (Edgard Blucher, Săo Paulo, 1994).
- [2] T. Hogg and B.A. Huberman, Proceedings of the National Academy of Sciences, 81 (1984).
- [3] A.T. Winfree, The Geometry of Biological Time (Springer-Verlag, Nova Iorque, 1980).
- [4] M.J. Feigenbaum, J. Stat. Phys. 19, 25 (1978).
- [5] R.V. Buskirk and C. Jeffries, Phys. Rev. A 31, 3332 (1985).
- [6] R.M. May, Nature 261, 459 (1976).
- [7] Programa Maple 8, Maplesoft Inc., USA.
Datas de Publicação
-
Publicação nesta coleção
24 Abr 2006 -
Data do Fascículo
Jun 2005
Histórico
-
Recebido
14 Set 2004 -
Aceito
16 Dez 2004