跳到主要内容

九州大学 システム情報科学府 情報理工学専攻・電気電子工学専攻 2023年8月実施 解析学・微積分

Author

Casablanca, Miyake

Description

(1) 積分

I=0x5exp(x4)dxI = \int_{0}^{\infty}x^5\exp(-x^4)dx

を計算せよ。ただし, exp(x2)dx=π\int_{-\infty}^{\infty}\exp(-x^2)dx = \sqrt{\pi} を証明なしに用してよい。

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

dydx+y=xsinhx\frac{dy}{dx} + y = x\sinh x

(3) 複素関数 f(z)=1z4+1f(z) = \frac{1}{z^4 + 1} を考える。次の各問いに答えよ。

  • (a) f(z)f(z) の極をすべて求めよ。

  • (b) 下図に示す半円 CC に沿った複素積分 Cf(z)dz\oint_{C}f(z)dz を求めよ。ただし, R>1R > 1 とする。

题目描述

  1. 计算广义积分
    I=0x5ex4dx.I=\int_0^\infty x^5e^{-x^4}\,dx.
    可以不加证明地使用 ex2dx=π\int_{-\infty}^\infty e^{-x^2}\,dx=\sqrt\pi
  2. 求微分方程
    dydx+y=xsinhx\frac{dy}{dx}+y=x\sinh x
    的通解。
  3. 对复函数
    f(z)=1z4+1,f(z)=\frac1{z^4+1},
    求全部极点;再对上图所示、半径 R>1R>1 的半圆闭路 CC,计算 Cf(z)dz\oint_Cf(z)\,dz

考点

  • 广义积分:通过幂次代换把积分化为 Gamma/高斯型积分并计算。
  • 一阶线性微分方程:使用积分因子,并处理含双曲正弦的非齐次项。
  • 留数定理:求 z4+1z^4+1 的全部零点,判断半圆轮廓内的极点并汇总留数。

Kai - By Casablanca

(1)

(ex2dx)2=πx5ex4dx=x24dex4=14x2ex4+14ex4dx20x5ex4dx=140et2dt=π/8\begin{aligned} % (\int_{-\infty}^{\infty}e^{-x^2}dx)^2 &= \int e^{-x^2}dx \int e^{-y^2}dy = \int e^{-(x^2+y^2)dxdy} = \pi \\ \left(\int_{-\infty}^{\infty}e^{-x^2}dx \right)^2 &= \pi \\ \int x^5 e^{-x^4}dx &= \int-\frac{x^2}{4}de^{-x^4} = -\frac{1}{4}x^2e^{-x^4} + \frac{1}{4}\int e^{-x^4}dx^2 \\ \int_0^{\infty}x^5 e^{-x^4}dx &= \frac{1}{4}\int_0^{\infty}e^{-t^2}dt = \sqrt{\pi}/8 \end{aligned}

(2)

dydx+y=x(exex)/2ex(y+y)=12x(e2x1)(exy)=12xe2x12x(12xe2x12x)dx=14xe2x18e2xx24+C,\begin{aligned} \frac{dy}{dx} + y &= x(e^x - e^{-x})/2 \\ e^x(y' + y) &= \frac{1}{2}x(e^{2x} - 1) \\ (e^x y)' &= \frac{1}{2}xe^{2x} - \frac{1}{2}x \\ \int (\frac{1}{2}xe^{2x} - \frac{1}{2}x)dx &= \frac{1}{4}xe^{2x} - \frac{1}{8}e^{2x} - \frac{x^2}{4} + C, \end{aligned}
y=14xex18ex14x2ex+Cexy = \frac{1}{4}xe^x - \frac{1}{8}e^x - \frac{1}{4}x^2e^{-x} + Ce^{-x}

where CC is a constant.

(3)

(a)

Consider z4=1z^4 = -1,

we get

z1=22+22i,z2=2222i,z3=22+22i,z4=2222iz_1 = \frac{\sqrt{2}}{2} + \frac{\sqrt{2}}{2}i , z_2 = \frac{\sqrt{2}}{2} - \frac{\sqrt{2}}{2}i , z_3 = -\frac{\sqrt{2}}{2} + \frac{\sqrt{2}}{2}i , z_4 = -\frac{\sqrt{2}}{2} - \frac{\sqrt{2}}{2}i

and these are the poles

(b)

Cf(z)dz=2πiRes[f(z),22+22i]+2πiRes[f(z),22+22i]=22π\begin{aligned} \oint_Cf(z)dz &= 2\pi i \text{Res}[f(z),\frac{\sqrt{2}}{2} + \frac{\sqrt{2}}{2}i] + 2\pi i \text{Res}[f(z),\frac{\sqrt{2}}{2} + \frac{\sqrt{2}}{2}i] \\ &= \frac{\sqrt{2}}{2}\pi \end{aligned}

Kai - By Miyake

(1)

exp(x2)\exp(-x^2) は偶関数であるから、

exp(x2)dx=π \begin{align} \int_{- \infty}^\infty \exp \left( -x^2 \right) dx = \sqrt{\pi} \end{align}

より、

0exp(x2)dx=π2 \begin{align} \int_0^\infty \exp \left( -x^2 \right) dx = \frac{\sqrt{\pi}}{2} \tag{A} \end{align}

がわかる。

y=x2y=x^2 とおくと、 dy=2xdxdy=2xdx であり、次のように計算できる:

I=120y2exp(y2)dy=120y(12exp(y2))dy(  は yによる微分)=14[yexp(y2)]0+140exp(y2)dy=π8.( 式 (A) ) \begin{align} I &= \frac{1}{2} \int_0^\infty y^2 \exp \left( - y^2 \right) dy \\ &= \frac{1}{2} \int_0^\infty y \left( - \frac{1}{2} \exp \left( - y^2 \right) \right)' dy & \left( \ ' \text{ は } y { による微分 } \right) \\ &= - \frac{1}{4} \left[ y \exp \left( - y^2 \right) \right]_0^\infty + \frac{1}{4} \int_0^\infty \exp \left( - y^2 \right) dy \\ &= \frac{\pi}{8} . & ( \because \text{ 式 (A) } ) \end{align}

(2)

dydx+y=0 \begin{align} \frac{dy}{dx} + y = 0 \end{align}

の一般解は

y=Cex(C は積分定数 ) \begin{align} y &= Ce^{-x} & ( C \text{ は積分定数 } ) \end{align}

である。このことを考慮して、与えられた微分方程式

dydx+y=xsinhx \begin{align} \frac{dy}{dx} + y = x \sinh x \end{align}

y=zexy=ze^{-x} ( zzxx の関数 ) を代入すると、

dzdxexzex+zex=xsinhxdzdxex=xexex2dzdx=12(xe2xx)  z=12(xe2xx)dx=12(12xe2x14e2xx22)+C(C は積分定数 ) \begin{align} \frac{dz}{dx} e^{-x} - z e^{-x} + z e^{-x} &= x \sinh x \\ \frac{dz}{dx} e^{-x} &= x \frac{e^x - e^{-x}}{2} \\ \frac{dz}{dx} &= \frac{1}{2} \left( x e^{2x} - x \right) \\ \therefore \ \ z &= \frac{1}{2} \int \left( x e^{2x} - x \right) dx \\ &= \frac{1}{2} \left( \frac{1}{2} x e^{2x} - \frac{1}{4} e^{2x} - \frac{x^2}{2} \right) + C & ( C \text{ は積分定数 } ) \end{align}

を得るので、求める一般解は

y=14xex18ex14x2ex+Cex(C は積分定数 ) \begin{align} y &= \frac{1}{4} x e^x - \frac{1}{8} e^x - \frac{1}{4} x^2 e^{-x} + C e^{-x} & ( C \text{ は積分定数 } ) \end{align}

であることがわかる。

(3)

(a)

z4+1=0z^4+1=0z=reiθz= r e^{i \theta} ( r>0,θr \gt 0, \theta は実数) を代入すると

r4e4iθ=1 \begin{align} r^4 e^{4i \theta} = -1 \end{align}

であり、

r=1,  4θ=π+2π 整数  \begin{align} r = 1 , \ \ 4 \theta = \pi + 2 \pi \cdot \text{ 整数 } \end{align}

となる。よって、 f(z)f(z)

z=eπ/4,e3π/4,e5π/4,e7π/4=1+i2,1i2,1i2,1+i2 \begin{align} z &= e^{\pi/4}, e^{3\pi/4}, e^{5\pi/4}, e^{7\pi/4} \\ &= \frac{1+i}{\sqrt{2}}, \frac{1-i}{\sqrt{2}}, \frac{-1-i}{\sqrt{2}}, \frac{-1+i}{\sqrt{2}} \end{align}

に1位の極をもつ。

(b)

(a) で求めた4つの極は

z0=eπ/4=1+i2,z1=e3π/4=1+i2,z2=e5π/4=1i2,z3=e7π/4=1i2 \begin{align} z_0 &= e^{\pi/4} = \frac{1+i}{\sqrt{2}}, \\ z_1 &= e^{3\pi/4} = \frac{-1+i}{\sqrt{2}}, \\ z_2 &= e^{5\pi/4} = \frac{-1-i}{\sqrt{2}}, \\ z_3 &= e^{7\pi/4} = \frac{1-i}{\sqrt{2}} \end{align}

であるが、このうち CC の内側にあるのは z0,z1z_0, z_1 である。

z0z1=22=2,z0z2=2+2i2=2(1+i),z0z3=2i2=2i,z1z0=22=2,z1z2=2i2=2i,z1z3=2+2i2=2(1+i) \begin{align} z_0 - z_1 &= \frac{2}{\sqrt{2}} = \sqrt{2} ,\\ z_0 - z_2 &= \frac{2+2i}{\sqrt{2}} = \sqrt{2} (1+i) ,\\ z_0 - z_3 &= \frac{2i}{\sqrt{2}} = \sqrt{2} i ,\\ z_1 - z_0 &= \frac{-2}{\sqrt{2}} = - \sqrt{2} ,\\ z_1 - z_2 &= \frac{2i}{\sqrt{2}} = \sqrt{2} i ,\\ z_1 - z_3 &= \frac{-2+2i}{\sqrt{2}} = \sqrt{2} (-1+i) \end{align}

なので、 z0z_0 における留数は

limzz0(zz0)f(z)=1(z0z1)(z0z2)(z0z3)=122(1+i)=142(1+i) \begin{align} \lim_{z \to z_0} (z-z_0) f(z) &= \frac{1}{(z_0-z_1)(z_0-z_2)(z_0-z_3)} \\ &= \frac{1}{2 \sqrt{2} (-1+i)} \\ &= - \frac{1}{4 \sqrt{2}} (1+i) \end{align}

であり、 z1z_1 における留数は

limzz1(zz1)f(z)=1(z1z0)(z1z2)(z1z3)=122(1i)=142(1i) \begin{align} \lim_{z \to z_1} (z-z_1) f(z) &= \frac{1}{(z_1-z_0)(z_1-z_2)(z_1-z_3)} \\ &= \frac{1}{2 \sqrt{2} (1-i)} \\ &= \frac{1}{4 \sqrt{2}} (1-i) \end{align}

である。よって、留数定理より

Cf(z)dz=2πi(142(1+i)+142(1i))=π2 \begin{align} \oint_C f(z) dz &= 2 \pi i \left( - \frac{1}{4 \sqrt{2}} (1+i) + \frac{1}{4 \sqrt{2}} (1-i) \right) \\ &= \frac{\pi}{\sqrt{2}} \end{align}

がわかる。