• 検索結果がありません。

SynchronizationPhenomenainTwoRingsofCoupledvanderPolOscillators 2 つのリング状結合ファンデアポール発振器の同期現象

N/A
N/A
Protected

Academic year: 2021

シェア "SynchronizationPhenomenainTwoRingsofCoupledvanderPolOscillators 2 つのリング状結合ファンデアポール発振器の同期現象"

Copied!
4
0
0

読み込み中.... (全文を見る)

全文

(1)

社団法人 電子情報通信学会

THE INSTITUTE OF ELECTRONICS,

INFORMATION AND COMMUNICATION ENGINEERS

信学技報

TECHNICAL REPORT OF IEICE.

2 つのリング状結合ファンデアポール発振器の同期現象

成相 大貴 TranMinh Hai 上手 洋子 西尾 芳文

徳島大学工学部 〒 770–8506 徳島県徳島市南常三島 2–1 E-mail: {nariai,minhhai,uwate,nishio}@ee.tokushima-u.ac.jp

あらまし

本研究では、2 種類のリング状結合発振器を抵抗で結合した新たな結合発振器システムを提案する。1 つ目 の回路は 3 個のファンデアポール発振器を抵抗でリング状に結合させた。2 つ目の回路は 3 個の van der Pol 発振器の 構造を変えた発振器をインダクタを介して抵抗で落とすタイプでリング状に結合させた。私たちはこれら2つのリン グ状結合発振器を抵抗で結合させ、その抵抗の結合強度に着目して研究を行った。結合強度を変えることで同期現象 がどのように変化するのかを調査した。

キーワード

van der Pol 発振器、同期現象

Synchronization Phenomena

in Two Rings of Coupled van der Pol Oscillators

Daiki NARIAI , Tran MINH HAI , Yoko UWATE , and Yoshifumi NISHIO

Electrical and Electronic Engineering, Tokushima University, 2-1 Minamijosanjima, Tokushima, 770-8506 Japan E-mail: {nariai,minhhai,uwate,nishio}@ee.tokushima-u.ac.jp

Abstract Synchronization phenomena have been studied in various fields since long time ago. In this study, I propose a novel coupled oscillatory system such as two rings of oscillators coupled by resistors. The first ring con- sists of three van der Pol oscillators connected by resistors. The second ring consists of three oscillators connected by inductors and resistors. I focus on coupling strength of coupled oscillatory networks consisting of two kinds of oscillators. I investigate how to change synchronization phenomena by changing the coupling strength.

Key words van der Pol oscillator, Synchronization

1. ま え が き

同期現象は、非線形システムにおける最も代表的な現象であ り、自然界においてもホタルの点滅や心臓の鼓動などと同期現 象によって発生する現象が多く存在する。そのため古くから同 期現象は、工学や生物学などの様々な分野において研究されて いる。その中でも、

van der Pol

発振器の同期現象は、周波数 を変えることで様々な自然現象に近似することができる。さら に

van der Pol

発振器の回路構造はシンプルで

van der Pol

発 振器の結合系は扱いやすく、多くの研究者が

van der Pol

発振 器の様々なネットワークを提案している

[1]

[3]

。しかし

1

種 類のネットワークでの同期現象は今まで研究されてきたが、

2

種類のネットワークを持つ研究の報告は少ない。自然界におい ても

1

つのネットワークだけで成り立っているものは少なく、

様々なネットワーク感の相互作用によって成り立っている。私 たちは異なる同期現象が発生するように結合した

2

種類の

van

der Pol

発振器のネットワークを用いて研究を行う。

本研究では、結合構造が異なる

2

種類のリング状結合

van der Pol

発振器を抵抗で結合して同期現象の観測を行った。

1

つ目 のリングは、

3

個の

van der Pol

発振器を抵抗で結合したリン グ状結合

van der Pol

発振器を使用する。

2

つ目のリングは、

3

個の

van der Pol

発振器の構造を少し変え、インダクタを介し

て抵抗で落とす結合をしたリング状結合

van der Pol

発振器を 使用する。この

2

つのリング状結合

van der Pol

発振器間を抵 抗で結合させ、その抵抗の結合強度を変えた場合のの同期現象 の調査を行う。

2. システムモデル

図1に本研究で使用する

van der Pol

発振器の回路構造を示 す。この回路は負性抵抗、インンダクタ、キャパシタから構成 されている。

2

に提案する

van der Pol

発振器のネットワークモデルを

— 1 —

(2)

示す。

3

個の

VDP

を抵抗で結合した同相同期が発生する

1

つ 目のリングと、

3

個の

NC

をインダクタを介して抵抗で落とす 結合をした

3

相同期が発生する

2

つ目のリングを

3

つの抵抗

R

1

, R

2

, R

3で結合する。

R

1

, R

2

, R

3の値を変化させた場合、

同相同期が発生するネットワークと

3

相同期が発生するネット ワークを結合させたネットワークの同期現象がどのように変化 するのかを観測する。

C L

i gn

v Cn

i Ln

VDP

1 van der Pol

発振器.

circuit 2

circuit 3

circuit 6 circuit 4

circuit 1 VDP

VDP

VDP R

R2

R3

2L 2L 2L

2L 2L

2L i1a

i1b i3a

i3b i2a i2b

i4a i4b

i5b

i6a i6b i5a ir4

ir5

ir3 ir1

ir2

ir6

R R

R1

circuit 5

NC NC

NC R'

R'

R'

i

gn

v

Cn C NC

2

システムモデル.

この回路の

1

つ目のリングの回路方程式を式

(1)

に示す。

 

 

 

 

 

 

 

−i

gn

i

cn

i

Ln

= i

an

+ i

bn

+ i

rn

L

didtLn

= v

n

v

n

v

i

= (i

an

i

b(i)

)R v

n

v

j

= (i

bn

i

a(j)

)R v

n

v

n+3

= (i

rn

i

r(n+3)

)R

3

(1)

ここで、

n = 1, 2, 3

i = 2, 3, 1

j = 3, 1, 2

とする。

この回路の

2

つ目のリングの回路方程式を式

(2)

に示す。

 

 

 

 

 

 

 

 

−i

gn

i

cn

= i

an

+ i

bn

+ i

rn

v

n

2L

didtan

= (i

an

+ i

b(i)

)R

v

n

2L

didtbn

= −(i

bn

+ i

a(j)

)R

v

n

v

n−3

= (i

rn

i

r(n−3)

)R

3

(2)

ここで、

n = 4, 5, 6

i = 5, 6, 4

j = 6, 4, 5

とする。

非線形抵抗の

I V

特性を式

(3)

に示す

:

i

gn

= g

1

v

n

+ g

3

v

3n

(3)

また、変数変換に使用するパラメータを以下のように置くこ とによって

t =

LGτ, v

n

= √

g1

3g3

x

n

, i

n

=

g1C 3g3L

y

n

, ε = g

1

C

L

, α =

R1

L

C

, β = R

C

L

, γ

n

=

R1

n

L

C

.

(4)

1

つ目のリングの正規化方程式式

(5)

2

つ目のリングの正規化 方程式

(6)

が得られる。

 

 

 

 

 

 

˙

x

n

= ε(x

n

x

3n

) y

n

γ

n

(x

n

x

n+3

) +α(−x

n

+ x

i

+ x

j

)

˙

y

n

= x

n

.

(n = 1, 2, 3)

(5)

 

 

 

 

 

 

 

 

 

˙

x

n

= ε(x

n

x

3n

) y

an

y

bn

+ γ

n−3

(x

n

x

n−3

)

˙

y

an

= x

n

β(y

an

+ y

b(i)

)

˙

y

bn

= x

n

β(y

bn

+ y

a(j)

)

(n = 4, 5, 6)

(6)

パラメータ

ε

は非線形強度を意味し、

α, β,γ

nは抵抗

R

、抵抗

R

、抵抗

R

nの結合強度を意味する。

3. シミュレーション結果

提案するシステムモデルの

γ

nの値を変化させた時の同期現 象の変化に着目して研究を行った。他の結合強度は

ε = 0.05, α = 0.05, β = 0.05

に固定する。

γ

1

= γ

2

= γ

3

= 0.02

の時の発振器間の位相差を図

3

、図

4

に 示す。図

3

、図

4

の結果は初期値が異なる。これら

2

つの結果 により、初期値によって同期現象が変化することが確認できた ので、結合強度

γ

1

, γ

2

, γ

3のうち

1

つまたは

2

つの結合強度を 弱め、初期値によって変化しないように同期現象を制御する。

γ

1

= 0.001, γ

2

= 0.001, γ

3

= 0.02

の場合の発振器間の位相 差を図

5

に示す。初期値に関係なく、

1

つ目のリングでは同相 同期が観測できる。図

6

3

つの

NC

の電圧の時系列波形を示 す。この時系列波形と発振器間の位相差により、

2

つ目のリン グでは

3

相同期が発生する。これらより、

3

つのうち

2

つの結 合強度を弱めた場合は

1

つ目のリングでは同相同期、

2

つ目の リングでは

3

相同期が発生する。

γ

1

= 0.02, γ

2

= 0.005, γ

3

= 0.02

の場合の発振器間の位相差 を図

7

に示す。初期値に関係なく、

1

つ目のリングの

VDP

では それぞれ同相同期が発生し、

2

つ目のリングの

NC

では

circuit 4

circuit 6

でも同相同期が発生する。この結果より、

3

つの うち

1

つの結合強度を弱めた場合は

1

つ目のリングでは同相同

— 2 —

(3)

x 3 x 3

x 1

x 1

x 2 x 2

x 4 x 5 x 6 x 3 x 2

x 1 x 5

x 5 x 4

x 4 x 6

x 6

3

発振器間の位相差

(

γ1

=

γ2

=

γ3

= 0.02 ).

x 2

x 2 x 1

x 1 x 3

x 3

x 1 x 2 x 3 x 6 x 5

x 4

x 6 x 6 x 4

x 4 x 5

x 5

4

発振器間の位相差

(

γ1

=

γ2

=

γ3

= 0.02 ).

x 2

x 2 x 3

x 3 x 1

x 1

x 1 x 2 x 3 x 6 x 5

x 4

x 4 x 4

x 5

x 5 x 6

x 6

5

発振器間の位相差

(

γ1

= 0.001, γ

2

= 0.001, γ

3

= 0.02 ).

x n

time

6

電圧の時系列波形.

期、

2

つ目のリングの

2

つの

NC

間で同相同期が発生する。

結合強度を変えることによって、同じシステムモデルで様々 な同期現象を観測することができた。これらの結果より結合強 度を決めることによって同期現象を制御することができる。さ らに結合強度を変えることによって、

2

つ目のリングの同期現 象が変化することが分かった。

x 2 x 2 x 3

x 3 x 1

x 1

x 4 x 5 x 6

x 3 x 2

x 1 x 5

x 5 x 6

x 6 x 4

x 4

7

発振器間の位相差

(

γ1

= 0.02, γ

2

= 0.005, γ

3

= 0.02 ).

次に、

2

つ目のリングの位相差と結合強度の関係に着目して 調査する。図

8

に結合強度

γ

2と

circuit 4

circuit 5

circuit 5

circuit 6

circuit 6

circuit 4

の位相差の関係を表した グラフを示す。

γ

2の値を

0.000

から

0.030

まで

0.001

間隔で変 化させた場合の

2

つ目のリングの位相差を観測する。この時の

γ

2以外の結合強度

γ

1

= 0.02, γ

3

= 0.02

、初期値は一定の値に 固定する。

γ

2の値が

0.02

より小さい

(

他の

2

つの結合強度より小さい

)

場合について考える。

γ

2が

0.000

から

0.012

まで

circuit 4

circuit 5

circuit 5

circuit 6

の位相差は逆相同期が発生す る。

0.12

を超えると、

circuit 4

circuit 5

circuit 5

circuit 6

は位相差が

60

付近にしだいに近づく。

circuit 6

circuit 4

においては、

γ

2が

0.02

より少しでも小さければ必ず同相同期 が発生する。

次に

γ

2の値が

0.02

より大きい

(

他の

2

つの結合強度より大 きい

)

場合について考える。

γ

2が大きくなるに従って、

circuit 4

circuit 5

circuit 5

circuit 6

の位相差は

60

circuit 6

circuit 4

の位相差は

120

に近づく。

8 2

つ目のリングにおける位相差と結合強度の関係.

9,10

γ

1

= 0.02, γ

2

= 0.03, γ

3

= 0.02

時の

1

つ目のリン グ、

2

つ目のリングの時系列波形を示す。図

9

より、

1

つ目の

— 3 —

(4)

リングの発振器では同相同期が発生する。図

10

より、

2

つ目の リングでは

3

相同期が

1

つ目のリングの同相同期に寄せられた

60

ずれの同期現象が観測できた。

time x n

9 γ

1

= 0.02, γ

2

= 0.03, γ

3

= 0.02

時の

1

つ目のリングの時系列 波形.

time x n

10 γ

1

= 0.02, γ

2

= 0.03, γ

3

= 0.02

時の

2

つ目のリングの時系列 波形.

4. ま と め

本研究では提案した

2

つのリング状結合

van der Pol

発振器 における、結合強度を変えることにより同期現象の調査、初期 値によって変わらないようにする制御を行った。

3

つの結合強 度が全て等しい時は同期現象が初期値によって変化する。しか し

3

つのうち

1

つの結合強度を弱くした場合では、初期値に関 係なく

1

つ目のリング状結合発振器では同相同期が発生し、

2

つ目のリング状結合発振器では

2

つの

van der Pol

発振器で同 相同期がおこる。

3

つのうち

2

つの結合強度を弱くした場合で は、初期値に関係なく

1

つ目のリング状結合発振器では同相同 期が発生し、

2

つ目のリング状結合発振器では

3

相同期または

3

相同期が同相同期に寄せられた同期現象が観測することがで きる。

今後の展望として、提案システムモデルの回路実験を行いた い。また、リング状ではなく様々な構造における

2

種類のネッ トワークにおける同期現象を観測していきたい。

文 献

[1] M. A. Barr´ on, M. Sen, “Synchronization of four coupled van der Pol oscillators" ,Nonlinear Dynamics, pp. 357-367, June 2009.

[2] Y. Uwate, Y. Nishio, R. Stoop, “Synchronization in Three Coupled van der Pol Oscillators with Different Coupling Strength", Proceedings of RISP International Workshop on Nonlinear Circuits, Communications and Signal Processing (NCSP’10), pp. 109-112, Mar. 2010.

[3] K. Matsumura, T. Nagai, Y. Uwate, Y. Nishio, “Analy- sis of Synchronization Phenomenon in Coupled Oscillator Chains", Proceedings of IEEE International Symposium on Circuits and Systems (ISCAS’12), pp. 620-623, May 2012.

[4] H. Kume, Y. Uwate, Y. Nishio, “Synchronization Phenom- ena in Different Sizes of Rings of Coupled Oscillators", Pro- ceedings of International Symposium on Nonlinear Theory and its Applications (NOLTA’14), pp. 369-372, Sep. 2014.

— 4 —

図 1 van der Pol 発振器.
図 8 2 つ目のリングにおける位相差と結合強度の関係.

参照

関連したドキュメント

Several Synchronization Modes in Coupled Multi-State Chaotic Circuits 和田 昌浩 1. Masahiro Wada 西尾

Ogucht, “Stability and Synchronization in Networks of Time-delayed Coupled Nonlinear Systems” (¡Special Is- sue¿Network Dynamics and its Applications), Institute of Systems, Control

Abstract In this paper, an approach for clustering on complex networks by using synchronization phenomena of coupled chaotic circuits is proposed.. The proposed method is based

Abstract In this study, we observe synchronization on Hindmarsh Rose model which produce spiking-bursting behavior like real biological neurons.. There is phenomenon

[1] K.Ago, Y.Uwate, Y.Nishio, “Investigation of Synchro- nization in Coupled Chaotic Circuit Network with Local Bridge”, IEEE Workshop on Nonlinear Networks December 12-13, 2014.

Nishio, “Investigation of Synchro- nization in Coupled Chaotic Circuit Network with Local Bridge”, IEEE Workshop on Nonlinear Networks December 12-13, 2014.

Nishio,“Chaos Propagation in a Ring of Coupled Circuits Generating Chaotic and Three-Periodic Attractor“ Proceedings of IEEE Asia Pacific Conference on Circuits and Systems

We confirm that the coupled chaotic circuits network are formed several clusters which are defined by using chaos synchronization.. Finally, we apply this proposed chaotic