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,NC1 きゅうで、あなのない(単連結たんれんけつな)領域りょういきなら 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-π))

る。yyt=π連続れんぞくであり、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=0anxn代入だいにゅうしておなべき比較ひかくすると、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+20

となる。判定はんていによりすべての 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
タブを全て閉じる