markdown
一般二階線型常微分方程式の基礎理論md ab8c2a1
lecture/math/differential-equations/general-second-order-linear-odes-overview.lecture.n.md
Download PDF

一般二階線型常微分方程式いっぱんにかいせんけいじょうびぶんほうていしきの基礎理論きそりろん

mathdifferential-equationssecond-orderlecture

1導入どうにゅう

この講義こうぎでは、二階線型にかいせんけいという語ごが定数係数ていすうけいすうを意味いみしないことを明確めいかくにし、一般論いっぱんろんと特殊解法とくしゅかいほうを分離ぶんりする。

2対象たいしょうとする方程式ほうていしき

一般いっぱんの二階線型にかいせんけいの標準形ひょうじゅんけいは

y′'+P(x)y′+Q(x)y=R(x)

である。I を区間くかん、x0∈I とし、P,Q,R:I→K は連続れんぞくであると仮定かていする。ただし、ここで K は R または C を表あらわす。この条件じょうけんのもとで、任意にんいの a,b∈K に対たいし、初期条件しょきじょうけん y(x0)=a,y′(x0)=b を満みたす解かいが区間くかん I 全体ぜんたいで一意いちいに存在そんざいする。

3一般理論いっぱんりろんと特殊解法とくしゅかいほうの区別くべつ

特性方程式とくせいほうていしきは二階線型にかいせんけいの一般解法いっぱんかいほうではない。係数けいすうが定数ていすうであり、erx の微分びぶんで形かたちが保存ほぞんされる場合ばあいに適用てきようできる特殊とくしゅな解法かいほうである。この区別くべつをあらかじめ明確めいかくにしなければ、変数係数へんすうけいすうの方程式ほうていしきに誤あやまって特性方程式とくせいほうていしきを適用てきようすることになる。

4線型理論せんけいりろんの骨格こっかく

一般いっぱんの二階線型にかいせんけいでは、まず

L[y]=y′'+P(x)y′+Q(x)y

と定義ていぎする。ここで L は関数かんすうを関数かんすうへ送おくる線型作用素せんけいさようそlinear operatorである。実際じっさい、u,v を 2 回かい微分可能びぶんかのうな関数かんすう、α,β を定数ていすうとすると、微分びぶんの線型性せんけいせいより

L[αu+βv]=αL[u]+βL[v]

である。

data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md

同次方程式どうじほうていしき

y′'+P(x)y′+Q(x)y=0

の解集合かいしゅうごうは、加法かほうと定数倍ていすうばいで閉とじるため線型空間せんけいくうかんになる。ただし、閉とじていることだけでは 2 次元じげんとは言いえない。2 次元じげんであることは、初期値問題しょきちもんだいの存在そんざい・一意性いちいせいから従したがう。

基準点きじゅんてん x0 を 1 つ固定こていし、同次どうじ解かい y に

Φ(y)=(y(x0),y′(x0))

を対応たいおうさせる。Φ は線型写像せんけいしゃぞうである。係数けいすう P,Q が連続れんぞくであるというこのページの仮定かていのもとでは、任意にんいの初期値しょきち (a,b) に対たいして

y(x0)=a,y′(x0)=b

を満みたす同次どうじ方程式ほうていしきの解かいが一意いちいに存在そんざいする。したがって Φ は同次どうじ解空間かいくうかんから K2 への同型どうけいであり、同次どうじ解空間かいくうかんは 2 次元じげんである。

したがって独立どくりつな 2 解かい y1,y2 が得えられれば、

yh=C1y1+C2y2

が同次解どうじかいの一般形いっぱんけいである。

非同次方程式ひどうじほうていしきでは、特殊解とくしゅかい yp を 1 つ得えれば

y=yh+yp

が一般解いっぱんかいである。この分解ぶんかいは線型性せんけいせいに由来ゆらいする。

厳密げんみつには、非ひ同次どうじ方程式ほうていしきは

L[y]=R

である。L[yp]=R を満みたす特殊解とくしゅかい yp が 1 つあるとする。h∈kerL なら

L[yp+h]=L[yp]+L[h]=R

なので yp+h は解かいである。逆ぎゃくに y も L[y]=R の解かいなら

L[y-yp]=L[y]-L[yp]=R-R=0

であるから、y-yp∈kerL である。したがって非ひ同次どうじ解かい全体ぜんたいは

yp+kerL

である。ここでの kerL が、同次どうじ方程式ほうていしきの解空間かいくうかんである。

5一次連立系いちじれんりつけいとしての観点かんてん

一般いっぱんの二階線型にかいせんけいは、一階連立系いっかいれんりつけいへ変換へんかんすると理論りろんの構造こうぞうが明確めいかくになる。x1=y,x2=y′ と設定せっていすると、

(x1x2)′=(01-Q(x)-P(x))(x1x2)+(0R(x))

を得える。係数けいすうが連続れんぞくなら、この一階連立系いっかいれんりつけいの標準理論ひょうじゅんりろんにより初期値しょきち y(x0),y′(x0) から解かいが一意いちいに決定けっていされる。二階にかいで初期条件しょきじょうけんが 2 本ほん必要ひつようになる理由りゆうは、状態じょうたいが y と y′ の 2 成分せいぶんを持もつからである。

この観点かんてんは、定数係数ていすうけいすうの特性方程式とくせいほうていしきと連立系れんりつけいの固有値こゆうちが同おなじ情報じょうほうを表あらわしていることも説明せつめいする。定数係数ていすうけいすうなら行列ぎょうれつが定数ていすうとなり、固有値こゆうちにより指数関数しすうかんすうの時間発展じかんはってんへ分解ぶんかいできる。

一次連立系いちじれんりつけいへの変換へんかん、行列指数関数ぎょうれつしすうかんすう、固有値こゆうちとの対応たいおうは、後続講義こうぞくこうぎで詳述しょうじゅつする。

data/lecture/math/differential-equations/first-order-systems-and-matrix-exponentials.lecture.n.md

6Wronskian と独立性どくりつせい

独立どくりつな 2 解かい y1,y2 を確認かくにんするとき、Wronskian

W(y1,y2)=y1y2′-y1′y2

を用もちいる。W(x0)≠0 なら、係数けいすうが連続れんぞくな区間くかんでは y1,y2 は基本解系きほんかいけいを構成こうせいする。実際じっさい、Abel の恒等式こうとうしき

W(x)=W(x0)exp(-∫x0xP(s)ds)

により、Wronskian は区間内くかんないの他ほかの点てんでも 0 にならない。逆ぎゃくに、2 解かいが線型従属せんけいじゅうぞくなら Wronskian は恒等的こうとうてきに 0 である。したがって Wronskian は 2 解かいの線型独立性せんけいどくりつせいを判定はんていし、同次解空間どうじかいくうかんが 2 次元じげんであるという既知きちの事実じじつと併あわせて、その 2 解かいが基本解系きほんかいけいを構成こうせいすることを確認かくにんする。

7解法かいほうの選択基準せんたくきじゅん

状況じょうきょう候補こうほとなる方法ほうほう理由りゆう
定数係数ていすうけいすう特性方程式とくせいほうていしきerx の形かたちが保存ほぞんされる
定数係数ていすうけいすうの非同次方程式ひどうじほうていしきで、外力項がいりょくこうを含ふくむ有限次元関数空間ゆうげんじげんかんすうくうかんが微分びぶんで不変ふへん未定係数法みていけいすうほう有限次元ゆうげんじげんの試行空間しこうくうかんを選択せんたくできる。共鳴時きょうめいじには同次解どうじかいとの重複ちょうふくを除のぞくよう試行空間しこうくうかんを調整ちょうせいする
基本解系きほんかいけいが既知きちで外力項がいりょくこうが連続れんぞく定数変化法ていすうへんかほう同次解どうじかいから特殊解とくしゅかいを構成こうせいできる
変数係数へんすうけいすうで初等解しょとうかいが困難こんなん級数解法きゅうすうかいほう係数比較けいすうひかくで局所解きょくしょかいを構成こうせいする

8具体例ぐたいれい

y′'+xy′+y=0

は二階線型にかいせんけいであるが、定数係数ていすうけいすうではない。したがって r2+xr+1=0 のような特性方程式とくせいほうていしきを作成さくせいしてはならない。x が残存ざんぞんするため、erx の代入だいにゅうは代数方程式だいすうほうていしきを生成せいせいしない。

一方いっぽう、y′'-3y′+2y=0 では係数けいすうが定数ていすうである。erx を代入だいにゅうすると r2-3r+2=0 を得える。この対比たいひにより、特性方程式法とくせいほうていしきほうの本質ほんしつが二階にかいではなく定数係数ていすうけいすうにあることが判明はんめいする。

9適用範囲てきようはんいと限界げんかい

ここでの一般論いっぱんろんは線型性せんけいせいと、区間全体くかんぜんたいにおける正規化後せいきかごの係数けいすう P,Q および外力項がいりょくこう R の連続性れんぞくせいに基礎きそを置おく。正規化前せいきかまえの方程式ほうていしきで最高階係数さいこうかいけいすうが 0 になる場合ばあいには、その係数けいすうが 0 でない区間くかんへ制限せいげんしてから標準形ひょうじゅんけいへ変形へんけいする必要ひつようがある。非線型ひせんけいの二階方程式にかいほうていしきでは、解空間かいくうかんの線型構造せんけいこうぞうも重合じゅうごうの原理げんりも一般いっぱんには成立せいりつしない。

10演習えんしゅうリンク

data/exercise/math/differential-equations/second-order-linear-constant-coefficient-odes.exercise.n.md data/exercise/math/differential-equations/second-order-linear-odes.exercise.n.md

11関連かんれんリンク

data/lecture/math/differential-equations/second-order-linear-constant-coefficient-odes.lecture.n.md data/lecture/math/differential-equations/variation-of-parameters-and-wronskian.lecture.n.md data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md

12前まえの講義こうぎ

data/lecture/math/differential-equations/second-order-linear-constant-coefficient-odes.lecture.n.md

13次つぎの講義こうぎ

data/lecture/math/differential-equations/extensions-of-second-order-linear-odes.lecture.n.md

Foundations of General Second-Order Linear Ordinary Differential Equations

1Introduction

This lecture distinguishes general second-order linear equations from their constant-coefficient special case and separates structural results from specialized solution methods.

2Equations Under Consideration

The standard form of a general second-order linear equation is

y′'+P(x)y′+Q(x)y=R(x).

Let I be an interval, let x0∈I, and assume that P,Q,R:I→K are continuous, where K denotes either R or C. For every a,b∈K, these assumptions guarantee a unique solution on all of I satisfying y(x0)=a and y′(x0)=b.

3General Theory and Specialized Methods

The characteristic equation is not a general method for every second-order linear equation. It applies when the coefficients are constant, because differentiation preserves the form of erx. This distinction prevents the erroneous application of a characteristic polynomial to a variable-coefficient equation.

4Structure of the Linear Theory

Define

L[y]=y′'+P(x)y′+Q(x)y.

The operator L maps functions to functions and is linear: for twice differentiable functions u,v and constants α,β∈K,

L[αu+βv]=αL[u]+βL[v].
data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md

The solution set of the homogeneous equation L[y]=0 is a vector space. Closure under addition and scalar multiplication alone does not establish its dimension. Fixing x0∈I, define

Φ(y)=(y(x0),y′(x0)).

Existence and uniqueness for the initial-value problem show that Φ:kerL→K2 is a linear isomorphism. Consequently, kerL is two-dimensional. If y1,y2 are linearly independent solutions, every homogeneous solution has the form

yh=C1y1+C2y2.

For the nonhomogeneous equation L[y]=R, suppose that one particular solution yp is known. If h∈kerL, then

L[yp+h]=L[yp]+L[h]=R.

Conversely, if L[y]=R, then L[y-yp]=0. Hence the full solution set is the affine space

yp+kerL.

This identity is the precise meaning of the familiar decomposition y=yh+yp.

5First-Order System Formulation

Set x1=y and x2=y′. The equation is equivalent to

(x1x2)′=(01-Q(x)-P(x))(x1x2)+(0R(x)).

Continuity of the coefficient matrix and forcing term yields existence and uniqueness on I by the standard theory of linear first-order systems. Two initial conditions are required because the state has the two components y and y′. In the constant-coefficient case, the matrix is constant, and its eigenvalues encode the same information as the roots of the characteristic polynomial.

The subsequent lecture on first-order systems develops this conversion, the matrix exponential, and the correspondence with eigenvalues in detail.

data/lecture/math/differential-equations/first-order-systems-and-matrix-exponentials.lecture.n.md

6Wronskian and Linear Independence

For two homogeneous solutions, define the Wronskian by

W(y1,y2)(x)=y1(x)y2′(x)-y1′(x)y2(x).

If W(y1,y2)(x0)≠0 at one point of I, then y1,y2 form a fundamental system on I. Indeed, Abel's identity gives

W(x)=W(x0)exp(-∫x0xP(s)ds),

so the Wronskian cannot vanish elsewhere on the interval. Conversely, linear dependence forces the Wronskian to vanish identically. Thus the Wronskian tests whether two solutions are linearly independent; together with the previously established two-dimensionality of the homogeneous solution space, this test determines whether they form a fundamental system.

7Criteria for Selecting a Solution Method

SituationCandidate methodRationale
Constant coefficientsCharacteristic-equation methodDifferentiation preserves the form of erx
Constant-coefficient nonhomogeneous equation for which a finite-dimensional function space containing the forcing term is invariant under differentiationMethod of undetermined coefficientsA finite-dimensional trial space is available; under resonance, it must be adjusted to remove overlap with the homogeneous solutions
General continuous forcing with a known fundamental systemVariation of parametersConstructs a particular solution from the homogeneous solutions
Variable coefficients for which elementary solutions are unavailablePower-series or Frobenius methodConstructs local solutions through coefficient comparison

8Examples

The equation

y′'+xy′+y=0

is second-order and linear but does not have constant coefficients. Substitution of erx leaves dependence on x and therefore does not produce an algebraic equation; an expression such as r2+xr+1=0 is not a characteristic equation.

By contrast, y′'-3y′+2y=0 has constant coefficients. Substitution of erx gives r2-3r+2=0. Thus, applicability of the characteristic-equation method depends on constant coefficients, not merely on the order of the equation.

9Scope and Limitations

The vector-space and superposition results depend on linearity. They do not generally hold for nonlinear second-order equations. The global assertion on I also depends on continuity of the normalized coefficients P,Q and forcing term R throughout that interval; if the leading coefficient of an unnormalized equation vanishes, the equation must first be restricted to intervals on which normalization is valid.

10Exercises

data/exercise/math/differential-equations/second-order-linear-constant-coefficient-odes.exercise.n.md data/exercise/math/differential-equations/second-order-linear-odes.exercise.n.md

11Related Lectures

data/lecture/math/differential-equations/variation-of-parameters-and-wronskian.lecture.n.md data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md

12Previous Lecture

data/lecture/math/differential-equations/second-order-linear-constant-coefficient-odes.lecture.n.md

13Next Lecture

data/lecture/math/differential-equations/extensions-of-second-order-linear-odes.lecture.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
タブを全て閉じる