= 2019/03/28 =

What we did

 * センサーの対角化を進める。

=== 数学的背景 ===
 * 進める前に，センサーの対角化について数学的にどうなるかをまとめておこう。
 * ゴールは，「各自由度が分離できる」「自由度ごとにどれだけ振れるかの大きさを把握する」こと
 * 後述するが，２つ目の目標は，絶対的な大きさまで把握しようとすると難しい。
 * TMS-VISは剛体の６自由度だが，まず簡単のために水平方向の３自由度(Longitudinal, Transverse, Yaw)で考える。
 * センサーの幾何学的配置からある程度自由度が分離できているとして，その自由度の基底を(e_Lg, e_Tg, e_Yg)^tとする。このセンサーで得た各自由度のスペクトルを(L_1, T_1, Y_1)
 * 目標の理想的な，真の自由度の基底を(e_Li, e_Ti, e_Yi)^tとする。真の自由度で測定できた場合の各自由度のスペクトルを(L_0, T_0, Y_0)
 * 対角化のイメージとしては，これら(L_1, T_1, Y_1)と(L_0, T_0, Y_0)との変換行列を求めることである。
{{attachment:matrix.png}}

 * ここで，幾何学的配置からのセンサーは真の自由度に近いとするという近似をする。つまり，上の行列の対角成分は1であるとする。
 * また，真の自由度からのスペクトルは測定できないが，それに対して何かしら定数倍がかかっているものは測定できるとする。

{{attachment:diagimage.png}}

 * ここで，f_DOFは各自由度の共振周波数である。
 * よって各共振周波数のピークの比で構成されていることがわかる。
 * この行列の逆行列を求めれば良い。

=== 手順 ===
 1. アクチュエータの対角化行列を求めて，E2Cに挿入。
(19/04/12 アクチュエータの対角化行列乗せてなかったので追記)
今回，以下のようなE2Cの対角化行列Dを使用している。

{{attachment:act_diag.png||width=600}}
 2. 自由度ごとにwhite-noiseを入れて励起自由度のPSDとその他の自由度のPSDとの比をとる。
 3. 励起自由度とその他の自由度との位相伝達関数から符号を決める。~0度なら+, ~180度なら-
 4. PSDの比から上の行列を作って，逆行列を作る。