スロットアンテナの理論と設計 — バビネの原理から放射特性を理解する

金属板にナイフで細長い溝(スロット)を切り込むと、そこから電磁波が放射される — 直感的には不思議に思えるかもしれません。金属は電磁波を遮蔽するものなのに、なぜ穴を開けただけで「アンテナ」として機能するのでしょうか?

この問いに答える鍵がバビネの原理(Babinet’s principle)です。バビネの原理を理解すると、スロットアンテナはダイポールアンテナの「相補的な双子」であることがわかります。そして、ダイポールアンテナの性質がわかっていれば、スロットアンテナの放射パターンやインピーダンスも自動的に導けるのです。

スロットアンテナの理論を理解することは、以下の場面で大きな力を発揮します。

  • 航空機のアンテナ: 航空機の機体表面にスロットを切ることで、空気抵抗を増やさずにアンテナを設置できます。機体そのものが導体板の役割を果たすため、スロットアンテナは航空宇宙分野で広く使われています
  • 導波管スロットアレイ: 導波管の壁面に複数のスロットを配列することで、高利得・低損失のアンテナアレイが構成できます。レーダーや衛星通信の送受信アンテナとして不可欠です
  • マイクロストリップアンテナの理解: パッチアンテナの放射メカニズムはスロットモデルで理解できるため、スロットアンテナの理論はパッチアンテナ設計の基礎にもなります

本記事の内容

  • スロットアンテナの基本構造と動作の直感的理解
  • バビネの原理の物理的意味と数学的定式化
  • スロットダイポールとワイヤダイポールの等価関係
  • インピーダンスの相補関係式の導出
  • 導波管スロットアレイの設計原理
  • Pythonによるダイポールとスロットの放射パターン比較

前提知識

この記事を読む前に、以下の記事を読んでおくと理解が深まります。

スロットアンテナとは

基本構造

スロットアンテナの最もシンプルな形態は、無限に広い完全導体板に切り込んだ細長い開口(スロット)です。スロットの長さを $L$、幅を $W$($W \ll L$)とします。このスロットの中央部に高周波電源を接続し、スロットの短辺方向に電圧を加えると、スロット内に電界が生じ、電磁波が空間に放射されます。

日常的なアナロジーで考えると、スロットアンテナは「音叉の振動」に似ています。音叉の隙間(スリット)から音波が放射されるように、金属板のスロットから電磁波が放射されるのです。

なぜスロットから電磁波が放射されるのか

無限導体板に電磁波が入射すると、板の表面に誘導電流が流れ、入射波を打ち消す散乱波が生じます。これが金属板による遮蔽の原理です。

ところが、スロットがあると、スロットを横切る方向の電流が流れることができません。電流は遮断され、スロットの両縁に沿って迂回します。この電流の不連続(ギャップ)が、スロット内に電界を生み出します。スロット内の電界は時間的に振動し、マクスウェル方程式に従って磁界を誘導し、電磁波として空間に放射されるのです。

つまり、スロットアンテナの放射は導体板上の電流を「切断」することによって生じると理解できます。ダイポールアンテナが「電流を流す」ことで放射するのに対し、スロットアンテナは「電流を遮断する」ことで放射する — この対照性が、次に解説するバビネの原理の核心です。

バビネの原理

光学におけるバビネの原理

バビネの原理はもともと光学(スカラー回折理論)で知られた原理です。光学版のバビネの原理は次のように述べられます。

「不透明なスクリーンによる回折パターンと、そのスクリーンの相補的なスクリーン(元のスクリーンの開口部と遮蔽部を入れ替えたもの)による回折パターンの和は、スクリーンがない場合の光の振幅と一致する。」

つまり、ある障害物による回折と、その障害物の「ネガ像」による回折は、互いに補い合って元の波を再現するということです。

電磁気学版のバビネの原理

電磁気学では、バビネの原理はベクトル場に拡張されます。スカラー版との本質的な違いは、電界と磁界の役割が入れ替わることです。

いま、無限導体板に開口(スロット)$S$ があるスクリーンを考えます。このスクリーンの「相補的スクリーン」とは、もとのスクリーンの開口部分が導体で、導体部分が開口になったものです。

電磁気学版のバビネの原理は、次のように述べられます。

「元のスクリーンを通過する電磁波の電界 $\bm{E}_1$ と、相補的スクリーンを通過する電磁波の磁界 $\bm{H}_2$ は、次の関係で結ばれる。」

$$ \bm{E}_1 + \eta_0 \bm{H}_2 = \bm{E}_0 $$

ここで $\bm{E}_0$ はスクリーンがない場合(自由空間)の電界、$\eta_0 = \sqrt{\mu_0/\epsilon_0} \approx 377\,\Omega$ は自由空間の波動インピーダンスです。

この式が意味するのは、元のスクリーンの開口を通過する電界は、相補的スクリーンの開口を通過する磁界(を $\eta_0$ 倍したもの)と足すと、自由空間の電界になるということです。

電界と磁界の入れ替わり

バビネの原理の最も重要な帰結は、相補的な2つの構造では電界と磁界の役割が入れ替わるということです。

ワイヤダイポールアンテナが電流(電荷の流れ)によって電界を放射するのに対し、その相補構造であるスロットアンテナは等価磁流(磁荷の流れ)によって磁界を放射します。放射パターンの形状は同じですが、電界と磁界の偏波面が90°回転します。

この電界と磁界の交換は、単に形式的なものではなく、実用的に重要な帰結をもたらします。特に、インピーダンスの相補関係と偏波の回転は、アンテナ設計に直接影響する性質です。

では、この原理をスロットとダイポールの等価関係に具体的に適用してみましょう。

スロットダイポールとワイヤダイポールの等価関係

相補的構造

長さ $L$、幅 $W$ のスロットアンテナと、長さ $L$、直径 $d \approx W$ のワイヤダイポールアンテナを考えます。スロットの形状がダイポールの「ネガ像」になっている — つまり、導体の部分と空間の部分が入れ替わっている — とき、これらは相補的な構造です。

放射パターンの比較

バビネの原理により、スロットアンテナの放射パターンはダイポールアンテナの放射パターンと形状は同一ですが、偏波面が90°回転します。

ダイポールアンテナの電界は、ダイポールの軸を含む面(E面)内で最大となり、その面に直交する面(H面)では放射パターンが無指向性です。

スロットアンテナでは、この関係が反転します。スロットの長軸に垂直な面がE面となり、長軸を含む面がH面になります。具体的には:

  • ワイヤダイポール(z軸方向): $E_\theta$ 成分が支配的(θ偏波)
  • スロットダイポール(z軸方向): $E_\phi$ 成分が支配的(φ偏波)

つまり、ダイポールの放射パターンを90°回転させ、電界と磁界を入れ替えたものがスロットの放射パターンです。

等価磁流モデル

スロットアンテナの解析をさらに進めるために、等価磁流の概念を導入します。

ダイポールアンテナでは、導体上を電流 $\bm{J}$ が流れ、それが電磁波を放射します。スロットアンテナでは、スロット内の電界を等価な磁流 $\bm{M}$ に置き換えて解析します。

スロット内の電界が $\bm{E}_s$ であるとき、等価磁流は

$$ \bm{M} = -2\hat{\bm{n}} \times \bm{E}_s $$

で定義されます。ここで $\hat{\bm{n}}$ は導体板の法線ベクトル、係数2は無限導体板による鏡像効果を考慮したものです。

この等価磁流を用いると、スロットの放射問題は磁流源からの放射問題に帰着されます。磁流源の放射は、マクスウェル方程式の双対性(duality)により、電流源の放射問題と全く同じ数学構造で解けます。ただし、電界と磁界の役割が入れ替わります。

放射パターンの等価関係を理解したところで、次にインピーダンスの相補関係を導出しましょう。

インピーダンスの相補関係

ブッカーの関係式

バビネの原理から導かれる最も実用的な結果の一つが、相補的アンテナのインピーダンスの関係式です。この関係式はブッカーの関係式(Booker’s relation)と呼ばれます。

スロットアンテナの入力インピーダンスを $Z_s$、その相補構造であるワイヤダイポールアンテナの入力インピーダンスを $Z_d$ とすると

$$ Z_s \cdot Z_d = \frac{\eta_0^2}{4} $$

が成り立ちます。ここで $\eta_0 \approx 377\,\Omega$ は自由空間の波動インピーダンスです。

ブッカーの関係式の導出

ブッカーの関係式を導出しましょう。バビネの原理から、相補的な2つのスクリーンの透過波の電界と磁界には次の関係が成り立ちます。

自由空間を伝搬する平面波を入射波とします。入射波の電界を $\bm{E}_0$、磁界を $\bm{H}_0$ とすると、自由空間の波動インピーダンスの関係から

$$ |\bm{E}_0| = \eta_0 |\bm{H}_0| $$

です。

スロットスクリーンの透過電界を $\bm{E}_s$、相補スクリーン(ダイポール的構造)の透過磁界を $\bm{H}_d$ とすると、バビネの原理より

$$ \bm{E}_s = \eta_0 \bm{H}_d $$

の対応関係があります(スクリーンがない場合の入射波を基準として)。

次に、アンテナとして給電する場合を考えます。スロットアンテナの給電電圧を $V_s$、流れ込む等価電流を $I_s$ とし、ダイポールアンテナの給電電圧を $V_d$、電流を $I_d$ とします。

バビネの原理の双対性から、スロットの電圧がダイポールの電流に対応し、スロットの電流がダイポールの電圧に対応します。

$$ V_s \leftrightarrow \eta_0 I_d / 2, \quad I_s \leftrightarrow V_d / (2\eta_0) \cdot \eta_0 = V_d / 2 $$

ここで係数 $1/2$ はスロットが片面にのみ放射する場合と両面に放射する場合の補正です(無限導体板のスロットは両面放射)。

インピーダンスの定義 $Z = V / I$ を用いると

$$ Z_s = \frac{V_s}{I_s}, \quad Z_d = \frac{V_d}{I_d} $$

上の対応関係を代入すると

$$ Z_s \cdot Z_d = \frac{V_s}{I_s} \cdot \frac{V_d}{I_d} = \frac{\eta_0^2}{4} $$

これがブッカーの関係式です。

具体的な数値例

半波長ダイポールアンテナの入力インピーダンスは $Z_d = 73 + j42.5\,\Omega$ です(共振長よりわずかに長い場合)。共振長に調整すると $Z_d \approx 73\,\Omega$(純抵抗)となります。

ブッカーの関係式から、対応する半波長スロットのインピーダンスは

$$ Z_s = \frac{\eta_0^2}{4 Z_d} = \frac{(377)^2}{4 \times 73} \approx \frac{142129}{292} \approx 487\,\Omega $$

となります。ダイポールのインピーダンスが73 $\Omega$ であるのに対し、スロットは487 $\Omega$ と大幅に高いことがわかります。

この高いインピーダンスは、50 $\Omega$ 系への直接整合が困難であることを意味します。実用的なスロットアンテナでは、スロットの位置をずらす(オフセット給電)、スロットの形状を変える(ボウタイ型スロット)、整合回路を用いるなどの手法でインピーダンスを調整します。

自己相補構造のインピーダンス

特に興味深いのは自己相補構造の場合です。自己相補構造とは、元の構造と相補構造が完全に同一(回転や反転で一致)のものを言います。

自己相補構造では $Z_s = Z_d$ なので

$$ Z^2 = \frac{\eta_0^2}{4} $$

$$ Z = \frac{\eta_0}{2} = \frac{377}{2} \approx 188.5\,\Omega $$

この結果は周波数に依存しません。これは螺旋アンテナやボウタイアンテナなどの自己相補構造が広帯域で一定のインピーダンスを持つことの理論的根拠です。

インピーダンスの理論を理解したところで、次に実用的に非常に重要な導波管スロットアレイについて解説しましょう。

導波管スロットアレイ

導波管壁面のスロット

導波管スロットアレイは、矩形導波管の壁面にスロットを配列したアンテナです。導波管が給電線路と放射素子の支持構造を兼ねるため、非常にコンパクトで低損失なアレイアンテナを実現できます。レーダーや衛星通信で広く使われています。

矩形導波管のTE$_{10}$モードの管内波長を $\lambda_g$ とします。管内波長は自由空間の波長 $\lambda_0$ と導波管の幅広面の寸法 $a$ を用いて

$$ \lambda_g = \frac{\lambda_0}{\sqrt{1 – \left(\frac{\lambda_0}{2a}\right)^2}} $$

で与えられます。$\lambda_g > \lambda_0$ であることに注意してください。

スロットの種類と位置

導波管壁面のスロットは、その位置と方向によって異なる励振特性を持ちます。

幅広面縦スロット(broad wall longitudinal slot): 導波管の幅広面上に、導波管の長手方向に沿ってスロットを切ります。スロットを導波管の中心線からオフセットさせることで、管壁電流を切断し、放射が生じます。オフセット量が大きいほど強く励振されます。

スロットの中心線からのオフセット $x$ における規格化コンダクタンスは、近似的に

$$ \frac{g}{g_0} \approx \sin^2\left(\frac{\pi x}{a}\right) $$

で与えられます。ここで $a$ は導波管の幅広面の寸法、$g_0$ は最大コンダクタンスです。

狭壁面傾斜スロット(narrow wall inclined slot): 導波管の狭壁面に斜めにスロットを切ります。傾斜角が大きいほど強く励振されます。航空機レーダーなどで使われます。

アレイの設計

導波管スロットアレイの設計では、各スロットの位相を揃えることが重要です。

TE$_{10}$モードの壁面電流は $\lambda_g / 2$ ごとに方向が反転します。したがって、$\lambda_g / 2$ 間隔でスロットを配列すると、隣接するスロットは逆位相で励振されます。

同位相で放射するためには、隣接するスロットを導波管中心線に対して交互に反対側にオフセットさせます。これにより、電流の方向反転とオフセットの反転が打ち消し合い、すべてのスロットが同位相で放射されます。

スロット間隔 $d = \lambda_g / 2$ のアレイのアレイファクタは

$$ \text{AF}(\theta) = \sum_{n=0}^{N-1} w_n \, e^{j n k_0 d \sin\theta} $$

で表されます。ここで $w_n$ は各スロットの重み(励振振幅)、$k_0 = 2\pi/\lambda_0$ は自由空間の波数です。

全スロットを等振幅($w_n = 1$)で励振した場合のアレイファクタは

$$ \text{AF}(\theta) = \frac{\sin\left(\frac{N k_0 d \sin\theta}{2}\right)}{\sin\left(\frac{k_0 d \sin\theta}{2}\right)} $$

となります。この式の最大値は $\theta = 0$(ブロードサイド方向)で $N$ となり、アレイの利得はスロット数 $N$ に比例して増加します。

導波管スロットアレイの理論を押さえたところで、Pythonを使ってダイポールとスロットの放射パターンを比較し、バビネの原理を目で確認しましょう。

Pythonによるダイポールとスロットの放射パターン比較

バビネの原理から、ダイポールとスロットの放射パターンは形状が同じで偏波面が90°回転する関係にあります。これをPythonで可視化します。

import numpy as np
import matplotlib.pyplot as plt

# --- 半波長ダイポールの放射パターン ---
theta = np.linspace(0.001, np.pi, 500)  # 0からπまで(0とπの特異点を避ける)

# 半波長ダイポールのE面パターン
# F(theta) = cos(pi/2 * cos(theta)) / sin(theta)
F_dipole = np.abs(np.cos(np.pi / 2 * np.cos(theta)) / np.sin(theta))
F_dipole /= np.max(F_dipole)  # 正規化

# スロットアンテナのパターン(バビネの原理より同一形状)
F_slot = F_dipole.copy()

fig, axes = plt.subplots(1, 3, figsize=(18, 5))

# (a) ダイポールのE面パターン
ax = axes[0]
theta_deg = np.degrees(theta)
F_dB = 20 * np.log10(F_dipole + 1e-10)
F_dB = np.clip(F_dB, -30, 0)
ax.plot(theta_deg, F_dB, 'c-', linewidth=2, label='Dipole E-plane')
ax.set_xlabel('Theta [degrees]', fontsize=12)
ax.set_ylabel('Normalized Pattern [dB]', fontsize=12)
ax.set_title('Half-wave Dipole\nE-plane Pattern', fontsize=13)
ax.set_xlim(0, 180)
ax.set_ylim(-30, 3)
ax.grid(True, alpha=0.3)
ax.axhline(y=-3, color='r', linestyle='--', alpha=0.5, label='-3 dB')
ax.legend()

# (b) スロットのH面パターン(ダイポールのE面と同じ形状)
ax = axes[1]
ax.plot(theta_deg, F_dB, 'm-', linewidth=2, label='Slot H-plane')
ax.set_xlabel('Theta [degrees]', fontsize=12)
ax.set_ylabel('Normalized Pattern [dB]', fontsize=12)
ax.set_title('Half-wave Slot\nH-plane Pattern', fontsize=13)
ax.set_xlim(0, 180)
ax.set_ylim(-30, 3)
ax.grid(True, alpha=0.3)
ax.axhline(y=-3, color='r', linestyle='--', alpha=0.5, label='-3 dB')
ax.legend()

# (c) 極座標での比較(両方重ねて表示)
ax = plt.subplot(133, projection='polar')
theta_full = np.concatenate([theta, 2*np.pi - theta[::-1]])
F_full = np.concatenate([F_dipole, F_dipole[::-1]])
ax.plot(theta_full, F_full, 'c-', linewidth=2, label='Dipole E / Slot H')

# ダイポールH面 = スロットE面(無指向性)
F_omni = np.ones_like(theta_full)
ax.plot(theta_full, F_omni, 'm--', linewidth=1.5, alpha=0.7, label='Dipole H / Slot E')

ax.set_title('Pattern Comparison\n(Polar)', fontsize=13, pad=20)
ax.legend(loc='lower right', fontsize=9)

plt.tight_layout()
plt.show()

上のグラフから、バビネの原理に基づくダイポールとスロットの関係が明確に確認できます。

  1. ダイポールのE面パターン(左)とスロットのH面パターン(中央)は完全に同一の形状です。これはバビネの原理の直接的な帰結であり、放射パターンの「形」は相補構造で保存されることを示しています。
  2. 極座標表示(右)では、8の字パターン(ダイポールE面 = スロットH面)と無指向性パターン(ダイポールH面 = スロットE面)の対応が見て取れます。偏波面が90°回転する関係が視覚的に理解できます。
  3. メインローブ方向($\theta = 90°$、すなわちアンテナ軸に垂直な方向)で最大放射となるドーナツ型のパターンは、ダイポールもスロットも共通です。

Pythonによるインピーダンスの相補関係の検証

ブッカーの関係式 $Z_s \cdot Z_d = \eta_0^2/4$ を周波数に対して検証します。

import numpy as np
import matplotlib.pyplot as plt

# --- ダイポールのインピーダンス(近似モデル) ---
def dipole_impedance(freq, L=0.15):
    """
    半波長ダイポールのインピーダンスの近似計算。
    共振周波数付近のみ有効。

    Parameters
    ----------
    freq : array  周波数 [Hz]
    L : float     ダイポール長 [m]
    """
    c = 3e8
    lam = c / freq
    kL = 2 * np.pi * L / lam  # 電気的長さ k*L

    # 近似式: 共振時 R ≈ 73 Ohm, X は共振からのずれに比例
    R = 73 * np.ones_like(freq)
    # リアクタンスの近似(共振からのずれに対して線形近似)
    f0 = c / (2 * L)  # 共振周波数
    X = 120 * np.pi * (freq - f0) / f0 * 0.5  # 近似

    return R + 1j * X

# 周波数範囲
L_dipole = 0.15  # 15 cm → f0 = 1 GHz
c = 3e8
f0 = c / (2 * L_dipole)
freq = np.linspace(0.7 * f0, 1.3 * f0, 500)

# ダイポールインピーダンス
Z_d = dipole_impedance(freq, L_dipole)

# ブッカーの関係式によるスロットインピーダンス
eta0 = 377  # Ohm
Z_s = eta0**2 / (4 * Z_d)

fig, axes = plt.subplots(2, 2, figsize=(14, 10))

# (a) ダイポールの抵抗とリアクタンス
ax = axes[0, 0]
ax.plot(freq / 1e9, np.real(Z_d), 'c-', linewidth=2, label='$R_d$')
ax.plot(freq / 1e9, np.imag(Z_d), 'c--', linewidth=2, label='$X_d$')
ax.set_xlabel('Frequency [GHz]', fontsize=12)
ax.set_ylabel('Impedance [$\\Omega$]', fontsize=12)
ax.set_title('Dipole Impedance', fontsize=13)
ax.legend(fontsize=11)
ax.grid(True, alpha=0.3)
ax.axhline(y=0, color='gray', linewidth=0.5)

# (b) スロットの抵抗とリアクタンス
ax = axes[0, 1]
ax.plot(freq / 1e9, np.real(Z_s), 'm-', linewidth=2, label='$R_s$')
ax.plot(freq / 1e9, np.imag(Z_s), 'm--', linewidth=2, label='$X_s$')
ax.set_xlabel('Frequency [GHz]', fontsize=12)
ax.set_ylabel('Impedance [$\\Omega$]', fontsize=12)
ax.set_title('Slot Impedance (from Booker)', fontsize=13)
ax.legend(fontsize=11)
ax.grid(True, alpha=0.3)
ax.axhline(y=0, color='gray', linewidth=0.5)

# (c) Z_s * Z_d の積
ax = axes[1, 0]
product = Z_s * Z_d
ax.plot(freq / 1e9, np.real(product), 'g-', linewidth=2, label='Re$(Z_s Z_d)$')
ax.plot(freq / 1e9, np.imag(product), 'g--', linewidth=2, label='Im$(Z_s Z_d)$')
ax.axhline(y=eta0**2/4, color='r', linestyle=':', linewidth=2,
           label=f'$\\eta_0^2/4 = {eta0**2/4:.0f}\\,\\Omega^2$')
ax.set_xlabel('Frequency [GHz]', fontsize=12)
ax.set_ylabel('$Z_s \\cdot Z_d$ [$\\Omega^2$]', fontsize=12)
ax.set_title('Verification: $Z_s \\cdot Z_d = \\eta_0^2/4$', fontsize=13)
ax.legend(fontsize=10)
ax.grid(True, alpha=0.3)

# (d) VSWRの比較(50 Ohm基準 vs 188 Ohm基準)
ax = axes[1, 1]
Z0_50 = 50
Z0_188 = eta0 / 2

Gamma_d = (Z_d - Z0_50) / (Z_d + Z0_50)
VSWR_d = (1 + np.abs(Gamma_d)) / (1 - np.abs(Gamma_d))

Gamma_s = (Z_s - Z0_188) / (Z_s + Z0_188)
VSWR_s = (1 + np.abs(Gamma_s)) / (1 - np.abs(Gamma_s))

ax.plot(freq / 1e9, VSWR_d, 'c-', linewidth=2, label='Dipole (50 $\\Omega$)')
ax.plot(freq / 1e9, VSWR_s, 'm-', linewidth=2, label='Slot (188 $\\Omega$)')
ax.axhline(y=2, color='r', linestyle='--', alpha=0.5, label='VSWR = 2')
ax.set_xlabel('Frequency [GHz]', fontsize=12)
ax.set_ylabel('VSWR', fontsize=12)
ax.set_title('VSWR Comparison', fontsize=13)
ax.set_ylim(1, 5)
ax.legend(fontsize=11)
ax.grid(True, alpha=0.3)

plt.tight_layout()
plt.show()

この4つのグラフから、ブッカーの関係式の妥当性と実用的な意味が明確に確認できます。

  1. ダイポールのインピーダンス(左上): 共振周波数(1 GHz)で抵抗が約73 $\Omega$、リアクタンスがゼロとなる典型的な共振特性が見られます。共振点から離れるとリアクタンスが急激に増大します。
  2. スロットのインピーダンス(右上): ブッカーの関係式から計算されたスロットのインピーダンスは、共振時に約487 $\Omega$ と高い値を示します。ダイポールとは「逆数的」な関係にあることがわかります。
  3. $Z_s \cdot Z_d$ の検証(左下): 積の実部が全周波数で $\eta_0^2/4 \approx 35525\,\Omega^2$ に一致し、虚部がゼロです。ブッカーの関係式が厳密に成り立つことが確認できます。
  4. VSWR比較(右下): ダイポールを50 $\Omega$ 系に、スロットを188 $\Omega$(自己相補インピーダンス)系に接続した場合、インピーダンス帯域幅の傾向は類似しています。

Pythonによる導波管スロットアレイの放射パターン

最後に、導波管スロットアレイのアレイファクタを計算し、スロット数による指向性の変化を可視化します。

import numpy as np
import matplotlib.pyplot as plt

def array_factor(theta, N, d, wavelength):
    """
    等間隔リニアアレイのアレイファクタを計算する。

    Parameters
    ----------
    theta : array  角度 [rad]
    N : int         素子数
    d : float       素子間隔 [m]
    wavelength : float  波長 [m]
    """
    k = 2 * np.pi / wavelength
    psi = k * d * np.sin(theta)

    # AF = sin(N*psi/2) / sin(psi/2)
    numerator = np.sin(N * psi / 2)
    denominator = np.sin(psi / 2)

    # 0/0 の特異点を処理
    AF = np.where(np.abs(denominator) < 1e-10,
                  N * np.ones_like(theta),
                  numerator / denominator)
    return np.abs(AF) / N  # 正規化

# パラメータ
wavelength = 0.03  # 30 mm (10 GHz)
d = 0.5 * wavelength  # 半波長間隔
theta = np.linspace(-np.pi/2, np.pi/2, 1000)

# 異なるスロット数でのアレイファクタ
N_values = [4, 8, 16, 32]
colors = ['#FF6B6B', '#4ECDC4', '#45B7D1', '#96CEB4']

fig, axes = plt.subplots(1, 2, figsize=(14, 6))

# (a) 直交座標表示
ax = axes[0]
for N, color in zip(N_values, colors):
    AF = array_factor(theta, N, d, wavelength)
    AF_dB = 20 * np.log10(AF + 1e-10)
    AF_dB = np.clip(AF_dB, -40, 0)
    ax.plot(np.degrees(theta), AF_dB, color=color, linewidth=1.5,
            label=f'N = {N} slots')

ax.set_xlabel('Theta [degrees]', fontsize=12)
ax.set_ylabel('Array Factor [dB]', fontsize=12)
ax.set_title('Waveguide Slot Array: Array Factor', fontsize=13)
ax.set_xlim(-90, 90)
ax.set_ylim(-40, 3)
ax.legend(fontsize=11)
ax.grid(True, alpha=0.3)
ax.axhline(y=-3, color='r', linestyle='--', alpha=0.3)

# (b) 極座標表示(N=16の例)
ax = plt.subplot(122, projection='polar')
for N, color in zip(N_values, colors):
    AF = array_factor(theta, N, d, wavelength)
    ax.plot(theta, AF, color=color, linewidth=1.5, label=f'N={N}')

ax.set_title('Array Factor (Polar)', fontsize=13, pad=20)
ax.legend(loc='lower right', fontsize=9)

plt.tight_layout()
plt.show()

導波管スロットアレイの放射パターンから、以下の重要な特性が読み取れます。

  1. スロット数 $N$ が増えるほどメインビームが狭くなることが明確に見て取れます。$N = 4$ では半値幅が約25°ですが、$N = 32$ では約3°まで狭まります。ビーム幅は概ね $\lambda / (Nd)$ に比例して狭くなります。
  2. サイドローブが現れることも確認できます。等振幅励振の場合、第1サイドローブは約-13 dBです。実用的なレーダーアンテナでは、テイラー分布やチェビシェフ分布などの振幅テーパーを適用してサイドローブを-20〜-40 dBに抑制します。
  3. グレーティングローブは発生していないことがわかります。これは素子間隔が $d = \lambda/2$ であるため、可視領域($|\theta| < 90°$)内にグレーティングローブが生じないためです。

実用上の注意点

有限導体板の影響

理論上のスロットアンテナは無限導体板を仮定していますが、実際の導体板は有限サイズです。有限導体板では、板の端で電流が回折し、裏面にも放射が漏れます。その結果、放射パターンにリップルが生じたり、前後比(F/B比)が低下したりします。

一般的な目安として、導体板の各辺はスロット長の2倍以上が望ましいとされています。

スロットの形状バリエーション

基本的なストレートスロットに加えて、以下のような形状バリエーションが実用上使われます。

ボウタイスロット: スロットの中央を広げた形状。帯域幅が広がります。

リングスロット: 環状のスロット。円偏波を発生させることが可能です。

折り返しスロット: スロットをU字型やH字型に折り返した形状。小型化が可能です。

テーパードスロット(ビバルディアンテナ): スロット幅を指数関数的に広げた形状。非常に広い帯域と高い指向性を実現する進行波アンテナです。

まとめ

本記事では、スロットアンテナの理論と設計について解説しました。

  • バビネの原理: 相補的スクリーンの透過波において、電界と磁界の役割が入れ替わります。これにより、スロットアンテナはダイポールアンテナの「相補的な双子」として理解できます
  • 放射パターンの等価関係: スロットのH面パターンはダイポールのE面パターンと同一形状であり、偏波面が90°回転します
  • ブッカーの関係式: 相補的アンテナのインピーダンスの積は $Z_s \cdot Z_d = \eta_0^2/4$ を満たします。半波長スロットのインピーダンスは約487 $\Omega$ です
  • 自己相補構造: $Z = \eta_0/2 \approx 188.5\,\Omega$ の周波数非依存インピーダンスが得られます
  • 導波管スロットアレイ: 導波管壁面のスロットを配列することで、高利得・低損失のアレイアンテナが構成でき、レーダーや衛星通信で広く使われています

バビネの原理は、一見すると抽象的に見えますが、「導体と開口を入れ替えると、電界と磁界も入れ替わる」という明快な物理を表現しています。この原理を理解することで、スロットアンテナだけでなく、パッチアンテナやFSS(周波数選択表面)など、開口を持つ構造全般の放射メカニズムを見通せるようになります。

次のステップとして、以下の記事も参考にしてください。