跳到主要内容

九州大学 システム情報科学府 共通 2012年8月実施 微分方程

Author

思齐塾, 祭音Myyura

Description

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

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

(2) y2y3y=ex+e3x+cosxy'' - 2y' - 3y = e^x + e^{3x} + \cos x

(3) y4y+7y6y+2y=0y'''' - 4y''' + 7y'' - 6y' + 2y = 0

题目描述

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

y=3x+y5x3y5.y'=\frac{3x+y-5}{x-3y-5}.
y2y3y=ex+e3x+cosx.y''-2y'-3y=e^x+e^{3x}+\cos x.
y4y+7y6y+2y=0.y''''-4y'''+7y''-6y'+2y=0.

Kai

(1) y=3x+y5x3y5y' = \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+k5X3Y+h3k5\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+k5=03h + k - 5 = 0 and h3k5=0h - 3k - 5 = 0 .

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

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

Then dYdX=3X+YX3Y\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+vXX3vX=3+v13vv + X\frac{dv}{dX} = \frac{3X + vX}{X - 3vX} = \frac{3 + v}{1 - 3v}

XdvdX=3+v13vv=3+vv+3v213v=3+3v213vX\frac{dv}{dX} = \frac{3 + v}{1 - 3v} - v = \frac{3 + v - v + 3v^2}{1 - 3v} = \frac{3 + 3v^2}{1 - 3v}

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

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

1311+v2dvv1+v2dv=lnX+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)=lnX+C1\frac{1}{3} \arctan(v) - \frac{1}{2} \ln(1 + v^2) = \ln|X| + C_1

13arctan(YX)12ln(1+(YX)2)=lnX+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)=lnX+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)=lnX+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)+lnX=lnX+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+1x2)12ln((x2)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

(2) y2y3y=ex+e3x+cosxy'' - 2y' - 3y = e^x + e^{3x} + \cos x

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

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

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

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

For exe^x , try AexAe^x . Then Aex2Aex3Aex=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 cosx\cos x , try Ccosx+DsinxC\cos x + D\sin x . Then (CcosxDsinx)2(Csinx+Dcosx)3(Ccosx+Dsinx)=cosx(-C\cos x - D\sin x) - 2(-C\sin x + D\cos x) - 3(C\cos x + D\sin x) = \cos x .

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

The general solution is y=c1e3x+c2ex14ex+14xe3x15cosx110sinxy = 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) y4y+7y6y+2y=0y'''' - 4y''' + 7y'' - 6y' + 2y = 0

The characteristic equation is r44r3+7r26r+2=0r^4 - 4r^3 + 7r^2 - 6r + 2 = 0 .

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

(r44r3+7r26r+2)÷(r22r+1)=r22r+2(r^4 - 4r^3 + 7r^2 - 6r + 2) \div (r^2 - 2r + 1) = r^2 - 2r + 2 .

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

Therefore the general solution is y=c1ex+c2xex+c3excosx+c4exsinxy = c_1 e^x + c_2 xe^x + c_3 e^x \cos x + c_4 e^x \sin x .