制御と振動の数学/第一類/演算子法の誕生/演算子法の完成

§1

いまや我々は演算子法を完成させ得る段階に達した． $p$ が微分で，${\frac {1}{p}}$ が積分を表すというのであるから， 微分と積分の関係を表す次の二つの公式，

(1.10)

$\int _{0}^{t}{\frac {dx(t)}{dt}}dt=x(t)-x(0)$

${\frac {d}{dt}}\int _{0}^{t}x(t)dt=x(t)$

がまず頭に浮かぶ．初期値 $x(0)$ を含んでいる前者を採用するのが妥当であろう．

そこで，すでに述べたように，

$\int _{0}^{t}f(\tau )d\tau ={\frac {1}{p}}f(t)$

によって式 (1.10) を記号 $p$ を用いて表せば，

(1.11)

${\frac {1}{p}}x'(t)=x(t)-x(0)$

ここに

$x'(t):={\frac {dx(t)}{dt}}$

となる．いま $p$ を普通の数と同じように取り扱ってよいものとすれば，$p$ を両辺に掛けて，

(1.12)

$x'(t)=px(t)-px(0)$

なる関係式を得る．ここではもはや $p$ は微分という意味をもたない． もし $p={\frac {d}{dt}}$ ならば $px(0)=0$ となってしまうからである． もちろん $x(0)\equiv 0$ のときは，$p$ は微分と考えてなんら差し支えはないであろう． 式 (1.12) のように変形しておいたのは， 先に $x'$ を $px$ で置き換えたように，今度は $x'$ を $px-px(0)$ で置き換えようという下心である．

ここで式 (1.12) の興味深い応用を示そう． この式に $x=e^{at}$ を代入すると，

$ae^{at}=pe^{at}-p$

これを $e^{at}$ について解くと

(1.13)

$e^{at}={\frac {p}{p-a}}$

となる． この公式は実質的にはすでにこの式およびこの式にて得られている． この式を $a$ について $n$ 回微分すると，

$t^{n}e^{at}={\frac {n!\ p}{(p-a)^{n+1}}}$

よって，

(1.14)

${\frac {p}{(p-a)^{n+1}}}={\frac {t^{n}}{n!}}e^{at}$

を得る．特に $a=0$ とおくと，

(1.14a)

${\frac {1}{p^{n}}}={\frac {t^{n}}{n!}}$

となる．これもすでに得られた式 (1.9) と一致する．

例7$\quad$

${\frac {dx}{dt}}+x=e^{-t},x(0)=x_{0}$

を解け．

これは例5 で $a=-1$ とおいたものである． 前節のやり方では解決できなかったことを思い出してもらいたい． この式を $p$ で表すと，

$px-px_{0}+x={\frac {p}{p+1}}$

これを $x$ について解くと

$x={\frac {p}{p+1}}x_{0}+{\frac {p}{(p+1)^{2}}}$

ここで公式 (1.13), (1.14) を用いて $t$ の関数に戻すと，

$x(t)=e~{-t}x_{0}+te^{-t}$

例8$\quad$

これが正しい解であることを確かめよ．

解答例$\quad$

$x={\frac {p}{p+1}}x_{0}+{\frac {p}{(p+1)^{2}}}$ のとき，

$x'(t)=-x_{0}e^{-t}-te^{-t}+e^{-t}$

$\therefore x'(t)+x(t)=-x_{0}e^{-t}-te^{-t}+e^{-t}+x_{0}e^{-t}+te^{-t}=e^{-t}.$

また

$x(0)=e^{-t}x_{0}+te^{-t}|_{t=0}=x_{0}\quad (\because e^{0}=1.)$

例9$\quad$

${\frac {dx}{dt}}+x=t^{2}+3t,x(0)=x_{0}$

を解け．

式(1.12), 式(1.9)を用いて，$p$ の式に書き換えると，

$px-px_{0}+x={\frac {2}{p^{2}}}+{\frac {3}{p}}$

$x$ について解くと，

$x={\frac {p}{p+1}}x_{0}+{\frac {2+3p}{(p+1)p^{2}}}$

となる．これを時間関数に戻すにあたって，(1.13), 式(1.9) が使いやすいように、部分分数に分解すると，

$x={\frac {p}{p+1}}x_{0}+{\frac {p}{p+1}}+{\frac {2}{p^{2}}}+{\frac {1}{p}}-1$

となる．これは

$x(t)=e^{-t}x_{0}+e^{-t}+t^{2}+t-1$

に他ならない．

例10$\quad$

これが正しい解であることを確かめよ．

解答例$\quad$

$x(t)=e^{-t}x_{0}+e^{^{-}t}+t^{2}+t-1$ のとき，

$x'(t)=-x_{0}e^{-t}-e^{-t}+2t+1$．

ゆえに

$x(t)+x'(t)=(e^{-t}x_{0}+e^{^{-}t}+t^{2}+t-1)+(-x_{0}e^{-t}-e^{-t}+2t+1)=t^{2}+3t$．

また

$x(0)=e^{0}\cdot x_{0}+e^{0}-1=x_{0}$．

以上より $x(t)=e^{-t}x_{0}+e^{-t}+t^{2}+t-1$ は ${\frac {dx}{dt}}+x=t^{2}+3t,x(0)=x_{0}$ のひとつ．

例11$\quad$

${\frac {dx}{dt}}+x=e^{t},x(0)=x_{0}$ を解け．

解答例$\quad$

$px-px_{0}+x=e^{t}$

$(p+1)x=px_{0}+{\frac {p}{p-1}}$

$x={\frac {p}{p+1}}x_{0}+{\frac {p}{(p+1)(p-1)}}$

$={\frac {p}{p+1}}x_{0}+{\frac {1}{2}}\left({\frac {1}{p+1}}+{\frac {1}{p-1}}\right)$

$={\frac {p}{p+1}}x_{0}+{\frac {1}{2}}\left({\frac {1}{p+1}}-1+{\frac {1}{p-1}}+1\right)$

$={\frac {p}{p+1}}x_{0}+{\frac {1}{2}}\left({\frac {p}{p-1}}-{\frac {p}{p+1}}\right)$

$=x_{0}e^{-t}+{\frac {1}{2}}\left(e^{t}-e^{-t}\right)$

験算をする．

${\frac {dx}{dt}}=-x_{0}e^{-t}+{\frac {1}{2}}\left(e^{t}+e^{-t}\right)$

$\therefore {\frac {dx}{dt}}+x=-x_{0}e^{-t}+{\frac {1}{2}}\left(e^{t}+e^{-t}\right)$$+x_{0}e^{-t}+{\frac {1}{2}}\left(e^{t}-e^{-t}\right)$

$=e^{t}$

で与方程式を満たす．

§2

さらに 2 階微分に対する公式を導いてみよう． 式 (1.12) を 2 度用いると，

$x''=px'(t)-px'(0)=p\{px(t)-px(0)\}-px'(0)$

となる．よって

(1.15)

$x''=p^{2}x(t)-p^{2}x(0)-px'(0)$

を得る．ここでも $x(0)=x'(0)=0$ のときは $x''=p^{2}x$ となり，$p$ が微分を表すと考えてよいことを示している．

式 (1.15) において $x(t)=\sin \beta t$ とおけば，

$-\beta ^{2}\sin \beta t=p^{2}\sin \beta t-p\beta$

(1.16)

$\therefore \sin \beta t={\frac {\beta p}{p^{2}+\beta ^{2}}}$

を得る．また，

$\sin t={\frac {p}{p^{2}+1}}={\frac {1}{p}}{\frac {1}{(1+{\frac {1}{p^{2}}})}}$$={\frac {1}{p}}-{\frac {1}{p^{3}}}+{\frac {1}{p^{5}}}-{\frac {1}{p^{7}}}+\cdots$

$t-{\frac {t^{3}}{3!}}+{\frac {t^{5}}{5!}}-{\frac {t^{7}}{7!}}+\cdots$

これは $\sin t$ の Taylor 展開である．

例12$\quad$

(1)

(1.17)

$\cos \beta t={\frac {p^{2}}{p^{2}+\beta ^{2}}}$

(2)

$\cos t=1-{\frac {t^{2}}{2!}}+{\frac {t^{4}}{4!}}-{\frac {t^{6}}{6!}}+\cdots$

を示せ．

解答例$\quad$

式 (1.15) に $x(t)=\cos \beta t$ を代入して，

$-\beta ^{2}\cos \beta t=p^{2}\cos \beta t-p^{2}$

$(p^{2}+\beta ^{2})\cos \beta t=p^{2}$

$\therefore \cos \beta t={\frac {p^{2}}{p^{2}+\beta ^{2}}}$

すなわち

$\cos t={\frac {1}{1-(-{\frac {1}{p^{2}}})}}=1-{\frac {1}{p^{2}}}+{\frac {1}{p^{4}}}-{\frac {1}{p^{6}}}+\cdots$

$=1-{\frac {t^{2}}{2!}}+{\frac {t^{4}}{4!}}-{\frac {t^{6}}{6!}}+\cdots$

$\diamondsuit$

これらの結果（式(1.16)，式 (1.17) を用いると

$e^{it}={\frac {p}{p-i}}={\frac {p(p+i)}{(p^{2}+1)}}={\frac {p^{2}}{p^{2}+1}}+i{\frac {p}{p^{2}+1}}$

$=\cos t+i\sin t$

を得る．ここに $i$ は虚数単位である．また，

$\cos t={\frac {p^{2}}{p^{2}-i^{2}}}={\frac {1}{2}}\left({\frac {p}{p-i}}+{\frac {p}{p+i}}\right)$

$={\frac {e^{it}+e^{-it}}{2}}$

同様に，

$\sin t={\frac {p}{p^{2}-i^{2}}}={\frac {1}{2i}}\left({\frac {p}{p-i}}-{\frac {p}{p+i}}\right)$

$={\frac {e^{it}-e^{-it}}{2i}}$

を得る．これらは有名な Euler の公式 である．

例13$\quad$

${\frac {d^{2}x}{dt^{2}}}+x=t,x(0)=1,x'(0)=0$

を解け．

公式 (1.15) 等を用いて，$p$ の式で表すと，

$p^{2}x-p^{2}+x={\frac {1}{p}}$

$p^{2}x+x=p^{2}+{\frac {1}{p}}$

$x={\frac {p^{2}}{p^{2}+1}}+{\frac {1}{p(p^{2}+1)}}$

$={\frac {p^{2}}{p^{2}+1}}-{\frac {p}{p^{2}+1}}+{\frac {1}{p}}$

ここで公式(1.16), (1.17) を用いて $t$ の関数に戻すと，

$x(t)=\cos t-\sin t+t$

$\diamondsuit$

例14$\quad$ 例12 の解の正しいことを確かめよ．

解答例$\quad$

$x'(t)=-\sin t-\cos t+1$

$x''(t)=-\cos t+\sin t$

$x''(t)+x=-\cos t+\sin t+\cos t-\sin t+t=t,x(0)=1-0+0=1,x'(0)=-0-1+1=0$

よって

$x(t)=\cos t-\sin t+t$

は与方程式の解のひとつ．

$\diamondsuit$

例15$\quad$

${\frac {d^{2}x}{dt^{2}}}+4{\frac {dx}{dt}}+3x=5e^{2t},x(0)=2,x'(0)=1$

を解け．

解答例$\quad$

$p^{2}x-p^{2}\cdot 2-p+4px-4p\cdot 2+3x={\frac {5p}{p-2}}$

$(p^{2}+4p+3)x=2p^{2}+9p+{\frac {5p}{p-2}}$

$x={\frac {2p^{2}}{(p+1)(p+3)}}+{\frac {9p}{(p+1)(p+3)}}+{\frac {5p}{(p+1)(p+3)(p-2)}}$

第一項

${\frac {2p^{2}}{(p+1)(p+3)}}={\frac {Ap}{p+1}}+{\frac {Bp}{p+3}}$

とおいて

${\frac {Ap}{p+1}}+{\frac {Bp}{p+3}}={\frac {1}{(p+1)(p+3)}}\cdot \left\{(A+B)p^{2}+(3A+B)p\right\}$

より

$A+B=2,\quad 3A+B=0\quad \therefore A=-1,\quad B=3$

したがって

${\frac {2p^{2}}{(p+1)(p+3)}}={\frac {-p}{p+1}}+{\frac {3p}{p+3}}$

同様に第二項は

${\frac {9p}{(p+1)(p+3)}}={\frac {{\frac {9}{2}}p}{p+1}}+{\frac {-{\frac {9}{2}}p}{p+3}}$

第三項は

${\frac {5p}{(p+1)(p+3)(p-2)}}={\frac {-{\frac {5}{6}}p}{p+1}}+{\frac {{\frac {1}{2}}p}{p+3}}+{\frac {{\frac {1}{3}}p}{p-2}}$

ゆえに

$x={\frac {(-1+{\frac {9}{2}}-{\frac {5}{6}})p}{p+1}}+{\frac {(3-{\frac {9}{2}}+{\frac {1}{2}})p}{p+3}}+{\frac {{\frac {1}{3}}p}{p-2}}$

$={\frac {{\frac {8}{3}}p}{p+1}}+{\frac {-p}{p+3}}+{\frac {{\frac {1}{3}}p}{p-2}}$

$={\frac {8}{3}}e^{-t}-e^{-3t}+{\frac {1}{3}}e^{2t}$

$\diamondsuit$