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=0、y(1)=1 から A=1 である。したがって

y=x

である。

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

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

n 階線型同次かいせんけいどうじ ODE の基底解きていかいを y1,…,yn とする。同一どういつの点てん x0 で y,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
タブを全て閉じる