バックナンバーはこちら。
https://www.simulationroom999.com/blog/compare-matlabpythonscilabjulia2-backnumber/
はじめに
正規方程式を導出するまでの説明。
今回は二次形式の微分(勾配)を実際の多項式に適用してみる。
登場人物
博識フクロウのフクさん
![指差しフクロウ](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/指差しフクロウ.png)
イラストACにて公開の「kino_k」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=iKciwKA9&area=1
エンジニア歴8年の太郎くん
![技術者太郎](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/技術者01アップ.png)
イラストACにて公開の「しのみ」さんのイラストを使用しています。
https://www.ac-illust.com/main/profile.php?id=uCKphAW2&area=1
ロードマップ【再掲】
![太郎くん](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/技術者02アップ.png)
まずはロードマップの再掲で現在位置の確認。
二次形式の微分のところだね。
![正規方程式に至る道、二次形式、対称行列、二次形式の微分、グラム行列、二乗和誤差、正規方程式](https://www.simulationroom999.com/blog/wp-content/uploads/2023/03/02_正規方程式に至る道-1024x313.png)
![フクさん](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/指差しフクロウ.png)
今回は、具体的な多項式に対して、二次形式の微分を求める。
使用する多項式は以下としよう。
\(
f(x,y)=3x^2+2y^2+5xy
\)
∇による偏導関数
![フクさん](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/指差しフクロウ.png)
まずは\(\nabla\)による偏導関数。
\(
\nabla f(x,y)=
\begin{bmatrix}
\displaystyle\frac{\partial f(x,y)}{\partial x}\\
\displaystyle\frac{\partial f(x,y)}{\partial y}\\
\end{bmatrix}=
\begin{bmatrix}
6x+5y\\
5x+4y
\end{bmatrix}
\)
![フクさん](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/お休みフクロウ.png)
上記を元に以下が成立すればOK
\(
2AX=
\begin{bmatrix}
6x+5y\\
5x+4y
\end{bmatrix}
\)
普通に手計算
![太郎くん](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/技術者01アップ.png)
あとは、これを各ツール、各言語で確認するって感じか。
![フクさん](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/指差しフクロウ.png)
そうだね。
と言っても、実は手計算でも確認できちゃうけどね。
\(
A=
\begin{bmatrix}
3 & 5/2\\
5/2 & 2
\end{bmatrix}
\)
\(
2AX=
2
\begin{bmatrix}
3 & 5/2\\
5/2 & 2
\end{bmatrix}
\begin{bmatrix}
x\\
y
\end{bmatrix}=
\begin{bmatrix}
6x+5y\\
5x+4y
\end{bmatrix}
\)
![太郎くん](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/技術者03アップ.png)
つまり、ツールで計算させるまでもない!!
![フクさん](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/お休みフクロウ.png)
まぁ、ここはツールに慣れることを、目的とした作業になるかもね。
まとめ
![フクさん](https://www.simulationroom999.com/blog/wp-content/uploads/2020/05/指差しフクロウ.png)
まとめだよ。
- 具体的な二次形式の多項式に対して微分。
- ∇による微分結果確認。
- 二次形式の微分の公式による結果確認。
- ツールで計算させるまでもないが、一応やっておく。
バックナンバーはこちら。
コメント