射影作用素
表示
(直交射影から転送)
線型代数学および函数解析学における射影作用素あるいは単に射影︵しゃえい、英: projection︶とは、いわゆる射影︵投影︶を一般化した概念である。有限次元ベクトル空間 Vの場合は、V 上の線型変換 P: V→ Vであって、冪等律 P2= Pを満たすものを言う。ベクトル vの像 Pvを vの射影という。射影作用素はベクトル空間 Vを U⊕W と直和分解したときに、V の元 v= u+ w(u ∈ U, w∈ W) を uに写すような変換である。ベクトル空間の次元が無限次元の場合には、連続性を考慮しなければならない。例えばヒルベルト空間 における射影作用素とは、 上の有界線型作用素 であって、冪等律 P2= Pを満たすものを言う。このときさらに自己共役性 P∗ = Pを持つときには直交射影︵ちょっこうしゃえい、英: orthogonal projection︶という[1]。直交射影のことを単に射影と呼ぶこともある[2]。
この定義は抽象的ではあるが、投影図法の考え方を一般化し、定式化したものになっている。 幾何学的対象上の射影の影響は、その対象の各点における射影の影響を調べることでわかる。
平易な例
[編集]直交射影
[編集]
例えば、三次元空間 R3の点 (x, y, z) を点 (x, y, 0) へ写す写像は xy-平面の上への射影である。この写像は行列
によって表現される。実際、この行列 Pの任意のベクトルへの作用は
となり、これが射影を定めること︵つまり P= P2を満たすこと︶は
なる計算によって確かめられる。
斜交射影
[編集]
︵定義は後述するが︶直交でない︵斜交︶射影の簡単な例として
を挙げることができる。行列の積の定義に従って計算すれば
故に Pが実際に射影となることが分かる。
この射影 Pが直交射影となるのは α = 0 のときであり、かつそのときに限る。
分類
[編集]
以下、本節において考えるベクトル空間はすべて有限次元であるものと仮定する︵この場合、射影の連続性などを気にしなくともすむ︶。
変換 Tは kに沿った mの上への射影である。T の値域は mで あり、T の零空間は kに等しい。
本項冒頭の導入文で述べたとおり、射影 Pは冪等律すなわち P2= Pを満たすような線型変換である。
もととなるベクトル空間を Wとする。W の部分線型空間 Uおよび Vが、それぞれ Pの値域および零空間︵核︶であるものと仮定すると、基本的な性質として
●P は U上に恒等作用素 Iとして作用する。つまり、
●直和分解 W= U⊕ Vが成立する。すなわち、W の各ベクトル xは Uの元 uと Vの元 vを用いて x= u+ vなる形に一意的に表される。これには
とすればよい。
などが成り立つことがわかる。射影の値域と核は互いに﹁相補的﹂なもので、P と Q= I − P も同じく﹁相補的﹂である。すなわち、作用素 Qもやはり射影を定め、Q の値域は Pの核、Q の核は Pの値域となる。逆もまた然り。
このとき Pを︵核︶V に沿った︵値域︶U の上への射影と言い、また Qを Uに沿った Vの上への射影と呼ぶ。
ベクトル空間の、部分空間の直和への分解は一般には一意的でない。従って、部分空間 Vが与えられたとき、その値域︵もしくは核︶が Vとなるような射影は一般に複数存在しうる。
射影のスペクトルが {0, 1} に含まれることは
から分かる。射影の固有値となれるのは 0 および1に限られるが、それらに対応する固有空間は射影の核および値域に他ならない。
自明でない射影は最小多項式が となり、これは相異なる一次因子の積となっているから、P は対角化可能である。
直交射影
[編集]
考えているベクトル空間に内積が定義されていれば、直交性や︵線型作用素の自己共軛性︶といったような内積に付随するさまざまな概念を用いることができるようになる。直交射影は、値域 Uと核 Vとが互いに直交する部分空間になっているような射影をいう。射影が直交射影であるための必要十分条件は、それが自己共軛であること、即ち実ベクトル空間の場合には、ある直交基底に関する表現行列 Pが対称行列︵P = PT︶であり、複素ベクトル空間の場合には、表現行列 Pがエルミート行列︵P = (P*)T)︶となることである。実際に、x, yが射影の定義域に属するベクトルのとき、 Px∈ U, y− Py∈ Vであり、かつ を正定値内積として
が成り立つから、Px と y− Pyとが任意の x, yに関して互いに直交するのは、P = PTP︵これは P= PTかつ P= P2に同値︶のときであり、かつそのときに限る[3]。
直線の上への直交射影の場合が最も簡単であろう。直線上の単位ベクトル uをとれば、当該の射影は
で与えられる。この作用素は uを変えないし、また uに直交する全てのベクトルを零化する。このことは、u を含むどんな直線の上への射影についても正しい[4]。これを見るのに簡単な方法は、勝手なベクトル xを直線上の成分︵つまり射影されたベクトルを考える︶とそれに垂直な成分との和
と考えることである。これに射影を施せば、平行なベクトル同士の内積と垂直なベクトル同士の内積の性質から
を得る。
この等式は任意次元の部分空間の上への直交射影にも拡張することができる。u1, ..., ukを部分空間 Uの正規直交基底とし、各列ベクトルが u1, ..., ukになっている k-次正方行列を Aと書けば、所期の射影が
で表される[5]。これは内積を使えば
と書くこともできる。行列 ATは Uの直交成分が消える部分等距変換であり、A は Uを考えている全体空間へ埋め込む等長変換になっている。従って PAの値域は Aの終空間 (final space) であり、また ATAが U上の恒等変換であることは明らかである。
上記の議論で正規直交条件は落とすこともできる。即ち、u1, …, ukを︵必ずしも正規直交でない︶基底とし、それらを列ベクトルに持つ行列を Aと書けば、求める射影は
と書ける[6]。この場合も行列 Aは Uの全体空間への埋め込みになっているが、しかし一般にはもはや等距変換ではない。ここで行列 (ATA)−1 はノルムを回復する﹁正規化因子﹂である。実際、階数1の作用素 uuTは ‖u‖ ≠ 1 のとき射影にならないが、これをuTu = ‖u‖2 で割って得られる u(uTu)−1uT は uで張られる部分空間の上への射影になる。
この射影の値域となるベクトル空間が︵基底ではなくて︶枠 (frame) で張られているとき︵つまり生成元の数が次元の値よりも大きいとき︶には、上記の公式は
という形になる。ここで はムーア・ペンローズ擬似逆行列を表す。このような場合には、射影作用素を構成する方法は無数にあり、これはその無数の可能性のうちの一つに過ぎないことに注意すべきである。
あるいは、行列 が正則で ATB= 0︵つまり、B は Aの零空間行列︶のときには
が成り立つ。直交条件を強めて、正則行列 Wに対して ATWB= ATWTB= 0 が成り立つものとすれば、
が成立する。
これらの公式は︵転置行列を随伴行列に取り替えれば︶複素内積空間でも成立する。
斜交射影
[編集]
直交射影でないような射影のことを、斜交射影と呼ぶこともある。直交射影ほど頻繁ではないが、この種の射影は二次元に描画された空間図形を表すのにも用いられる。
斜交射影はその値域と核によって定まり、与えられた値域と核を持つ射影の行列表現の式は次のように求められる。まず射影の値域の基底を成すベクトルを u1, …, ukとし、それらを列ベクトルとして並べた n× k行列を Aと書く。射影の値域と核とは互いに補空間になっているから、核の次元は n − k である。従って、射影の核の直交補空間の次元は kであり、v1, …, vkがその基底を成すものとして、それらを並べた行列を Bと書く。このとき、当該の射影は
によって定まる。この公式を、上で直交射影に対してやったように拡張することもできる[7]。
標準形
[編集]
体上の d-次元ベクトル空間上の射影 P= P2は、その最小多項式が x2 − x で相異なる一次因子の積に分解されるから、対角化可能である。従って、適当な基底を選べば Pは、r を Pの階数として
なる形に表すことができる。ここで、Ir は r-次単位行列、0d−r は次数 d − r の零行列である。複素ベクトル空間で内積を持つ場合には、適当な正規直交基底を選んで、P の表現行列を
なる形にすることができる[8]。ただし、σ1 ≥ σ2 ≥ … ≥ σk > 0 とする。また、k, s, mは整数で、実数 σi は一意に定まる。2k + s+ m= dであることに注意せよ。このときの、Im ⊕ 0s なる因子は、その上に Pが直交射影として作用する最大の不変空間に対応しており︵故に P自体が直交射影となるのは k= 0 のとき、かつそのときに限る︶、かつ σi-ブロックが Pの斜交成分に対応している。
ノルム空間上の射影作用素
[編集]
考えるベクトル空間 Xが︵有限次元とは限らない︶ノルム空間のとき、︵有限次元の場合には関係ないが︶解析学的なことも考えないといけないので、ここでは Xはバナッハ空間であることを仮定する。
先に述べた代数的な概念の多くはこの文脈においても有効である。例えば、互いに補空間となるような部分空間への Xの直和分解が与えられればやはり射影が定まるし、逆に射影からそのような直和分解が得られる。実際、X が直和分解 X= U⊕ Vを持つとき、P(u + v) = uで定義される作用素はやはり値域 Uおよび核 Vの射影である︵P2 = Pは明らかである︶。一方 Pが X上の射影、即ち P2= Pを満たすならば (I − P)2 = (I − P) は容易に確かめられ、即ち (I − P) もまた射影となる。関係式 I= P+ (I − P) から Xが Ran(P) ⊕ Ran(I − P) なる直和に分解されることが従う。
しかし、有限次元の場合とは対照的に、射影は一般に連続とは限らない。実際、X の部分空間 Uがノルムの定める位相に関して閉でないときは Uの上への射影は連続でない。同じことだが、連続な射影 Pの値域は必ず閉部分空間でなければならない。更には、連続射影の︵実は一般の連続線型作用素の︶核は閉部分空間である。従って、連続射影 Pは Xの互いに補空間となる閉部分空間の直和への分解 X= Ran(P) ⊕ Ker(P) = Ran(P) ⊕ Ran(I − P) を与える。
逆は、適当な仮定を追加すれば成り立つ。U を Xの閉部分空間とすると、X = U⊕ Vとなる閉部分空間 Vが存在する場合に限り、値域が U, 核が Vとなる射影 Pは連続である。これは閉グラフ定理から従う。即ち、xn → xかつ Pxn→ yとするとき、Px = yが示されればよい。U が閉で、{Pxn} ⊂ Uだから yは Uに属し、Py = yが成り立つ。また、xn − Pxn= (I − P)xn → x− yである。このとき、V は閉で {(I − P)xn} ⊂ Vだったから、x − y∈ V即ち P(x − y) = Px− Py= Px− y= 0 を得て、主張が示される。
今の議論では U, Vがともに閉であるという仮定が効いているが、閉部分空間 Uが与えられたときにその閉補空間 Vの存在は一般には保証されない。ただし、ヒルベルト空間では直交補空間をとることで常にそれができる。バナッハ空間の場合には、一次元部分空間が常に閉補空間を持つことが、ハーン・バナッハの定理から直ちに従う。実際、U を uが張る一次元部分空間とすると、ハーン・バナッハから、有界線型汎函数 φ で φ(u) = 1 なるものがとれる。このとき、作用素 P(x) := φ(x)u は P2= Pを満足し、射影となる。φ の有界性から Pの連続性が出るから、従って Ker(P) = Ran(I − P) が Uの閉補空間となる。
そうは言うものの、開写像定理により、バナッハ空間上の任意の連続射影は開写像であることが言える。
応用およびさらに進んだ議論
[編集]
射影︵直交射影とその他︶は、線形代数の問題でのいくつかの計算アルゴリズムにおいて、重要な役割を果たす。
●QR分解 (ハウスホルダー変換とグラム・シュミットの正規直交化法を参照)
●特異値分解
●ヘッセンベルク形式 への変換 (多くの固有値計算アルゴリズムでの最初の処理)
●線形回帰
上で述べたように、射影というのは冪等作用素の特別なものであり、解析学的には直交射影は特性函数の非可換な一般化になっている。可測集合の特性函数を考えることから測度論が始まったように、冪等作用素は︵例えば半単純多元環などの︶分類にも用いられ、それゆえ想像のつくとおり、射影作用素も作用素環論の文脈で極めて頻繁に用いられる。特に、フォン・ノイマン環はその射影の成す完備束によって生成される。
物理への応用
[編集]
量子論では,ある条件を満たす状態の全体は状態空間の部分空間と考えることができるので,量子力学的な命題と部分空間,すなわち射影演算子とを対応させることができる(量子論理).
統計力学では、運動の粗視化を射影演算子を使って定式化する方法(射影演算子の方法)がある.
分子対称性、分子振動、格子振動、結晶の波動関数では、任意の関数からある対称性に従う関数のみを作りたい時に、射影演算子が用いられる。たとえば射影演算子を用いれば、既約表現の表現行列からその基底関数︵基準振動、基準モードなど︶を求めることができる。
一般化
[編集]
より一般に、ノルム空間の間の写像 T: V→ Wが与えられたとき、同じようにこれが核の直交補空間上の等距写像となることを要求することができる。その は等距であり、特に全射でなければならない。直交射影の場合というのは Wが Vの部分空間であるときである。リーマン幾何学においてこのことはリーマン沈め込みの定義に使われている。
注釈
[編集]- ^ Reed & Simon 1980, p. 187.
- ^ Reed & Simon 1980, p. 188.
- ^ Meyer 2000, p. 433.
- ^ Meyer 2000, p. 431.
- ^ Meyer 2000, equation 5.13.4.
- ^ Meyer 2000, equation 5.13.3.
- ^ Meyer 2000, equation 7.10.39.
- ^ Doković, D. Ž. (August 1991). “Unitary similarity of projectors”. Aequationes Mathematicae 42 (1): 220–224. doi:10.1007/BF01818492 .
参考文献
[編集]
●Dunford, N.; Schwartz, J. T. (1958). Linear Operators, Part I: General Theory. Interscience. Zbl 0084.10402
●Meyer, Carl D. (2000). Matrix Analysis and Applied Linear Algebra. Society for Industrial and Applied Mathematics. ISBN 978-0-89871-454-8. Zbl 0962.15001
●Reed, Michael; Simon, Barry (1980). Methods of modern mathematical physics I: Functional analysis (Rev. and enl. ed.). Academic Press. ISBN 0-12-585050-6. MR0751959. Zbl 0459.46001
●齋藤, 正彦﹃線型代数入門﹄︵初版︶東京大学出版会︿基礎数学1﹀、1966年。ISBN 978-4-13-062001-7。
●黒田, 成俊﹃関数解析﹄共立出版株式会社︿共立数学講座15﹀、1980年。ISBN 978-4-320-01106-9。
関連項目
[編集]- Centering matrix, which is an example of a projection matrix.
- 直交化
- 不変部分空間
- 跡
- ダイクストラの射影アルゴリズム (Dykstra's projection algorithm): to compute the projection onto an intersection of sets
外部リンク
[編集]- MIT Linear Algebra Lecture on Projection Matrices at Google Video, from MIT OpenCourseWare
- Planar Geometric Projections Tutorial - a simple-to-follow tutorial explaining the different types of planar geometric projections.