Upgrade to Pro
— share decks privately, control downloads, hide ads and more …
Speaker Deck
Features
Speaker Deck
PRO
Sign in
Sign up for free
Search
Search
数値シミュレーション / Python 13
Search
kaityo256
PRO
January 14, 2020
Education
5
3.6k
数値シミュレーション / Python 13
プログラム基礎同演習13
kaityo256
PRO
January 14, 2020
Tweet
Share
More Decks by kaityo256
See All by kaityo256
論文の読み方 / How to survey
kaityo256
PRO
175
120k
リンゴゲームと貧富の差 / Origin of the disparity of wealth
kaityo256
PRO
12
13k
渡辺研Slackの使い方 / Slack Local Rule
kaityo256
PRO
8
7.5k
時間の矢について / Time's arrow
kaityo256
PRO
12
17k
t-SNEをざっくりと理解 / Overview of t-SNE
kaityo256
PRO
2
520
未定義動作でFizz Buzz / Undefined Fizz Buzz
kaityo256
PRO
1
720
卒論の書き方 / Happy Writing
kaityo256
PRO
29
19k
再帰呼び出し / Python Recursion
kaityo256
PRO
0
1.5k
分子動力学法(3) 実装と高速化の詳細 / Simulation 07
kaityo256
PRO
0
640
Other Decks in Education
See All in Education
パフォーマンス・チューニング入門
oracle4engineer
PRO
2
600
自己紹介 / who-am-i
yasulab
2
3.1k
Earthquake and Disaster Prevention Information for UTokyo International Students
utokyoissr2360
0
630
MonacaEducation導入手引き2024年第1版
asial_edu
1
150
Data Presentation - Lecture 5 - Information Visualisation (4019538FNR)
signer
PRO
0
1.8k
Human Perception and Colour Theory - Lecture 2 - Information Visualisation (4019538FNR)
signer
PRO
0
1.7k
Best Wedding day perfume
vicjon
0
260
Avoin jakaminen ja Creative Commons -lisenssit
matleenalaakso
0
1.1k
人生の転機からチャンスを掴む「シュロスバーグの4Sモデル」/4s-models
yuko_yokouchi
2
720
WordPressを教える人のための視点と考え方
crebowinfo
0
230
千葉県印西市立・原山小学校における新たな学び「情報探究の時間」実践報告』
codeforeveryone
1
700
生成AIを活用できる大学教職員になる-基本と実践-
gmoriki
0
300
Featured
See All Featured
Fashionably flexible responsive web design (full day workshop)
malarkey
398
65k
Testing 201, or: Great Expectations
jmmastey
28
6.3k
Optimising Largest Contentful Paint
csswizardry
8
2.4k
Creatively Recalculating Your Daily Design Routine
revolveconf
210
11k
Atom: Resistance is Futile
akmur
259
25k
Automating Front-end Workflow
addyosmani
1356
200k
We Have a Design System, Now What?
morganepeng
43
6.7k
Intergalactic Javascript Robots from Outer Space
tanoku
266
26k
Art, The Web, and Tiny UX
lynnandtonic
289
19k
No one is an island. Learnings from fostering a developers community.
thoeni
16
2.1k
Web development in the modern age
philhawksworth
202
10k
JavaScript: Past, Present, and Future - NDC Porto 2020
reverentgeek
40
4.4k
Transcript
1 数値シミュレーション プログラミング基礎同演習 慶應義塾大学理工学部物理情報工学科 渡辺 2019/1/14
2 Keep派の正解確率は1/3 Keep派とChange派で 場合は尽きている Change派の正解確率は2/3 「Change派正解」は、「Keep派正解」の余事象 Change派の方が得
3 道が確率pで通行可能な時、左から右に通過できる確率は?
4 p 通行不能 通行可能 ある確率(0.5)を境に「ほぼ通行不能」から「ほぼ通行可能」に急激に変化する この現象を相転移、相転移する確率pを臨界点と呼ぶ
5 ・CDを買うとポスターがランダムにもらえる ・44種類全て揃えると特別なイベントに招待 ・全て揃えるのにCDを何枚買う必要があるか?
6 はじめはサクサク、後になるほど渋くなる CD購入枚数 新しいポスターを得る確率 最初の方の試行で「これは楽勝だ」と誤認させる→沼へ 100枚購入した時点で コンプ率90%程度 10枚に一度しか ポスターが増えない
7 途中であきらめるとそれまでのコストが無駄に 100枚購入した時点で、「あれ?」と気が付く しかし、手元には「ほぼそろっている」ポスター 今諦めたら、100枚購入した投資が無駄になる 沼へ ※ サンクコスト (撤退時に回収できない費用)効果
8 Q. コンプガチャは違法にすべきなの? A. 我々がどのような世界を目指すかに依る 嘘がなければ高金利は許されるか? どこまでを「自己責任」とするか? ギャンブルは禁止すべき?制限すべき? スマホのn年縛りは? 答えはない。我々が相談し、決めるしかない
9 数値シミュレーション • 差分化 • 拡散方程式 • 反応拡散方程式
10 あるルールに従う系の振舞を 再現・予想すること それをコンピュータ上でやるのがコンピュータシミュレーション
11 訓練のためのシミュレーション フライトシミュレータ ドライブシミュレータ 娯楽のためのシミュレーション シミュレーションゲーム (戦略、経営等) ロールプレイングゲーム 予測のためのシミュレーション 天気予報・自動車の設計・創薬
12 物理とは 我々が住むこの世界を理解・記述する学問 この世界のルール この世界は微分方程式で記述されている これを支配方程式(Governing Equation)と呼ぶ
13 古典的な物体の運動を記述する = ∙ = , ∙ = 0 ×
= − , × − = ニュートンの運動方程式 = 1 + 流体の流れを記述する ナビエ・ストークス方程式 電場や磁場の振る舞いを記述する マクスウェル方程式 例:天体の運動、弾道計算 例:天気予報 例:スマートフォンの開発等
14 = 2 2 拡散方程式 物質やエネルギーなどが広がって (拡散して)いく様子を表現する方程式
15 = 2 2 拡散方程式 温度の時間変化を表す ある点での温度が上がるか下がるか? あがるならどれくらい上がるか?
16 = 2 2 拡散方程式 ある地点での二階微分の値 = () 2 2
> 0 2 2 < 0 下に凸なら正 = () 上に凸なら負
17 = 2 2 拡散方程式 まわりの平均との差を表す 下に凸=自分は平均より低い 平均 平均 上に凸=自分は平均より高い
18 拡散方程式 まわりを見て 自分が平均以下なら増えようとする 自分が平均以上なら減ろうとする 自分 自分は平均より低い 平均 = 2
2
19 支配方程式 = 知りたい現象を記述する微分方程式 これを解けば未来がわかる ほとんどの微分方程式は厳密に解くことができない 数値的に近似解を求める コンピュータシミュレーション
20 コンピュータは 離散的値しか扱えない この世界は連続的 計算機が扱えるように連続的な値をとびとびの値にすることを離散化と呼ぶ
21 = − 2 2 離散化には空間の離散化と時間の離散化がある 時間 空間
22 空間の離散化 時間の離散化 拡大するとピクセルに 静止画像を高速コマ送り 我々が計算機を通して目にするものは離散化されている
23 連続的な世界 離散的な世界 この領域全体の物理量を この点での値で代表させる
24 微分を差分で近似すること(離散化の一種) + ℎ = + ℎ′ + ℎ2 テイラー展開を一次まで考える
二次以上を無視する ≈ + ℎ − () ℎ ′ について解く 微分が差分で近似された
25 O Time t () 時刻tにおける傾き O Time t ()
≈ + −() t+dt 時間変化=現在と少し未来の差
26 + ℎ = + ℎ′ + ℎ2 2 ′′()
テイラー展開を二次まで考える − ℎ = − ℎ′ + ℎ2 2 ′′() 2 2 ≈ + ℎ − 2 + ( − ℎ) ℎ2 ′ を消去 二次の中央差分 二階微分
27 = − 2 2 ある場所の時間変化量は まわりの平均との差をへらそうとする この地点での次のステップでの値を 周りの値をみて決める 上記の操作をすべての地点について繰り返すと
次のステップ(少し未来)での「世界」がわかる 「ステップ」を繰り返せば、遠い未来の世界がわかる
28 Belousov-Zhabotinsky (BZ) 反応 いくつかの溶液を混ぜる 色が周期的に変わっていく 周期的な化学反応が拡散と組み合わさると?
29 二種類の化学物質uとvがお互いに反応しながら拡散する式 = ∆ − 2 + (1 − )
= ∆ + 2 − + 拡散 反応 チューリング・パターン 反応と拡散がおりなす不思議な模様