部分積分って何だ?

部分積分って何だ?

皆さまお疲れ様です。

前々回のラプラス変換の記事で繰越となっていました『部分積分』について、今回は書いていきます。なお、前々回分はコチラです。

新制度以降の23年分の過渡現象をラプラス変換で解き終わりました。過渡現象問題について思うところを書いていきます。

部分積分とは

2つの関数を一方は積分しつつ、もう一方は微分するという少し変則的な微積分の公式です。

\begin{align}\int^a_b f(x)g(x)dx=[f(x)G(x)]^a_b-\int^a_b f'(x)G(x)dx\\ \end{align}

$G(x)$は$g(x)$を積分した『原始関数』というものです。$g(x)=x$とすると、$G(x)=\frac{1}{2}x^2$です。

導出方法

$f(x)$と$g(x)$のどちらを積分して微分すれば良いかが分かなくなることもありますので、簡単に導出する方法を紹介します。

ベースとなるのは2関数微分です。

\begin{align}(f\cdot g)’=f’\cdot g+f\cdot g’\\ \end{align}

簡単のために$(x)$を省略しています。
\begin{align}f\cdot g’=(f\cdot g)’-f’\cdot g\\ \end{align}

と移項して、両辺を積分します。

\begin{align}\int^a_b f\cdot g’dt=\int^a_b (f\cdot g)’dt-\int^a_b f’\cdot gdt\\ \end{align}

右辺第1項では積分と微分が打ち消し合います。また、$g$→$g’$の微分関係を$G$→$g$に1段階シフトさせると、

\begin{align}\int^a_b f\cdot gdx=[f\cdot G]^a_b-\int^a_b f’\cdot Gdx\\ \end{align}

となります。

どういうときに使うのか

電験ではなかなか使う機会が少ない公式です。積分しても形が変わらない$e^x$が入っている関数で、例えば$xe^x$を積分するときに使用します。

\begin{align}\int^a_b xe^xdx&=[xe^x]^a_b-\int^a_b e^xdx\\
&=(ae^a-be^b)-[e^x]^a_b\\
&=(ae^a-be^b)-(e^a-e^b)\\
&=(a-1)e^a-(b-1)e^b\\
\end{align}

というようにです。

ラプラス変換における部分積分

ラプラス変換の問題を解いていて部分積分を使うことになるとは思いませんでした。新制度(平成7年度)以降の23年間で21回過渡現象の問題が出題された中で、部分積分を使用したのは平成16年度の1回のみでした。

グラフの概形を5種類答えるだけの形式で、最後のグラフの概形を答えるためだけに部分積分以外の計算も含めて、延々とルーズリーフの表裏2枚くらいにわたって計算しました。過渡現象以外の他の問題を考えると明らかに捨てた方が良いので、本番であれば部分積分などせず、適当に答えを選んでおしまいです。

妥協も必要です笑

ラプラス変換の計算の中でどのように登場してきたかを話し始めると長くなりますので省略しますが、計算の最後の最後で

\begin{align}\int^t_0 e^{a\tau}sinb\tau d\tau\\ \end{align}

を計算する必要が出てきました。

ラプラス変換の解答集ではきちんと書いています。

これからが本当の地獄だ…

実は部分積分をよく使う関数として、$e^x$と三角関数を組み合わせた関数があります。今回がそのパターンです。計算を延々としていきますがお付き合い下さい。

\begin{align}&\int^t_0 e^{a\tau}sinb\tau d\tau\\
&=[\frac{1}{a}e^{a\tau}sinb\tau]^t_0-\int^t_0 \frac{1}{a}e^{a\tau}bcosb\tau d\tau\\
&=\frac{1}{a}e^{at}sinbt-([\frac{1}{a^2}e^{a\tau}bcosb\tau]^t_0+\int^t_0 \frac{1}{a^2}e^{a\tau}b^2sinb\tau d\tau)\\
&=\frac{1}{a}e^{at}sinbt-\frac{b}{a^2}(e^{at}cosbt-1)-\frac{b^2}{a^2}\int^t_0 e^{a\tau}sinb\tau d\tau\\
\end{align}

$\int^t_0 e^{a\tau}sinb\tau d\tau$からスタートしてまた$\int^t_0 e^{a\tau}sinb\tau d\tau$が登場してきましたので、等式のように扱って以下のような処理をします。

\begin{align}\int^t_0 e^{a\tau}sinb\tau d\tau&=\frac{1}{a}e^{at}sinbt-\frac{b}{a^2}(e^{at}cosbt-1)-\frac{b^2}{a^2}\int^t_0 e^{a\tau}sinb\tau d\tau\\
(1+\frac{b^2}{a^2})\int^t_0 e^{a\tau}sinb\tau d\tau&=\frac{1}{a}e^{at}sinbt-\frac{b}{a^2}(e^{at}cosbt-1)\\
\int^t_0 e^{a\tau}sinb\tau d\tau&=\frac{1}{1+\frac{b^2}{a^2}}[\frac{1}{a}e^{at}sinbt-\frac{b}{a^2}(e^{at}cosbt-1)]\\
&=\frac{1}{a(a^2+b^2)}[ae^{at}sinbt-b(e^{at}cosbt-1)]\\
\end{align}

こういう計算をしていると、度々途中で何をしているのか分からなくなります笑 そういうといきはもう機械的に計算するしかありません。

まとめ

面倒な計算をしてきましたが、部分積分は電験では知っておく必要が無いかと思います。もちろん電験1種でも使用した記憶はありません。

これを覚えるくらいならば、部分分数分解とかを完璧にした方がまだ実用的ですね笑

それでは次回!