Gan's blog - フォーラム

メイン > > 電気回路理論おもちゃ箱 > 分布定数回路の過渡現象演習問題

(1) 2 3 »

投稿者	スレッド
webadm	投稿日時: 2013-5-11 21:03
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	分布定数回路の過渡現象演習問題いよいよ最後の章の演習問題に取り組むことにしよう。

webadm	投稿日時: 2013-5-12 20:14
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	波動方程式最初は波動方程式を解くことから次の２階線型偏微分放映式を波動方程式という。この方程式を解き、その結果を用いて無損失線路の電圧、電流分布の一般解を求めよ。 $\begin{eqnarray}<br />\frac{\partial^2 e}{\partial x^2}=\frac{1}{v^2}\frac{\partial^2 e}{\partial t^2}<br />\end{eqnarray}$ というもの。学校でも物理で波動方程式を解く講義を必ず受けるが、これも実は電信方程式の特殊なケースにすぎず、歴史的に早い時期に解くことができた易しい問題。最初から電信方程式に取り組むと全てのケースを解くのに２０年とかかかった人も居るぐらいだから。一番単純なケースだけやってお茶を濁すというぐらいしか学校では無理だからだ。波動方程式の一般解はd'Alembertが最初に導いたことでd'Alembertの解として知られているが、d'Alembertの方法はなぞるのと長いので、今日的にはそれ以外の色々な簡潔な導出方法があって、そちらを使ってd'Alembertの解を得るのが普通だ。ここでは著者の解とも市販の物理数学のテキストにも無いHeavisideの演算子法をつかって同じ結果を得てみることにする。いずれの方法でも時間と距離それぞれに関して常微分方程式を解くだけである。問題の偏微分方程式を以下の用に作用素方程式に書き直す $\begin{eqnarray}<br />\frac{\partial^2}{\partial x^2}&=&p^2,\ \ \ \ \frac{\partial^2}{\partial t^2}=q^2\\<br />\left(p^2-\frac{1}{v^2}q^2\right)e&=&0<br />\end{eqnarray}$ 次に時間で２重積分して初期条件を与えるポテンシャル項を出現させる $\begin{eqnarray}<br />\frac{1}{q}\left(\frac{1}{q}\left(p^2-\frac{1}{v^2}q^2\right)e+K_0\right)+K_1&=&0<br />\end{eqnarray}$ これをeに関して時間、距離に関して順番に演算子法で解くと $\begin{eqnarray}<br />e\left(x,t\right)&=&\frac{\left( {K}_{1}\,{q}^{2}+{K}_{0}\,q\right) \,{v}^{2}}{q^2-{p}^{2}\,{v}^{2}}=K_1\,v^2\frac{q^2}{q^2-p^2\,v^2}+\frac{K_0\,v^{\cancel{2}}}{p\,\cancel{v}}\frac{q\,p\,v}{q^2-p^2\,v^2}\\<br />&=&K_1\,v^2\,\cosh\,p\,v\,t+\frac{K_0\,v}{p}\sinh\,p\,v\,t\\<br />&=&K_1\,v^2\frac{e^{p\,v\,t}+e^{-p\,v\,t}}{2}+\frac{K_0\,v}{p}\frac{e^{p\,v\,t}-e^{-p\,v\,t}}{2}\\<br />&=&\frac{1}{2}\left(K_1\,v+\frac{K_0}{p}\right)v\,e^{p\,v\,t}+\frac{1}{2}\left(K_1\,v-\frac{K_0}{p}\right)v\,e^{-p\,v\,t}\\<br />&=&\frac{1}{2}\left(K_1\,v\,H\left(x-v\,t\right)+\int_{0}^{x}K_0\,H\left(z-v\,t\right)dz\right)v+\frac{1}{2}\left(K_1\,v\,H\left(x+v\,t\right)-\int_{0}^{x}K_0\,H\left(z+v\,t\right)dz\right)v\\<br />&=&e_1\left(x-v\,t\right)+e_2\left(x+v\,t\right)<br />\end{eqnarary}$ ということになる。これがd'Alembertの解である。e1,e2は初期条件によって決まる任意の距離と時間に依存する一価関数であり、互いに時間とともに速度vで空間座標軸を逆方向に移動することがわかる。さて題意では電流分布についての一般解も求めている。電流分布に関しても電信方程式はまったく同じだから、 $\begin{eqnarray}<br />i\left(x,t\right)&=&i_1\left(x-v\,t\right)+i_2\left(x+v\,t\right)<br />\end{eqnarray}$ と表される。これらの電圧分布と電流分布の一般解を電信方程式の基礎方程式に代入すると $\begin{eqnarray}<br />-\frac{\partial}{\partial x}\left(e_1\left(x-v\,t\right)+e_2\left(x+v\,t\right)\right)&=&L\frac{\partial}{\partial t}\left(i_1\left(x-v\,t\right)+i_2\left(x+v\,t\right)\right)+R\left(i_1\left(x-v\,t\right)+i_2\left(x+v\,t\right)\right)\\<br />-\frac{\partial}{\partial x}\left(i_1\left(x-v\,t\right)+i_2\left(x+v\,t\right)\right)&=&C\frac{\partial}{\partial t}\left(e_1\left(x-v\,t\right)+e_2\left(x+v\,t\right)\right)+G\left(e_1\left(x-v\,t\right)+e_2\left(x+v\,t\right)\right)<br />\end{eqnarray}$ 題意では無損失線路の場合なので、R=G=0としてよく $\begin{eqnarray}<br />\eta&=&x-v\,t,\ \ \ \ \xi=x+v\,t\\<br />-\frac{\partial}{\partial \eta}e_1\left(\eta\right)-\frac{\partial}{\partial \xi}e_2\left(\xi\right)&=&L\,v\left(-\frac{\partial}{\partial \eta}i_1\left(\eta\right)+\frac{\partial}{\partial \xi}i_2\left(\xi\right)\right)\\<br />-\frac{\partial}{\partial \eta}i_1\left(\eta\right)-\frac{\partial}{\partial \xi}i_2\left(\xi\right)&=&C\,v\left(-\frac{\partial}{\partial \eta}e_1\left(\eta\right)+\frac{\partial}{\partial \xi}e_2\left(\xi\right)\right)<br />\end{eqnarray}$ ここでvは、無損失線路での波の伝搬速度なので $\begin{eqnarray}<br />v&=&\frac{1}{\sqrt{L\,C}}<br />\end{eqnarray}$ を代入すると $\begin{eqnarray}<br />-\frac{\partial}{\partial \eta}e_1\left(\eta\right)-\frac{\partial}{\partial \xi}e_2\left(\xi\right)&=&\sqrt{\frac{L}{C}}\left(-\frac{\partial}{\partial \eta}i_1\left(\eta\right)+\frac{\partial}{\partial \xi}i_2\left(\xi\right)\right)\\<br />-\frac{\partial}{\partial \eta}i_1\left(\eta\right)-\frac{\partial}{\partial \xi}i_2\left(\xi\right)&=&\sqrt{\frac{C}{L}}\left(-\frac{\partial}{\partial \eta}e_1\left(\eta\right)+\frac{\partial}{\partial \xi}e_2\left(\xi\right)\right)<br />\end{eqnarray}$ また無損失線路の特性インピーダンスZ0は $\begin{eqnarray}<br />Z_0&=&\sqrt{\frac{L}{C}}<br />\end{eqnarray}$ であることから、 $\begin{eqnarray}<br />-\frac{\partial}{\partial \eta}e_1\left(\eta\right)-\frac{\partial}{\partial \xi}e_2\left(\xi\right)&=&Z_0\left(-\frac{\partial}{\partial \eta}i_1\left(\eta\right)+\frac{\partial}{\partial \xi}i_2\left(\xi\right)\right)\\<br />-\frac{\partial}{\partial \eta}i_1\left(\eta\right)-\frac{\partial}{\partial \xi}i_2\left(\xi\right)&=&\frac{1}{Z_0}\left(-\frac{\partial}{\partial \eta}e_1\left(\eta\right)+\frac{\partial}{\partial \xi}e_2\left(\xi\right)\right)<br />\end{eqnarray}$ という関係が成り立つ。第二式の両辺にZ0を乗じて、項を消去すると $\begin{eqnarray}<br />\frac{\partial}{\partial \eta}e_1\left(\eta\right)&=&Z_0\frac{\partial}{\partial \eta}i_1\left(\eta\right)\\<br />\frac{\partial}{\partial \xi}e_2\left(\xi\right)&=&-Z\0\frac{\partial}{\partial \xi}i_2\left(\xi\right)<br />\end{eqnarray}$ 両辺を積分して、変数変換を元に戻すと $\begin{eqnarray}<br />e_1\left(x-v\,t\right)&=&Z_0\,i_1\left(x-v\,t\right)\\<br />e_2\left(x+v\,t\right)&=&-Z_0\,i_2\left(x+v\,t\right)<br />\end{eqnarray}$ という結果が得られる。これを元の電流分布の式に代入すると $\begin{eqnarray}<br />i\left(x,t\right)&=&\frac{1}{Z_0}\left(e_1\left(x-v\,t\right)-e_2\left(x+v\,t\right)\right)<br />\end{eqnarray}$ ということになる。著者はこの関係を一言で済ませているが、それだともやもやして納得いかないという場合に几帳面に示すとこうことになる。こうやって常々基礎計算力を磨いておくとぼけ防止になるし、普段仕事上でのリスクを予見する能力を維持することにもなるので技術者として実益もあり悪いことではない。 P.S d'Alembertの解は２対あって、ひとつは上に上げたような時間を固定して空間軸上を見るのと、空間を固定して時間軸上の変化を見る場合である。後者の場合、e1(t-x/v),e2(t+x/v)のような２つの関数の和で表されることになる。いずれにせよ波動方程式は宇宙空間のように無損失な媒体を波が伝わる一番簡単なケースなので、それ以外のケースになるといきなり難しくなることが想像に難くない。

webadm	投稿日時: 2013-5-18 21:57
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	無歪み線路前問はよく考えたら無損失線路と題すべきだったかもしれない。電信方程式の解の分類としては無損失線路が一番簡単で、次に少し一般化に近く、無損失線路のケースに近いのが無歪み線路。無歪み線路の条件はHeavisideによってR/L=G/Cであることが提示された。これの有難味は、更に一般化した電信方程式の解と比べることによって明らかになるが、それはまた別の機会に。無歪み線路の電圧、電流分布の一般式を導出せよ。というもの。電信方程式の基礎方程式から始めよう $\begin{eqnarray}<br />\begin{bmatrix}\frac{\partial}{\partial x}&L\frac{\partial}{\partial t}+R\cr C\frac{\partial}{\partial t}+G& \frac{\partial}{\partial x}\end{bmatrix}\begin{bmatrix}e\left(x,t\right)\cr i\left(x,t\right)\end{bmatrix}&=&0<br />\end{eqnarray}$ これをHeavisideの演算子とベクトルで書き直すと $\begin{eqnarray}<br />\frac{\partial}{\partial x}&=&p,\ \ \ \ \ \frac{\partial}{\partial t}=q\\<br />{\bf U}&=&\begin{bmatrix}e\left(x,t\right)\cr i\left(x,t\right)\end{bmatrix}\\<br />\begin{bmatrix}p&L\,q+R\cr C\,q+G&p\end{bmatrix}{\bf U}&=&0<br />\end{eqnarray}$ 従って自明な解U＝0以外の解が存在するには $\begin{eqnarray}<br />\left\|\begin{array}p&L\,q+R\cr C\,q+G&p\end{array}\right\|&=&p^2-\left(L\,q+R\right)\left(C\,q+G\right)=p^2-L\,C\left(q+\frac{R}{L}\right)\left(q+\frac{G}{C}\right)\\<br />&=&p^2-L\,C\left(\left(q+\frac{\frac{R}{L}+\frac{G}{C}}{2}\right)^2-\left(\frac{\frac{R}{L}-\frac{G}{C}}{2}\right)^2\right)=0<br />\end{eqnarray}$ が成り立つことが必要十分条件である。ここで無損失線路の条件R/L=G/Cを適用すると $\begin{eqnarray}<br />\alpha&=&\frac{\frac{R}{L}+\frac{G}{C}}{2}\\<br />p^2-L\,C\left(q+\alpha\right)^2&=&0\\<br />\frac{1}{L\,C}p^2-\left(q+\alpha\right)^2&=&0<br />\end{eqnarray}$ ということになる。従って、解は以下の式を満足することになる $\begin{eqnarray}<br />\left(\frac{1}{L\,C}p^2-\left(q+\alpha\right)^2\right){\bf U}&=&0\\<br />\end{eqnarray}$ これを時間tで二重積分して空間に関するベクトルポテンシャル、K0,K1を出現させると $\begin{eqnarray}<br />\frac{1}{q}\left(\frac{1}{q}\left(\frac{1}{\sqrt{L\,C}}p^2-\left(q+\alpha\right)^2\right){\bf U}+{\bf K_0}\right)+{\bf K_1}&=&0<br />\end{eqnarray}$ これをUに関して解くと $\begin{eqnarray}<br />{\bf U}&=&\frac{{\bf K_1}q^2+{\bf K_0}q}{\left(q+\alpha\right)^2-\frac{1}{L\,C}p^2}=\frac{{\bf K_1}q^2}{\left(q+\alpha\right)^2-\frac{1}{L\,C}p^2}+\frac{1}{\alpha}\frac{{\bf K_0}q\alpha}{\left(q+\alpha\right)^2-\frac{1}{L\,C}p^2}\\<br />&=&e^{-\alpha\,t}{\bf K_1}\cosh\frac{1}{\sqrt{L\,C}}p\,t+\frac{e^{-\alpha\,t}}{\alpha}{\bf K_0}\sinh\frac{1}{\sqrt{L\,C}}p\,t\\<br />&=&e^{-\alpha\,t}{\bf K_1}\frac{e^{\frac{1}{\sqrt{L\,C}}p\,t}-e^{-\frac{1}{\sqrt{L\,C}}p\,t}}{2}+\frac{e^{-\alpha\,t}}{\alpha}{\bf K_0}\frac{e^{\frac{1}{\sqrt{L\,C}}p\,t}-e^{-\frac{1}{\sqrt{L\,C}}p\,t}}{2}\\<br />&=&\frac{e^{-\alpha\,t}}{2}\left({\bf K_1}+\frac{1}{\alpha}{\bf K_0}\right)e^{\frac{1}{\sqrt{L\,C}}p\,t}+\frac{e^{-\alpha\,t}}{2}\left({\bf K_1}-\frac{1}{\alpha}{\bf K_0}\right)e^{-\frac{1}{\sqrt{L\,C}}p\,t}\\<br />&=&\frac{e^{-\alpha\,t}}{2}\left({\bf K_1}+\frac{1}{\alpha}{\bf K_0}\right)H\left(x-\frac{1}{\sqrt{L\,C}}t\right)+\frac{e^{-\alpha\,t}}{2}\left({\bf K_1}-\frac{1}{\alpha}{\bf K_0}\right)H\left(x+\frac{1}{\sqrt{L\,C}}t\right)\\<br />&=&e^{-\alpha\,t}\left({\bf U_1}+{\bf U_2}\right)\\<br />{\bf U_1}&=&\begin{bmatrix}e_1\left(x-\frac{1}{\sqrt{L\,C}}t\right)\cr i_1\left(x-\frac{1}{\sqrt{L\,C}}t\right)\end{bmatrix}\\<br />{\bf U_1}&=&\begin{bmatrix}e_1\left(x+\frac{1}{\sqrt{L\,C}}t\right)\cr i_2\left(x+\frac{1}{\sqrt{L\,C}}t\right)\end{bmatrix}<br />\end{eqnarray}$ ということになる。前問と同様に電流分布は電圧の進行波と反射波の重ね合わせの関係から $\begin{eqnarray}<br />{\bf U}&=&e^{-\alpha\,t}\begin{bmatrix}e_1\left(x-\frac{1}{\sqrt{L\,C}}t\right)+e_2\left(x-\frac{1}{\sqrt{L\,C}}t\right)\cr \frac{1}{Z_0}\left(e_1\left(x-\frac{1}{\sqrt{L\,C}}t\right)-e_2\left(x+\frac{1}{\sqrt{L\,C}}t\right)\right)\end{bmatrix}\\<br />&=&e^{-\alpha\,t}\begin{bmatrix}e_1\left(x-\frac{1}{\sqrt{L\,C}}t\right)+e_2\left(x-\frac{1}{\sqrt{L\,C}}t\right)\cr \sqrt{\frac{C}{L}}\left(e_1\left(x-\frac{1}{\sqrt{L\,C}}t\right)-e_2\left(x+\frac{1}{\sqrt{L\,C}}t\right)\right)\end{bmatrix}<br />\end{eqnarray}$ ということになる。無歪み線路は無損失線路と同様に実数値の特性インピーダンスを持つことは以下に示す通り $\begin{eqnarray}<br />Z_0&=&\frac{e\left(x,t\right)}{i\left(x,t\right)}=\sqrt{\frac{L\frac{\partial}{\partial t}+R}{C\frac{\partial}{\partial t}+G}}=\sqrt{\frac{L\left(\frac{\partial}{\partial t}+\frac{R}{L}\right)}{C\left(\frac{\partial}{\partial t}+\frac{G}{C}\right)}}\\<br />&=&\sqrt{\frac{L}{C}}\ \ \ \ \left(\frac{R}{L}=\frac{G}{C}\right)<br />\end{eqnarray}$ P.S 著者は電分布の一般解を電信方程式に代入して前問のように電流分布の解を導いているが、そうしてもよい。特性インピーダンスが実数値なのは無損失線路と無歪み線路だけで、それ以外の一般の線路では複素数になる（周波数特性がフラットでない）点だけは憶えておく必要がある。これが信号が遠くまで歪みなく届くには無損失線路（宇宙空間）か無歪み線路（同軸ケーブル）のどちらかでなければならない理由である。これが深刻な技術的な問題となったのは電話の時代になってからである。それ以降このことが常識と化した。近代で身近なところではGHzのクロック周波数で動作する超LSI内の配線は必然的に半導体シリコン内で無歪み線路を形成する必要があるので、半導体ウェハー表面からどのくらい掘った（エッチングした）位置でその特性（R/L=G/C)が得られるか見極める必要がある。実際に回路をエッチングして特性を確認するわけだが、そう何度も試し掘りして試行錯誤的に見つけるわけにもいかない。事前に有る程度近い位置を解析的に求める必要がある。そこのところは半導体メーカーの超極秘のノウハウである。

webadm	投稿日時: 2013-5-19 4:19
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	無限長無損失線路これまでの問題は初期値問題だったが、今度は境界値問題無限長無損失線路の送電端にes(t)=Em e^-bt cosωtを加えるとき、任意の点xの電圧、電流はどうなるか。というもの。これも理論のときに直流電圧を加えた場合を解いてしまったが、今回は交流電圧を加えた場合となる。無損失線路(R=G=0)の基礎方程式をHeavisideの演算子とベクトルで表すと $\begin{eqnarray}<br />\frac{\partial}{\partial x}&=&p,\ \ \ \ \frac{\partial}{\partial t}=q\\<br />{\bf U}&=&\begin{bmatrix}e\left(x,t\right)\cr i\left(x,t\right)\end{bmatrix}\\<br />\begin{bmatrix}p&L\,q\cr C\,q&p\end{bmatrix}{\bf U}&=&0<br />\end{eqnarray}$ ということになる。従って特性方程式は $\begin{eqnarray}<br />\left\|\begin{array}p&L\,q\cr C\,q&p\end{array}\right\|&=&p^2-L\,C\,q^2=0<br />\left(p^2-L\,C\,q^2\right){\bf U}&=&0<br />\end{eqnarray}$ ということになる。これを距離で二重積分しベクトルポテンシャルを出現させると $\begin{eqnarray}<br />\frac{1}{p}\left(\frac{1}{p}\left(p^2-L\,C\,q^2\right){\bf U}+{\bf K_0}\right)+{\bf K_1}&=&0<br />\end{eqnarray}$ これをUに関して解くと $\begin{eqnarray}<br />{\bf U}&=&-\frac{{\bf K_1}p^2+{\bf K_0}p}{p^2-q^2\,L\,C}=-{\bf K_1}\left(\frac{p^2}{p^2-q^2\,L\,C}\right)-\frac{{\bf K_0}}{q\,\sqrt{L\,C}}\left(\frac{p\,q\,\sqrt{L\,C}}{p^2-q^2\,L\,C}\right)\\<br />&=&-{\bf K_1}\cosh\,x\,q\,\sqrt{L\,C}-\frac{{\bf K_0}}{q\,\sqrt{L\,C}}\sinh\,x\,q\,\sqrt{L\,C}<br />\end{eqnarray}$ また距離に関して解いただけの途中だが、この段階で距離に関して以下の境界条件を与えることができる。 $\begin{eqnarray}<br />\left.{\bf U}\right\|_{x=0}&=&\begin{bmatrix}E_m\,e^{-b\,t}\,\cos\omega\,t\cr \frac{1}{Z_0}E_m\,e^{-b\,t}\,\cos\omega\,t\end{bmatrix}\\<br />\left.{\bf U}\right\|_{x=x_0}&=&\begin{bmatrix}0\cr 0\end{bmatrix}<br />\end{eqnarray}$ ここでx0を送電端から十分遠い位置とする。これを先の途中解に適用するとベクトルポテンシャルに関する連立方程式が得られる $\begin{eqnarray}<br />\begin{bmatrix}E_m\,e^{-b\,t}\,\cos\omega\,t\cr \frac{1}{Z_0}E_m\,e^{-b\,t}\,\cos\omega\,t\end{bmatrix}&=&-{\bf K_1}<br />\begin{bmatrix}0\cr 0\end{bmatrix}&=&-{\bf K_1}\cosh\,x_0\,q\,\sqrt{L\,C}-\frac{{\bf K_0}}{q\,\sqrt{L\,C}}\sinh\,x_0\,q\,\sqrt{L\,C}<br />\end{eqnarray}$ 第一式を第二式に代入してK0を解くと $\begin{eqnarray}<br />{\bf K_1}&=&\begin{bmatrix}-E_m\,e^{-b\,t}\,\cos\omega\,t\cr -\frac{1}{Z_0}E_m\,e^{-b\,t}\,\cos\omega\,t\end{bmatrix}\\<br />{\bf K_0}&=&\begin{bmatrix}q\,\sqrt{L\,C}\coth\left(x_0\,q\,\sqrt{L\,C}\right)E_m\,e^{-b\,t}\cos\omega\,t\cr \frac{q\,\sqrt{L\,C}\coth\left(x_0\,q\,\sqrt{L\,C}\right)E_m\,e^{-b\,t}\cos\omega\,t}{Z_0}\frac{}{}\end{bmatrix}<br />\end{eqnarray}$ ここでx0→∞とするとK0は $\begin{eqnarray}<br />{\bf K_0}&=&\begin{bmatrix}q\,\sqrt{L\,C}E_m\,e^{-b\,t}^cos\omega\,t\cr \frac{q\,\sqrt{L\,C}E_m\,e^{-b\,t}\cos\omega\,t}{Z_0}\end{bmatrix}<br />\end{eqnarray}$ ということになる。これらベクトルポテンシャルを途中解に代入して完全に解くと $\begin{eqnarray}<br />{\bf U}&=&\begin{bmatrix}\cosh\left(x\,q\,\sqrt{L\,C}\right)E_m\,e^{-b\,t}\cos\omega\,t-\sinh\left(x\,q\,\sqrt{L\,C}\right)E_m\,e^{-b\,t}\cos\omega\,t\cr \frac{1}{Z_0}\left(\cosh\left(x\,q\,\sqrt{L\,C}\right)E_m\,e^{-b\,t}\cos\omega\,t-\sinh\left(x\,q\,\sqrt{L\,C}\right)E_m\,e^{-b\,t}\cos\omega\,t\right)\end{bmatrix}\\<br />&=&\begin{bmatrix}e^{-x\,q\,\sqrt{L\,C}}E_m\,e^{-b\,t}\cos\omega\,t\cr \frac{1}{Z_0}e^{-x\,q\,\sqrt{L\,C}}E_m\,e^{-b\,t}\cos\omega\,t\end{bmatrix}\\<br />&=&\begin{case}<br />\begin{bmatrix}0\cr 0\end{bmatrix}\ \ \ \ \left(t-x\sqrt{L\,C}\lt 0\right)\\<br />\begin{bmatrix}E_m\,e^{-b\,\left(t-x\,sqrt{L\,C}\right)}\cos\omega\left(t-x\,\sqrt{L\,C}\right)\cr \frac{1}{Z_0}E_m\,e^{-b\left(t-x\sqrt{L\,C}\right)}\cos\omega\left(t-x\,\sqrt{L\,C}\right)\end{bmatrix}\ \ \ \ \left(t-x\sqrt{L\,C}\ge 0\right)<br />\end{case}<br />\end{eqnarray}$ ということになる。著者は信号がまだ伝わっていない空間での解を書いていないが、著者の解は既に信号が伝わった以降の解である点を注意しておく。 Z0は無損失線路の特性インピーダンスで $\begin{eqnarray}<br />Z_0&=&\sqrt{\frac{L}{C}}<br />\end{eqnarray}$ である。 P.S 最後に残った時間推移演算子は時間に関する項に全て作用することに注意しなければならない。なので本来はベクトルポテンシャルK0,K1は最初の段階で演算子関数の右側に現れなければならない。順序交換すると誤った結果を導くことになる。Maximaを使うとそんなことは気にもとめてくれないので、勝手に演算順序を交換してしまうので注意。時間推移演算子が最後に残ったが、それについては共立の「数学公式」のHeavisideの演算子の章には何故か見あたらない。幸いにして手元にあるCourant & Hilbertの"Methods of Mathematical Physics"のheavisideの演算子法の解説の中に書いてある。p524に引用： 5) A very important operator which is not rational, the exponential operator, is introduces for constant h by the definition $\begin{eqnarray}<br />\left(14\right)\ \ \ \ e^{-h\,p}f\left(t\right)&=&f\left(t-h\right)<br />\end{eqnarray}$ 理論のところでこれが出てきた時に安易に"Taylorの定理により"と書いた憶えがあるが、それはHeaviside流のこじつけであって、Taylorの定理そのものは解析関数の近似級数展開で級数展開そのものが関数と同値であることを意味しているわけではない。 Courant & Hilbert本には先ほどの引用した直後に以下のように補足されている。引用： which is suggested by the Taylor series for $exp(-hp)$ and f(t-h). However, this plausibility consideration cannot provide a satisfactory justification, because the definition should not depend on the analytic character of the function f(t). Our definition is justified rather by the fact that it implies the relations $exp(-hp)exp(-kp)f(t)=exp(-(h+k)p)f(t)$ and $(d/dh)exp(-hp)f(t)=-pexp(-hp)f(t)$ the latter is equivalent to the equation $(d/dh)f(t-h)=-(d/dt)f(t-h).$

webadm	投稿日時: 2013-5-19 23:57
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	無限長無歪み線路次は境界値問題で、無限長無歪み線路に関するもの無限長無歪み線路の送電端にes(t)=Em sinωtを加えるとき、任意の点xの電圧、電流はどうなるか。というもの。基礎方程式をHeavisideの演算子とベクトルで表すと $\begin{eqnarray}<br />\frac{\partial}{\partial x}&=&p,\ \ \ \ \frac{\partial}{\partial t}=q\\<br />{\bf U}&=&\begin{bmatrix}e\left(x,t\right)\cr i\left(x,t\right)\end{bmatrix}\\<br />\begin{bmatrix}p&L\,q+R\cr C\,q+G&p\end{bmatrix}{\bf U}&=&0<br />\end{eqnarray}$ ということになる。これの特性方程式は $\begin{eqnarray}<br />\left\|\begin{array}p&L\,q+R\cr C\,q+G&p\end{array}\right\|&=&p^2-\left(L\,q+R\right)\left(C\,q+G\right)=p^2-L\,C\left(q+\frac{R}{L}\right)\left(q+\frac{G}{C}\right)\\<br />&=&p^2-L\,C\left(\left(q+\frac{1}{2}\left(\frac{R}{L}+\frac{G}{C}\right)\right)^2-\left(\frac{1}{2}\left(\frac{R}{L}-\frac{G}{C}\right)\right)^2\right)\\<br />&=&p^2-L\,C\left(\left(q+\alpha\right)^2-\beta^2\right)\\<br />&=&0\\<br />\alpha&=&\frac{1}{2}\left(\frac{R}{L}+\frac{G}{C}\right),\ \ \ \ \beta=\frac{1}{2}\left(\frac{R}{L}-\frac{G}{C}\right)<br />\end{eqnarray}$ 従って無歪み線路(R/L=G/C)の場合 $\begin{eqnarray}<br />\left(p^2-L\,C\left(q+\alpha\right)^2\right){\bf U}&=&0<br />\end{eqnarray}$ ということになる。前問と同様に距離で二重積分してベクトルポテンシャルを出現させると $\begin{eqnarray}<br />\frac{1}{p}\left(\frac{1}{p}\left(p^2-L\,C\left(q+\alpha\right)^2\right){\bf U}+{\bf K_0}\right)+{\bf K_1}&=&0<br />\end{eqnarray}$ これをUに関して解くと $\begin{eqnarray}<br />{\bf U}&=&-\frac{{\bf K_1}p^2+{\bf K_0}p}{p^2-L\,C\left(q+\alpha\right)^2}=-\left(\frac{p^2}{p^2-L\,C\left(q+\alpha\right)^2}\right){\bf K_1}-\left(\frac{p\,\sqrt{L\,C}\left(q+\alpha\right)}{p^2-L\,C\left(q+\alpha\right)^2}\right)\frac{{\bf K_0}}{\sqrt{L\,C}\left(q+\alpha\right)}\\<br />&=&-{\bf K_1}\cosh\left(x\,\sqrt{L\,C}\left(q+\alpha\right)\right)-\frac{{\bf K_0}}{\sqrt{L\,C}\left(q+\alpha\right)}\sinh\left(x\sqrt{L\,C}\left(q+\alpha\right)\right)<br />\end{eqnarray}$ ここで以下の無限長線路の境界条件を与える $\begin{eqnarray}<br />\left.{\bf U}\right\|_{x=0}&=&\begin{bmatrix}E_m\,\sin\omega\,t\cr \frac{1}{Z_0}E_m\,\sin\omega\,t\end{bmatrix}\\<br />\left.{\bf U}\right\|_{x=x_0}&=&\begin{bmatrix}0\cr 0\end{bmatrix}<br />\end{eqnarray}$ ここでx0は送電端より十分遠い距離とする。これを途中の解に適用すると $\begin{eqnarray}<br />\begin{bmatrix}E_m\,\sin\omega\,t\cr \frac{1}{Z_0}E_m\,\sin\omega\,t\end{bmatrix}&=&-{\bf K_1}\\<br />\begin{bmatrix}0\cr 0\end{bmatrix}&=&-{\bf K_1}\cosh\left(x_0\,\sqrt{L\,C}\left(q+\alpha\right)\right)-\frac{{\bf K_0}}{\sqrt{L\,C}\left(q+\alpha\right)}\sinh\left(x_0\sqrt{L\,C}\left(q+\alpha\right)\right)\\<br />\end{eqnarray}$ 第一式を第二式に代入してx0→∞とすると $\begin{eqnarray}<br />{\bf K_0}&=&\begin{bmatrix}\sqrt{L\,C}\left(q+\alpha\right)\coth\left(x_0\sqrt{L\,C}\left(q+\alpha\right)\right)E_m\,\sin\omega\,t\cr \frac{1}{Z_0}\sqrt{L\,C}\left(q+\alpha\right)\coth\left(x_0\sqrt{L\,C}\left(q+\alpha\right)\right)E_m\,\sin\omega\,t\end{bmatrix}\\<br />&=&\begin{bmatrix}\sqrt{L\,C}\left(q+\alpha\right)E_m\,\sin\omega\,t\cr \frac{1}{Z_0}\sqrt{L\,C}\left(q+\alpha\right)E_m\,\sin\omega\,t\end{bmatrix}\ \ \ \left(x_0\to\infty\right)\\<br />{\bf K_1}&=&\begin{bmatrix}-E_m\,\sin\omega\,t\cr -\frac{1}{Z_0}E_m\,\sin\omega\,t\end{bmatrix}<br />\end{eqnarray}$ ということになる。これを途中解に代入して完全に解くと $\begin{eqnarray}<br />{\bf U}&=&\begin{bmatrix}e^{-x\,\sqrt{L\,C}\alpha}e^{-x\,\sqrt{L\,C}q}E_m\,\sin\omega\,t\cr \frac{1}{Z_0}e^{-x\,\sqrt{L\,C}\alpha}e^{-x\,\sqrt{L\,C}q}E_m\,\sin\omega\,t\end{bmatrix}\\<br />&=&\begin{case}<br />\begin{bmatrix}0\cr 0\end{bmatrix}\ \ \ \ \left(t\lt x\sqrt{L\,C}\right)\\<br />\begin{bmatrix}e^{-x\,\sqrt{L\,C}\alpha}E_m\,\sin\omega\left(t-x\sqrt{L\,C}\right)\cr \frac{1}{Z_0}e^{-x\,\sqrt{L\,C}\alpha}E_m\,\sin\omega\left(t-x\sqrt{L\,C}\right)\end{bmatrix}\ \ \ \ \left(t\ge x\sqrt{L\,C}\right)<br />\end{case}<br />\end{eqnarray}$ ということになる。これでも十分だが更に整理すると $\begin{eqnarray}<br />\sqrt{L\,C}\alpha&=&\sqrt{L\,C}\frac{1}{2}\left(\frac{R}{L}+\frac{G}{C}\right)=\frac{1}{2}\left(\sqrt{\frac{C}{L}R^2}+\sqrt{\frac{L}{C}G^2}\right)=\frac{1}{2}\left(\sqrt{\cancel{C}\,R\frac{G}{\cancel{C}}}+\sqrt{\cancel{L}\,G\frac{R}{\cancel{L}}}\right)\ \ \ \left(\frac{R}{L}=\frac{G}{C}\right)\\<br />&=&\sqrt{R\,G}<br />\end{eqnarray}$ であることから $\begin{eqnarray}<br />{\bf U}&=&\begin{case}<br />\begin{bmatrix}0\cr 0\end{bmatrix}\ \ \ \ \left(t\lt x\sqrt{L\,C}\right)\\<br />\begin{bmatrix}e^{-x\,\sqrt{R\,G}}E_m\,\sin\omega\left(t-x\sqrt{L\,C}\right)\cr \frac{1}{Z_0}e^{-x\,\sqrt{R\,G}}E_m\,\sin\omega\left(t-x\sqrt{L\,C}\right)\end{bmatrix}\ \ \ \ \left(t\ge x\sqrt{L\,C}\right)<br />\end{case}<br />\end{eqnarray}$ としても良い。ここで特性インピーダンスZ0は無歪み線路では実数値で $\begin{eqnarray}<br />Z_0&=&\sqrt{\frac{L}{C}}<br />\end{eqnarray}$ である。

webadm	投稿日時: 2013-5-20 6:55
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	無限長RC線路これも理論の時にやってしまったが、境界値問題である t=0で直流電圧Eを加えるとき無限長線路の電圧、電流分布はどうなるか。またE=1[V],R=3[Ω/km],C=0.2[uF/km]のときの電圧、電流を図示せよ。というもの。まずは電圧、電流の解を求める必要がある。 RC線路の基礎方程式をHeavisideの演算子とベクトルで表すと $\begin{eqnarray}<br />\frac{\partial}{\partial x}&=&p,\ \ \ \ \frac{\partial}{\partial t}=q\\<br />{\bf U}&=&\begin{bmatrix}e\left(x,t\right)\cr i\left(x,t\right)\end{bmatrix}\\<br />\begin{bmatrix}p&R\cr C\,q&p\end{bmatrix}{\bf U}&=&0<br />\end{eqnarray}$ ということになる。これの特性方程式は $\begin{eqnarray}<br />\left\|\begin{array}p&R\cr C\,q&p\end{array}\right\|&=&p^2-R\,C\,q=0\\<br />\left(p^2-R\,C\,q\right){\bf U}&=&0<br />\end{eqnarray}$ ということになる。これは熱方程式である。距離で２重積分してベクトルポテンシャル項を出現させると $\begin{eqnarray}<br />\frac{1}{p}\left(\frac{1}{p}\left(p^2-R\,C\,q\right){\bf U}+{\bf K_0}\right)+{\bf K_1}&=&0<br />\end{eqnaray}$ これをUについて解くと $\begin{eqnarray}<br />{\bf U}&=&-\frac{{\bf K_1}p^2+{\bf K_0}p}{p^2-R\,C\,q}=-\left(\frac{p^2}{p^2-R\,C\,q}\right){\bf K_1}-\left(\frac{p\,\sqrt{R\,C\,q}}{p^2-R\,C\,q}\right)\frac{{\bf K_0}}{\sqrt{R\,C\,q}}\\<br />&=&-{\bf K_1}\cosh\,x\sqrt{R\,C\,q}-\frac{{\bf K_0}}{\sqrt{R\,C\,q}}\sinh\,x\sqrt{R\,C\,q}<br />\end{eqnarray}$ ということになる。ここで以下の無限長線路の境界条件を与えると $\begin{eqnarray}<br />\left.{\bf U}\right\|_{x=0}&=&-{\bf K_1}=\begin{bmatrix}E\cr \frac{E}{Z_0}\end{bmatrix}\\<br />\left.{\bf U}\right\|_{x=x_0}&=&-{\bf K_1}\cosh\,x_0\,\sqrt{R\,C\,q}-\frac{{\bf K_0}}{\sqrt{R\,C\,q}}\sinh\,x_0\,\sqrt{R\,C\,q}=\begin{bmatrix}0\cr 0\end{bmatrix}<br />\end{eqnarray}$ 第一式を第二式に代入してK0を解くと $\begin{eqnarray}<br />{\bf K_0}&=&\begin{bmatrix}\sqrt{R\,C\,q}E\,\coth\,x_0\,\sqrt{R\,C\,q}\cr \frac{1}{Z_0}\sqrt{R\,C\,q}E\,\coth\,x_0\,\sqrt{R\,C\,q}\end{bmatrix}\\<br />&=&\begin{bmatrix}\sqrt{R\,C\,q}E\cr \frac{E}{Z_0}\sqrt{R\,C\,q}\end{bmatrix}\ \ \ \ \left(x_0\to \infty\right)<br />{\bf K_1}&=&\begin{bmatrix}-E\cr -\frac{E}{Z_0}\end{bmatrix}<br />\end{eqnarray}$ これを途中解に代入して完全に解くと $\begin{eqnarray}<br />{\bf U}&=&\begin{bmatrix}E\,\left(\cosh\,x\,\sqrt{R\,C\,q}-\sinh\,x\,\sqrt{R\,C\,q}\right)\cr \frac{E}{Z_0}\left(\cosh\,x\,\sqrt{R\,C\,q}-\sinh\,x\,\sqrt{R\,C\,q}\right)\end{bmatrix}\\<br />&=&\begin{bmatrix}E\left(\frac{\cancel{e^{x\sqrt{R\,C\,q}}+e^{-x\sqrt{R\,C\,q}}}}{2}-E\frac{\cancel{e^{x\sqrt{R\,C\,q}}}-e^{-x\sqrt{R\,C\,q}}}{2}\right)\cr E\left(\sqrt{\frac{C\,q}{R}}\frac{\cancel{e^{x\sqrt{R\,C\,q}}}+e^{-x\sqrt{R\,C\,q}}}{2}-\sqrt{\frac{C\,q}{R}}\frac{\cancel{e^{x\sqrt{R\,C\,q}}}-e^{-x\sqrt{R\,C\,q}}}{2}\right)\end{bmatrix}\\<br />&=&\begin{bmatrix}Ee^{-x\,\sqrt{R\,C\,q}}\cr E\sqrt{\frac{C\,q}{R}}e^{-x\sqrt{R\,C\,q}}\end{bmatrix}\\<br />&=&\begin{bmatrix}E\left(1-erf\left(\frac{x\sqrt{R\,C}}{2\,\sqrt{t}}\right)\right)\cr E\sqrt{\frac{C}{R\,\pi\,t}}e^{-\frac{x^2\,R\,C}{4\,t}}\end{bmatrix}\\<br />Z_0&=\sqrt{\frac{R}{C\,q}}<br />\end{eqnarray}$ ということになる。 RC線路の場合には線路が抵抗Rで直列に接続されているので、信号はすぐさま抵抗分圧でどんなに遠い受電端でも即時に伝わるという矛盾があるため、波が伝わらないという時空間が存在しない。著者は定性的に電圧と電流の式を調べているが、面倒なのでMaximaでプロットしてしまおう。電圧分布を時間と距離に関してプロットすると電流分布を時間と距離に関してプロットするとということになる。著者はついでに電流のピーク値を求める式を導出している。

webadm	投稿日時: 2013-5-20 9:24
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	損失のある無限長線路次も境界値問題で、更に一般化した線路損失のあるケース G=0の無限長線路にt=0で直流電圧Eをくわえたときの電圧、電流分布を求めよ。というもの。 R＞０、L＞0、C＞0、G=0の線路の基礎方程式をHeavisideの演算子とベクトルで表すと $\begin{eqnarray}<br />\frac{\partial}{\partial x}&=&p,\ \ \ \ \frac{\partial}{\partial t}=q\\<br />{\bf U}&=&\begin{bmatrix}e\left(x,t\right)\cr i\left(x,t\right)\end{bmatrix}\\<br />\begin{bmatrix}p&L\,q+R\cr C\,q&p\end{bmatrix}{\bf U}&=&0<br />\end{eqnarray}$ ということになる。これの特性方程式は $\begin{eqnarray}<br />\left\|\begin{array}p&L\,q+R\cr C\,q&p\end{array}\right\|&=&p^2-C\,q\left(L\,q+R\right)=p^2-L\,C\,q\left(q+\frac{R}{L}\right)=0\\<br />\left(p^2-L\,C\,q\left(q+\frac{R}{L}\right)\right){\bf U}&=&0<br />\end{eqnarray}$ 距離で二重積分してベクトルポテンシャルを出現させると $\begin{eqnarray}<br />\frac{1}{p}\left(\frac{1}{p}\left(p^2-L\,C\,q\left(q+\frac{R}{L}\right)\right){\bf U}+{\bf K_0}\right)+{\bf K_1}&=&0<br />\end{eqnarray}$ これをUに関して解くと $\begin{eqnarray}<br />{\bf U}&=&-\frac{{\bf K_1}p^2+{\bf K_0}p}{p^2-C\,L\,q\left(q+\frac{R}{L}\right)}=-\left(\frac{p^2}{p^2-C\,L\,q\left(q+\frac{R}{L}\right)}\right){\bf K_1}-\left(\frac{p\,\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}}{p^2-C\,L\,q\left(q+\frac{R}{L}\right)}\right)\frac{{\bf K_0}}{\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}}\\<br />&=&-{\bf K_1}\cosh\,x\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}-\frac{{\bf K_0}}{\sqrt{C\,Lq\left(q+\frac{R}{L}\right)}}\sinh\,x\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}<br />\end{eqnarray}$ ということになる。ここで以下の無限長線路の境界条件を与えると $\begin{eqnarray}<br />\left.{\bf U}\right\|_{x=0}&=&-{\bf K_1}=\begin{bmatrix}E\cr \frac{E}{Z_0}\end{bmatrix}\\<br />\left.{\bf U}\right\|_{x=x_0}&=&-{\bf K_1}\cosh\,x_0\,\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}-\frac{{\bf K_0}}{\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}}\sinh\,x_0\,\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}=\begin{bmatrix}0\cr 0\end{bmatrix}<br />\end{eqnarray}$ 第一式を第二式に代入してK0を解くと $\begin{eqnarray}<br />{\bf K_0}&=&\begin{bmatrix}E\,\sqrt{C\,L\,q\,\left(q+\frac{R}{L}\right)}\coth\,x_0\,\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}\cr \frac{E}{Z_0}\coth\,x_0\,\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}\end{bmatrix}\\<br />&=&\begin{bmatrix}E\,\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}\cr \frac{E}{Z_0}\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}\end{bmatrix}\ \ \ \ \left(x_0\to \infty\right)<br />{\bf K_1}&=&\begin{bmatrix}-E\cr -\frac{E}{Z_0}\end{bmatrix}<br />\end{eqnarray}$ ということになる。これを途中解に代入して完全に解くと $\begin{eqnarray}<br />{\bf U}&=&\begin{bmatrix}E\left(\cosh\,x\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}-\sinh\,x\,\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}\right)\cr \frac{E}{Z_0}\left(\cosh\,x\,\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}-\sinh\,x\,\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}\right)\end{bmatrix}\\<br />&=&\begin{bmatrix}E^{-x\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}}\cr E\sqrt{\frac{C}{L}}\sqrt{\frac{q}{q+\frac{R}{L}}}e^{-x\sqrt{C\,L\,q\left(q+\frac{R}{L}\right)}}\end{bmatrix}\\<br />&=&\begin{bmatrix}E\,e^{-x\sqrt{C\,L}\sqrt{\left(q+\frac{R}{2\,L}\right)^2-\left(\frac{R}{2\,L}\right)^2}}\cr E\sqrt{\frac{C}{L}}\frac{q}{\sqrt{\left(q+\frac{R}{2\,L}\right)^2-\left(\frac{R}{2\,L}\right)^2}}e^{-x\sqrt{C\,L}\sqrt{\left(q+\frac{R}{2\,L}\right)^2-\left(\frac{R}{2\,L}\right)^2}}\end{bmatrix}\\<br />&=&\begin{bmatrix}E\,e^{-\frac{R}{2\,L}t}e^{-x\sqrt{C\,L}\sqrt{q^2-\left(\frac{R}{2\,L}\right)^2}}\cr E\sqrt{\frac{C}{L}}e^{-\frac{R}{2\,L}t}\frac{q}{\sqrt{q^2-\left(\frac{R}{2\,L}\right)^2}}e^{-x\sqrt{C\,L}\sqrt{q^2-\left(\frac{R}{2\,L}\right)^2}}\end{bmatrix}\\<br />&=&\begin{case}<br />\begin{bmatrix}0\cr 0\end{bmatrix}\ \ \ \ \left(t-x\sqrt{C\,L}\lt 0\right)\\<br />\begin{bmatrix}E\left(e^{-\frac{R}{2\,L\,\sqrt{C\,L}}x}+\frac{R}{2\,L\,\sqrt{C\,L}}x\int_{x\,\sqrt{C\,L}}^{t}\frac{e^{-\frac{R}{2\,L}\tau}I_1\left(\frac{R}{2\,L}\sqrt{\tau^2-\left(x\sqrt{C\,L}\right)^2}\right)}{\sqrt{\tau^2-\left(x\sqrt{C\,L}\right)^2}}d\tau\right)\cr E\sqrt{\frac{C}{L}}e^{-\frac{R}{2\,L}t}I_0\left(\frac{R}{2\,L}\sqrt{t^2-\left(x\sqrt{C\,L}\right)^2}\right)\end{bmatrix}\ \ \ \ \left(t-x\sqrt{C\,L}\ge 0\right)\\<br />\end{case}<br />Z_0&=&\sqrt{\frac{L\,q+R}{C\,q}}=\sqrt{\frac{L}{C}}\sqrt{\frac{q+\frac{R}{L}}{q}}<br />\end{eqnarray}$ ということになる。 I0,I1はそれぞれ変形Bessel関数である $\begin{eqnarray}<br />I_n\left(z\right)&=&i^{n}J_n\left(i\,z\right)<br />\end{eqnarray}$ ちなみに著者の最後の電圧分布の解は、距離xではなく時間tで微分するように記述されているなど間違っていることに注意。例え距離で微分したとしても得られる解は最初に与えた境界条件を満足しない。それに本来なら電流分布の解を距離xで微分しなければならないが、それは(1/s)を全体に適用した後に行わなければならない。著者はもしかして演算子法を知らないか、知っていてもMiksinskyや吉田の著書を読んだことがないのだろうと想像する。さすがにこの解は電圧、電流分布の計算は要求していない。数値計算が難しいので著者も正解を出せないからだろう。これは読者の課題としよう( ´∀`)

webadm	投稿日時: 2013-5-21 7:38
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	続：無損失線路次は再び無損失の問題無損失線路の電圧、電流をつぎの初期条件のもとに解け。 $e\left(0,t\right)=e_S\left(t\right)\ \ \ \ i\left(0,t\right)=i_S\left(t\right)$ というもの。以前の問題の結果から解は以下の特性方程式を満足する。 $\begin{eqnarray}<br />\frac{\partial}{\partial x}&=&p,\ \ \ \ \frac{\partial}{\partial t}=q\\<br />{\bf U}&=&\begin{bmatrix}e\left(x,t\right)\cr i\left(x,t\right)\end{bmatrix}\\<br />\left(p^2-L\,C\,q^2\right){\bf U}&=&0<br />\end{eqnarray}$ これを例によって距離で二重積分してベクトルポテンシャルを出現させると $\begin{eqnarray}<br />\frac{1}{p}\left(\frac{1}{p}\left(p^2-L\,C\,q^2\right){\bf U}+{\bf K_0}\right)+{\bf K_1}&=&0<br />\end{eqnarray}$ これをUについて解くと $\begin{eqnarray}<br />{\bf U}&=&-\frac{{\bf K_1}p^2+{\bf K_0}p}{p^2-L\,C\,q^2}=-\left(\frac{p^2}{p^2-L\,C\,q^2}\right){\bf K_1}-\left(\frac{p\,\sqrt{L\,C}q}{p^2-L\,C\,q^2}\right)\frac{{\bf K_0}}{\sqrt{L\,C}q}\\<br />&=&-{\bf K_1}\cosh\,x\sqrt{L\,C}q-\frac{{\bf K_0}}{\sqrt{L\,C}q}\sinh\,x\sqrt{L\,C}q\\<br />&=&-{\bf K_1}\frac{e^{x\sqrt{L\,C}q}+e^{-x\sqrt{L\,C}q}}{2}-\frac{{\bf K_0}}{\sqrt{L\,C}q}\frac{e^{x\sqrt{L\,C}q}-e^{-x\sqrt{L\,C}q}}{2}\\<br />&=&-\frac{1}{2}\left({\bf K_1}+\frac{{\bf K_0}}{\sqrt{L\,C}q}\right)e^{x\sqrt{L\,C}q}-\frac{1}{2}\left({\bf K_1}-\frac{{\bf K_0}}{\sqrt{L\,C}q}\right)e^{-x\sqrt{L\,C}q}<br />\end{eqnarray}$ 従って題意の初期条件を与えると $\begin{eqnarray}<br />\left.{\bf U}\right\|_{x=0}&=&-{\bf K_1}=\begin{bmatrix}e_S\left(t\right)\cr i_S\left(t\right)\end{bmatrix}<br />\end{eqnarray}$ これでベクトルポテンシャルK1は解けたが、K0がまだ残っている。どうすんだこっから（´Д｀；）とりあえず代入して整理してみると $\begin{eqnarray}<br />{\bf U}&=&\frac{1}{2}\begin{bmatrix}e_S\left(t+x\sqrt{L\,C}\right)+e_S\left(t-x\sqrt{L\,C}\right)\cr i_S\left(t+x\sqrt{L\,C}\right)+i_S\left(t-x\sqrt{L\,C}\right)\end{bmatrix}-\frac{1}{2\sqrt{L\,C}q}\left(e^{x\sqrt{L\,C}q}-e^{-x\sqrt{L\,C}q}\right){\bf K_0}<br />\end{eqnarray}$ ということになる。さていやなことに未定積分項K0が残ってしまった。 K0を決定するにはx≠0での境界条件が必要になる。ここでどうして著者は未定項を含まない解を得ているのか謎である。境界条件がx=0に関してだけならどうしても完全には解けないはずである。著者の解を見ると、予めe(x,0)=i(x,0)=e'(x,0)=i'(x,0)=0という条件で導いたLaplace変換の結果を用いている。つまり題意以外の初期条件があるのである。これは明らかに詐欺だろう。この時点でもうまじめに解く気が失せてしまった。著者の隠していた初期条件を用いて著者と同じ結果が得られるかどうか確認するのは読者の課題としよう( ´∀`) 問題文でのみ与えられたx=0の初期条件を満足するならK0は任意に決めてよいことになる。幸いにしてK0=0としても、題意の初期条件を満足するので、それを解とすることにする。K0≠0でx=0の初期条件を満足する解は無数にあると思われる。それを全部調べるわけにはいかない。 es(t)=sin(t),L=C=1としてプロットするとということになる。これはwave equationで画像検索すると似たようなグラフが出てくるが、時空境界線より先には波はまだ届いてない領域があるので、そこは平坦ではければならないが、そうなってないものがほとんどである。どんぶらこっこどんぶらこっこと1/2振幅の２つの波が重なって進んでいく様子がよくわかる。学生時代はコンピューターとか教室になかったから、黒板の上で数式から波を想像するしかなかったが、よほどの想像力豊かか数学の天才でないかぎり、このグラフを思い浮かべるのは無理だろう。やさしい時代になった。 P.S 最初に送電端(x=0)での電圧と電流は一次独立ではないと思ったが、良く考えると送電端に反射波が到達している場合には、電圧と電流は一次独立となる。これは送電端の終端条件にもよるが、すくなくとも反射波の電圧が加わった時点で終端回路に電流が流れ、それは送電端に加わる電圧によって流れる電流を打ち消したり増大させることになる。そうすると必ずしも送電端の電圧と電流が一次独立ではないとは言えない。もっと良く考えれば、無限に長い線路の途中の任意の点を送電端としても問題としては成り立つ。x=0と定めた線路上からxが増大する方向は線路定数が均一だとすればよい。逆にx<0の空間に線路があってもよいし、実際の電源が印可されているのはずっと手前でもよいし、x=0より手前の線路の特性はどうなっていても構わないということになる。問題は長い線路上のx=0とした点からx>0の均一な特性の線路上の電圧と電流である。演算子法で導出した途中解を更に無理矢理完全に特と、第二項は積分項になることが想像に難くない。1/pの演算子関数があるので時間に関する積分になるのである。同時にexp(±x√LCq)演算子関数があるので、時間シフトが作用することも予想される。これらを総合すると、解は一次元のStokesの定理（特殊な場合のGreenの定理やGaussの定理を一般化したもの）に対応することが予想される。StokesというとWillam Thomsonの終生の友人で数学者だが、その業績を学ぶ機会は普通の人はまずない。でもHeavisideの演算子法を操ると意識するしないに関わらずその結果が見え隠れするのは実に興味深い。

webadm	投稿日時: 2013-5-24 15:12
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	続：無歪み線路次は前問の無歪み線路バージョン無歪み線路の電圧、電流を次ぎの初期条件のもとに解け。 $\begin{eqnarray}<br />e\left(0,t\right)&=&e_S\left(t\right),\ \ \ \ i\left(0,t\right)=i_S\left(t\right)<br />\end{eqnarray}$ まじめに解こうかどうしようか悩むところ。無歪み線路の基礎方程式をHeavisideの演算子とベクトルで表すと $\begin{eqnarray}<br />\frac{\partial}{\partial x}&=&p,\ \ \ \ \frac{\partial}{\partial t}=q\\<br />{\bf U}&=&\begin{bmatrix}e\left(x,t\right)\cr i\left(x,t\right)\end{bmatrix}\\<br />\begin{bmatrix}p&L\,q+R\cr C\,q+G&p\end{bmatrix}{\bf U}&=&0<br />\end{eqnarray}$ これの特性方程式は $\begin{eqnarray}<br />\left\|\begin{array}p&L\,q+R\cr C\,q+G&p\end{array}\right\|&=&p^2-\left(L\,q+R\right)\left(C\,q+G\right)=p^2-L\,C\left(q+\frac{R}{L}\right)\left(q+\frac{G}{C}\right)\\<br />&=&p^2-L\,C\left((q+\frac{1}{2}\left(\frac{R}{L}+\frac{G}{C}\right)^2-\left(\frac{\frac{R}{L}-\frac{G}{C}}{2}\right)^2right)\\<br />&=&p^2-L\,C\left(\left(q+\alpha\right)^2-\beta^2\right)\\<br />&=&p^2-L\,C\left(q+\alpha\right)^2\ \ \ \ \left(\beta=0\right)<br />&=&0\\<br />\left(p^2-L\,C\left(q+\alpha\right)^2\right){\bf U}&=&0\\<br />\alpha&=&\frac{1}{2}\left(\frac{R}{L}+\frac{G}{C}\right),\ \ \ \ \beta=\frac{1}{2}\left(\frac{R}{L}-\frac{G}{C}\right)<br />\end{eqnarray}$ 例によって距離xで二重積分してベクトルポテンシャルを出現させると $\begin{eqnarray}<br />\frac{1}{p}\left(\frac{1}{p}\left(p^2-L\,C\left(q+\alpha\right)^2\right){\bf U}+{\bf K_1}\right)+{\bf K_0}&=&0<br />\end{eqnarray}$ これをUについて解くと $\begin{eqnarray}<br />{\bf U}&=&-\frac{{\bf K_1}p^2+{\bf K_0}p}{p^2-L\,C\,\left(q+\alpha\right)^2}=-\left(\frac{p^2}{p^2-L\,C\left(q+\alpha\right)^2}\right){\bf K_1}-\left(\frac{p\,\sqrt{L\,C}\left(q+\alpha\right)}{p^2-L\,C\left(q+\alpha\right)^2}\right)\frac{{\bf K_0}}{\sqrt{L\,C}\left(q+\alpha\right)}\\<br />&=&-{\bf K_1}\cosh\,x\sqrt{L\,C}\left(q+\alpha\right)-\frac{{\bf K_0}}{\sqrt{L\,C}\left(q+\alpha\right)}\sinh\,x\sqrt{L\,C}\left(q+\alpha\right)\\<br />\end{eqnarray}$ ということになる。これにx=0での初期条件を適用すると $\begin{eqnarray}<br />\left.{\bf U}\right\|_{x=0}&=&-{\bf K_1}=\begin{bmatrix}e_S\left(t\right)\cr i_S\left(t\right)\end{bmatrix}\\<br />{\bf K_1}&=&-\begin{bmatrix}e_S\left(t\right)\cr i_S\left(t\right)\end{bmatrix}<br />\end{eqnarray}$ とK1が決まる。これを元の途中解に代入すると $\begin{eqnarray}<br />{\bf U}&=&\begin{bmatrix}e_S\left(t\right)\cr i_S\left(t\right)\end{bmatrix}\cosh\,x\sqrt{L\,C}\left(q+\alpha\right)-\frac{{\bf K_0}}{\sqrt{L\,C}\left(q+\alpha\right)}\sinh\,x\sqrt{L\,C}\left(q+\alpha\right)\\<br />&=&\begin{bmatrix}e_S\left(t\right)\cr i_S\left(t\right)\end{bmatrix}\frac{e^{x\sqrt{L\,C}\left(q+\alpha\right)}+e^{-x\sqrt{L\,C}\left(q+\alpha\right)}}{2}-\frac{{\bf K_0}}{\sqrt{L\,C}\left(q+\alpha\right)}\sinh\,x\sqrt{L\,C}\left(q+\alpha\right)\\<br />&=&\begin{bmatrix}e^{x\sqrt{L\,C}\alpha}e_S\left(t+x\sqrt{L\,C}\right)+e^{-x\sqrt{L\,C}\alpha}e_S\left(t-x\sqrt{L\,C}\right)\cr e^{x\sqrt{L\,C}\alpha}i_S\left(t+x\sqrt{L\,C}\right)+e^{-x\sqrt{L\,C}\alpha}i_S\left(t-x\sqrt{L\,C}\right)\end{bmatrix}-\frac{{\bf K_0}}{\sqrt{L\,C}\left(q+\alpha\right)}\sinh\,x\sqrt{L\,C}\left(q+\alpha\right)\\<br />&=&\frac{1}{2}\begin{bmatrix}e^{x\sqrt{R\,G}}e_S\left(t+x\sqrt{L\,C}\right)+e^{-x\sqrt{R\,G}}e_S\left(t-x\sqrt{L\,C}\right)\cr e^{x\sqrt{R\,G}}i_S\left(t+x\sqrt{L\,C}\right)+e^{-x\sqrt{R\,G}}i_S\left(t-x\sqrt{L\,C}\right)\end{bmatrix}-\frac{{\bf K_0}}{\sqrt{L\,C}\left(q+\alpha\right)}\sinh\,x\sqrt{L\,C}\left(q+\alpha\right)<br />\end{eqnarray}$ ということになる。境界条件がもうひとつ与えられないとK0が決定できないが。これで解けないと投げ出すのは早い。K0=0の場合の解は確かに与えられた初期条件を満足する。しかしK0≠0の解が存在する可能性が依然として残る。これはいわばDiophantus方程式を実数領域に拡張したような問題で、与えられたx=0の初期条件を満たすK0がひつと以上存在することは明らかである。K0≠0の解が有限個なのかそれとも無数に存在するかは現時点ではわからない。またその解を全てプロットすることもできない。著者はいかにも問題文から完全な解を導いてみせているが、それは予めe(x,0)=i(x,0)=e'(x,0)=i'(x,0)=0という初期条件を与えて導いたLaplace変換の結果を用いているためである。だったらその条件を問題文に含めないと詐欺だろう。暗黙の条件としてそれが常識だというならこちらの完敗というしかない。著者の暗黙の初期条件を使って完全な解を導くのは読者の課題としよう( ´∀`) Heavisideの演算子法でさえ、無数の他の解の存在を示している。Laplace変換ではそれがまったく示すことができないのだと思ってしまう。無損失線路の場合は送電端での源信号がそのまま空間を伝わるのに対して無歪み線路の場合は、進行波は遠くに行くにつれ減衰し、逆に反射波は送電端に近づくほど減衰し、遠ざかるほど増大する。無限長線路の場合受電端の条件を与えない限り、送電端から遠ざかるほど無限に増大する源信号と同じ形の反射波が解の中に含まれる。 K0=0とした解が題意で与えられた初期条件を満たすことから、解であることは明らかである。従って解は $\begin{eqnarray}<br />{\bf U}&=&\frac{1}{2}\begin{bmatrix}e^{x\sqrt{R\,G}}e_S\left(t+x\sqrt{L\,C}\right)+e^{-x\sqrt{R\,G}}e_S\left(t-x\sqrt{L\,C}\right)\cr e^{x\sqrt{R\,G}}i_S\left(t+x\sqrt{L\,C}\right)+e^{-x\sqrt{R\,G}}i_S\left(t-x\sqrt{L\,C}\right)\end{bmatrix}<br />\end{eqnarray}$ ちなみに著者の解には誤植があるので注意しておく。Laplace変換だと無限長線路の解しか得られないことになる（送電端から距離が遠ざかるほど信号が減衰していく）。Laplace変換は積分変換なので、右半平面上での積分が収束しないと扱えない。なので無限遠点に近づくほど値が小さくなる解しか扱えないわけである。解がたったそれだけではないことは明らかである。 eS(t)=sin(t),L=C=G=B=1としてプロットすると（ただし縦軸範囲を-1～1に制限）ということになる。これでわかるのはx=0のプロットが題意で与えられた通り源信号eS(t)そのものであることからこれが解であることは確かである。しかし、無限に距離と共に増大する反射波との合成のため、xが増大するとともに振幅が増大していくという現実にはあり得ないかもしれない解となっている。もしかしたら宇宙はこんな式で表されるのかもしれない、遠く離れるほど増大していくし。いや現実に存在する、Tesla Coilがそれだ。Tesla Coilは１次コイルで励起された二次コイルの送電端からずっと離れた受電端にだけ大気の絶縁破壊を起こす程の大電圧が発生する不思議なあれである。ひょっとしてNikola Teslaはこの解の存在を知ってTesla Coilを思いついたのかもしれない。普通に行儀のよい解き方をするとこの解は得られない。もしかしたらNikola TeslaもHeavisideの演算子法を使ったのかもしれない。もっと縦軸の範囲を広くして巨視的な観点で解を眺めると全容が見えてくる。縦軸のスケールを制限しないと、送電点側は受電端側に比べると振幅が借り着無く小さくほとんど無視できる電圧しか発生していない。線路14m長のところで送電端に給電される電圧の200万倍近い高電圧が発生している。Tesla Towerが完成していればそれ以上の規模だったかもしれない。もっともこれはL,Cがあり得ない程超高透磁率、超高誘電率の線路を使った場合の極端な例だが。実際に利用可能な線路を使って計算してみるのは読者の課題としよう( ´∀`) 以降の問題でちゃんと送電端と受電端の境界条件を与える問題が登場することを期待する。 P.S 著者と同じ解は得られなかったけれども、Tesla Coilの解を偶然に見つけてしまった。怪我の功名である。 Tesla Coilは明らかにここで見いだした解によって送電端から離れれば離れるほど電圧振幅が大きくなる定在波を生じる分布定数回路であることは確か。単純に巻き線比を極端に大きくした感応コイルではない。感応コイルは密結合トランスなので、二次側の両端子間に高電圧が発生する。Tesla Coilはそれとはまったく異なり、送電端は送電電圧以上発生せず、送電端から離れれば離れるほど高い電圧が発生するという不思議な性質も、この解によって説明される。学校でこの解を教えないのは、きっと第二のNikola Teslaを生み出さないようにするために意図的に仕組んだカリキュラムだからかもしれない。とにかく初期状態が線路の電圧と電流が全て0にしてしまえば、絶対にTesla Coilの解は出てこない。それと計算を簡単にするために,L=C=R=ω=1としたが、これもTesla Coilの解に必要な共振条件を満たす。計算機のなかった時代Nikola Tesraもきっとそうやって試し算をしてみたら、想定外のモンスター解が現れてびっくりしたかもしれない。これは実際に作って試すしかないと思ったに違いない。いや～電気ってほんと面白いですね。 Tesla Coilをまだ知らない諸君は、Googleで画像検索してみたまえ。絶対魅了されるはずだ。芥川龍之介も「歯車」で青い閃光になりたいと夢見たと書いてある。

webadm	投稿日時: 2013-6-8 20:35
Webmaster 登録日: 2004-11-7 居住地: 投稿: 3088	続々：無損失線路次はようやく有限長の線路の問題受電端を開放にした長さlの無損失線路にt=0で直流電圧Eを加えた。電圧、電流を求めよ。というもの。これは理論の時に最初に考えたモデル。線路は無損失なのだから線路上で電力が消費されることがなく、加えられた信号は受電端の方へ一定の速度で伝わっていく。半無限長線路の場合には永遠に受電端に到達しないので無限の彼方まで信号が伝わっていくだけだった。今度は有限長で線路の終端が存在し、信号の波はそこから先へは伝わることができず、かつ終端で電力が消費されるわけでもないので、行き場を失った波は反射率1で送電端の方へ逆戻りすると考えられる。それ以降は送電端から供給される進行波と受電端でUターンした反射波が線路上で重なっていき、それが送電端まで到達する。送電端に到達した反射波は、それより先に伝わることができない。かつ線路は無損失なので、送電端で電力が消費されることもなく、終端条件によって進行波と同じ方向に逆戻りすることになる。ということで波が送電端から受電端に届く所用時間間隔で線路上に行ったり来たりする反射波と進行波が重なり合う現象が発生することが予想される。タイミングチャートで描くとということになる。進行波は送電端に電圧Eが加えられている限り途切れることはないのと、反射波も進行波が途切れない限り途切れることなく発生する。それがUターン毎に時間遅れで重なり合うことになる。ことは簡単ではない。さてこれを数式で表すことができるのだろうか？　そこがこの問題の中心ということになる。ストラテジーとしては・電信方程式の基礎方程式から無損失線路の方程式を導出・演算子法でx=0とx=lにおける境界条件を与えて方程式を解くということになるが、そう簡単にいかない予感がする。題意では送電端に加えられた電圧と、受電端が開放されているということしか条件が与えられていない。前者は送電端の電圧条件を与え、後者は受電端の電流条件を与えるが、送電端の電流条件と受電端の電圧条件は与えられていない。とりあえず基礎方程式からベクトルとHeaviside演算子を使った式を導出してみよう無損失線路なので $\begin{eqnarray}<br />\frac{\partial}{\partial x}&=&p,\ \ \ \ \frac{\partial}{\partial t}=q\\<br />{\bf U}&=&\begin{bmatrix}e\left(x,t\right)\cr i\left(x,t\right)\end{bmatrix}\\<br />\begin{bmatrix}p&L\,q\cr C\,q&p\end{bmatrix}{\bf U}&=&0\\<br />\end{eqnarray}$ ということになる。これの特性方程式から $\begin{eqnarray}<br />\left\|\begin{array}p&L\,q\cr C\,q&p\end{array}\right\|&=&p^2-L\,C\,q^2=0\\<br />\left(p^2-L\,C\,q^2\right){\bf U}&=&0<br />\end{eqnarray}$ ということになる。これを例によって距離で二重積分してUについて解くと $\begin{eqnarray}<br />\frac{1}{p}\left(\frac{1}{p}\left(p^2-L\,C\,q^2\right){\bf U}+{\bf K_0}\right)+{\bf K_1}&=&0\\<br />{\bf U}&=&-\frac{{\bf K_1}p^2+{\bf K_0}p}{p^2-L\,C\,q^2}=-{\bf K_1}\frac{p^2}{p^2-L\,C\,q^2}-\frac{{\bf K_0}}{\sqrt{L\,C}q}\frac{p\sqrt{L\,C}q}{p^2-L\,C\,q^2}\\<br />&=&-{\bf K_1}\cosh\,x\sqrt{L\,C}q-\frac{{\bf K_0}}{\sqrt{L\,C}q}\sinh\,x\sqrt{L\,C}q<br />\end{eqnarray}$ ということになる。これをx=0とx=lにおける境界条件を与えると $\begin{eqnarray}<br />\left.{\bf U}\right\|_{x=0}&=&-{\bf K_1}=\begin{bmatrix}E\cr i\left(0,t\right)\end{bmatrix}\\<br />{\bf K_1}&=&-\begin{bmatrix}E\cr i\left(0,t\right)\end{bmatrix}\\<br />\left.{\bf U}\right\|_{x=\ell}&=&-{\bf K_1}\cosh\,\ell\sqrt{L\,C}q-\frac{{\bf K_0}}{\sqrt{L\,C}q}\sinh\,\ell\sqrt{L\,C}q=\begin{bmatrix}e\left(\ell,t\right)\cr 0\end{bmatrix}\\<br />{\bf K_0}&=&-\frac{\sqrt{L\,C}q}{\sinh\,\ell\sqrt{L\,C}q}\left({\bf K_1}\cosh\,\ell\sqrt{L\,C}q+\begin{bmatrix}e\left(\ell,t\right)\cr 0\end{bmatrix}\right)=\frac{\sqrt{L\,C}q}{\sinh\,\ell\sqrt{L\,C}q}\left(\begin{bmatrix}E\cr i\left(0,t\right)\end{bmatrix}\cosh\,\ell\sqrt{L\,C}q-\begin{bmatrix}e\left(\ell,t\right)\cr 0\end{bmatrix}\right)\\<br />&=&\begin{bmatrix}E\sqrt{L\,C}q\coth\,\ell\sqrt{L\,C}q-\frac{\sqrt{L\,C}q}{\sinh\,\ell\sqrt{L\,C}q}e\left(\ell,t\right)\cr i\left(0,t\right)\sqrt{L\,C}q\,\coth\,\ell\sqrt{L\,C}q\end{bmatrix}<br />\end{eqnarray}$ ということになる。ポテンシャルK0,K1は解けたが、K0の中に未知関数e(l,t),i(0,t)が現れてしまっている。どうすんだこれ（´Д｀；）とりあえずK0,K1を途中解に代入して整理すると $\begin{eqnarray}<br />{\bf U}&=&\begin{bmatrix}E\cr i\left(0,t\right)\end{bmatrix}\cosh\,x\sqrt{L\,C}q-\frac{1}{\cancel{\sqrt{L\,C}q}}\begin{bmatrix}E\cancel{\sqrt{L\,C}q}\coth\,\ell\sqrt{L\,C}q-\frac{\cancel{\sqrt{L\,C}q}}{\sinh\,\ell\sqrt{L\,C}q}e\left(\ell,t\right)\cr i\left(0,t\right)\cancel{\sqrt{L\,C}q}\,\coth\,\ell\sqrt{L\,C}q\end{bmatrix}\sinh\,x\sqrt{L\,C}q\\<br />&=&\begin{bmatrix}E\,\left(\cosh\,x\sqrt{L\,C}q-\frac{1}{\sinh\,\ell\sqrt{L\,C}q}\left(\cosh\,\ell\sqrt{L\,C}q-\frac{e\left(\ell,t\right)}{E}\right)\sinh\,x\sqrt{L\,C}q\right)\cr i\left(0,t\right)\left(\cosh\,x\sqrt{L\,C}q-\frac{\cosh\,\ell\sqrt{L\,C}q}{\sinh\,\ell\sqrt{L\,C}q}\,\sinh\,x\sqrt{L\,C}q\right)\end{bmatrix}\\<br />&=&\begin{bmatrix}\frac{E}{\sinh\,\ell\sqrt{L\,C}q}\left(\sinh\,\ell\sqrt{L\,C}q\,\cosh\,x\sqrt{L\,C}q-\left(\cosh\,\ell\sqrt{L\,C}q-\frac{e\left(\ell,t\right)}{E}\right)\sinh\,x\sqrt{L\,C}q\right)\cr \frac{i\left(0,t\right)}{\sinh\,\ell\sqrt{L\,C}q}\left(\sinh\,\ell\sqrt{L\,C}q\,\cosh\,x\sqrt{L\,C}q-\cosh\,\ell\sqrt{L\,C}q\,\sinh\,x\sqrt{L\,C}q\right)\end{bmatrix}\\<br />&=&\begin{bmatrix}E\frac{\sinh\left(\ell-x\right)\sqrt{L\,C}q}{\sinh\,\ell\sqrt{L\,C}q}+e\left(\ell,t\right)\frac{\sinh\,x\sqrt{L\,C}q}{\sinh\,\ell\sqrt{L\,C}q}\cr i\left(0,t\right)\frac{\sinh\left(\ell-x\right)\sqrt{L\,C}q}{\sinh\,\ell\sqrt{L\,C}q}\end{bmatrix}<br />\end{eqnarray}$ ということになる。この結果から判ることは線路の電圧分布が送信端の源信号Eと受電端の反射信号源e(l,t)のそれぞれの信号端からの距離と時間に依存すること、電流分布は受電端が開放のため常に受電端から反射する電流は0であるため送電端の電流と受電端からの距離と時間のみに依存するということである。さてこれを実関数に変換するのはテクニックが要る。分子と分母の両方に時間シフト演算子があるからだ。どうすんだこれ（´Д｀；）だんだんとどの本にも書いてない式が現れてくると、孤高の世界に足を踏み入れたことを感じる。後にも先にもここに足を踏みいれたのは自分しかいないのだと思うと心細くなる。Heavisideはこれを２０年間耐えたばかりか、その結果を公表したばかりに大変な災難に巻き込まれることになったのだから、それに比べればまだましである。とにかく今は自分を信じるしかない。上の式を指数関数に展開して整理すると $\begin{eqnarray}<br />{\bf U}&=&\begin{bmatrix}E\frac{e^{\left(\ell-x\right)\sqrt{L\,C}q}-e^{-\left(\ell-x\right)\sqrt{L\,C}q}}{e^{\ell\sqrt{L\,C}q}-e^{-\ell\sqrt{L\,C}q}}+e\left(\ell,t\right)\frac{e^{x\sqrt{L\,C}q}-e^{-x\sqrt{L\,C}q}}{e^{\ell\sqrt{L\,C}q}-e^{-\ell\sqrt{L\,C}q}}\cr i\left(0,t\right)\frac{e^{\left(\ell-x\right)\sqrt{L\,C}q}-e^{-\left(\ell-x\right)\sqrt{L\,C}q}}{e^{\ell\sqrt{L\,C}q}-e^{-\ell\sqrt{L\,C}q}}\end{bmatrix}\\<br />&=&\begin{bmatrix}E\frac{e^{-x\sqrt{L\,C}q}-e^{-\left(2\,\ell-x\right)\sqrt{L\,C}q}}{1-e^{-2\,\ell\sqrt{L\,C}q}}+e\left(\ell,t\right)\frac{e^{-\left(\ell-x\right)\sqrt{L\,C}q}-e^{-\left(\ell+x\right)\sqrt{L\,C}q}}{1-e^{-2\,\ell\sqrt{L\,C}q}}\cr i\left(0,t\right)\frac{e^{-x\sqrt{L\,C}q}-e^{-\left(2\,\ell-x\right)\sqrt{L\,C}q}}{1-e^{-2\,\ell\sqrt{L\,C}q}}\end{bmatrix}\\<br />&=&\begin{bmatrix}\left(E\left(e^{-x\sqrt{L\,C}q}-e^{-\left(2\,\ell-x\right)\sqrt{L\,C}q}\right)+e\left(\ell,t\right)\left(e^{-\left(\ell-x\right)\sqrt{L\,C}q}-e^{-\left(\ell+x\right)\sqrt{L\,C}q}\right)\right)\sum_{n=0}^{\infty}e^{-2\,n\,\ell\sqrt{L\,C}q}\cr i\left(0,t\right)\left(e^{-x\sqrt{L\,C}q}-e^{-\left(2\,\ell-x\right)\sqrt{L\,C}q}\right)\sum_{n=0}^{\infty}e^{-2\,n\,\ell\sqrt{L\,C}q}\end{bmatrix}\\<br />&=&\begin{bmatrix}\sum_{n=0}^{\infty}\left(E\,\left(u\left(t-\left(2\,n\,\ell+x\right)\sqrt{L\,C}\right)-u\left(t-\left(2\,\left(n+1\right)\ell-x\right)\sqrt{L\,C}\right)\right)+e\left(\ell,t-\left(\left(2\,n+1\right)\ell-x\right)\sqrt{L\,C}\right)-e\left(\ell,t-\left(\left(2n+1\right)\ell+x\right)\sqrt{L\,C}\right)\cr \sum_{n=0}^{\infty}\left(i\left(0,t-\left(2\,n\,\ell+x\right)\sqrt{L\,C}\right)-i\left(0,t-\left(2\left(n+1\right)\ell-x\right)\sqrt{L\,C}\right)\end{bmatrix}<br />\end{eqnarray}$ ということになる。良くみるとHeavisideが導出した電信方程式の一般解に良く似ている。演算子法を使うと終端の２つの信号源から伝搬する信号に関する無限級数で表されることになる。これは著者の解とは大分見かけが違う。またしてもHeavisideが歩んだ２０年の小径を横切った感じがする。１９世紀にHeavisideは確かにこの辺を通り過ぎたのだ。上の解を最初に描いたジグザグ図に当てはめてみると。ということになる。これをグラフにプロットするのもやっかいである。 0＜x＜lとして、線路上の点xにおける電圧をプロットしてみるとということになる。線路は無損失なので、このパターンが延々に繰り返されることになる。電流について同様にジグザグ図を描いてみると i(0,0)=E/Z0として電流をプロットしてみるとということになる。解が再帰的にe(l,t)やi(0,t)を参照しているので時間を遡った境界点の電圧と電流値に依存する。時間と距離だけ与えられてぱっと計算するというわけにはいかず、最低でもt=0から１周期分は計算しておく必要がある。周期関数なので、 $\begin{eqnarray}<br />e\left(x,t\right)&=&e\left(x,t+\frac{4\,\ell}{c}\right)\\<br />i\left(x,t\right)&=&i\left(x,t+\frac{4\,\ell}{c}\right)<br />\end{eqnarray}$ ということになる。再帰的でない定式化は読者の課題としよう( ´∀`) P.S 演算子法で解くと、式の上ではまるで時間軸を逆方向に移動する波があるようにも解釈できる。d'Alembertの解には二組あるというのを前にも書いたが、ひとつは空間上（x軸上）を移動する進行波と反射波、もう一組が時間軸上を移動する進行波と反射波の解である。１９世紀には時間軸上を逆戻りする波などナンセンスとされたが、２０世紀に入って量子力学によってそれは必須となった。

(1) 2 3 »

スレッド表示 | 新しいものから

前のトピック | 次のトピック | トップ

投稿するにはまず登録を