バックナンバーはこちら。
https://www.simulationroom999.com/blog/compare-matlabpythonscilabjulia2-backnumber/
はじめに
正規方程式を導出するまでの説明。
今回は二次形式の続き。
登場人物
博識フクロウのフクさん

イラストACにて公開の「kino_k」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=iKciwKA9&area=1
エンジニア歴8年の太郎くん

イラストACにて公開の「しのみ」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=uCKphAW2&area=1
ロードマップ【再掲】

まずはロードマップを再掲。


今回は二次形式の話の続きで、行列表現と多項式表現の可逆性維持の話になる。
二次形式の行列表現と多項式表現

前回、二次形式を行列表現できたが、
これをもう少し具体的に書いてみる。
\(ax^2+bxy+cy^2\)の場合
\(
\begin{bmatrix}
x \\
y
\end{bmatrix}^T
\begin{bmatrix}
a && c/2 \\
c/2 && b
\end{bmatrix}
\begin{bmatrix}
x \\
y
\end{bmatrix}=
ax^2+by^2+cxy
\)
\(ax^2+by^2+cz^2+dxy+eyz+fxz\)の場合
\(
\begin{bmatrix}
x \\
y \\
z
\end{bmatrix}^T
\begin{bmatrix}
a && d/2 && f/2 \\
d/2 && b && e/2 \\
f/2 && e/2 && c
\end{bmatrix}
\begin{bmatrix}
x \\
y \\
z
\end{bmatrix}=
ax^2+by^2+cz^2+dxy+eyz+fxz
\)

対称行列にする意味ってあるの?

多項式から行列にする場合、対称行列という制約を設けないと、相互変換できないってのはあるな。

なるほど。
制約があるから可逆性が保証出来て、結果として表現が違うだけで同じものと言えるってわけか。

そうそう。
今後について

で、ここまでで二次形式の行列表現、多項式表現が等しいことを証明したわけだが、
これが本当に同じなのかを各ツール、各言語で確認するのを挟もうと思う。

そうだねー。
なんか、ここまででも結構ごちゃごちゃした感じがするもんねー。

そこで、一点問題がある!

ずっと問題だらけやー!

現状の行列表現のままだと、行列を利用した一括変換ができない。
まぁfor文使えば解決ではあるのだけど、
折角だから、ベクトル、行列であるが故の利点を享受しておきたい。
というわけで、この話を次回やってから、実際に各ツール、各言語で計算させてみる。
って流れになると思う。

もう好きにしてくれ・・・。
まとめ

- 二次形式の行列表現と多項式表現の関係性を示した。
- 行列Aが対称行列になることを制約とすると行列表現と多項式表現に可逆性を付加することができる。
- これに伴い、等しいものとして扱うことが可能となる。
- 実際に行列表現と多項式表現が等しいかを各ツール、各言語で確認したいが、その前にベクトル、行列のまま演算する仕掛けを考える必要がある。
バックナンバーはこちら。
コメント