markdown
ODE の方法選択と次のトラックmd c506bab
lecture/math/differential-equations/laplace-transforms-series-and-modeling.lecture.n.md
Download PDF

ODE の方法選択ほうほうせんたくと次つぎのトラック

date2026-07-14document_iddoc_496f2f825a99d777606d1aab5ca84caadescriptionODEで求める成果と方程式の構造を先に判定し、明示解・局所級数・定性解析・数値解から次の数学トラックまでを選ぶ出口講義。prerequisitesステップ関数・デルタ分布・因果的畳み込み / べき級数解法と Frobenius 法 / 相平面と安定性 / 物理モデルとPDEへの橋渡しtype講義content_typelecturestatusactiverelateddata/lecture/math/differential-equations/analytic-numerical-and-qualitative-ode-methods.lecture.n.md / data/lecture/math/differential-equations/classifying-first-order-odes.lecture.n.md / data/lecture/math/differential-equations/first-order-systems-and-matrix-exponentials.lecture.n.md / data/lecture/math/differential-equations/linearization-and-eigenvalue-stability.lecture.n.md / data/lecture/math/differential-equations/step-functions-delta-functions-and-convolution.lecture.n.md / data/lecture/math/differential-equations/power-series-and-frobenius-methods.lecture.n.md / data/lecture/math/differential-equations/physical-models-and-transition-to-pdes.lecture.n.md / data/lecture/math/differential-equations/variation-of-parameters-and-wronskian.lecture.n.md / data/lecture/math/differential-equations/phase-planes-and-stability.lecture.n.md / data/lecture/math/differential-equations/beyond-euler-runge-kutta-and-implicit-methods.lecture.n.md / data/lecture/math/analysis/introduction-to-laplace-transform.lecture.n.md / data/lecture/math/analysis/introduction-to-fourier-transform.lecture.n.md / data/lecture/math/partial-differential-equations/introduction-to-pdes.lecture.n.md / data/lecture/math/partial-differential-equations/initial-and-boundary-value-problems.lecture.n.md / data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md / data/lecture/math/partial-differential-equations/separation-of-variables-and-fourier-series.lecture.n.md / data/lecture/math/partial-differential-equations/fourier-transforms-and-pdes.lecture.n.md / data/exercise/math/differential-equations/ode-method-selection.exercise.n.md
mathdifferential-equationstransformsmodelinglecture

1導入どうにゅう

この講義こうぎでは、公式こうしきの名前なまえから解法かいほうを選えらぶのではなく、まず「何なにを求もとめるか」と「方程式ほうていしきがどの構造こうぞうを持もつか」を確認かくにんし、その組合くみあわせから候補こうほを絞しぼることを説明せつめいする。

明示式めいじしきが必要ひつようなのか、局所解きょくしょかいを構成こうせいしたいのか、安定性あんていせいや長時間挙動ちょうじかんきょどうを知しりたいのか、数値すうちを得えたいのかで、適切てきせつな出口でぐちは変かわる。一ひとつの方法ほうほうが使つかえることと、それが最短さいたん・最良さいりょうであることも同おなじではない。

2選択前せんたくまえに固定こていする契約けいやく

  1. 未知関数みちかんすう、独立変数どくりつへんすう、領域りょういきを決きめる。
  2. 線型せんけいか非線型ひせんけいか、係数けいすうは定数ていすうか変数へんすうかを調しらべる。
  3. 右辺うへんや入力にゅうりょくの形かたち、初期しょき・境界きょうかいデータの位置いちを確認かくにんする。
  4. 明示解めいじかい、局所展開きょくしょてんかい、定性的性質ていせいてきせいしつ、数値近似すうちきんじのどれを成果せいかとするかを決きめる。
  5. 適用条件てきようじょうけんと解かいの意味いみを確認かくにんし、外はずれたときの次候補つぎこうほを用意よういする。

3目的もくてきと構造こうぞうから候補こうほを選えらぶ

目的もくてきと構造こうぞう第一候補だいいちこうほ適用条件てきようじょうけんと次候補つぎこうほ
変数分離形へんすうぶんりけいの一階方程式いっかいほうていしき分離積分ぶんりせきぶん割わる前まえに定数解ていすうかいを確認かくにんし、初期値しょきちを代入だいにゅうした後あとで解かいの最大存在区間さいだいそんざいくかんを調しらべる
一階線型方程式いっかいせんけいほうていしき積分因子せきぶんいんし標準形ひょうじゅんけい y′+p(x)y=q(x) に直なおし、p,q の連続区間れんぞくくかんごとに解とく
完全微分形かんぜんびぶんけいポテンシャルの復元ふくげんM,N が C1 級きゅうで、穴あなのない(単連結たんれんけつな)領域りょういきなら My=Nx を確認かくにんする。穴あなのある領域りょういきでは大域たいいきポテンシャルが存在そんざいしない場合ばあいがある
一階方程式いっかいほうていしきの型かたがまだ不明ふめい一階方程式いっかいほうていしきの分類表ぶんるいひょう変数分離へんすうぶんり・線型せんけい・完全形かんぜんけい・同次形どうじけい・Bernoulli 形けいを順じゅんに照合しょうごうする
定数係数線型ていすうけいすうせんけいのスカラー高階斉次方程式こうかいせいじほうていしき特性方程式とくせいほうていしき根こんの重複ちょうふくと実数じっすう・複素数ふくそすうの設定せっていを確認かくにんする
定数係数線型方程式ていすうけいすうせんけいほうていしきで、非斉次項ひせいじこうが多項式たこうしき・指数しすう・正弦余弦せいげんよげんなど微分びぶんで閉とじる有限和ゆうげんわ未定係数法みていけいすうほう共鳴きょうめいを補正ほせいする。一般いっぱんの右辺うへんなら定数変化法ていすうへんかほうや畳たたみ込こみを考かんがえる
定数係数線型連立系ていすうけいすうせんけいれんりつけい行列指数関数ぎょうれつしすうかんすう・Jordan 構造こうぞう斉次系せいじけいは eA(t-t0)、強制系きょうせいけいは定数変化公式ていすうへんかこうしきと行列畳ぎょうれつたたみ込こみを使つかう
片側時間かたがわじかんの定数係数線型初期値問題ていすうけいすうせんけいしょきちもんだいで、切替きりかえ・遅延ちえん・衝撃入力しょうげきにゅうりょくを含ふくむラプラス変換へんかん通常入力つうじょうにゅうりょくでは区分的連続性くぶんてきれんぞくせい・指数位数しすういすうと共通収束半平面きょうつうしゅうそくはんへいめんを確認かくにんし、デルタ入力にゅうりょくは分布ぶんぷとして扱あつかう
解析的係数かいせきてきけいすうを持もつ通常点つうじょうてんの局所解きょくしょかいべき級数法きゅうすうほう最高階係数さいこうかいけいすうが展開中心てんかいちゅうしんで非零ひれいであり、正規化後せいきかごの係数けいすうと右辺うへんが解析的かいせきてきであることを確認かくにんする
正則特異点せいそくとくいてんの局所解きょくしょかいFrobenius 法ほう指標方程式しひょうほうていしきの根こんと対数項たいすうこうの要否ようひを確認かくにんする。不規則特異点ふきそくとくいてんでは一般いっぱんに直接適用ちょくせつてきようできない
線型系せんけいけいの原点げんてんの安定性あんていせい固有値こゆうち・Jordan 構造こうぞうすべての実部じつぶが負ふなら大域漸近安定たいいきぜんきんあんていである。すべての実部じつぶが 0 以下いかで虚軸上きょじくじょうの block が大おおきさ 1 なら Lyapunov 安定あんてい、正実部せいじつぶまたは虚軸上きょじくじょうの非自明ひじめい block があれば不安定ふあんていである
C1 級非線型系きゅうひせんけいけいの平衡点近傍へいこうてんきんぼう線型化せんけいか・固有値こゆうちすべての実部じつぶが負ふなら局所漸近安定きょくしょぜんきんあんてい、正せいの実部じつぶが一ひとつでもあれば不安定ふあんていである。正せいの実部じつぶがなく実部じつぶ 0 を含ふくむときは、非線型項ひせんけいこうや別べつの道具どうぐが必要ひつようになる
2 次元自律系じげんじりつけいの軌道配置きどうはいち相平面そうへいめん平衡点へいこうてん、流ながれの向むき、軌道きどうの形かたちを明示解めいじかいなしに調しらべられる。大域的結論たいいきてきけつろんには追加理論ついかりろんが必要ひつようになることがある
明示式めいじしきが困難こんなんで数値すうちが必要ひつようRunge-Kutta 法ほう・陰解法いんかいほう評価区間ひょうかくかんでの解かいの存在そんざい・一意性いちいせいと、次数じすうに必要ひつような滑なめらかさを確認かくにんする。陽的方法ようてきほうほうでは安定性あんていせいのため極端きょくたんに小ちいさい刻きざみ幅はばが必要ひつようになる硬かたい問題もんだいstiff problemでは陰解法いんかいほうを検討けんとうし、各段階かくだんかいの代数方程式だいすうほうていしきも所定精度しょていせいどまで解とく

この表ひょうは排他的はいたてきな分類ぶんるいではない。同おなじ問題もんだいに複数ふくすうの方法ほうほうを使つかい、明示式めいじしきを数値解すうちかいの検証けんしょうに、定性解析ていせいかいせきを計算範囲けいさんはんいの設計せっけいに使つかうこともある。

既存きそんの方法選択講義ほうほうせんたくこうぎは、閉形式へいけいしきに固執こしつせず解析かいせき・数値すうち・定性ていせいのどの成果せいかを求もとめるかを決きめる。一階方程式いっかいほうていしきの分類講義ぶんるいこうぎは、その後あとに具体的ぐたいてきな型かたを判定はんていする。この出口講義でぐちこうぎは、それらを高階方程式こうかいほうていしき・連立系れんりつけい・変換へんかん・級数きゅうすう・PDE まで含ふくむ最終経路さいしゅうけいろへまとめる。

data/lecture/math/differential-equations/analytic-numerical-and-qualitative-ode-methods.lecture.n.md data/lecture/math/differential-equations/classifying-first-order-odes.lecture.n.md

線型化せんけいかが与あたえるのは平衡点近傍へいこうてんきんぼうの局所結論きょくしょけつろんである。遠方えんぽうの初期値しょきち、大域吸引性たいいききゅういんせい、有限時間発散ゆうげんじかんはっさんには別べつの解析かいせきが必要ひつようになる。また、局所誤差きょくしょごさ・大域誤差たいいきごさ・絶対安定性領域ぜったいあんていせいりょういきは別々べつべつの概念がいねんであり、刻きざみ幅制御はばせいぎょだけで自動的じどうてきにすべてが保証ほしょうされるわけではない。

4ラプラス変換へんかんが閉とじる例れい

ここでは H(t-π) を、t<π で 0、t>π で 1 となる Heaviside 階段関数かいだんかんすうとする。t=π での値あたいは以下いかの連続解れんぞくかいに影響えいきょうしない。初期値問題しょきちもんだい

y′'+y=H(t-π),y(0)=0,y′(0)=0

を考かんがえる。片側かたがわラプラス変換へんかんを Y(s)=[PARSE ERROR: Undefined("Command(\"mathcal\")")]L{y}(s) とする。ここでは y,y′ が各有限区間かくゆうげんくかんで絶対連続ぜったいれんぞくで、必要ひつような導関数どうかんすうと入力にゅうりょくが区分的連続くぶんてきれんぞくかつ指数位数しすういすうを持もち、共通きょうつうの収束半平面しゅうそくはんへいめんで計算けいさんできるとする。このとき

[PARSE ERROR: Undefined("Command(\"mathcal\")")]L{y′}=sY-y(0),[PARSE ERROR: Undefined("Command(\"mathcal\")")]L{y′'}=s2Y-sy(0)-y′(0)

および [PARSE ERROR: Undefined("Command(\"mathcal\")")]L{H(t-π)}=e-πs/s より

Y(s)=e-πss(s2+1)=e-πs(1s-ss2+1).

遅延公式ちえんこうしきで逆変換ぎゃくへんかんすると

y(t)=H(t-π)(1-cos(t-π))

を得える。y と y′ は t=π で連続れんぞくであり、t<π と t>π の各区間かくくかんで方程式ほうていしきを直接確認ちょくせつかくにんできる。デルタ入力にゅうりょくは通常関数つうじょうかんすうについての上うえの仮定かていの外そとにあり、分布ぶんぷとしての拡張かくちょうを使つかう。

data/lecture/math/differential-equations/step-functions-delta-functions-and-convolution.lecture.n.md data/lecture/math/analysis/introduction-to-laplace-transform.lecture.n.md

5級数法きゅうすうほうが閉とじる例れい

y′'+xy′+y=0

は係数けいすうが全平面ぜんへいめんで解析的かいせきてきであり、すべての有限点ゆうげんてんが通常点つうじょうてんである。y=∑n=0∞anxn を代入だいにゅうして同おなじ冪べきを比較ひかくすると、n[PARSE ERROR: Undefined("Command(\"ge\")")]0 で

(n+2)(n+1)an+2+(n+1)an=0,an+2=-ann+2

を得える。a0,a1 を自由じゆうに選えらぶと偶数列ぐうすうれつと奇数列きすうれつがそれぞれ決きまり、二階線型方程式にかいせんけいほうていしきの 2 自由度じゆうどが現あらわれる。零れいでない偶数ぐうすう・奇数部分級数きすうぶぶんきゅうすうでは、2 項こうおきの項こうの絶対値ぜったいちの比ひは

|an+2xn+2anxn|=|x|2n+2⟶0

となる。比ひの判定はんていによりすべての x で収束しゅうそくするので、この例れいの級数解きゅうすうかいの収束半径しゅうそくはんけいは無限大むげんだいである。

展開中心てんかいちゅうしんが通常点つうじょうてんならべき級数きゅうすう、特異点とくいてんのうち正則特異点せいそくとくいてんなら Frobenius 法ほうが候補こうほになる。不規則特異点ふきそくとくいてんではこの二択にたくを機械的きかいてきに使つかえない。

data/lecture/math/differential-equations/power-series-and-frobenius-methods.lecture.n.md

6Fourier 級数きゅうすう・Fourier 変換へんかん・PDE を分わける

  • 有界区間ゆうかいくかんの線型境界値問題せんけいきょうかいちもんだいで、境界条件きょうかいじょうけんを満みたす固有関数こゆうかんすうが互たがいに直交ちょっこうし、対象たいしょうの関数かんすうを展開てんかいできるなら、固有関数展開こゆうかんすうてんかいが候補こうほになる。たとえば標準的ひょうじゅんてきな二階作用素にかいさようそ -X′' と X(0)=X(L)=0 では、固有関数こゆうかんすう sin(nπx/L) による正弦級数せいげんきゅうすうが現あらわれる。非斉次境界条件ひせいじきょうかいじょうけんは、まず既知関数きちかんすうを引ひいて斉次化せいじかすることがある。
  • 全空間ぜんくうかんの定数係数問題ていすうけいすうもんだいでは Fourier 変換へんかんが候補こうほになる。入口講義いりぐちこうぎでは、まず何回なんかいでも微分びぶんでき、関数かんすうとすべての導関数どうかんすうが無限遠むげんえんでどの多項式たこうしきの逆数ぎゃくすうよりも速はやく 0 へ近ちかづく Schwartz 関数かんすうを安全あんぜんな適用範囲てきようはんいとする。
  • 空間分布くうかんぶんぷを連続場れんぞくばとして残のこすモデルでは、局所保存則きょくしょほぞんそくと構成則こうせいそくから PDE を導みちびく。方程式ほうていしきだけでなく、領域りょういき・初期境界しょききょうかいデータ・解かいのクラスまで指定していして問題もんだいになる。

非線型問題ひせんけいもんだいでも Fourier 表現ひょうげんを使つかうことはあるが、積せきが周波数間しゅうはすうかんの畳たたみ込こみとなり、成分せいぶんは独立どくりつに対角化たいかくかされない。不規則領域ふきそくりょういきでは単純たんじゅんな三角関数さんかくかんすうではなく、領域りょういきと境界条件きょうかいじょうけんに応おうじた固有関数こゆうかんすうや数値的方法すうちてきほうほうが必要ひつようになる。

集中定数しゅうちゅうていすう ODE と連続場れんぞくば PDE のモデル選択せんたくは物理ぶつりモデルの橋渡はしわたし講義こうぎが担当たんとうする。ここでは導出どうしゅつを繰くり返かえさず、PDE トラックの履修順りしゅうじゅんだけを示しめす。

data/lecture/math/differential-equations/physical-models-and-transition-to-pdes.lecture.n.md data/lecture/math/partial-differential-equations/introduction-to-pdes.lecture.n.md data/lecture/math/partial-differential-equations/initial-and-boundary-value-problems.lecture.n.md data/lecture/math/partial-differential-equations/classification-of-second-order-linear-pdes.lecture.n.md data/lecture/math/vector-calculus/laplacian-and-physical-applications.lecture.n.md data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md

上うえの入口いりぐちと代表三方程式だいひょうさんほうていしきを終おえた後あとに、有界区間ゆうかいくかんの変数分離へんすうぶんりと全空間ぜんくうかんの Fourier 変換へんかんへ進すすむ。

data/lecture/math/analysis/introduction-to-fourier-transform.lecture.n.md data/lecture/math/partial-differential-equations/separation-of-variables-and-fourier-series.lecture.n.md data/lecture/math/partial-differential-equations/fourier-transforms-and-pdes.lecture.n.md

7明示解以外めいじかいいがいの出口でぐち

定性解析ていせいかいせきは「式しきで解とけないから妥協だきょうする」方法ほうほうではない。安定性あんていせい、周期軌道しゅうききどう、長時間挙動ちょうじかんきょどうのように、明示式めいじしきとは別べつの問といへ直接答ちょくせつこたえる。数値解法すうちかいほうも値あたいを出だすだけではなく、局所誤差きょくしょごさ・大域誤差たいいきごさ・安定性領域あんていせいりょういきを管理かんりして初はじめて方法ほうほうになる。

data/lecture/math/differential-equations/phase-planes-and-stability.lecture.n.md data/lecture/math/differential-equations/beyond-euler-runge-kutta-and-implicit-methods.lecture.n.md

8演習えんしゅうと確認かくにん

横断的おうだんてきな方法選択ほうほうせんたくは、次つぎの専用演習せんようえんしゅうで「成果せいか・構造こうぞう・第一候補だいいちこうほ・仮定かてい・次候補つぎこうほ」を一組ひとくみとして確認かくにんする。

data/exercise/math/differential-equations/ode-method-selection.exercise.n.md
raw .n.md をコピー
loc をコピー (filepath:line ~ line)
copy share link
copy encoded share link
path をコピー
copy share link
copy encoded share link
copy share link
copy encoded share link
タブを全て閉じる