Master Stability Function

Machine Translated by Google
Il s'agit du manuscrit accepté mis à disposition via CHORUS. L'article a été publié comme
suit :
Fonctions de stabilité maîtrisées pour une synchronisation complète, intracouche et

intercouche dans les réseaux multiplex d'oscillateurs Rössler couplés Longkun Tang, Xiaoqun
Wu, Jinhu Lü, Junan Lu et
Raissa M. D'Souza Phys. Rév. E 99, 012304 — Publié le 3 janvier 2019 DOI :
10.1103/PhysRevE.99.012304
Fonctions de stabilité maîtrisées pour complet, intracouche et intercouche

synchronisation dans les réseaux multiplex de Rössler couplés
Longkun Tang,1, Xiaoqun Wu,2, 3 Jinhu Lü,4 Junan Lu,2 et Raissa M. D'Souza3
1Laboratoire clé de sciences informatiques de l'Université de la province du Fujian,
École des sciences mathématiques, Université Huaqiao, Quanzhou 362021, Chine.
2École de Mathématiques et Statistiques,
Université de Wuhan, Wuhan 430072, Chine.
3Département d'informatique, Université de Californie, Davis CA 95616, ÉtatsUnis.

4
Institut des sciences des systèmes, Académie des mathématiques et des sciences des systèmes,
Académie chinoise des sciences, Pékin 100190, Chine.
Abstrait
Les phénomènes de synchronisation présentent un grand intérêt dans toutes les disciplines et un intérêt croissant
dans un contexte de réseau multiplex. Pour le réseau multiplex d'oscillateurs Rössler couplés, nous
montrent comment la Master Stability Function, un cadre célèbre pour analyser la synchronisation sur
un seul réseau, peut être étendu à certaines classes de réseaux multiplex avec différents intra
fonctions de couplage couche et intercouche. Nous dérivons trois équations principales de stabilité qui déterminent
respectivement les régions nécessaires de synchronisation complète, de synchronisation intracouche et
synchronisation intercouches. Nous calculons ces trois régions explicitement pour le cas d'un système à deux couches
réseau d'oscillateurs de Rössler et montrent que le chevauchement des régions détermine le type de
synchronisation réalisée. En particulier, si la fonction de couplage inter ou intracouche est telle que
la région de synchronisation intercouche ou intracouche est vide, une synchronisation complète ne peut pas
être atteint quelle que soit la force de couplage. De plus, pour toute structure de réseau, le
l'apparition de synchronisation intracouche et intercouche dépend principalement des fonctions de couplage
de nœuds au sein d’une couche et entre les couches, respectivement. Notre analyse mathématique exige que
les supraLaplaciens intra et intercouches font la navette. Mais nous montrons que ce n’est qu’un moyen suffisant, et non
condition nécessaire et que les résultats puissent être appliqués de manière plus générale.
Motsclés : Réseau multiplex ; fonction de stabilité principale ; synchronisation complète ; synchronisation intracouche
nisation; synchronisation intercouches ; région synchronisée.
* À qui la correspondance doit être adressée : [email protected]
1
INTRODUCTION
La synchronisation dans un réseau d'éléments connectés est essentielle au bon fonctionnement
d'une grande variété de systèmes naturels et artificiels, des réseaux cérébraux à l'énergie électrique
grilles. Cela a stimulé un grand nombre de recherches sur les propriétés de synchronisation.
liens de réseaux complexes, avec des topologies de petit monde, sans échelle et autres types de topologies [1–15].
Pourtant, de nombreux phénomènes de synchronisation, comme dans les réseaux électriques, n’impliquent pas de péché.
gle réseau de manière isolée mais s'appuient sur la synchronisation complète d'un ensemble de plus petits
réseaux. Et plus généralement, audelà des réseaux uniques, on comprend désormais que
les interactions entre les réseaux sont de plus en plus importantes et que les interactions peuvent avoir un impact sur la
processus dynamiques [16–21]. Un paradigme qui capture de nombreux acteurs interdépendants du monde réel
Les réseaux sont ceux des réseaux multiplexés. Ici, le même ensemble de nœuds existe dans plusieurs couches
des réseaux, où chaque couche représente un type d'interaction différent, l'état interne de
les nœuds correspondants dans chaque couche peuvent être distincts et le modèle de connectivité entre
les nœuds de chaque couche peuvent être distincts [22, 23]. À titre d'exemple, considérons le système social en ligne
d'un ensemble d'individus. Ils peuvent interagir sur Twitter ou sur Facebook ou sur Linkedin ou
sur une combinaison des trois, et chaque couche peut avoir son propre modèle de connectivité, mais
il y a généralement une influence qui se propage entre eux [24]. Compte tenu de la nécessité d'étudier la dynamique
processus physiques sur des réseaux complexes en couches et la large applicabilité de la synchronisation,
nous étudions ici les phénomènes de synchronisation sur les réseaux multiplex, un domaine qui a attiré
une attention croissante ces dernières années.
L'une des méthodes les plus importantes pour étudier la synchronisation réseau sur des réseaux uniques
est la méthode de la fonction de stabilité principale (MSF) proposée par Pecora et Carroll [25]. Comme
établi via l'approche MSF, si un réseau peut ou non réaliser la synchronisation
est déterminé non seulement par la structure du réseau, mais aussi par la dynamique nodale et par
la fonction de couplage interne qui décrit les interactions entre les différents composants.
éléments des vecteurs d’état des nœuds connectés [26–28]. En d’autres termes, la dynamique nodale,
la topologie du réseau et la fonction de couplage interne sont trois éléments fondamentaux dans l'étude
synchronisation du réseau. Les deux derniers font l’objet de la plus grande attention, et le premier est
généralement défini comme un système chaotique spécifique. Tels que les circuits de Lorenz, Chen, Chua et
Systèmes Rössler, etc. Ici, le système chaotique de Rössler est sélectionné comme nœud du réseau
dynamique due au fait que le système peut être mis en œuvre par des circuits et appliqué à
2
communication sécurisée. Plus important encore, nous nous concentrons sur les différentes fonctions de couplage interne
au sein et entre les couches dans le cadre du réseau multiplex, ainsi que les différents intracouches
topologies.
Les études actuelles sur les phénomènes de synchronisation dans les réseaux multiplex analysent un multiplex
réseau comme un seul grand réseau composite dont la topologie est décrite par un supra
Matrice laplacienne. Cela nécessite que la fonction de couplage interne soit la même, quel que soit le
si les nœuds sont liés par un bord intracouche ou intercouche et le framework MSF peut
donc être directement appliqué. Les valeurs propres de ce supraLaplacien sont ensuite utilisées pour analyser le
stabilité de l'état de synchronisation complète dans les réseaux multiplex. Par exemple, Solé
Ribalta et coll. [29] ont étudié les propriétés spectrales du Laplacien des réseaux multiplex,
et discuté de la synchronisabilité via le rapport propre de la matrice laplacienne. Aguirre et coll.
Al. [30] ont étudié l'impact du degré de nœud connecteur sur la synchronisabilité de deux réseaux en étoile avec une
liaison intercouches et ont montré que connecter les nœuds de haut degré (bas degré) de chaque réseau est la stratégie
la plus (la moins) efficace pour réaliser la synchronisation. Xu et coll. [31] ont étudié la synchronisabilité des réseaux à
deux couches pour trois couplages spécifiques
modèles, et déterminé qu'il existe une valeur optimale du couplage intercouches
force pour maximiser la synchronisation complète dans les réseaux à deux couches qu’ils analysent.
Li et coll. [32] ont étudié la synchronisabilité d'un réseau duplex composé de deux étoiles
réseaux à deux liens intercouches en donnant une expression analytique contenant le plus grand
et les plus petites valeurs propres non nulles de la matrice laplacienne, le poids du lien, ainsi que le
taille du réseau.
En 2012, Sorrentino et al. [33, 34] ont considéré un modèle « d’hyperréseau » innovant
constitué d’un ensemble de N nœuds qui interagissent via plusieurs types de fonctions de couplage. Note
le contraste avec un réseau multiplex, où un ensemble de N nœuds existe sur chacun des M
couches distinctes (pour un total de M × N nœuds), et chaque nœud peut être dans un état différent dans chaque couche.
(Voir par exemple, Fig. 1.) Dans le modèle « hyperréseau », il n'y a que N nœuds
au total et chaque nœud ne peut être que dans un seul état à un moment donné. À ce titre, l’accent est mis sur
synchronisation complète et trois situations sont trouvées où la topologie du réseau est telle
qu'on peut découpler les effets des fonctions d'interaction de la structure des réseaux
et appliquer l'approche MSF [33, 34]. Très récemment, del Genio et al. étendu cette analyse à un éventail plus large de
scénarios, toujours en utilisant une approche MSF [35], et montrer comment
le modèle « hyperréseau » de [33, 34] est équivalent à un réseau où les nœuds ont de nombreux
3
différents types d’interactions (ou « couches » d’interaction). Bien que ces travaux considèrent que
les nœuds peuvent interagir les uns avec les autres via différentes fonctions de couplage, ils ne capturent pas
la richesse des phénomènes pouvant se produire dans les réseaux multiplex tels que les réseaux intracouches et
synchronisation intercouches.
Jusqu'à présent, seules des études limitées se sont concentrées sur la synchronisation intracouche et intercouche.
Par exemple, Gambuzza et. al [36] ont analysé la synchronisation d'une population d'oscillateurs
couplés indirectement via un milieu inhomogène. Le système est formalisé en termes
d'un réseau à deux couches, où la couche supérieure est composée d'oscillateurs déconnectés, et la
la couche inférieure est constituée d'oscillateurs couplés selon une topologie donnée, et chaque nœud de
la couche supérieure est reliée à son homologue de la couche inférieure. Par simulations numériques,
ils ont montré l'apparition d'une synchronisation intracouche sans cohérence intercouche, ce qui
c'estàdire un état dans lequel les nœuds d'une couche sont synchronisés entre eux sans être
synchronisés avec ceux de l’autre couche. Peu de temps après, SevillaEscoboza et. Al
[37] ont étudié la synchronisation intercouches dans un réseau duplex de couches identiques,
et a montré qu'il existe des cas où chaque nœud d'une couche donnée peut se synchroniser
avec sa réplique dans l'autre couche, que la synchronisation intracouche soit ou non
se produit. Ces découvertes dans des systèmes spécifiques fournissent des bases utiles pour élucider une approche plus
approche fondamentale de l'analyse des phénomènes de synchronisation dans les réseaux multiplexes. Dans
En fait, comme nous le montrerons ici, les équations maîtresses de stabilité peuvent être dérivées pour systématiquement
prédire quand la synchronisation intra et intercouches est simultanément prise en charge et quand
ils ne sont pas simultanément pris en charge pour certaines classes de réseaux multiplex.
Sur la base des motivations cidessus, nous développons ici une méthode Master Stability Function
qui capture une caractéristique essentielle des réseaux multiplex, à savoir le couplage intercouches
Cette fonction peut être distincte de la fonction de couplage intracouche. Ainsi, contrairement aux précédents
Grâce à ces approches, nous pouvons analyser différents types de comportements cohérents, y compris des
synchronisation, synchronisation intracouche et synchronisation intercouches en multiplex
réseaux, cependant, nous sommes limités à certaines classes de topologies. En particulier, nous dé
Rive l'équation maîtresse de stabilité pour un réseau multiplex où le supralaplacien de
les connexions intracouches et celle des connexions intercouches font la navette, telles que définies en détail
cidessous. Nous dérivons en outre deux formes réduites de l'équation maîtresse de stabilité correspondant
aux seules interactions intercouches ou intracouches. Nous montrons ensuite comment trois différents nécessaires
les régions de synchronisation peuvent être calculées à partir du MSF des trois maîtres de stabilité
4
équations. Enfin, nous montrons comment appliquer explicitement le MSF multiplex en analysant un modèle spécifique.
exemple concret d'un réseau à deux couches d'oscillateurs de Rössler avec une topologie intracouche identique
structures et connexions intercouches un à un. Pour une applicabilité plus large de ce multi
approche plex MSF, nous illustrons en outre que les trois équations principales de stabilité peuvent toujours
être utilisé pour prédire la zone de synchronisation pour certaines classes de réseaux multiplex avec
supraLaplaciens non commutatifs.
II. UN CADRE DE FONCTION DE STABILITÉ MAÎTRE POUR LES CLASSES DE
RÉSEAUX MULTIPLEXES
A. Un modèle de réseau multiplex
Nous considérons un réseau multiplex constitué de M couches constituées chacune de N nœuds.

(k) (k) X (k) (k) un
L'état du ième nœud dans la kème couche est spécifié par x = (x ,X
∙∙∙
je je1 , je2 , dans le ),
vecteur d'état à m dimensions. L’évolution du système multiplex complet peut s’écrire comme suit :
N M
(k) (k) (je)
x˙i (k) = f(x je )−c je
je H(xj (k) )−d dklΓ(x je ), je = 1, 2, ∙ ∙ ∙ , N; k = 1, 2, ∙ ∙ ∙ , M, (1)
j=1 l=1
(k)
où x˙i (k) = f(x je ) (je = 1, 2, ∙ ∙ ∙ , N; k = 1, 2, ∙ ∙ ∙ , M) décrit la dynamique isolée
pour le ième nœud de la kème couche, et f(∙) : R m → R m est une fonction vectorielle bien définie,
H(∙) : R m → R m et c sont la fonction de couplage interne et la force de couplage pour les nœuds
au sein de chaque couche, respectivement, et Γ(∙) : R m → R m et d sont la fonction de couplage interne
et la force de couplage pour les nœuds à travers les couches, respectivement. Pour plus de simplicité et de clarté, ici
on laisse H(x) = Hx et Γ(x) = Γx, à savoir les fonctions de couplage entre nœuds sont linéaires
(on peut donc aussi appeler matrices de couplage internes H et Γ). De plus, l'accouplement intérieur
la matrice des nœuds au sein d'une couche, H, est identique pour toutes les couches et la matrice de couplage interne
pour les nœuds répartis sur deux couches, Γ est le même pour toutes les paires de couches.
Éléments l (k)
je décrire la matrice laplacienne des nœuds dans la kème couche. Explicitement, si
(k)
le ième nœud est connecté au jème nœud dans la kème couche, l je = −1, sinon
(k) (k) =−
N (k)
je
je = 0, et l ii
je
j=1 ij , pour je, j = 1, 2, ∙ ∙ ∙ , N et k = 1, 2, ∙ ∙ ∙ , M. De même, si un
le nœud de la kème couche est connecté à sa réplique dans la lème couche, dkl = −1, sinon
M
dkl = 0 et dkk = − l=1 dkl, pour k, l = 1, 2, ∙ ∙ ∙ , M.
Par souci de simplicité, notons
5
(k) (1)
x1
(k) f(x1 ) X
_
˜f(x (1))
(k) (k) (2)
x2
_ f(x 2 ) X ˜f(x (2))
x (k) = .. , ˜f(x (k) ) = .. , X= .. , F(x) = .. ,
. . . .
(k) (M)
XN (k) f(xN ) X ˜f(x (M) )
alors l'évolution du réseau multiplex (équation 1) peut être réécrite comme
L je
x˙ = F(x) − c(L H)x − d(L C)x, (2)
L je
où L représente le supralaplacien des connexions intracouches et L pour le supra
L (1)
L
L (2)
Laplacien des connexions intercouches. En détail, L = Ml=1L (k) =
..
.
L (M)
je je
et moi =L DANS . Voici l'opération de somme directe, IN est la matrice d'identité N par N,
(k) est l'opération produit de Kronecker, L (k) = (l je )N×N est la matrice laplacienne des nœuds
je
dans la kième couche, et L = (dkl)M×M représente la matrice laplacienne intercouches. Plus
des détails sur les supraLaplaciens et les modèles de réseau multiplex peuvent être trouvés dans les réf. [18, 23, 29,
31, 32] et les références qui y figurent.
B. Trois équations maîtresses de stabilité
La méthode de la fonction de stabilité principale [25] est l'une des méthodes les plus importantes pour étudier
stabilité des systèmes identiques couplés synchronisés. Il simplifie un système en réseau à grande échelle.
à un système de taille de nœud via la diagonalisation et le découplage, tant que le couplage interne
les fonctions pour toutes les paires de nœuds sont identiques. Ainsi, déterminer si un réseau peut atteindre
la synchronisation peut être transformée en déterminant si tous les modes caractéristiques du réseau
tomber dans les régions synchronisées correspondantes. Dans ce qui suit, nous établirons un mas
cadre de stabilité ter pour les réseaux multiplex avec intercouche et intracouche non identiques
fonctions de couplage interne.
Selon l’idée du cadre de stabilité principal [25], pour étudier la synchronisation du réseau, nous pouvons linéariser
l’équation dynamique (2) à 1M 1N s, où s est un état synchrone du réseau satisfaisant s˙ = f(s ) et 1M désigne
un Mdimensionnel
6
vecteur dont toutes les entrées sont 1. On obtient ainsi l’équation variationnelle suivante :
L je
˙ξ = [IM×N Df(s) − c(L H) − d(L Γ)]ξ, (3)
où ξ = x − 1M 1N s et IM×N est la matrice identité d'ordre M × N.

L Moi et L
Supposons que L sont des matrices symétriques et satisfont L LL I = L IL L . Après
diagonalisation et découplage (voir l'annexe B pour plus de détails), nous obtenons l'équation principale de stabilité du
multiplex pour un système décrit par l'équation. (1):
y˙ = [Df(s) − αH − βΓ]y, (4)
L Moi et L
où α = cλ, β = dµ, λ et µ sont les valeurs propres de L respectivement, et satisfaire
2 + m = 0.
2
minutes
Puisque cette équation peut être un système variable dans le temps, en particulier si s(t) est fonction de
temps, ses valeurs propres peuvent ne pas être utiles pour déterminer la stabilité. Par conséquent, le plus grand
Exposant de Lyapunov (LLE) de l'équation. (4) est utilisé à la place, qui est une fonction de α et β, noté σ (α, β) et appelé
fonction de stabilité principale multiplex pour l'équation. (1). Veuillez consulter l'annexe
A pour plus d'informations sur les exposants de Lyaponuv.
Lorsque λ = 0 et µ = 0, il n'y a pas de couplages intercouches quel que soit d, pour d arbitrairement
choisi dans [0, +∞), et Eq. (4) se réduit à
y = [Df(s) − αH]y, (5)
Il est clair que l’équation. (5) devient exactement l'équation principale de stabilité de chaque réseau intracouche
indépendant (pas de couplages intercouches).
De même, lorsque λ = 0 et µ = 0, on peut obtenir l'équation suivante
y˙ = [Df(s) − βΓ]y, (6)
quelle que soit la force de couplage c, pour c arbitrairement choisi dans [0, +∞). Éq. (6) devient exactement l'équation
principale de stabilité pour chaque réseau intercouche indépendant (pas de couplages intracouches).
Pour un réseau monocouche, une condition nécessaire pour que le collecteur de synchronisation soit
stable est que le plus grand exposant de Lyapunov σ(α) de l'équation. (5) inférieur à zéro [38]. Par analogie avec une
seule couche, pour l'équation principale de stabilité du multiplex (4), σ(α, β) < 0 est une valeur nécessaire
condition de stabilité du collecteur de synchronisation dans un réseau multiplex.
7
Il convient de noter en particulier le cas où les fonctions de couplage intra et intercouches
sont identiques. Ici H = Γ, et Eq. (4) se transforme en y˙ = [Df(s) − γH]y, (avec γ = α + β).
Il s'agit de l'équation principale de stabilité pour le réseau composite unique correspondant où
le modèle de couplage interne (fonction) entre deux nœuds quelconques est identique, et un seul supra
Le laplacien peut décrire sa topologie. C’estàdire que l’équation maîtresse de stabilité du
le réseau composite unique est un cas particulier de l'équation. (4).

L
L'hypothèse selon laquelle L Moi et L sont symétriques et satisfont aux conditions de commutativité.
tion est une condition importante pour le découplage du système et restreint notre approche de
s'appliquant à l'ensemble de la classe des réseaux multiplex. Mais cette hypothèse peut être assouplie, car
L
c'est seulement une condition suffisante mais non nécessaire. Considérons d’abord le cas où L
je
et moi sont commutatifs mais non symétriques. Comme le montre l'Annexe E, le même maître
les équations de stabilité (4) à (6) (qui correspondent respectivement aux équations (E10), (E9) et (E6) peuvent
être dérivé à condition que le réseau multiplex ait une topologie intracouche identique sur
L
chaque couche et que la matrice laplacienne intracouche L et le Laplacien intercouche
matrix L je
peut être diagonalisable et avoir de réelles valeurs propres. Il y a des cours importants
de réseaux du monde réel qui correspondent à ce paradigme, comme le système de référence en fonctionnement continu.
l’infrastructure d’information géospatiale des stations de détection (CORS) [39–41] discutée en détail dans
Annexe F.
L
Considérons ensuite le cas où L Moi et L ne faites pas la navette. Comme le montre la simulation
résultats avec un réseau de Rössler à deux couches avec des supraLaplaciens non commutatifs, les trois
Les équations principales de stabilité (4) à (6) peuvent toujours être utilisées pour prédire le comportement de synchronisation du réseau.
iors. En particulier, pour les réseaux duplex, si la topologie du réseau est différente sur chaque couche,
mais il existe un couplage pondéré identique un à un des nœuds entre les couches, nous pouvons prédire
synchronisation complète et synchronisation intracouche. Si la topologie sur chaque couche est
identique, mais le couplage pondéré un à un n'est pas identique, nous pouvons prédire une
synchronisation et synchronisation intercouches. (Voir la partie simulation pour plus de détails.)
C. Régions synchronisées
En utilisant les équations maîtresses de stabilité multiplex développées cidessus, nous pouvons analyser trois
types de comportements de synchronisation : synchronisation complète, synchronisation intracouche
et synchronisation intercouches. Ici, nous définissons les régions qui prennent en charge chaque comportement et
8
dans les sections suivantes et dans le matériel supplémentaire, nous montrons que c'est le chevauchement de ces éléments
régions qui déterminent le type de modèle de synchronisation affiché par un réseau multiplex.
Pour le réseau multiplex complet, à partir de l'équation principale de stabilité du multiplex (4), nous pouvons
calculer la région
Rα,β = {(α, β)|σ(α, β) < 0},
qui est appelée région synchronisée commune (qui prend en charge la synchronisation complète du réseau).
Chaque fois que σ(α, β) < 0, perturbations transversales à la variété de synchronisation
s'éteignent et le réseau est dit synchronisable.
De l’équation. (5), nous pouvons obtenir la région pour la synchronisation intracouche. La région dépend
uniquement sur la valeur du paramètre α, mais pour permettre ultérieurement une comparaison sur l'ensemble du paramètre
espace, nous incluons explicitement le paramètre β dans la définition de la région,
Intra
R = {(une, b)| σ(a) < 0},
a, b
où LLE(α) est le plus grand exposant de Lyapunov pour l'équation principale de stabilité (5). De même, à partir de
l’équation. (6), nous obtenons la région de synchronisation intercouches
Inter
R = {(une, b)|σ(b) < 0}.
a, b
Nous appelons ces régions dans l’espace des paramètres de (α ≥ 0, β ≥ 0) la synchronisation correspondante.
régions nisées par rapport à α et β.
Lorsque les structures topologiques du réseau sont spécifiées, nous pouvons déterminer λ et µ (les
L
valeurs propres de L Moi et L ) directement, et alors les régions Rα,β, RIntra et RInter peuvent être α,β
un B
paramétré simplement en termes de forces de couplage c et d, notées Rc,d, RIntra RInter c,d . Nous et
CD
appelons ces régions les régions synchronisées correspondantes par rapport aux cou
plings c et d.
¨
III. RÉSEAU À DEUX COUCHES D'OSCILLATEURS ROSSLER
Avec le développement du cadre MSF multiplex, nous analysons maintenant plus en profondeur un problème spécifique
exemple d'un réseau à deux couches d'oscillateurs Rössler et calculer les différents types de
régions synchronisées.
9
Le célèbre oscillateur chaotique de Rössler est décrit comme
x˙ = −y − z,
y˙ = x + ay, (7)
z˙ = z(x − c) + b ,
où a = b = 0,2 et c = 9. C'est la fonction f dans le réseau multiplex Eq. (1). Que
c'estàdire que l'état de chaque nœud du réseau est un vecteur tridimensionnel avec chaque composant
évoluant par l'équation. (7). Pour les matrices de couplage interne H et Γ, nous considérons la famille de choix (où i, j = 1, 2, 3),
3×3
qui correspondent à la forme la plus simple Iij qui représente une matrice dont
R (i, j) l'élément est un et les autres éléments sont nuls. La topologie intercouches est définie comme étant une
connexion à un, c'estàdire que chaque nœud d'une couche est connecté à un nœud homologue
dans l'autre couche.
Ensuite, nous donnons d'abord un aperçu de la synchronisation intracouche et intercouche, puis
calculer les régions paramétrées de synchronisation pour l'intracouche générale (inconnue)
topologies et montrer comment déterminer les régions synchronisées après avoir spécifié les
topologies de couches dans la dernière soussection.
A. Synchronisation intracouche et intercouche
Il est bien connu qu'une synchronisation complète signifie que tous les nœuds d'un réseau arrivent à
un état identique. Mais pour les réseaux multiplex, il est également très important d'étudier les réseaux intracouches.
synchronisation et synchronisation intercouches. Comme le montre la figure 1, la synchronisation intracouche
La nisation signifie que tous les nœuds de chaque couche atteignent un état identique, tandis que les nœuds intercouches
la synchronisation signifie que chaque nœud d'une couche atteint le même état que ses homologues dans
d'autres couches.
B. Régions synchronisées pour les topologies intracouches inconnues
Les régions de synchronisation calculées à partir du multiplex MSF sont paramétrées par α
et β, et ne nécessitent donc pas que la topologie inter et intracouche soit spécifiée. Figure 2
montre les régions synchronisées telles que paramétrées par (α, β) pour un réseau multiplex à deux couches
d'oscillateurs Rössler avec une topologie arbitraire pour différentes combinaisons d'intercouches et
matrices de couplage intracouche H et Γ. Ici, l'ombrage vert (gris) représente les régions
dix
(un) (b)
d d
c c
FIGUE. 1. Représentation schématique de (a) la synchronisation intracouche et (b) la synchronisation intercouche
nisation, dans un réseau multiplex de deux couches.
Rα,β tel qu'obtenu à partir de l'équation principale de stabilité (4). Les régions RIntra telles un
qu'obtenues
B à partir de l'Eq.
(5) sont entourés de lignes bleues en pointillés et les régions RInter obtenues à partir de l'équation.
un B
(6) entouré de
lignes pointillées rouges.
La synchronisation se produit dans la région lorsque le critère MSF est négatif, c'estàdire
lorsque σ(α, β) < 0. Ainsi, à partir de la figure 2, nous pouvons facilement obtenir la région synchronisée conjointe :
Rα,β ≈ {(α, β)| 0,2 < α + β < 4,6} pour H = I11 et Γ = I11,
Rα,β ≈ {(α, β)| 0,23 < α < 4,3, β ≥ 0} pour H = I11 et Γ = I13,
un β
Ra,β ≈ {(α, β)| + 0,2 0,18 > 1, β > h(α)} pour H = I11 et Γ = I22
Rα,β ≈ {(α, β)| β > 0,2, α ≥ 0} pour H = I13 et Γ = I22,
3 − 0,012a 2
+ 0,12a − 0,35.
où h(α) = −10−5α 4 + 0,00057α
En particulier, en laissant β = 0 dans Rα,β, nous avons l'intervalle Rα = (0,2, 4,6) pour H = I11, Rα =
pour H = I13 et Rα = (0,18, ∞) pour H = I22. De même, en laissant α = 0 dans Rα,β, nous avons Rβ = (0,2, 4,6)
pour Γ = I11, Rβ = pour Γ = I13 et Rβ = (0,18, ∞) pour Γ = I22. Ici, les intervalles Rα {α| σ(α) < 0} et Rβ {β|
σ(β) < 0} peut également être obtenu à partir des équations. (5) et (6), respectivement.
Par conséquent, pour trois types de modèles de couplage, c'estàdire H = I11 et tout Γ, H = I13 et
tout Γ, et H = I22 et tout Γ, nous obtenons la région suivante pour la synchronisation intracouche,
11
5 dix
4.6 (un) (b)
4 8
3 6
b
b
2 4
1 2
0,2
0
0,2 1 2 3 4 4,6 5 0,2 1 2 3 4 4,6 5 6
un un
2 2
(c) (d)
1,5
b
1 1
0,5
0,2
0
0,2 2 4.6 6 8 0 2 4 6 8 dix
un un
FIGUE. 2. Les régions synchronisées par rapport à α et β, Rα,β peintes en vert (gris),
Entrer
un B
entouré par les lignes bleues en pointillés, et RInter un B
délimité par des lignes pointillées rouges. Ici
l'oscillateur Rössler est pris comme dynamique nodale, et la matrice de couplage intracouche H et la
les matrices de couplage intercouches Γ sont choisies comme suit : (a) H = I11, Γ = I11, (b) H = I11, Γ = I13,
(c) H = I11, Γ = I22, (d) H = I13, Γ = I22.
respectivement.
R.Intra = {(une, b)| 0,2 < α < 4,6, β ≥ 0}, RIntra = ,

un B un B
et
R.Intra = {(une, b)| 0,18 < α < +∞, β ≥ 0}.

un B
De manière analogue, on peut obtenir RInter

un B en remplaçant α par β, et H par Γ dans le RIntra cidessus un B .
Comme le montre la figure 2, pour H = I11 et Γ = I11, les régions
R.Intra = {(une, b)| 0,2 < α < 4,6, β ≥ 0} et R

Inter
= {(une, b)| 0,2 < β < 4,6, α ≥ 0},
un B un B
12
qui sont entourés respectivement par des lignes pointillées bleues et rouges.
Pour H = I11 et Γ = I13,
R. Intra = {(une, b)| 0,2 < α < 4,6, β ≥ 0} et R

Inter
= ,
un B un B
où RIntra est entouré de lignes bleues en pointillés.

un B
Pour H = I11 et Γ = I22, RIntra est la partie délimitée par les lignes bleues en pointillés, et
un B
RInter est celui audessus de la ligne pointillée rouge. Pour être exact,
un B
R. Intra = {(une, b)| 0,2 < α < 4,6, β ≥ 0} et R

Inter
= {(une, b)| 0,18 < β < +∞, α ≥ 0}.
un B un B
Pour H = I13 et Γ = I22, RIntra = , et RInter = {(une, b)| 0,18 < β < +∞, α ≥ 0},
un B un B
qui est audessus de la ligne pointillée rouge.
D'une manière générale, un réseau multiplex avec une topologie spécifiée peut réaliser des
synchronisation lorsque tous les modes caractéristiques du réseau non nuls, y compris ceux du
Les Laplaciens intracouches et intercouches tombent dans la région synchronisée. Pour un réseau à deux couches
fonctionnent avec des topologies intracouches identiques, notre analyse théorique (voir Annexe E) approfondit
montre qu'un réseau duplex peut réaliser une synchronisation complète lorsque tous les
et RInterα, β . Par conséquent, accordez

les modes caractéristiques se situent à l’intersection de Rα,β, RIntra un B
En ce qui concerne la région qui se chevauche, on peut déterminer si le réseau atteint un niveau complet
synchronisation ou non après avoir spécifié la topologie. Cependant, que se passetil lorsque tous
les modes caractéristiques non nuls ne tombent pas dans l'intersection ? D'autres simulations montrent
que dans ce cas, le réseau pourrait prendre en charge d'autres comportements dynamiques cohérents, tels que
synchronisation intracouche ou intercouche.
C. Régions synchronisées avec des topologies intracouches données
Pour pousser l'analyse plus loin, il faut préciser la topologie du Rössler à deux couches
réseau d'oscillateurs. Pour plus de simplicité, supposons que les deux couches ont le même intracouche
topologie, et chaque nœud d'une couche est connecté à sa réplique dans l'autre couche. Considérer
que chaque couche est un réseau en étoile composé de 5 nœuds. Ensuite, la matrice laplacienne intracouche
13
est
4 −1 −1 −1 −1
−1 1 0 0 0
L= −1 0 1 00 ,
−1 0 0 1 0
−1 0 0 0 1
L0
et la matrice supralaplacienne intracouche est L L = Le Laplacien intercouche
0 litre
.
1 −1
matrix L je = et la matrice supralaplacienne intercouche L je
=L je
I5. Il
−1 1
,
est facile de vérifier que L LL I = L IL L , L et les valeurs caractéristiques de L Moi et L sont λ =
0, 0, 1, 1, 1, 1, 1, 1, 5, 5 et µ = 0, 0, 0, 0, 0, 2, 2, 2, 2, 2, respectivement.
Nous pouvons calculer directement les valeurs propres λ et µ et paramétrer les valeurs propres synchronisées.
régimes par les forces de couplage c et d (plutôt que les plus générales α et β) pour toutes les
différentes combinaisons des matrices de couplage internes H et Γ. (Voir l'Annexe B pour plus
de détails sur la transformation de Rα,β en Rc.d.) Par conséquent, pour H = I11 et Γ = I11, la région
par rapport aux paramètres c et d est
Rc,d ≈ {(c, d)| 0,2 < c + 2d, c + 0,4d < 0,92}.
De même, pour H = I11 et Γ = I13, alors
Rc,d ≈ {(c, d)| 0,23 < c < 0,86, d ≥ 0} ;
pour H = I11 et Γ = I22, alors
c d
Rc,d ≈ {(c, d)| + 0,2
0,09
> 1, d > h(c)}
où 12
4 3
(−625 ∙ 10−5c + 0,07125c − 0,2c
2
+ 0,6c − 0,35) ; et pour H = I13 et Γ = I22, alors Rc,d
≈ {(c, d)| d > 0,1, c ≥ 0}. Ces régions, Rc, d, sont représentées sur les Fig. 36 par les lignes pleines
dans les panneaux (c) pour les différents choix de H et Γ considérés.
Pour tester nos prédictions théoriques, nous résolvons ensuite numériquement le duplex de Rössler en réseau
système et identifier les régions de paramètres qui prennent en charge les trois comportements cohérents différents.
iors : synchronisation complète, synchronisation intracouche et synchronisation intercouche.
14
0,2 5
(un) 4.6 (b)
0 4
3
−0,2
b
2
−0,4
1
−0,6
0,2
0,2 2 3 4,6 6 8 9 10 0,2 1 2 3 4 4,6 5
un B un
2.5 2.5
(c) (d)
2 2
1,5 1,5
d
d
1
1 1 1
3
3
0,5 0,5
0,1 2 0,1 2
0 0,2 0,4 0,6 0,8 1 0 0,2 0,4 0,6 0,8 1
c c
FIGUE. 3. Régions synchronisées sur le réseau pour H = I11 et Γ = I11. (a) l'intervalle synchronisé de
le réseau Rössler intracouche/intercouche indépendant par rapport à α/β ; (b) le synchronisé
région par rapport à α et β pour les réseaux de Rössler ; (c) la région synchronisée par rapport à
couplages c et d pour un duplex de Rössler constitué de deux couches en étoile avec une couche intermédiaire unàun
Connexions; (d) des zones de synchronisation numérique par rapport aux couplages c et d, dans lesquelles le
La région marron (gris foncé) représente la zone de synchronisation complète, le jaune (gris moyen)
est pour la synchronisation intracouche, et le cyan (gris clair) est la synchronisation intercouche et le
la région blanche représente la nonsynchronisation.
Nous quantifions que le système a atteint le type de comportement spécifique via la synchronisation.
erreurs de tion telles que définies à l'Annexe D. En limitant les valeurs de ces erreurs, nous développons
trois fonctions d'indicateur différentes, qui identifient que le système a atteint un niveau macroscopique
ordre de la forme : Id = 3 lorsque le réseau atteint une synchronisation complète, Id = 2 pour
synchronisation intracouche, Id = 1 pour la synchronisation intercouche et Id = 0 pour aucun des
cas cidessus. Voir l'annexe D pour plus de détails.
15
dix
0,2 (un) (b)
8
0
6
−0,2
b
4
−0,4
H=I11
2
−0,6 C=I13
0
0,2 1 2 3 4,6 6 7 8 9 10 0,2 1 2 3 4 4,6 5
un B un
3 3
(c) (d)
2.5 2.5
2 2
d
d
1,5 2 1,5
2
1 1
0,5 0,5
0 0
0 0,23 0,5 0,86 1 1.2 0 0,2 0,4 0,6 0,8 1
c c
La figure 3 montre les régions synchronisées du réseau pour le réseau en étoile à deux couches de Rössler.
oscillateurs pour le scénario H = I11 et Γ = I11. En détail, le panneau (a) affiche la synchronisation.
intervalles nisés des réseaux Rössler intracouches et intercouches indépendants par rapport
à α ou β, qui peut être calculé à partir des équations principales de stabilité (5) et (6) (sans
considération de d ou c), respectivement. Puisque H = Γ, les deux intervalles se chevauchent. Panneau (b)
donne la région synchronisée par rapport à α et β pour ce réseau de Rössler calculé
16
à partir de l’équation principale de stabilité (4). Le panneau (c) montre la région synchronisée en fonction de la
force de couplage intra et intercouches c et d. Le panneau (d) montre le nombre numérique
fonction indicatrice calculée (c'estàdire les valeurs calculées numériquement de l'erreur de synchronisation
telles que classées dans l'équation (D4) donnée à l'annexe D) par rapport aux couplages c et d pour
ce réseau Rössler duplex. Ici, la zone marron (gris foncé) marquée « 3 » représente
synchronisation complète, la zone jaune (gris moyen) étiquetée « 2 » est destinée à l'intracouche
synchronisation, l'étiquette cyan (gris clair) avec « 1 » est pour la synchronisation intercouches, et
la région blanche représente les cas contraires.
D'autres choix pour les fonctions de couplage H et Γ sont présentés dans les figures 4 à 6 pour ce même
réseau en étoile à deux couches d'oscillateurs Rössler. Les résultats sont analogues à ceux de la figure 3. Il
Il convient de noter que dans les panneaux (c) de toutes ces figures, les régions de couche complète intracouche
et la synchronisation intercouches prédite par le multiplex MSF Eqs. (4) à (6), représenté par
les régions marron (gris foncé), jaune (gris moyen) et cyan (gris clair), respectivement,
peut capturer tous les comportements présentés par des simulations numériques directes présentées dans des panneaux
(d).
Nous montrons ensuite comment ces zones distinctes peuvent être déterminées à partir des trois régions : Rc,d,
Entrer
c,d ,
et RInter dérivé des équations. (4), (5) et (6). En effet, les intersections
CD
des régions déterminent le type de comportement cohérent qui est stable. Plus précisément, le
l'intersection des trois régions détermine une synchronisation complète, l'intersection
de Rc,d et RIntra CD détermine la synchronisation intracouche et l'intersection de Rc,d et
RInter
CD
détermine la synchronisation intercouches.
Par exemple, pour le cas où H = I11 et Γ = I11, la région synchronisée Rc,d =

=
{(c, ré)| c + 2d > 0,2, c + 0,4d < 0,92}, la région synchronisée intracouche RIntra CD
{(c, ré)| 0,2 < c < 0,92, d ≥ 0} et la région synchronisée intercouches RInter CD
= {(c, ré)|c ≥
0, 0,1 < d < 2,3 }. L'intersection de ces trois parties est {(c, d)| c > 0,2, d > 0,1, c+0,4d <
0,92}, comme indiqué par le numéro « 3 » dans le panneau (c) de la figure 3, qui coïncide essentiellement avec le
Zone de synchronisation complète calculée numériquement en couleur marron (gris foncé) dans le panneau
(d). De plus, la simple synchronisation intracouche (sans synchronisation intercouche)
la zone en jaune (gris moyen) dans le panneau (d) coïncide avec la région étiquetée « 2 » dans le panneau
(c) : Rc,d ∩ RIntra − RInter = {(c, ré)| c > 0,2, 0 ≤ d < 0,1, c + 0,4d < 0,92}, et le simple
CD CD
zone de synchronisation intercouches (sans synchronisation intracouche) en cyan (gris clair)
s'accorde bien avec la région étiquetée « 1 » dans le panneau (c) : Rc,d ∩ RInter Entrez = {(c, ré)| 0 ≤
CD CD
17
0,3 2
(un) (b)
0 1,5
b
1
−0,5 H=I11
C=I22 0,5
−1 0
0,18 1 2 3 4,6 6 7 8 9 10 0,2 2 4.6 6 8
un B un
1,5 1,5
(c) (d)
1 3
1 1
1
1
d
d
3
0,5 0,5
0,1 0,1 2
0 0
2 0,2 0,92 2 3 4 0,2 0,92 2 3 4
c c
c < 0,2, d > 0,1, c + 0,4d < 0,92}. Des observations similaires peuvent être obtenues dans les panneaux (c) et
(d) des Fig. 46.
En d’autres termes, la zone réelle pour une synchronisation complète est déterminée par l’interc’està
section de Rc,d, RIntra

CD
et RInter
c,d , dire Rc,d∩RIntra CD
∩RInter
c,d . De plus, le simple intracouche
la zone de synchronisation est déterminée par l'intersection de la région synchronisée et de l'intracouche
région synchronisée soustrayant la partie synchronisée intercouches, c'estàdire Rc,d∩RIntraCD −RInter

c,d .
18
2
(un) (b)
0,2
−0,2
b
1
−0,4
H=I13
−0,6
C=I22
−0,8
0,2
−1
0,18 1 2 3 4 5 6 7 8 9 10 0 2 4 6 8 dix
un B un
2.5 2.5
(c) (d)
2 2
1,5 1,5 1
1
d
d
1 1
0,5 0,5
0,1 0,1
0 1 2 3 4 0 1 2 3 4
c c
La simple zone de synchronisation intercouches est déterminée par Rc,d ∩ RInter Entrez .
CD CD
De plus, lorsque la dynamique nodale et les structures de réseau sont données, RIntra CD
et RInterCD
sont principalement déterminés par les matrices de couplage internes des nœuds intracouches (H) et
les nœuds intercouches (Γ) respectivement, et Rc,d est déterminé par les deux. En particulier, si le
la matrice de couplage intercouches Γ rend la région synchronisée intercouches RInter CD

vide, alors
le réseau multiplex ne peut pas réaliser la synchronisation intercouches, ce qui entraîne l'échec de
19
2.5 2.5
(un) (b)
2 2
1,5 1,5
d
d
1 1
1 1
3 3
0,5 0,5
0,1 2 0,1 2
0,04 0,2 0,4 0,6 0,921 0,04 0,2 0,4 0,6 0,8 0,921
c c
2 2
(c) (d)
1,5
1
1
d
d
1 3 1 1 3
0,5
0,1 2 0,1 2
0,04 0,2 0,4 0,6 0,8 1 1.2 0,04 0,2 0,4 0,6 0,92 1.2
c c
FIGUE. 7. Régions synchronisées en réseau pour les réseaux Rössler composés de deux monocouches entièrement
réseaux connectés avec différents H et Γ, H = I11, Γ = I11 pour (a) et (b), et H = I11, Γ = I22
(c) et (d). (a) et (c) sont les régions synchronisées autour de c et d ; (b) et (d) sont respectivement
zones de synchronisation numérique correspondantes, dans lesquelles la région marron (gris foncé) représente
zone de synchronisation complète, le jaune (gris moyen) est pour la synchronisation intracouche, et le
le cyan (gris clair) correspond à la synchronisation intercouches et la région blanche représente la nonsynchronisation.
synchronisation complète, comme le montre la figure 4. Si la matrice de couplage intracouche H fait
la région synchronisée intracouche RIntra CD

vide, alors le réseau multiplex ne peut pas atteindre
synchronisation intracouche, ce qui conduit également à l'échec de la synchronisation complète, comme indiqué
sur la figure 6.
Afin de vérifier les résultats précédents sur une topologie multiplex différente, nous considérons un
réseau duplex composé de deux couches réseau entièrement connectées avec une couche intermédiaire unàun
Connexions. Les résultats présentés dans les Fig. Les figures 7 et 8 illustrent à nouveau les observations cidessus.
20
3 3
(un) (b)
2.5 2.5
2 2
2
2
d
d
1,5 1,5
1 1
0,5 0,5
0
0 0,05 0,5 0,86 1 1.2 0 0,2 0,4 0,6 0,8 1
c c
4 4
(c) (d)
1
3 3
1
d
d
2 2
1 1
0,1 0,1
0 0,2 0,4 0,6 0,8 1 0 0,2 0,4 0,6 0,8 1
c c
FIGUE. 8. Régions synchronisées en réseau pour les réseaux Rössler composés de deux monocouches entièrement
réseaux connectés avec différents H et Γ, H = I11, Γ = I13 pour (a) et (b), et H = I13, Γ = I22
(c) et (d). (a) et (c) sont les régions synchronisées autour de c et d ; (b) et (d) sont respectivement
zones de synchronisation numérique correspondantes, dans lesquelles la région marron (gris foncé) signifie simplement
synchronisation complète, la région jaune (gris moyen) signifie une simple synchronisation intracouche,
la région cyan (gris clair) signifie la synchronisation intercouches et la région bleue signifie la nonsynchronisation entre les couches.
synchronisation.
¨
IV. RÉSEAU ROSSLER À DEUX COUCHES AVEC SUPRA
LAPLACIENS
Jusqu'à présent, nous avons analysé le cas des supraLaplaciens commutatifs avec lesquels nous dérivons
les trois équations principales de stabilité (4)(6). Cependant, la condition de commutativité restreint
notre approche s'applique à l'ensemble de la classe des réseaux multiplex, ce n'est qu'une solution suffisante mais
ce n'est pas une condition nécessaire et peut être assoupli.
21
2.5 2.5
(un) (b)
2 2
d 1,5 1,5
d
1
1 1 1
3
3
0,5 0,5
0,1 2 2
0
0 0,2 0,4 0,6 0,8 1 0 0,2 0,4 0,6 0,8 1
c c
2.5 2.5
(c) (d)
2 2
1,5 1,5
1
d
d
1 1
3
3
0,5 0,5
0 2 0 2
0 0,2 0,4 0,6 0,8 1 0 0,2 0,4 0,6 0,8 1
c c
FIGUE. 9. Le cas des matrices supralaplaciennes non commutatives avec différentes topologies intracouches et
poids de couplage intercouches identiques. Régions synchronisées sur le réseau calculées à partir des équations de stabilité principales
(a) et les zones de synchronisation numérique (b)(c)(d) pour H = I11 et Γ = I11. La première couche du
le réseau duplex est de type étoile, la deuxième couche est celle générée à partir du type étoile avec 1 (b), 2 (c)
et 3 (d) bords supplémentaires, respectivement. Le couplage un à un entre les couches est identique.
Nous considérons ici deux cas de supraLaplaciens non commutatifs pour montrer que le
trois équations principales de stabilité peuvent toujours être utilisées pour prédire la synchronisation du réseau
comportements. l'un est un réseau duplex qui a une topologie différente sur chaque couche et un à un
couplage pondéré identique des nœuds entre les couches. L'autre est un réseau duplex qui
a une topologie identique sur chaque couche et un couplage pondéré un à un non identique des nœuds
entre les couches.
Pour le premier cas, considérons des réseaux duplex spécifiques avec 5 nœuds sur chaque couche et un
couplage à un de nœuds entre couches, où une couche est de type étoile et l'autre
est du type étoile avec 1, 2 ou 3 arêtes supplémentaires. Dans ce cas, il est facile de vérifier que
L1 0
matrices supralaplaciennes intra et intercouches L
L = Moi et L = L je
DANS ne pas
0 L2
22
2 2
(un) (b)
1,5 1,5
3
1
3 1
d
d
1 1 1
0,5 0,5
0,1 2
0 0 2
0,2 0,92 2 3 0 1 2 3
c c
2 2
(c) (d)
1,5 1,5
1
3 1
3
d
d
1 1
0,5 0,5
0 2 0 2
0 1 2 3 0 1 2 3
c c
FIGUE. 10. Le cas des matrices supralaplaciennes non commutatives avec différentes topologies intracouches et
poids de couplage intercouches identiques. Régions synchronisées sur le réseau calculées à partir des équations de stabilité principales
(a) et les zones de synchronisation numérique (b)(c)(d) pour H = I11 et Γ = I22. La première couche du
le réseau duplex est de type étoile, la deuxième couche est celle générée à partir du type étoile avec 1 (b), 2 (c)
et 3 (d) bords supplémentaires, respectivement. Le couplage un à un entre les couches est identique.
1 −1
Je fais la navette. Ici L = et les plus petites valeurs propres non nulles des deux intra
−1 1
,
les matrices lapaciennes de couche L1 et L2 sont égales et leurs plus grandes valeurs propres sont également égales, c'estàdire
λ2 = 1 et λN = 5.
Pour le deuxième cas, considérons des réseaux duplex spécifiques qui ont un type d'étoile ou un réseau identique.
topologie entièrement connectée sur chaque couche et couplage unàun pondéré non identique entre
entre couches, la matrice supralapalacienne intercouches étant L je = L je diag{2, 1, 1, 1, 1},

1 −1
où ici L je =
Dans ce cas, Lje la plus petite valeur propre non nulle de µ2 = 2 et sa
−1 1
.
plus grande valeur propre µN = 4. Il est facile de vérifier que LL et moije ne faites pas la navette.
Les figures 9 à 12 montrent les résultats pour les deux classes différentes de réseaux duplex cidessus avec
différentes combinaisons de H et Γ. Nous constatons toujours que les régions superposées obtenues à partir de
23
les trois équations principales de stabilité coïncident étroitement avec les zones calculées numériquement pour
les trois différents types de comportements synchrones. Plus précisément, la zone réelle pour compléter
la synchronisation est déterminée par Rc,d ∩ RIntra CD ∩ RInter pour les deux classes de noncommutatif
c,d supraLaplaciens. Pour les réseaux duplex avec différentes topologies intracouches (la première classe),
la zone de synchronisation intracouche est déterminée par Rc,d∩RIntra c,d . Pour les réseaux duplex avec
couplage unàun pondéré non identique (la deuxième classe), la synchronisation intercouches
l’aire est déterminée par Rc,d ∩ RInter c,d . Ces résultats mettent en lumière les faits importants qui
la différence des topologies intracouches peut conduire à la modification de l'intercouche réelle
les régions synchronisées et les poids de couplage unàun intercouches non identiques peuvent conduire à
le changement des régions synchronisées intracouches réelles.
En d’autres termes, même si ici les matrices supralaplaciennes inter et intracouches ne
ne font pas la navette, les trois régions synchronisées prédisent toujours les zones réelles pour une synchronisation complète.
chronisation et synchronisation intracouche, ou pour une synchronisation complète et intercouche
synchronisation. Par conséquent, la condition de commutation n’est pas nécessaire pour nos résultats, elle
n’est suffisant que pour notre analyse théorique. Notamment pour le cas de différents intracouches
topologies, on peut appliquer ces trois régions synchronisées pour prédire les zones réelles pour
synchronisation complète et synchronisation intracouche. Dans quelle mesure l'observation
s’applique reste une question ouverte.
¨
V. UN RÉSEAU À TROIS COUCHES D'OSCILLATEURS ROSSLER
Considérons ici un réseau à trois couches d'oscillateurs Rössler avec une topologie interne identique
(comme la structure entièrement connectée) et le couplage de type chaîne entre les couches. C'est le
L
matrice supralaplacienne intracouche L = I3 L, et la matrice supralaplacienne intercouches
1 −1 0
L je
=
L =L DANS , Je où L −1 2 −1 .
0 −1 1
Comme le montrent les fig. 13 et 14, les trois régions synchronisées calculées à partir des trois
Les équations maîtresses de stabilité peuvent également prédire les zones réelles des trois phases synchrones.
comportements, vérifiant en outre que nos résultats peuvent s’appliquer audelà des réseaux duplex.
24
2 2
(un) (b)
1,5 1,5
d
1 1
1
1
0,5 3 0,5
3
0 2 0
0 0,2 0,4 0,6 0,8 1 0 0,2 0,4 0,6 0,8 1
c c
2 2
(c) (d)
1,5 1,5
1
1
3
d
d
1 3 1 1
0,5 0,5
0 2 0
0 1 2 3 0 1 2 3
c c
FIGUE. 11. Le cas des matrices supralaplaciennes non commutatives. Régions synchronisées en réseau à partir de
équations principales de stabilité (à gauche) et zones de synchronisation numérique (à droite) pour les réseaux de Rössler avec
topologies intracouches de type étoile identiques et poids de couplage unàun non identiques entre les couches. Ici
je
la matrice supralaplacienne intercouche L = [1 1 ; −1 1] diag{2, 1, 1, 1, 1}, H = I11 et Γ = I11 pour
(a)(b), et H = I11 et Γ = I22 pour (c)(d).
VI. DISCUSSION
En résumé, nous développons un cadre de fonction maître de stabilité qui capture un élément essentiel
caractéristique des réseaux multiplex, que les fonctions de couplage intracouche et intercouche peuvent
être distinct. Nous définissons ici une matrice supralaplacienne distincte pour les connexions intracouches,
L
noté L , et un pour les connexions intercouches, noté L
je
L. _ Si L et moi je
faire la navette, le
Le réseau multiplex peut être facilement découplé et ainsi les modes caractéristiques du réseau intra
Les couches du Laplacien sont séparées de celles du Laplacien intercouches. (Notez cette commutation
La condition est une condition suffisante mais non nécessaire pour notre analyse théorique. Voir Sec.
IV pour plus de détails.) Nous pouvons ensuite développer une équation principale de stabilité multiplexe, Eq. (4), à
établir la région nécessaire pour une synchronisation complète. Dans la limite de l'absence d'intercalaire
le couplage du multiplex MSF se réduit à une équation maîtresse de stabilité pour chaque
25
2 2
(un) (b)
1,5 1,5
d
1 1
1
1
0,5 3 0,5 3
0,1 2
0 0
0,04 0,2 0,4 0,6 0,8 1 0,04 0,2 0,4 0,6 0,8 1
c c
2 2
(c) (d)
1,5 1,5
1 3
d
d
3 1 1
1
0,5 0,5
0 2 0
0,04 0,5 1 1,5 0 0,5 1 1,5
c c
FIGUE. 12. Le cas des matrices supralaplaciennes non commutatives. Régions synchronisées en réseau à partir de
équations principales de stabilité (à gauche) et zones de synchronisation numérique (à droite) pour les réseaux de Rössler avec
topologies intracouches entièrement connectées identiques et poids de couplage unàun non identiques entre les couches.
je
Ici la matrice supralaplacienne intercouche L = [1 1 ; −1 1] diag{2, 1, 1, 1, 1}, H = I11 et Γ = I11 pour
(a)(b), et H = I11 et Γ = I22 pour (c)(d).
couche, Eq. (5), nous permettant de calculer la région nécessaire à la synchronisation intracouche.
Dans la limite de l'absence de couplage intracouche, le multiplex MSF se réduit à une stabilité maître
équation pour chaque réseau intercouche indépendant, Eq. (6), nous permettant de calculer le
région nécessaire à la synchronisation intercouches.
Pour utiliser explicitement le framework MSF multiplex, il faut spécifier f(∙) (c'estàdire le
dynamique nodale) et les fonctions de couplage inter et intracouches (c'estàdire respectivement H et Γ).
Nous considérons spécifiquement un réseau à deux couches d'oscillateurs de Rössler et diverses formes
de H et Γ. Nous constatons que les différents types de comportements cohérents observés dans le réseau
sont déterminés par les intersections des trois régions nécessaires décrivant une synchronisation complète.
chronisation, synchronisation intracouche et synchronisation intercouche. Compte tenu d'un
topologie du réseau, ces régions peuvent ensuite être paramétrées par les couches intra et intercouches.
forces de couplage (c'estàdire c et d respectivement). La synchronisation complète est stable lorsque
26
2 2
(un) (b)
1,5 1,5
d
1
d
1
1 3 1 3
0,5 0,5
0,2 0,2
2 2
0 0
0,04 0,2 0,4 0,6 0,8 1 0,04 0,2 0,4 0,6 0,8 1
c c
2 2
(c) (d)
1,5 1,5
3 1 3
1
d
d
1 1
0,5 0,5
0,2 0,2
2 2
0 0
0,04 0,5 0,92 1,5 0,04 0,5 0,92 1,5
c c
FIGUE. 13. Dans le cas d'un réseau à trois couches entièrement connecté, la liaison intercouches est un lien (couche I — couche
II—couche III) Régions synchronisées en réseau à partir des équations principales de stabilité (à gauche) et de la synchronisation numérique.
zones de nisation (à droite) pour les réseaux de Rössler avec différentes combinaisons de H et Γ. (a)(c) H = I11 et
Γ = I11, (b)(d) H = I11 et Γ = I22.
c et d se situent tous deux dans le chevauchement des trois régions. La synchronisation intracouche est stable
lorsque c et d tombent tous deux dans le chevauchement de la région synchronisée commune et de la couche intra
région synchronisée. La synchronisation intercouches est stable lorsque c et d tombent dans le
chevauchement de la région synchronisée commune et de la région synchronisée intercouches.
Pour une dynamique nodale de réseau donnée, la région synchronisée conjointe est principalement déterminée
par les deux matrices de couplage internes H et Γ. De même, la région synchronisée intracouche est
principalement déterminé par la matrice de couplage intracouche H et la synchronisation intercouche
région par la matrice de couplage intercouches Γ. Par conséquent, en plus de la dynamique nodale, le
La fonction de couplage interne est un facteur essentiel pour déterminer quel type de synchronisation
le réseau arrivera. Si H est sous une forme telle que la région synchronisée intracouche
est vide, la synchronisation intracouche est instable quelle que soit la taille de l'espace intracouche
la force de couplage est. De même, si Γ est sous une forme telle que la région synchronisée intercouches
27
3 3
(un) (b)
2.5 2.5
2 2
2
2
d
d
1,5 1,5
1 1
0,5 0,5
0 0
00h05 0,5 0,86 1 1.2 0 0,2 0,4 0,6 0,8 1 1.2
c c
3 3
(c) (d)
2.5 2.5
1
2 2
1
d
d
1,5 1,5
1 1
0,5 0,5
0,2
0
0 0,2 0,4 0,6 0,8 1 1.2 0 0,2 0,4 0,6 0,8 1 1.2
c c
FIGUE. 14. Dans le cas d'un réseau à trois couches entièrement connecté, la liaison intercouches est un lien (couche I — couche
II—couche III) Régions synchronisées en réseau à partir des équations principales de stabilité (à gauche) et de la synchronisation numérique.
zones de nisation (à droite) pour les réseaux de Rössler avec différentes combinaisons de H et Γ. (a)(c) H = I11 et
Γ = I13, (b)(d) H = I13 et Γ = I22.
est vide, la synchronisation intercouches est instable quelle que soit la taille de l'espace intercouches.
la force de couplage est. Dans les deux cas, la synchronisation complète ne se produira pas, quel que soit le
la force de couplage.
Ici, nous avons étudié théoriquement et numériquement des réseaux duplex spécifiques de
Oscillateurs Rössler où les deux couches ont la même structure topologique. Notre approche
peut être appliqué aux réseaux multiplex avec différents choix pour la dynamique nodale interne,
différentes fonctions de couplage inter et intracouches, et davantage de couches.
Comme ce travail introduit une approche systématique pour analyser les modèles de synchronisation
dans les réseaux multiplex, l'accent est ici mis sur le cas le plus simple des réseaux multiplex où
la matrice supraLaplacienne des connexions intracouches est commutative avec celle des
connexions intercouches. Notre cadre est également valable à condition que le réseau multiplex
a une topologie intracouche qui est identique sur chaque couche et que à la fois le Lapla
28
L
matrice cienne L et la matrice laplacienne intercouches L
je
peut être diagonalisable et avoir un réel
valeurs propres. Nous vérifions numériquement dans la section IV que les équations maîtresses de stabilité dérivées
ici peut s'appliquer à une classe plus large de réseaux multiplex avec supra
Laplaciens, mais nous ne pouvons prédire que la région de synchronisation complète et intracouche
synchronisation, ou la région de synchronisation complète et de synchronisation intercouches,
et nous ne pouvons pas prédire simultanément le chevauchement de ces trois comportements de synchronisation.
iors. Établir les conditions minimales exactes dans lesquelles notre cadre peut être appliqué
reste une question ouverte importante.
REMERCIEMENTS
Ce travail est soutenu en partie par le Programme national clé de recherche et de développement
de Chine dans le cadre de la subvention 2016YFB0800401, en partie par la National Natural Science Foundation
de Chine au titre des subventions 61573004, 11501221, 61573262, 61621003 et 61532020, en partie par le
Bureau de recherche de l'armée américaine dans le cadre de l'Initiative de recherche universitaire multidisciplinaire Award No.
W911NF1310340 et accord de coopération W911NF0920053, en partie par le Programme
pour les excellents talents du nouveau siècle à l'Université de la province du Fujian en 2016 et en partie par le
Programme de promotion pour les enseignants jeunes et d'âge moyen en recherche scientifique et technologique
de l'Université Huaqiao (ZQNYX301).
Annexe A : Exposants de Lyapunov
L'exposant de Lyapunov mesure la contraction exponentielle ou le taux d'expansion de l'intérieur.
perturbations finisimales. Pour un système dynamique à temps continu à n dimensions
x˙ = G(x) (A1)
Lyapunov Les exposants sont déterminés par l'équation linéarisée par rapport à la référence
trajectoire s(t)
U˙ = J(s(t))U
avec la condition initiale U(0), où J est la matrice jacobienne de G, et s(t) satisfait l’équation (A1).
Soit vi(0)(i = 1, 2, ∙ ∙ ∙ , n) est le vecteur orthonormé de U(0). Les exposants de Lyapunov sont
29
défini comme suit

1
σi = lim ln ||U(t)vi(0)||,
t→∞ t
et le plus grand est appelé le plus grand exposant de Lyapunov (qui est supérieur à zéro pour les systèmes
chaotiques). Il joue un rôle clé dans l'analyse de la stabilité des systèmes contrôlés. On peut ajuster le
paramètre (ici la force de couplage) de telle sorte que le plus grand exposant de Lyapunov
est inférieur à zéro, et ainsi le système est contrôlé selon la trajectoire souhaitée.
Annexe B : Découplage du système de réseau multiplex
Supposons que les matrices supralaplaciennes L L Moi et L sont des matrices symétriques et satisfont
L LL I = L IL L , alors il existe une matrice inversible P telle que
P −1L LP = diag{λ1, ∙ ∙ ∙ , λM, λM+1, ∙ ∙ ∙ , λM×N },
P −1L IP = diag{µ1, ∙ ∙ ∙ , µM, µM+1, ∙ ∙ ∙ , µM×N },
λ1 = ∙ ∙ ∙ = λM < λM+1 ≤ ∙ ∙ ∙ ≤ λM×N , µk ≥ 0 (k = 1, 2, ∙ ∙ ∙ où 0 = , M × N),
et diag{υ1, ∙ ∙ ∙ , υM} désigne une matrice diagonale dont le jième élément diagonal est υj (j = , M). 1, 2, ∙ ∙ ∙
En désignant un nouveau vecteur η = [η 1 , η 2 , M×N ] = (P Im) −1ξ, on peut transformer le , η

∙∙∙
équation variationnelle (3) en
η˙ = [IM×N Df(s) − c(diag{λ1, ∙ ∙ ∙ , λM×N } H) − d(diag{µ1, ∙ ∙ ∙ , µM×N } Γ)]η. (B1)
Cela donne en outre
η˙k = [Df(s) − cλkH − dµkΓ]ηk, k = 1, 2, ∙ ∙ ∙ , M × N. (B2)
Ici, ηk représente le mode de perturbation dans l'espace propre généralisé associé à
λk et µk. Un critère pour que la variété de synchronisation soit (asymptotiquement) stable est que tous les
exposants transversaux de Lyapunov de l'équation variationnelle (B2) soient strictement négatifs.
Clairement, ces exposants de Lyapunov dépendent de la dynamique des nœuds f(∙), du réseau intra
et les forces de couplage intercouches c et d, ainsi que les matrices de couplage H et Γ. Par conséquent,
nous pouvons obtenir les trois équations principales de stabilité : Eqs. (4), (5) et (6).
30
Annexe C : Calcul des régions synchronisées R c,d
On peut calculer trois régions synchronisées en ce qui concerne les paramètres α et β : R,
RIntra et RInter des équations. (4), (5) et (6), respectivement. De plus, lorsque le réseau
topologies sont données, on peut calculer directement les valeurs caractéristiques du supraLaplacien
matrices et paramétrer ces régions en termes de c et d, puisque α = cλ et β = dµ.
Par exemple, lorsque H = I11 et Γ = I11, les modes caractéristiques non nuls α = cλ et
β = cµ devrait se situer dans Ra,β = {(α, β)| 0,2 < α + β < 4,6}, et par conséquent la région avec
par rapport aux paramètres c et d est
Rc,d = {(c, d)| 0,2 < c + 2d, c + 0,4d < 0,92}.
Pour d'autres combinaisons de H et Γ, les régions synchronisées par rapport aux paramètres c
et d peut être obtenu de la même manière.
Annexe D : Erreurs de synchronisation et fonction de l'indicateur.
Pour mesurer l'étendue de la synchronisation intracouche, intercouche et complète, nous introduisons
déduisez les indices suivants :

N
(k) 1 (k)
ETIntra(t) = x je (k) (t), k = 1, 2, ∙ ∙ ∙ (t) − x ,M (D1)
N je = 1
où ∙ et x (k) (t) est l'état moyen de tous les nœuds du kième est un opérateur de norme,
(k)
couche au temps t. Ainsi E Intra(t) est l'erreur de synchronisation des nœuds de la kème couche au moment
t, à savoir l’erreur de synchronisation intracouche.
De même, l'erreur de synchronisation intercouches est définie comme

N M
1 (k)
EInter(t) = x
je (t) − xi(t), (D2)
MN je = 1 k=1
et l'erreur de synchronisation complète est définie comme

M N
1 (k)
E(t) = (t) − x(t), x je (D3)
NouveauMexique
k=1 je = 1
où xi(t) est l'état moyen du nœud i dans chaque couche et de ses homologues dans les autres couches,
et x(t) est celui de tous les nœuds du réseau multiplex.
31
Avec ces définitions, nous utilisons la fonction indicatrice suivante pour représenter
synchronisation, synchronisation intracouche et synchronisation intercouche :
(k)
3, EInter(t) < et E Intra(t) < pour tout t > T0, (k)
2, EInter(t) ≥ et E Intra(t) < pour tout t > T0,

Identifiant =
(k) (D4)
1, EInter(t) < et E Intra(t) ≥ pour tout t > T0,
0, autre.
Ici, T0 est une valeur de seuil temporel et est un seuil donné pour les erreurs de synchronisation.
Dans les simulations, = 1,0 × 10−2 , et T0 = 0,8Ttotal (Ttotal est le temps d'évolution total).
Il est évident que le réseau atteint une synchronisation complète lorsque Id = 3, intracouche
synchronisation lorsque Id = 2, synchronisation intercouche lorsque Id = 1, et aucun des
cidessus lorsque Id = 0.
Annexe E : Analyse théorique pour le cas de synchronisation complète
Considérons un réseau duplex composé de deux sousréseaux avec la même topologie interne
et une connectivité intercouche un à un entre les nœuds. L'évolution dynamique peut être
écrit comme :
N
x˙i = f(xi) − c lijHxj − d Γ(axe − non),

j=1
N je = 1, 2, ∙ ∙ ∙ , N. (E1)
y˙i = f(yi) − c lijHyj − d Γ(byi − bxi),
j=1
je
une −une
Ici la matrice laplacienne intercouches L =
(a et b sont des réels non négatifs
−bb
2 2
constantes satisfaisant un +b = 0), indiquant que l'échange d'informations entre les couches
est asymétrique et pondéré lorsque a = b. Notez que le réseau duplex (E1) avec a = b = 1 a
été discutée dans le texte principal. Les matrices supralaplaciennes intracouches et intercouches,
L je je
L = I2 L et L =L IN satisfait la condition commutative, où L = (lij )N×N est
la matrice laplacienne intracouche et Im est une matrice identité d'ordre m.
Ensuite, nous expliquerons théoriquement comment les modèles de synchronisation observés nécessitent le
chevauchement des différentes régions de synchronisation sous deux aspects : la synchronisation intracouche
équations de stabilité de chronisation dérivées ensuite dans la Sec. E 1 et la synchronisation intercouches
équations de stabilité dérivées de la Sec. E2.
32
1. Équations de stabilité de synchronisation intracouche
Soit sx(t) et sy(t) les états synchrones intracouches de la couche x et de la couche y,
respectivement, qui sont dominés par les équations suivantes :
s˙x = f(sx) − ad Γ(sx − sy),

(E2)
s˙y = f(sy) − bd Γ(sy − sx).
La linéarisation du réseau duplex (E1) aux états synchrones intracouches sx et sy donne
˙ N
δxi = Df(sx)δxi − à Γ(δxi − δyi) − c lijHδxj ,
˙
j=1
N
je = 1, 2, ∙ ∙ ∙ , N. (E3)
δyi = Df(sy)δyi − bd Γ(δyi − δxi) − c lijHδyj ,
j=1
Df(x) 0
Notons δzi = (δx T
je , δy T ) T , Df (sx, sy) = δz = (δz T
1, δz 2T , ∙∙∙
, δz N
T T
) ,
je
0 Df(s) ,
puis l'équation. (E3) peut être réécrit comme
˙
je
δz = IN [Df (sx, sy) − d (L Γ)]δz − c(L (I2 H))δz. (E4)
Puisque la matrice laplacienne L = (lij ) est symétrique (en supposant que les liens au sein de chaque couche sont
non orienté), il existe une matrice inversible P telle que
l1
l2
P −1LP = ,
..
.
λN
−1 T T T T
ici 0 = λ1 < λ2 ≤ λ3 ≤ ∙ ∙ ∙ ≤ λN . Soit ξ = (P I2m) δz, et ξ = (ξ , XN ) ,
∙∙∙
1 fois2 ,
nous avons
˙ξj = [Df (sx, sy) − d (L

je
Γ)]ξj − cλk(I2 H)ξj , j = 1, 2, ∙ ∙ ∙ , N.
En négligeant l'indice j, on peut obtenir la forme générale de l'équation maîtresse de stabilité
pour la synchronisation intracouche :
je
η˙ = [Df (sx, sy) − d (L Γ)]η − α(I2 H)η, (E5)
où α = cλ, et λ est toute valeur propre non nulle de la matrice laplacienne L.
33
Si le réseau duplex (E1) atteint une synchronisation complète, cela signifie que les deux réseaux intra
les états synchrones de couche sx et sy convergent vers le même état s dominé par l'état isolé
système nodal : s˙ = f(s). Alors l’équation variationnelle de (E2) en s
˙δs = [Df(s) − (a + b)dΓ]δs
une −une
devrait être stable. Puisque les valeurs propres de L
je = sont µ = 0, a + b et β =
−bb
µd où µ est la valeur propre non nulle, l'équation variationnelle cidessus peut être en conséquence
transformé sous la forme générale suivante :
˙δs = [Df(s) − βΓ]δs, (E6)
qui est stable pour β RInter

b et toute valeur de force de couplage α. Cela donne
R.Inter = {β |σ(β) < 0},

b
où σ(β) est le plus grand exposant de Lyapunov de l'équation (E6). Pour plus de commodité, nous incluons
le paramètre α dans RInterb , et obtenir
R.Inter = {(α, β)|σ(β) < 0, α ≥ 0}. (E7)

un B
Nous appelons RInter

un B la région synchronisée intercouches par rapport à α et β.
Simultanément, lorsque le réseau duplex atteint une synchronisation complète, l'équation. (E5)
est stable en s, ce qui signifie que l'équation suivante est stable à l'origine :
je
η˙ = [I2 Df(s) − d (L Γ)]η − α(I2 H)η. (E8)
une −une
Diagonaliser la matrice L
je = et faire une simple transformation linéaire,
−bb
,
on peut obtenir les équations découplées de (E4) :
˙ζ1 = [Df(s) − αH]ζ1, (E9)
et
˙ζ2 = [Df(s) − αH − (a + b)dΓ]ζ2.
Semblable à l’argument cidessus, remplacez (a + b)d par β, et le deuxième des éléments découplés
les équations se transforment en
˙ζ2 = [Df(s) − αH − βΓ]ζ2. (E10)
34
L’équation (E9) est stable lorsque α RIntra un B

{α |σ(α) < 0, β ≥ 0}, ici σ(α) est le plus grand
Exposant de Lyapunov de l'équation. (E9) avec paramètre α. De même, l’équation (E10) est stable lorsque
(α, β) Rα,β {(α, β)|σ(α, β) < 0}, où σ(α, β) est le plus grand exposant de Lyapunov de
Éq. (E10). Pour plus de commodité, nous appelons RIntra

un B , RInter et Rα,β l'intracouche, intercouche
un B
et des régions synchronisées conjointes, respectivement. Étant donné une topologie de réseau intracouche spécifiée,
Entrer
un B , RInter et Rα,β peuvent être paramétrés par c (la force de couplage intracouche) et d
α,β
(la force de couplage intercouches), notée RIntra c,d , RInter et Rc,d, respectivement.
CD
En résumé, pour atteindre une synchronisation complète dans le réseau duplex (E1), il faut
que trois équations de stabilité de synchronisation, (E6), (E9) et (E10), sont simultanément
écurie. Ainsi, les modes caractéristiques intracouches α = λc et la caractéristique intercouches
les modes β = µd doivent tomber dans le chevauchement de RIntra un B , RInter et Rα,β, soit RIntra ∩
un B un B
RInter ∩ Rα,β. Cela indique que les forces de couplage intracouche et intercouche doivent diminuer
α,β
dans le chevauchement de RIntra

c,d , RInter et Rc,d, c'estàdire RIntra ∩ RInter ∩ Rc,d, lorsque la couche intracouche
c, cc, d CD
je
topologie du réseau L et mode de liaison intercouches L sont précisées.
2. Équations de stabilité de synchronisation intercouches
À la sec. E 1 nous sommes partis des équations de stabilité de synchronisation intracouche, nous pouvons
analyser également ce problème à partir de l’approche de synchronisation intercouches. Notons si(t)(i =
1, 2, ∙ ∙ ∙ , N) comme états synchrones intercouches, qui sont dominés par les éléments suivants
équations :
N
s˙i = f(si) − c lijHsj , je = 1, 2, ∙ ∙ ∙ ,N (E11)

j=1
Linéariser le réseau duplex (E1) aux états synchrones intercouches si , alors on obtient
˙ N
δxi = Df(si)δxi − c lijHδxj − ad Γ(δxi − δyi),
j=1
˙ N je = 1, 2, ∙ ∙ ∙ , N. (E12)
δyi = Df(si)δyi − c lijHδyj − bd Γ(δyi − δxi),
j=1
Notant δzi = δxi − δyi , on obtient de l’équation (E12) que
˙ N
δzi = Df(si)δzi − c lijHδzj − (a + b)d Γδzi , je = 1, 2, ∙ ∙ ∙ ,N

j=1
35
T
Soit δZ = (δz 1,
δz 2T , ∙∙∙ T T
, δz N ) , nous obtenons ainsi l'équation maîtresse de stabilité suivante pour
synchronisation intercouches :
˙
δZ = DF(s1, s2, ∙ ∙ ∙ , sN )δZ − c(L H)δZ − (a + b)d(IN Γ)δZ, (E13)
où
Df(s1)
Df(s2)
DF(s1, s2, ∙ ∙ ∙ , sN ) = .
..
.
Df(sN )
Il convient de noter que l'équation (E11) dominant l'état synchrone intercouches si(t) peut
être linéarisé en s(t) comme
N
˙
δsi = Df(s)δsi − c lijHδsj , je = 1, 2, ∙ ∙ ∙ , N (E14)
j=1
et peut donc être réécrit comme
˙
δS = [IN Df(s) − c(L H)]δS, (E15)
T T T T
où δS = (δs , δs , ∙∙∙
, δs ) . Diagonalisons la matrice laplacienne L = (lij ), et on peut
obtenir les équations découplées
˙ξk = [Df(s) − cλkH]ξk, k = 2, 3, ∙ ∙ ∙ , N, (E16)
où 0 = λ1 < λ2 ≤ λ3 ≤ λN sont les valeurs propres de L, et la forme générale de (E16) est
η˙ = [Df(s) − αH]η. (E17)
Maintenant, si le réseau duplex (E1) réalise une synchronisation complète, cela signifie que
si(i = 1, 2, ∙ ∙ ∙ , N) converge vers un état synchrone s. Il faut donc exiger que
Éq. (E17) est stable. Évidemment, l’équation. (E17) est stable à l’origine lorsque α RIntra un B {une |s(une) <
0, β ≥ 0} (σ(α) est le plus grand exposant de Lyapunov de l'équation (E17) ). Compte tenu de l'intracouche
structure et liaison intercouches, RIntra un B peut être paramétré par les forces de couplage
c et d, notés RIntra c,d .
Simultanément, remplaçons si(i = 1, 2, ∙ ∙ ∙ , N) dans l'équation. (E13), il y a
˙
δZ = [IN Df(s) − c(L H)]δZ − (a + b)d(IN Γ)δZ. (E18)
36
En diagonalisant la matrice L et en effectuant une simple transformation linéaire, on peut obtenir le
équations découplées suivantes :
˙ζ1 = [Df(s) − (a + b)dΓ]ζ1,
˙ζ2 = [Df(s) − αH − (a + b)dΓ]ζ2.
Semblable à la méthode de traitement dans (I), remplacez (a + b)d par β, les équations découplées cidessus
peut être transformé en
˙ζ1 = [Df(s) − βG]ζ1, (E19)
˙ζ2 = [Df(s) − αH − βΓ]ζ2. (E20)
L’équation (E19) est stable lorsque β RInter un B

{β |σ(β) < 0, α ≥ 0}, ici σ(β) est le plus grand
Exposant de Lyapunov de l'équation. (E19). De même, l’équation. (E20) est stable lorsque (α, β) Rα,β
{(α, β)|σ(α, β) < 0}, et σ(α, β) est le plus grand exposant de Lyapunov de l'équation. (E20). Pour
commodité, nous appelons RIntra RInter et Rα,β les intracouches, intercouches et articulaires synchronisés
une,b , une,b
régions, respectivement. Étant donné une topologie de réseau intracouche spécifiée, RIntraα,β , α,βRInter et Rα,β
peuvent être paramétrés par c (la force de couplage intracouche) et d (la force de couplage intercouche
force), noté RIntra c,d , c,d

RInter et Rc,d, respectivement.
En résumé, pour obtenir une synchronisation complète en réseau duplex (E1), il faut
que les équations. (E17), (E19) et (E20) sont simultanément stables. Cela fonctionne lorsque (α, β)
Entrer
un B
∩RInter
α,β ∩Rα,β, ou lorsque (c, d) RIntraCD ∩RInter
c,d ∩Rc,d après la topologie intracouche et
les moyens de liaison intercouches sont définitivement donnés. Il vérifie à nouveau que la synchronisation complète
se produit dans le chevauchement de RIntra RInter et Rα,β.

une,b , une,b
Annexe F : Un exemple concret, le système CORS
Le réseau de stations de référence en fonctionnement continu (CORS) [39] comprend
de stations de référence du système mondial de navigation par satellite (GNSS) fonctionnant en permanence, un
réseau de communication et centres de données. Grâce à l'observation continue par satellite GNSS
Lites et traitement des mesures GNSS, le système CORS est largement appliqué à divers domaines
tels que : positionnement tridimensionnel, navigation et synchronisation à différents niveaux de précision,
suivi et détermination de l'orbite des satellites, maintien du cadre de référence de la terre,
recherche géodynamique telle que les tremblements de terre et les mouvements des plaques, ainsi que le niveau de la mer, l'ionosphère
37
et surveillance de la vapeur d'eau. Le système CORS est un système essentiel d’information géospatiale.
infrastructure avec de nombreux pays et régions du monde, tels que la Chine, l'Amérique et l'Europe
et l'Australie ayant mis en place des systèmes CORS.
FIGUE. 15. Un schéma des systèmes CORS du système chinois de navigation par satellite BeiDou. La couche cidessous est
la couche physique composée de certains récepteurs GNSS, et la couche supérieure est le traitement des données correspondant
couche.
Le système CORS peut être considéré comme un réseau multicouche composé d'une couche physique
et une couche de données, comme le montre la figure 15. La couche physique se compose de N récepteurs, qui suivent le
mêmes m satellites, pour recevoir les données de positionnement. Parce qu'il existe des horloges biaisées entre
récepteurs, ces récepteurs doivent ajuster leurs horloges et obtenir une synchronisation temporelle pour
améliorer la précision du positionnement [41, 42]. C'estàdire la synchronisation intracouche avec
le respect du temps est requis dans la couche physique. Toutefois, le traitement des données n'est pas
effectué par les récepteurs, mais plutôt par les unités de traitement de données dans les centres de données [40]. Ces
unités de traitement de données relient respectivement leurs récepteurs pour former une couche de données virtuelle
(couche supérieure sur la figure 15) et un réseau de hiérarchie numérique synchrone (SDH). Chaque unité de la couche de données
doit réaliser la synchronisation des données pour un positionnement précis, ce qui implique que la couche de données
doit réaliser une synchronisation intracouche.
[1] DJ Watts et SH Strogatz, Nature 393, 440 (1998).
[2] M. Barahona et LM Pecora, Phys. Le révérend Lett. 89, 054101 (2002).
[3] H. Hong, MY Choi et BJ Kim, Phys. Rév.E 65, 026139 (2002).
38
[4] MEJ Newman, SIAM Rév. 45, 167 (2003).
[5] J. Lu, X. Yu, G. Chen et D. Cheng, IEEE Trans. Système de circuits. I 51, 787 (2004).
[6] J. Lu et G. Chen, IEEE Trans. Automate. Contr. 50, 841 (2005).
[7] J. Zhou, JA Lu et J. Lü, IEEE Trans. Automate. Contr. 51, 652 (2006).
[8] CW Wu, Synchronisation dans un réseau complexe de systèmes dynamiques non linéaires (World Sci
entic, Singapour, 2007) pp. 51123.
[9] S. Boccalettia, V. Latorab, Y. Morenod, M. Chavezf et D.U. Hwang, Phys. Rep. 424, 175
(2006).
[10] A. Arenas, A. D´ıazGuilera, J. Kurths, Y. Morenob et C. Zhou, Phys. Rep. 469, 93 (2008).
[11] Y. Chen, J. Lu, F. Han et X. Yu, Syst. Contr. Lett. 60, 517 (2011).
[12] L. Huang, YC Lai et RA Gatenby, Chaos 18, 013101 (2008).
[13] L. Donetti, PI Hurtado et MA Muñoz, Phys. Le révérend Lett. 95, 188701 (2005).
[14] LM Pecora, F. Sorrentino, AM Hagerstrom, TE Murphy et R. Roy, Nat. Commun.
5, 4079 (2014).
[15] L. Tang, JA Lu et G. Chen, Chaos 22, 023121 (2012).
[16] CD Brummitt, RM D'Souza et EA Leicht, Proc. Natl Acad. Sci. ÉtatsUnis 109, E680
(2012).
[17] F. Radicchi et A. Arenas, Nat. Phys. 9, 717 (2013).
[18] S. Goomez, A. DíazGuilera, JGG. nes, CJ PérezVicente, Y. Moreno et A. Arenas,
Phys. Le révérend Lett. 110, 028701 (2013).
[19] MD Domenico, A. Sol'eRibalta, S. Gomez et A. Arenas, Proc. Natl. Acad. Sci. CERF
111, 8351 (2014).
[20] G. D'Agostino et A. Scala, Réseaux de réseaux : la dernière frontière de la complexité
(Springer, Berlin, 2014).
[21] T. VallesCatala, FA Massucci, R. Guimera et M. SalesPardo, Phys. tour. X6, 011036
(2016).
[22] M. Kivel a, A. Sands, M. Barthelemy, JP Gleeson, Y. Moreno et MA Porter, J. Complex
Réseau. 2, 203 (2014).
[23] S. Boccaletti, G. Bianconi, R. Criado, CD Genio, J. GómezGardenes, M. Romance,
I. SendiñaNadal, Z. Wang et M. Zanin, Phys. Rep.544, 1 (2014).
[24] O. Yagan, D. Qian, J. Zhang et D. Cochran, IEEE Journal on Selected Areas in Commu
39
communications 31, 1038 (2013).
[25] LM Pecora et TL Carroll, Phys. Le révérend Lett. 80, 2109 (1998).
[26] L. Tang, JA Lu, J. Lu et X. Yu, Int. J. Bifurcat. Chaos 22, 1250282 (2012).
[27] L. Tang, JA Lu, J. Lu et X. Wu, Int. J. Bifurcat. Chaos 24, 1450011 (2014).
[28] L. Tang, X. Wu, J. Lu et JA Lu, Chaos 25, 033101 (2012).
[29] A. SoléRibalta, MD Domenico, NE Kouvaris, A. D´ıazGuilera, S. Gómez et A. Arenas,
Phys. Tour. E88, 032807 (2013).
[30] J. Aguirre, R. SevillaEscoboza, R. Gutiérrez, D. Papo et JM Buld'u, Phys. Rev. Lett.
112, 248701 (2014).
[31] M. Xu, J. Zhou, JA Lu et X. Wu, Eur. Phys. JB 88, 240 (2015).
[32] Y. Li, X. Wu, JA Lu et J. Lü, IEEE Trans. Circuits Sys. II : Mémoires express 63, 206
(2015).
[33] F. Sorrentino, New J. Phys. 14, 033035 (2012).
[34] D. Irving et F. Sorrentino, Phys. Rév.E 86, 056102 (2012).
[35] CID Genio, J. GómezGardeñes, I. Bonamassa et S. Boccaletti, Science Advances 2,
Non. 11 (2016).
[36] LV Gambuzza, M. Frasca et JGG. nda, EPL 110, 20010 (2015).
[37] R. SevillaEscoboza, IS na Nadal, I. Leyva, R. Guti'errez, JM Buld'u et S. Boccaletti,
Chaos 26, 065304 (2016).
[38] GA Leonov et NV Kuznetsov, Int. J. Bifurcat. Chaos 17, 1079 (2007).
[39] C. Rizos, GPS Solut. 11, 151 (2007).
[40] G. Blewitt, dans GPS for Geodesy, édité par P. Teunissen et A. Kleusberg (Springer, Berlin
Heidelberg, 1998) pp.
[41] H. Liu, RF Zhang, JN Liu et M. Zhang, Sci. Technologie chinoise. Sci. 59, 9 (2016).
[42] D. Kim et RB Langley, NavigationAlexandria 49, 205 (2002).
40

Master Stability Function

Transféré par

Droits d'auteur :

Formats disponibles

Master Stability Function

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Master Stability Function

Transféré par

Droits d'auteur :

Formats disponibles

Machine Translated by Google

Fonctions de stabilité maîtrisées pour une synchronisation complète, intracouche et

Fonctions de stabilité maîtrisées pour complet, intra­couche et inter­couche

1Laboratoire clé de sciences informatiques de l'Université de la province du Fujian,

École des sciences mathématiques, Université Huaqiao, Quanzhou 362021, Chine.

2École de Mathématiques et Statistiques,

Université de Wuhan, Wuhan 430072, Chine.

3Département d'informatique, Université de Californie, Davis CA 95616, États­Unis.

Académie chinoise des sciences, Pékin 100190, Chine.

respectivement les régions nécessaires de synchronisation complète, de synchronisation intra­couche et

l'apparition de synchronisation intra­couche et inter­couche dépend principalement des fonctions de couplage

nisation; synchronisation inter­couches ; région synchronisée.

* À qui la correspondance doit être adressée : [email protected]

La synchronisation dans un réseau d'éléments connectés est essentielle au bon fonctionnement

une attention croissante ces dernières années.

établi via l'approche MSF, si un réseau peut ou non réaliser la synchronisation

et discuté de la synchronisabilité via le rapport propre de la matrice laplacienne. Aguirre et coll.

deux couches pour trois couplages spécifiques

modèles, et déterminé qu'il existe une valeur optimale du couplage inter­couches

scénarios, toujours en utilisant une approche MSF [35], et montrer comment

couplés indirectement via un milieu inhomogène. Le système est formalisé en termes

synchronisation, synchronisation intra­couche et synchronisation inter­couches en multiplex

Rive l'équation maîtresse de stabilité pour un réseau multiplex où le supra­laplacien de

supra­Laplaciens non commutatifs.

II. UN CADRE DE FONCTION DE STABILITÉ MAÎTRE POUR LES CLASSES DE

A. Un modèle de réseau multiplex

Nous considérons un réseau multiplex constitué de M couches constituées chacune de N nœuds.

je je1 , je2 , dans le ),

au sein de chaque couche, respectivement, et Γ(∙) : R m → R m et d sont la fonction de couplage interne

Par souci de simplicité, notons

31, 32] et les références qui y figurent.

B. Trois équations maîtresses de stabilité

fonctions de couplage interne.

où ξ = x − 1M 1N s et IM×N est la matrice identité d'ordre M × N.

multiplex pour un système décrit par l'équation. (1):

y˙ = [Df(s) − αH − βΓ]y, (4)

A pour plus d'informations sur les exposants de Lyaponuv.

choisi dans [0, +∞), et Eq. (4) se réduit à

y = [Df(s) − αH]y, (5)

indépendant (pas de couplages inter­couches).

De même, lorsque λ = 0 et µ = 0, on peut obtenir l'équation suivante

y˙ = [Df(s) − βΓ]y, (6)

condition de stabilité du collecteur de synchronisation dans un réseau multiplex.

Il convient de noter en particulier le cas où les fonctions de couplage intra­ et inter­couches

Il s'agit de l'équation principale de stabilité pour le réseau composite unique correspondant où

Le laplacien peut décrire sa topologie. C’est­à­dire que l’équation maîtresse de stabilité du

le réseau composite unique est un cas particulier de l'équation. (4).

synchronisation complète et synchronisation intra­couche. Si la topologie sur chaque couche est

synchronisation et synchronisation inter­couches. (Voir la partie simulation pour plus de détails.)

types de comportements de synchronisation : synchronisation complète, synchronisation intra­couche

Rα,β = {(α, β)|σ(α, β) < 0},

Chaque fois que σ(α, β) < 0, perturbations transversales à la variété de synchronisation

s'éteignent et le réseau est dit synchronisable.

espace, nous incluons explicitement le paramètre β dans la définition de la région,

l’équation. (6), nous obtenons la région de synchronisation inter­couches

régions nisées par rapport à α et β.

Le célèbre oscillateur chaotique de Rössler est décrit comme

où a = b = 0,2 et c = 9. C'est la fonction f dans le réseau multiplex Eq. (1). Que

Fonctions de stabilité maîtrisées pour complet, intracouche et intercouche

3Département d'informatique, Université de Californie, Davis CA 95616, ÉtatsUnis.

respectivement les régions nécessaires de synchronisation complète, de synchronisation intracouche et

l'apparition de synchronisation intracouche et intercouche dépend principalement des fonctions de couplage

nisation; synchronisation intercouches ; région synchronisée.

modèles, et déterminé qu'il existe une valeur optimale du couplage intercouches

synchronisation, synchronisation intracouche et synchronisation intercouches en multiplex

Rive l'équation maîtresse de stabilité pour un réseau multiplex où le supralaplacien de

supraLaplaciens non commutatifs.

indépendant (pas de couplages intercouches).

Il convient de noter en particulier le cas où les fonctions de couplage intra et intercouches

Le laplacien peut décrire sa topologie. C’estàdire que l’équation maîtresse de stabilité du

synchronisation complète et synchronisation intracouche. Si la topologie sur chaque couche est

synchronisation et synchronisation intercouches. (Voir la partie simulation pour plus de détails.)

types de comportements de synchronisation : synchronisation complète, synchronisation intracouche

l’équation. (6), nous obtenons la région de synchronisation intercouches

Ensuite, nous donnons d'abord un aperçu de la synchronisation intracouche et intercouche, puis

calculer les régions paramétrées de synchronisation pour l'intracouche générale (inconnue)

topologies de couches dans la dernière soussection.

A. Synchronisation intracouche et intercouche

synchronisation et synchronisation intercouches. Comme le montre la figure 1, la synchronisation intracouche

B. Régions synchronisées pour les topologies intracouches inconnues

FIGUE. 1. Représentation schématique de (a) la synchronisation intracouche et (b) la synchronisation intercouche

qui est audessus de la ligne pointillée rouge.

et RInterα, β . Par conséquent, accordez

synchronisation intracouche ou intercouche.

C. Régions synchronisées avec des topologies intracouches données

iors : synchronisation complète, synchronisation intracouche et synchronisation intercouche.

le réseau R¨ossler intracouche/intercouche indépendant par rapport à α/β ; (b) le synchronisé

la région blanche représente la nonsynchronisation.

synchronisation intracouche, Id = 1 pour la synchronisation intercouche et Id = 0 pour aucun des

cas cidessus. Voir l'annexe D pour plus de détails.

le réseau R¨ossler intracouche/intercouche indépendant par rapport à α/β ; (b) le synchronisé

la région blanche représente la nonsynchronisation.

force de couplage intra et intercouches c et d. Le panneau (d) montre le nombre numérique

de Rc,d et RIntra CD détermine la synchronisation intracouche et l'intersection de Rc,d et

(d). De plus, la simple synchronisation intracouche (sans synchronisation intercouche)

zone de synchronisation intercouches (sans synchronisation intracouche) en cyan (gris clair)

le réseau R¨ossler intracouche/intercouche indépendant par rapport à α/β ; (b) le synchronisé

la région blanche représente la nonsynchronisation.

(d) des Fig. 46.

région synchronisée soustrayant la partie synchronisée intercouches, c'estàdire Rc,d∩RIntraCD −RInter

le réseau R¨ossler intracouche/intercouche indépendant par rapport à α/β ; (b) le synchronisé

la région blanche représente la nonsynchronisation.

la matrice de couplage intercouches Γ rend la région synchronisée intercouches RInter CD

synchronisation complète, comme le montre la figure 4. Si la matrice de couplage intracouche H fait

la région synchronisée intracouche RIntra CD

entre couches, la matrice supralapalacienne intercouches étant L je = L je diag{2, 1, 1, 1, 1},

la différence des topologies intracouches peut conduire à la modification de l'intercouche réelle

le changement des régions synchronisées intracouches réelles.

En d’autres termes, même si ici les matrices supralaplaciennes inter et intracouches ne

chronisation et synchronisation intracouche, ou pour une synchronisation complète et intercouche

synchronisation complète et synchronisation intracouche. Dans quelle mesure l'observation

région nécessaire à la synchronisation intercouches.

chronisation, synchronisation intracouche et synchronisation intercouche. Compte tenu d'un

forces de couplage (c'estàdire c et d respectivement). La synchronisation complète est stable lorsque

région synchronisée. La synchronisation intercouches est stable lorsque c et d tombent dans le