方向場ほうこうば・Euler 法ほう・誤差ごさと安定性あんていせいの入口いりぐち

date2026-05-26description方向場と Euler 法を、解析解が得られない場合の挙動把握・誤差評価・数値安定性の入口として整理する。prerequisites微分方程式の入口 / 初期値問題の存在・一意性type講義statusactive

mathdifferential-equationsnumerical-methodslecture

導入どうにゅう

このページの核心かくしんは、解析解かいせきかいが得えられない場合ばあいでも、方向場ほうこうばで解かいの傾向けいこうを把握はあくし、Euler 法ほうで近似解きんじかいを構成こうせいし、誤差ごさと安定性あんていせいを別概念べつがいねんとして確認かくにんすることである。

何なにを解とくページか

標準形ひょうじゅんけいは

y^{'} = f (x, y), y (x_{0}) = y_{0}

である。方向場ほうこうばDirection field は、各点かくてん $(x, y)$ に傾かたむき $f (x, y)$ を配置はいちした図式ずしきである。Euler 法ほうEuler method は、現在点げんざいてんの傾かたむきで短みじかい直線近似ちょくせんきんじを行おこなう数値解法すうちかいほうである。

なぜこの方針ほうしんを選えらぶのか

微分びぶんは局所的きょくしょてきな一次近似いちじきんじを与あたえる。したがって、現在点げんざいてんで $y^{'} = f (x, y)$ を評価ひょうかすれば、短みじかい区間くかんでは $y$ がどちらへ変化へんかするかを推定すいていできる。方向場ほうこうばは式しきを解とかずに全体傾向ぜんたいけいこうを確認かくにんする道具どうぐであり、Euler 法ほうはその局所近似きょくしょきんじを反復はんぷくする方法ほうほうである。

解法かいほうの流ながれ

初期値しょきち $y (x_{0}) = y_{0}$ と刻きざみ幅はば $h$ を設定せっていする。Euler 法ほうでは

x_{n + 1} = x_{n} + h, y_{n + 1} = y_{n} + h f (x_{n}, y_{n})

と更新こうしんする。これは Taylor 展開てんかい

y (x + h) = y (x) + h y^{'} (x) + O (h^{2})

の一次項いちじこうで切きる近似きんじである。したがって局所打切きょくしょうちきり誤差ごさは $O (h^{2})$ であり、有限区間ゆうげんくかんで反復はんぷくした大域誤差たいいきごさは標準的ひょうじゅんてきに $O (h)$ である。

具体例ぐたいれい

増加ぞうかする例れい

y^{'} = y, y (0) = 1

で $h = 0.1$ とすると、Euler 法ほうは $y_{n + 1} = 1.1 y_{n}$ を与あたえる。したがって $y_{1} = 1.1, y_{2} = 1.21$ である。正確解せいかくかい $e^{x}$ と比較ひかくすると、刻きざみ幅はばを縮小しゅくしょうすると近似きんじが改善かいぜんする。

安定性あんていせいが問題もんだいになる例れい

y^{'} = a y, a < 0

の真しんの解かいは減衰げんすいする。Euler 法ほうでは

y_{n + 1} = (1 + a h) y_{n}

である。数値解すうちかいが減衰げんすいするには $| 1 + a h | < 1$ が必要ひつようである。たとえば $a = - 10$ 、 $h = 0.3$ では $1 + a h = - 2$ となり、真しんの解かいは減衰げんすいするにもかかわらず数値解すうちかいは振動しんどうしながら増大ぞうだいする。

精度せいどと安定性あんていせいの区別くべつ

精度せいどは、刻きざみ幅はばを小ちいさくしたとき真しんの解かいへどの速はやさで接近せっきんするかを表あらわす。安定性あんていせいは、反復はんぷくが誤差ごさを増幅ぞうふくするか抑制よくせいするかを表あらわす。高精度こうせいどな方法ほうほうでも、問題もんだいと刻きざみ幅はばの組合くみあわせによって不安定ふあんていになることがある。

陰的方法いんてきほうほうへの接続せつぞく

Forward Euler 法ほうは現在点げんざいてんの傾かたむき $f (x_{n}, y_{n})$ を使用しようする。一方いっぽう、Backward Euler 法ほうは次点じてんの傾かたむき $f (x_{n + 1}, y_{n + 1})$ を使用しようし、

y_{n + 1} = y_{n} + h f (x_{n + 1}, y_{n + 1})

を解とく。未知量みちりょう $y_{n + 1}$ が右辺うへんにも出現しゅつげんするため、各段階かくだんかいで方程式ほうていしきを解とく計算負担けいさんふたんが生しょうじる。しかし減衰問題げんすいもんだいでは、安定性あんていせいが改善かいぜんされることがある。

モデル方程式ほうていしき $y^{'} = a y$ に適用てきようすると、Backward Euler 法ほうは

y_{n + 1} = \frac{1}{1 - a h} y_{n}

を与あたえる。 $a < 0$ なら $| 1 / (1 - a h) | < 1$ が広ひろい刻きざみ幅はばで成立せいりつする。この比較ひかくは、数値解法すうちかいほうの選択せんたくでは精度せいどだけでなく安定性あんていせいを確認かくにんする必要ひつようがあることを示しめす。

方向場ほうこうばで判定はんていできること

方向場ほうこうばは厳密解げんみつかいを生成せいせいしないが、解析前かいせきまえの診断しんだんに有効ゆうこうである。

確認項目かくにんこうもく	方向場ほうこうばから得えられる情報じょうほう
平衡解へいこうかい	傾かたむきが 0 になる水平線すいへいせん
増減ぞうげん	傾かたむきの符号ふごう
安定性あんていせい	近傍きんぼうの矢印やじるしが平衡解へいこうかいへ向むかうか
数値計算すうちけいさんの危険きけん	急激きゅうげきな変化へんかや傾かたむきの大おおきい領域りょういき

どこまで成なり立たつか

方向場ほうこうばは解曲線かいきょくせんの候補こうほを可視化かしかする道具どうぐであり、1 本ほんの厳密解げんみつかいを与あたえるものではない。Euler 法ほうは入口いりぐちとして重要じゅうようだが、実用計算じつようけいさんでは Runge-Kutta 法ほう、Backward Euler 法ほう、適応刻てきおうきざみ幅はばなどを検討けんとうする。

次つぎに参照さんしょうするページ

data/lecture/math/differential-equations/Euler法の次に何が来るか：Runge-Kuttaと陰的方法の入口-講義.n.md

演習えんしゅうリンク

data/exercise/math/differential-equations/存在一意性と数値解法-基本演習.n.md