\def\Mtt#1#2#3#4{\left(\begin{array}{rr}#1 \\ #3\end{array}\right)} \def\Mto#1#2{\left(\begin{array}{r}#1 \\ #2\end{array}\right)} \def\Mss#1#2#3#4#5#6#7#8#9{\left(\begin{array}{rrr}#1 \\ #4 \\ #7 \end{array}\right)} \def\Mqqa#1#2#3#4#5#6#7#8{\left(\begin{array}{rrrr}#1 \\ #5 \\} \def\Mqqb#1#2#3#4#5#6#7#8{#1 \\ #5\end{array}\right)} \def\sth{\sin\theta}\def\cth{\cos\theta}\def\Sth{S\theta}\def\Cth{C\theta} \def\sph{\sin\phi}\def\cph{\cos\phi}\def\Sph{S\phi}\def\Cph{C\phi} \def\sal{\sin\alpha}\def\cal{\cos\alpha}\def\Sal{S\alpha}\def\Cal{C\alpha} \def\sps{\sin\psi}\def\cps{\cos\psi}\def\Sps{S\psi}\def\Cps{C\psi}

東北学院大学工学部でＩＴ教育の試みの一つとして配布した KNOPPIX EduTG は、ＣＤ一枚から立ち上がる、OS(Linux)と様々なアプリケーションプログラムをセットにしてドイツで開発されたKNOPPIXをつくばで日本語化したものに、さらにＩＴ教育委員会の要望を組み込んで開発されたものです。
そのため、一般的なKNOPPIXには含まれていないものの工学部で使いそうなツールも搭載してあります。

このページでは、せっかく配布されているKNOPPIXをロボット工学の助けに使わない手はないということで、活用法をまとめていきます。

なお、このページの例は、

家でブロードバンドルータ等を使っている場合はKNOPPIXを立ち上げれば、たぶんMozillaでホームページを見ることができるので、みながら＆コピペする。
大学でやるばあいは(指導教員他の許可取ること)、となりのWindowsパソコンと見比べながら。
USBメモリにごそっとダウンロードしておいてあとで KNOPPIXにもっていって

等で試すことができます。 (紙に書くとたいへんでしょう。＜敢えて印刷は避ける方向)

主要アプリケーションの一つであるOpenOfficeに含まれる表計算ソフト、Calcは簡単なロボット制御シミュレーションができます。
表計算ソフトでシミュレーションをやるときの、基本的な手順です(筆者独断)。

まず、各状態変数の関係式をリストアップする。
横方向に各状態変数をおいていく。
状態方程式のように、「状態変数の微分値＝」で式が与えられるときは数値積分が必要なので、状態変数そのものと、その微分値の枠を取る。
縦方向に時間軸をとる。左端の列は時刻。
単純な変数間の関係式をまず時刻０の行に入力する。
システムへの入力値は適当に設定しておく。
積分していくような変数は、初期値を書きこねおく。
入力した１行を次の行にコピーする。
時刻＝１行上の時刻＋Δt
変数＝一行上の変数＋変数の微分値×Δt
の関係を入力。
以上を下方向にだぁ〜っとコピーして、入力を適宜設定して得られた結果を吟味する。

Δtは特定のセル(枠)に書いておいて、「$行$列」という固定枠表記で数式中で使います。すると、そこを変えるだけで計算への時間ステップの影響が確認できます。同じように制御ゲインなども固定枠に書いておきます。実際に、この車輪移動ロボットのページのアニメーションなどは、OpenOfficeのCalcでまずは動作確認および軌道計算や車輪の回転などを確認して、それを３次元ＣＧソフトのデータに結合してつくっています。そのなかから、円弧走行、

のデータを見てみます。

変数の名前を適当につけていたため、ちょっと分かりにくいかも知れませんが、

VR,VR: 左右車輪の進行速度
DR,DL: VR,VLの積分値＝両車輪の移動経路長
V,omega: 車両の平均速度および角速度
Theta,Theta': 角速度を積分したものと、その中間値。ロボットの姿勢を示す。
X,Y: 車両の位置。VをTheta'方向に積算。

となっています。また、$D$3にはロボットの車幅(ωの計算に必要)、$D$4には時間刻がはいっています。
このシートでは、VR,VLを適宜設定するとその際の軌道が計算されます。 X,Yをグラフにすれば、実際の軌跡がわかるわけです。なお、グラフにするときには、縦軸と横軸が同じ縮尺になるように気をつけましょう。違う縮尺だと、不要な勘違いを起こす危険性があります。
ＣＧ映像を作るためには、これに加えて車輪の回転角度や旋回中心の表示などいろいろ追加で計算しています。

この講義では特に前半、やたらとcos/sinのつまった行列の演算を扱いました。これらはある程度は手計算でできますが、コンピュータがやってくれるならまかせるにこしたことはありません(もちろん、手を動かすことは重要)。

もっとも、本来コンピュータは数値計算は得意ですが、数式をそのまま処理するということは容易ではありません。数式の解釈や展開、単純化は簡単ではないのです。しかし、KNOPPIX EduTGに含まれているMaxima(xmaxima)というソフトウエアが数式処理をしてくれますので、この講義で出てきた演算を例につかってみます。

XMaximaは起動後、メニュー→数学演習→xmaximaで起動します。下半分には簡単な使い方がでます。青い表示の数式をクリックするとその演算が実行されます。解説文は英語ですが、微積の使い方などはすぐにわかるとおもいます。ここではそれらは飛ばして、実際に行列演算をやってみます。参考：Maximaのマニュアル＠狸穴プチテストで出題した問題、

Ｑ：
θ1およびθ2で回転させる行列

に対して、

を計算し、結果を論じよ。

Ａ：
計算については

の結果を論じることは無理としても、計算だけやらせてみます。

まず、行列を定義します。Maximaでの行列は

A:matrix([a,b],[c,d]);        行列Aを定義
B:matrix([e,f],[g,h]);        行列Bを定義
A+B;                          A+Bを演算
A.B;                          行列積 AB
A*B;                          注意これは各要素ごとに乗じる:使えない

という感じで演算します。注意点は行列積でうっかり＊をつかうと違う結果になることです。「．」をつかいます。

さて、上の例題の通りの定義は

A:matrix([cos(th1),-sin(th1)],[sin(th1),cos(th1)]);
B:matrix([cos(th2),-sin(th2)],[sin(th2),cos(th2)]);
A.B;

です。実際にXmaximaの上の枠で入力してEnterすると妥当な結果が得られます(最後の「;」に注意。忘れるとちゃんと動かない)。
さて、この結果、積和公式でまとめて欲しいものです。それを行うには

TRIGREDUCE(%);

と入力します。「%」は直前の結果を使うという意味です。手を抜かないでやるならば、「TRIGREDUCE(A.B);」とします(TRIGREDUCEは積和公式などで単純化を試みる)。ここまでの実行例を下に示します。
次は、オイラー角の回転行列です。座標変換のページから抜粋。

上記オイラー角での回転を数式で考えてみましょう。

まず、Ｚ軸回りの回転(0→1)の角度をφ(phi, ファイ)として、

となります。ただし、

と置き換えています(この先面倒なので)。
つぎにＸ軸回り(1→2)の角度をθとして、

を得ます(同じ置き換えパターン)。最後のＺ軸回り(2→3)の角度をψ(psi, プサイ)として、

となります。

角度(φ, θ, ψ)で表される回転を総合して得られる回転行列は、

となります。

これをXMaximaにかけてみます。

R01:matrix([AC,-AS,0],[AS,AC,0],[0,0,1]);
R12:matrix([1,0,0],[0,BC,-BS],[0,BS,BC]);
R23:matrix([CC, -CS,0],[CS,CC,0],[0,0,1]);
R01.R12.R23;

ただし、簡単にするために、φ、θ、ψをそれぞれ、A,B,Cとおき、しかも、cos(θ)をBCとおいています。これは小細工で、maximaがこちらの意図とは関係なく、アルファベット順に変数を並べてしまうため、角度順に並ぶように敢えてこういう表記をしました。
結果は、右図のようになりました。良く見比べると、式が一致していることが確認できます。 (xmaximaがあっていたのではなく、このページがあっていた、といったほうがいいかも。もちろん、xmaximaの使い方が正しかったとも言えます。) 次はほとんど同じですが、レポートに出したロールピッチヨーの方です。

Ｑ：
Ｚ-Ｙ-Ｘの順番で回転させるロールピッチヨー角回転変換の変換行列

を求めよ。

Ａ：

検算例：
３行：

１行：

２行・３行：

まずは、同じように

R01:matrix([AC,-AS,0],[AS,AC,0],[0,0,1]);
R12:matrix([BC,0,BS],[0,1,0],[-BS,0,BC]);
R23:matrix([1,0,0],[0,CC,-CS],[0,CS,CC]);
R01.R12.R23;

をしてみます。

答えはちゃんとあいます。

次に、せっかくなので検算も行ってみます。
検算するには、最終的に１になるか０になるかをチェックしたいので、ACなどの表記ではだめです。面倒ですが、ちゃんとcos/sinで打ち込みます。

R01:matrix([cos(th1),-sin(th1),0],[sin(th1),cos(th1),0],[0,0,1]);
R12:matrix([cos(th2),0,sin(th2)],[0,1,0],[-sin(th2),0,cos(th2)]);
R23:matrix([1,0,0],[0,cos(th3),-sin(th3)],[0,sin(th3),cos(th3)]);
R03:R01.R12.R23;
ROW(R03,1).ROW(R03,1);
TRIGSIMP(%);
TRIGSIMP(ROW(R03,1).ROW(R03,2));

結果：

ここでは、行列R03をR01,R12,R23の積としてまず定義しています。結果は前と同じ意味ですが、大きくなり過ぎたので列ごとに分解されて表示されたようです。
次に、ROWで行列の行を抜いてきます。ROW(R03,1)でR03の１行目がとれます。これを二つとって内積をとると、回転行列なので１になるはずです。ただ内積をとるとごちゃっと出てきますが、簡略化の処理、TRIGSIMPを使用すると１になります。異なる行(ここでは１と２)を取ってきて内積をとった場合は０になります。 ROWのかわりに列をとるCOLを使っても、同列なら１、異列なら０になります。

数式処理ソフトは使い方さえ間違わなければ強力です。ですが、うちまちがったり、ちょっと勘違いしただけで当然結果はかわりますので、過信は禁物です。必ず妥当性をチェックしながら進みましょう。また、これに頼ると中間結果に注意を払わなくなってしまいます。やっぱりレポートは手で計算して、最後にこれでチェックするくらいがいいでしょう。