markdown
初期値問題と境界値問題md d61bd6b
lecture/math/differential-equations/initial-and-boundary-value-problems.lecture.n.md
Download PDF

初期値問題しょきちもんだい境界値問題きょうかいちもんだい

mathdifferential-equationsinitial-value-problemboundary-value-problemlecture

1導入どうにゅう

この講義こうぎでは、微分方程式びぶんほうていしき条件じょうけんわせて解集合かいしゅうごう決定けっていする方法ほうほう説明せつめいする。微分方程式びぶんほうていしき積分せきぶんすると、ふつう任意定数にんいていすうふく関数族かんすうぞくられる。条件じょうけんは、その関数族かんすうぞく制約せいやくすが、条件じょうけん整合性せいごうせいによってはかい存在そんざいしない場合ばあい複数ふくすうかいのこ場合ばあいもある。

たとえば

y=2x

積分せきぶんすると、

y=x2+C

である。このままでは Cまらない。y(0)=1 という条件じょうけんくわわると C=1 になり、かいy=x2+1まる。

2初期値問題しょきちもんだい

初期値問題しょきちもんだいinitial value problemとは、ある基準点きじゅんてん未知関数みちかんすうやその導関数どうかんすうあたい指定していする問題もんだいである。たとえば

y=f(x,y),y(x0)=y0

一階いっかい初期値問題しょきちもんだいである。

二階にかいなら、たとえば

y'=F(x,y,y),y(x0)=y0,y(x0)=v0

のように、関数かんすうあたい速度そくどにあたる導関数どうかんすうあたい指定していすることがおおい。

時間発展じかんはってんあらわ問題もんだいでは、初期値しょきち現在げんざい状態じょうたい指定していし、そこから将来しょうらい状態じょうたい決定けっていする。ただし、独立変数どくりつへんすう時間じかんでない場合ばあいにも、おな数学的形式すうがくてきけいしきもちいる。

3境界値問題きょうかいちもんだい

境界値問題きょうかいちもんだいboundary value problemとは、区間くかん端点たんてん領域りょういき境界きょうかい条件じょうけん指定していする問題もんだいである。たとえば

y'=0,y(0)=0,y(1)=1

境界値問題きょうかいちもんだいである。

境界値問題きょうかいちもんだいかんがえる理由りゆうは、ぼう温度分布おんどぶんぷげんかたち定常状態ていじょうじょうたいのように、両端りょうたん境界きょうかい条件じょうけん全体ぜんたいかたち拘束こうそくする現象げんしょうがあるからである。

4命題めいだい1:条件じょうけん任意定数にんいていすう制約せいやくする

仮定かていとして、微分方程式びぶんほうていしき一般解いっぱんかいが 1 任意定数にんいていすう Cふく

y=ϕ(x,C)

あたえられているとする。条件じょうけん y(x0)=y0すと、許容きょようされる C

ϕ(x0,C)=y0

たすあたい制約せいやくされる。この方程式ほうていしき一意解いちいかいつときにかぎり、C一意いちい決定けっていされる。

4.1証明しょうめい

かいy=ϕ(x,C)けるという仮定かていから、基準点きじゅんてん x0 でのあたいϕ(x0,C) である。条件じょうけん y(x0)=y0たすには、

ϕ(x0,C)=y0

必要ひつようである。ぎゃくに、この方程式ほうていしきたす Cえらべば、

y(x0)=ϕ(x0,C)=y0

となる。したがって条件じょうけん任意定数にんいていすう制約せいやくする方程式ほうていしきとして機能きのうする。この方程式ほうていしきかいがなければ、提示ていじされた一般解族いっぱんかいぞくには条件じょうけんたす要素ようそがない。複数ふくすうかいがあれば、複数ふくすう候補こうほのこる。

この命題めいだいは、一般解族いっぱんかいぞく条件じょうけん操作そうさを、任意定数にんいていすうたいする方程式ほうていしき帰着きちゃくできることをしめす。

5れい1:一階いっかい初期値問題しょきちもんだい

y=-y,y(0)=1

かんがえる。まず方程式ほうていしきだけをくと、

y=Ce-x

である。ここで y(0)=1使つかうと、

1=Ce0=C

だから

y=e-x

である。

検算けんざんする。y=e-x なら y=-e-x=-y であり、また y(0)=1 である。したがって方程式ほうていしき条件じょうけん両方りょうほうたす。

6れい2:二階にかい境界値問題きょうかいちもんだい

y'=0,y(0)=0,y(1)=1

かんがえる。y'=0 を 1 かい積分せきぶんすると y=A、さらに積分せきぶんすると

y=Ax+B

である。境界条件きょうかいじょうけん y(0)=0 から B=0y(1)=1 から A=1 である。したがって

y=x

である。

初期値問題しょきちもんだいでも境界値問題きょうかいちもんだいでも、基本的きほんてき操作そうさおなじである。まず一般解いっぱんかいもとめ、条件じょうけん代入だいにゅうして任意定数にんいていすう制約せいやくする。

n 任意定数にんいていすうがある場合ばあいでも、n 条件じょうけんがあれば自動的じどうてきまるわけではない。一般いっぱんには、条件じょうけん代入だいにゅうによって任意定数にんいていすうかんする方程式系ほうていしきけいられ、その可解性かかいせい一意性いちいせい確認かくにんする必要ひつようがある。線型同次せんけいどうじ ODE の解空間かいくうかん線型独立せんけいどくりつかい基底解きていかいbasis solutionという。線型せんけい ODE の一般解いっぱんかい同次部分どうじぶぶん基底解きていかい線型結合せんけいけつごうあらわし、線型せんけい条件じょうけん場合ばあいには、この方程式系ほうていしきけい連立一次方程式れんりついちじほうていしきとなる。その係数行列けいすうぎょうれつ正則せいそくであれば、任意定数にんいていすう一意いちい決定けっていされる。

n 階線型同次かいせんけいどうじ ODE の基底解きていかいy1,,yn とする。同一どういつてん x0y,y,,y(n-1)あたい指定していする標準的ひょうじゅんてき初期条件しょきじょうけんでは、積分定数せきぶんていすうかんする係数行列けいすうぎょうれつ

(yj(i-1)(x0))i,j=1n

である。この行列式ぎょうれつしきx0 における Wronskian W(y1,,yn)(x0) という。Wronskian が 0 でなければ、この行列ぎょうれつ正則せいそくであり、積分定数せきぶんていすう一意いちい決定けっていされる。一方いっぽうことなる端点たんてん境界条件きょうかいじょうけん一般いっぱん線型条件せんけいじょうけんでは、別途べっと構成こうせいされる条件行列じょうけんぎょうれつ正則性せいそくせい判定はんていする必要ひつようがあり、その行列式ぎょうれつしき一般いっぱんには Wronskian ではない。

data/lecture/math/linear-operator/affine-spaces-and-condition-determination.lecture.n.md data/lecture/math/differential-equations/variation-of-parameters-and-wronskian.lecture.n.md

7初期値問題しょきちもんだい境界値問題きょうかいちもんだいちが

項目こうもく初期値問題しょきちもんだい境界値問題きょうかいちもんだい
条件じょうけん位置いち1 てんあたい導関数どうかんすう区間くかん端点たんてん境界きょうかい
典型的てんけいてき解釈かいしゃく指定していされた状態じょうたいから発展はってん決定けっていする境界きょうかいのデータによって分布ぶんぷ拘束こうそくする
代表例だいひょうれい人口増加じんこうぞうか冷却れいきゃく運動うんどう定常温度ていじょうおんどげん固有振動こゆうしんどう
注意点ちゅういてん存在そんざい一意性いちいせいには仮定かてい必要ひつようかいがない、または複数ふくすうある場合ばあいがある

ここで固有振動こゆうしんどうという言葉ことばてくるが、このページではくわしくあつかわない。境界条件きょうかいじょうけんたす特別とくべつ振動しんどうだけがゆるされる問題もんだいがあり、そのさき固有値こゆうち固有関数こゆうかんすうまなぶ。

8注意ちゅうい条件じょうけん本数ほんすうだけでは十分じゅうぶんでない

二階にかい微分方程式びぶんほうていしきでは、一般解いっぱんかい通常つうじょう 2 任意定数にんいていすうあらわれる。そのため 2 条件じょうけん必要ひつようになることがおおい。しかし、条件じょうけんが 2 あればかならず 1 ほんかい存在そんざいする、というわけではない。

たとえば

y'=0

かんがえる。1 かい積分せきぶんすると y=A、もう 1 かい積分せきぶんすると

y=Ax+B

である。ここで境界きょうかいでの導関数どうかんすう指定していする条件じょうけん

y(0)=0,y(1)=1

す。y=A だから、y(0)=0 から A=0 である。一方いっぽうy(1)=1 から A=1 でなければならない。これは矛盾むじゅんする。したがって、この境界値問題きょうかいちもんだいにはかいがない。

一方いっぽう境界条件きょうかいじょうけん

y(0)=0,y(1)=0

なら、A=0 だけがまり、B任意にんいである。したがって

y=B

がすべてかいであり、かい無限個むげんこある。

このように、条件じょうけん本数ほんすう階数かいすうっていても、境界条件きょうかいじょうけん位置いち方程式ほうていしき性質せいしつによって、かいがない場合ばあい複数ふくすうある場合ばあいこる。

この問題もんだいは、存在そんざい一意性いちいせい理論りろんにおける主要しゅよう分析対象ぶんせきたいしょうである。

data/lecture/math/differential-equations/existence-uniqueness-and-lipschitz-condition.lecture.n.md

9かた手順てじゅん

  1. 未知関数みちかんすう独立変数どくりつへんすう確認かくにんする。
  2. 階数かいすう確認かくにんする。
  3. 条件じょうけんおなてんあつまっているか、端点たんてん境界きょうかいかれているかを確認かくにんする。
  4. 一般解いっぱんかいられる場合ばあいは、それをもとめる。
  5. 条件じょうけん代入だいにゅうして任意定数にんいていすうたいする制約方程式せいやくほうていしきき、許容きょようされる定数ていすう決定けっていする。
  6. 最後さいごもと方程式ほうていしき条件じょうけん代入だいにゅうして検算けんざんする。

10つぎ講義こうぎ

つぎに、一階微分方程式いっかいびぶんほうていしき方程式ほうていしき形式けいしきおうじて分類ぶんるいし、適用可能てきようかのう解法かいほう判定はんていする。

data/lecture/math/differential-equations/classifying-first-order-odes.lecture.n.md

11演習えんしゅうリンク

data/exercise/math/differential-equations/introduction-to-odes-and-direct-integration.exercise.n.md data/exercise/math/differential-equations/existence-uniqueness-and-numerical-methods.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
タブを全て閉じる