物理のかぎしっぽ 記事ソース/多変数の微分関数の逆変換 のバックアップ(No.1)

記事ソース/多変数の微分関数の逆変換

これはrst2hooktailの記事ソース保存・変換用です(詳細).

コンバート

最近コンバートされた結果: HTMLPDFTeX

公開・更新メニュー ▼▲

記事ソースの内容

=========================================================================
多変数の微分関数の逆変換
=========================================================================

皆さんは、極座標つまり $(x,y)$ と $(r , \theta)$ の変換行列(ヤコビアン)を見て [*]_ 、

.. [*]  $r= \sqrt{x^2+y^2}, \theta = \tan^{-1}(y/x)$ ですね。

<tex>
\dfrac{dr}{dx} \neq (\dfrac{dx}{dr})^{-1} \tag{##}
</tex>

であることに戸惑った経験はありませんか?

そういうことができるのは、どんな時なのかということについて調べてみました。

具体例(2次元極座標)
=========================

<tex>
x = r \cos \theta, y = \sin \theta \tag{##}
</tex>

<tex>
r= \sqrt{x^2+y^2}, \theta = \tan^{-1}(y/x) \tag{##}
</tex>

の微分形式を考えてみます。

<tex>
\begin{pmatrix}
d r \\
d \theta
\end{pmatrix}
&=
\begin{pmatrix}
\dfrac{dr}{dx} & \dfrac{dr}{dy} \\
\dfrac{d \theta}{dx} & \dfrac{d \theta}{dy}
\end{pmatrix}
\begin{pmatrix}
d x \\
d y
\end{pmatrix} \\
&=
\begin{pmatrix}
\dfrac{x}{\sqrt{x^2+y^2}} & \dfrac{y}{\sqrt{x^2+y^2}} \\
\dfrac{-y}{x^2+y^2} & \dfrac{x}{x^2+y^2}
\end{pmatrix}
\begin{pmatrix}
d x \\
d y
\end{pmatrix} \\
&= 
\begin{pmatrix}
\cos \theta & \sin \theta \\
-\sin \theta /r & \cos \theta /r
\end{pmatrix}
\begin{pmatrix}
d x \\
d y
\end{pmatrix} \tag{##}
</tex>

そして逆変換は、

<tex>
\begin{pmatrix}
d x \\
d y
\end{pmatrix}
&=
\begin{pmatrix}
\dfrac{dx}{dr} & \dfrac{dx}{d \theta} \\
\dfrac{dy}{dr} & \dfrac{dy}{d \theta}
\end{pmatrix}
\begin{pmatrix}
d r \\
d \theta
\end{pmatrix} \\
&=
\begin{pmatrix}
\cos \theta & -r \sin \theta \\
\sin \theta & r \cos \theta
\end{pmatrix}
\begin{pmatrix}
d r \\
d \theta
\end{pmatrix} \tag{##}
</tex>

ですね。確かに、

<tex>
\dfrac{dr}{dx} = \cos \theta \neq (\cos \theta)^{-1} (\dfrac{dx}{dr})^{-1} \tag{##}
</tex>

<tex>
\dfrac{d \theta}{dy} = \cos \theta/r \neq (r \cos \theta)^{-1} = (\dfrac{dy}{d \theta })^{-1} \tag{##}
</tex>

と逆関数の微分法は成り立っていないようです。

一般論
==========================-

ここで、一般的な場合に拡張して関係を調べてみましょう。変数同士の変換行列のランクを落とさないと仮定して、
二対二組 $(a,b),(x,y)$ の変数間の変換を考えます。

<tex>
\begin{pmatrix}
da \\
db
\end{pmatrix}
=
\begin{pmatrix}
\dfrac{da}{dx} & \dfrac{da}{dy} \\
\dfrac{db}{dx} & \dfrac{db}{dy}
\end{pmatrix}
\begin{pmatrix}
dx \\
dy
\end{pmatrix} \tag{##}
</tex>

という関係が成立していたとすると、仮定より、
上式の行列は逆を持ちます。すると、

<tex>
\begin{pmatrix}
dx \\
dy
\end{pmatrix} 
&=
\dfrac{1}{\dfrac{da}{dx}\dfrac{db}{dy}-\dfrac{db}{dx}\dfrac{da}{dy}}
\begin{pmatrix}
\dfrac{db}{dy} & -\dfrac{da}{dy} \\
-\dfrac{db}{dx} & \dfrac{da}{dx}
\end{pmatrix}
\begin{pmatrix}
da \\
db
\end{pmatrix} \\
&\equiv J^{-1} \begin{pmatrix}
\dfrac{db}{dy} & -\dfrac{da}{dy} \\
-\dfrac{db}{dx} & \dfrac{da}{dx}
\end{pmatrix}
\begin{pmatrix}
da \\
db
\end{pmatrix} \tag{##}
</tex>

となります。ここで記号 $\equiv$ はこれでヤコビアン $J$ を定義するという意味です。

これと、

<tex>
\begin{pmatrix}
dx \\
dy
\end{pmatrix} 
=
\begin{pmatrix}
\dfrac{dx}{da} & \dfrac{dx}{db} \\
\dfrac{dy}{da} & \dfrac{dy}{db}
\end{pmatrix}
\begin{pmatrix}
da \\
db
\end{pmatrix} \tag{##}
</tex>

と比較します。すると、

<tex>
J^{-1} \begin{pmatrix}
\dfrac{db}{dy} & -\dfrac{da}{dy} \\
-\dfrac{db}{dx} & \dfrac{da}{dx}
\end{pmatrix}
&=
\begin{pmatrix}
\dfrac{dx}{da} & \dfrac{dx}{db} \\
\dfrac{dy}{da} & \dfrac{dy}{db}
\end{pmatrix} \tag{##}
</tex>

という関係が成立します。

再び極座標
================================

さて、得られた結果の検証をしてみましょう。 $(a,b,x,y) \to (r,\theta,x,y) $ とすると、

まず、

<tex>
J = 1/r(\cos^2 \theta + \sin^2 \theta) = 1/r \tag{##}
</tex>

となります。




@@author:クロメル@@
@@accept:2011-06-24@@
@@category:熱力学@@
@@id:invOfDifOfMPF@@
トップ   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Modified by 物理のかぎプロジェクト PukiWiki 1.4.6 Copyright © 2001-2005 PukiWiki Developers Team. License is GPL.
Based on "PukiWiki" 1.3 by yu-ji Powered by PHP 5.3.29 HTML convert time to 0.010 sec.