4.電信方程式の一般解(ダランベールの解)

4.1 はじめに

ここからは式(3.8)の電信方程式を解いていきます。電信方程式は簡単には解かせてくれない方程式なので、変数分離法を使って単一正弦波についての解を求めたりするのが一般的です。変数分離と難しく聞こえるものの、実は一番やさしい解法で、いわゆる電気回路演算で用いられるjωなどの複素演算(フェーザ法)に対応させることと全く同じです。

他にはラプラス変換や直接シミュレーションする方法などもありますが、ここでは、求積解である「ダランベールの解」と呼ばれる一般解を求めていきます。

ちなみに、偏微分方程式は解けるほうが珍しいようで、解き方については難解な教科書が山ほど出ているほどです。

4.1 無損失の電信方程式

\begin{eqnarray} \pdiffB{2}{v}{x} &=& LC \pdiffB{2}{v}{t} + (LG+CR) \pdiffA{v}{t} + RGi \tag{3.8a} \\ \pdiffB{2}{i}{x} &=& LC \pdiffB{2}{i}{t} + (LG+CR) \pdiffA{i}{t} + RGi \tag{3.8b} \end{eqnarray}

電信方程式(3.8)は、偏微分方程式であるため、常微分方程式と異なり解くのが非常に厄介です。そのままの形では解くことができません。しかしながらR=G=0としたときに限り特別に求積法と呼ばれる解法を使うことができます。 \begin{eqnarray} \pdiffB{2}{v}{x} &=& LC \pdiffB{2}{v}{t} \tag{4.1a} \\ \pdiffB{2}{i}{x} &=& LC \pdiffB{2}{i}{t} \tag{4.1b} \end{eqnarray}

この式(4.1)は一般に無損失線路の方程式と呼ばれています。エネルギ損失の原因となる抵抗が一切ないからです。オリジナルの式(3.7)よりも随分とシンプルになりました。電圧と電流それぞれの式は全く同じ形の方程式なので、どちらか片方を解けば、直ちに他の解も分かります。ここでは電圧を解くことにします。

4.2 特別な解法(ダランベールの解法)

4.2.1 解くための準備事項

ここからは微分方程式を解くにあたって、唐突な手法がたくさん飛び出します。怪しげに思わないで、まぁそういうものだと思ってください。まず裏技その1は方程式の係数LCを以下の形に変換することから始まります。 \[ u=\frac{1}{\sqrt{LC}} \tag{4.2} \]

すると式(4.1a)の無損失電圧式は \[ \pdiffB{2}{v}{x} = \frac{1}{u^2} \pdiffB{2}{v}{t} \tag{4.3} \\ \] となります。

どうしてわざわざこんなことをするのか疑問になると思いますが、こうすると上手に解けてしまうという理由だけです。ちょっとだけタネ明かしをすると、uは波の波形の速度を表しているので、このような変換をしています。

次の裏技その2も突然に出てきます。今度は変数の変換です。ここでは二つの変数、距離xと時間tを、新しく考えた変数ξ(クシー)とη(イ−タ)に置換してしまうのです。 \begin{eqnarray} \xi &=& x-ut\\ \eta &=& x+ut \tag{4.4} \end{eqnarray}

ここで変数変換後の電圧関数v(x,t)は \[ v(x,t) \Rightarrow v(\xi,\eta) = v \left( \frac{\xi + \eta}{2} , \frac{-\xi + \eta}{2u} \right) \tag{4.5} \\ \] と表現することができます。

4.2.2 変数置換後の左辺

式(4.3)の左辺・・・すなわち距離xで偏微分している項・・・にξとηに変換された式(4.5)を代入するため、xで偏微分すると、(補足1 \begin{eqnarray} \pdiffA{}{x} v(\xi,\eta) &=& \pdiffA{v}{\xi} \pdiffA{\xi}{x} + \pdiffA{v}{\eta} \pdiffA{\eta}{x} \\ &=& \pdiffA{v}{\xi} + \pdiffA{v}{\eta} \tag{4.6} \end{eqnarray} ここでは、 \begin{eqnarray} \pdiffA{\xi}{x} &=& \pdiffA{}{x} (x-ut) = 1\\ \pdiffA{\eta}{x} &=& \pdiffA{}{x} (x+ut) = 1 \end{eqnarray} となっていることがポイントです。

式(4.6)をもう一度xで偏微分すると、(補足2 \begin{eqnarray} \pdiffB{2}{}{x} v(\xi,\eta) &=& \pdiffA{}{\xi} \left( \pdiffA{v}{\xi}+\pdiffA{v}{\eta} \right) \pdiffA{\xi}{x} + \pdiffA{}{\eta} \left( \pdiffA{v}{\xi} + \pdiffA{v}{\eta} \right) \pdiffA{\eta}{x} \\ &=& \pdiffA{}{\xi} \left( \pdiffA{v}{\xi}+\pdiffA{v}{\eta} \right) + \pdiffA{}{\eta} \left( \pdiffA{v}{\xi} + \pdiffA{v}{\eta} \right) \\ &=& \left( \pdiffB{2}{v}{\xi} + \pdiffC{v}{\eta}{\xi} \right) +\left( \pdiffC{v}{\eta}{\xi} + \pdiffB{2}{v}{\eta} \right) \\ &=& \pdiffB{2}{v}{\xi} + 2 \pdiffC{v}{\eta}{\xi} + \pdiffB{2}{v}{\eta} \\ &=& \pdiffB{2}{v}{\xi} + 2 \pdiffC{v}{\xi}{\eta} + \pdiffB{2}{v}{\eta} \end{eqnarray} . よって、式(4.3)の左辺は、 \[ \pdiffB{2}{v}{x} = \pdiffB{2}{v}{\xi} + 2 \pdiffC{v}{\xi}{\eta} + \pdiffB{2}{v}{\eta} \tag{4.7} \] とまとめられます。xで偏微分したのに、xが式中に現れず完全に消えていることが面白いところです。

4.2.3 変数置換後の右辺

次に、式(4.3)の右辺項を変数置換し、時間tで偏微分すると、 \begin{eqnarray} \pdiffA{}{t} v(\xi,\eta) &=& \pdiffA{v}{\xi} \pdiffA{\xi}{t} + \pdiffA{v}{\eta} \pdiffA{\eta}{t} \\ &=& -u \pdiffA{v}{\xi} + u \pdiffA{v}{\eta} \tag{4.8} \end{eqnarray} ここでは、 \begin{eqnarray} \pdiffA{\xi}{t} &= \pdiffA{}{x} (x-ut) =& -u \\ \pdiffA{\eta}{t} &= \pdiffA{}{x} (x+ut) =& u \end{eqnarray} と、定数項になっていることがポイントです。

さらに(4.8)をもういちど時間tで微分することで、式(4.3)の右辺項が求められます。 \begin{eqnarray} \pdiffB{2}{}{t} v(\xi,\eta) &=& \pdiffA{ }{\xi} \left( -u \pdiffA{v}{\xi} + u \pdiffA{v}{\eta} \right) \pdiffA{\xi}{t} \\ & & \hspace{3em} + \pdiffA{}{\eta} \left( -u \pdiffA{v}{\xi} + u \pdiffA{v}{\eta} \right) \pdiffA{\eta}{t} \\ &=& u^2 \pdiffA{}{\xi} \left( \pdiffA{v}{\xi} - \pdiffA{v}{\eta} \right) \\ & & \hspace{3em} + u^2 \pdiffA{}{\eta} \left( -\pdiffA{v}{\xi} + \pdiffA{v}{\eta} \right) \\ &=& u^2 \left( \pdiffB{2}{v}{\xi} - \pdiffC{v}{\eta}{\xi} \right) \\ & & \hspace{3em} + u^2 \left( -\pdiffC{v}{\xi}{\eta} + \pdiffB{2}{v}{\eta} \right) \\ &=& u^2 \left( \pdiffB{2}{v}{\xi} - 2 \pdiffC{v}{\xi}{\eta} + \pdiffB{2}{v}{\eta} \right) \end{eqnarray} よって、 \[ \pdiffB{2}{}{t} v(\xi,\eta) = u^2 \left( \pdiffB{2}{v}{\xi} - 2 \pdiffC{v}{\xi}{\eta} + \pdiffB{2}{v}{\eta} \right) \tag{4.9} \] とまとめることができ、今度も時間tで偏微分したにもかかわらず、やはりtが式から消えています。

4.2.4 変数置換後の両辺

これで、式(4.3) \[ \pdiffB{2}{v}{x} = \frac{1}{u^2} \pdiffB{2}{v}{t} \tag{4.3} \] の変数置換された後の両辺が求められたので、式(4.7)と式(4.9)を(4.3)に代入して、 \begin{eqnarray} \pdiffB{2}{v}{\xi} + 2 \pdiffC{v}{\xi}{\eta} + \pdiffB{2}{v}{\eta} &=& \frac{1}{u^2} u^2 \left( \pdiffB{2}{v}{\xi} - 2 \pdiffC{v}{\xi}{\eta} + \pdiffB{2}{v}{\eta} \right) \\ &=& \pdiffB{2}{v}{\xi} - 2 \pdiffC{v}{\xi}{\eta} + \pdiffB{2}{v}{\eta} \tag{4.10} \end{eqnarray}

さらに両辺を整理すれば、 \[ \pdiffC{v}{\xi}{\eta} =0 \tag{4.11} \] という、重要な関係式が得られます。

4.2.5 両辺の積分による一般解

式(4.11)をηで積分します \[ \int \pdiffC{v}{\xi}{\eta} d \eta =\pdiffA{v}{\xi} = K(\xi) \tag{4.12} \] ここで、Kは積分定数です。ここでは、ξのみの関数をηで偏微分したときに0となること、すなわち \[ \pdiffA{K(\xi)}{\eta}=0 \] であることを利用して、積分定数とみなすわけです。

さらに式(4.12)をξで積分すると電圧vが出てきます \[ v(\xi,\eta) = \int K(\xi) \ d \xi +g(\eta) \tag{4.13} \] ここで、\( g(\eta) \)は積分定数であり、 \[ \pdiffA{g(\eta)}{\xi}=0 \] を利用して、積分定数とみなしています。

式(4.13)の積分項は、\( \xi \)のみの関数なので、まとめて\( f(\xi) \)と表すことができます。すると式(4.13)は、 \[ v(\xi,\eta) = f(\xi) + g(\eta) \tag{4.14} \] とまとめることできて、最終的に \[ v(x,t) = f(x-ut) + g(x+ut) \tag{4.15} \] が求めるべき方程式となります。

この解は、ダランベールの解(d'Alembert's solution)と呼ばれます。電気分野以外では、楽器の弦などの一次元方向の振動現象を表すときの方程式です。つまり、分布定数回路は弦楽器の物理振動を電気振動に置き換えたような振る舞いをする回路として考えることができるのです。

4.A 補足

4.A1 合成関数の偏微分法

2変数関数\( z=z(\xi,\eta) \)があり、\( \xi = \xi(x,t), \ \eta=\eta(x,t) \)と、それぞれの変数がx,tの合成関数になっているとします。関数zが連続で偏導関数をもつとき、zはxまたはtで偏微分することができて、 \begin{eqnarray} \pdiffA{z}{x} &=& \pdiffA{z}{\xi} \pdiffA{\xi}{x} + \pdiffA{z}{\eta} \pdiffA{\eta}{x} \\ \pdiffA{z}{t} &=& \pdiffA{z}{t} \pdiffA{\xi}{t} + \pdiffA{z}{\eta} \pdiffA{\eta}{t} \end{eqnarray} が成立します。

この公式はチェインルール(chain rule)と呼ばれていて、合成関数の常微分とは異なりますので注意が必要です。

4.A2 偏微分順序の変更

前ページでも紹介しましたが、2変数関数\( f=f(x,y) \) がC2級関数(=2階の偏導関数がすべて存在して連続の時)であれば、x->yの順番で偏微分した関数 \( f_{xy} \)と、y->xの順番で偏微分した関数\( f_{yx} \)は、一致します。式で書けば、 \[ \frac{\partial^2 f}{\partial y \partial x} = \frac{\partial^2 f}{\partial x \partial y} \]

電気回路の世界では、不連続な関数というものをあまり扱わないので、計算においては偏微分の順序変更を頻繁に行っています。