2012年2月8日水曜日

行列の掛け算の定理

 
 
大学への数学Ⅲ&Cの勉強
行列と連立1次方程式

【解説】
行列の掛け算はどういう演算に応用でき、どういう演算には応用できないかは、行列の掛け算の定理により定められます。
以下で、この行列の掛け算の定理を証明します。

 
(条件1)ベクトルDが、以下の式1と式2によりベクトルFに変換されること。その式の係数が行列Aを成すこと。

式1と式2の演算の係数を添え字を付けて区別して式3で表します。
この式3は、式3’とあらわすと便利です。
この表し方はアインシュタインが提案した表記方法です(アインシャタインの縮約記法)。
この式3’は、変数で表した添え字の変数名が同じ場合に、その変数名の添え字のあらゆる場合(この場合は1と2のみ)について和をとることにした式です。

(条件2)ベクトルFが、以下の式4と式5によりベクトルGに変換されること。その式の係数が行列Bを成すこと。

式4と式5の演算の係数を添え字を付けて式6で表します。
この式6は、アインシュタインの縮約記法の式6’であらわせます。

式6’に式3’を代入すると式7が得られます。
式7は、更に式8に変形できます(式8はアインシュタインの縮約記法であらわした式です)。
式8’のように、演算Aと演算Bを合成した演算Cを考えます。
式8から、演算Cをあらわす行列Cの要素の値を計算する式9が得られます。
こうして、演算AとBの行列の要素を用いた式9で演算Cをあらわす行列Cの要素が計算できることがわかりました。

この演算Cを、式10のように、演算Aと演算Bの積(掛け算)と定義します。
そして、その演算Aと演算Bの掛け算である演算Cの行列の要素は式9’で計算できます。
これが、行列の掛け算の定理です。

 すなわち、行列の掛け算の定理を、以下のようにあらわすことができます。
(条件1)ベクトルを式1~2でベクトルに変換する演算Aの行列Aが定義され、
(条件2)ベクトルを式4~5でベクトルに変換する演算Bの行列Bが定義されている場合に以下のことが成り立つ。
(結果)
 第1のベクトルを演算Aで変換して第2のベクトルを得、更に、第2のベクトルを演算Bで変換して第3のベクトルを得る場合、
その2つの演算は、第1のベクトルを第3のベクトルに変換する演算Cにまとめることができる。
演算Cの行列Cは、式9’で計算することで求めることができる。
その演算Cを、演算Aに演算Bを掛け算した演算であると定義する。
演算Cの行列の要素を計算する式9’の計算を、行列の掛け算と定義する。

【計算の形】
 なお、この掛け算の計算方式をベクトルAとBの内積にあてはめると、以下の形の行列の計算になる。

 この計算は、横ベクトルAと縦ベクトルBの積であらわせる。
 一方、アインシュタインの縮約記法では、ベクトルを縦ベクトルと横ベクトルの形に区別しないで内積の計算を表現できる。

 
 

リンク:
追加講:三角形の面積と行列式
高校数学の目次

2012年2月7日火曜日

回転変換の行列の計算で三角関数の加法定理を得る

 

大学への数学Ⅲ&Cの勉強
行列と連立1次方程式

【解説】

(X,Y)座標であらわされるベクトルを、左回りに角度θ回転させて、(X’,Y’)座標であらわされるベクトルにする変換は以下の式になります。

X’=(cosθ)X-(sinθ)Y
Y’=(sinθ)X+(cosθ)Y
この変換は次式のように行列であらわせます。

θ=α+β
の場合、
(1)角度αの回転変換の行列(回転行列)は、行列をベクトルに掛け算してベクトルを得る式であらわされている。
(2)角度βの回転変換の行列(回転行列)も、行列をベクトルに掛け算してベクトルを得る式であらわされている。
(3)そのため、以下の、行列の掛け算の定理が成り立つ。
すなわち、
(A)第1のベクトルに角度αの回転変換をして第2のベクトルを得、
(B)その第2のベクトルに角度βの回転変換を行って第3のベクトルを得る変換は、
(C)1つの行列であらわせる。
その行列は、
角度αの回転行列に角度βの回転行列を行列の掛け算計算で得た行列と等しい。

また、第1のベクトルに角度αの回転変換をして得た第2のベクトルに角度βの回転変換を行って第3のベクトルを得る変換は、θ=α+βの回転変換と同じである。

その回転行列を計算して比較すると、上式のように、三角関数の加法定理が得られます。

(補足)
なお、
(1)角度αの回転行列と角度βの回転行列を掛け合わせた計算計算結果に加法定理を適用すると角度(α+β)の回転行列になるので、
(2)角度αの回転に角度βの回転を重ねると角度(α+β)の回転になる
と説明されることがあるが、
その説明は本末転倒です。

角度αの回転に角度βの回転を重ねると角度(α+β)の回転になるのは自明のことであって、
(1)その自明な事実に、
(2)2つの演算を重ね合わせた演算の行列は行列の掛け算によって導出できる行列の掛け算の定理を適用することで、
(3)加法定理という関係式が導き出されます。

定理とその定理から導出できる結果とは、その関係を逆転させて、結果から大元の定理を導き出すこともできます。
そのため、どの事実が自明で主となる事実であり、どの事実が従属的に分かる関係であるかを的確に判断する数学的センスが大切です。

(蛇足)
ここで、2θの回転行列とその2分の1の回転行列の要素の間には以下の関係があることを覚えておくと便利です。



 

リンク:
追加講:三角形の面積と行列式
高校数学の目次

2012年2月6日月曜日

座標の回転変換の式のおぼえかた


《目次》行列演算からベクトル解析まで

【解説】

(X,Y)の座標であらわすベクトルを、左回りに角度θ回転させて、(X’,Y’)の座標であらわすベクトルにする変換は以下の式になります。

X’=(cosθ)X-(sinθ)Y
Y’=(sinθ)X+(cosθ)Y

この変換の式のおぼえ方は、上の図を思い描いて、図から、

X座標成分のみがあるベクトルを回転させたベクトルのY’座標が(sinθ)・Xであること。
Y座標のみがあるベクトルを回転させたベクトルのX’座標が(-sinθ)・Yであること。

を想像する。
次に、それ以外は、係数がcosθの項を加えることで、回転変換の式ができあがります。

(補足)
三角関数の加法定理を知らないと回転行列が導き出せないという誤解が流通しているようですが、そのようなことはありません。
X座標成分のみがあるベクトルを回転させたベクトルのY’座標が(sinθ)・Xである。
X座標成分のみがあるベクトルを回転させたベクトルのX’座標が(cosθ)・Xである。

そのため、回転行列Aの要素は、
A11=(cosθ)
A21=(sinθ)
Y座標成分のみがあるベクトルを回転させたベクトルのX’座標が(-sinθ)・Yである。
Y座標成分のみがあるベクトルを回転させたベクトルのY’座標が(cosθ)・Yである。

そのため、回転行列Aの要素は、
A12=(-sinθ)
A22=(cosθ)
こうして、加法定理を使わずに、回転行列の全ての要素がわかります。


この回転行列が全てのベクトルを左回りに角度θ回転することの証明:
(1)全てのベクトルはX方向のベクトルとY方向のベクトルに分解されます。
(2)そして、その各々のベクトルが同じ角度θ回転します。
(3)そのため、その成分を合成したベクトルもまた角度θ回転します。 
(証明おわり)



リンク:
複素数平面での座標回転を応用した例
追加講:三角形の面積と行列式
高校数学の目次

行列演算からベクトル解析まで目次

2012年2月5日日曜日

行列の固有値が重根の場合の固有ベクトル




大学への数学Ⅲ&Cの勉強
行列と連立1次方程式

【解説】
2行2列の行列の場合は、固有値λは、以下の行列式を0にするλを計算することで求められます。
εsp(As1-λEs1)(Ap2-λEp2)=0
ここで、行列Espは単位行列です。
また、そのλの解(固有値)α、β毎に、固有ベクトルがあります。

すなわち、行列は、固有値毎のあるベクトル(固有ベクトル)を固有値倍にします。

一方λの解が重根である場合、すなわち、α=βの場合は、単位行列に比例する行列でなければ、
固有値αに対応する固有ベクトルは1つしか無い。
(固有ベクトルが2つ無ければ、その行列は対角化できない異常な行列です。その行列は、対角化できないが、ここをクリックした先のサイトのように三角化はできる。)

(注意)単位行列に比例する行列や0行列の場合は、例外的に、全方向のベクトルが固有ベクトルであり、固有値の重複した数と同じ2つの固有ベクトルを持つことに注意。
3行3列の行列において、固有値が重複し、かつ、その重複した数と同じ数の固有ベクトルを持つ行列の例は、ここをクリック先のページの行列Bを参照。


《第1の例》
考え易くするために、λ=α=β=0となる行列Aspを考える。
そのような場合の行列Aspの例を下図に示す。

εsp(As1-λEs1)(Ap2-λEp2)
=(A11-λ)(A22-λ)-A12A21
=(1-λ)(-1-λ)+1=λ2=0

この行列の余因子行列Cptは、上図のように求められる。
よって、この行列Aspは、余因子行列の列ベクトル

をλ=0倍の0ベクトルに変換する。この列ベクトル

が固有ベクトルである。


一方、行列Aspは、ベクトル

を、

に変換する。

すなわち、固有値0に対する固有ベクトルに平行な、固有ベクトルの(-a-b)倍の、ベクトルに変換する。
すなわち、列ベクトル(a,b)に対して、列ベクトル(-1,1)に垂直な列ベクトル(-1,-1)との内積(-a-b)を計算して、その値を列ベクトル(-1,1)に掛け算した列ベクトルを求める、という変換を行なうのです。

行列Aspが任意の列ベクトル

を固有ベクトルに平行なベクトルに変換した後で再度行列Aspで変換すれば、それは0ベクトルになる。


同様に、λ=α=β=0となる2行2列の行列Aspは、どの行列も、
その行列で任意のベクトルを変換すれば、固有値が0の固有ベクトルに平行なベクトルに変換し、
2重に変換すれば、任意のベクトルを0ベクトルに変換する。
任意のベクトルを0ベクトルに変換する元の行列は0行列である。すなわち、
AspApt=0st
である。

【第1の証明】
このことの証明を簡単に行なうには、ケイリー・ハミルトンの定理を先に証明しておいて、ケイリー・ハミルトンの定理を使って、これを証明するのが良いように考える。

(証明開始)
2行2列のケイリー・ハミルトンの定理は、固有値αとβに関して、
(Asp-αEsp)(Apt-βEpt)=0st
である。
このα=β=0の場合は、
AspApt=0st
になる。
また、Aspの余因子行列Cmtの列ベクトルをCtとすると、それは、固有値α=β=0の固有ベクトルであるから、
AspCp=0s
また、任意のベクトルDtに関して、
AspAptDt=0s
∴
AptDtは、行列Aspで0ベクトルに変換される固有ベクトルCpに平行なベクトルである。
(証明おわり)

【第2の証明】
(証明開始)
重根の場合、すなわち、α=β=0となる場合に、
1つの固有ベクトル
(C11, C21)=(-1,1)≡Cp 
が得られ、その固有ベクトルに対して、
(Amp-α・Emp)Cp=0m  (式1)
になる。
次に、行列(Amp-α・Emp)による任意のベクトルDpの変換結果のベクトルは、以下の式2であらわせる。すなわち、 
(1)ベクトルDmがベクトルCmと線形独立な場合は、ベクトルDmと固有ベクトルCmを用いて、 
(Amp-α・Emp)Dp=gDm+hCm  (式2)
とあらわせる。
(2)また、ベクトルDmがベクトルCmと線形独立で無くベクトルCmの定数倍の場合は、g=0、h=0とあらわせて、やはり式2であらわせる。

(仮定)ここで、g≠0と仮定すると、

式1により、
(Amp-α・Emp)(h/g)Cp=0m  (式3) 
が得られる。
この式3を式2に足し合わせると、
(Amp-α・Emp)(Dp+(h/g)Cp)=gDm+hCm ,
Amp(Dp+(h/g)Cp)=(α+g)・(Dm+(h/g)Cm) (式4)
式4は、ベクトル(Dp+(h/g)Cp)が行列Ampにより(α+g)倍に変換されることをあらわしているので、
α+gが固有値であることになる。
一方、固有値は重根であって、αのみであるので、
g=0
になる。
しかし、g≠0と仮定していたので、これと矛盾する。
よって、g≠0とした最初の仮定が不適である。
∴ g=0

任意のベクトルFに関する式2は、
(Amp-α・Emp)Dp=hCm
∴
(Amp-α・Emp)Dtは、行列Aspで0ベクトルに変換される固有ベクトルCpに平行なベクトルである。
(証明おわり)

《第2の例》
固有値が1の重根を持つ行列A:

は、固有ベクトル

を1つのみ持つ。
その固有ベクトルと同じ固有ベクトルを持つ行列B:

は、行列Aとの積が交換可能(可換)な関係があり、以下の式が成り立つ。

固有ベクトルを共有する行列同士の積が交換可能であるという関係は一般的に成り立つ。




リンク:
固有ベクトルが同じ行列
追加講:三角形の面積と行列式
高校数学の目次

ベクトルの座標を回転させる行列の固有値と固有ベクトル

 
 
大学への数学Ⅲ&Cの勉強
行列と連立1次方程式

【解説】
座標値(x,y)であらわされるベクトルを、左回りに角度θ回転させて、座標値(x’,y’)であらわされるベクトルにする変換をあらわすと以下の式になります。
x’=(cosθ)x-(sinθ)y
y’=(sinθ)x+(cosθ)y
この変換は、以下のように行列Aspであらわされています。
X’s=AspXp 
ただし、この表記方法はアインシュタインの行列の積の表記方法です。通常の行列表現は以下のものです。

そして、ベクトル及び行列の各要素は以下の通りです。
x’1=x’
x’2=y’
x1=x
x2=y
A11=cosθ
A12=-sinθ
A21=sinθ
A22=cosθ
この行列Aspは回転変換をあらわしますが、この行列は、あるベクトル(固有ベクトル)をλ倍にする変換でもある。
λは、以下の行列式を0にするλ(それを固有値と呼ぶ)を計算することで求められます。
εsp(As1-λEs1)(Ap2-λEp2)=0
ここで、行列Espは単位行列です。
また、そのλの解(固有値)α、β毎に、
(1)λの解α(=固有値)については、
行列(Asp-αEsp)の余因子行列Cptを計算し、その余因子行列Cptの列だけを取り出した列ベクトルCp1(又はCp2)が、
(Asp-αEsp)Cp1=0s
(0sは0ベクトル)
AspCp1=α・EspCp1=α・Cs1
を満足する。Cp1がα倍になる。
Cp1が固有ベクトルである。
Cp2もα倍になる固有ベクトルであるが、Cp1に平行なベクトルであるので、Cp1を求めるのもCp2を求めるのも同じことである。

(2)λの解β(=固有値)については、
行列(Asp-βEsp)の余因子行列C’ptを計算し、その余因子行列C’ptの列だけを取り出した列ベクトルC’p1(又はC’p2)が、
AspC’p1=β・C’s1
を満足する。
C’p1(又はC’p2)が固有ベクトルである。

(3)固有値を具体的に求める
εsp(As1-λEs1)(Ap2-λEp2)=0
(cosθ-λ)(cosθ-λ)+(sinθ)(sinθ)=0
λ2-2(cosθ)λ+1=0
(λ-cosθ)2=-sin2θ
(λ-cosθ)=±i・sinθ
λ=cosθ±i・sinθ
α=cosθ+i・sinθ
β=cosθ-i・sinθ

(4)この固有値αに対する行列(Asp-αEsp)の余因子行列Cptを計算する。
A11-αE11=-i・sinθ
A12-αE12=-sinθ
A21-αE21=sinθ
A22-αE22=-i・sinθ
余因子行列Cptは、
C11=-i・sinθ
C12=sinθ
C21=-sinθ
C22=-i・sinθ
よって、固有ベクトルは、
(-i・sinθ, -sinθ)
 sinθ≠0のとき、
-sinθで割り算して、もっと簡単にした固有ベクトルは、
(i, 1)
(注意)sinθ=0の場合、この割り算計算は(0で割り算するので)無効です。sinθ=0の場合、すなわち、回転角度が0の場合の行列Aspの固有ベクトルが(i, 1)であるかかどうかが定まりません。実際、回転角度が0の場合の行列Aspは単位行列Espになりますが、その固有ベクトルは、あらゆる方向のベクトル全てが固有ベクトルになります。(注意おわり)


この固有ベクトルに左から行列Aspを掛け算すると、
(i・cosθ-sinθ, i・sinθ+cosθ)
=(i・sinθ+cosθ)(i, 1)=α(i, 1)
 行列Aspは、実数のベクトルだけを考えると回転変換しかしないのに、複素数で考えると、固有ベクトルを固有値倍に変換する。

(5)固有値βに対する行列(Asp-βEsp)の余因子行列Cptを計算する。
A11-βE11=i・sinθ
A12-βE12=-sinθ
A21-βE21=sinθ
A22-βE22=i・sinθ
余因子行列Cptは、
C11=i・sinθ
C12=sinθ
C21=-sinθ
C22=i・sinθ
よって、固有ベクトルは、
(i・sinθ, -sinθ)
sinθ≠0のとき、
sinθで割り算して、もっと簡単にした固有ベクトルは、
(i, -1)
この固有ベクトルに左から行列Aspを掛け算すると、
(i・cosθ+sinθ, i・sinθ-cosθ)
=(-i・sinθ+cosθ)(i, -1)=β(i, -1)
 このように、固有ベクトルを固有値倍に変換する。




リンク:
追加講:三角形の面積と行列式
高校数学の目次

2012年2月3日金曜日

行列式が0になる行列は、余因子行列の列ベクトルを0ベクトルに変換する

 
 
大学への数学Ⅲ&Cの勉強
行列と連立1次方程式

行列Ampに対する余因子行列Bstを、以下の式で定義します。

(1)
2行2列の行列Ampの余因子行列Bstは、
エディントンの行列式の計算記号εmpを使って、アインシュタインの縮約記法であらわして、
B1m=εmpAp2
B2p=εmpAm1
で計算します。
(補足)これらの式は、行列式εmpAm1Ap2の計算式から、余因子行列Bstに対応して、行列Aのs列目の項を除去した式です。
この計算に利用する行列
εmp
は、
ε12=1,
ε21=-1
であり、それ以外の
εmp=0
です。
余因子行列Bstは、具体的には、
B11=A22
B12=-A12
B21=-A21
B22=A11

こうして計算した余因子行列Bstを行列Ampに掛け算すると、
行列Ampの行列式=Δとすると、
B1mAm1=εmpAm1Ap2=Δ
B1mAm2=εmpAm2Ap2=0
B2pAp1=εmpAm1Ap1=0
B2pAp2=εmpAm1Ap2=Δ
∴
BstAtp=Δ・Esp
このEspは、単位行列。

 (2)
3行3列の行列Ampの余因子行列Bstは、
エディントンの行列式の計算記号εmpsを使って、アインシュタインの縮約記法であらわして、
B1m=εmpsAp2As3
B2p=εmpsAm1As3
B3s=εmpsAm1Ap2
で計算します。
(補足)これらの式は、行列式εmpsAm1Ap2As3の計算式から、余因子行列Bstに対応して、行列Aのs列目の項を除去した式です。
この計算に利用する行列
εmpr
は、
ε123=ε231=ε312=1
ε213=ε321=ε132=-1
であり、それ以外の
εmpr=0
です。

こうして計算した余因子行列Bstを行列Ampに掛け算すると、
3行3列の行列Ampの行列式=Δとすると、
B1mAm1=εmpsAm1Ap2As3=Δ
B1mAm2=εmpsAm2Ap2As3=0
・・・
結局、
BsmAmp=Δ・Esp
このEspは、単位行列。

(3)
このように、余因子行列Bsmは、
行列Ampの左から掛け算することで、Δ・Esp
となる行列です。

(4)
また、行列Ampの右から掛け算する場合も、同じことになります。
そうなるのは、行列式が以下の式でもあらわせることに由来します。
つまり、2行2列の行列式Δは、
Δ=εmpAm1Ap2
Δ=εmpA1mA2p
ともあらわせ、
また、3行3列の行列式Δは、
Δ=εmpsAm1Ap2As3
Δ=εmpsA1mA2pA3s
ともあらわせることに由来します。
(この証明は省略)

以下では、余因子行列Bsmの右側からの掛け算でも、この関係を満足することを確認するために、2行2列の行列を例に、直接計算してみます。
B11=A22
B12=-A12
B21=-A21
B22=A11
この余因子行列Btmの左から行列Astを掛け算すると、
A1tBt1=A11B11+A12B21=A11A22-A12A21=Δ
A1tBt2=A11B12+A12B22=-A11A12+A12A11=0
A2tBt1=A21B11+A22B21=A21A11-A22A21=0
A2tBt2=A21B12+A22B22=-A21A12+A22A11=Δ
∴
AstBtm=Δ・Esm
このように、余因子行列Btmは、
行列Astの右から掛け算しても左から掛け算しても、Δ・Esm
となる行列です。

(5)
ベクトルBt1に対して、
(ベクトルBt1は、(B11, B21)です。)
AstBt1=Δ・Es1
そのため、行列Astの行列式Δ=0の場合、
行列Astは、ベクトルBt1を0ベクトルに変換します。
このベクトルはこの行列の固有値0に対する固有ベクトルです。
また、ベクトルBt2に対しても、
AstBt2=Δ・Es2
そのため、行列Astの行列式Δ=0の場合、
行列Astは、ベクトルBt2も0ベクトルに変換します。

2行2列の行列Astの行列式Δ=0の場合、
ベクトルBt1とベクトルBt2は、同じ方向を向いた平行なベクトルです(証明は省略)。そのため、そのうち一方だけで(ただし、それが0ベクトルで無いベクトルを選ぶこと)、0ベクトルに変換されるベクトルの代表にすることができます。

このように、行列Astの行列式Δ=0の場合、
その余因子行列の列ベクトルBt1(又はBt2)は、行列Astが0ベクトルに変換します。

しかし、行列Aが、
1,1,1
0,0,0
0,0,0
のような行列の場合、
その余因子行列Bは、全ての成分が0である0行列になる。
その行列Aは、縦ベクトル(1,-1,0)や縦ベクトル(1,0,ー1)を0ベクトルに変換する。
しかし、その縦ベクトルを余因子行列Bの各列ベクトル(0ベクトル)の合成によって作ることができない。そのように、行列Aによって0ベクトルに変換される縦ベクトルのうち、余因子行列Bの縦ベクトルから作れないベクトルもあるので注意すること。




リンク:
追加講:三角形の面積と行列式
高校数学の目次