跳到主要内容

東京大学 工学系研究科 2025年8月実施 数学 第1問

Author​

祭音Myyura (co-authored with GPT 5.6 SOL)

Description​

I​

次の微分方程式の一般解 y(x)y(x) を求めよ。

y′′+6y′+8y−3e−x=0.y''+6y'+8y-3e^{-x}=0.
x2y′′−2y−4x3=0.x^2y''-2y-4x^3=0.

II​

微分方程式

x(2x3−1)y′+2x2+y−2xy2=0x(2x^3-1)y'+2x^2+y-2xy^2=0

を考える。

  1. 特解が y=A0+A1x+A2x2y=A_0+A_1x+A_2x^2 の形で与えられるとき、A0,A1,A2A_0,A_1,A_2 を求めよ。
    1. の特解を用いて一般解を求めよ。

题目描述​

I. 分别求下列二阶线性常微分方程的一般解:

y′′+6y′+8y−3e−x=0,y''+6y'+8y-3e^{-x}=0,

以及

x2y′′−2y−4x3=0.x^2y''-2y-4x^3=0.

II. 对一阶非线性方程

x(2x3−1)y′+2x2+y−2xy2=0,x(2x^3-1)y'+2x^2+y-2xy^2=0,

先假设存在二次多项式特解 y=A0+A1x+A2x2y=A_0+A_1x+A_2x^2,求三个系数;再利用该特解求一般解。

Kai​

I.1​

同次方程式の特性方程式は

λ2+6λ+8=(λ+2)(λ+4)=0\lambda^2+6\lambda+8=(\lambda+2)(\lambda+4)=0

である。右辺を 3e−x3e^{-x} とみて特解を yp=Ae−xy_p=Ae^{-x} とおくと、

(1−6+8)Ae−x=3Ae−x=3e−x(1-6+8)Ae^{-x}=3Ae^{-x}=3e^{-x}

より A=1A=1 である。したがって、

y=C1e−2x+C2e−4x+e−x\boxed{y=C_1e^{-2x}+C_2e^{-4x}+e^{-x}}

を得る。

I.2​

x≠0x\ne0 の区間で Cauchy-Euler 型方程式として解く。同次解を y=xmy=x^m とおくと、

m(m−1)−2=(m−2)(m+1)=0,m(m-1)-2=(m-2)(m+1)=0,

したがって m=2,−1m=2,-1 である。また特解を yp=Ax3y_p=Ax^3 とおけば、

x2(6Ax)−2Ax3=4Ax3=4x3x^2(6Ax)-2Ax^3=4Ax^3=4x^3

より A=1A=1 である。ゆえに、

y=C1x2+C2x+x3(x≠0)\boxed{y=C_1x^2+\frac{C_2}{x}+x^3}\qquad(x\ne0)

となる。

II.1​

y=A0+A1x+A2x2y=A_0+A_1x+A_2x^2 を代入して xx の各次数の係数を比較すると、

0=(−2A22+4A2)x5+(−4A1A2+2A1)x4+(−4A0A2−2A12)x3+(−4A0A1−A2+2)x2−2A02x+A0.\begin{aligned} 0={}&(-2A_2^2+4A_2)x^5+(-4A_1A_2+2A_1)x^4\\ &+(-4A_0A_2-2A_1^2)x^3+(-4A_0A_1-A_2+2)x^2\\ &-2A_0^2x+A_0. \end{aligned}

これを満たす係数は

A0=0,A1=0,A2=2\boxed{A_0=0,\qquad A_1=0,\qquad A_2=2}

であり、特解は

yp=2x2y_p=2x^2

である。

II.2​

Riccati 方程式を

y′=22x3−1y2−1x(2x3−1)y−2x2x3−1y'=\frac{2}{2x^3-1}y^2-\frac{1}{x(2x^3-1)}y-\frac{2x}{2x^3-1}

と書き、既知の特解を用いて

y=2x2+1uy=2x^2+\frac1u

とおく。代入して整理すると、uu は線形方程式

u′+8x3−1x(2x3−1)u=−22x3−1u'+\frac{8x^3-1}{x(2x^3-1)}u=-\frac{2}{2x^3-1}

を満たす。ここで

8x3−1x(2x3−1)=1x+6x22x3−1\frac{8x^3-1}{x(2x^3-1)} =\frac1x+\frac{6x^2}{2x^3-1}

なので、積分因子は

μ(x)=x(2x3−1)\mu(x)=x(2x^3-1)

である。したがって、

ddx[x(2x3−1)u]=−2x\frac{d}{dx}\left[x(2x^3-1)u\right]=-2x

より

u=C−x2x(2x3−1).u=\frac{C-x^2}{x(2x^3-1)}.

よって一般の一パラメータ解は

y=2x2+x(2x3−1)C−x2=x(2Cx−1)C−x2\boxed{ y=2x^2+\frac{x(2x^3-1)}{C-x^2} =\frac{x(2Cx-1)}{C-x^2} }

である。また、変数変換で除外された特解 y=2x2y=2x^2 自身も元の方程式を満たす。以上の式は分母および微分項の係数が零にならない区間ごとに考える。

Reference​