2021 Computer Simulation

数物科学科 2021年後期計算機シミュレーション (木3-4) [情報基盤セ3F301]

担当者(所属)：高橋信介
(理工)： taka_lec耳hirosaki-u.ac.jp ←耳を@に置換えて下さい

【レポート提出状況】

自由提出レポート課題や提出締切についてはTeamsのファイルタブにアップしてあります。

第01回 Oct.07	第02回 Oct.14	第03回 Oct.21	第04回 Oct.28
第05回 Nov.04	第06回 Nov.11	第07回 Nov.18	第08回 Nov.25
第09回 Dec.02	第10回 Dec.09	第11回 Dec.16	第12回 Dec.23
第13回 Jan.13	第14回 Jan.20	第15回 Jan.27

Linux, Ubuntu
login, logout, gedit, gfortran, gnuplot, e-mail
ls, time, cd, pwd, cat
プリントpage1～2
fortranをチョット復習
この授業の受講を希望する場合は taka_lec 耳 hirosaki-u.ac.jp 宛にメールを送信．
gnuplotをちょっと練習

次回のプリント(page11～22)を配布しました

落下運動の運動方程式を解く(2階の微分方程式を解く)
次のファイルをダウンロード（或いはこれをコピペ）して利用します。

空気抵抗を無視した場合の放物運動(1次元)
数値解と解析解を比較する
速度の1乗に比例する空気抵抗を考慮した場合の放物運動(1次元)
空気抵抗を考慮する場合
sp140.f は sp131.f を読み込み，修正・編集・保存して利用します.
gnuplotを繰り返し練習しておきましょう。

空気抵抗を考慮しない放物運動(2次元)

速度の2乗に比例する空気抵抗を考慮した場合の放物運動(2次元)
　Euler Richardson法
（プログラムは何か所か消してある部分があります）

sp200.f
sp210.f
sp211.f
sp200.plt
sp200.eps
未提出・訂正レポートは遅れても必ず提出すること.
追加や修正を提出する場合は，どこを，追加や修正したか
を必ず説明し，その部分だけを送信すること.

単振動，減衰振動，強制振動

計算の方針を詳しく説明する必要あり。

自分の直感

ただしk=0, k2=0とすればそれぞれ抵抗無しの計算として比較できますね

レポート提出の際は提出のルールを守る(プリントpage2参照)

数百行数千行もメールに載せて送信しない

単位のある数値には単位を必ず付ける。

未提出・訂正レポートは遅れても必ず提出すること.
追加や修正を提出する場合は，どこを，追加や修正したか
を必ず説明し，その部分及びその関連部分だけを送信すること.

RC回路（電気回路の過渡現象）
キルヒホッフの法則を使います

LRC回路
- sp500.f (sp401.fのコピペ編集必要)
- sp510.f (sp401.fのコピペ編集必要)
- sp500.plt
- sp501.plt
- sp510.plt
- sp510.eps

第10回 Dec.09
ここから確率論的シミュレーション
一様乱数の使い方を復習しておきましょう

一様乱数
度数分布と確率分布
- sp610.f
- sp610.plt

第11回 Dec.16　
遅れています。今日は「一様乱数」「度数分布と確率分布」プリント30ページまででした。

【レポート課題３】を出しました。
【配布資料】

「一様乱数の分布」のばらつきを調べる。

第12回 Dec.23
一様乱数を使って他の確率関数に従う乱数をつくる

指数乱数（逆関数法）
正規乱数

第13回 Jan.13, 2022
実習室で対面授業
観測データの解析１

演習
- 070109.DAT 観測データ(20000events)
- sp800.f (ASデータファイルを読む．時間差計算)
- sp800.plt (飛来時間差の図)
- sp805.f (1時間毎の飛来数の推移)
- sp805.plt (1時間毎の飛来数の推移の図)
- sp810.f (飛来時間差の分布)
- sp810.plt (飛来時間差の分布 ===> 指数分布)

第14回 Jan.20
まとめ（Web授業）

1/13に配布したプリントはpdfファイルをTeamsにアップしてあります。
私有パソコンにfortranをインストールする動画及びpdfファイルをTeamsにアップしてあります。参考にしてください。
9分18秒付近から再生してください。

第15回 Jan.27(休講)　→　Feb.03

期末試験（Web）

第16回 Feb.3　→　Feb.05 (Web)

【レポート提出状況】
r 印：レポ－ト受信，_印：未着信，/印：不備有 or/and 説明不足，
G印：Good and/or がんばった，O印：OK(良以上)，o印：可，x印：再提出,
V印：Very good and/or がんばった，

A：Abstract（rep3-1, rep3-2, rep3-3 全体の要旨），
p：計算処理の方針(policy)の説明 and/or 修正・作成したprogram及びその説明，
k：計算結果(考察に必要な部分及びその前後の部分のみとその説明)，
f ：考察に必要な図(.PNG形式のファイル)とその説明，
K：考察(←これが大切，文章は良く推敲する)

Feb.13, 2022. 18:00　（これ以降更新しません。以降の受信はメールで返信します）
自由提出レポート課題や提出締切についてはTeamsのファイルタブにアップしてあります。

rep#	0	rep1				rep2				rep3-1,2,3				_				_
出題	1	10/21				11/18				12/16				/				/_
締切	_	10/28, 10:20				11/25, 10:20				12/23, 10:20				/ , :				/_ , :
pkf	m	04	08	_	_	p	k	f	K	A,p	k	f	K	p	k	f	K	p	k	f	K
a05	O	O	o	_	_	r	r	o	o	_,_	_	_	_	_	_	_	_	_	_	_	_
a18	O	O	O	_	_	_	_	_	_	_,_	_	_	_	_	_	_	_	_	_	_	_
a46	O	O	_	_	_	_	_	_	_	_,_	_	_	_	_	_	_	_	_	_	_	_
a66	O	O	O	_	_	V	V	V	V	V,V	V	V	V	_	_	_	_	_	_	_	_
a67	/	O	o	_	_	G	G	G	G	V,V	V	V	V	_	_	_	_	_	_	_	_
a74	/	O	o	_	_	_	O	O	O	r12,r12	r12	r12	r12	_	_	_	_	_	_	_	_
a77	O	O	O	_	_	_	_	_	_	_,_	_	_	_	_	_	_	_	_	_	_	_
a79	O	O	O	_	_	r	r	r	r	_,_	_	_	_	_	_	_	_	_	_	_	_
b03	O	x	_	_	_	_	O	O	O	r1,r1	r1	_	r1	_	_	_	_	_	_	_	_
b41	O	O	_	_	_	O	O	O	O	O1r23,O1r23	O1r23	_	O1r23	_	_	_	_	_	_	_	_
b42	O	x	xx	_	_	_	_	_	_	_,_	_	_	_	_	_	_	_	_	_	_	_
b73	O	O	xx	_	_	O	O	O	O	r1,r1	r1	_	r1	_	_	_	_	_	_	_	_
b75	_	_	_	_	_	_	_	_	_	_,_	_	_	_	_	_	_	_	_	_	_	_
c02	O	O	O	_	_	_	_	_	_	_,_	_	_	_	_	_	_	_	_	_	_	_
c33	O	O	O	_	_	r	O	O	O	O1r2,O1r2	O1r2	r2	O1r2	_	_	_	_	_	_	_	_
pkf	m	04	08	_	_	p	k	f	K	A,p	k	f	K	p	k	f	K	p	k	f	K
rep#	0	rep1				rep2				rep3				_				_

taka_lec 耳 hirosaki-u.ac.jp