The eighth homework-Routes to Chaos: Period Doubling




1. Abstract



We have seen that at low driving forces the damped, nonlinear pendulum exhibits simple oscillatory motion, while at high drive it can be chaotic. However, when we consider this question again, the question must come out, that is, exactly how does the transition from simple to chaotic behavior take place? This is the main problem that we will deal with in the following content. Of course, we also will disscuss when the driving force is larger than 1.2, what will happen then?



I'm so sorry that the vpython didn't work on my computer. So, I have to use upperclassman's GIF.



Question 3.18

Calculate Poincare section for the pendulum as it undergoes the period-doubling route to chaos. Plot \omega versus \theta, with one point plotted for each drive cycle, as in Figure 3.9. Do this for F_D =1.4, 1.44, 1.465, using the other parameters as given in connection with Figure 3.10. You should find that after removing the points corresponding to the initial transient the attractor in the period-1 regime will contain only a single point. Likewise, if the behavior is period n, the attractor will contain n discrete points.


Question 3.20

Calculate the bifurcation diagrams for the pendulum in the vicinity of F_D=1.35 to 1.5. Make a magnified plot of the diagram (as compared to Figure 3.11) and obtain an estimate of the Feigenbaum \delta parameter.




2. Background and Introduction



At last homework, we have ever taked about the phenomenon-chaos and we have disscussed a example of Poincare section-strange attractor. When the system is chaos, assume that



Quickly, it's Poincare section just like below.



But now, we will introduce more examples of Poincare section.

Bifurcation diagram is very helpful to analyze the transition to chaos. It can show us lines for \theta as a function of drive amplitude, which was constructed in the following manner. For each value of F_D we have calculated \theta as a function of time. After waiting for 300 driving periods so that the initial transients have decayed away, we plotted \theta at times that were in phase with the driving force as a function of F_D. Here we plotted points up to the 428th driving period.


After that, we can calculate the Feigenbaum \delta parameter through following formula



In theory, when n levels off to infinity, \delta=4.669.




3. Body Content and Conclusion



1. Period Doubling with different values of the drive amplitude.



In the last time, we just gave the results when F_D is smaller than 1.2, now, I'll give results for \theta as a fuction of time for our pendulum for several different values of the drive amplitude.
Here Click the Code!

Different drive amplitude

Surprisingly, when F_D=1.2, the system is in chaotic state. When F_D=1.35, the system become well-aligned again, its period is the same as the drive period. When F_D=1.44, its period is the twice as the drive period. When F_D=1.465, its period is the fourth times as the drive period.

  • Conclusion: If we were to increase the driven amplitude further, the period would double again as the pendulum would switch to a motion that has a period eight times that of the drive. The period-doubling cascade would continue if the drive were increased further.


2.Question 3.18



If the time, which we plot, satisfys conditions below, from above figure, there must be some points instead of lines.
Click the Code


F_D=1.4

F_D=1.44

F_D=1.465

  • Analysis: As the same as 1., when F_D=1.4, there is just one point. When F_D=1.44, there are two points. When F_D=1.465, there are four points. Be careful, my each figure has two points that are used to mark conditions!
  • Conclusion: from figure, it quite easy find that after removing the points corresponding to the initial transient the attractor in the period-1 regime will contain only a single point. Likewise, if the behavior is period n, the attractor will contain n discrete points.

Click the Code

Combination Comparision

  • In order to understand it more directly, I just give a figure that we plot just in the drive period. And above four diagrams were ploted when F_D=1.2, F_D=1.4, F_D=1.44 and F_D=1.465.

In fact, if we didn't set the values of each axis, one point will be a line, why is it happens?


Calculational Error
  • As a matter of fact, computer can't give accurate \pi, for example, it may give 3.141592654. Thus, in a microscopic regime, it may not a point, just like what I showed you above.


3. Question 3.20



Bifurcation diagram is a quite good method to tell us the transition to chaos. I spend a flood of time on operating the program, it need a long time to calculate. If you find that it seems that my code can't work, don't worried or surprised about, just wait for the results. It is about several minutes to several hours. Of course, it depend on your accuracy.
Click the Code

Bifurcation Diagram

  • Conclusion: For each value of F_D we have some values of \theta, just like question 1. and queation 2., it is obvious that when F_D=1.4, there just is one value of \theta, when F_D=1.44, there are two branchs. There may are more branchs when F_D is larger.

In this part, the other probelm is to find the Feigenbaum \delta. It requires high accuracy. So, I spend almost an hour to get the result.


More Accurate Bifurcation Diagram

A Magnified plot of the diagram1

A Magnified plot of the diagram2
A Magnified plot of the diagram3
  • Make some magnifield plots of the diagram, in this way, we can find the points that from period 2^(n-1) to period 2^(n).
  • Just like figures above, we can get some some branch-points, such as, F_D1=1.4228, F_D2=1.45841, F_D=1.47505, F_D=1.47616, F_D=1.476425, F_D=1.476481, F_D=1.476493.
  • Calculate them, we can obtain 2.14, 14.991, 4.189, 4.732, 4.667. If we continue this process, we will get the Feigenbaum \delta approximate to 4.669.



4. Reference and Acknowledgement

[1] Computational Physics (Second Edition), Nicholas J. Giordano, Hisao Nakanishi;
[2] Upperclassman-Yuqiao Wu
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 201,312评论 5 473
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 84,578评论 2 377
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 148,337评论 0 333
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,134评论 1 272
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,161评论 5 363
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,303评论 1 280
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,761评论 3 393
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,421评论 0 256
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,609评论 1 295
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,450评论 2 317
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,504评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,194评论 3 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,760评论 3 303
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,836评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,066评论 1 257
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,612评论 2 348
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,178评论 2 341

推荐阅读更多精彩内容

  • **2014真题Directions:Read the following text. Choose the be...
    又是夜半惊坐起阅读 9,311评论 0 23
  • 上一章 下一章 首页 瓜多不再是新毅的粉丝,新毅很想抚去瓜多的伤痛,但瓜多总是有意无意的退避着,这世间最难愈合的...
    GentleC阅读 214评论 0 3
  • 记得女朋友当初给我到单位送饭,让我很吃惊。因为在我眼中,她有些懒且任性,便好奇问她为何要给我大老远送饭。她说想趁爱...
    托尔西奇阅读 218评论 0 0
  • 愿我们还能看见"黄河之水天上来,奔流到海不复还。万里长城永不倒,千里黄河水滔滔"的壮景。我希望更多人学习黄河的伟大...
    emmmmmm哦阅读 123评论 0 0
  • 在最近的一段时间里,我渐渐地对阅读有了兴趣。 当我可以安静下来读一本书时,仿佛觉得自己在与作者对话;瞬间感觉回到了...
    冰磊阅读 255评论 0 0