先週は

非斉次方程式の解+斉次方程式の解=非斉次方程式の解

非斉次方程式
(An(x)(ddx)n+An1(x)(ddx)n1++A1(x)ddx+A0(x))y=C(x)
と、上の式でC(x)=0とした斉次方程式
(An(x)(ddx)n+An1(x)(ddx)n1++A1(x)ddx+A0(x))y=0
を考える。非斉次方程式の解としてy1(x)を1つ、斉次方程式の解をとしてy0(x)を1つ、それぞれ見つけたとする。y0(x)+y1(x)もまた、非斉次方程式の解である。

というところまでやったので、今週はまず、これの簡単な例をやってみよう。

ddxy=x+y

という非斉次線型微分方程式を解きたい。これは「変数分離できる形」にはなってない。

まずは試行錯誤で解を探す。たとえばy=ax+bが解になるだろうか、ということを考えてみる。代入してみると、

addxy=x+ax+by0=(1+a)x+ba

となるから、a=1,b=1が解となる。これでy=x1という解が見つかったわけである。ここで「バンザイ、解が見つかった」と終わってはいけない。なぜならこの解は「特別なある1つの解」であって、全ての解を求めていないのである。

関数y=x1は上のグラフであり、この線の上という(全x-y平面から見たらほんとに狭い)範囲の上での「解」を求めたに過ぎない。前に述べたように、このような解を「特解」と呼ぶ。我々が求めたいのは全x-y平面を埋め尽くす、「一般解」である。

非斉次になっているのはxという項のせいだから、これを消して

ddxy=y

という斉次方程式を作ってみると、この解は

y=Cex

である。非斉次方程式の解は特解にこの「斉次方程式の一般解」を足せば作ることができる。すなわち、

y=x1+Cex

が「一般解」なのである。

それで「一般解」としてOKなのはどうしてですか?
前に話した、「n階微分方程式の階はちょうどn個の未定パラメータを含む」というのが効いてきます。つまりここでCという未定のパラメータ(積分定数ですが)があるので、これで全ての解が表せてます。「全て」というところを納得するには、x=0の時にy=C1になる、ということを考えるといいです。特解y=x1の方はx=0の時y=1と固定されていまっているけど、「一般解」の方はy=C1Cは任意なので、どのようなyの値に対しても対応する解がちゃんとあります。
下のグラフで言うと、y=x1+Cexは、グラフの全平面を覆い尽くすことができる、ということです。

グラフを上の図に示した。グラフはC0.5ずつ変えた線を示しているが、もちろん線と線の隙間にもちゃんと線があり、全平面を埋め尽くしている。どのような初期値(x0,y0)から出発しても、この微分方程式に従うその後の変化がわかることになる。

重ねあわせの原理のおかげで

(斉次方程式の一般解)+(非斉次方程式の特解)=(斉次方程式の一般解)

という関係が成立するおかげで、このような計算ができる。

ここでやったことは以下のように考えてもよい。まず特解y=x1を見つけたから、「実際の解は特解に近い形をしているだろう」と推測し、「とりあえず特解に未知の関数Yを足したものが解だろう」とあたりをつけて、y=x1+Yと置いてみる。これを元の微分方程式に代入すれば、

ddx(x1+Y)=xx1+Yy1+ddxY=1+Y

となるから、後はddxY=Yという微分方程式を解けばよい。