#author("2021-01-19T14:00:55+09:00","","") [[MechanicalUniverse]] * R2年度 [#zbdee2bd] ** 1月18日の授業をオンデマンド課題で行います。 [#v807dcaf] 交流回路で交流電源、抵抗、コイルをつないだ回路のオシロスコープに表示されるグラフを作成する。最終的な抵抗とコイルの電圧降下を表示するグラフを作成するが、縦軸、横軸のないグラフの書き方を示すので、縦軸、横軸を作成せよ。 まずは簡単なグラフの書き方。 import numpy as np import matplotlib.pyplot as plt t = np.linspace(0, 2.) I = np.sin(2*np.pi*t) plt.plot(t, I) plt.show() これだと数値に意味が与えられていない。そこで2サイクル分だけ表示する。標準的な単位で示した物理量を入れたプログラムにする。 import numpy as np import matplotlib.pyplot as plt R = 100 I = 0.01 freq= 50. L=0.1 x = np.linspace(0, 2.) t = x /freq omega= 2* np.pi*freq VR = R*I*np.sin(2*np.pi*x) VL = omega * L * I *np.cos(2*np.pi*x) plt.plot(x, VR) plt.plot(x, VL) plt.show() 上に示したプログラムでは東日本の50Hzを仮定している。抵抗の振幅が1Vになるようにしている。numpyをつかうので自動的に高速に計算できる数値のタイプに変換されているようである。つぎに横軸を時間軸に変えるグラフへの変更を示すので、これに縦軸、横軸がそれぞれ電圧、時間であることがわかるようにせよ。 plt.plot(t, VR) plt.plot(t, VL) さらに電源、抵抗、コンデンサーの回路について、抵抗、コンデンサーの電圧降下のぐらふの書き方を示すので同様に縦軸、横軸について改善すること。 import numpy as np import matplotlib.pyplot as plt R = 100 I = 0.01 freq= 50. C = 1.e-4 x = np.linspace(0, 2.) t = x /freq omega= 2* np.pi*freq VR = R*I*np.sin(2*np.pi*x) VC = -I /(omega*C)*np.cos(2*np.pi*x) plt.plot(t, VR) plt.plot(t, VC) plt.show() * H29年度,H30年度 [#k0b9a79d] ** 原子物理学 [#h6eda756] 蛍光灯を見ると黄色に二本線が見えるが、この図のどれにあたるか。5460Aは緑。黄色はもう少し波長が長い。特に5770Aは2個の電子で作るスピンが変化する。でも意外に明るい。光は軌道角運動量が1だけ変化して出てくる。 http://physwiki.apps01.yorku.ca/index.php?title=Main_Page/PHYS_3220/Excitation_Potentials ただしこの蛍光管はかなりまれ。見える輝線がほぼ水銀由来というもの。消えかかった印刷を読むとpanasonic のFL18SS・Wというものらしい。50cmの直管蛍光灯ということなので、あったのが不思議なものだった。WでなくEX-Nなんかはきっと蛍光物質の輝線がいっぱい出てくる。今どきの蛍光塗料のせいでみんなこれしかない。そうなると、このあたりを見ながら水銀の輝線を確認してもよいかも。殺菌灯では水銀の線がよく見える。 http://seppina.cocolog-nifty.com/blog/2016/07/post-463a.html 水銀は2電子系だが、1電子系で有名なのはナトリウム。ただこのランプで高圧だと輝線は見えないらしい。D線は高倍率で見ると2本見える。分離は4A。 ** 6.5-7 電位、キャパシター、誘電体 [#pe18b6ed] H30年度は教科書をそのまま見ていきます。 誘電体の面白い性質として圧電素子がある、「圧電素子 原理」で検索 http://www.jiyuken.dicis.net/thunderbolt/playwithspark.html #ref(電位・キャパシター・誘電体2.pptx) ** 6.3 電場 [#nab64eaa] 琉球大学の前野昌弘さんのページにflashやjavaのプログラムがアップされている。 http://irobutsu.a.la9.jp/mybook/ykwkrEM/sim.html 興味ある人は本も買ってみるとよいかも。図書館に本はあります。 ** 5.5 熱力学第2法則 [#x6088afb] 教科書の二つの箱に分布する粒子の確率分布は2項分布になる。標準偏差はNpqの平方根。pが1/2の時、2σ/NはNの平方根となる。2項分布をgnuplotで表示させてみるには以下のページの式をつかう。階乗のlogをとると大きな数字にならないという説明で指数関数で戻してやっている。少数でも階乗もどきの計算をさせるのにガンマ関数を使っている。 #ref(binomial.gnu) http://agw.hatenablog.jp/entry/20111212/1323754951 ** 5.2 熱の移動 [#q8532da9] 熱放射についてプランクの法則をgnuplotで描画する。数字が大きくなるのを避けるため、プランク定数と光速の積を2πで割ったものを197とする電子ボルトとナノメートルの表記にした後、ボルツマン定数も電子ボルトにしてつじつまを合わせる。その後マイクロメートルに変換したり、温度を千度単位にしたりして、以下のプログラムを得た。プロットしてみよう。 #ref(plank.plt) #ref(plank2.plt) * H19年度 物理学II [#vf4b413b] 使用したパワーポイントファイル。ただし学内からのみアクセス可能。 ** 波・熱 [#v3cba66e] / 平衡と弾性&ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/balanceElastic.ppt); / 流体 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/fluid.ppt); / 振動 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/vibration.ppt); / 波動I &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/waveI.ppt); / 波動II &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/waveII.ppt); / 温度・熱・第一法則 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/firstLaw.ppt); / 気体分子運動論 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/MolecularM.ppt); / エントロピーと熱力学第2法則 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/irreversible.ppt); / ** 電磁気学 [#nda1d6d4] / 電荷&ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/denjiki/denka.ppt); / 電場 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/denjiki/denba.ppt); / ガウスの法則 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/denjiki/gauss.ppt);/ 電位 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/denjiki/denni.ppt); / 電気容量 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/denjiki/Capacitor.ppt); / 電流と抵抗 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/denjiki/currentRegistance.ppt); / 回路 &ref(http://natsci.kyokyo-u.ac.jp/~takasima/butsuIIppt/denjiki/Circuit.ppt); /