跳到主要内容

九州大学 システム情報科学府 情報学専攻・情報知能工学専攻・電気電子工学専攻 共通 2013年度入学 数学 問2(微分方程式)

Author​

思齐塾, 祭音Myyura

Description​

出典:九州大学公式問題(保存版)。

次の微分方程式の一般解を求めよ.なお、 y′y' は関数 y(x)y(x) の xx に関する1階導関数を表している。

(1) y′=3x+y−5x−3y−5y' = \frac{3x + y - 5}{x - 3y - 5}

(2) y′′−2y′−3y=ex+e3x+cos⁡xy'' - 2y' - 3y = e^x + e^{3x} + \cos x

(3) y′′′′−4y′′′+7y′′−6y′+2y=0y'''' - 4y''' + 7y'' - 6y' + 2y = 0

题目描述​

求下列三个微分方程的通解,其中 y′y' 表示函数 y(x)y(x) 关于 xx 的一阶导数,y′′,y′′′,y′′′′y'',y''',y'''' 依次表示更高阶导数:

y′=3x+y−5x−3y−5.y'=\frac{3x+y-5}{x-3y-5}.
y′′−2y′−3y=ex+e3x+cos⁡x.y''-2y'-3y=e^x+e^{3x}+\cos x.
y′′′′−4y′′′+7y′′−6y′+2y=0.y''''-4y'''+7y''-6y'+2y=0.

Kai​

(1) y′=3x+y−5x−3y−5y' = \frac{3x + y - 5}{x - 3y - 5}​

Let x=X+h,y=Y+kx = X + h, y = Y + k . Then y′=dYdXy' = \frac{dY}{dX} .

dYdX=3(X+h)+(Y+k)−5(X+h)−3(Y+k)−5=3X+Y+3h+k−5X−3Y+h−3k−5\frac{dY}{dX} = \frac{3(X+h) + (Y+k) - 5}{(X+h) - 3(Y+k) - 5} = \frac{3X + Y + 3h + k - 5}{X - 3Y + h - 3k - 5}

Choose hh and kk such that 3h+k−5=03h + k - 5 = 0 and h−3k−5=0h - 3k - 5 = 0 .

Solving these equations gives h=2,k=−1h = 2, k = -1 .

So x=X+2,y=Y−1x = X + 2, y = Y - 1 .

Then dYdX=3X+YX−3Y\frac{dY}{dX} = \frac{3X + Y}{X - 3Y} .

Let Y=vXY = vX . Then dYdX=v+XdvdX\frac{dY}{dX} = v + X\frac{dv}{dX} .

v+XdvdX=3X+vXX−3vX=3+v1−3vv + X\frac{dv}{dX} = \frac{3X + vX}{X - 3vX} = \frac{3 + v}{1 - 3v}

XdvdX=3+v1−3v−v=3+v−v+3v21−3v=3+3v21−3vX\frac{dv}{dX} = \frac{3 + v}{1 - 3v} - v = \frac{3 + v - v + 3v^2}{1 - 3v} = \frac{3 + 3v^2}{1 - 3v}

1−3v3+3v2dv=dXX\frac{1 - 3v}{3 + 3v^2}dv = \frac{dX}{X}

∫1−3v3(1+v2)dv=∫dXX\int \frac{1 - 3v}{3(1 + v^2)}dv = \int \frac{dX}{X}

13∫11+v2dv−∫v1+v2dv=ln⁡∣X∣+C1\frac{1}{3} \int \frac{1}{1 + v^2}dv - \int \frac{v}{1 + v^2}dv = \ln|X| + C_1

13arctan⁡(v)−12ln⁡(1+v2)=ln⁡∣X∣+C1\frac{1}{3} \arctan(v) - \frac{1}{2} \ln(1 + v^2) = \ln|X| + C_1

13arctan⁡(YX)−12ln⁡(1+(YX)2)=ln⁡∣X∣+C1\frac{1}{3} \arctan(\frac{Y}{X}) - \frac{1}{2} \ln(1 + (\frac{Y}{X})^2) = \ln|X| + C_1

13arctan⁡(YX)−12ln⁡(X2+Y2X2)=ln⁡∣X∣+C1\frac{1}{3} \arctan(\frac{Y}{X}) - \frac{1}{2} \ln(\frac{X^2 + Y^2}{X^2}) = \ln|X| + C_1

13arctan⁡(YX)−12ln⁡(X2+Y2)+12ln⁡(X2)=ln⁡∣X∣+C1\frac{1}{3} \arctan(\frac{Y}{X}) - \frac{1}{2} \ln(X^2 + Y^2) + \frac{1}{2} \ln(X^2) = \ln|X| + C_1

13arctan⁡(YX)−12ln⁡(X2+Y2)+ln⁡∣X∣=ln⁡∣X∣+C1\frac{1}{3} \arctan(\frac{Y}{X}) - \frac{1}{2} \ln(X^2 + Y^2) + \ln|X| = \ln|X| + C_1

13arctan⁡(YX)−12ln⁡(X2+Y2)=C1\frac{1}{3} \arctan(\frac{Y}{X}) - \frac{1}{2} \ln(X^2 + Y^2) = C_1

13arctan⁡(y+1x−2)−12ln⁡((x−2)2+(y+1)2)=C1\frac{1}{3} \arctan(\frac{y+1}{x-2}) - \frac{1}{2} \ln((x-2)^2 + (y+1)^2) = C_1

This arctangent expression is local to intervals with x≠2x\ne2. More generally, choose a continuous argument θ=Arg⁡((x−2)+i(y+1))\theta=\operatorname{Arg}((x-2)+i(y+1)) along the solution and write

θ3−12log⁡((x−2)2+(y+1)2)=C1.\frac{\theta}{3}-\frac12\log\bigl((x-2)^2+(y+1)^2\bigr)=C_1.

This also describes regular solutions crossing x=2x=2 with y≠−1y\ne-1. The original differential equation requires x−3y−5≠0x-3y-5\ne0; the point (2,−1)(2,-1) is excluded.

(2) y′′−2y′−3y=ex+e3x+cos⁡xy'' - 2y' - 3y = e^x + e^{3x} + \cos x​

The homogeneous equation is y′′−2y′−3y=0y'' - 2y' - 3y = 0 .

The characteristic equation is r2−2r−3=0r^2 - 2r - 3 = 0 .

(r−3)(r+1)=0(r - 3)(r + 1) = 0 . So r=3,−1r = 3, -1 .

The homogeneous solution is yh=c1e3x+c2e−xy_h = c_1 e^{3x} + c_2 e^{-x} .

For exe^x , try AexAe^x . Then Aex−2Aex−3Aex=exAe^x - 2Ae^x - 3Ae^x = e^x . So −4A=1-4A = 1 , and A=−14A = -\frac{1}{4} .

For e3xe^{3x} , try Bxe3xBxe^{3x} . Then B(6e3x+9xe3x)−2B(e3x+3xe3x)−3Bxe3x=e3xB(6e^{3x} + 9xe^{3x}) - 2B(e^{3x} + 3xe^{3x}) - 3Bxe^{3x} = e^{3x} . So 4B=14B = 1 , and B=14B = \frac{1}{4} .

For cos⁡x\cos x , try Ccos⁡x+Dsin⁡xC\cos x + D\sin x . Then (−Ccos⁡x−Dsin⁡x)−2(−Csin⁡x+Dcos⁡x)−3(Ccos⁡x+Dsin⁡x)=cos⁡x(-C\cos x - D\sin x) - 2(-C\sin x + D\cos x) - 3(C\cos x + D\sin x) = \cos x .

−4C−2D=1,2C−4D=0-4C - 2D = 1, 2C - 4D = 0 . So C=2DC = 2D . −8D−2D=1-8D - 2D = 1 . So D=−110D = -\frac{1}{10} and C=−15C = -\frac{1}{5} .

The general solution is y=c1e3x+c2e−x−14ex+14xe3x−15cos⁡x−110sin⁡xy = c_1 e^{3x} + c_2 e^{-x} - \frac{1}{4} e^x + \frac{1}{4}xe^{3x} - \frac{1}{5}\cos x - \frac{1}{10} \sin x .

(3) y′′′′−4y′′′+7y′′−6y′+2y=0y'''' - 4y''' + 7y'' - 6y' + 2y = 0​

The characteristic equation is r4−4r3+7r2−6r+2=0r^4 - 4r^3 + 7r^2 - 6r + 2 = 0 .

By inspection, r=1r = 1 is a root twice. So (r−1)2=r2−2r+1(r-1)^2 = r^2 - 2r + 1 is a factor.

(r4−4r3+7r2−6r+2)÷(r2−2r+1)=r2−2r+2(r^4 - 4r^3 + 7r^2 - 6r + 2) \div (r^2 - 2r + 1) = r^2 - 2r + 2 .

So r2−2r+2=0r^2 - 2r + 2 = 0 . r=2±4−82=1±ir = \frac{2 \pm \sqrt{4 - 8}}{2} = 1 \pm i .

Therefore the general solution is y=c1ex+c2xex+c3excos⁡x+c4exsin⁡xy = c_1 e^x + c_2 xe^x + c_3 e^x \cos x + c_4 e^x \sin x .