質量行列

印刷用ページはサポート対象外です。表示エラーが発生する可能性があります。ブラウザーのブックマークを更新し、印刷にはブラウザーの印刷機能を使用してください。

質量行列（しつりょうぎょうれつ）は解析力学では、対称行列 M であり、ある系の一般化座標系qの時間微分 ${\dot {q}}$ とその系の運動エネルギー T との関係を次式で表す。

T={\frac {1}{2}}\mathbf {\dot {q}} ^{\textsf {T}}\mathbf {M} \mathbf {\dot {q}}

ここで、 $\mathbf {\dot {q}} ^{\textsf {T}}$ は、ベクトル $\mathbf {\dot {q}}$ の転置を表す。^[1]この方程式は、質量 $m$

$m$ と速度vを持つ粒子の運動エネルギーの公式に似ている。すなわち

T={\frac {1}{2}}m|\mathbf {v} |^{2}={\frac {1}{2}}\mathbf {v} \cdot m\mathbf {v}

これは、システムの各粒子の位置をqで表すことによって得られる。

一般に、質量行列Mは状態qに依存し、したがって時間とともに変化する。

ラグランジュ力学は、常微分方程式(実際には連立微分方程式のシステム)を生成する。これは、システム内のすべての粒子の位置を完全に定義する一般化座標の任意のベクトルによってシステムの進展を記述する。上記の運動エネルギー式は、すべての粒子の全運動エネルギーを表す方程式の1つの項である。

例

1つの空間次元における質量のシステム。

q={\begin{bmatrix}x_{1}&x_{2}\end{bmatrix}}^{\textsf {T}}

。

T=\sum _{i=1}^{2}{\frac {1}{2}}m_{i}{\dot {x_{i}}}^{2}

T={\frac {1}{2}}{\dot {q}}^{\textsf {T}}M{\dot {q}}

M={\begin{bmatrix}m_{1}&0\\0&m_{2}\end{bmatrix}}

M=\operatorname {diag} \left[m_{1}I_{n_{1}},\,m_{2}I_{n_{2}},\,\ldots ,\,m_{N}I_{n_{N}}\right]

_niが_I×N '_I nは単位行列、または：

M={\begin{bmatrix}m_{1}&\cdots &0&0&\cdots &0&\cdots &0&\cdots &0\\\vdots &\ddots &\vdots &\vdots &\ddots &\vdots &\ddots &\vdots &\ddots &\vdots \\0&\cdots &m_{1}&0&\cdots &0&\cdots &0&\cdots &0\\0&\cdots &0&m_{2}&\cdots &0&\cdots &0&\cdots &0\\\vdots &\ddots &\vdots &\vdots &\ddots &\vdots &\ddots &\vdots &\ddots &\vdots \\0&\cdots &0&0&\cdots &m_{2}&\cdots &0&\cdots &0\\\vdots &\ddots &\vdots &\vdots &\ddots &\vdots &\ddots &\vdots &\ddots &\vdots \\0&\cdots &0&0&\cdots &0&\cdots &m_{N}&\cdots &0\\\vdots &\ddots &\vdots &\vdots &\ddots &\vdots &\ddots &\vdots &\ddots &\vdots \\0&\cdots &0&0&\cdots &0&\cdots &0&\cdots &m_{N}\\\end{bmatrix}}

q={\begin{bmatrix}x&y&\alpha \end{bmatrix}}

{\begin{aligned}x_{1}&=(x,y)+R(\cos \alpha ,\sin \alpha )&v_{1}&=\left({\dot {x}},{\dot {y}}\right)+R{\dot {\alpha }}(-\sin \alpha ,\cos \alpha )\\x_{2}&=(x,y)-R(\cos \alpha ,\sin \alpha )&v_{2}&=\left({\dot {x}},{\dot {y}}\right)-R{\dot {\alpha }}(-\sin \alpha ,\cos \alpha )\end{aligned}}

2T=m{\dot {x}}^{2}+m{\dot {y}}^{2}+mR^{2}{\dot {\alpha }}^{2}-2Rd\sin(\alpha ){\dot {x}}{\dot {\alpha }}+2Rd\cos(\alpha ){\dot {y}}{\dot {\alpha }}

$m=m_{1}+m_{2}$ と $d=m_{1}-m_{2}$ 。この式は、次のように行列形式で記述でき

T={\frac {1}{2}}{\dot {q}}^{\textsf {T}}M{\dot {q}}

M={\begin{bmatrix}m&0&-Rd\sin \alpha \\0&m&Rd\cos \alpha \\-Rd\sin \alpha &Rd\cos \alpha &R^{2}m\end{bmatrix}}

^ Mathematical methods for physics and engineering, K.F. Riley, M.P. Hobson, S.J. Bence, Cambridge University Press, 2010, ISBN 978-0-521-86153-3