物理のかぎしっぽ 記事ソース/一階常微分方程式の基礎知識

記事ソース/一階常微分方程式の基礎知識

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

コンバート

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

公開・更新メニュー ▼▲

記事ソースの内容

============================================================
一階常微分方程式の基礎知識
============================================================

特定の形の常微分方程式を解く時に知っていると便利なことを書きます。

その形とは、 $t$ の関数について、

<tex>
\frac{df(t)}{dt}=g(t)f(t) \tag{##}
</tex>

というものです。

これは、両辺を $f(t)$ で割ってやると、

<tex>
\frac{\dot{f}}{f}=g(t) \tag{##}
</tex>

となり、 $C$ を積分定数として両辺の $t$ での不定積分を取って、

<tex>
ln|f|= \int g(t) dt + C \tag{##}
</tex>

よって、両辺を指数関数の肩に乗せれば、

<tex>
|f|=exp(\int g(t) dt +C)=C^\prime exp(\int g(t) dt) \tag{##}
</tex>

となります。

特に $g(t)$ が $t$ の冪関数の時、 $C$ を定数として、

<tex>
\dot{f} = C t^a f  \tag{##}
</tex>

この時には、

<tex>
f(t)=
\begin{cases}
exp(\dfrac{C t^{a+1}}{a+1}) \ \ \ \ (when \ \   a \neq -1)  \\
t^{C} \ \ \ \ (when \ \  a = -1)
\end{cases} \tag{##}
</tex>

となります。

ここで、 

<tex>
\lim_{a \to -1} t^{a+1} \to 1\ \ , \ \ \lim_{a \to -1 \pm 0} \frac{1}{a+1} = \pm \infty
</tex>

より、

<tex>
\lim_{a \to -1 \pm 0 }exp^{\frac{C t^{a+1}}{a+1}} = \infty \ \ \mathrm{or}\ \ 0 \tag{##}
</tex> 

となり、極限が $t^C$ にならないことに注意してください。

ちなみに、 $\dot{f}=(t^2 - 2t +3)f$ など、微分して一般の多項式が出てくるものが欲しければ、

例えば、この例だと、

<tex>
\dot{f} &= \frac{d}{dt} \left( exp(t^3/3) \times exp(-2t^2/2) \times exp(3t) \right) \\
&= \frac{d}{dt}exp(t^3/3-t^2+3t) \\
&= (t^2 -2t +3)f(t) \tag{##}
</tex>

の様に、積を作ることで和が作られます。

二階の微分方程式については、目下勉強中です (^^; 。

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

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