逆行列を求めるプログラムプログラムで行列演算をしていると逆行列が必要になることがしばしばあります。 2×2の逆行列なら簡単に求まるからいいんですが、それより大きな行列になると公式らしい公式が見当たらないので仕方なくLU分解のルーチンを組む、なんてことを私は何回か経験しました。 ただ、そのLU分解のルーチンが結構厄介で、ピボット操作などもちゃんとルーチンに組み込んでおかないとゼロ除算とかが発生していまい、本当は逆行列が存在するにもかかわらず算出できない、なんてことになってしまいます。 そもそもこのテの数値演算のアルゴリズムは数式上では正確でもコンピュータに計算させると大きな誤差が発生してしまうので困ったもんです。 ごちゃごちゃしたプログラムを組まなくても逆行列が出せないかなぁ、と思ってちょっと考察してみました。
