計(jì)算物理第9次作業(yè)

Name: 賀一珺
Student Number: 2014302290002

Question

  • 3.30 Investigate the Lyapunov exponent of the stadium billiard for several values of \alpha. You can do this qualitatively by examing the behavior for only one set of initial conditions for each value of \alpha you consider, or more quantitatively by averaging over a range of initial conditions for each value of \alpha.

  • 3.31 Study the behavior for other types of tables. One interesting possibility is a square table with a circular interior wall located either in the center, or slightly off-center. Another possibility is an elliptical table.

Abstract

In last homework I investigated the process to chaos, using the bifurcation diagram, which is quite stimulating. Now I will get closer to the world of chaos. I will show you two totally different systems which will also lead you to the real world which is full of chaos. They are Lorenz model which comes from weather problems and the billiard problem which lays in our daily life. There are also some new tricks that you must pay attention to when doing your progremming and I will give you a brief introduction to those small tricks. Let's prolong our legs kicking tricks together! I hate tricks! Also, I will show you the trejectory and phase space of the billiard problem which in my opinion will certainly makes you upsurge!

Background

Lorenz system[1]

The Lorenz system is a system of ordinary differential equations (the Lorenz equations, note it is not Lorentz) first studied by Edward Lorenz. It is notable for having chaotic solutions for certain parameter values and initial conditions. In particular, the Lorenz attractor is a set of chaotic solutions of the Lorenz system which, when plotted, resemble a butterfly or figure eight.

Intermittent_Lorenz_Attractor_-_Chaoscope

The model is a system of three ordinary differential equations now known as the Lorenz equations:


Here x, y, and z make up the system state, t is time, and σ, ρ, β are the system parameters. The Lorenz equations also arise in simplified models for lasers, dynamos, thermosyphons, brushless DC motors, electric circuits, chemical reactions and forward osmosis.

Billiard ball[2]

A billiard ball is a small, hard ball used in cue sports, such as carom billiards, pool, and snooker. The number, type, diameter, color, and pattern of the balls differ depending upon the specific game being played. Various particular ball properties such as hardness, friction coefficient and resilience are important to accuracy.

Billiards_balls

Crutial points

  • The system gives you a path to chaos. Thus you have to minimize every possible error. When handeling the billiard problem, you have to take care of boundary carefully. Since we devided the total time into many small parts, the volue on the boundary are actually skipped. What should we do is to add code and approxe the boundary to ecrease the error.
  • The hardest problem I meet in this homework is how to progrem the boundary condition correctly. The solotion could be approached by adjust the variables in the situation of simple circle.

Plotting

Before the begining

Lorenz Model

I tried out Lorenz model in our textbook. Here is initial condition:
Time step = 0.0001
σ = 10
b = 8/3
r = 25

Lorenz Model z versus time
3D trejectory

I adjusted the volue of r and plotted again. In this picture the volue of r are 25, 10, 5, retrospectively.

z versus time

It can be derived from this picture that when r = 25 the system is in the state of chaos actually.

The Billiard Problem

I also tired out some simple billiard problems in this section, which include a billiard on a square table and a billiard on a circle table.

Square table

Initial Condition:
Time step = 0.01
x = 0.2
y = 0
vx = 1
vy = 0.6

trajectory of a billiard on a table
phase space plot
Circular table

Initial condition:
x = 0.2
y = 0
vx = 1
vy = 0.6

trajectory of a billiard on a table
phase space

It's beautiful, isn't it? Now let's begin our work.

Problem 3.30

Plot the trajectory on four possible table first. Note that I used the trick mensioned before here. The initial condition is the same as last plot.

the trajectory on four possible table
phase space

Here is the seperation versus time when alpha is 0:

alpha = 0

Here is the seperation versus time when alpha is 0.001:

alpha = 0.001

Here is the seperation versus time when alpha is 0.01:

alpha = 0.01

Here is the seperation versus time when alpha is 0.1:

alpha = 0.1

Problem 3.31

Square table with a circular interior wall
trajectory of a billiard on a table
phase space
Elliptical table
trajectory of a billiard on a table
Phase space

A 3D version by vpython

Circular table

Conclution

I investigated two more systems in this homework, but there are still some problems. The first is that the phase space of 3.30 is actually wrong. However, it's hard for me to figure out the crutial points. Also, I got some stimulating results in this homework. These two sytems exhibit all charicters of chaos.

Acknowledgement

[1]wikipedia
[2]Shuxu Hu(胡墅緒) for that he helped me using vpython to work out 3D diagrem
[3]Shan Tan(譚善) for his idea on velocity calculating near boundary
[4]Prof. Cai

How to contact me

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末境蜕,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子胖腾,更是在濱河造成了極大的恐慌,老刑警劉巖,帶你破解...
    沈念sama閱讀 216,651評(píng)論 6 501
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異烹玉,居然都是意外死亡卖宠,警方通過(guò)查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,468評(píng)論 3 392
  • 文/潘曉璐 我一進(jìn)店門(mén)坪稽,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)曼玩,“玉大人鳞骤,你說(shuō)我怎么就攤上這事∈蚺校” “怎么了豫尽?”我有些...
    開(kāi)封第一講書(shū)人閱讀 162,931評(píng)論 0 353
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)顷帖。 經(jīng)常有香客問(wèn)我美旧,道長(zhǎng)贬墩,這世上最難降的妖魔是什么榴嗅? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 58,218評(píng)論 1 292
  • 正文 為了忘掉前任,我火速辦了婚禮陶舞,結(jié)果婚禮上录肯,老公的妹妹穿的比我還像新娘。我一直安慰自己吊说,他們只是感情好论咏,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,234評(píng)論 6 388
  • 文/花漫 我一把揭開(kāi)白布。 她就那樣靜靜地躺著颁井,像睡著了一般厅贪。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上雅宾,一...
    開(kāi)封第一講書(shū)人閱讀 51,198評(píng)論 1 299
  • 那天养涮,我揣著相機(jī)與錄音,去河邊找鬼眉抬。 笑死贯吓,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的蜀变。 我是一名探鬼主播悄谐,決...
    沈念sama閱讀 40,084評(píng)論 3 418
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼库北!你這毒婦竟也來(lái)了爬舰?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 38,926評(píng)論 0 274
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤寒瓦,失蹤者是張志新(化名)和其女友劉穎情屹,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體杂腰,經(jīng)...
    沈念sama閱讀 45,341評(píng)論 1 311
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡垃你,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,563評(píng)論 2 333
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片惜颇。...
    茶點(diǎn)故事閱讀 39,731評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡皆刺,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出官还,到底是詐尸還是另有隱情芹橡,我是刑警寧澤毒坛,帶...
    沈念sama閱讀 35,430評(píng)論 5 343
  • 正文 年R本政府宣布望伦,位于F島的核電站,受9級(jí)特大地震影響煎殷,放射性物質(zhì)發(fā)生泄漏屯伞。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,036評(píng)論 3 326
  • 文/蒙蒙 一豪直、第九天 我趴在偏房一處隱蔽的房頂上張望劣摇。 院中可真熱鬧,春花似錦弓乙、人聲如沸末融。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 31,676評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)勾习。三九已至,卻和暖如春懈玻,著一層夾襖步出監(jiān)牢的瞬間巧婶,已是汗流浹背。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 32,829評(píng)論 1 269
  • 我被黑心中介騙來(lái)泰國(guó)打工涂乌, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留艺栈,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 47,743評(píng)論 2 368
  • 正文 我出身青樓湾盒,卻偏偏與公主長(zhǎng)得像湿右,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子罚勾,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,629評(píng)論 2 354

推薦閱讀更多精彩內(nèi)容