周期性边界条件下最近邻Kuramoto模型的“迟滞”现象

周期性边界条件下最近邻Kuramoto模型的“迟滞”现象 周期性边界条件下最近邻 Kuramoto 模型的 "迟滞"现象 王念川1，肖井华2 (1. 北京邮电大学理学院，北京 100876; 5 2. 北京邮电大学，北京 100876) 摘要:Kuramoto 模型自从 20 世纪中页提出以后，由于其包含简单而丰富的动力学现象，人 们对该模型的研究从未间断。本论文主要研究在最近邻带周期性边界条件的 Kuromoto 模型 下的同步问题。研究发现，由于振子初始相位的随机性，导致结果中振子频率的多态现象。 又利用序参量时序图，来验证结果的正确性。通过剥离态算法，本文得到分叉图中“迟滞现 10 象”。 关键词:非线性;同步;迟滞现象 中图分类号:[O415.6] 15 Phenomenon of arrearage in the nearest neighbor Kuramoto model with periodic boundary condition WANG Nianchuan1, XIAO Jinghua2 (1. Beijing University of Posts and Telecommunication, School of science, Beijing 100876; 2. Beijing University of Posts and Telecommunication, Beijing 100876) Abstract: In this paper,we introduce about synchronization and the progress about the research of 20 it .Here we mainly study about the problem of synchronization in nearest neighbor Kuramoto model with periodic boundary condition.From our result ,there are multistable behavior about the frequency of the oscillators as the difference of initial condition.Then we use the order parameter to verify our result.finally,through stripped algorithm,we can get the arrearage phenomenon in the graph. 25 Keywords: nonlinear; synchronization; arrearage phenomenon 0 引言 同步是自然界的一种基本现象[1]。它通常是指:至少在两个振动系统相位之间的一种协 调一致现象。关于同步的最早研究可以追溯到 1673 年惠更斯(C.Huygens)关于耦合单摆的 30 同步现象的观察。实际上，若干个耦合单元之间通过相互作用达到同步的现象在许多领域中 屡见不鲜。例如，对人体的心肺功能的研究 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 明，人的心跳和呼吸频率是在若干个有理的比 例的同步(锁相)之间变化，尽管他们的自然频率非常不同;在生态学上的一个典型的例子 就是马来西亚岛屿上成千上万只萤火虫在夜晚同步地闪动荧光，每个只虫子(个体)尽管各 有差异，但他们之间通过某种方式耦合确可以达到同步状态;在化学方面，人们对化学波等 35 的研究也发现许多与同步有关的现象;在光学中，人们发现若得到高功率的激光相干输入， 通过耦合使各激光器达到同步状态为最佳。上述几个简单例子说明同步现象不仅在自然界广 泛存在，而且在实际应用中是非常重要的。 1 Kuramoto 模型 在 20 世纪 70、80 年代，Kuramoto 一直致力于集体同步问题的研究，并提出了一个可 40 以很好地解释集体同步行为的模型[2]。由具有任意自由频率的相振子通过相位差的正弦值进 作者简介:王念川，男，硕士研究生，主要研究方向:非线性动力学。E-mail: w394832880@126.com 通信联系人:肖井华，(1965-)，男，教授，混沌控制及其在通信中的应用。E-mail: jhxiao@bupt.edu.cn -1- 行耦合而组成的系统，被称为 Kuramoto 模型。振子在无相互作用时以其自身频率,，i, 振动， 有全局耦合的相互作用时系统的运动方程为: N d,i 式(1) , N , , ，i ， ,i ), i , 1, 2, sin(, j dt N K j ,1 K 代表耦合强度， N 为系统尺寸，,i 为振子 i 的相位，是一个绕着极限环运动 其中，45 并随时间增长的变量， ,i 表征振子的瞬时频率，且振子间的相互作用函数为周期函数 sin(, ) 。振子间由于相互作用它们离开各自的固有频率。如果耦合强度足够大，所有振子 将达到同步，它们的平均频率相等，也就是说，此时振子的相位存在着固定的关系，出现锁 相(锁频)的现象。存在一个临界的耦合强度 Kc ，当 K , K c 时，所有振子没有达到相同 50 步，当 K , K c 时，所有振子达到锁相的状态[3]。 我们可以很容易理解，如果系统振子本身的频率是固定的，那么振子在不同的排列方式 下，达到同步时所需的临界耦合强度肯定是不同的。本文所讨论的模型为有周期边界条件的 最近邻耦合 Kuramoto 振子系统，方程如下: N K 式(2) ,i , ，i ， ,i ) ， sin(,i 1 ,i )]. [sin(, i，1 3 j,1 随着耦合强度的增加，系统会逐渐达到同步。在由非同步到同步的过程中，存在一个临 55 界的耦合强度 Kc ，当 K , K c 时，部分振子会先达到同步;当 K , K c 时，所有振子的频率 都相等，即锁频。在实验模拟中，我们认为振子的平均频率相等时，振子就达到了同步。 1.1 不同参数结构下系统的临界耦合强度值 我们令振子的自身频率固定，当变换振子的位置时，系统的参数结构与所有振子达到锁 60 相时的临界耦合强度之间存在一定的关系;系统尺寸为 N ，每个振子 i(i , 1, 2, , N ) 的自 , N ) 。所有振子在一个环上排列，按照式(2)，最近邻的振子间有 身频率为 ，i (i , 1, 2, 直接的强度为 K 的耦合相互作用。不同的振子排列方式，所有振子达到锁相时所需的临界 耦合强度是不尽相同的[4]。由 数学 数学高考答题卡模板高考数学答题卡模板三年级数学混合运算测试卷数学作业设计案例新人教版八年级上数学教学计划 知识，我们可以算出 N 个振子的所有不相同的排列数目 N 5 为 PN / 2 N 。以 5 个振子为例，不相同的排列数目为: (1 / 2)P5 / 5 , 12 。 65 图 1 振子的排列方式 (a)振子的排列方式为 12345，(b) 振子的排列方式 23451,(c)振子的排列方式为 15432。 从环上最上面的振子顺时针编号，振子的这三种排列方式对应的结构是相同的，(b)经过顺时针旋转(c)经过 对称变换都可以得到(a) 70 -2- 图 2. N , 5个 Kuramoto 振子组成一维耦合环系统在 12 种不同的振子排列方式下的分岔树及临界耦合强度 值 K c ，初始自身频率相对固定。 75 2 随机初始条件下的相同步 当我们将每个振子的初始相位随机时，对于排列 14253，,得到的分岔图如下: 图 3.各个振子的频率随耦合强度的增加的分岔图。横轴为耦合强度 K，纵轴为每个振子的平均频率 F。 可以看到，在 K>1.5 之后，还有两处地方出现了不同步的现象。而以前一直认为，当所 80 有振子的平均频率都相等后，系统处于锁相状态，也就是随着耦合强度的增加，所有振子的 平均频率会一直相等。 同样，我们遍历 6 个振子，发现在排列方式为 135264 时，得到的分岔图中也有类似的 现象。 -3- 85 图 4.随机初始条件下 6 个振子的分岔图 也就是说在同步之后，这里也出现了不同步的现象。而且，在随机初始相位的情况下。 振子的分岔图也不再像固定初始相位时那么光滑。 90 序参量 R[5]的表达式为: N 1 i , i , l , re e N l , 1 式(3) 当系统频率同步的时候，则每个振子的频率相等，则在相同时内每个振子相位改变相同， 所以序参量是不变的，对与上面的等式就只是左边的相位, 发生变化。序参量的大小表征 相位同步的程度。当所有振子的相位相等时，r=1。而实际上，如果所有振子的初始相位不 95 相同的话，无论 K 为多大，经过时间 t 后，相位不可能都相等，r 都不可能等于 1。所以序 参量只能是随着 K 的增加无限趋近于 1.当 K>Kc 时，序参量时不变。下面的图，K