制御と振動の数学/第一類/演算子法の誕生/ある思いつきのソースを表示

==§1==

まず 1 階の微分方程式

{{制御と振動の数学/equation|<math>\frac{dx}{dt} + ax = f(t)</math>}}

を解くことを考えよう．ここに <math>a</math> は定数とする．いま，

{{制御と振動の数学/equation|<math>p := \frac{d}{dt}</math><ref>
記号 <math>:=</math> は左辺を右辺で定義することを表す．
</ref>
|tag=(1.1)|label=eq:1.1}}

とおくと，式[[制御と振動の数学/第一類/演算子法の誕生/ある思いつき#eq1.1|(1.1)]]は，

{{制御と振動の数学/equation|<math>px + ax = f(t)</math>}}

すなわち

{{制御と振動の数学/equation|<math>(p + a)x = f(t)</math>}}

と書くことことができる．ここで <math>p</math> を普通の数のように考えることができるものと仮定して，<math>p + a</math> で割ると，

{{制御と振動の数学/equation|<math>x = \frac{f(t)}{p + a}</math>}}

となる．そこでなんらかの方法で割り算が実行できて，簡単に解が求まれば好都合である．しばらく厳密な論理性の追求はあとまわしにし，
正しい解が得られればよいと割り切って，自由奔放に考えていこう．

<!-- ex:001:start-->
<div id="ex:1">
<strong>例1</strong><math>\quad</math>
{{制御と振動の数学/equation|<math>\frac{dx}{dt} + x = t^2 + 3t</math>|tag=(1.2)|label=eq:1.2}}
<!-- ex:001:end-->

この場合，<math>p</math> を用いて <math>x</math> を表すと，

{{制御と振動の数学/equation|<math>x = \frac{t^2 + 3t}{p + 1}</math>}}

となる．<math>p</math> が微分であることを念頭において，小学校以来おなじみの割り算を筆算で実行しよう．


<math>
\begin{array}{lcl}
       &    &  t^2  & t  & -1 \\
\hline
1 + p  & )  &  t^2  & 3t & \\
       &    &  t^2  & 2t & \\
\hline
       &    &       & t  & \\
       &    &       & t  &  +1 \\
\hline
       &    &       &    &  -1 \\
       &    &       &    &  -1 \\
\hline
       &    &       &    &   0 \\
\end{array}
</math>

{{制御と振動の数学/equation|<math>(t^2 + 3t) / (1 + p) = (t^2 + t -1) \cdots\cdots 0</math>}}

と割り切れ，商
{{制御と振動の数学/equation|<math>x = t^2 + t -1</math>|tag=(1.3)|label=eq:1.3}}

を得る．上の割り算の手順は次の通りである．まず <math>t^2</math> を <math>1</math> で割ると <math>t^2</math> が立つ．
そこで，<math>t^2</math> に <math>1 + p</math> をかけて，

{{制御と振動の数学/equation|<math>(1 + p)t^2 = t^2 + 2t</math>}}

と計算する．<math>p</math> は微分だから、<math>p\cdot t^2 = 2t</math>．被除数 <math>t^2 + 3t</math> から上式をを引くと，

{{制御と振動の数学/equation|<math>(t^2 + 3t) - (t^2 + 2t) = t</math>}}

となる．次に <math>t</math> を <math>1</math> で割ると <math>t</math> が立つ．これに <math>(1 + p)</math> をかけ，<math>t + 1</math> を得る．
<math>t</math> からこれを引いて <math>-1</math> を得る．<math>\cdots</math> という通常の手順に従えばよいのである．
このようにして求めた商，式[[制御と振動の数学/第一類/演算子法の誕生/ある思いつき#eq1.3|(1.3)]]が
式[[制御と振動の数学/第一類/演算子法の誕生/ある思いつき#eq1.2|(1.2)]]の解である保証はない．それゆえ吟味が必要である．

{{制御と振動の数学/equation|<math>x = t^2 + t - 1</math>}}
{{制御と振動の数学/equation|<math>x' = 2t + 1</math>}}
{{制御と振動の数学/equation|<math>x + x'= t^2 + 3t</math>}}

よって間違いなく，式[[制御と振動の数学/第一類/演算子法の誕生/ある思いつき#eq1.2|(1.2)]]の解である．初期値は <math>x(0) = -1</math> となっている．

<!-- ex:002:start-->
<div id="ex:2">
<strong>例2</strong><math>\quad</math>
{{制御と振動の数学/equation|<math>2\frac{d^2 x}{dt^2} + 3\frac{dx}{dt} + x = t^3 + 9t^2</math>}}

を上の方法で解き，それが解であることを確かめよ．

<strong>解答例</strong><math>\quad</math>

<math>
\begin{array}{lcl}
              &    &  t^3  &      & -12t & +36 \\
\hline
1 + 3p + 2p^2 & )  &  t^3  & +9t^2 & \\
              &    &  t^3  & +9t   & +12t \\
\hline
              &    &       &       & -12t \\
              &    &       &       & -12t & -36 \\
\hline
              &    &       &       &      & +36 \\
              &    &       &       &      & +36 \\
\hline
              &    &       &       &      &   0 \\
\end{array}
</math>


<math>(t^3 + 9t^2) / (1 + 3p + 2p^2) = t^3 \cdots\cdots -12t</math>

<math>-12t / (1 + 3p + 2p^2) = -12t \cdots\cdots 36</math>

<math>36 / (1 + 3p + 2p^2) = 36 \cdots\cdots 0</math>

<math>\therefore x = t^3 -12t + 36</math><br /><br />
次に検算を実施する。<br />
<math>x = t^3 - 12t + 36</math> のとき、<br />
:<math>x' = 3t^2 -12</math>
:<math>x'' = 6t</math>
:<math>\therefore 2x'' + 3x' + x = 2\cdot 6t + 3(3t^2 - 12) + (t^3 - 12t + 36)</math>
:<math>= {\cancel{12t}} + 9t^2 {\cancel{-36}} + t^3 {\cancel{-12t}} +{\cancel{36}}</math>
:<math>= t^3 + 9t^2</math>
すなわち <math>x = t^3 -12t + 36</math><br />は<br /><math>2\frac{d^2 x}{dt^2} + 3\frac{dx}{dt} + x = t^3 + 9t^2</math><br /> の解のひとつ．


<!-- ex:002:end -->

==§2==

[[制御と振動の数学/第一類/演算子法の誕生/ある思いつき#ex:1|例1]]では，除数が <math>1 + p</math> であった．定数 <math>1</math> があるため，割り算が簡単になった．
そこで定数項の欠けた次の例を考えてみよう．

<!-- ex:003:start-->
<div id="ex:3">
<strong>例3</strong><math>\quad</math>
{{制御と振動の数学/equation|<math>\frac{d^2x}{dt^2} + \frac{dx}{dt} = t^2 + 3t</math>}}
<!-- ex:003:end-->


<math>p = \frac{d}{dt}, \quad p^2 = \frac{d^2}{dt^2}</math> とおけば，

{{制御と振動の数学/equation|<math>(p^2 + p)x = t^2 + 3t</math>}}

となるから，割り算を実行すると，


<math>
\begin{array}{lcl}
         &    &  \frac{t^3}{3}  & \frac{t^2}{2} & -t     \\
\hline
p + p^2  & )  &  t^2            & 3t            &        \\
         &    &  t^2            & 2t            &        \\
\hline
         &    &                 & t             &        \\
         &    &                 & t             &  \ 1   \\
\hline
         &    &                 &               &  -1    \\
         &    &                 &               &  -1    \\
\hline
         &    &                 &               &  \ \ 0 \\
\end{array}
</math>


<math>(t^2 + 3t) / (p + p^2) = \frac{t^3}{3} \cdots\cdots t</math>

<math>t / (p + p^2) = \frac{t^2}{2} \cdots\cdots -1</math>

<math>-1 / (p + p^2) = -t \cdots\cdots 0</math>

<math>x = \frac{t^3}{3} + \frac{t^2}{2} - t</math>

となる．手順は次のとおりである．<math>t^2</math> を <math>p</math> で割るというのは，<math>p</math> を掛けると <math>t^2</math> となる式を求めるということであるから，
<math>p</math> が微分演算であることを思い出すと，答えとして <math>\frac{t^3}{3}</math> を得る．次に，

{{制御と振動の数学/equation|<math>(p + p^2)\frac{t^3}{3} = t^2 + 2t</math>}}

を被除数 <math>t^2 + 3t</math> から引き，<math>t</math> を得る．<math>t</math> を <math>p</math> で割ると，上と同様に考えて，<math>\frac{t^2}{2}</math> となる．
以下同様．

<strong>験算</strong>

{{制御と振動の数学/equation|<math>x = \frac{t^3}{3} + \frac{t^2}{2} - t</math>}}
{{制御と振動の数学/equation|<math>x' = t^2 + t - 1</math>}}
{{制御と振動の数学/equation|<math>x'' = 2t + 1</math>}}
{{制御と振動の数学/equation|<math>x' + x'' = t^2 + 3t</math>}}

今度もうまくいった．初期値は <math>x(0) = 0, x'(0) = -1</math> である．

==§3==

[[制御と振動の数学/第一類/演算子法の誕生/ある思いつき#ex:3|例3]]の計算の中で，<math>p</math> で割るという演算を抜き出してみると，

{{制御と振動の数学/equation|<math>\frac{1}{p}\cdot t^2 = \frac{t^3}{3}</math>}}
{{制御と振動の数学/equation|<math>\frac{1}{p}\cdot t = \frac{t^2}{2}</math>}}
{{制御と振動の数学/equation|<math>\cdots\cdots</math>}}

である．これらの事実から次のことが分かる．<math>p</math> で割るということは，積分する（あるいは原始関数を求める）ということにほかならぬ．つまり，

{{制御と振動の数学/equation|<math>\frac{1}{p} := \int dt</math>|tag=(1.4)|label=eq:1.4}}

を意味することが分かる．

このことを念頭において，[[制御と振動の数学/第一類/演算子法の誕生/ある思いつき#ex:3|例3]]の計算を次のように少し修正してみよう．
積分を先にすませておくのである．

{{制御と振動の数学/equation|<math>x = \frac{t^2 + 3t}{p + p^2} = \frac{1}{1 + p} \left\{ \frac{1}{p} (t^2 + 3t) \right\}</math>}}
{{制御と振動の数学/equation|<math>= \frac{1}{1 + p} \cdot \left\{ \frac{t^3}{3} + \frac{3t^2}{2} \right\}</math>}}

と変形しておいてから，割り算を実行する．

<math>
\begin{array}{lcl}
       &    &  \frac{t^3}{3}  & \frac{t^2}{2}  & -t     &  \ \ 1   \\
\hline
1 + p  & )  &  \frac{t^3}{3}  & \frac{3t^2}{2} &        &          \\
       &    &  t^2            & t^2            &        &          \\
\hline
       &    &                 & \frac{t^2}{2}  &        &          \\
       &    &                 & \frac{t^2}{2}  &  \ \ t &          \\
\hline
       &    &                 &               &  -t     &          \\
       &    &                 &               &  -t     & -1       \\
\hline
       &    &                 &               &         &  \ \ \ 1 \\
       &    &                 &               &         &  \ \ \ 1 \\
\hline
       &    &                 &               &         &  \ \ \ 0 
\end{array}
</math>



<math>\left( \frac{t^3}{3} + \frac{3t^2}{2} \right) / (1 + p) = \frac{t^3}{3} \cdots\cdots \frac{t^2}{2}</math>

<math>\left( \frac{t^2}{2} \right) / (1 + p) = \frac{t^2}{2} \cdots\cdots -t</math>

<math>-t / (1 + p) = -t \cdots\cdots 1</math>

<math>1 / (1 + p) = 1 \cdots\cdots 0</math>

よって，

{{制御と振動の数学/equation|<math>x = \frac{t^3}{3} + \frac{t^2}{2} - t + 1</math>}}

を得る．これは[[制御と振動の数学/第一類/演算子法の誕生/ある思いつき#ex:3|前に述べたもの]]と定数 <math>+1</math> だけ異なっているが，やはり解であることは験算するまでもなく明らかである．
ただし，初期値は <math>x(0) = 1, x'(0) = -1</math> である．定数の差は不定積分に伴う積分定数のとり方に依存する．同じ解が欲しければ，

{{制御と振動の数学/equation|<math>\frac{1}{p}(t^2 + 3t) = \frac{t^3}{3} + \frac{3t^2}{2} - 1</math>}}

としておけばよい．積分定数の選び方によって解が異なるということは，欠陥ではなく長所である．色々と異なる初期値を持った解が得られることを示唆しているからである．
我々の前途に光明を投じてくれているのである．

<!-- ex:004:start-->
<div id="ex4">
<strong>例4</strong><math>\quad</math>
{{制御と振動の数学/equation|<math>a_n\frac{d^n x}{dt^n} + a_{n-1}\frac{d^{n-1} x}{dt^{n-1}} + a_{n-2}\frac{d^{n-2} x}{dt^{n-2}} + \cdots + a_0x = f(t)</math>}}

を上述の方法で解いたとき，常に割り切れ，その積が解となっていることを示せ．
<!-- ex:004:end-->



<references />

[[カテゴリ:微分方程式]]