応用解析 ・ 第 2 回 求積法 2 dy = ay dx y(x0 ) = y0 (1) (2) という問題を解いた。(1) を満たす関数は y(x) = Ceax , で、さらに (2) を満たす関数は y(x) = y0 ea(x−x0 ) C∈R であることがわかった。 ... C のこと ... この C を含んだ解全体(無限個の解) ... C を定めた(一つの)解 ... 一般解に含まれる積分定数をどのように選んでも表現できない解 注: m 階の常微分方程式の場合、 回積分して解くので、m 個の積分定数が出てきて、それを定めるのに 個の条件が必要になる。 以下は、次の常微分方程式の解き方(求積法)を学ぶ: • 変数分離形 • 同次形 • 1 階線形 • 2 階線形(定数係数) 2.1 変数分離形 (separation of variables) dy = f (x)g(y) dx 次のように解く。 1 dy g(y) dx dG(y(x)) dx = f (x) (g(y) 6= 0 のとき) = f (x) (関数 1/g の原始関数を G(y) = ∫ 最後に、g(y) 6= 0 と仮定したことによりのぞかれた解がないか調べる。 1 1/g(y) dy とおいた) 例. バクテリアの増殖をより正確に表すロジスティック方程式を解く。 dy = (a − by)y dx y 6= a/b と y 6= 0 のとき、以下の式は同値である。 dy 1 = (a − by)y dx ( ) 1 1 b dy − = a y by − a dx ∫ ∫ ) b 1( 1 dy − dy = a y by − a ) 1( log |y| − log |by − a| = a y = log by − a y = by − a y(1 − bC 00 eax ) = y(x) = • y 6= 0 の仮定について: (* 変数を分離する) 1 1 ∫ (* それぞれの辺を積分する) dx x+C ax + C 0 (* 変形する) 0 0 ±eax+C = ±eC eax = C 00 eax −aC 00 eax aC 00 , 00 bC − e−ax (C 00 ∈ R, C 00 6= 0) C 00 6= 0 00 という定数関数が方程式をみたす。この解は、上式において C = とおけば実現される。 00 という定数関数が方程式をみたす。この解は C → • y 6= a/b の仮定について:同様に、 とすれば実現される。 時刻 x = 0 でバクテリアの量は y(0) = y0 だったとする。この初期条件を用いて、一般解の中の積分定数 C 00 を 定め、特殊解を求めることができる。つまり、 一般解に代入すれば、 y(x) = a b + ( ya0 − b)e−ax という解を得る。 シンプルなモデル du/dt = au の解と比べるために、それぞれのモデルの解のグラフを描いてみるとよい。 2 2.2 同次形 (homogeneous) (y) dy =f dx x 解き方: 1. 新しい未知関数を z(x) = により導入する。 2. z が満たす微分方程式を導く。 y(x) = z(x) · x であるから、 y 0 (x) z0 z 0 (x) · x + z(x) = f (z(x)) f (z) − z = x = 3. 得られた方程式を z について解く(変数分離形になっているので解ける)。 4. z(x) に y(x)/x を戻し、y(x) を x を用いて表す。 5. ステップ 3. で方程式を解くときに f (z) − z で割るので、f (z) − z = 0 の根に関係する解を調べる。つまり、 ζ が f (z) − z = 0 の根ならば、 y(x) = も解であることに注意する。 例. y0 = ( y )2 x + 2y x 1. y = xz とおく 2. y 0 = z + xz 0 なので、z に対する方程式は z + xz 0 = z 2 + 2z ⇒ z0 = z2 + z x 3. 得られた方程式を z について解く z = C 0x z+1 (C 0 = ±eC , C 0 ∈ R, C 0 6= 0) 4. z を y/x で置き換え、さらに書き直せば、 y/x = C 0x y/x + 1 ⇒ y(x) = C 0 x2 1 − C 0x 5. f (z) − z = z 2 − z = 0 を満たすのは z = 0 と z = −1 であり、これはそれぞれ y(x) = 0 と y(x) = −x とい う解に対応しているが、これらは一般解で C 0 = 0, C 0 = ∞ とおいて得られる特殊解である。 3
© Copyright 2024 ExpyDoc