ニュートンの運動方程式
こんにちは
今回は少し簡単な記事です。
今回のモチベーションは、「なぜ運動方程式は2階の微分方程式なの?」ということです。
高校に入って物理を始める時に最初に物体の運動について勉強すると思います。
その中で学ぶのがニュートン力学(Newtonian Mechanics)です。
そこに出てくる法則が3つあります。復習で書きますね。(ニュートン力学 - Wikipedia参照)
まず前提の定義としてどんな物体であっても大きさは無視したただの点であり、その点にすべての質量が集まっているものとする。またこの点を質点という。
運動の第1法則(慣性の法則)
質点は力が作用しない限り、 静止または等速直線運動をする
質点の運動量の時間的な変化は、質点に加えられた力の大きさに比例し、方向は与えられた力の方向に依存する
運動の第3法則(作用・反作用の法則)
二つの質点があり、その間で相互作用があるときに質点1から質点2への力と質点2から質点1への力は大きさが等しく、方向が反対である
ニュートンが当時書いた『プリンキピア』の中では実は運動の第1法則と第2法則しか書いていないらしいですが。笑
ということで第1法則、第2法則を考えましょう。
ニュートンの当時の気持ちになってみましょう。
力を加えると質点の運動は変わりそうだ。ここで運動とは位置の時間変化であるので、まずは1階微分である速度が変化することを考えてみましょう。すなわち、比例係数として
となります。ここで一般化するために、力を次のように時間に依存した量として考えていきます。ただし、は定ベクトルとする。
この関係を先ほどの式に代入する。
この式での問題点は速度の式において、質点が静止している時は運動の第1法則が成り立つが、等速直線運動をしているときには成り立たないということです。すなわち、がある一定の有限値を持つ場合はがある有限値を持つ必要がある。しかし、第1法則は力が作用しない限り成り立つものなので矛盾する。
では、次に2階微分である加速度が変化することを考えてみましょう。比例係数はとします。
これにより先ほどの問題点が解決され、物質の運動を表現できそうです。
では、次に3階微分である躍度(加加速度)が変化することを考えてみましょう。比例係数はとします。
今度は力がなかったとしても速度が一定の有限値にならないという問題が生じてしまっています。ということは3階微分では表せないということです。同様に4階以上でも表せません。
つまり運動方程式は位置の時間に関する2階微分である加速度は力を加えることで変化するという方程式です。
今回はここがゴールなのでここまでにします!
2017-12-30
積分と和分
こんにちは
では、今日も少し数学的なお話をしますね。
積分と和分、すなわちとについて話します。。
高校でなんかやったなって人もいるかもしれませんが、もっとわかりやすく話しますね。当然この関係は物理をやる時にもかなり大事になります。
では、上の図で考えていきます。今回は少し単純化してお話していくのでご了承ください。
まず、ある関数において区間]における面積を求めよう。
こんなことを言われれば次の積分によって面積を求めると思います。
ここまでは高校の最初の積分の授業で行うと思います。では、図の右側のように長方形を区切ったようなモデルで考えてみましょうか。この長方形の長さは長方形の左下の座標のときの座標とします。
この時に一つ一つの長方形の幅をとします。この時に長方形の左下の座標がにある一つの長方形の面積は次のようになることは簡単にわかります。
ここまでついてこれていますか?わからなければコメントくださいね!
では先に進めます。
では、長方形の数が個あるときのを考えましょう。今回は単純化のためはすべて一定の大きさとします。
ここでは何を表すかといういうと、を表し、一つ一つの幅という意味で追加しています。例えば、仮に今回の長方形2つを1つの長方形とするならとなります。
そうすると、今回の面積を和分(すなわち和)を求めると次のようになります。
また、は次のようになります。ただし、今回、としています。
よって、面積は次のようになります。
高校数学では、この和分のときのを実は省略しています。理由は単純でというように差がである整数のみを扱うからです。すなわち、だからです。
そういうことで、高校数学風に書き直すと次のようになります。
これは余談ですので、元の標識でもう少し進めます。では、最後にこの長方形の幅が無限に小さくなった場合、すなわち飛び飛びでなく連続的な長方形の和を考えます。実はこれが本当の積分というものです。幅が無限に小さいのでとすれば良いので次の関係式が成り立ちます。
ちなみにここまでは普通に理解できるじゃんって思う方も多いかもしれません。ただ、この表式を次のように書かれると一気に謎になる人が多いので最後に書いておきます。
こう見ると積分と和分って同じ式のように見えてきませんか?見えてきたら今回の目的は達成です!!
ちなみに、とは微小な差を表すので次のようになっています。
だから微小面積も次のようになります。
最後にまとめてみると、はの連続版ということです。これって結局、連続と離散の関係になっていますね。だから積分と和分は連続と離散の関係と覚えておきましょう。そうすると実は微分と差分も一緒じゃんとわかりますよ。近いうちに話しますね!では!
2018-12-29
整数の直交
こんにちは
今日は簡単な記事書きますね。
下のものを見て何かわかりますか?
いやいやいや、1188と4165は直交(図形なら垂直と言います)じゃないよね?そんな感想だったら嬉しいです。笑
では、この1188と4165が直交ってどういうこと?
今日は
タイトルにあるように互いに素の表示についてお話しします。
「量子力学〜状態のブラケット表示①〜_の中でもベクトルに触れましたが、ベクトルだと直交していることが簡単にわかりますね。すなわち内積が0になればベクトルは直交していると言えます。一方、今回の1188と4165はなぜ?という疑問が生まれます。
では、数を扱うのでよくやることとして素因数分解をしてみましょう。
となっています。これだけではまだわかりにくいのある整数を素因数分解したら次のようになることを考えましょう。
このときにを次のようにベクトル表示しましょう。
とすると、1188や4165はどうなるのかというと次のようになります。
この2つのベクトルの内積を考えると、
となります。すなわち、この2つの整数1188と4165は直交しているということがわかります。また、このように2つの数の最小公約数が1であるときにその2数は「互いに素」であると言います。
つまり、互いに素である2数は直交していて、と書けるのです。
以上です。面白いので友達に自慢してみてください。
2018-12-29
量子力学〜状態のブラケット表示①〜
こんにちは
とうとう量子力学での記事も3つ目となりますね。
ここまでどうだったでしょうか。ぜひ参考にしてくださいね!
前回までは、粒子の状態は関数で表してきました。おそらく大学の最初の頃は関数の方がわかりやすいと思います。しかし、関数以外の状態の表し方が実は存在します。それがディラック(P. A. M. Dirac)により導入されたブラケット表示というベクトル表示です。
ブラケット(bracket)という単語を辞書で調べて見るとわかりますが、角括弧(<>)という意味です。
なので、状態をのように表示します。また、この表示での名前を紹介すると、をケットベクトル、をブラベクトルと言います。覚えやすさを考えれば、<>がブラケットと英語でいうので、"<"がブラ、">"がケットという風に覚えてしまうのが簡単でしょう。
では、関数で親しんできた状態の表示がなぜベクトルで表示できるの?と疑問に思った方もいると思います。
その疑問を解消するためにもここで内積というものについて振り返ってみたいです。
内積
ベクトルの内積についてちょっと振り返ってみてください。
の2つのベクトルがあったときにその内積は
となります。
正直これが何を表すのかわかっていない方もチラホラいたりするので、もう少し意味について考えますね。
結論から言ってしまうと内積は正射影をしていることと同じです。これじゃわからないのでもう少し砕くとあるベクトル(図の棒)に光を当ててそのベクトルの影を映し出すことに近いです。 一番わかりやすい例は仕事だと思います。
通常、物理で仕事は力×距離で習うと思います。上の図だとだということですね。一方で、下の図だととおそらく多くの高校生は習います。これが意味するのは力ベクトル*1の方向成分と距離をかけたものです。もう少しわかりやすくするために単位ベクトルというものを導入します。単位ベクトルとは長さが1のベクトルです。なぜ単位というのかは、例えば距離の単位はと言われれば1 m、質量の単位なら1 kgと皆さん答えると思います。よっぽどのひねくれ者でなければ、2 mや8 kgを単位としている人はいませんよね。では、距離方向の単位ベクトルは距離で割り算すればよいのでとなります。よって、この単位ベクトルと力の内積は次のようになります。
すなわち、力の距離方向成分が出たというわけです。イメージつかめたでしょうか?内積というのは結局ある方向成分への正射影であるということです。仕事はただここで出たと距離を掛け合わせれば良いので同じ結果が得られます。(当然ですが笑)
では、話を戻しましょう。
なんとなく内積は射影なのかとお分かり頂けたと思います。では、ブラケットでは何をしたいかを話します。という表示があったときにこれはというブラベクトルとというケットベクトルの内積を表します。少しずつ気づいた方も出てきたかと思いますが、状態ベクトルを実空間の位置に正射影するとそれは状態関数になるということなのです。
余談ですが、僕たちが見える状態というのはほとんど場合、実空間への射影だから見えるのですね。もしかしたら別の側面が別の空間には見えているかもしれません。
今日はここまでにしますが、話をまとめると粒子の状態は関数()だけでなくベクトルでも表示できるということ、内積とは正射影であるということ、状態ベクトルの実空間()への正射影が状態関数であるということです。
では、また!
2017-12-28
*1:大学ではベクトルをボールド(太字)で書くことが多い
量子力学〜井戸型ポテンシャル中の粒子〜
こんにちは
前回は「量子力学〜波動関数と波動方程式 - 導入〜」ということで波動関数と波動方程式について紹介しました。
今回は大学の量子力学の授業でおそらくほぼ絶対に触れる井戸型ポテンシャル中の電子についてのお話をできればなと思います!
井戸型ポテンシャルと聞くと少し苦手意識を持つかもしれませんが、簡単に言ってしまえば無限に高い壁に囲まれた空間にいる人を考えていることとあまり変わりません。また、一番想像しやすいのは地面だと思います。地面はポテンシャルが無限大の空間であり、自分たちは地面より下に簡単にいくことができません。この観点をなんとなく持つと理解が深まると思います。
では、パターンに分けて紹介します。
1. 一次元の井戸型ポテンシャル中の粒子
では、まず上のような井戸型ポテンシャル中の粒子を考えます。
ポテンシャルとは上でも言いましたが壁です。粒子は壁に囲まれていて外には出られない状態にあるということが図に表されています。
では、この時の粒子の状態を考えていきます。今回は定常状態(時間に依存しない)を考え、且つ1次元を考えているのでとしています。
その準備としてポテンシャルを数式で表すと以下のようになります。
このポテンシャルの場でSchrödinger方程式を考えてみましょう。
ただし、である。この式を丁寧に展開して、状態をを求めていきますね。
を用いると次のように展開できる。
ただし、である。
(i)の範囲
の範囲を考える。
と置くと次のようになる。
これは単振動の運動方程式と同じ形をしているので、一般解は次のようになる。
は係数である。
この係数を考えるために、一旦の範囲の計算はここまでにして次の範囲を先に考えていこう。
(ii)の範囲
この範囲ではまずとして議論を進めていく。
Schrödinger方程式は以下のようになる。
と置くと次のようになる。今ここででなく、と書いたのはが今回十分に大きい場合を考えているためである。また、を実数にするためでもある。すると方程式は次のように変形されていく。
この方程式の解は次のようになる。
は係数である。
この解を考える時に想像してほしいことは、高いポテンシャルの中に入り込んだとしても無限に遠くなると粒子は存在確率が0になるということである。すなわち、となるはずである。よって、状態は次のようになる。
また、今回はのときにとなるので、状態は改めて書くと以下のようになる。
すなわち、無限大の障壁の中では粒子が存在しないということを表す。
では、(i)の範囲の続きに戻ろう。改めてをここに記す。
関数のでの境界条件を考えると次のようになる。
あと少しなので頑張ろう。さらに上2式をBについてまとめると以下のようになる。
すなわち、となる。実部と虚部に分けるとこの解がどうなるかがわかりやすいだろう。
ということで、関数は次のようになる。
前回「量子力学〜波動関数と波動方程式 - 導入〜」でもお話しましたが、粒子がとる状態ごとの確率はで表せます。今回で言えば、1次元の線上の場所ごとの状態の確率を示します。例えば、なら存在がOO%といった具合です。また、その確率を全体で積分したもの(全体で足し合わせたもの)は確率なので1にいなります。ということで、今回の状態関数についても2乗して積分してその係数を求めてみましょう。
が奇数(odd)のときを考える。
が偶数(even)のときを考える。
ゆえに、(i)の範囲で状態は次のようになる。
これを実際に図示してほしい。どうなるだろう。見た目としては、境界で節になるような波であることがわかると思う。この形を覚えておくことでテストなどでもイメージがしやすくなるだろう。また、この境界に定在波の節がくることを波のピン留めという。また、エネルギーも簡単なので計算してみると次のように量子化していることがわかる。
一方テストなどではポテンシャルの範囲が変わってくることがちょくちょくある。もしもが次のようになっていたら、平行移動を考えてほしい。
すなわち、状態関数を状態関数に変えてみると答えが出るだろう。
2. 三次元の井戸型ポテンシャル中の粒子
これはイメージと結果だけ書きますね。
イメージとしては長さが、幅が、高さがの直方体の空間を想像してほしい。その空間から粒子が抜けられたないとすると粒子は直方体中に収まりますね。やることは実は一次元の井戸と変わりません。ただ、イメージをして欲しかったというだけです。またのそれぞれの方向が独立になっていることを考えましょう。すなわち、がに依存したりと互いに依存することはないものとします。その時、モデルは一次元量子井戸が3方向で存在すると思えば問題ないです!ということは、状態関数に対するシュレディンガー方程式は次のようになります。
このときにが独立であることを利用するとのように変数分離することができる。よって、式を変形していこう。
はポテンシャルのそれぞれの方向の和、はエネルギーのそれぞれの方向の和である。また、今回は方向のポテンシャルはには依存しないものとする。同様にもそれぞれやに依存しないものとする。エネルギーも同様に考えるが、座標には依存しないものとする。ただし、方向のエネルギーをそれぞれとする。
この式は恒等式であり、は独立であったのでそれぞれの方向で常に0にならなければならない。ゆえに3つの方程式が得られる。
少し変形する。
すなわち、1次元井戸型ポテンシャルを全方向でとけばいいだけということがわかる。
今日はここまでしよう。まぁテストなどの参考にぜひ!
2017-12-28
量子力学〜波動関数と波動方程式 - 導入〜
こんにちは
たまには真面目に物理の勉強の話をしようと思います! 量子力学について本日はお話しますね。 極めて簡単なところから話すのでつまらない人もいるかもしれませんがよろしくお願いします。
正直、今回ははてなブログでTeXを使いたいというのがモチベーションですが。
まず量子力学の勉強を始めると最初に出会うのが以下のようなシュレディンガー方程式(Schrödinger方程式)という方程式だと思います。
正直これは何を言ってるかわかりませんね。そもそもこの方程式は何を表すのでしょうか。
結論から言ってしまうとこの方程式は粒子の状態を示す方程式です。すなわち、この方程式を解くと粒子の状態がわかるということになります。また、状態の2乗であるは存在確率を示します。
では、状態と言われてもピンとこないと思います。なので、人で考えるとわかりやすいと思います。
例えば、Aさんが東京駅(位置:)に9:00:00(時間:)にいるときの状態を表すのがAさんの状態関数となり、Aさんはの確率でその状態になっていことがわかります。
つまり、粒子に話を戻すと、粒子がいつどこにいてどんな状態かをを表すのがこの状態関数であり、その状態を記述する方程式がシュレディンガー方程式ということです。また、その状態関数の2乗は粒子がその状態である確率を示すということです。
また、状態関数のことを波動関数、シュレディンガー方程式のことを波動方程式と言います。
わかりましたか?シュレディンガー方程式を解くということは粒子がどの時間にどこにいるかを知ることと同じことです。
次に上に書いたシュレディンガー方程式でおそらく気になってしまったのは、やの上にあるハットは何を示すかということだろうと思います。
これは量子力学で大切になる演算子と言われるものです。演算子とは演算とあるように、ある状態へなにかしらの演算させるという意味で演算子なのです。特にOO演算子はある状態を演算することでOOという物理量を算出するためにあります。ある状態(状態関数)の運動量が知りたいなら運動量演算子を、位置が知りたいのなら位置演算子を、エネルギーが知りたいのならエネルギー演算子を状態関数に作用させればいいわけです。
余談ですが、さっきのAさんの状態で行くとAさんの体温を知りたいので体温演算子というものをとすると、を計算することでAさんの体温がわかるということです。
では、シュレディンガー方程式に話を戻したいので、もう一度先の方程式を書きますね。
この方程式で、左辺の状態関数の前の部分を特にハミルトニアン(Hamiltonian)といいます。今回はそれぞれの演算子の具体的な形について最後に触れて終わりにしますね。詳細は今後別の機会に書くかもしれません。
正直まずはこの形は覚えてしまいたいところです。最後によく見る形にシュレディンガー方程式を書き直します。
今日は波動関数、波動方程式、演算子などについて触れてきましたがそれぞれが何かわかりましたか?
まとめると、波動関数は粒子の状態を表す関数、波動方程式はその状態を記述する方程式、演算子は状態からある物理量を取り出すために波動関数に作用するものです。
うだうだ書いてきましたが結局はTeXの練習をしたかっただけです。笑
2017-12-23
RaspberryPi3の初期設定の仕方
こんにちは
pythonを始めておおよそ1週間が経ち、Lineに仮想通貨の値段が通知されるプログラムを作成しました。
ただ、次にしたいと思ったのが、この通知を自分のノートPCが付いていない時にも実行すること。
ということでRaberry Pi 3 Model Bを買ってみました。
そもそもRaberry Pi 3 Model Bとは・・・?
「Wikipedia : Raspberry Pi参照」
要するに教育用で作られた小さなコンピュータである。
実際に中にはOpenOfficeやMathematicaなど理系の心を惹くようなアプリ、さらにはPython、Javaなどの言語ができます!
調べるともっといろんなことができそうですね。(以下のサイト参照)
Raspberry Pi 3にOpenCV 3.1のインストール | TomoSoft
ではでは実際に行ったことを忘備録として。
準備
付属品も含め、自分は以下の商品を買っています。
・Raspberry Pi3 Model B ボード&ケースセット
・電源(2.4 A)
・USBキーボード、USBマウス
・micro SDカード(class4を買ったけど、classは10の方がおすすめ)
初期設定
では、買ってみたものを使ってRaspberry Pi3を初期設定していこう。
下を参照して進めた。ただ、RaspberryPiのバージョンが変わっているので書き直します。
作業1
SD Card Fomatterをインストールし、SDカードをフォーマット(
FAT)!
作業2
Rasberry PiのホームページでNOOBSのzipファイルをダウンロード
zipファイル解凍し、SDカードのディレクトリ(SDカードのフォルダと考えてOK)の直下に解凍後のフォルダ群をコピーする!
作業3
今フォルダ群を保存したSDカードをRasberry Pi本体に挿れ、さらに本体にUSBマウス、USBキーボード、HDMIケーブルを挿れる。
HDMIケーブルはPCやTVのディスプレイに挿入できます。
その後、電源を本体に挿すと、本体の電源が入る。
このときの動作は
・Rasberry Pi本体の電源あたりが赤く光っており、且つ緑の点滅がある
・ディスプレイが虹色になり、その後OSの選択画面が出てくるので、Raspbian(RECOMMENDED)を選択し、左上のInstallをクリックする。Warningが出てくるが「yes」で大丈夫。
・インストールが終わったら以下の動作へ。
作業4
あとは下のサイトを参照して、wifiやsshの設定を行って完了!
https://qiita.com/susu37/items/142f92b6a8dce6e38737
このやり方でできなかったら、気兼ねなくコメントしてください!
2017-12-23