8.3 有理函數之瑕積分 (Improper Integrals of Rational Functions)

本節目的在於利用留數定理計算微積分中的有理函數的瑕積分之值。一般而言瑕積分分成以下三類:

0f(x)dx=limb0bf(x)dx0f(x)dx=limaa0f(x)dxf(x)dx=limaacf(x)dx+limbcbf(x)dx, cR.\begin{align} &\int_0^\infty f(x) dx = \lim_{b\to\infty} \int_0^b f(x) dx \tag{第一型} \\ &\int_{-\infty}^0 f(x) dx = \lim_{a\to-\infty} \int_a^0 f(x) dx \tag{第二型} \\ &\int_{-\infty}^\infty f(x) dx = \lim_{a\to-\infty} \int_a^c f(x)dx+\lim_{b\to\infty} \int_c^b f(x) dx, ~c\in\mathbb{R}. \tag{第三型} \\ \end{align}

因此只要上述極限值存在,則稱對應的瑕積分存在。若(第三型)的瑕積分存在,則其值亦可計算如下:

f(x)dx=limRRRf(x)dx.\int_{-\infty}^\infty f(x)dx=\lim_{R\to\infty}\int_{-R}^R f(x)dx.

範例8.3-1. 考慮瑕積分 2xx2+1dx\displaystyle \int_{-\infty}^\infty \frac{2x}{x^2+1} dx,由於 2xx2+1dx=ln(x2+1)+C, CR\displaystyle \int \frac{2x}{x^2+1}dx = \ln(x^2+1)+C,~C\in\mathbb{R},得知

02xx2+1dx=limb0b2xx2+1dx=limbln(x2+1)0b=limbln(b2+1)=,\int_0^\infty \frac{2x}{x^2+1}dx = \lim_{b\to\infty} \int_0^b \frac{2x}{x^2+1}dx=\lim_{b\to\infty} \ln(x^2+1)\bigg|_{0}^{b}=\lim_{b\to\infty} \ln(b^2+1)=\infty,

亦即 2xx2+1dx\displaystyle \int_{-\infty}^\infty \frac{2x}{x^2+1} dx 不存在,然而

limRRR2xx2+1dx=limRln(x2+1)RR=0.\lim_{R\to\infty} \int_{-R}^R \frac{2x}{x^2+1}dx=\lim_{R\to\infty} \ln(x^2+1)\bigg|_{-R}^{R}=0.

亦即 limRRRf(x)dx\displaystyle \lim\limits_{R\to\infty}\int_{-R}^R f(x)dx 之值存在,不代表瑕積分 f(x)dx\displaystyle \int_{-\infty}^\infty f(x) dx 存在。

有理函數積分

定義8.3-1.fC(,)f\in C(-\infty,\infty),則定義 Cauchy 主值 (principal value, P.V.)為

P.V.f(x)dx=limRRRf(x)dx.P.V.\int_{-\infty}^\infty f(x)dx=\lim_{R\to\infty} \int_{-R}^R f(x)dx.

範例8.3-2. P.V.xdx=limRRRxdx=limRx22RR=0\displaystyle P.V. \int_{-\infty}^\infty x dx = \lim\limits_{R\to\infty} \int_{-R}^R x dx = \lim\limits_{R\to\infty} \frac{x^2}{2}\bigg|_{-R}^{R} = 0 ,但 xdx\displaystyle \int_{-\infty}^\infty xdx 不存在,此因 limb0bxdx=\displaystyle\lim\limits_{b\to\infty} \int_0^b x dx = \infty

範例8.3-3. 計算 P.V.1x2+1dx\displaystyle P.V. \int_{-\infty}^{\infty}\frac{1}{x^2+1}dx 之值。

定理8.3-1.f(z)=P(z)Q(z)\displaystyle f(z) = \frac{P(z)}{Q(z)},其中 deg(Q)deg(P)+2\deg(Q)\ge \deg(P)+2 。若 Q(x)0, xRQ(x)\neq 0,~\forall x\in\mathbb{R} ,則

P.V.f(x)dx=2πij=1kRes[PQ,zj] \displaystyle P.V. \int_{-\infty}^{\infty} f(x)dx = 2\pi i \sum\limits_{j=1}^k \text{Res}\left[\frac{P}{Q},z_j\right],

其中 zj, j=1,2,,kz_j,~j=1,2,\ldots,kff 落在上半平面之極點,即 Im(zj)0\text{Im}(z_j)\ge 0

範例8.3-4. 計算  dx(1+x2)(4+x2)\displaystyle \int_{-\infty}^{\infty} ~\frac{dx}{(1+x^2)(4+x^2)} 之值。

範例8.3-5. 計算 dx(4+x2)4\displaystyle \int_{-\infty}^{\infty} \frac{dx}{(4+x^2)^4} 之值。

定理8.3-1的作法,可以應用到更廣義的區域形式,如下例所示:

範例8.3-6. 計算 I=0dx1+x3I=\displaystyle \int_{0}^{\infty}\frac{dx}{1+x^3} 之值。

範例8.3-7.I=P.V.eax1+ebxdx, 0<a<bI=\displaystyle P.V.\int_{-\infty}^{\infty}\frac{e^{a x}}{1+e^{b x}}dx,~0<a<b ,證明

I=11e2πiablimRΓ1+Γ2+Γ3+Γ4eaz1+ebzdzI=\frac{1}{1-e^{2\pi i\frac{a}{b}}} \lim\limits_{R\to\infty}\int_{\Gamma_1+\Gamma_2+\Gamma_3+\Gamma_4} \frac{e^{a z}}{1+e^{b z}}dz

以及 II 之值,其中 Γ1, Γ2, Γ3\Gamma_1,~\Gamma_2,~\Gamma_3Γ4\Gamma_4 如下圖所示。

如若 0<a<10<a<1,則由範例8.3-7P.V.eax1+exdx=πsin(aπ)\displaystyle P.V. \int_{-\infty}^{\infty}\frac{e^{a x}}{1+e^x}dx=\frac{\pi}{\sin(a\pi)}。試問若欲計算 P.V.eaxex1dx\displaystyle P.V. \int_{-\infty}^{\infty}\frac{e^{a x}}{e^x-1}dx ,則圍線積分要如何取才行?

應用在級數求和

p(x)p(x) 為二階(含)以上實係數多項式且無整數零點,即 p(k)0, kZp(k)\neq 0, ~\forall k\in\mathbb{Z}。欲計算 k=01p(k)\sum\limits_{k=0}^\infty \frac{1}{p(k)}k=0(1)kp(k)\sum\limits_{k=0}^\infty \frac{(-1)^k}{p(k)} 之值,則可以從 cotπz\cot \pi zcscπz\csc\pi z 之極點來著手,亦即可考慮函數 F(z)=πcotπzp(z)F(z)=\frac{\pi \cot \pi z}{p(z)}的函數,則此函數之極點除了 p(z)p(z) 的零點外,尚有 cotπz\cot\pi z 之極點,即 sinπz=0\sin \pi z=0,可得 πz=kπ\pi z=k \pi,即 z=k=0,±1,±2,z=k=0, \pm 1, \pm 2,\ldots

以下先考慮一般形的級數和 n=0f(n)\sum\limits_{n=0}^\infty f(n),取以下的圍線 CNC_N

圖8.3.1 圍線 CNC_N 圖形,其四個邊從底邊順時間數,分別為路徑 Γ1, Γ2, Γ3, Γ4\Gamma_1,~\Gamma_2,~\Gamma_3,~\Gamma_4

CN=Γ1+Γ2+Γ3+Γ4,C_N=\Gamma_1+\Gamma_2+\Gamma_3+\Gamma_4,

其中

Γ1:z=(N+12)(ti), t[1,1],Γ2:z=(N+12)(1+ti), t[1,1],Γ3:z=(N+12)(t+i), t[1,1],Γ4:z=(N+12)(1ti), t[1,1].\begin{align*} \Gamma_1:z=(N+\frac12)(t-i),~t\in[-1,1],\quad & \Gamma_2:z=(N+\frac12)(1+ti),~t\in[-1,1],\\ \Gamma_3:z=(N+\frac12)(-t+i),~t\in[-1,1],\quad & \Gamma_4:z=(N+\frac12)(-1-ti),~t\in[-1,1]. \end{align*}

針對此圍線,函數 FF 之積分為

CNF(z)dz=2πik=NNRes[F,k]+2πij=1rNRes[F,zpj]\oint_{C_N} F(z) dz = 2\pi i \sum_{k=-N}^N \text{Res}[F,k]+2\pi i \sum_{j=1}^{r_N} \text{Res}[F,z_{p_j}]

其中 zpjz_{p_j}F(z)F(z) 在圍線 CN\mathcal{C}_N 內之第 jj 個零點,總數為 rNr_N

引理8.3-1. CNC_N 如圖8.3.1所示,則於此圍線上 cot(πz)<A|\cot(\pi z)|<A 其中 AA 為常數。

定理8.3-2.f(z)f(z) 延如圖8.3.1的 CNC_N 路徑滿足 f(z)M/zm, m>1|f(z)|\le M/|z|^m, ~m>1MM 為與 NN 無關之常數,則

n=f(n)=j=1kRes[f(z)πcot(πz),zpj],\sum_{n=-\infty}^\infty f(n)=-\sum_{j=1}^{k} \text{Res}[f(z) \pi \cot(\pi z),z_{p_j}],

其中 zpjz_{p_j}f(z)f(z) 在複數平面的第 jj 個極點,j=1,2,,kj=1,2,\ldots,k

範例8.3-8. 證明 n=1n2+a2=πacothπa\displaystyle \sum\limits_{n=-\infty}^\infty \frac{1}{n^2+a^2}=\frac{\pi}{a}\coth \pi a 其中 a>0a>0

推論8.3-1.f(z)=P(z)Q(z)f(z)=\frac{P(z)}{Q(z)},其中 P(z)P(z)Q(z)Q(z) 為多項式滿足 deg(Q)deg(P)+2\deg(Q)\ge \deg(P)+2 的條件,且 f(z)f(z) 之極點為 zp1, zp2, ,zpkz_{p_1},~z_{p_2},~\ldots,z_{p_k},並沒有落在整數上的極點,則

n=P(n)Q(n)=j=1kRes[πcot(πz)P(z)Q(z),zpj].\sum_{n=-\infty}^\infty\frac{P(n)}{Q(n)}=-\sum_{j=1}^k \text{Res}[\frac{\pi \cot(\pi z)P(z)}{Q(z)},z_{p_j}].

推論8.3-1 之特例如下:

推論8.3-2.f(z)=1p(z)f(z)=\frac{1}{p(z)},其中 p(z)p(z) 大於或等於二階多項式條件,且 p(z)p(z) 之零點為 zp1, zp2, ,zpkz_{p_1},~z_{p_2},~\ldots,z_{p_k},並且沒有落在整數點上,則

n=1p(n)=j=1kRes[πcot(πz)p(z),zpj].\sum_{n=-\infty}^\infty\frac{1}{p(n)}=-\sum_{j=1}^k \text{Res}[\frac{\pi \cot(\pi z)}{p(z)},z_{p_j}].

範例8.3.3. 證明 n=11n2=π6\displaystyle\sum\limits_{n=1}^\infty \frac{1}{n^2}=\frac{\pi}{6}

沿用上面的討論,留數定理可以應用在下面類型的無窮級數計算,歸納如下:

級數方法
n=f(n)\sum\limits_{n=-\infty}^\infty f(n)使用 CNf(z)πcot(πz)dz\displaystyle\oint_{C_N} f(z)\pi \cot(\pi z) dz ,其中 CNC_N 是以 (n+12)(±1±i)(n+\frac12)(\pm 1 \pm i) 為頂點正方形
n=(1)nf(n)\sum\limits_{n=-\infty}^\infty (-1)^nf(n)使用 CNf(z)πcsc(πz)dz\displaystyle\oint_{C_N} f(z)\pi \csc(\pi z) dz ,其中 CNC_N 是以 (n+12)(±1±i)(n+\frac12)(\pm 1 \pm i) 為頂點正方形
n=f(n+12)\sum\limits_{n=-\infty}^\infty f(n+\frac12)使用 CNf(z)πtan(πz)dz\displaystyle\oint_{C_N} f(z)\pi \tan(\pi z) dz ,其中 CNC_N 是以 n(±1±i)n(\pm 1 \pm i) 為頂點正方形
n=(1)nf(n+12)\sum\limits_{n=-\infty}^\infty (-1)^n f(n+\frac12)使用 CNf(z)πsec(πz)dz\displaystyle\oint_{C_N} f(z)\pi \sec(\pi z) dz ,其中 CNC_N 是以 n(±1±i)n(\pm 1 \pm i) 為頂點正方形

習題

  1. 證明 n=11n2+a2=π2acothaπ12a2\displaystyle \sum\limits_{n=1}^\infty \frac{1}{n^2+a^2}=\frac{\pi}{2a}\coth a\pi-\frac{1}{2a^2} 其中 a>0a>0
  1. 證明 n=11n4=π490\displaystyle\sum\limits_{n=1}^\infty \frac{1}{n^4} = \frac{\pi^4}{90}
  1. 證明 n=(1)nn2=π212\displaystyle\sum\limits_{n=-\infty}^\infty \frac{(-1)^n}{n^2} = -\frac{\pi^2} {12}
  1. 證明 n=1(1)nn4=7π4720\displaystyle\sum\limits_{n=1}^\infty \frac{(-1)^n}{n^4} = \frac{7\pi^4}{720}
  1. 證明 n=1(1)n(a+n)2=π2cot(πa)csc(πa)\displaystyle\sum\limits_{n=1}^\infty \frac{(-1)^n}{(a+n)^2} = \pi^2 \cot(\pi a)\csc(\pi a)

參考文獻

[SOCV] Spiegel, M. R., Lipschutz, S., Schiller, J. .J., Spellman, D., Schaum’s Outline of Complex Variables, 2nd ed., McGraw-Hill, Inc, USA, 2009.