部分積分
| class = plainlist | titlestyle = padding-bottom:0.25em; | pretitle = Part of a series of articles about | title = 解析学 | listtitlestyle = text-align:center; | liststyle = border-top:1px solid #aaa;padding-top:0.15em;border-bottom:1px solid #aaa; | expanded =
| abovestyle = padding:0.15em 0.25em 0.3em;font-weight:normal; | above =
テンプレート:Endflatlistテンプレート:Startflatlist
| list2name = differential | list2titlestyle = display:block;margin-top:0.65em; | list2title = テンプレート:Bigger | list2 =
定義 |
---|
テンプレート:Endflatlist |
概念 |
テンプレート:Endflatlist |
法則と恒等式 |
テンプレート:Endflatlist |
| list3name = integral | list3title = テンプレート:Bigger | list3 =
定義 |
---|
Integration by |
| list4name = series | list4title = テンプレート:Bigger | list4 =
収束判定法 |
---|
| list5name = vector | list5title = テンプレート:Bigger | list5 =
定理 |
---|
| list6name = multivariable | list6title = テンプレート:Bigger | list6 =
Formalisms |
---|
定義 |
| list7name = specialized | list7title = テンプレート:Bigger | list7 = テンプレート:Startflatlist
}}部分積分(ぶぶんせきぶん、英: Integration by parts)とは、微分積分学・解析学における関数の積の積分に関する定理であり、積の積分をより計算が容易な積分に変形するために頻繁に使われる手法である。
具体的には、2つの微分可能な関数 u(x)、v(x)、区間 a ≤ x ≤ b に対して成り立つ以下のような関係式を指す[1]。
- [math]\int_a^b u(x)v'(x)\,dx = \left[u(x)v(x)\right]_a^b - \int_a^b u'(x)v(x)\,dx [/math]
不定積分の場合であれば、同様に以下の関係式が成り立つ。
- [math]\int u(x) v'(x)\, dx = u(x) v(x) - \int u'(x) v(x)\, dx\![/math]
またはより簡潔に
- [math]\int u \, dv=uv-\int v \, du.\![/math]
と表記される。ここで du と dv は x の関数 u, v の微分、即ち
- [math]du=u'(x)dx, \quad dv=v'(x)dx [/math]
である。
Contents
導出
上記の定理は以下にように導出される。 u(x) と v(x) がともに微分可能な関数であるとき、積の微分法則(ライプニッツ則)より
- [math]\frac{d}{dx}\left(u(x)v(x)\right) = v(x) \frac{d}{dx}\left(u(x)\right) + u(x) \frac{d}{dx}\left(v(x)\right)\![/math]
両辺を区間 a ≤ x ≤ b で x に関して積分して
- [math]\int_a^b \frac{d}{dx}\left(u(x)v(x)\right)\,dx = \int_a^b u'(x)v(x)\,dx + \int_a^b u(x)v'(x)\,dx [/math]
ここで微分積分学の基本定理より、
- [math]\int_a^b \frac{d}{dx}\left(u(x)v(x)\right)\,dx = \left[u(x)v(x)\right]_a^b[/math]
であるから、
- [math]\left[u(x)v(x)\right]_a^b = \int_a^b u'(x)v(x)\,dx + \int_a^b u(x)v'(x)\,dx [/math]
即ち以下の部分積分の公式を得る。
- [math]\int_a^b u(x)v'(x)\,dx = \left[u(x)v(x)\right]_a^b - \int_a^b u'(x)v(x)\,dx [/math]
不定積分の場合も同様に導出出来る。
ここで左辺の ∫uv′dx は v′ (v の 導関数) を含んでいるから、まず v (v′ の 原始関数)を見つける必要があり、次いで部分積分の公式を適用し、積分 ∫vu′dx を計算する。
(具体的な計算例は後述)
視覚的な解釈
パラメーター t によって (x, y) = (f(t), g(t)) で表された曲線を定義する。この曲線が局所的に一対一対応であると仮定すると、
- [math]x(y) = f(g^{-1}(y))[/math]
- [math]y(x) = g(f^{-1}(x))[/math]
青色の領域の面積は、
- [math]A_1=\int_{y_1}^{y_2}x(y)dy[/math]
同様に赤色の領域の面積は、
- [math]A_2=\int_{x_1}^{x_2}y(x)dx[/math]
にそれぞれ対応する。
A1 と A2 を足し合わせた領域全体は、大きい方の長方形の面積 x2y2 から小さい方の長方形の面積 x1y1 を除いたものに等しい。
- [math]\overbrace{\int_{y_1}^{y_2}x(y)dy}^{A_1}+\overbrace{\int_{x_1}^{x_2}y(x)dx}^{A_2}=\biggl.x_iy_i\biggl|_{i=1}^{i=2}[/math]
近傍で曲線が滑らかであれば、これは不定積分に一般化できる。
- [math]\int xdy + \int y dx = xy[/math]
変形して、
- [math]\int xdy = xy - \int y dx[/math]
つまり部分積分は、青色の領域の面積が領域全体の面積と赤色の領域の面積から導かれることに相当すると考える事が出来る。
またこのように可視化することにより、関数 f(x) の積分が分かっている時に逆関数 f−1(x) の積分が部分積分で求められることが理解出来る。実際、関数 x(y) と y(x) は逆関数の関係にあり、積分 ∫x dy は ∫y dx が分かっていれば上記のようにして計算可能である。
部分積分を用いた積分計算
基本方針
部分積分は機械的に積分を求められる方法ではなく、むしろある程度の試行錯誤を要する場合がある。基本的な方針は、ある一つの関数が与えられた時に、それを部分積分公式に当てはめて変形した場合に出現する積分項がもとの積分よりも計算が容易になるように、その関数を2つの関数の積 u(x)v(x) に分割するというものである[2]。下記の式は良い分割の方法を探すのに役立つであろう。
- [math]\int uv \,dx = u \int v \,dx - \int \left ( u' \int v \,dx \right )\,dx\![/math]
右辺で u は微分されて、逆に v は積分されていることに注意。即ち、微分された時に単純な形になる関数を u に、また積分された時に単純な形になる関数を v に選択するのが良いことが分かる。簡単な例として以下の積分を考えてみると、
- [math]\int \frac{\ln x}{x^2}\,dx\![/math]
ln x を微分すると 1/x であることから、(ln x) を u の部分として選択し、また 1/x2 の不定積分が −1/x であることから、(1/x2) を v の部分として選択する。すると公式により、
- [math]\int \frac{\ln x}{x^2}\,dx = -\frac{\ln x}{x} - \int \biggl( \frac{1}{x} \biggr) \biggl( -\frac{1}{x} \biggr) \, dx\![/math]
−1/x2 の不定積分は 1/x となる。
別の例として、u' (∫v dx) の積が約分により簡単な形になるように u と v を選択することもある。例えば次の例では、
- [math]\int\sec^2x\ln|\sin x|dx[/math]
u(x) = ln |sin x| として、また v(x) = sec2x とすると、 u を微分すると 合成関数の微分法により 1/ tan x となり、 v を積分すると tan x となる。したがって公式により、
- [math]\int\sec^2x\ln|\sin x|dx=\tan x\ln|\sin x|-\int\tan x\frac{1}{\tan x}dx[/math]
被積分関数は 1 となり、積分して x となる。積が簡単な形になる組み合わせを探すにはある程度の試行錯誤が必要なことがある。
その他いくつかのテクニックを以下の例で示す。
例
- 多項式と三角関数
- [math]I=\int x\cos (x) \,dx\,[/math]
この積分を計算するには、
- [math]u = x \Rightarrow d u = dx[/math]
- [math]dv = \cos(x)\,dx \Rightarrow v = \int\cos(x)\,dx = \sin x[/math]
とすると、
- [math] \begin{align} \int x\cos (x) \,dx & = \int u \, dv \\ & = uv - \int v \, du \\ & = x\sin (x) - \int \sin (x) \,dx \\ & = x\sin (x) + \cos (x) + C, \end{align} \![/math]
ここで C は積分定数である。
下記の式における x のより高位の累乗では、
- [math]\int x^n e^x \,dx,\,\int x^n\sin (x) \,dx,\,\int x^n\cos (x) \,dx\,[/math]
部分積分を繰り返し使って同様に計算出来る。1回部分積分を適用する度に x の指数が1ずつ下がる。
- 指数関数と三角関数
部分積分の仕組みを考えるためによく使われる例として、
- [math]I=\int e^{x} \cos (x) \,dx [/math]
を計算する。ここでは、部分積分を2回行う。最初に
- [math]u = \cos(x) \Rightarrow du = -\sin(x)\,dx[/math]
- [math]dv = e^x \, dx \Rightarrow v = \int e^x \,dx = e^x[/math]
とすると、
- [math]\int e^{x} \cos (x) \,dx = e^{x} \cos (x) + \int e^{x} \sin (x) \,dx.\![/math]
となる。残った積分項に対して再度部分積分を行う。
- [math]u = \sin(x) \Rightarrow du = \cos(x)\, dx[/math]
- [math]dv = e^x \,dx \Rightarrow v = \int e^x \,dx = e^x[/math]
として、
- [math]\int e^x \sin (x) \,dx = e^x \sin (x) - \int e^x \cos (x) \,dx [/math]
これらを組み合わせて、
- [math]\int e^x \cos (x) \,dx = e^x \cos (x) + e^x \sin (x) - \int e^x \cos (x) \,dx [/math]
同じ積分項が等式の両辺に出現しているので
- [math]2 \int e^{x} \cos (x) \,dx = e^{x} ( \sin (x) + \cos (x) ) + C\![/math]
と変形出来て、
- [math]\int e^x \cos (x) \,dx = {e^x ( \sin (x) + \cos (x) ) \over 2} + C_2\![/math]
となる。C と C2 ( = C/2) は積分定数である。
[math]\int \sec^3 x \, dx[/math] のような積分も同様の方法を使って計算出来る。
- 関数に形式的に1を掛ける
更によく知られた例を挙げる。被積分関数を1とそれ自身の積と考えて部分積分を行う方法である。これは、被積分関数の導関数が分かっていて、更にその導関数に xを乗じた関数の積分が計算可能な場合に有効である。
最初の例として ∫ ln(x) dx を考える。これを以下のように1と自身の積として考えて、
- [math]I=\int \ln (x) \cdot 1 \,dx.\![/math]
次のようにおくと、
- [math]u = \ln(x) \Rightarrow du = \frac{dx}{x}[/math]
- [math]dv = dx \Rightarrow v = x\ [/math]
以下のように計算出来る[3]。
- [math] \begin{align} \int \ln (x) \,dx & = x \ln (x) - \int \frac{x}{x} \,dx \\ & = x \ln (x) - \int 1 \,dx \\ & = x \ln (x) - x + C \end{align} [/math]
次の例として arctan(x) の積分を考える。
- [math]I=\int \arctan (x) \,dx [/math]
これを以下のように書き換える。
- [math]\int \arctan (x) \cdot 1 \,dx [/math]
次のようにおくと、
- [math]u = \arctan(x) \Rightarrow du = \frac{dx}{1 + x^2}[/math]
- [math]dv =dx \Rightarrow v = x[/math]
以下のように計算出来る。
- [math] \begin{align} \int \arctan (x) \,dx & = x \arctan (x) - \int \frac{x}{1 + x^2} \,dx \\[8pt] & = x \arctan (x) - {1 \over 2} \ln \left( 1 + x^2 \right) + C \end{align} [/math]
ここでは逆関数の微分法を使用した。
部分積分の再帰的適用
部分積分を ∫ v du に対して再帰的に適用することにより、次の公式を得る。
- [math]\int uv = u v_1 - u' v_2 + u'' v_3 - \cdots + (-1)^{n-1}\ u^{(n-1)} \ v_{n} + (-1)^n \int{u^{(n)}v_{n}}.\![/math]
ここで、u′ は u の1次導関数、u′′ は2次導関数であり、u(n) は n 次導関数を表す。[math]v_{n}[/math] は以下のように定義される。
- [math]v_{n+1}(x)=\int\! \int\ \cdots \int v \ (dx)^{n+1}.\![/math]
上記の式は、uv1 から開始して1つ目の項は順に微分して行き、2つ目の項は積分して行けば計算出来る(同時に符号を反転しながらであるが)。特に、u(k + 1) があるk + 1 で 0 になる時には u(k) の項までで終了するため、便利な公式である。
拡張
多因子への拡張
(積の微分法則の一般化も参照のこと)
3つの関数 u(x), v(x), w(x) の積の微分法則に対して積分を行うと、同様に以下のような結果を得る。
- [math]\int_a^b u v \, dw = u v w - \int_a^b u w \, dv - \int_a^b v w \, du[/math]
一般的に n個の関数の積の場合は、
- [math]\frac{d}{dx} \left(\prod_{i=1}^n u_i(x) \right)= \sum_{j=1}^n \prod_{i\neq j}^n u_i(x) \frac{du_j(x)}{dx} [/math]
即ち、
- [math] \Bigl[ \prod_{i=1}^n u_i(x) \Bigr]_a^b = \sum_{j=1}^n \int_a^b \prod_{i\neq j}^n u_i(x) \, du_j(x) [/math]
ここで右辺の積は、同じ項で微分を取った関数を除く全ての関数の積を取るものとする。
スティルチェス積分
リーマン=スティルチェス積分(またはスティルチェス積分)とは、スティルチェスによるリーマン積分の拡張である。
リーマン=スティルチェス積分に関しても、被積分関数 f および積分関数 g に対して部分積分公式が
- [math]\int_a^b f\,dg=f(b)g(b)-f(a)g(a)-\int_a^b g\,df[/math]
なる形で成り立つ。
また、リーマン=スティルチェス積分および(狭義の)ルベーグ積分の一般化であるルベーグ=スティルチェス積分(またはルベーグ=ラドン積分)に対しても、以下の形で部分積分公式が定式化される。
2つの有界変動関数 U, V に対して U または V のいずれかが連続、若しくは U および V がともに正常("regular")となるような点では、
- [math]\int_a^b U\,dV+\int_a^b V\,dU=U(b+)V(b+)-U(a-)V(a-),\quad(a\lt b)[/math]
が成立する。
詳細はリーマン=スティルチェス積分およびルベーグ=スティルチェス積分を参照。
高次元への拡張
部分積分を高次元の場合に対して拡張することが出来る。
[math]\Omega \subset \mathbb{R}^n[/math] を区分的に滑らかな境界 [math]\Gamma[/math] を持つ有界な開集合とし、[math]\mathbf{n}[/math] を [math]\Gamma[/math] への外向き単位面法線ベクトル、[math]u[/math] と [math]\mathbf{v}[/math] をそれぞれ [math]\Omega[/math] の閉包において連続微分可能な関数およびベクトル値関数として定義する。
この時、 [math]\nabla \cdot (u \mathbf{v}) = \nabla u \cdot \mathbf{v} + u \nabla \cdot \mathbf{v}[/math] に対してガウスの発散定理を適用すると、
- [math] \int_\Omega \nabla \cdot (u \mathbf{v}) \, d\Omega = \int_\Omega ( \nabla u \cdot \mathbf{v} + u \nabla \cdot \mathbf{v}) \, d\Omega =\int_{\Gamma} (u \mathbf{v}) \cdot \mathbf{n}\, d\Gamma[/math]
であるから、以下の部分積分公式が得られる。
- [math] \int_{\Omega} \nabla u \cdot \mathbf{v}\, d\Omega = \int_{\Gamma} u (\mathbf{v}\cdot \mathbf{n})\, d\Gamma - \int_\Omega u\, \nabla\cdot\mathbf{v}\, d\Omega[/math]
また、[math]\mathbf{v}=\nabla v[/math],[math]v\in C^2(\bar{\Omega} )[/math] なる [math]v[/math] で表される時、
- [math] \int_{\Omega} \nabla u \cdot \nabla v\, d\Omega = \int_{\Gamma} u\, \nabla v\cdot \mathbf{n}\, d\Gamma - \int_\Omega u\, \nabla^2 v\, d\Omega[/math]
となり、グリーンの第一恒等式(Green's first identity)が得られる。
同様に、任意の階数の微分可能テンソル場 [math]\boldsymbol{F}[/math] と [math]\boldsymbol{G}[/math] に対して、発散定理より以下の部分積分公式が導かれる。
- [math] \int_{\Omega} \boldsymbol{F}\otimes\boldsymbol{\nabla}\boldsymbol{G}\, d\Omega = \int_{\Gamma} \mathbf{n}\otimes(\boldsymbol{F}\otimes\boldsymbol{G})\, d\Gamma - \int_{\Omega} \boldsymbol{G}\otimes\boldsymbol{\nabla}\boldsymbol{F}\, d\Omega [/math]
ここで [math]\otimes[/math] はテンソル積を表す。[math]\boldsymbol{F}[/math] が恒等テンソルに等しい時は、発散定理の式を得る。
- [math] \int_{\Omega}\boldsymbol{\nabla}\boldsymbol{G}\, d\Omega = \int_{\Gamma} \mathbf{n}\otimes\boldsymbol{G}\, d\Gamma \, [/math]
添字表記で表すと以下のようになる。
- [math] \int_{\Omega} F_{ijk....}\,G_{lmn...,p}\, d\Omega = \int_{\Gamma} n_p\,F_{ijk...}\,G_{lmn...}\, d\Gamma - \int_{\Omega} G_{lmn...}\,F_{ijk...,p}\, d\Omega \, [/math]
ここで[math]\boldsymbol{F}[/math] と [math]\boldsymbol{G}[/math] がともに2階のテンソルであるような特殊な場合を考え、1つの添字の縮約を取ると、
- [math] \int_{\Omega} F_{ij}\,G_{pj,p}\, d\Omega = \int_{\Gamma} n_p\,F_{ij}\,G_{pj}\, d\Gamma - \int_{\Omega} G_{pj}\,F_{ij,p}\, d\Omega \, [/math]
即ち
- [math] \int_{\Omega} \boldsymbol{F}\cdot(\boldsymbol{\nabla}\cdot\boldsymbol{G})\, d\Omega = \int_{\Gamma} \mathbf{n}\cdot(\boldsymbol{G}\cdot\boldsymbol{F}^T)\, d\Gamma - \int_{\Omega} (\boldsymbol{\nabla}\boldsymbol{F}):\boldsymbol{G}^T\, d\Omega \, [/math]
となる。
応用
部分積分の解析学におけるいくつかの応用例を挙げる。
特殊関数
ガンマ関数は広義積分を用いて定義される特殊関数である。部分積分を使うと、これが階乗の拡張になっていることが分かる。
- [math]\begin{align} \Gamma(z) & = \int_0^\infty d\lambda e^{-\lambda} \lambda^{z-1} \\ & = - \int_0^\infty d\left(e^{-\lambda}\right) \lambda^{z-1} \\ & = - \left[e^{-\lambda}\lambda^{z-1}\right]_0^\infty + \int_0^\infty d\left(\lambda^{z-1}\right) e^{-\lambda} \\ & = 0 + \int_0^\infty d\lambda\left(z-1\right) \lambda^{z-2} e^{-\lambda} \\ & = (z-1)\Gamma(z-1) \\ \end{align} [/math]
このようにして、以下のよく知られた等式が得られる。
- [math]\Gamma(z) = (z-1)\Gamma(z-1)\,[/math]
整数 z に対してこの公式を繰り返し適用することで階乗が得られる。
- [math]\Gamma(z+1) = z![/math]
調和解析
調和解析、特にフーリエ解析における部分積分の応用例を挙げる。よく知られた例として、関数のフーリエ変換の収束が、関数の滑らかさに依存していることを示すものである。
- 導関数のフーリエ変換
f が k 回連続微分可能であり、更に k 次までの導関数が無限大で 0 に収束する時、そのフーリエ変換は以下の関係式を満たす。
- [math](\mathcal{F}f^{(k)})(\xi) = (2\pi i\xi)^k \mathcal{F}f(\xi)[/math]
ここで f (k) は f の k 次導関数を表す。
導関数のフーリエ変換に対して部分積分を適用すると、以下の結果を得る。
- [math]\begin{align} (\mathcal{F}f')(\xi) &= \int_{-\infty}^\infty e^{-2\pi iy\xi} f'(y)\,dy \\ &=\left[e^{-2\pi iy\xi} f(y)\right]_{-\infty}^\infty - \int_{-\infty}^\infty (-2\pi i\xi e^{-2\pi iy\xi}) f(y)\,dy \\ &=2\pi i\xi \int_{-\infty}^\infty e^{-2\pi iy\xi} f(y)\,dy \\ &=2\pi i\xi \mathcal{F}f(\xi). \end{align}[/math]
この結果を繰り返し適用することによって、一般の k に対する結果が得られる。同様の手法は導関数のラプラス変換を求める際にも利用出来る。
- フーリエ変換の収束
上記の結果により、f と f (k) が積分可能ならば、
- [math]\vert\mathcal{F}f(\xi)\vert \leq \frac{I(f)}{1+\vert 2\pi\xi\vert^k}[/math], ただし [math]I(f):=\int_{-\infty}^\infty\Bigl(\vert f(y)\vert + \vert f^{(k)}(y)\vert\Bigr) dy[/math].
言い換えると、f がこれらの条件を満足するならば、そのフーリエ変換は無限大で高々 1/|ξ|k のオーダーで収束するということである。特に、k ≥ 2 ならばフーリエ変換は積分可能である。
証明にはフーリエ変換の定義から直ちに得られる次の関係を用いる。
- [math]\vert\mathcal{F}f(\xi)\vert \leq \int_{-\infty}^\infty \vert f(y) \vert \,dy[/math]
節の冒頭で述べたのと同様の考え方により、次の結果が得られる。
- [math]\vert(2\pi i\xi)^k \mathcal{F}f(\xi)\vert \leq \int_{-\infty}^\infty \vert f^{(k)}(y) \vert \,dy[/math]
この2つの不等式を片々加えて 1 + |2πξ|k で除することにより上記の結果が得られる。
作用素論
作用素論における部分積分の利用例の1つとして、-Δ ( Δ は ラプラス作用素) が L2 において正値作用素であるということが挙げられる(Lp 空間を参照)。
f が滑らかでコンパクトな台を持つならば、部分積分を用いることにより以下の結果を得る。
- [math]\begin{align} \langle -\Delta f, f \rangle_{L^2} &= -\int_{-\infty}^\infty f''(x)\overline{f(x)}\,dx \\ &=-\left[f'(x)\overline{f(x)}\right]_{-\infty}^\infty + \int_{-\infty}^\infty f'(x)\overline{f'(x)}\,dx \\ &=\int_{-\infty}^\infty \vert f'(x)\vert^2\,dx \geq 0. \end{align}[/math]
その他の応用
- 弱微分、およびシュワルツの分布関数の定式化
- スツルム=リウヴィル理論における境界値問題
- 変分法によるオイラー=ラグランジュ方程式の導出
脚注
- ↑ Konrad Königsberger: Analysis 1. Springer-Verlag, Berlin u. a., 2004, ISBN 3-540-41282-4, 202.
- ↑ Yvonne Stry: Mathematik kompakt: für Ingenieure und Informatiker. 3., bearb. Auflage, Springer-Verlag, 2010, ISBN 3642111912, 314.
- ↑ Otto Forster: Analysis Band 1: Differential- und Integralrechnung einer Veränderlichen. Vieweg-Verlag, 8. Aufl. 2006, ISBN 3-528-67224-2, 210.
参考文献
- Evans, Lawrence C. (1998). Partial Differential Equations. Providence, Rhode Island: American Mathematical Society. ISBN 0-8218-0772-2.
- Arbogast, Todd; Jerry Bona (2005). Methods of Applied Mathematics (Pdu).
関連項目
外部リンク
- Weisstein, Eric W. “Integration by Parts”. MathWorld(英語). Template:Cite webの呼び出しエラー:引数 accessdate は必須です。