5.8 RWG基底関数を用いたポテンシャル積分について
RWG基底関数を用いたベクトルポテンシャルおよびスカラポテンシャルにおいて，観測点と波源が一致する特異点を含めた積分について説明する$^\dagger$．
$\dagger$ D. R. Wilton, S. M. Rao, A. W. Glisson, D. H. Schaubert, O. M. Al-Bundak, and C. M. Butler,
“Potential integrals for uniform and linear source distributions on polygonal and polyhedral domains,”
IEEE Trans. Antennas Propagat. , vol.32, no.3, pp.276-281 (1984).
三角形領域$T_q$おける電流3成分$(i=1,2,3)$から領域$T_p$の中心点（位置ベクトル$\VEC{r}_p^{(c)}$）への寄与をまとめて考えると，
これらの積分は次のようになる．
\begin{eqnarray}
\VEC{A}_{pq,i} &=& \frac{\mu}{4\pi} \iint _{T_q} \left( \frac{l_{q,i}}{2A_q} \VECi{\rho}_{q,i} \right) \frac{e^{-jkR_p}}{R_p} dS'
\ \ \ \ \ \ (i=1,2,3) \\
\Phi_{pq,i} &=& \mp \frac{1}{4\pi j \omega \epsilon} \iint _{T_q} \left( \frac{l_{q,i}}{A_q} \right) \frac{e^{-jkR_p}}{R_p} dS'
\ \ \ \ \ \ (i=1,2,3)
\end{eqnarray}
ここで，
\begin{eqnarray}
R_p &=& \big| \VEC{r}_p^{(c)} - \VEC{r}' \big|
\\
\VECi{\rho}_{q,i} &=& \pm (\VEC{r}' - \VEC{v}_{q,i}) \ \ \ \ \ \ (i=1,2,3)
\end{eqnarray}
ただし，$ \VEC{v}_{q,i}$および$l_{q,i} \ (i=1,2,3)$は，三角形領域$T_q$の頂点の位置ベクトルおよびその対辺の長さ，
$A_q$はその面積を示す．また，上式の$\pm$の符号は，領域$T_q$が$T^{(\pm)}$のどちらかで対応して決まる．
各辺$i$（$i=1,2,3$）に関する計算に用いる変数の定義
簡単のため，$i=1$を考え，観測点を一般化して$\VEC{r}$とおく．また，$T_q$の符号は$(+)$について考え，各変数の添字を略すとベクトルポテンシャル$\VEC{A}$は，
\begin{gather}
\VEC{A} = \frac{\mu}{4\pi} \frac{l}{2A} \iint _T \VECi{\rho}_m^{(+)} \frac{e^{-jkR}}{R} dS'
\end{gather}
上式の積分項 $\VEC{I}$ を，次のような項に分けて計算していく．
\begin{eqnarray}
\VEC{I}
&=& \iint _T \VECi{\rho}_m^{(+)} \frac{e^{-jkR}}{R} dS'
\nonumber \\
&=& \iint _T (\VECi{\rho}' - \VECi{\rho}_v ) \frac{e^{-jkR}}{R} dS'
\nonumber \\
&=& \iint _T (\VECi{\rho}' - \VECi{\rho}) \frac{e^{-jkR}}{R} dS'
+ (\VECi{\rho} - \VECi{\rho}_v ) \iint _T \frac{e^{-jkR}}{R} dS'
\end{eqnarray}
ただし，$\VECi{\rho}$，$\VECi{\rho}'$，$\VECi{\rho}_v$は，
観測点の位置ベクトル$\VEC{r}$，電流源の位置ベクトル$\VEC{r}'$，
共有する辺に対向する三角形の頂点$\VEC{v}_m^{(+)}$を各々面$S$に投影したベクトルを示す．
RWG基底関数に関わる変数の定義
さらに，
\begin{gather}
\frac{e^{-jkR}}{R} = \left( \frac{e^{-jkR}}{R} -\frac{1}{R} \right) + \frac{1}{R} \\
\lim _{R\to 0} \left( \frac{e^{-jkR}}{R} -\frac{1}{R} \right) = -jk
\end{gather}
を考慮して，
\begin{eqnarray}
\VEC{I}
&=& \iint _T (\VECi{\rho}' - \VECi{\rho}) \frac{e^{-jkR}-1}{R} dS'
+ \iint _T \frac{\VECi{\rho}' - \VECi{\rho}}{R} dS'
\nonumber \\
&&+ (\VECi{\rho} - \VECi{\rho}_v ) \iint _T \frac{e^{-jkR}-1}{R} dS'
+ (\VECi{\rho} - \VECi{\rho}_v ) \iint _T \frac{dS'}{R}
\label{eq:I}
\end{eqnarray}
上式の第1項および第3項は特異点を持たないため，容易に数値積分できる．
3角形領域の積分について
まず，三角形領域$T_q$における波源の位置ベクトル$\VEC{r}'$は，
\begin{gather}
\VEC{r}' = \VEC{v}_{q,1} + \xi _q ( \VEC{v}_{q,2} - \VEC{v}_{q,1}) + \eta _q ( \VEC{v}_{q,3} - \VEC{v}_{q,1} )
\ \ \ \ \ (0 \leq \xi _q, \eta _q \leq 1)
\end{gather}
簡単のため，添字$q$は省略して，
\begin{eqnarray}
\VEC{r}'
&=& \VEC{v}_1 + \xi (\VEC{v}_2 - \VEC{v}_1 ) + \eta (\VEC{v}_3 - \VEC{v}_1)
\ \ \ \ \ (0 \leq \xi, \eta \leq 1)
\nonumber \\
&=& (1-\xi -\eta ) \VEC{v}_1 + \xi \VEC{v}_2 + \eta \VEC{v}_3
\end{eqnarray}
これより，面積要素$dS'$は，
\begin{eqnarray}
dS'
&=& \left| \frac{\partial \VEC{r}'}{\partial \xi} \times \frac{\partial \VEC{r}'}{\partial \eta} \right| d\xi d\eta
\nonumber \\
&=& \big| (\VEC{v}_2 - \VEC{v}_1 ) \times (\VEC{v}_3 - \VEC{v}_1 ) \big| d\xi d\eta
\nonumber \\
&=& 2 A_q d\xi d\eta
\end{eqnarray}
よって，三角形領域$T_q$の積分は，次のような2重積分の形で表される．
\begin{gather}
\iint _{T_q} g(\VEC{r}') dS'
= 2 A_q \int _{\xi=0}^1 \int _{\eta=0}^{1-\xi} g \Big( (1-\xi -\eta ) \VEC{v}_1 + \xi \VEC{v}_2 + \eta \VEC{v}_3 \Big) d\eta d\xi
\end{gather}
3角形領域の面積分に関わる変数の定義
特異点を含む積分について
式\eqref{eq:I}の第2項および第4項の積分$I_2$，$I_4$
\begin{gather}
I_2 = \iint _T \frac{\VECi{\rho}' - \VECi{\rho}}{R} dS', \ \ \ \ \ \
I_4 = \iint _T \frac{dS'}{R}
\end{gather}
については，周回積分への変換と特異点のふるまいを考慮した積分を行う．まず，準備として，導体面に沿う2次元演算子
$\nabla _s$
を作用させたときの計算を行う．
\begin{align}
&\VEC{P} = P \VEC{u}_p \equiv \VECi{\rho}' - \VECi{\rho}
\\
&P = |\VECi{\rho}' - \VECi{\rho} |
\\
&R = |\VEC{r} - \VEC{r}' |
= |d\VEC{n} + \VECi{\rho} - \VECi{\rho}' |
= |d\VEC{n} - P \VEC{u}_p |
= \sqrt{d^2 + P^2}
\\
&\nabla' _s R
= \frac{\partial R}{\partial P} \VEC{u}_p = \frac{P}{R} \VEC{u}_p
= \frac{\VEC{P}}{R}
= \frac{\VECi{\rho}' - \VECi{\rho}}{R}
\end{align}
これより，積分$I_2$に対して，特異点$P=0$の点を含む微小領域を$s_\epsilon $として，次のように領域$s_\epsilon $とそれ以外に分けて積分し，それから極限を求める．
\begin{eqnarray}
I_2
&=& \iint _S \frac{\VECi{\rho}' - \VECi{\rho}}{R} dS'
\nonumber \\
&=& \lim _{s_\epsilon \to 0} \iint _{S-s_\epsilon} \nabla _s' R dS'
+ \lim _{s_\epsilon \to 0} \iint _{s_\epsilon} \frac{\VECi{\rho}' - \VECi{\rho}}{R} dS'
\label{eq:I2}
\end{eqnarray}
また，
\begin{gather}
\frac{R}{P} = \frac{\sqrt{d^2+P^2}}{P}
= \sqrt{\frac{d^2}{P^2} +1}
\end{gather}
\begin{eqnarray}
\frac{\partial}{\partial P} \left( \frac{R}{P} \right)
&=& \frac{\partial}{\partial P} \left( \sqrt{\frac{d^2}{P^2} +1} \right)
\nonumber \\
&=& -\frac{d^2}{P^3} \frac{P}{\sqrt{d^2+P^2}}
= -\frac{d^2}{P^2} \frac{1}{R}
\end{eqnarray}
\begin{eqnarray}
\nabla'_s \cdot \left( \frac{R}{P} \VEC{u}_\rho \right)
&=& \nabla'_s \frac{R}{P} \cdot \VEC{u}_\rho + \frac{R}{P} \nabla'_s \cdot \VEC{u}_\rho
\nonumber \\
&=& \frac{\partial}{\partial P} \left( \frac{R}{P} \right) + \frac{R}{P} \frac{1}{P}
\nonumber \\
&=& -\frac{d^2}{P^2} \frac{1}{R} + \frac{R}{P^2}
\nonumber \\
&=& \frac{-d^2+R^2}{P^2 R}
= \frac{P^2}{P^2 R}
= \frac{1}{R}
\end{eqnarray}
同様にして，積分$I_4$も特異点$P=0$の点を含む微小領域を$s_\epsilon $として，
\begin{eqnarray}
I_4
&=& \iint _T \frac{dS'}{R}
\nonumber \\
&=& \lim _{s_\epsilon \to 0} \iint _{S-s_\epsilon } \nabla _s' \cdot \left( \frac{R\VEC{P}}{P^2} \right) dS'
+ \lim _{s_\epsilon \to 0} \iint _{s_\epsilon} \frac{dS'}{R}
\label{eq:I4}
\end{eqnarray}
と分けて考えることにする．
周回積分への変換
式\eqref{eq:I4}の第1項は$P=0$の点を除いた積分であり，次のようにガウスの2次元発散定理により周回積分に変換できる．
\begin{eqnarray}
&&\lim _{\epsilon \to 0} \iint _{S-s_\epsilon } \nabla _s' \cdot \left( \frac{R\VEC{P}}{P^2} \right) dS'
\nonumber \\
&=& \lim _{\epsilon \to 0} \oint _{C-c_\epsilon} \frac{R\VEC{P}}{P^2} \cdot d\VEC{u}'
\nonumber \\
&=& \oint _{C} \frac{R}{P^2} \big( \VEC{P} \cdot \VEC{u}' \big) dl'
+ \lim _{\epsilon \to 0} \oint _{-c_\epsilon} \frac{R}{P} \big( \VEC{u}_p \cdot \VEC{u}' \big) dl'
\end{eqnarray}
ただし，周回積分路$C$は面$S$の周辺，
$c_\epsilon$は面$s_\epsilon$の周辺の各々閉じた経路を示す．また，ベクトル線積分要素
$d\VEC{u}' (= \VEC{u}' dl')$の方向は，面$S$，$s_\epsilon$上で積分路$C$，$c_\epsilon$の法線方向にとる．上式の第1項は，積分路$C$が3角形領域の辺$i$に沿う直線の積分路$L_i$となるので，
\begin{gather}
\oint _{C} \frac{R}{P^2} \big( \VEC{P} \cdot \VEC{u}' \big) dl'
= \sum _i \big( \VEC{P}^{(0)} \cdot \VEC{u}_i' \big) \int _{L_i} \frac{R}{P^2} dl'
\end{gather}
ここで，
\begin{eqnarray}
P^2 &=& (P^{(0)})^2 + l'^2
\\
R^2 &=& d^2+(P^{(0)})^2 + l'^2 = d^2+P^2
\end{eqnarray}
より，
\begin{eqnarray}
\frac{R}{P^2}
&=& \frac{R}{P^2} \frac{R}{R}
= \frac{P^2+d^2}{P^2R}
\nonumber \\
&=& \frac{1}{R} + \frac{d^2}{P^2R}
\nonumber \\
&=& \frac{1}{\sqrt{d^2+(P^{(0)})^2 + l'^2}} + \frac{d^2}{\big\{ (P^{(0)})^2 + l'^2 \big\} ^2 \sqrt{d^2+(P^{(0)})^2 + l'^2}}
\end{eqnarray}
辺$i$（積分路$L_i$）の両端の座標成分を$l=l_i^{(-)}, l_i^{(+)}$とすると，
\begin{eqnarray}
\int _{L_i} \frac{R}{P^2} dl'
&=& \int _{L_i} \left( \frac{1}{R} + \frac{d^2}{P^2R} \right) dl'
\nonumber \\
&=& \int _{l_i^{(-)}}^{l_i^{(+)}} \frac{dl'}{\sqrt{d^2+\big( P_i^{(0)}\big) ^2 + l'^2}}
\nonumber \\
&&+ d^2 \int _{l_i^{(-)}}^{l_i^{(+)}} \frac{dl'}{\left\{ \big( P_i^{(0)}\big) ^2 + l'^2 \right\} \sqrt{d^2+\big( P_i^{(0)}\big) ^2 + l'^2}}
\label{eq:intRP2}
\end{eqnarray}
辺$i$に沿う線積分に関わる変数の定義（$i=1$のとき）
上式は，置換積分（$t \equiv x + \sqrt{x^2 + a^2}$）より得られる不定積分公式（積分定数$C$は省略）
\begin{gather}
\int \frac{dx}{\sqrt{x^2+a^2}} = \ln \big| x + \sqrt{x^2 + a^2} \big|
\end{gather}
を用い，
\begin{gather}
R_i^{(\pm)} = \sqrt{d^2+\big( P_i^{(0)} \big) ^2 + \left( l_i^{(\pm)} \right) ^2 }
\end{gather}
とおいて，式\eqref{eq:intRP2}の第1項は，
\begin{eqnarray}
&&\int _{l_i^{(-)}}^{l_i^{(+)}} \frac{dl'}{\sqrt{d^2+\big( P_i^{(0)}\big) ^2 + l'^2}}
\nonumber \\
&=& \left[ \ln \left| l'+ \sqrt{d^2+\big( P_i^{(0)}\big) ^2 + l'^2 } \right| \right] _{l_i^{(-)}}^{l_i^{(+)}}
\nonumber \\
&=& \ln \left| l_i^{(+)} + \sqrt{d^2+\big( P_i^{(0)}\big) ^2 + \left( l_i^{(+)} \right) ^2 } \right|
\nonumber \\
&&- \ln \left| l_i^{(-)} + \sqrt{d^2+\big( P_i^{(0)}\big) ^2 + \left( l_i^{(-)} \right) ^2 } \right|
\nonumber \\
&=& \ln \left| l_i^{(+)} + R_i^{(+)} \right| - \ln \left| l_i^{(-)} + R_i^{(-)} \right|
\nonumber \\
&=& \ln \left| \frac{l_i^{(+)} + R_i^{(+)}}{l_i^{(-)} + R_i^{(-)}} \right|
\end{eqnarray}
また，置換積分（$|q| x \equiv p \sqrt{x^2+p^2+q^2} \ \tan t $）より得られる不定積分公式（積分定数$C$は省略）
\begin{gather}
\int \frac{dx}{\big( x^2 + p^2 \big) \sqrt{x^2+p^2+q^2}}
= \frac{1}{|pq|} \tan ^{-1} \left( \frac{|q|x}{p \sqrt{x^2+p^2+q^2}} \right)
\end{gather}
を用いると，式\eqref{eq:intRP2}の第2項は，
\begin{eqnarray}
&&d^2 \int _{l'=l_i^{(-)}}^{l_i^{(+)}} \frac{dl'}{\left\{ \big( P_i^{(0)} \big) ^2 + l'^2 \right\}
\sqrt{d^2+\big( P_i^{(0)} \big) ^2 + l'^2}}
\nonumber \\
&=& \frac{d^2}{P_i^{(0)} |d|} \left[ \tan ^{-1} \left( \frac{|d| \ l'}{P_i^{(0)} \sqrt{d^2+\big( P_i^{(0)} \big) ^2+l'^2}} \right) \right] _{l'=l_i^{(-)}}^{l_i^{(+)}}
\nonumber \\
&=& \frac{|d|}{P_i^{(0)}} \left\{ \tan ^{-1} \left( \frac{|d|l_i^{(+)}}{P_i^{(0)} R_i^{(+)}} \right)
- \tan ^{-1} \left( \frac{|d|l_i^{(-)}}{P_i^{(0)} R_i^{(-)}} \right) \right\}
\end{eqnarray}
積分路$-c_\epsilon$においては，$\VEC{u}_p \cdot \VEC{u}' = -1$，$dl' = Pd\phi$，$R = \sqrt{d^2+\epsilon ^2}$より，
\begin{eqnarray}
\lim _{\epsilon \to 0} \oint _{-c_\epsilon} \frac{R}{P} \big( \VEC{u}_p \cdot \VEC{u}' \big) dl'
&=& \lim _{\epsilon \to 0} \int _0^{\alpha} \frac{R}{P} (-1) P d\phi
\nonumber \\
&=& \lim _{\epsilon \to 0}(-R) \int _0^{\alpha} d\phi
\nonumber \\
&=& \lim _{\epsilon \to 0} -\sqrt{d^2+\epsilon ^2} \cdot \alpha
\nonumber \\
&=& -\alpha |d|
\end{eqnarray}
一方，式\eqref{eq:I4}の$I_4$の第2項は，$P=0$の点を座標原点にしたローカルな極座標系を考え，積分範囲$s_\epsilon$を半径$\epsilon $，平面角$\alpha$ [rad]の円形領域にとり，
$t \equiv d^2 + P^2$の置換積分を実行すると，次のようになる．
\begin{eqnarray}
\lim _{\epsilon \to 0} \iint _{s_\epsilon} \frac{dS'}{R}
&=& \lim _{\epsilon \to 0} \int _{\phi =0}^\alpha \int _{P=0}^\epsilon \frac{PdP d\phi}{\sqrt{d^2+P^2}}
\nonumber \\
&=& \lim _{\epsilon \to 0} \Big[ \phi \Big] _0^\alpha \int \frac{1}{\sqrt{t}} \frac{dt}{2}
\nonumber \\
&=& \lim _{\epsilon \to 0} \alpha \Big[ \sqrt{t} \Big]
\nonumber \\
&=& \lim _{\epsilon \to 0} \alpha \Big[ \sqrt{d^2+P^2} \Big] _0^\epsilon
\nonumber \\
&=& \lim _{\epsilon \to 0} \alpha \big( \sqrt{d^2+\epsilon^2}-|d| \big)
= 0
\end{eqnarray}
よって，
\begin{eqnarray}
I_4 &=& \iint _T \frac{dS'}{R}
\nonumber \\
&=& - \alpha |d| + \sum _i \big( \VEC{u}_{P0,i} \cdot \VEC{u}_i' \big)
\nonumber \\
&& \cdot \left[ P_i^{(0)} \ln \left| \frac{l_i^{(+)} + R_i^{(+)}}{l_i^{(-)} + R_i^{(-)}} \right| \right.
\nonumber \\
&&\left. + |d| \left\{ \tan ^{-1} \left( \frac{|d|l_i^{(+)}}{P_i^{(0)} R_i^{(+)}} \right)
- \tan ^{-1} \left( \frac{|d|l_i^{(-)}}{P_i^{(0)} R_i^{(-)}} \right) \right\} \right]
\end{eqnarray}
ただし，
\begin{gather}
\alpha = \sum _i \big( \VEC{u}_{P0,i} \cdot \VEC{u}_i' \big)
\left\{ \tan ^{-1} \left( \frac{l_i^{(+)}}{P_i^{(0)}} \right)
- \tan ^{-1} \left( \frac{l_i^{(-)}}{P_i^{(0)}} \right) \right\}
\end{gather}
また，
\begin{eqnarray}
\beta &\equiv& \tan ^{-1} \frac{l_i^{(\pm)}}{P_i^{(0)}}
\\
\gamma &\equiv& \tan ^{-1} \frac{|d| l_i^{(\pm)}}{P_i^{(0)} R_i^{(\pm)}}
\end{eqnarray}
とおき，正接の加法定理
\begin{gather}
\tan (\beta -\gamma )= \frac{\tan \beta -\tan \gamma}{1+\tan \beta \tan \gamma}
\end{gather}
を用いると，
\begin{eqnarray}
\big( l_i^{(\pm)} \big) ^2 &=& \big( P_i^{(\pm)} \big) ^2 - \big( P_i^{(0)} \big) ^2
\\
\big( P_i^{(\pm)} \big) ^2 &=& \big( R_i^{(\pm)} \big) ^2 - \big| d \big| ^2
\\
\big( P_i^{(0)} \big) ^2 &=& \big( R_i^{(0)} \big) ^2 - \big| d \big| ^2
\end{eqnarray}
より，次の関係が得られる（導出省略）．
\begin{eqnarray}
&&\tan ^{-1} \left( \frac{l_i^{(\pm)}}{P_i^{(0)}} \right)
- \tan ^{-1} \left( \frac{|d| l_i^{(\pm)}}{P_i^{(0)} R_i^{(\pm)}} \right)
\nonumber \\
&=& \tan ^{-1} \left( \frac{P_i^{(0)} l_i^{(\pm)}}{\big( R_i^{(0)} \big) ^2 + |d| R_i^{(\pm)}} \right)
\end{eqnarray}
これより，
\begin{eqnarray}
I_4 &=& \iint _T \frac{dS'}{R}
\nonumber \\
&=& \sum _i \big( \VEC{u}_{P0,i} \cdot \VEC{u}_i' \big)
\left[ P_i^{(0)} \ln \left| \frac{l_i^{(+)} + R_i^{(+)}}{l_i^{(-)} + R_i^{(-)}} \right| \right.
\nonumber \\
&&- |d| \left\{ \tan ^{-1} \left( \frac{P_i^{(0)} l_i^{(+)}}{ \big( R_i^{(0)} \big) ^2 + |d| R_i^{(+)}} \right) \right.
\nonumber \\
&&\left. \left. - \tan ^{-1} \left( \frac{P_i^{(0)} l_i^{(-)}}{ \big( R_i^{(0)} \big) ^2 + |d| R_i^{(-)}} \right) \right\} \right]
\end{eqnarray}
一方，積分$I_2$についても，式\eqref{eq:I2}の第1項は2次元勾配定理より，周回積分（積分路$C$および$c_\epsilon$）に変換でき，第2項は$I_4$の計算と同様にしてゼロになる（導出省略）．よって，
\begin{eqnarray}
I_2
&=& \lim _{s_\epsilon \to 0} \iint _{S-s_\epsilon} \nabla _s' R dS'
\nonumber \\
&=& \lim _{\epsilon \to 0} \oint _{C-c_\epsilon} R \VEC{u}' dl'
\nonumber \\
&=& \sum _i \VEC{u}_i' \int _{L_i} R dl'
+ \lim _{\epsilon \to 0} \int _{-c_\epsilon} R \VEC{u}' dl'
\nonumber \\
&=& \sum _i \VEC{u}_i' \int _{l_i^{(-)}}^{l_i^{(+)}} \sqrt{d^2+\big( P_i^{(0)}\big) ^2 + l'^2 } \ dl'
+ \lim _{\epsilon \to 0} R \int _0^\alpha \VEC{u}' \epsilon d\phi
\end{eqnarray}
上式の第2項はゼロ，
第1項は置換積分（$t \equiv x + \sqrt{x^2 + a^2}$）より得られる不定積分公式（積分定数$C$は省略）
\begin{gather}
\int \sqrt{x^2+a^2} \ dx = \frac{1}{2} \left( x \sqrt{x^2 + a^2} + a^2 \ln \big| x+\sqrt{x^2 + a^2} \big| \right)
\end{gather}
より，
\begin{eqnarray}
I_2
&=& \sum _i \VEC{u}_i' \frac{1}{2}
\left[ l' \sqrt{d^2+\big( P_i^{(0)}\big) ^2 + l'^2 } \right.
\nonumber \\
&&\left.
+ \left\{ d^2+\big( P_i^{(0)}\big) ^2 \right\} \ln \left| l' + \sqrt{d^2+\big( P_i^{(0)}\big) ^2 + l'^2 } \right|
\right] _{l_i^{(-)}}^{l_i^{(+)}}
\nonumber \\
&=& \frac{1}{2} \sum _i \VEC{u}_i'
\Big[ \left\{ l^{(+)} R_i^{(+)} + \big( R_i^{(0)} \big) ^2 \ln \big| l_i^{(+)} + R_i^{(+)} \big| \right\}
\nonumber \\
&&- \left\{ l^{(-)} R_i^{(-)} + \big( R_i^{(0)} \big) ^2 \ln \big| l_i^{(-)} + R_i^{(-)} \big| \right\} \Big]
\nonumber \\
&=& \frac{1}{2} \sum _i \VEC{u}_i' \left\{ \big( R_i^{(0)} \big) ^2
\ln \left| \frac{l_i^{(+)} + R_i^{(+)}}{l_i^{(-)} + R_i^{(-)}} \right|
+ l_i^{(+)} R_i^{(+)} - l_i^{(-)} R_i^{(-)} \right\}
\end{eqnarray}
前のページに戻る
「電磁波工学特論」へ戻る
