MATLAB,Python,Scilab,Julia比較 第2章 その37【二次形式の微分①】

MATLAB,Python,Scilab,Julia比較 第2章 その37【二次形式の微分①】 数値計算
MATLAB,Python,Scilab,Julia比較 第2章 その37【二次形式の微分①】

バックナンバーはこちら。
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

二次形式の微分の位置づけ

フクさん
フクさん

今回からは二次形式の微分の話に突入する。
ロードマップを再掲しておくので、位置づけを確認しておくと良いだろう。

正規方程式に至る道、二次形式、対称行列、二次形式の微分、グラム行列、二乗和誤差、正規方程式
太郎くん
太郎くん

前回まででやってきた、二次形式と対称行列を利用する感じなのかな?

フクさん
フクさん

そうだね。
二次形式の係数部を対称行列と定義しているため、計算結果がシンプルになる。

太郎くん
太郎くん

どうシンプルになるのかさっぱりだけど、
そこは話を聞いてみればわかるかな。

説明手順

フクさん
フクさん

とりあえず、説明手順を書き出しておこう。

  • \(\nabla\)(ナブラ)について
  • 二次形式の微分(勾配)
  • 適当な多項式に当てはめてみる
  • 各ツール、各言語で確認してみる。
太郎くん
太郎くん

\(\nabla\)(ナブラ)とか謎の言葉が・・・。

∇(ナブラ)

フクさん
フクさん

まず、\(\nabla\)(ナブラ)について。

太郎くん
太郎くん

なんか不穏な感じしかしない。
\(\nabla\)の記号は見たことあるけど。

フクさん
フクさん

そんな変なものじゃない。
ベクトルに対しての偏微分だ。

太郎くん
太郎くん

十分以上に変な話だーーーー!!

フクさん
フクさん

いやいや、偏微分はなんとなくわかったでしょ?!

太郎くん
太郎くん

まぁ、特定の軸に沿っての微分だから、むしろシンプルになったものだとは聞いた。

フクさん
フクさん

ベクトルに対する偏微分は、ベクトルの各要素に対して個別に微分をするってだけで、
これもやはり複雑な事象をシンプルにしてくれる考え方だ。

太郎くん
太郎くん

まぁ一個一個に着目すればシンプルになるってことか・・・。

∇の定義

フクさん
フクさん

とりあえず\(\nabla\)の定義を確認しておこう。

\(
\nabla=
\begin{bmatrix}
\partial / \partial x_1 \\
\vdots \\
\partial / \partial x_n \\
\end{bmatrix}
\)

フクさん
フクさん

これを\(x_1,\dots, x_n \)のベクトルを持った関数に適用すると以下になる。

\(
\nabla f(x_1,\dots, x_n)=
\begin{bmatrix}
\displaystyle \frac{\partial f(x_1,\dots, x_n)}{\partial x_1} \\
\vdots \\
\displaystyle \frac{\partial f(x_1,\dots, x_n)}{\partial x_n} \\
\end{bmatrix}=
\mathrm{grad} f
\)

太郎くん
太郎くん

ようわからんが、\(x_1,\dots, x_n \)の各要素単位で微分をしたいんだなーってことしかわからんな。

フクさん
フクさん

ここではその程度の認識でOKだ。
あと、勾配を示す、\(\mathrm{grad} f\)もあるから、この点も覚えておくと良いだろう。

まとめ

フクさん
フクさん

まとめだよ。

  • 二次形式の微分についての話へ突入。
  • ∇(ナブラ)について説明。
    • ベクトルに対しての偏微分。
    • 各要素に対しての微分を行うだけなので、複雑な概念ではない。

バックナンバーはこちら。

コメント

タイトルとURLをコピーしました