回転前後の座標が分かっている場合の回転中心座標を求める
剛体の移動を追跡する実験をしたときに、その剛体の2点に測定点を貼り付けてカメラ等で動画撮影する事がよくあります。この時、その剛体の特性がよく分かってない場合、その回転中心(もし拘束条件がなければ重心)はどのように求められるか?という解法について、とある事情で必要になったのでメモがてら書いておきます。
具体的に言うと、下図のような2点p0, p1が、ある点c周りに回転してP0, P1になった時、cの座標と回転角θを求めたいということです。
結果は以下の通りになりました。
以下証明など
続きを読むヤコビ反復法が収束するための必要十分条件(2次元の場合)
ヤコビ法(ヤコビの反復法)を使って線形方程式
の求解する時の収束条件について、世間一般では「係数行列Aが対角優位であること」ぐらいしか書いていない。しかしこれは十分条件であって、つまり、対角優位でなくても収束することがある。実際、
という係数行列は、2行目で対角成分2<非対角成分3なので対角優位ではないが、実際には収束する。
では必要十分条件はなんだろうか?一般の場合は難しそうなのでとりあえず2次元で考えてみたところ導出できたのでメモがてら。
結論だけ先に書くと、2次元でヤコビ反復法が収束するための必要十分条件は、対角成分の積の絶対値が非対角成分の積の絶対値より大きいことだった。つまり
に対して、であれば、ヤコビ反復法は収束するということである。
以下証明(というか考え方)。
続きを読む