部分積分

積分公式を一番よく覚えているのは大学入試直前ではないでしょうか. 大学生以上になると授業での演習量が減るのでどんどん忘れて行きます. 授業の最後にたまに演習問題があることがあるんですが, 部分積分がさっぱり分からなくなっていて問題が解けなかったので復習しました.

部分積分の公式

つぎの積分公式

\int f'(x)g(x)\,dx=f(x)g(x)-\int f(x)g'(x)\,dx

を部分積分といいます.2つの関数のかけ算になっている被積分関数の1つ(この場合は f(x) ) だけを部分的に積分するので部分積分といいます.定積分の場合は

\int_a^b f'(x)g(x)\,dx=\left[f(x)g(x)\right]_a^b-\int_a^b f(x)g'(x)\,dx

となります.それぞれの公式の証明は教科書に書いてあると思いますが,部分積分の公式は覚えておくべきです. これだけではよく分からないので例題をみて行きましょう.

例題1. (e^x)x の積分

つぎの積分を実行せよ.

\int e^xx\,dx

被積分関数がかけ算になっています.

\int e^x\,dx=e^x ,\quad \int x\,dx=\frac{1}{2}x^2

ということは簡単に分かりますから, e^xx が別々なら積分は簡単です. しかしそれらの掛け算だとどうしていいか分かりません. こういうときに部分積分を思い出します.

\int e^xx\,dx=\int (e^x)'x\,dx

ですから,部分積分の公式を使うと

\int e^xx\,dx &= \int (e^x)'x\,dx\\ &= e^xx-\int e^x(x)'\,dx\\ &= e^xx-\int e^x\,dx\\ &= e^xx-e^x+C

という具合に積分が実行できます.補足ですが,逆に

\int e^xx\,dx &= \int e^x\left(\frac{1}{2}x^2\right)'\,dx\\ &= e^x\frac{1}{2}x^2 - \int (e^x)'\frac{1}{2}x^2\,dx\\ &= e^x\frac{1}{2}x^2 - \int e^x\frac{1}{2}x^2\,dx\\ &= \cdots

としてしまうといつまでたっても計算が終わりません. 2つの関数のうちどちらを (f(x))' の形にすればいいかは, 何度か部分積分の問題を解いていったら自然に分かると思います [*]

[*]もちろん,掛け算の形になっているからといって, いつでも部分積分が通用するわけではありません. だから積分って難しいんですよね….

例題2. x sin x の積分

つぎの積分を実行せよ.

\int x\sin x\,dx

微分して \sin x になるのは -\cos x ですから,

\int x\sin x\,dx=\int x(-\cos x)'\,dx

として部分積分すればいいです.

\int x\sin x\,dx &= \int x(-\cos x)'\,dx\\ &= x(-\cos x)-\int (x)'(-\cos x)\,dx\\ &= -x\cos x+\int \cos x\,dx\\ &= -x\cos x+\sin x+C

例題3. log x の積分

つぎの積分を実行せよ.

\int \log x\,dx

これは被積分関数が掛け算になっていません.一見簡単そうにみえますが, \log x の積分公式はないのですぐには解けません.そこで部分積分の公式をつかうために,

\int \log x\,dx = \int 1\cdot\log x\,dx

というふうに無理矢理掛け算にしておきます. x で微分して1になるのは x ですから, 部分積分の公式からつぎのように積分計算ができます.

\int \log x\,dx &= \int (x)'\cdot\log x\,dx\\ &= x\log x-\int x(\log x)'\,dx\\ &= x\log x-\int x\frac{1}{x}\,dx\\ &= x\log x-\int 1\,dx\\ &= x\log x-x+C

例題4. 少し複雑な積分

1種類の単原子分子からなる気体のマクスウェルの速度分布関数は次式で表される.

f(v)=4\pi\left(\frac{m}{2\pi k_BT}\right)^{3/2}v^2\exp\left(-\frac{mv^2}{2k_BT}\right)

この気体分子の速度の平均値 <x> は,速度に上式を掛け, v 空間で積分することで求められる. <x> を計算せよ.

この問題が例の,崎間が授業中解けなかったってやつです(他の人たちは解けてた).まず,問題文の通りにやると

<x> &= \int_0^{\infty}4\pi\left(\frac{m}{2\pi k_BT}\right)^{3/2}v^2 \exp\left(-\frac{mv^2}{2k_BT}\right)v\,dv\\    &= 4\pi\left(\frac{m}{2\pi k_BT}\right)^{3/2} \int_0^{\infty}v^3\exp\left(-\frac{mv^2}{2k_BT}\right)\,dv

となります.このままではやりにくいので v^2=x と置いて積分変数を置換します.

\frac{d}{dv}v^2=\frac{dx}{dv} \to dv=\frac{1}{2v}dx

より

\int_0^{\infty}v^3\exp\left(-\frac{mv^2}{2k_BT}\right)\,dv &= \int_0^{\infty}v^3\exp\left(-\frac{mv^2}{2k_BT}\right)\,\frac{1}{2v}dx\\ &= \frac{1}{2}\int_0^{\infty}v^2\exp\left(-\frac{mv^2}{2k_BT}\right)\,dx\\ &= \frac{1}{2}\int_0^{\infty}x\exp\left(-\frac{mx}{2k_BT}\right)\,dx

です.ここで部分積分を行います,ってうわーなんだこれは.メゲずに黙々とやりましょう.

\frac{1}{2}\int_0^{\infty}&x\exp\left(-\frac{mx}{2k_BT}\right)\,dx\\ &= \frac{1}{2}\int_0^{\infty}x\left(-\frac{2k_BT}{m} \exp\left(-\frac{mx}{2k_BT}\right)\right)'\,dx\\ &= \frac{1}{2}\left[x\left(-\frac{2k_BT}{m}\exp\left(-\frac{mx}{2k_BT}\right)\right)\right]_0^{\infty}   -\frac{1}{2}\int_0^{\infty}(x)'\left(-\frac{2k_BT}{m}\exp\left(-\frac{mx}{2k_BT}\right)\right)\,dx\\ &= -\frac{k_BT}{m}\left[x\exp\left(-\frac{mx}{2k_BT}\right)\right]_0^{\infty} +\frac{k_BT}{m}\int_0^{\infty}\exp\left(-\frac{mx}{2k_BT}\right)\,dx\\ &= 0+\frac{k_BT}{m}\left[-\frac{2k_BT}{m} \exp\left(-\frac{mx}{2k_BT}\right)\right]_0^{\infty}\\ &= -\frac{2{k_B}^2T^2}{m^2}(0-1)\\ &= \frac{2{k_B}^2T^2}{m^2}

なんとかスッキリしました. exp(-\infty)=0 になるのはいいですよね.

これでやっと,求めるべき平均速度 <x>

<x> &= 4\pi\left(\frac{m}{2\pi k_BT}\right)^{3/2}\frac{2{k_B}^2T^2}{m^2}\\    &= 4\pi\frac{m}{2\pi k_BT}\sqrt{\frac{m}{2\pi k_BT}}\frac{2{k_B}^2T^2}{m^2}\\    &= \frac{4k_BT}{m}\sqrt{\frac{m}{2\pi k_BT}}\\    &= \sqrt{\frac{4^2{k_B}^2T^2}{m^2}\frac{m}{2\pi k_BT}}\\    &= \sqrt{\frac{8k_BT}{\pi m}}

と求まります.やたらと繁雑でしたが,計算間違いをしないようにしましょう (この記事を書くとき僕は何回も計算を間違えてました).