物理のかぎしっぽ 記事ソース/グリーン関数と逆行列

記事ソース/グリーン関数と逆行列

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

コンバート

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

公開・更新メニュー ▼▲

記事ソースの内容

============================================================
グリーン関数と逆行列
============================================================

グリーン関数と逆行列の類似性について書きます。



グリーン関数
======================

演算子 $L$ (ただし $L$ は、 $x$ についての演算子)
と既知関数 $f(x)$ 、未知関数 $\psi(x)$ 
があり以下の関係を満たすとします。

<tex>
L\psi(x)=f(x) \tag{##}
</tex>

ここで、グリーン関数 $G(x-x^\prime)$ を
次の性質を持つ関数として、
定義します [*]_ 。

<tex>
LG(x-x^\prime)=-\delta(x-x^\prime) \tag{##}
</tex>

.. [*] その求め方はここでは書きません。普通の演算子 $L$ に対してそんな関数が存在するということだけを知っておいてください。

これを使って、形式的に次のように書きます。

<tex>
G(x-x^\prime)=- L^{-1}\delta(x-x^\prime) \tag{##}
</tex>

すると、結局未知関数 $\psi(x)$ は次のように求まります。

<tex>
\psi(x) 
&= - \int L^{-1} \{ \delta(x-x^\prime) \} f(x^\prime) dx^\prime \\
&= \int G(x-x^\prime) f(x^\prime) dx^\prime \tag{##}
</tex>

ここで、式 $(1)$ を確認しておきましょう。
式 $(4)$ の両辺に $L$ を作用させます。
すると、 $L$ は $x$ のみに関わるので、

<tex>
L \psi(x) 
&= - \int L \{ L^{-1} \{ \delta(x-x^\prime) \} \} f(x^\prime) dx^\prime \\
&= \int \delta(x-x^\prime) f(x^\prime) dx^\prime \\
&= f(x) \tag{##}
</tex>

となり、確かに式 $(1)$ が成立していることが分かります。

逆行列
===================

一方、有限次元の逆行列を持つ行列 $A$ 、
既知ベクトル $\bm{f}$ 、
未知ベクトル $\bm{x}$ について、
次のような方程式を考えます。

<tex>
A\bm{x}=\bm{f} \tag{##}
</tex>

ここで、グリーン関数の行列版とでも言えるような、
次のベクトル群 $\bm{g}_i$ を導入します。

<tex>
A\bm{g}_1
=
\begin{pmatrix}
1 \\
0 \\
\vdots
\end{pmatrix} 
\equiv \bm{d}_1
\tag{##}
</tex>

<tex>
A\bm{g}_2
=
\begin{pmatrix}
0 \\
1 \\
\vdots
\end{pmatrix} 
\equiv \bm{d}_2
\tag{##} 
</tex>

すると、

<tex>
\bm{f} =\sum_i f_i A \bm{g}_i \tag{##}
</tex>

ですから、

<tex>
\bm{x} &= A^{-1} \bm{f} \\ 
&= A^{-1} \sum_i f_i A \bm{g}_i \\
&= \sum_i f_i \bm{g}_i \\
&= \sum_i \bm{g}_i f_i \tag{##}
</tex>

と表せます。ここで、 $f_i$ はベクトル $\bm{f}$ の第 $i$ 成分
です。
ここで、比較のため、式 $(4)$ をもう一度書きます。

<tex>
\psi(x) = \int G(x-x^\prime) f(x^\prime) dx^\prime \tag{4}
</tex>


きれいに二つが対応しているのが見てとれますね。
というのは、つまり、δ関数 $\bm{d}_i$ に逆作用素(逆行列) $A^{-1}$ を掛けると、
グリーン関数 $\bm{g}_i$ になります。一方、非斉次項の $f(x)$ に対応するのが、 $f_i$ なわけです。

それでは、今日はこの辺で。

@@author:クロメル@@
@@accept:2010-04-20@@
@@category:物理数学@@
@@id:greenFunctionAndMatrix@@
トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新の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.009 sec.