Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

MISE EN PLACE D'UN CIRCUIT CHUA POUR DÉMONTRER

BIFURCATIONS ET ÉTRANGES ATTRACTEURS DANS UNE CLASSE

Ahmad Byagowi1 et Witold Kinsner1,2,3

1
Département de génie électrique et informatique
Université du Manitoba, Winnipeg, MB, Canada R3T 5V6
et 2L'Institut des Sciences Mathématiques Industrielles
1
et Laboratoires de recherche en télécommunications, TRLabs
{[email protected] | [email protected]}

Résumé : Cet article décrit la conception et la mise en œuvre 1. INTRODUCTION

d'un circuit Chua à double défilement pour démontrer le chaos
dans les systèmes dynamiques aux étudiants d'un cours 1.1 Utilité des démonstrations en classe
d'études supérieures afin d'améliorer leur visualisation et leur
compréhension des attracteurs étranges et des arbres de Enseigner aux étudiants en génie du premier cycle et des cycles
bifurcation de Feigenbaum. supérieurs dans un environnement de classe implique généralement
L'enseignement des systèmes dynamiques (c'est­à­dire un certain nombre de techniques, notamment des présentations orales et
des systèmes non linéaires pouvant présenter du chaos) est des discussions, du matériel visuel (images et vidéo) pour
souvent considéré comme difficile en raison de la modélisation augmenter les idées difficiles, des devoirs écrits, des quiz, des
mathématique impliquée et de l'inclusion du quatrième état tests et des examens, ainsi que des laboratoires pratiques. Les
d'attracteur étrange, en plus de la stabilité ponctuelle démonstrations en classe d'objets existants (puces, processeurs,
traditionnelle, de la stabilité cyclique et de la stabilité toroïdale, cartes, dispositifs électromécaniques, sous­systèmes et
comme trouvé dans les systèmes dynamiques. Un cours systèmes) peuvent également être très efficaces. Les
d'études supérieures est offert à l'Université du Manitoba depuis démonstrations développées spécifiquement pour le cours par
de nombreuses années pour fournir à la fois (i) une théorie les participants au cours sont particulièrement importantes, car
unifiée des dimensions fractales, ainsi que de nombreuses ils comprennent les défis liés à l'acquisition de matériel difficile.
mises en œuvre pratiques d'algorithmes pour calculer les Dans le cas de systèmes dynamiques (c'est­à­dire des
dimensions fractales, y compris le spectre des dimensions systèmes dynamiques non linéaires pouvant présenter du chaos),
Rényi requis pour la caractérisation. des attracteurs étranges en utilisant l'analyse
le matériel multifractale.
enseigné nécessite une sophistication considérable de
Leon Chua a développé un circuit non linéaire simple la part des étudiants, et les explications par les techniques
capable de produire une riche collection de phénomènes traditionnelles peuvent ne pas être adéquates [Kins12]. Par
dynamiques, allant des points fixes aux points de cycle, aux conséquent, l'utilisation de démonstrations d'exemples réels de
bifurcations standards (doublement de période), à d'autres systèmes dynamiques peut faciliter la compréhension des
routes standards vers le chaos et au chaos lui­même. La raison concepts sous­jacents, en particulier lorsque les démonstrations
pour laquelle ce circuit spécifique a été choisi comme outil de peuvent se comporter de manière entièrement contrôlée.
démonstration de classe est triple : (i) le circuit possède un Il existe de nombreux exemples de systèmes dynamiques
modèle analytique et peut être simulé, (ii) le circuit peut être qui peuvent être utilisés à des fins de démonstration. Cependant,
mis en œuvre à l'aide de composants disponibles dans le seuls quelques systèmes non linéaires satisfont simultanément
commerce, et (iii) les signaux du circuit peuvent être acquis aux trois exigences clés suivantes : (i) ils peuvent être décrits
sans affecter ni altérer de manière significative son fonctionnement.par un modèle mathématique approprié (par exemple, des
Cet article décrit le l'architecture, équations différentielles), peuvent être modélisés sur un
la mise en œuvre, la vérification et les tests du système Chua, ordinateur, et (iii) peuvent être applicable en laboratoire pour
ainsi qu'une analyse des données obtenues au cours de la des démonstrations en classe et des recherches. De tels
phase actuelle de développement. Bien qu'il existe de exemples doivent également satisfaire à des exigences
nombreuses implémentations possibles du circuit de Chua, supplémentaires, notamment : (i) être peu coûteux à mettre en
notre implémentation présente plusieurs caractéristiques de œuvre avec des composants commerciaux disponibles sur
conception innovantes pour la rendre plus applicable et étagère (COTS), (ii) faciles à entretenir, (iii) faciles à modifier,
améliorer l'apprentissage des élèves en classe. (iv) durables dans de multiples utilisations. , et doit être facile à
surveiller sans affecter la qualité du signal mesuré (les systèmes
Mots­clés : Chaos ; visualisation du chaos ; diagramme de dynamiques sont sensibles aux conditions initiales et donc aux sondes des capteurs
bifurcation ; Circuit Chua ; Potentiomètre numérique multi­prises.

CEEA Conf. 2012 ; Document ChuaChaos­Full­v40.docx

096 Winnipeg, MB; 17­20 juin 2012 – 1 sur 10 – 6 avril 2012
Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

2. CONTEXTE DU CIRCUIT CHUA Troisièmement, étant donné que le comportement chaotique se produit dans
les boucles de tension et de courant électriques du circuit, le circuit peut être

2.1 Le circuit Chua comme système dynamique détecté (surveillé) à l'aide d'instruments de laboratoire tels qu'un oscilloscope, ou
il peut être capturé à l'aide d'instruments d'acquisition de données standard à

Vers 1983, Leon Chua a développé un circuit non linéaire haute impédance. . Étant donné que la résistance est le seul composant qui a une

(aujourd'hui connu sous le nom de circuit Chua) capable de relation indépendante du temps entre la tension et le courant (tracé VI), la

démontrer une riche collection de dynamiques. modification de la valeur de la résistance dans le circuit Chua peut être considérée

phénomènes, allant des points fixes aux points de cycle, les comme

bifurcations standards (doublement de période), d'autres routes

standards vers le chaos et le chaos lui­même (par exemple, la meilleure façon de modifier le comportement du circuit.
[Chua07], [Chua92], [ChPW95], [Chua06], [PaCh89], [ Kapi02],
[Ogor97], [Deva92], [Ott93], [AlSY96], [PeJS92], [Strog94],
[BrHT10], [MaSi05], [OtSY94], [Kins12]). Quelques années plus
tard, Chua, Komuro et Matsumoto ont prouvé rigoureusement que
l'attracteur Chua est effectivement chaotique.
[ChKM86]. De plus, le circuit peut répondre à toutes les exigences
décrites dans la section précédente.
Comme le montre la figure 1, le circuit Chua se compose de
trois composants de stockage d'énergie (deux condensateurs C1
et C2, une inductance L), une résistance R (ou de manière
équivalente conductance G = 1/R) et une résistance non linéaire
rN ou conductance. gN tel que la diode Chua représentée sur la figure 2.
Cet élément non linéaire fournit l'énergie nécessaire au circuit. Notez
Fig. 2 Résistance négative dans le circuit Chua.
que nous utilisons des lettres majuscules pour désigner des valeurs
(De [Kins12])
constantes et des lettres minuscules pour désigner des variables et
des valeurs qui ne sont pas constantes dans le temps. Le circuit
2.2 Modèle de circuits
comporte deux boucles : la boucle non linéaire avec C1 (à gauche) et
la boucle purement réactive avec C2 (à droite).
Le circuit peut être modélisé comme un circuit RLC du troisième ordre avec
Les deux boucles sont couplées par la résistance R. Lorsque R
quatre éléments linéaires et un élément non linéaire. À partir de la loi de Kirchhoff,
est très grand, le couplage est faible, et il est fort lorsque R
on peut écrire l'ensemble suivant de trois équations du premier ordre [Kapi02, Sec.
est petite. La résistance non linéaire rN et le couplage entre les
deux boucles sont en grande partie responsables du comportement
6.4], [Kins12; Ch1, Partie3]
extraordinaire de ce circuit simple.

" dvC1 1 1
$
= (vC2 ! dans) C1 !

g(vC1 )
dt RC1 C1
$

dvC2 1 1
= (vC1 ! dans) C2 +
$$

# je L (1)
dt RC2 C2
$

langue
$
= !vC2
dt
%
$

Fig. 1 Circuit Chua avec deux condensateurs (C1 et C2), une où la conductance non linéaire est approximée par un
inductance (L), une résistance (R) et une résistance non linéaire fonction linéaire antisymétrique par morceaux en trois segments
(RN). (De [Kins12])
(Fig. 2), de pentes Ga (dans la partie médiane) et Gb

Le circuit satisfait aux trois exigences. Premièrement, le circuit (les sections extérieures) et les points d'arrêt à Bp = ±1
peut être modélisé à l’aide de techniques de circuit standard.
1
() # %
g(x* ! vg ) = Gb x *+ Ga " Gb x *+Bp
$ " x *"Bp &
Deuxièmement, le circuit est implémentable, maintenable et 2

peut être construit à l’aide de composants COTS. Le comportement (2)

du circuit dépend exclusivement des valeurs de ses composants.
La modification du circuit peut être réalisée soit en utilisant Ce système dynamique autonome (c'est­à­dire un système qui
différents composants à valeurs discrètes, soit en utilisant des évolue grâce à sa dynamique naturelle, sans aucune
composants variables.

CEEA Conf. 2012 ; Document 096 ChuaChaos­Full­v40.docx

Winnipeg, MB; 17­20 juin 2012 – 2 sur 10 – 6 avril 2012
Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

injection de signal externe dans le circuit) peut être davantage 3. CONCEPTION DU MATÉRIEL DU CIRCUIT CHUA
simplifié sous la forme sans dimension suivante
3.1 Conceptions alternatives
$ dx
=!
( ) y " x " h(x) Plusieurs implémentations alternatives pour le circuit Chua ont
dt
&

&
été envisagées pour la démonstration du cours fractal et chaos
toi
(3) [Kins12]. Le circuit sélectionné offre une flexibilité permettant de
&

% =y"x+z
&
dt démontrer différents états du circuit Chua et peut fournir un contrôle
détaillé du comportement du circuit.
langue
#et
&

="
&

'
dt
En utilisant un seul potentiomètre numérique multiprise au lieu de
où la résistance constante utilisée dans le circuit Chua, un mode de
réalisation flexible du circuit a été mis en œuvre. Le potentiomètre
$
vC1 vC2 je
L numérique peut être réglé soit manuellement par l'utilisateur, soit
&
X= y= z=
& Pb Pb BpG automatiquement par un microcontrôleur. Le changement de valeur
&

tG du potentiomètre étant automatisé, différents schémas de contrôle

&
! = peuvent être mis en œuvre afin d'ajuster le comportement du circuit
&

C2 Chua. De plus, en modifiant la valeur de la résistance à plusieurs

% (4)
Ici Go reprises, la transition des états (doublement de période) dans le circuit
&

une = b=
&

g g Chua peut être observée clairement et l'ensemble du diagramme de

&

bifurcation peut être démontré.

&

" = C2
#=
C2
LG2
&

C1
'
et 3.2 La diode Chua
1
h(x) = bx + un ! b " x +1 ! x !1 $%
()# (5) L'élément non linéaire rN n'existant pas en tant que composant
2
passif, il doit être implémenté à l'aide d'éléments actifs tels que
des amplificateurs opérationnels (amplis opérationnels). Comme
Notez que x est lié à la tension sur C1, y se rapporte à la tension sur
C2 et z représente le courant circulant dans l'inducteur. Les paramètres
il existe de nombreuses implémentations du circuit Chua, il existe
critiques sont les rapports de conductance de la résistance non linéaire
de nombreuses implémentations alternatives de l'élément non
dans linéaire (également connu sous le nom de diode Chua).
La figure 3 montre le schéma d'une implémentation matérielle
sa conductance médiane Ga à la conductance de couplage G, et
populaire du circuit Chua avec une implémentation à double ampli
la conductance externe Gb à G. Les autres paramètres α et β opérationnel de la diode Chua non linéaire.
relient respectivement la valeur du condensateur C2 à C1 et à L.

Le circuit affiche des classes de comportement distinctes pour

différentes combinaisons de α, β, a et b [Kins12]. Par exemple, le
chaos complet à double défilement se développe pour α = 10,0, β
= 14,87, a = –1,27 et b = –0,68.

2.3 Autres circuits chaotiques

En plus du circuit Chua original, il existe de nombreuses

variantes. Par exemple, des treillis de ce type
les circuits produisent un comportement plus complexe.
Fig. 3. Implémentation matérielle utilisée pour le circuit Chua.
Il existe également de nombreux circuits chaotiques simples similaires.
Par exemple, Sprott a introduit des circuits chaotiques « jerk » (De [Kenn92])
décrit par une seule équation différentielle ordinaire autonome du
troisième ordre, avec diverses non­linéarités, souvent représentées La diode Chua a été implémentée par de nombreux auteurs
par des approximations linéaires par morceaux pour faciliter une (par exemple, [ElKe00], [Cros03]), mais un circuit robuste a été
implémentation électronique robuste avec des diodes et des développé pour la première fois par Kennedy [Kenn92], [KeRS00].
amplificateurs opérationnels [Spro00a], [Spro00b]. Ces circuits sont Une bonne discussion sur les amplis opérationnels liés à de tels
également évolutifs sur une large gamme de fréquences. circuits peut être trouvée, par exemple, dans [HoHi89, Ch. 4],
[SeSm09], [NaSe03].

CEEA Conf. 2012 ; Document 096 ChuaChaos­Full­v40.docx

Winnipeg, MB; 17­20 juin 2012 – 3 sur 10 – 6 avril 2012
Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

potentiomètre numérique pour faire varier le comportement du circuit.

3.3 Alimentation isolée Un potentiomètre numérique multi­tap (1024­tap) (MAX5483
de MAXIM) est sélectionné [Maxi10]. Un schéma fonctionnel de ce
La mise en œuvre du circuit Chua illustré à la figure 3 est basée composant est présenté à la figure 5.
sur une puce à double ampli opérationnel (TL082CN) qui
nécessite une double alimentation de ±15 V. De plus, le circuit Chua 3.5 Microcontrôleur
produit une large gamme de fréquences qui peuvent interférer avec
les composants numériques utilisés dans le circuit. Par exemple, le Le circuit Chua de base ne nécessite aucun ordinateur, ni
potentiomètre numérique est basé sur des circuits électroniques pour microprocesseur, ni microcontrôleur. Notre circuit Chua automatisé
imiter la variabilité de l'impédance et nécessite des conditions de décrit dans cet article nécessite un microcontrôleur tel que le
tension strictes de la part d'une alimentation. Ainsi, la robustesse du PIC16F876A [Micr03]. Ce contrôleur a été sélectionné en raison de la
circuit pourrait être augmentée en isolant l'alimentation du potentiomètre grande expérience et du soutien disponible à l'Université du Manitoba.
numérique et de son circuit de commande. Cela a été réalisé par un De nombreux autres microcontrôleurs intégrés pourraient également
convertisseur DC­DC qui convertit une seule alimentation de 5 V en être utilisés pour exécuter les fonctions requises.
±15 V pour fournir la tension aux amplificateurs opérationnels. Le
DCP0205­15DP a été sélectionné à cet effet et son schéma fonctionnel L'une des fonctions est de contrôler le potentiomètre numérique.
est présenté sur la figure 4 [Burr99]. Le pot reçoit des commandes via le
bus d'interface périphérique série (SPI), selon un protocole spécifique.
Le microcontrôleur envoie les valeurs de la résistance requise du
potentiomètre numérique entre deux points de consigne pendant une
durée souhaitée, produisant ainsi une forme d'onde en dents de scie,
comme le montre la figure 6.

Figure 4. Convertisseur CC­CC DDCP0205­15DP.

3.4 Potentiomètre numérique multiprise

Fig. 6. Signal en dents de scie utilisé pour faire varier la valeur

du potentiomètre.

La forme d'onde en dents de scie est implémentée en utilisant

code logiciel sur le microcontrôleur PIC16F876A. Le microcontrôleur
augmente la valeur de son registre. Lorsque la valeur du registre
atteint le point de consigne haut, il place la valeur du point de consigne
bas sur le registre. L'incrémentation se répète afin de fournir le
périodique
forme d'onde dans une plage requise par le mode de fonctionnement
spécifique du circuit Chua. Les deux points de consigne sont trouvés
afin de régler le circuit Chua variant entre p0 et p∞. L'implémentation
du générateur de forme d'onde en dents de scie se fait d'abord en
langage C, et
5. Potentiomètre numérique MAX5483.
puis le code est compilé et programmé dans le microcontrôleur, à
(Depuis [Maxi10])
l'aide du progiciel fourni pour les microcontrôleurs PIC.

Le comportement du circuit Chua dépend de toutes les valeurs de

ses composants réactifs, passifs et actifs. Afin de démontrer le
3.6 Implémentation de la maquette
diagramme de bifurcation sur un oscilloscope, le comportement du
circuit Chua doit se répéter progressivement et à plusieurs reprises
Afin de faciliter la mise en œuvre, une maquette a été sélectionnée
du premier point de bifurcation p0 au dernier point de bifurcation p∞ . pour la réalisation du circuit Chua et de la section contrôleur numérique.
Dans ce travail, la résistance de couplage R est remplacée par une
résistance multiprise

CEEA Conf. 2012 ; Document 096 ChuaChaos­Full­v40.docx

Winnipeg, MB; 17­20 juin 2012 – 4 sur 10 – 6 avril 2012
Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

unité. Cette démonstration de bifurcation constitue le premier objectif

principal de cette expérience.
La figure 8 montre une capture d'écran d'un diagramme de bifurcation
répété sur un oscilloscope numérique (Tektronix
TDS1012B en utilisant sa fonction de capture d'écran intégrée).
Afin d'obtenir une résolution plus élevée du diagramme de bifurcation,
la période de rampe de résistance du
le potentiomètre numérique est réglé sur 10 secondes. En d’autres termes,
Fig. 7. Implémentation matérielle de la conception. la résistance du circuit Chua est effacée du point de consigne bas vers le
point de consigne haut en 10 secondes.
La figure 7 montre le circuit Chua sur le côté gauche de la
La figure 9 montre un diagramme de bifurcation de meilleure résolution
photographie et le microcontrôleur sur le côté droit. Le potentiomètre
afin d'obtenir plus de détails. Ceci a été réalisé en diminuant le paramètre
numérique se trouve au milieu de la maquette. Plusieurs indicateurs
de l'axe temporel de l'oscilloscope de 2,5 secondes à 1 seconde par unité
LED ont également été installés pour indiquer le fonctionnement du
de grille d'écran.
circuit. Les valeurs des composants utilisés pour le circuit Chua, à
l'exclusion de la mise en œuvre de la diode Chua, sont indiquées
dans le tableau 1.

Tableau 1. Valeur des composants utilisés dans le Chua

mise en œuvre des circuits.
Valeur du composant
Inductance L 18 mH

Condensateur C1 10 nF
Condensateur C2 100 nF
Résistance R 1 kΩ à 2 kΩ

4. RÉSULTATS EXPÉRIMENTAUX
4.1 Diagramme de bifurcation

Fig. 9. Un diagramme de bifurcation à plus haute résolution.

On voit sur la figure 9 que le côté gauche de la capture d'écran

représente une solution unique (1 période) qui augmente de manière
monotone jusqu'à la première bifurcation au milieu de l'écran. La bande
s'étend jusqu'aux trois quarts de l'écran où elle s'étend à la totalité de
l'écran. Le régime chaotique est indiqué par les points répartis sur la bande.

4.2 Affichage de l'attracteur Chua Strange en 2D

Le deuxième objectif principal de cette expérience est de démontrer

que le circuit Chua possède un attracteur étrange (ce nom a été inventé
par Ruelle et Takens en 1971 [RuTa71]). Si l'attracteur existe, le processus
Fig. 8. Diagramme de bifurcation capturé sur un oscilloscope
numérique (Tektronix TDS1012B). ne peut pas être aléatoire [Kins12]. Puisque le circuit Chua a trois degrés
de liberté, son espace d’états est tridimensionnel (3D). Une façon d'afficher
Dans tout cours sur les systèmes dynamiques, le diagramme de un attracteur étrange consiste à prendre une paire de variables d'état (par
bifurcation est l'un des outils de visualisation les plus importants. exemple, x et y) et à les transmettre aux deux déflecteurs (X, Y) d'un
S’il existe un diagramme de bifurcation, le processus sous­jacent peut oscilloscope.
conduire au chaos. Ce circuit a été conçu pour démontrer le diagramme
de bifurcation sur un oscilloscope à stockage ou via toute acquisition de Par exemple, la figure 10 montre une projection 2D du modèle 3D.
données appropriée. espace d'état du circuit Chua, alors qu'il est en mode attracteur étrange à
double défilement.

CEEA Conf. 2012 ; Document 096 ChuaChaos­Full­v40.docx

Winnipeg, MB; 17­20 juin 2012 – 5 sur 10 – 6 avril 2012
Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

Ce système d'acquisition de données enregistre 10 secondes de

données sur les trois canaux (x, y, z) du circuit Chua et stocke les
données sur un PC à l'aide de MatLab Simulink (MDL) [NaIn08]. Afin
de capturer la plus haute qualité des signaux à large bande générés
par le circuit Chua, la bande passante disponible la plus élevée pour
l'acquisition de données a été utilisée (50 000 échantillons par
seconde, avec une résolution de 8 bits). La longueur de l'enregistrement
est de 5 165 229 octets (4,92 Mio). À ce rythme, une carte SD de 4
Go peut stocker plus de 2,5 heures (8 130 s) de données.

La figure 12 montre une capture d'écran d'une fenêtre des

données capturées représentant les trois variables d'état du circuit
Chua. La capture d'écran a été obtenue à l'aide d'une fonction intégrée
de l'oscilloscope couleur à stockage numérique Tektronix TDS2014B.
Les trois canaux représentés sur la figure 12 sont la tension électrique
sur le condensateur C1 (Ch. 1), la tension sur le condensateur C2
(Ch. 2) et le courant électrique à travers l'inductance L (Ch. 3).

Fig. 10. Attracteur étrange à double défilement obtenu à partir du circuit

Chua (HP54603B).

On voit que l’étrange attracteur est en effet un double rouleau.

Bien que n'importe quelle paire de variables d'état puisse être utilisée
pour afficher l'attracteur étrange, les variables x et y
représentent les tensions sur les condensateurs qui peuvent être
facilement exploitées.

4.3 Affichage de l'attracteur Chua Strange en 3D

Puisqu'un oscilloscope ne peut pas afficher les trois dimensions

de l'attracteur réel, l'espace d'état 3D du circuit peut être reconstruit à
partir des trois dimensions x, y, z.
signaux, tels qu'acquis à partir du circuit Chua, à l'aide d'une unité
d'acquisition de données et d'un ordinateur. L'acquisition de données
sélectionnée pour cette expérience est la carte NI­DAQCard6062E de
National Instruments [NaIn05], comme le montre la Fig. 11. La carte Fig. 12. Données triple canal obtenues à partir du circuit Chua.
peut échantillonner à une vitesse allant jusqu'à 500 kilosamples par
seconde (kS/s) et 12 bits. résolution. Une boîte à bornes SCB068 est
également utilisée pour réduire le bruit et Les données ci­dessus (x, y, z) acquises à partir du circuit Chua
distorsion du signal. ont ensuite été affichées sous forme de tracé 3D, à l'aide de Mathworks
Maltab [Math12], et le tracé est illustré à la figure 13.

Fig. 11. Acquisition de données utilisée pour reconstruire l'espace

Fig. 13. Attracteur tridimensionnel à double étagement du circuit Chua.
d'état tridimensionnel.

CEEA Conf. 2012 ; Document 096 ChuaChaos­Full­v40.docx

Winnipeg, MB; 17­20 juin 2012 – 6 sur 10 – 6 avril 2012
Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

Notez que bien que le tracé ait été tourné par rapport à
l'affichage de l'oscilloscope de la figure 10, le caractère à double
défilement de l'objet reste inchangé. L'ordre des canaux a peut­être
été inversé lors du tracé des diagrammes.

4.3 Reconstruction d'un attracteur étrange

Le troisième objectif principal de cette expérience est de

reconstruire l’étrange attracteur à partir d’une seule série temporelle
acquise auprès du circuit Chua.

Reconstruction 2D : une projection 2D de l'étrange attracteur peut Fig. 15. Attracteur étrange 3D reconstruit à partir de vC1.
être reconstruite à partir d'un seul signal
observation, en traçant l'échantillon actuel du signal u[n+1] par Y a­t­il quelque chose de spécial dans la tension sur C1 qui
rapport à l'échantillon précédent du signal u[n]. nous permet de reconstruire l'étrange attracteur 3D avec une si
La figure 14 montre l'attracteur étrange 2D reconstruit à partir de la haute fidélité ? La réponse est clairement négative. Étant donné
tension du signal du condensateur C1 (vC1). que chaque série temporelle d'observation contient les mêmes
informations sur l'ensemble du système dynamique, l'attracteur
étrange pourrait être reconstruit à partir des deux autres observations,
0,8
soit la tension de C2, vC2, soit le courant traversant L, iL.
0,6

0,4
4.4 Reconstruction d'attracteurs inconnus
0,2

Dans le cas du circuit Chua, on sait que son espace des phases
CV1(n+1)

0
a trois dimensions. Ainsi, l’étrange attracteur Chua reconstruit doit
« vivre » en trois dimensions.
­0,2

On dit que la dimension d'encastrement de l'attracteur est DE = 3.

­0,4

Takens (par exemple, [Take81], [Take10]) a fourni un théorème

­0,6
remarquable qui stipule : un attracteur étrange, s'il existe, peut être
­0,8
reconstruit à partir d'une seule série temporelle, même si nous ne
­0,8 ­0,6 ­0,4 ­0,2 0 0,2 0,4 0,6 0,8
connaissons pas le nombre de dimensions d'intégration DE, = m et
CV1(n)
le décalage temporel τ entre les échantillons à impliquer dans la
Fig. 14. Attracteur étrange reconstruit en 2D à partir de VC1.
reconstruction. Pour clarifier cette affirmation, considérons une série
temporelle originale avec N
Ce résultat reconstruit à partir d'une seule série temporelle vC1
échantillons quantifiés désignés par
est remarquable car il ressemble à l'étrange attracteur affiché sur
l'oscilloscope (Fig. 10), généré par deux signaux distincts, vC1 et vC2.
{ } u[n] n = 1,2,3,...,N (6)
Puisque l’attracteur étrange a une forme non aléatoire et sans
ambiguïté, le processus sous­jacent ne peut pas être aléatoire. La
seule limitation du graphique 2D ci­dessus est que le double défilement L'étrange attracteur pourrait être reconstruit en plaçant ses points
se croise. successifs sur la trajectoire dans un pseudo­espace de phase à m
dimensions. Chaque point U[n] sur la trajectoire a m coordonnées
sélectionnées dans la série temporelle d'origine selon
Reconstruction 3D : la reconstruction de l'attracteur ci­dessus
Le processus peut être étendu à l'ensemble de l'attracteur étrange
3D en traçant les échantillons de signal avec un et deux retards, U[n] = (u[n],u[n + ! ]) ],u[n + !2],...,u[n + (m "1) !
u[n+2] contre u[n+1] contre u[n]. La figure 15 montre un tracé 3D de
l'attracteur étrange, également reconstruit à partir de vC1. Cet pour n = 1,2, 3,...,N ! (m !1) (7)
attracteur étrange 3D reconstruit est très similaire à l'attracteur
étrange assemblé à partir des trois données réelles capturées par Il existe plusieurs bonnes techniques pour trouver la dimension
l'acquisition de données.
d'intégration entière DE et le décalage temporel τ (par exemple,
[Addi97], [Will97], [SaYC91]). Puisque les deux paramètres

CEEA Conf. 2012 ; Document 096 ChuaChaos­Full­v40.docx

Winnipeg, MB; 17­20 juin 2012 – 7 sur 10 – 6 avril 2012
Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

doit être estimé à partir de la série chronologique elle­même, l’étrange papier. Le circuit Chua est facile à mettre en œuvre car il comporte
attracteur « vit » désormais dans un pseudo espace de phase. peu de composants et peut être construit à l’aide de composants
L'une des meilleures techniques pour trouver τ est le premier disponibles dans le commerce. Notre implémentation de circuit
critère d'information mutuelle minimale [Will97; Ch. 19]. Chua peut démontrer l’intégralité du diagramme de bifurcation et
Une bonne technique pour trouver DE est le faux critère de peut être utilisée pour afficher l’étrange attracteur.
distance du voisinage le plus proche [Addi97]. Le principe de cette Le potentiomètre multi­tap piloté par microcontrôleur permet d'effectuer
technique est que si la dimension de plongement du pseudo­ la démonstration rapidement et de manière répétée. Le temps gagné
espace des phases DP est trop petite, ou si DP < DE, les distances peut ensuite être consacré à la discussion des résultats.
entre les « faux » voisins changent. Cependant, lorsque DP ≥ DE,
les distances entre les « vrais » voisins cessent de changer. Cette implémentation peut également offrir la possibilité
d’automatiser les passages vers différents états du circuit Chua.

5. VÉRIFICATION DES RÉSULTATS

Remerciements
Pour vérifier que la fréquence d'échantillonnage de 50 kS/
s était suffisante pour ce circuit Chua, la stationnarité des trois Nous tenons à remercier le Département de génie électrique
signaux (vC1, vC2, iL) a été vérifiée (le premier et le deuxième et informatique de l'Université du Manitoba pour son soutien
instants uniquement), puis une analyse spectrale a été réalisée. financier partiel à ce projet.
pour établir la portée du signal et la fréquence de Nyquist a été
calculée. Les références
La faible stationnarité de chaque canal a été vérifiée sur toute
la longueur de l'enregistrement en calculant les deux premiers [Addi97] Paul S. Addison, Fractale et chaos : un
moments (moyenne et variance) pour trois fenêtres de tailles 100,
Cours illustré. Bristol, Royaume­Uni : Institute of Physics
1 000 et 10 000 échantillons, avec un chevauchement de 50 %.
Publishing, 1997, 256 pp. {ISBN 0­7503­0400­6 pbk}
Les plages de signaux ont également été calculées à partir
des données collectées. Les plages sont reproduites ici pour [AlSY96] Kathleen T. Alligood, Tim D. Sauer et
donner une idée des grandeurs attendues de tels circuits. James A. Yorke, Chaos. New York, NY : Springer­
Publishers, 1996, 603 pages. {ISBN 0­387­94677­2 pbk}
Tableau 2 : Plages de signaux. [BaBl09] Gregory L. Baker et James A. Blackburn, Le
pendule : une étude de cas en physique. Oxford,
Signal Valeur Min Valeur Max Crête­Crête
Royaume­Uni : Oxford University Press, 2009, 350 pages.
vC1 –0,7070 V 0,7827 V 1,4896 V
{ISBN : 978­0199557684}
vC2 –2,9316 V 3,1636 V 6,0952 V
Le ­41 mA 52 mA 92 mA [BlBa98] James A. Blackburn et Gregory L. Baker, "Une
comparaison des pendules chaotiques commerciaux", Am.
J. Phys., vol. 66, non. 9, pp. 821­830, septembre 1998.
Une analyse de Fourier a également été réalisée sur les trois Disponible à partir d'avril 2012 à partir de
canaux. Les spectres sont bien distincts, avec des coupures http://fraden.brandeis.edu/courses/phys39/chaos/Dae
vC1 = 5 700 Hz, vC2 = 3 200 Hz et iL = 3 000 Hz. Cela signifie dalon%20AJP.pdf
que la fréquence d'échantillonnage réelle utilisée dans l'expérience
[BrHT10] HW Broer, B. Hasselblatt et F. Takens (éd.),
était supérieure à la fréquence de Nyquist requise.
Handbook of Dynamical Systems, Vol. 3, Amsterdam :
Hollande du Nord/Elsevier, 2010. {978­0­
6. REMARQUES CONCLUSIONS
444­53141­4}
L'enseignement des systèmes dynamiques peut être amélioré [Burr99] Burr­Brown, miniature 2W isolé
en simulant et en démontrant le comportement de divers systèmes Convertisseurs DC/DC non régulés, DDCP0205­15DP
physiques. Très souvent dans le passé, les systèmes étaient Fiche de données. Tuscon, Arizona : Burr­Brown
mécaniques, comme les systèmes à pendule entraîné Corporation, 1999. Disponible le 30 mars 2012 auprès de
de Daedalon Corporation [SmBl06], Leybold Didactic, Pasco
www.datasheetcatalog.org/datasheets/120/193437_D
Scientific et TELAtomic, Inc. [BlBa98].
S.pdf
Même si les systèmes mécaniques étaient souvent coûteux et
difficiles à entretenir en classe, ils offraient une excellente [Chua92] Leon O. Chua, « La genèse du circuit de Chua »,
expérience pratique aux étudiants et aux chercheurs [BaBl09]. De Archives pour l'électronique et la technologie de transmission,
nombreux systèmes alternatifs sont électroniques, comme le vol. 46, n° 4, pp. 250­257, 1992. Disponible le 30 mars 2012
circuit Chua décrit dans ce auprès de

CEEA Conf. 2012 ; Document 096 ChuaChaos­Full­v40.docx

Winnipeg, MB; 17­20 juin 2012 – 8 sur 10 – 6 avril 2012
Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

www­ http://nonlinear.eecs.berkeley.edu/chaos/RobustOpA
inst.eecs.berkeley.edu/~ee129/sp09/handouts/Genesi mpRealizationOfChuaCircuit.pdf
sChuasCircuit.pdf [email protected]
[Chua06] Leon O. Chua, « Le circuit de Chua 10 ans plus tard », [Kins12] W. Kinsner, Ingénierie fractale et chaotique.
Journal international de théorie et d'applications des Notes de lecture. Winnipeg, Manitoba : Université du Manitoba,
circuits, vol. 22, non. 4, pp. 279­305, juillet­août 1994. {DOI : 2012, 941 p.
10.1002/cta.4490220404}. [MaSi05] Ganesh B. Maganti et Sahjendra N. Singh,
[Chua07] Leon O. Chua, « Circuit Chua », Scholarpedia, « Contrôle adaptatif modulaire du chaos dans le circuit de Chua »
tome 2, non. 10, pp, 1488, 2007. Journal international de bifurcation et de chaos en sciences
{doi:10.4249/scholarpedia.1488} Disponible le 30 mars 2012 appliquées et en ingénierie, vol. 15, non. 9, p.
sur 2973­2984, septembre 2005.
http://www.scholarpedia.org/article/Chua_circuit [Math12] MathWorks, MATLAB et Simulink. Natick, MA :
[ChKM86] Leon O. Chua, Motomasa Komuro et Takashi MathWorks, 2012. Disponible depuis le 30 mars 2012 auprès
Matsumoto, « La famille à double parchemin », de

Transactions IEEE sur les circuits et les systèmes, vol. 33, non. www.mathworks.com/

11, pp. 1072­1118, novembre 1986. [Maxi10] Maxime, 10 bits, MAX5481, MAX5482,
[ChPW95] Leon O. Chua, Ladislav Pivka et Chai Potentiomètres numériques non volatils à cône linéaire
Wah Wu, « Un circuit universel pour étudier les phénomènes MAX5483, MAX5484, fiche technique. Sunnyvale, Californie :
chaotiques », Phil. Trans. R. Soc. Londres. A, vol. 353, p. 65­84, Maxim Integrated Products, 2010. Disponible à compter du
16 octobre 1995. 30 mars 2012 auprès de
[Cros03] Michael Cross, « Le circuit de Chua », 2003. http://www.maxim­
Disponible à compter du 30 mars 2012 auprès de ic.com/datasheet/index.mvp/id/4730
http://www.cmp.caltech.edu/~mcc/chaos_new/Chua. [Micr03] Micropuce, microcontrôleur PIC16F876A,
HTML Fiche de données. Chandler, Arizona : Microchip Technology,
[Deva92] Robert L. Devaney, Un premier cours sur les systèmes 2003, 232 pages.
dynamiques chaotiques : théorie et expérience. Boston, MA : [NaIn05] National Semiconductor, DAQCard­6062E 12­
Addison­Wesley, 1992, 302 pp. {ISBN 0­201­ Spécifications de la famille PCMCIA (PC Card) d'E/S
55406­2} multifonctions bit. Austin, Texas : National Semiconductor,
[ElKe00] AS Elwakil et Michael P. Kennedy, « Implémentation décembre 2005. Disponible le 30 mars 2012 auprès de
améliorée de l'oscillateur chaotique de Chua à l'aide de www.ni.com/pdf/manuals/370724c.pdf
l'ampli­op à rétroaction de courant », IEEE Trans. Système de [NaIn08] National Instruments, Multisim (édition professionnelle,
circuits. I : Théorie fondamentale et applications, vol. 47, version d'essai de 30 jours). Austin, Texas : février 2008.
non. 1, pages 76 à 79, 2000. Disponible à partir du 30 mars 2012 auprès de
[HoHi89] Paul Horowitz et Winfield Hill, L'art de l'électronique. http://www.ni.com/multisim/
Cambridge, Royaume­Uni : Cambridge University Press, 1989 [NaSe03] National Semiconductor, amplificateur opérationnel
(2e éd.), 1 152 pages. {ISBN 978­ d'entrée et de sortie double rail­à­rail CMOS
0521370950} LMC6482AIN, fiche technique DS011713.
[Kapi02] Tomasz Kapitaniak, Chaos pour les ingénieurs. National Semiconductor, septembre 2003. Disponible le 30
New York, NY : Springer, 2002 (2e éd.), 142 p. mars 2012 auprès de
{ISBN 3­540­66574­9 pbk} http://nonlinear.eecs.berkeley.edu/chaos/LMC6482.p
[KeRS00] Michael Peter Kennedy, Riccardo Rovatti et Gianluca df

Setti, Chaotic Electronics in Telecommunications. [SmBl06] HJT Smith et James A. Blackburn, Guide de l'utilisateur du
Boca Raton, Floride : CRC Press, 2000, 445 pages {ISNBN pendule chaotique EM­52. Salem, MA : Daedalon Corporation,
0­8493­2348­7} 1998. 34 p. Disponible en avril 2012 auprès de
[Kenn92] Michael Peter Kennedy, « Ampli opérationnel robuste
réalisation du circuit de Chua », Frequenz, vol. 46, non. 3­4, pp. http://fraden.brandeis.edu/courses/phys39/chaos/
66­80, mars­avril 1992. Disponible le 30 mars 2012 auprès de EM­52%20Daedalon%20Manual.pdf

CEEA Conf. 2012 ; Document 096 ChuaChaos­Full­v40.docx

Winnipeg, MB; 17­20 juin 2012 – 9 sur 10 – 6 avril 2012
Machine Translated by Google

Proc. 2012 Association canadienne de formation en ingénierie (CEEA12) Conf.

[Spro00a] JC Sprott, « Systèmes chaotiques simples et Physique, vol. 20, pages 167 à 192, 1971 ; et vol. 23, pages
circuits », Am. J. Physique, vol. 68, non. 8, pp. 758­763, août 343­344, 1971.
2000. [SaYC91] Tim Sauer, James A. Yorke et Martin Casdagli,
[Spro00b] JC Sprott, « Une nouvelle classe de circuits chaotiques », « Embedology », Journal of Statistical Physics, vol. 65,
Lettres de physique A, vol. 266, p. 19­23, 14 février 2000. non. 3­4, pp. 579­616, 1991. {doi:10.1007/BF01053745}.
Disponible à compter du 30 mars 2012 auprès de

http://sprott.physics.wisc.edu/pubs/paper244.htm [SeSm09] Adel S. Sedra et Kenneth C. Smith,

[Ogor97] Maciej J. Ogorzalek, Chaos et complexité dans les circuits Circuits microélectroniques. Oxford, Royaume­Uni : Oxford
électroniques non linéaires. Singapour : World Scientific, University Press, 6e édition, 2009 (6e éd.), 1 456 pages.
1997, 276 pages. {ISBN981­02­2873­2} {ISBN­13 : 978­0195323030}
[Ott93] Edward Ott, Chaos dans les systèmes dynamiques. [Strog94] Steven H. Strogatz, Dynamique non linéaire et
Cambridge, Royaume­Uni : Cambridge University Press, 1993, Chaos : avec des applications à la physique, à la
385 pages. {ISNBN 0­521­43799­7 pbk} biologie, à la chimie et à l'ingénierie. Cambridge MA :
[OtSY94] Edward Ott, Tim Sauer et James A. Yorke (éd.), Faire Westview / Perseus Books Group, 1994, 498 pages.
face au chaos : analyse des données chaotiques et exploitation {ISBN 0­7382­0453­6 pbk}
des systèmes chaotiques. New York, New York : Wiley, 1994, [Take81] Floris Takens, « Détection d'attracteurs étranges en
418 pp. {ISNBN 0­471­02556­9} turbulence », dans DA Rand et L.­S. Young (éd.), Systèmes
[PaCh89] Thomas S. Parker et Léon O. Chua, dynamiques et turbulence, Notes de cours en mathématiques,
Algorithmes numériques pratiques pour les systèmes chaotiques. vol. 898 (Warwick 1980). New York, NY : Springer­Verlag.
New York, NY : Springer , 1989 , 346 pages . {ISBN0­ pp. 366­381, 1981.
387­96688­9} [Take10] F. Takens, « Théorie de la reconstruction et
[PeJS92] H.­O. Peitgen, H. Jurgens et D. Saupe, analyse de séries chronologiques non linéaires », dans [BrHT10],
Chaos et fractales : nouvelles frontières de la science. New York pp. 345­378, 2010.

(NY) : Springer Publishing House, 1992 (1ère éd.), 984 pages ; [Will97] Garnett P. Williams, Théorie du chaos apprivoisée.
2004 (2e éd.), 864 pages. {ISBN0­387­20229­3} Washington, DC : Joseph Henry Press, 1997, 499 p.
[RuTa71] F. Takens, F. et D. Ruelle, « Sur la nature de {ISBN0­309­06351­5}
turbulence », Communications en mathématiques

CEEA Conf. 2012 ; Document 096 ChuaChaos­Full­v40.docx

Winnipeg, MB; 17­20 juin 2012 – 10 sur 10 – 6 avril 2012