処理手順

直線状アンテナの放射インピーダンス

$$Z = \frac{V_0}{I(0)}$$

を計算する処理手順を以下に示す。 ただし、電流の展開項数N0は奇数とする。

1. 計算に必要なパラメータを読み込む:
   • $2L/\lambda_0$, $r_0/\lambda_0$,N0
   • $f_0$を基準周波数として、$f_s/f_0$から$f_e/f_0$まで、 $\Delta f/f_0$きざみでインピーダンスを求めることにする。
2. 周波数を$f_s/f_0$から$f_e/f_0$まで変えながら以下の処理を繰り返す。
   1. 計算に必要なパラメータを決定する: $kL$, $kr_0$, $k\Delta=2kL/(N0+1)$
   2. マトリックス要素$\ell_{mn}$を計算する:
      
      $$\mbox{\tt C(m,n)} \leftarrow \ell_{mn}$$
      
      

   3. 同様に$g_m$を計算する:
      
      $$\mbox{\tt C(m,N0+1)} \leftarrow g_m$$
      
      
      ここで、連立方程式を解くサブルーチンの都合で$g_m$の値を C(m,N0+1)に代入している。
   4. サブルールチンを利用して連立方程式を解く。

      C(m,N0+1) $, (m=1,2,3,\cdots\mbox{\tt N0})$に電流の 展開関数が返される。

   5. 放射インピーダンスを Z = 1/C(N0/2+1,N0+1) より求める。