院試hub

1階常微分方程式変数分離・線形・完全微分形の解き方と解の一意性

1階常微分方程式 y′=f(x,y)y' = f(x, y) は、式の型で解き方が決まる。変数分離形は両辺を積分し、線形方程式 y′+P(x)y=Q(x)y' + P(x)y = Q(x) は積分因子 e∫P dxe^{\int P\,dx} を掛け、完全微分形 M dx+N dy=0M\,dx + N\,dy = 0(My=NxM_y = N_x)はポテンシャル関数を求める。初期値問題の解がただ1つに決まるのは、ff が連続で yy についてリプシッツ条件を満たすときである。

微分方程式 の第1単元(全4単元)。

院試数学 基礎完成シリーズ 微分積分編

基礎ノートと類題を答案の形で解いた問題集(2,480円)。収録章と見本ページを商品ページで確認できる。過去問の本文は収録していない。

微分積分編の内容を見る

常微分方程式・一般解・初期値問題

定義

常微分方程式と解

1変数の未知関数 y(x)y(x) とその導関数を含む等式を常微分方程式といい、含まれる最高階の導関数の階数を方程式の階数という。1階の方程式を y′=f(x,y)y' = f(x, y) の形に書いたものを正規形という。

区間 II で微分可能な関数 y=φ(x)y = \varphi(x) が、II のすべての xx で φ′(x)=f(x,φ(x))\varphi'(x) = f(x, \varphi(x)) を満たすとき、φ\varphi を II 上の解という。任意定数を1つ含み、それを動かして得られる解の族を一般解、定数に値を入れた解を特殊解という。一般解の定数をどう選んでも得られない解を特異解という。

定義

初期値問題

方程式 y′=f(x,y)y' = f(x, y) と条件 y(x0)=y0y(x_0) = y_0 の組を初期値問題という。一般解の任意定数を初期条件で決めると答えになる。解は x0x_0 を含む区間でだけ定義されることがあり、その区間も答えの一部である。

解の存在と一意性の条件

定理・法則

解の存在と一意性(ピカール・リンデレーフの定理)

条件

f(x,y)f(x, y) は長方形 D={∣x−x0∣≤a, ∣y−y0∣≤b}D = \{\lvert x - x_0\rvert \le a, \ \lvert y - y_0\rvert \le b\} で連続で、yy についてリプシッツ条件を満たす。つまり定数 KK があって

∣f(x,y1)−f(x,y2)∣≤K∣y1−y2∣((x,y1),(x,y2)∈D)\lvert f(x, y_1) - f(x, y_2)\rvert \le K\lvert y_1 - y_2\rvert \quad ((x, y_1), (x, y_2) \in D)

が成り立つ。M=max⁡D∣f∣M = \max_D\lvert f\rvert、h=min⁡(a,b/M)h = \min(a, b/M) とする。

主張

初期値問題 y′=f(x,y)y' = f(x, y)、y(x0)=y0y(x_0) = y_0 の解が区間 ∣x−x0∣≤h\lvert x - x_0\rvert \le h にただ1つ存在する。

要点

逐次近似 yk+1(x)=y0+∫x0xf(t,yk(t)) dty_{k+1}(x) = y_0 + \int_{x_0}^{x} f(t, y_k(t))\,dt が一様収束することで示す。∂f/∂y\partial f/\partial y が DD で連続なら、平均値の定理からリプシッツ条件が成り立つ。ff が連続なだけなら解は存在する(ペアノの定理)が、一意とは限らない。

計算例

一意性が崩れる例

y′=3y2/3y' = 3y^{2/3}、y(0)=0y(0) = 0 の解を調べよ。

解答

y=0y = 0 は解である。y=x3y = x^3 も y′=3x2=3(x3)2/3y' = 3x^2 = 3(x^3)^{2/3} を満たし、y(0)=0y(0) = 0 なので解である。さらに任意の c≥0c \ge 0 について、x≤cx \le c で 00、x>cx > c で (x−c)3(x - c)^3 とつないだ関数も解になる。f(y)=3y2/3f(y) = 3y^{2/3} は連続だが、∂f/∂y=2y−1/3\partial f/\partial y = 2y^{-1/3} が y=0y = 0 の近くで有界でなく、リプシッツ条件を満たさない。

検算

y0≠0y_0 \neq 0 の初期値なら、y0y_0 の近くで ∂f/∂y\partial f/\partial y が連続なので解は局所的にただ1つである。

計算例

解が有限の xx で発散する例

y′=y2y' = y^2、y(0)=1y(0) = 1 を解け。

解答

変数を分離して dyy2=dx\frac{dy}{y^2} = dx、積分して −1y=x+C-\frac{1}{y} = x + C。y(0)=1y(0) = 1 から C=−1C = -1 で、y=11−xy = \frac{1}{1 - x}。解は x<1x < 1 でだけ定義され、x→1−0x \to 1 - 0 で +∞+\infty に発散する。

検算

y′=1(1−x)2=y2y' = \frac{1}{(1 - x)^2} = y^2。f(y)=y2f(y) = y^2 は全平面でなめらかでも、解が全区間で存在するとは限らない。一意性の定理が保証するのは x0x_0 の近くの区間だけである。

型ごとの解き方

1階常微分方程式の型と解法

変数分離形
y′=g(x)h(y) ⇒ ∫dyh(y)=∫g(x) dx+Cy' = g(x)h(y) \ \Rightarrow\ \int\frac{dy}{h(y)} = \int g(x)\,dx + C

h(y∗)=0h(y_*) = 0 となる定数 y=y∗y = y_* も解で、割り算で失われるので別に書く。

線形
y′+P(x)y=Q(x) ⇒ y=e−∫P dx(∫Q(x)e∫P dx dx+C)y' + P(x)y = Q(x) \ \Rightarrow\ y = e^{-\int P\,dx}\left(\int Q(x)e^{\int P\,dx}\,dx + C\right)

両辺に積分因子 e∫P dxe^{\int P\,dx} を掛けると、左辺が (ye∫P dx)′\left(ye^{\int P\,dx}\right)' になる。

同次形
y′=g ⁣(yx) ⇒ u=yx,xu′=g(u)−uy' = g\!\left(\frac{y}{x}\right) \ \Rightarrow\ u = \frac{y}{x}, \quad xu' = g(u) - u

uu についての変数分離形になる。

ベルヌーイ型
y′+P(x)y=Q(x)yn (n≠0,1) ⇒ z=y1−n,z′+(1−n)Pz=(1−n)Qy' + P(x)y = Q(x)y^n \ (n \neq 0, 1) \ \Rightarrow\ z = y^{1-n}, \quad z' + (1 - n)P z = (1 - n)Q

zz についての線形方程式になる。n>0n > 0 なら y=0y = 0 も解。

完全微分形
M dx+N dy=0,∂M∂y=∂N∂x ⇒ Φ(x,y)=C,Φx=M, Φy=NM\,dx + N\,dy = 0, \quad \frac{\partial M}{\partial y} = \frac{\partial N}{\partial x} \ \Rightarrow\ \Phi(x, y) = C, \quad \Phi_x = M, \ \Phi_y = N
積分因子
My−NxN=ρ(x)⇒μ=e∫ρ dx,Nx−MyM=σ(y)⇒μ=e∫σ dy\frac{M_y - N_x}{N} = \rho(x) \Rightarrow \mu = e^{\int\rho\,dx}, \qquad \frac{N_x - M_y}{M} = \sigma(y) \Rightarrow \mu = e^{\int\sigma\,dy}

μM dx+μN dy=0\mu M\,dx + \mu N\,dy = 0 が完全微分形になる。

定理・法則

完全微分形の判定

条件

M,NM, N は単連結な領域 DD(穴のない領域)で連続な1階偏導関数をもつ。

主張

M dx+N dyM\,dx + N\,dy がある関数 Φ\Phi の全微分 dΦd\Phi になるための必要十分条件は、DD で My=NxM_y = N_x が成り立つことである。このとき Φ(x,y)=∫(x0,y0)(x,y)(M dx+N dy)\Phi(x, y) = \int_{(x_0, y_0)}^{(x, y)}(M\,dx + N\,dy) は経路によらない。

要点

単連結でない領域では十分条件にならない。M=−yx2+y2M = \frac{-y}{x^2 + y^2}、N=xx2+y2N = \frac{x}{x^2 + y^2} は原点を除いた平面で My=NxM_y = N_x を満たすが、原点を回る閉曲線に沿った積分は 2π2\pi で、全平面からの1価のポテンシャルをもたない。

線形・完全微分形・ベルヌーイ型の計算例

計算例

線形方程式

y′+2xy=xy' + 2xy = x、y(0)=0y(0) = 0 を解け。

解答

積分因子は e∫2x dx=ex2e^{\int 2x\,dx} = e^{x^2}。両辺に掛けると

(yex2)′=xex2⟹yex2=12ex2+C\left(ye^{x^2}\right)' = xe^{x^2} \quad\Longrightarrow\quad ye^{x^2} = \frac{1}{2}e^{x^2} + C

よって y=12+Ce−x2y = \frac{1}{2} + Ce^{-x^2}。y(0)=0y(0) = 0 から C=−12C = -\frac{1}{2} で、y=12(1−e−x2)y = \frac{1}{2}\left(1 - e^{-x^2}\right)。

検算

y′=xe−x2y' = xe^{-x^2}、2xy=x−xe−x22xy = x - xe^{-x^2} で、和は xx。変数分離形 y′=x(1−2y)y' = x(1 - 2y) とみて解いても同じ式になる。

計算例

完全微分形

(2xy+3) dx+(x2+4y) dy=0(2xy + 3)\,dx + (x^2 + 4y)\,dy = 0 の一般解を求めよ。

解答

M=2xy+3M = 2xy + 3、N=x2+4yN = x^2 + 4y で My=2x=NxM_y = 2x = N_x なので完全微分形である。Φx=M\Phi_x = M を xx で積分して Φ=x2y+3x+h(y)\Phi = x^2y + 3x + h(y)。Φy=x2+h′(y)=N\Phi_y = x^2 + h'(y) = N から h′(y)=4yh'(y) = 4y、h=2y2h = 2y^2。

x2y+3x+2y2=Cx^2y + 3x + 2y^2 = C

検算

Φ=x2y+3x+2y2\Phi = x^2y + 3x + 2y^2 を偏微分すると Φx=2xy+3=M\Phi_x = 2xy + 3 = M、Φy=x2+4y=N\Phi_y = x^2 + 4y = N。

計算例

積分因子で完全微分形にする

y dx+2x dy=0y\,dx + 2x\,dy = 0(x>0x > 0)を解け。

解答

My=1M_y = 1、Nx=2N_x = 2 で完全微分形でない。Nx−MyM=1y\frac{N_x - M_y}{M} = \frac{1}{y} は yy だけの関数なので、μ=e∫dy/y=y\mu = e^{\int dy/y} = y。掛けると y2 dx+2xy dy=d(xy2)=0y^2\,dx + 2xy\,dy = d(xy^2) = 0 で、一般解は xy2=Cxy^2 = C。

検算

変数を分離すると dyy=−dx2x\frac{dy}{y} = -\frac{dx}{2x} から y=Cx−1/2y = Cx^{-1/2} となり、xy2=C2xy^2 = C^2 と同じ曲線族になる。

計算例

ベルヌーイ型

y′+y=y2y' + y = y^2 の一般解を求めよ。

解答

n=2n = 2 なので z=y−1z = y^{-1} とおく。z′=−y−2y′=−y−2(y2−y)=−1+zz' = -y^{-2}y' = -y^{-2}(y^2 - y) = -1 + z から z′−z=−1z' - z = -1。積分因子 e−xe^{-x} で解いて z=1+Cexz = 1 + Ce^{x}。よって y=11+Cexy = \frac{1}{1 + Ce^{x}}。z=y−1z = y^{-1} の置き換えで除いた y=0y = 0 も解である。

検算

y=11+Cexy = \frac{1}{1 + Ce^x} について y′=−Cex(1+Cex)2y' = -\frac{Ce^x}{(1 + Ce^x)^2}、y2−y=1−(1+Cex)(1+Cex)2=−Cex(1+Cex)2y^2 - y = \frac{1 - (1 + Ce^x)}{(1 + Ce^x)^2} = -\frac{Ce^x}{(1 + Ce^x)^2} で一致する。C=0C = 0 は定数解 y=1y = 1 にあたる。

答案でよくある誤り

  1. 1

    h(y)h(y) で割って変数分離し、h(y)=0h(y) = 0 の定数解を落とす。y′=xyy' = xy では y=0y = 0 も解で、一般解 y=Cex2/2y = Ce^{x^2/2} の C=0C = 0 に含めるか、別に書く。

  2. 2

    ∫dyy=log⁡y\int\frac{dy}{y} = \log y と書き、y<0y < 0 の解を失う。log⁡∣y∣\log\lvert y\rvert と書き、∣y∣=eCeG(x)\lvert y\rvert = e^{C}e^{G(x)} から y=±eCeG(x)y = \pm e^{C}e^{G(x)} を任意定数 AA にまとめる。

  3. 3

    線形方程式で y′+P(x)y=Q(x)y' + P(x)y = Q(x) の形に直す前に積分因子を作る。xy′+y=x2x y' + y = x^2 なら xx で割って P=1xP = \frac{1}{x} とする(x≠0x \neq 0 の区間で解く)。

  4. 4

    Φx=M\Phi_x = M を xx で積分したあと、積分定数を yy の関数 h(y)h(y) にしない。h(y)h(y) は Φy=N\Phi_y = N から決める。

  5. 5

    解の定義区間を書かない。y′=y2y' = y^2、y(0)=1y(0) = 1 の解 y=11−xy = \frac{1}{1 - x} は x<1x < 1 でだけ解であり、x>1x > 1 の枝は同じ初期値問題の解ではない。

微分方程式の前後の単元

このページの本文を数式のLaTeXのまま読むには Markdown版。最終更新 2026-11-16。誤りに気づいたら 運営者 へ連絡してほしい。

院試数学 基礎完成シリーズ 微分積分編

基礎ノートと類題を答案の形で解いた問題集(2,480円)。収録章と見本ページを商品ページで確認できる。過去問の本文は収録していない。

微分積分編の内容を見る