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

Initial Value Problems and Boundary Value Problems

1Introduction

This lecture explains how a differential equation and prescribed conditions jointly determine a solution set. Integrating a differential equation commonly produces a family of functions containing arbitrary constants. Conditions constrain that family, but incompatible or dependent conditions may leave no solution or multiple solutions.

For example, integrating

y′=2x

gives

y=x2+C.

The constant C is undetermined until a condition is imposed. The condition y(0)=1 yields C=1 and therefore selects the solution y=x2+1.

2Initial Value Problems

An initial value problem specifies the value of an unknown function and, when required, some of its derivatives at a common reference point. For example,

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

is a first-order initial value problem.

For a second-order equation, one commonly prescribes both the function and its first derivative:

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

When the independent variable represents time, these data may specify the current state from which an evolution is determined. The same mathematical form is also used when the independent variable is not time.

3Boundary Value Problems

A boundary value problem specifies conditions at the endpoints of an interval or on the boundary of a domain. For example,

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

is a boundary value problem.

Boundary value problems arise when endpoint or boundary data constrain an entire spatial profile, as in steady-state temperature distributions, the shapes of strings, and other equilibrium configurations.

4Proposition 1: Conditions Constrain Arbitrary Constants

Suppose that a general solution of a differential equation is given by

y=ϕ(x,C),

where C is an arbitrary constant. Imposing y(x0)=y0 restricts the admissible values of C to solutions of

ϕ(x0,C)=y0.

The constant C is determined uniquely only if this equation has a unique solution.

4.1Proof

The representation y=ϕ(x,C) gives y(x0)=ϕ(x0,C). Hence the condition y(x0)=y0 requires

ϕ(x0,C)=y0.

Conversely, any value of C satisfying this equation gives

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

Thus the condition becomes an equation constraining the arbitrary constant. If that equation has no solution, the conditioned problem has no solution within the stated general family; if it has multiple solutions, multiple candidates remain.

This proposition shows that imposing a condition on a general family can be reduced to solving an equation for its arbitrary constants.

5Example 1: A First-Order Initial Value Problem

Consider

y′=-y,y(0)=1.

Solving the differential equation alone gives

y=Ce-x.

The initial condition yields

1=Ce0=C,

so

y=e-x.

Verification is essential: y′=-e-x=-y and y(0)=1, so the function satisfies both the equation and the condition.

6Example 2: A Second-Order Boundary Value Problem

Consider

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

Integrating twice gives y′=A and then

y=Ax+B.

The boundary condition y(0)=0 gives B=0, while y(1)=1 gives A=1. Therefore,

y=x.

The same basic operation occurs in initial and boundary value problems: determine a general solution and substitute the conditions to constrain its arbitrary constants.

The presence of n arbitrary constants and n conditions does not by itself ensure unique determination. In general, substitution produces a system of equations for the constants, and its solvability and uniqueness must be checked. Linearly independent solutions that span the solution space of a homogeneous linear ODE are called basis solutions. When the homogeneous part of a general linear ODE solution is expressed as a linear combination of basis solutions and the imposed conditions are linear, the resulting system is linear. If its coefficient matrix is nonsingular, the arbitrary constants are uniquely determined.

Let y1,…,yn be a basis of solutions of an nth-order homogeneous linear ODE. For standard initial conditions specifying y,y′,…,y(n-1) at the same point x0, the coefficient matrix for the integration constants is

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

Its determinant is the Wronskian W(y1,…,yn)(x0). If the Wronskian is nonzero, the matrix is nonsingular and the integration constants are uniquely determined. For boundary conditions imposed at different endpoints or for general linear conditions, one must instead test the separately constructed condition matrix; its determinant is not generally a 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

7Comparison of Initial and Boundary Value Problems

ItemInitial value problemBoundary value problem
Location of dataFunction and derivative values at one reference pointData at endpoints or on a boundary
Typical interpretationDetermine an evolution from a specified stateConstrain a profile through boundary data
Representative applicationsPopulation growth, cooling, motionSteady temperature, strings, normal modes
Principal cautionExistence and uniqueness require hypothesesNo solution or multiple solutions may occur

The term normal mode appears in the table but is not developed here. In some problems, only particular oscillations satisfy the boundary conditions; their analysis leads to eigenvalues and eigenfunctions.

8Caution: The Number of Conditions Is Not Sufficient

A second-order differential equation commonly has a general solution containing two arbitrary constants, so two conditions are often required. Nevertheless, two conditions do not guarantee the existence of exactly one solution.

Consider

y′'=0.

Twice integrating gives

y=Ax+B,y′=A.

Now impose the derivative boundary conditions

y′(0)=0,y′(1)=1.

The first requires A=0, whereas the second requires A=1. These requirements are inconsistent, so the boundary value problem has no solution.

By contrast, if

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

then only A=0 is fixed and B remains arbitrary. Hence every constant function

y=B

is a solution, and the problem has infinitely many solutions.

Thus, even when the number of conditions matches the order, the placement of the conditions and the structure of the equation may lead to no solution or multiple solutions. Existence and uniqueness theory analyzes precisely these questions.

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

9Analysis Procedure

  1. Identify the unknown function and independent variable.
  2. Determine the order of the equation.
  3. Determine whether the conditions are imposed at one reference point or distributed across endpoints or a boundary.
  4. Determine a general solution when one is available.
  5. Substitute the conditions, solve the resulting constraint equations, and identify the admissible constants.
  6. Substitute the result into the original equation and all prescribed conditions to verify it.

10Next Lecture

The next lecture classifies first-order differential equations by form and explains how to identify applicable solution methods.

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

11Exercises

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
タブを全て閉じる