位置・速度・加速度と積分

位置 r(t) と速度 v(t) に時刻 t=t0 における初期条件 r(t0) , v(t0) を与えたとき, 加速度 a(t) を積分することで速度を, 速度を積分することで位置を求めることができる. (1)v=v(t0)+t0tadt,r=r(t0)+t0tvdt

位置 r ・速度 v ・加速度 a はそれぞれ, 時間微分をもちいて次のように書けるのであった.
位置・速度・加速度と微分 (2)v=drdta=dvdt=d2rdt2 すなわち, 位置 r(t) がわかっているならばその導関数 drdt を求めることで速度 v が定まり, 速度 v(t) がわかっているならばその導関数 dvdt を求めることで加速度 a が定まることを意味していた.

今回は逆に, 加速度がわかっているときに速度を, 速度がわかっているときに位置を求めるような手法について議論しよう. ここで大きなヒントとなるのは, va はそれぞれ rv を時間微分することで求められること, 数学的には微分と積分が互いに逆演算のようにみなせること, である.

下準備として積分法の解説を行うが, いわゆる区分求積法や置換積分などを知っている人は読み飛ばしてもらってかまわない.

下準備 : 積分法

詳しくは積分法も参照していただきたいが, 積分自体の説明と積分が微分と互いに逆演算のような関係にあることについて議論する.

定積分

いま, x を変数とした関数 f(x) が下図のように与えられているとする. この関数 f について, x=ax=b (>a) で区切られた区間 axbn 等分し, xk を次のように定義する. {x0=ax1=a+1banx2=a+2ban=xk=a+kban=xn1=a+(n1)banxn=b. また, (xkxk1)=ban という区間の幅を Δx と書くことにする.

このとき, Δx=banf(xk) の積 Sk=f(xk)Δx(k=1,2,,n) は底辺が Δx=ban で高さが f(xk) の長方形の符号付き面積に等しい[1]ここで符号付き面積といったのは, f(xk) が負の場合には Sk も負であり, その絶対値は図に示した長方形の面積に等しいからである. … Continue reading. そして, このような長方形を axb の全区間に渡って足し合わせた (3)S1+S2++Sn=f(x1)Δx+f(x2)Δx++f(xn)Δx=k=1nf(xk)Δx は上図に示したように, x=ax=b および f(x)x 軸とで囲まれた符号付き面積に近しい値となっている.

したがって, 分割数 n が非常に大きいという極限 n を適用したもの (4)limnk=1nf(xk)Δxaxb において, f(x)x 軸との間で囲まれる(符号付き)面積に一致することになる. 式(4)で定義されるような計算を定積分といい, 記号 (インテグラル)を用いて abf(x)dx:=limnk=1nf(xk)Δx=limnbank=1nf(xk) と表す.

定積分の性質

定積分は次の性質 (5)ab{kf(x)+lg(x)}dx=kabf(x)dx+labg(x)dxabf(x)dx=acf(x)dx+cbf(x)dx を満たす. これらの性質は定積分が面積に相当していることを鑑みると実に納得の行くものとなっている. さらに, 定積分の始点と終点を入れ換えると符号が反転し, baf(x)dx=abf(x)dx が成立する. この性質も上述の議論を追えば納得できるものであろう.

また, 定積分 abf(x)dx はある定数値であり, この式における x という文字は積分の計算結果にはあらわれてこない. したがって, x を別の文字に置き換えたとしても意味は変わらず, abf(x)dx=abf(x)dx=abf(t)dt=abf()d などが成立することには注意してほしい.

微分と積分の関係

関数の定積分を計算するにあたり, 毎回上記のような考察を行う必要はなく, 積分演算は微分の逆演算に相当していることを示しておこう.

いま, 変数 x の関数 f(x) を考え, 0xx の区間で f(x)x 軸とで囲まれた領域の符号付き面積は区間の終点の x 座標 x の関数であり, これを F(x) と定義しよう. F(x):=0xf(x)dx. 続いて, 定積分の範囲を微小な量 Δx だけ伸ばした F(x+Δx)F(x+Δx)=0x+Δxf(x)dx と計算することができる.

上図より, F(x+Δx)F(x) との差は近似的に f(x)Δx とみなすことができるので, (6)F(x+Δx)F(x)f(x)Δx f(x)F(x+Δx)F(x)Δx と書くことができる. ここで ほぼ等しいことを意味する記号である.

上式の右辺に対して Δx を限りなく小さくする極限操作 Δx0 を適用した limΔx0F(x+Δx)F(x)Δx は関数 F(x)導関数の定義そのものであり, 関数 fF の間には (7)dF(x)dx=f(x) が成立していることになる. つまり, 関数 F(x)f(x) を導関数に持つような関数なのである.

したがって, 関数 f を積分することで定義された F(x) について, F(x) を微分することで f(x) が求められることになり, 微分と積分は互いに逆演算のような関係にあることがわかる.

また, これで定積分を計算する術が得られたことになる. 実際, (8)abf(x)dx=a0f(x)dx+0bf(x)dx=0bf(x)dx0af(x)dx=F(b)F(a) と書くことができる. これが意味していることは, 関数 f(x)x=a から x=b までの定積分を計算したければ, f を導関数に持つような関数 F(x) を見つけて F(b)F(a) を計算せよ ということである[2] F のことを f原始関数などという. 原始関数については後述する..

この定積分計算を次のように簡略化して書くことにする. abf(x)dx=[F(x)]ab=F(b)F(a).

微分積分学の基本定理

(7)に示した ddx0xf(x)dx の導出過程を見返してみると, 定積分の始点は x=0 でない別の値 x=a を定積分の始点に選んでも ddxaxf(x)dx=ddxa0f(x)dx=const.+ddx0xf(x)dx=f(x) が成立する. ここで, 定積分 a0f(x)dx はとある定数値となるので, 微分するとゼロとなることをもちいた.

次の関係式 ddxaxf(x)dx=f(x)微分積分学の基本定理と呼ばれる.

原始関数と不定積分

ここでは, 式 (9)dF(x)dx=f(x) を満たすような F(x) は一つのみに定まるわけではないことを示しておこう.

まず, 式(9)を満たすような関数として F1 が見出されたとしよう. つまり, F1(x)dF1dx=f を満たすものである. この関数 F1 に定数 C を足した関数 F2=F1(x)+C も微分すると, dF2dx=f を満たすので, F2 も式(9)を満たす関数なのである.

以上より, 式(9)を満たすある関数に任意の定数 C を足した関数は全て式(9)を満たすことがわかる.

関数 f が与えられたとき, dFdx=f となるような関数 Ff原始関数と呼ぶ. また, f の原始関数を記号 f(x)dx であらわし, f不定積分という.

先に示したとおり, f の原始関数 F に定数 C を加えた関数 F+Cf の原始関数であり, f(x)dx=F(x)+C である. この任意の定数 C のことを積分定数という.

速度の時間積分

速度(1次元)の時間積分

1次元運動を行う点 P の位置を x(t) とすると, 速度 v(t) は位置の導関数 v(t)=dx(t)dt で与えられるのであった.

では, 速度 v が与えられたときに位置に関する情報を得るためにはどうするかといえば, 微分と積分が互いに逆演算のような関係であることを利用するのである. すなわち, 位置の時間微分が速度なので, 速度を時間で積分することで位置に関する情報が得られることを利用する.

まずは, 速度 v(t) がわかっているときに時刻 t=t1 から t=t2 の間の点 P の変位 x(t2)x(t1) を求めてみよう.

速度 v=dxdt の時刻 t1 から t2 の間の定積分は, t1t2v(t)dt=[x(t)]t1t2=x(t2)x(t1) である. この右辺は時刻 t1 から t2 の間の変位に他ならず, ある時間における変位は速度の(定)積分で計算可能であることが示された.

次に, 時刻 t0 における位置が x(t0) であることがわかっているとしよう. このとき, 時刻 t における位置 x(t) は速度の積分をもちいて次のように書くことができる. (10)x(t)=x(t)+x(t0)x(t0)=x(t0)+{x(t)x(t0)}=x(t0)+t0tv(t)dt. ここで, 右辺第1項は初期位置であり, 右辺第2項の積分計算は時刻 t0 から t までの間の変位である.

したがって, 位置 x(t0) と速度 v(t) が与えられたならば, あとは v の積分計算を行うことによって位置の時間変化を完全に記述することが可能となるのである. この x(t0)初期位置という. そして, この初期位置を設定するという行為は速度 v の(不定)積分で登場する任意定数 C を固定することに相当している.

vt グラフ

位置・速度・加速度の関係を, 縦軸を速度, 横軸を時刻とした vt グラフで考えてみよう.

vt グラフ上に与えられた曲線 v から微分・積分を用いて得られる情報は (11){dv(t)dt=a(t)t1t2v(t)dt=x(t2)x(t1) である. 第1式は v のある時刻の傾きが加速度に一致していることを意味し, 第2式はある時間内で v を積分するとその間の変位に一致していることを意味している.

上図には, t1t2v(t)dt が正となるような図を描いており, この場合の変位は正となる. 一方, t1t2v(t)dt が負の場合, すなわち, vt グラフの符号付き面積が負の場合には変位も負であることに注意してほしい.

また, vt グラフの符号付き面積からは変位しか読み取れないので, 初期位置が与えられてはじめて位置 x が完全に定まることも理解できるであろう.

1次元運動の総走行距離(道のり)

これまでに, vt グラフの符号付き面積が変位に一致することを議論した. しかし, 位置の移り変わりを追いかけていくにあたっては, 変位の他に総走行距離または道のりと呼ばれる量を考えることもできる.

変位がある時間幅の始点と終点でどれだけ位置が変わっているかという正味の移動距離を表す量であることに対し, 総走行距離(または道のり)は, 実際に移動した経路の距離を表す量である.

変位と総走行距離に違いが生じるのは, 物体が x 軸の負方向に運動する( v<0 である)瞬間が含まれている場合である. vt グラフで言えば, v<0 で符号付き面積の符号が負となる瞬間を含んでいるかどうかの違いである.

変位は, 正方向へ進むことでプラスの変位, 負方向へ進むことでマイナスの変位となる. そして, 最終的な変位はこれらを足し合わせることで得られるので, 正味の移動距離と方向しかわからないのである.

一方, 正方向に進んでも負方向に進んでも総走行距離は増加していくので, 総走行距離はどちらどれだけに進んだのかはわからず, 実際に移動した経路の長さという情報のみを含んでいる.

時刻 t1 から t2 の間の総走行距離 s は次式で与えられることになる. s=t1t2|v|dx=t1t2|dxdt|dx. 下に示した vt グラフには, ある速度 v(t) の曲線(黒の実線)を描いている. この vt グラフを用いて t1 から t3 の間の変位総走行距離を求める方法について考えよう. ただし, t1 から t2 の間は v>0 , t2 から t3 の間は v<0 とする.

まず, t1 から t3 の間の変位 Δx13 は次式で与えられる. (12)Δx13=t1t3v(t)dt=t1t2v(t)dt+t2t3v(t)dt. また, t1 から t3 の間の総走行距離(道のり) s13 は次式で与えられる. (13)s13=t1t2|v(t)|dt+t2t3|v(t)|dt.=t1t2v(t)dtt2t3v(t)dt. これは上図に示した |v(t)| の曲線(紫の鎖線)が時間軸と成す面積を求めていることに等しい.

速度(3次元)の時間積分

これまでに行なってきた, 速度を時間積分することで変位が得られるという議論は, x , y , z の各方向に対して成立する.

3次元空間内の点 P の座標を位置(ベクトル)を r=(x,y,z) であらわし, 速度(ベクトル) vv=drdt=(vx,vy,vz) であることがわかっているとする. このとき, 時刻 t1 から t2 の間の変位は (14)t1t2vdt=r(t2)r(t1) {t1t2vxdt=x(t2)x(t1)t1t2vydt=y(t2)y(t1)t1t2vzdt=z(t2)z(t1) と計算することができる.

したがって, 時刻 t0 における初期位置r(t0) であるとき, 時刻 t における位置 r(t) を次のように求めることができる. r(t)=r(t0)+t0tv(t)dt. これは, 初期位置 r(t0) と速度 v(t) が与えられたならば, v の時間積分を行うことによって時々刻々と変化する位置を完全に記述することが可能となる.

加速度の時間積分

加速度の時間積分によって速度が与えられることを議論しよう. なお, 以下の議論は速度の時間積分で変位が計算できたことと全く同じ構造となっている.

加速度(1次元)の時間積分

1次元運動を行う点 P の速度を v(t) とすると, 加速度 a(t) は速度の導関数 a(t)=dv(t)dt で与えられるのであった. したがって, 加速度を時間積分することである時間の間の速度の変化量を知ることができる.

加速度 a(t) がわかっているとき, 時刻 t=t1 から t=t2 の間の速度の変化 v(t2)v(t1) は次式のように v=dxdt の時刻 t1 から t2 の間の定積分で与えられる. t1t2a(t)dt=[v(t)]t1t2=v(t2)v(t1). また, 時刻 t0 における速度 v(t0) がわかっているとしよう. このとき, 加速度 a がわかっているならば, 時刻 t における速度 v(t) を次のように求めることができる. (15)v(t)=v(t0)+{v(t)v(t0)}=v(t0)+t0ta(t)dt. したがって, 速度 v(t0) と加速度 a が与えられたならば, あとは a の時間積分を行うことによって速度の時間変化を知ることが可能となる. この v(t0)初速度という. そして, この初速度を設定するという行為は加速度 a の(不定)積分で登場する任意定数 C を固定することに相当している.

at グラフ

速度・加速度の関係を, 縦軸を加速度, 横軸を時刻とした at グラフで考えてみよう.

加速度 a と速度 v について t1t2a(t)dt=v(t2)v(t1) という関係が成立していることから, at グラフ上に与えられた曲線 a(t) と時間軸および, ある時間幅で囲まれた符号付き面積を求めることにより, その区間における速度の正味の変化量を知ることができる.

また, at グラフの符号付き面積からは速度の変化しか読み取れないので, 初速度が与えられてはじめて v(t) が定まることも理解できるであろう.

加速度(3次元)の時間積分

これまでに行なってきた, 加速度を時間積分することで速度の正味の変化量が得られるという議論は, x , y , z の各方向に対して成立する.

3次元空間内の点 P の座標の速度(ベクトル)を v=(vx,vy,vz) であらわし, 加速度 aa=dvdt=(ax,ay,az) であることがわかっているとする. このとき, 時刻 t1 から t2 の間の変位は (16)t1t2adt=v(t2)v(t1) {t1t2axdt=vx(t2)vx(t1)t1t2aydt=vy(t2)vy(t1)t1t2azdt=vz(t2)vz(t1) と計算することができる.

したがって, 時刻 t0 における初速度v(t0) であるとき, 時刻 t における速度 v(t) を次のように求めることができる. v(t)=v(t0)+t0ta(t)dt. このように, 初速度 v(t0) と加速度 a(t) が与えられたならば, a の時間積分を行うことによって時々刻々と変化する速度を知ることができ, その速度を更に時間積分することで物体の位置が定まることになる.

経路に沿った積分

時刻 tA に点 rA=r(tA) いた点 P が3次元空間上を移動して時刻 tB に点 rB=r(tB) へ移動したとしよう. このとき, 時刻 tA から tB の間の変位 (rBrA) は, 時刻 tA から tB の間の速度 v時間積分 tAtBvdt=r(tB)r(tA) で与えられるのであった. このことを, 少し別の形で表現する手法を紹介しておこう.

証明は省略するが, x=x(t) の関数 f(x) の積分において成立する次の置換積分の公式 f(x)dx=f(x(t))dxdtdt が成立することが知られている. この式を速度の積分 tAtBv(t)dt の各成分に適用しよう. x 成分だけを具体的に書き出してみると, (17)tAtBvx(t)dt=tAtBdxdtdt=x(tA)x(tB)dx=x(tB)x(tA) が成立することになる. これは各成分に対して適用可能であるので, (18)tAtBv(t)dt=tAtBdr(t)dtdt=r(tA)r(tB)dr と言った具合にまとめて書くことにしよう. 以下では式(18)の意味についてもう少し物理的に考えてみよう.

P は3次元的な経路を通って移動するので, 今回 rA から rB へ移動する間に点 P が通過する経路を CAB と名付けよう.

この経路 CAB を下図のように n+1 個の点 r0=rA , r1 , , rk , , rn=rBn 個の線分にわけて, 点 rk から点 rk+1 への微小な変位を Δrk:=rkrk1 と定義しよう.

このような定義に従うと, 点 rA から点 rB への変位は微小変位 Δrk の総和 (19)k=1nΔrk=(r1r0)+(r2r1)++(rnrn1)=rnr0=rBrA で計算することができる.

ここで, 分割数 n を非常に大きくとる n という極限を上式に適用した limnk=1nΔrk において, Δrk は経路 CAB 上からわずかも外れない微小変位ベクトルとなり, これを経路 CAB に沿って順次足し合わせることを意味している. これはまさしく式(18)の意味していることに一致しており, rArBdr=limnk=1nΔrk と書くことができる.

以上の議論をまとめると, rA(tA) から rB(tB) への変位を計算する手法として, 物体の速度を時間で積分する手法微小な変位を積分する手法とがあることを紹介した.

物体の速度を時間で積分する手法 (20)r(tB)r(tA)=tAtBvdt=tAtBdrdtdt は時間間隔 Δt を非常に小さくし, その間の正味の移動量 vΔt を積み重ねることを意味している.

微小な変位を足し合わせる手法 (21)r(tB)r(tA)=rArBdr は経路にそった微小な変位 Δr を積み重ねることを意味している.

(20)と式(21)が同じ意味であることは, 置換積分の公式を認めれば明らかである. このようなことをくどくどと説明した理由は, 物理ではこのような議論の行き来が頻繁に行われるからである. これは, 積分時にどのような変数に注目したほうが議論対象の物理を語るのに楽[適切]か, と言うだけの話である. 積分時の変数の違いや細かな表現の違いに物怖じしないよう積分法の理解を深めておきたい.

脚注

脚注
1 ここで符号付き面積といったのは, f(xk) が負の場合には Sk も負であり, その絶対値は図に示した長方形の面積に等しいからである. また, いまは b>a を考えているが, より一般には ab の大小関係によって ban の符号が負になるなどの事情を含んでいる.
2 F のことを f原始関数などという. 原始関数については後述する.