制御と振動の数学/Laplace 変換/有理関数の原像/部分分数分解のソースを表示

演算子法を用いて微分方程式を解く際に，一番わずらわしい作業は，有理関数（分数式）の原像を求めることである．
今まで何回か述べてきたように，分数式をより簡単な分数式の和に分解しておいてから原像を求めるのである．
たとえば，
{{制御と振動の数学/equation|<math>\frac{1}{(s-1)(s+2)}=\frac{1}{3}\left( \frac{1}{s-1} - \frac{1}{s+2} \right) \sqsubset \frac{1}{3} (e^t - e^{-2t})</math>}}
のようになる．この操作を[[w:%E9%83%A8%E5%88%86%E5%88%86%E6%95%B0%E5%88%86%E8%A7%A3|部分分数分解]]という．この例は簡単であるが，次の分数式，
{{制御と振動の数学/equation|<math>R(s)=\frac{4}{(s^2 + 1)^2(s-1)^3}</math>|tag=(2.29)|label=eq:2.29}}
の場合はどうなるのであろうか．結論を先にいえば，係数を実数の範囲に制限しておくと，次のように分解できる．
{{制御と振動の数学/equation|<math>R(s)=\frac{A}{(s-1)^3}+\frac{B}{(s-1)^2}+\frac{C}{s-1}+\frac{Ds+E}{(s^2+1)^2}+\frac{Fs+G}{s^2+1}</math>}}
それをみるのは比較的簡単である．まず，
{{制御と振動の数学/equation|<math>R_1(s)=R(s)-\frac{A}{(s-1)^3}</math>}}
とおいて，<math>R_1(s)</math> が <math>R(s)</math> よりも簡単な分数式となるように <math>A</math> を決めればよい．
ここで「簡単な」とは分母の多項式が簡単になるという意味である．
{{制御と振動の数学/equation|<math>R_1(s)=\frac{4}{(s^2+1)^2(s-1)^3}-\frac{A}{(s-3)^3}=\frac{4-A(s^2+1)2}{(s^2+1)^2(s-3)^3}</math>}}
において分子が <math>(s-1)</math> を因数として持つように <math>A</math> を定めれば，分母分子が約されて簡略化が達成できる．
そのためには[[w:%E5%9B%A0%E6%95%B0%E5%AE%9A%E7%90%86|因数定理]]を用いて <math>A=1</math> と定めればよい<ref>
<math>(s-1)</math> で割り切ればよいから，<math>4-A(s^2+1)^2\bigg|_{s=1} = 4-A\cdot4 = 0 \therefore A = 1</math>．
</ref>．このとき，
{{制御と振動の数学/equation|<math>R_1(s)=\frac{-s^3-s^2-3s-3}{(s^2+1)^2(s-1)^2}</math>}}
となる．次に，
{{制御と振動の数学/equation|<math>R_2(s)=R_1(s) - \frac{B}{(s-1)^2}</math>}}
とおいて、上と同様な操作を繰り返すと <math>B = -2</math> と定まり<ref>
<math>-s^3-s^2-3s-3 - B(s^2+1)^2\bigg|_{s=1} = -8 - B\cdot 4 = 0 \therefore B = -2</math>．
</ref>，
{{制御と振動の数学/equation|<math>R_2(s)=\frac{2s^3 + s^2 + 4s + 1}{(s^2 + 1)^2(s-1)}</math>}}
となる．同様にして，
{{制御と振動の数学/equation|<math>R_3(s)=R_2(s) - \frac{C}{s-1}</math>}}
より，<math>C = 2</math><ref>
<math>2s^3+s^2+4s+1-C(s^2+1)^2\bigg|_{s=1} = 8 - 4C = 0 \therefore C = 2</math>． 
</ref>．よって，
{{制御と振動の数学/equation|<math>R_3(s) = \frac{-2x^3-3s+1}{(s^2+1)^2}</math>}}
と，遂には分母から因数 <math>(s-1)</math> が消えてしまう．次に最後の式の分子を
{{制御と振動の数学/equation|<math>-2s^3-3s+1=-2s(s^2+1)-s+1</math>}}
と変形し，
{{制御と振動の数学/equation|<math>R_3(s)=\frac{-s+1}{(s^2+1)^2}-\frac{2s}{s^2+1}</math>}}
を得る．よって [[制御と振動の数学/Laplace 変換/有理関数の原像/部分分数分解#eq:2.29|式(2.29)]] の部分分数展開が完成する．すなわち，
{{制御と振動の数学/equation|<math>\frac{4}{(s^2+1)^2(s-1)^3}=\frac{1}{(s-1)^3}-\frac{2}{(s-1)^2}+\frac{2}{s-1}-\frac{s-1}{(s^2+1)^2}-\frac{2s}{s^2+1}</math>}}

この例から，一般の場合を予測することは，そう難しくないであろう．一般の有理関数
{{制御と振動の数学/equation|<math>R(s):=\frac{b_0s^m + b_1s^{m-1} + b_2s^{m-2} + \cdots + bm}{s^n + a_1s^{n-1} + a_2s^{n-2} + \cdots + a_n}\quad (m<n)</math>|tag=(2.30)|label=eq:2.30}}
の場合の結果を述べる．

<div id="bubun-bunsu">
<strong>[部分分数定理]</strong>

<math>R(s)</math> を実係数の真の分数式<ref>
真の分数式とは，分子の次数が分母の次数より低い分数式（有理関数）のことをいう．
</ref>、<math>f(s)</math> をその分母とする．
{{制御と振動の数学/equation|<math>f(s):=s^n + a_1s^{n-1} + a_2s^{n-2} + \cdots + a_n</math>}}
は実数の範囲で必ず既約な 1 次式と 2 次式との積に因数分解できる<ref>
1 変数の多項式による方程式 <math>P(x)=\sum_{k}a_kx^k = 0, \quad a_k \in R</math> が虚根 <math>\alpha</math> を持つ，すなわち <math>P(\alpha) = 0</math> のとき，<br />
:<math>P(\bar{\alpha})=0, \quad</math> <math>\bar{\alpha}</math> は <math>\alpha</math> の共役複素数．…①<br />
①を示す．まず複素数 <math>\alpha, \beta</math> について，<br />
:<math>\overline{\alpha + \beta} = \bar{\alpha} + \bar{\beta}</math>…②
なぜならば，<math>\alpha = a + bi</math> とおくとき，<math>\bar{\alpha} = a - bi</math> で，
同じく <math>\beta = c + di</math> とおくと <math>\bar{\beta} = c - di</math>，<br />
:<math>\overline{\alpha + \beta} = \overline{(a + c) + (b + d)i}</math>
:<math>=(a + c) - (b + d)i = (a - bi) + (c - di) = \bar{\alpha} + \bar{\beta}</math>．
次に，<br />
:<math>\overline{\alpha \cdot \beta} = \bar{\alpha} \cdot \bar{\beta}</math>…③
なぜならば、同じく上記 <math>\alpha, \beta</math> について，
:<math>\overline{\alpha\cdot\beta} = \overline{(ac-bd)+(ad+bc)i}</math>
:<math>=(ac-bd)-(ad+bc)i = (a - bi)(c-di) = \bar{\alpha}\cdot\bar{\beta}</math>．
よって，<math>P(\alpha) = 0</math>より，
:<math>\overline{P(\alpha)} = \overline{0} = 0</math>
<math>\because</math> 実数 <math>r</math> の共役複素数は <math>r</math> そのものに等しい．…④<br />
②より共役複素数をとる操作と足し算は順序を入れ替えることができて，<br />
:<math>\overline{P(\alpha)} = \overline{\sum_k a_k\alpha^k} = \sum_k \overline{a_k\alpha^k} = 0</math>
③より共役複素数をとる操作と掛け算は順序を入れ替えることができて，<br />
:<math>\sum_k \overline{a_k\alpha^k} = \sum_k \overline{a_k}\cdot\overline{\alpha^k} = \sum_k a_k\cdot\overline{\alpha^k} = 0</math>
最後の変形は④による．
さらに③を適用して，<br />
:<math>\sum_k a_k\cdot\overline{\alpha^k} = \sum_k a_k\cdot(\overline{\alpha})^k = 0</math>
すなわち，<br />
:<math>\sum_k a_k\cdot(\overline{\alpha})^k = P(\bar{\alpha}) = 0</math>
ここに①を示せた.<br />
すなわち <math>P(x) = 0</math> が虚根を持つのなら，虚根 <math>\alpha</math> と <math>\bar{\alpha}</math> を同時に持つから，
:<math>P(x) = (x-\alpha)(x-\bar{\alpha})Q(x)</math>
と表せて，<math>(x-\alpha)(x-\bar{\alpha}) = x^2 - (\alpha + \bar{\alpha})x + \alpha\bar{\alpha}</math><br />
:<math>\alpha + \bar{\alpha} = a + bi + a - bi = 2a \in \mathbf{R}</math>．
:<math>\alpha\cdot\bar{\alpha} = (a + bi)(a - bi) = a^2 + b^2 \in \mathbf{R}</math>．
すなわち <math>(x-\alpha)(x-\bar{\alpha})</math> は実係数の二次式であり，<math>P(x)=0</math> の根は複素数空間内に存在することを前提として，<math>P(x)</math> は実係数の範囲で高々二次式までに因数分解される．<br />
 
</ref>から，それを，
{{制御と振動の数学/equation|<math>f(s) = \prod_{i=1}^{\mu}(s-\alpha_i)^{l_i}\cdot\prod_{j=1}^{\nu}(s^2+p_js+q_j)^{m_j}</math>}}
ここに，
{{制御と振動の数学/equation|<math>\sum_{i=1}^{\mu}l_i + 2\sum_{j=1}^{\nu}m_j = n</math>}}
とすれば，<math>R(s)</math> は，
{{制御と振動の数学/equation|<math>R(s)=\sum_{i=1}^{\mu}\left\{ \frac{A_{i1}}{s-\alpha_i} + \frac{A_{i2}}{(s-\alpha_i)^2} +\cdots+ \frac{A_{il_i}}{(s-\alpha_i)^{l_i}}\right\} + \sum_{j=1}^{\nu} \left\{ \frac{B_{j1}s + C_{j1}}{s^2 + p_js + q_j} + \frac{B_{j2}s+C_{j2}}{(s^2+p_js+q_j)^2} +\cdots+ \frac{B_{jm_{j}}s + C_{jm_j}}{(s^2+p_js+q_j)^{m_j}}\right\}</math>|tag=(2.31)|label=eq:2.31}}
と展開できる．ここに <math>A, B, C</math> 等は実定数である．これらの定数の個数はちょうど <math>n</math> で，<math>n</math> は分母の多項式 <math>f(s)</math> の次数に等しい．なお，この分解は一意に定まる．

<math>\diamondsuit</math>


この定理の証明は色々あるが，次のものが標準的である．詳細は [[制御と振動の数学/第一類/付録・部分分数展開|付録]] にゆずるが，ここでその概要を述べておこう．証明は 2 段に分かれる．

<div id="bubun-bunsu1">
<strong>[第一分解定理]</strong>

<math>f(s)</math> と <math>g(s)</math> を多項式とする．<math>f(s)</math> が互いに共通因子を持たない <math>f_1(s)</math> と <math>f_2(s)</math> の積であれば，
{{制御と振動の数学/equation|<math>\frac{g(s)}{f(s)} = \frac{g_1(s)}{f_1(s)} + \frac{g_2(s)}{f_2(s)}</math>}}
と分解される．左辺が真の分数式なら，右辺の二つの分数式もそれぞれ真の分数式とすることができる．

<math>\diamondsuit</math>

この定理を反復適用すれば，ひとまずは，
{{制御と振動の数学/equation|<math>R(s) = \sum_{i = 1}^{\mu}\frac{g_i(s)}{(s-\alpha_i)^{l_i}} + \sum_{j=1}^{\nu}\frac{h_j(s)}{(s^2 + p_js + q_j)^{m_j}}</math>|tag=(2.32)|label=eq:2.32}}
と分解できることは明らかである．

<strong>例</strong>

{{制御と振動の数学/equation|<math>\frac{4}{(s^2 + 1)^2(s-1)^3} = \frac{2s^2 - 6s + 5}{(s-1)^3} + \frac{-2s^3-3s+1}{(s^2+1)^2}</math>}}
第 2 段目の分解は次の定理による．

<div id="bubun-bunsu2">
<strong>[第二分解定理]</strong>

真の分数式 <math>\frac{g(s)}{[u(s)]^m}</math> は
{{制御と振動の数学/equation|<math>\frac{g(s)}{[u(s)]^m} = \frac{h_1(s)}{u(s)} + \frac{h_2(s)}{[u(s)]^2} + \frac{h_3(s)}{[u(s)]^3} + \cdots + \frac{h_m(s)}{[u(s)]^m}</math>}}
ここに，<math>u(s)</math> の次数 > <math>h_i(s)</math> の次数 <math>(i = 1, 2, \cdots, m)</math>
の形に分解できる．

<math>\diamondsuit</math>

式 [[制御と振動の数学/Laplace 変換/有理関数の原像/部分分数分解#eq:2.32|(2.32)]] の各項にこの定理を適用すれば，求める結果 [[制御と振動の数学/Laplace 変換/有理関数の原像/部分分数分解#eq:2.31|(2.31)]] を得る．

<strong>例</strong>

{{制御と振動の数学/equation|<math>\frac{2s^2-6s+5}{(s-1)^3} = \frac{2}{s-1} - \frac{2}{(s-1)^2} + \frac{1}{(s-1)^3}</math>}}
{{制御と振動の数学/equation|<math>\frac{-2s^3-3s+1}{(s^2+1)^2} = \frac{-2s}{s^2+1}+\frac{-s+1}{(s^2+1)^2}</math>}}
部分分数分解の原理は上述のとおりであるが，実用上は次の手法により求める．


<!-- ex:051:start-->
<div id="ex:51">
<strong>例51</strong><math>\quad</math>
{{制御と振動の数学/equation|<math>\frac{2}{s^4-2s^3+2s^2-2s+1}</math>}}
を部分分数に分解せよ．

<strong>解</strong>

分母を因数分解すると <math>(s^2+1)(s-1)^2</math> となるから，
{{制御と振動の数学/equation|<math>\frac{2}{(s-1)^2(s^2+1)} = \frac{A}{s-1}+\frac{B}{(s-1)^2}+\frac{Cs+D}{s^2 + 1}</math>}}
の形に分解できる．分母を払うと，
{{制御と振動の数学/equation|<math>2=A(s-1)(s^2+1)+B(s^2+1)+(Cs+D)(s-1)^2</math>}}
となるから，
{{制御と振動の数学/equation|<math>s=1</math> とおいて ; <math>2=2B, \therefore B = 1</math>}}
{{制御と振動の数学/equation|<math>s=0</math> とおくと ; <math>2=-A+B+D</math>}}
{{制御と振動の数学/equation|<math>s</math> の 2 次の係数を等置して ; <math>0 = -A + B + D -2C, \therefore C = 1</math>}}
{{制御と振動の数学/equation|<math>s</math> の 3 次の係数を等置して ; <math>0=A + C, \therefore A = -1</math>}}
{{制御と振動の数学/equation|最後に <math>D=0</math>}}
よって，
{{制御と振動の数学/equation|<math>\frac{2}{(s-1)^2(s^2+1)} = \frac{-1}{s-1}+\frac{1}{(s-1)^2}+\frac{s}{s^2 + 1}</math>}}
と分解される．
<math>\diamondsuit</math>
<!-- ex:051:end-->


<!-- ex:052:start-->
<div id="ex:52">
<strong>例52</strong><math>\quad</math>

次の分数式を部分分数に分解して原像を求めよ．
{{制御と振動の数学/equation|<math>\frac{s^2+2s}{(s^2+4s+5)^2}</math>}}

<strong>解答例</strong>

:<math>\frac{s^2+2s}{s^2 + 4s + 5} = \frac{As + B}{s^2 + 4s + 5} + \frac{Cs + D}{(s^2+4s+5)^2}</math>
:<math>s^2+2s = (s^2 + 4s + 5)(As + B) + Cs + D</math>
<math>\begin{cases}
A &= 0 \\
B + 4A &= 1 \\
4B + 5A + C &= 2 \\
5B + D &=0
\end{cases}</math><br />
よって，<br />
:<math>A =0</math>
:<math>\therefore B = 1</math>
:<math>4 + C = 2, \therefore C = -2</math>
:<math>5 + D = 0, \therefore D = -5</math>
:<math>\therefore \frac{s^2+2s}{s^2 + 4s + 5} = \frac{1}{s^2 + 4s + 5} + \frac{-2s-5}{(s^2+4s+5)^2}</math>
:<math>x \sqsupset X = \frac{1}{s^2 + 4s + 5} + \frac{-2s-5}{(s^2+4s+5)^2}</math>
とおくと，<br />
:<math>X = \frac{1}{(s+2)^2 + 1} + \frac{-2(s+2) -1}{\{(s+2)^2 + 1\}^2}</math>
よって [[制御と振動の数学/第一類/Laplace 変換/三角関数の Laplace 変換とその応用#eq:2.21|式(2.21)]] より
:<math>x = e^{-2t} \left\{ \sin t - \frac{t}{2}\sin t - \frac{1}{2}(\sin t - t\cos t) \right\}</math>
:<math>= \frac{e^{-2t}}{2} \left( \sin t - t\sin t + t\cos t \right)</math>

<math>\diamondsuit</math>
<!-- ex:052:end-->


<!-- ex:053:start-->
<div id="ex:53">
<strong>例53</strong><math>\quad</math>

次の分数式を部分分数に分解して原像を求めよ．
{{制御と振動の数学/equation|<math>\frac{1}{s(3s+1)^3}</math>}}

<strong>解答例</strong>

:<math>\frac{1}{s(3s+1)^3} = \frac{A}{s} + \frac{B}{3s+1} + \frac{C}{(3s+1)^2} + \frac{D}{(3s+1)^3}</math>
:<math>1 = A(3s+1)^3 + Bs(3s+1)^2 + Cs(3s+1) + Ds</math>
<math>s=0</math> を代入して，<math>A = 1</math>，<br />
<math>s = \frac{-1}{3}</math> を代入して，<math>D = -3</math>，<br />
以上から<br />
:<math>1 = (3s+1)^3 + Bs(3s+1)^2 + Cs(3s+1) -3s</math>
<math>s</math> で微分すると，<br />
:<math>0 = 3(3s+1)^2\cdot3 + B\{ (3s+1)^2 + s\cdot 2(3s + 1)\cdot 3\} + C(6s+1)-3</math>
整理して
:<math>0 = 3(3s+1)^2\cdot 3 + B(3s+1)(9s + 1) + C(6s + 1)-3</math>…①
<math>s=\frac{-1}{3}</math> を代入して <math>C = -3</math>．<br />
①に <math>C = -3</math> を代入した後で <math>s^0 = 1</math> の係数を比較して，<br />
:<math>0 = 9 + B -3 -3, \therefore B = -3</math>
すなわち，<br />
:<math>\frac{1}{s(3s+1)^3} = \frac{1}{s} + \frac{-3}{3s+1} + \frac{-3}{(3s+1)^2} + \frac{-3}{(3s+1)^3}</math>
:<math>=\frac{1}{s} + \frac{-1}{s + \frac{1}{3}} + \frac{\frac{-1}{3}}{ \left( s + \frac{1}{3} \right)^2} + \frac{\frac{-1}{9}}{\left(s+\frac{1}{3}\right)^3}</math>…②
第 4 項については，[[制御と振動の数学/第一類/Laplace 変換/f(t) の積分および微分の Laplace 変換#eq:2.8|式(2.8)]] より<br />
:<math>\mathcal{L}\left[\frac{t^2}{2}\right] = \frac{1}{s^3}</math> よって <math>\mathcal{L}\left[-\frac{1}{18}t^2\right] = \frac{-\frac{1}{9}}{s^3}</math>
以上より②の原像は<br />
:<math>1 + e^{-\frac{t}{3}} \left( -1 -\frac{1}{3}t - \frac{1}{18}t^2 \right)</math>

<math>\diamondsuit</math>
<!-- ex:053:end-->

[[カテゴリ:ラプラス変換]]