... 分解と呼ぶ1
これを Schmidt 分解、あるいは Gram-Schmidt の分解と呼ぶことを提案する本もある (シャトラン [1], 齋藤 [9])。 以下に見るように (1) の形の分解を得るアルゴリズムには 色々あるので、個人的には、 有名なアルゴリズムの名前に冠される ``Gram-Schmidt'' は適切ではないように 感じている。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... (と言われている2
øA.Björk, Solving linear least squares problems by Gram-Schmidt orthogonalization, BIT, Vol.7 (1967), pp.1-21.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... この本に書いてあることは今では世界の常識みたい3
Trefethen and Bau III の本の実験の mgs.m, clgs.m という名前の プログラムが世の中には一杯あるみたい。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... より4
これから、 相似変換で trace が不変であることが導かれる。 今必要なのはそれではないけれど。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
...$ 行列で$ 5
$B=(&beta#beta;_ij)$ が $j>i+1$ $$ $&beta#beta;_ij=0$ をみたすとき、$B$ は Hessenberg 行列という。
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... $ と仮定する)$ 6
MATLAB 風に書けば、 $B=A(1:n-1,1:n-1)$.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
... 以下のブロック上三角行列7
最初にこれを聞いたとき、Hessenberg 行列のことかと思ったが、 違う!
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.