直径 5 m のパラボラアンテナを 30 GHz で使うとします。このアンテナの放射パターンをきちんと測るには、送信アンテナをどれだけ離せばよいでしょうか。教科書の遠方界条件 $2D^2/\lambda$ に数値を入れると、答えは約 5 km です。5 km の電波暗室は存在しませんし、屋外に 5 km の測定レンジを確保するのも現実的ではありません。では、5 km より近い場所で測ったパターンは何が違うのでしょうか。まったくの別物なのでしょうか、それとも「少しぼけた」程度なのでしょうか。
この問いに答えるのが、アンテナ周辺の空間を距離で 3 つに切り分ける考え方です。アンテナのすぐそばはリアクティブ近傍界、遠く離れた場所はフラウンホーファー領域（遠方界）、そしてその中間に広がるのが本記事の主役、フレネル領域（放射近傍界）です。フレネル領域では電磁波はすでに「放射」として外向きに飛んでいるのに、放射パターンの形が測定距離ごとに変わってしまいます。ビームは太くなり、サイドローブは持ち上がり、ヌル（パターンの落ち込み）は埋まります。
この領域を理解しておくと、次のような場面で判断ができるようになります。
- アンテナ測定: 電波暗室のサイズが足りないとき、どれだけ誤差を許容すればよいか、あるいは近傍界測定＋NF-FF 変換に切り替えるべきかを定量的に判断できます
- レーダー・近距離無線: 至近距離の標的や、車載レーダーのバンパー越しの測定など、被測定物がフレネル領域に入る状況でのビーム形状の変化を予測できます
- 無線電力伝送・ビームフォーカシング: フレネル領域では「焦点」を作れるという性質を積極利用します。近傍界に電力を集中させる WPT や、ミリ波の近距離リンクの設計はこの領域の物理そのものです
- 衛星地球局の据付・校正: 大口径アンテナは実運用でも周囲の建物やレドームがフレネル領域に入るため、境界距離を知らないと校正結果を誤読します
本記事の内容
- アンテナ周辺の 3 領域（リアクティブ近傍界・フレネル領域・フラウンホーファー領域）の直感的な意味
- 開口面上の点から観測点までの距離のテイラー展開（省略なしの導出）
- 二次位相項から遠方界距離 $R = 2D^2/\lambda$ を導く（位相誤差 $\pi/8$ 基準）
- 三次位相項から内側境界 $R = 0.62\sqrt{D^3/\lambda}$ を導く
- フレネル領域で放射パターンがぼける理由 — 二次位相＝デフォーカスしたレンズ
- Python: 線状開口の放射積分を距離を変えて数値計算し、パターンの距離依存を実測
- Python: 第一サイドローブレベル・ヌル深さ・軸上利得損失の距離依存
- 「フレネル領域」と「フレネルゾーン」の用語の違いを図で整理
前提知識
この記事を読む前に、以下の記事を読んでおくと理解が深まります。
アンテナのまわりは 3 つの世界に分かれている
この図が本記事の全体像です。横軸はアンテナからの距離で、境界 $R_1 = 0.62\sqrt{D^3/\lambda}$ と $R_2 = 2D^2/\lambda$ で 3 つの帯に分かれています。各帯の中に描いた曲線は、実際にその距離で放射積分を解いて得た「そこで見えるパターン」で、いちばん近い帯（距離 $0.02R_2$）ではもはやビームの形すらなく開口の幅いっぱいに波打つ分布、中間の帯（$0.25R_2$）ではメインローブは見えるがヌルが埋まってサイドローブと地続きになった形、いちばん外の帯（$10R_2$）では教科書どおりの鋭いメインローブと深いヌルが並ぶ形になっています。距離が変わるとパターンの「形そのもの」が変わるというのが左と中央の帯の特徴で、それが消えて形が固定するのが右の帯だと押さえておいてください。
夜道で懐中電灯を点けたとき、レンズのすぐ前に手をかざすと、そこに見えるのは「ビーム」ではなくレンズの形そのものです。少し離すと、まだレンズの形を引きずった明るい四角が見えます。さらに離してようやく、私たちが「懐中電灯のビーム」と呼ぶ、周辺がなだらかに暗くなる円錐状の光が現れます。ここで重要なのは、離すにつれて明るい部分の「形」が変わっていくという点です。ある距離を超えると形は変わらなくなり、あとは離れた分だけ全体が大きく・暗くなるだけになります。
アンテナのまわりでもまったく同じことが起きています。電波は目に見えないので実感しにくいだけです。アンテナの開口（電波を放出している面）の大きさを $D$、波長を $\lambda$ とすると、距離 $R$ に応じて 3 つの世界が現れます。
リアクティブ近傍界（reactive near field） はアンテナのごく近くです。ここでは電磁エネルギーの大半が「放射」されずにアンテナのまわりに蓄えられ、半周期ごとにアンテナと空間の間を行き来しています。コンデンサやコイルに蓄えられるエネルギーと同じ性質のもので、正味の電力は外へ出ていきません。電界と磁界の位相がずれており、両者の比（波動インピーダンス）は自由空間の $377\,\Omega$ から大きく外れます。この領域に金属や誘電体を近づけると、アンテナの入力インピーダンス自体が変わってしまいます。スマートフォンを手で握るとアンテナ特性が変わるのは、手がリアクティブ近傍界に入るからです。
フレネル領域（Fresnel region）、別名 放射近傍界（radiating near field） はその外側です。ここではエネルギーはすでに外向きに流れており、電力は確かに放射されています。ところが、開口の各部分から届く波の位相が観測点でうまく揃わないため、パターンの形が距離とともに変わります。近いほどビームは太く、サイドローブとメインローブの区別は曖昧で、ヌルは埋まっています。この領域の解析には、球面波の距離を二次の項まで近似する「フレネル近似」が使えます。名前の由来はここにあります。
フラウンホーファー領域（Fraunhofer region）、すなわち 遠方界（far field） はさらに外側です。ここまで来ると、開口の各点から観測点へ向かう線がほぼ平行とみなせるようになり、位相差は「開口面上の位置 × 観測方向」という単純な一次式だけで決まります。その結果、放射パターンは開口面上の電磁界分布のフーリエ変換という美しい形に収束し、以後どれだけ離れても形は変わりません。変わるのは全体の大きさが $1/R$ で小さくなることだけです。私たちが「アンテナの放射パターン」と呼ぶものは、この領域でのパターンを指しています。
3 つの領域を分ける境界は、次の 2 本の式で与えられます。
$$ \begin{equation} R_1 = 0.62\sqrt{\frac{D^3}{\lambda}}, \qquad R_2 = \frac{2D^2}{\lambda} \end{equation} $$
$R < R_1$ がリアクティブ近傍界、$R_1 < R < R_2$ がフレネル領域、$R > R_2$ がフラウンホーファー領域です。この $0.62$ という半端な数字と、$2D^2/\lambda$ の $2$ という係数はどこから来たのでしょうか。どちらも「開口の端から届く波の位相が、中心から届く波に対してどれだけずれるか」という、たったひとつの基準から導かれます。次の節でその出発点となる放射積分を立て、距離をテイラー展開してみましょう。
放射積分と、距離のテイラー展開
舞台設定
図に描いたとおり、議論の全体はたった 2 本の線分の長さの差に集約されます。青い線が開口中心から観測点までの距離 $R$、赤い線が開口上の点 $x’$ から観測点までの距離 $r$ で、この $r – R$ がどれだけ $x’$ に依存するかがすべてを決めます。以下では厳密な $r = \sqrt{R^2 – 2Rx’\sin\theta + x’^2}$ を $x’/R$ で展開し、1 次項が遠方界パターンを作り、2 次項がフレネル領域の主犯となり、3 次項がフレネル近似の限界を決める、という順に読み解いていきます。
議論を見通しよくするため、まずは 1 次元の線状開口を考えます。長さ $D$ の開口が $x’$ 軸上に $-D/2 \le x’ \le D/2$ の範囲で置かれ、開口面は $z = 0$ の平面にあるとします。開口上の位置 $x’$ における電磁界の複素振幅を $A(x’)$ と書きます（一様励振なら $A(x’) = 1$ です）。
観測点 $P$ は、開口の中心から距離 $R$、法線（$z$ 軸）から測った角度 $\theta$ の方向にあるとします。デカルト座標で書けば
$$ P = (R\sin\theta,\; 0,\; R\cos\theta), \qquad \text{開口上の点} = (x’,\; 0,\; 0) $$
です。開口面上の微小要素からは球面波が出ており、観測点での寄与は $A(x’)\,e^{-jkr}/r$ に比例します。ここで $k = 2\pi/\lambda$ は波数、$r$ は開口上の点 $x’$ から観測点 $P$ までの距離です。全体の界はこれを開口全体で足し合わせた
$$ \begin{equation} E(\theta; R) = \int_{-D/2}^{D/2} A(x’) \frac{e^{-jkr(x’,\theta;R)}}{r(x’,\theta;R)}\, dx’ \end{equation} $$
で与えられます。これが放射積分です。近傍界も遠方界も、すべてこの一本の式に含まれています。領域を分けるのは物理が変わるからではなく、$r$ をどこまで近似できるかが変わるからにすぎません。
距離 $r$ を厳密に書く
まず $r$ を厳密に求めます。2 点間の距離なので
$$ r = \sqrt{(R\sin\theta – x’)^2 + (R\cos\theta)^2} $$
括弧を展開すると、$R^2\sin^2\theta – 2Rx’\sin\theta + x’^2 + R^2\cos^2\theta$ となり、$\sin^2\theta + \cos^2\theta = 1$ を使えば $R^2$ にまとまります。したがって
$$ \begin{equation} r = \sqrt{R^2 – 2Rx’\sin\theta + x’^2} \end{equation} $$
これは厳密な式で、近似は一切入っていません。この平方根が扱いにくいので、$x’/R$ が小さいことを利用して展開します。
二項展開を丁寧にたどる
$R$ を括り出して
$$ r = R\sqrt{1 + u}, \qquad u \equiv \frac{x’^2 – 2Rx’\sin\theta}{R^2} $$
と置きます。$u$ は $x’/R$ について 1 次と 2 次の項からなる小さい量です。$\sqrt{1+u}$ の二項展開（マクローリン展開）は
$$ \sqrt{1+u} = 1 + \frac{u}{2} – \frac{u^2}{8} + \frac{u^3}{16} – \cdots $$
です。ここが導出の要になります。$u$ 自体が $x’/R$ の 1 次項と 2 次項を含んでいるため、$u$ の 1 次・2 次・3 次のそれぞれが $x’$ の複数の次数に寄与します。$x’$ の次数ごとに項を集め直す必要があるのです。順に見ていきましょう。
まず $u/2$ の項に $R$ を掛けます。
$$ R \cdot \frac{u}{2} = \frac{x’^2 – 2Rx’\sin\theta}{2R} = -x’\sin\theta + \frac{x’^2}{2R} $$
$x’$ の 1 次項 $-x’\sin\theta$ と、2 次項 $x’^2/(2R)$ が出てきました。
次に $-u^2/8$ の項です。$u$ を 2 乗すると
$$ u^2 = \frac{(x’^2 – 2Rx’\sin\theta)^2}{R^4} = \frac{x’^4 – 4Rx’^3\sin\theta + 4R^2x’^2\sin^2\theta}{R^4} $$
これに $-R/8$ を掛けると
$$ -R\cdot\frac{u^2}{8} = -\frac{x’^4}{8R^3} + \frac{x’^3\sin\theta}{2R^2} – \frac{x’^2\sin^2\theta}{2R} $$
$x’$ の 2 次項 $-x’^2\sin^2\theta/(2R)$ と 3 次項 $x’^3\sin\theta/(2R^2)$ が現れました。ここで、先ほど $u/2$ から出た 2 次項 $x’^2/(2R)$ と合わせると
$$ \frac{x’^2}{2R} – \frac{x’^2\sin^2\theta}{2R} = \frac{x’^2(1 – \sin^2\theta)}{2R} = \frac{x’^2\cos^2\theta}{2R} $$
と、きれいに $\cos^2\theta$ にまとまります。三角関数の恒等式 $1 – \sin^2\theta = \cos^2\theta$ を使っただけですが、この $\cos^2\theta$ が後で「軸上（$\theta = 0$）で位相誤差が最大になる」という重要な結論を生みます。
最後に $u^3/16$ の項です。$x’$ の 3 次項だけを拾えばよいので、$u$ の 3 乗のうち $(-2Rx’\sin\theta)^3 = -8R^3x’^3\sin^3\theta$ の項だけを取り出します。
$$ R\cdot\frac{u^3}{16} \supset R \cdot \frac{-8R^3x’^3\sin^3\theta}{16R^6} = -\frac{x’^3\sin^3\theta}{2R^2} $$
先ほど $-u^2/8$ から出た 3 次項 $x’^3\sin\theta/(2R^2)$ と足し合わせると
$$ \frac{x’^3\sin\theta}{2R^2} – \frac{x’^3\sin^3\theta}{2R^2} = \frac{x’^3\sin\theta(1 – \sin^2\theta)}{2R^2} = \frac{x’^3\sin\theta\cos^2\theta}{2R^2} $$
となります。ここでも $1 – \sin^2\theta = \cos^2\theta$ が効いています。
展開の結論とその読み方
以上をまとめると、距離 $r$ は次のように展開されます。
$$ \begin{equation} r \simeq \underbrace{R}_{\text{共通位相}} \underbrace{- \, x’\sin\theta}_{\text{一次項}} + \underbrace{\frac{x’^2\cos^2\theta}{2R}}_{\text{二次項}} + \underbrace{\frac{x’^3\sin\theta\cos^2\theta}{2R^2}}_{\text{三次項}} + \cdots \end{equation} $$
この 4 つの項が、そのまま 3 領域の理論に対応しています。ひとつずつ意味を読み取る前に、この展開が本当に正しいかを厳密な距離から確かめておきましょう。
左は軸上（$\theta = 0$）で、厳密な $r$ から共通位相 $R$ と一次項 $-x’\sin\theta$ を引いた残りをプロットしたものです。どの距離でも $x’^2$ の放物線になり、$R = R_2$ の赤線は二次項の予測（黒破線）とぴったり重なって、開口端でちょうど $22.5^\circ$ に達しています。右は三次項が最大になる $\theta = 35.26^\circ$ で、さらに二次項も引いた残りを描いたもので、こちらは $x’^3$ の奇関数（原点対称）になっており、開口の左右で符号が反転する — つまりパターンを非対称に歪ませる誤差であることが目で確認できます。$R = R_2$（赤）ではこの三次項の残りはほぼゼロで、$0.05R_2$ まで近づいて初めて $22.5^\circ$ を超えます。
では 4 つの項の意味をひとつずつ読み取りましょう。
第 1 項 $R$ は $x’$ に依存しない共通の位相です。積分の外に $e^{-jkR}$ として括り出せるので、パターンの形には一切影響しません。全体の遅れと $1/R$ の減衰を生むだけです。
第 2 項 $-x’\sin\theta$ が遠方界パターンを作る主役です。この項だけを残して $r \simeq R – x’\sin\theta$（分母の $r$ は $R$ で近似）とすると、放射積分は
$$ E(\theta) \simeq \frac{e^{-jkR}}{R}\int_{-D/2}^{D/2} A(x’)\, e^{jkx’\sin\theta}\, dx’ $$
となります。これは開口分布 $A(x’)$ のフーリエ変換そのものです。$R$ が積分の中から完全に消えている点に注目してください。パターンの形が距離に依存しなくなるのは、まさにこの状況を指します。
第 3 項 $x’^2\cos^2\theta/(2R)$ が本記事の主役、二次位相項です。$1/R$ に比例するので、遠ざかれば消えますが近づくと効いてきます。この項が残っている状態が「フレネル領域」で、この項を残したまま行う近似をフレネル近似と呼びます。$x’^2$ に比例するということは、開口の中心から離れた場所ほど余計な位相遅れが加わるということです。これは球面波の波面の曲率そのもので、後で見るように「デフォーカスしたレンズ」と同じ効果を持ちます。
第 4 項 $x’^3\sin\theta\cos^2\theta/(2R^2)$ は三次位相項です。$1/R^2$ に比例するので、二次項よりずっと速く消えます。逆にいえば、この項が無視できなくなるほど近づくと、フレネル近似そのものが壊れます。$\sin\theta$ を含むので $\theta = 0$ ではゼロになり、$x’$ の奇数乗なので開口の左右で符号が反転する — つまりパターンを非対称に歪ませるタイプの誤差です。
ここまでで、3 つの領域を分ける物差しが「開口の端から届く波の位相がどれだけ余計にずれるか」であることが見えました。次は、この位相誤差にどれだけの上限を課すかを決め、境界距離を具体的に計算します。
遠方界距離 $2D^2/\lambda$ の導出
位相誤差 $\pi/8$ という基準
放射パターンを「遠方界とみなしてよい」と判断するには、どこかで線を引く必要があります。慣習的に使われる基準は
$$ \text{最大位相誤差} \le \frac{\pi}{8}\ \text{rad} = 22.5^\circ $$
です。この $\pi/8$ は経路長に直すと
$$ \Delta = \frac{\Delta\phi}{k} = \frac{\pi/8}{2\pi/\lambda} = \frac{\lambda}{16} $$
つまり $\lambda/16$ の経路差に相当します。この経路差が何と何の差なのかを図で押さえておきましょう。
青い破線が「遠方界ならこうなっているはず」という平面の等位相面、赤い曲線が観測点を中心とする本当の球面の等位相面です（縦方向は見やすさのために誇張してあります）。両者は開口中心で接し、開口端に向かうほど離れていきます。このずれが開口端からの波の余分な遠回り、すなわち経路差 $\Delta \simeq (D/2)^2/(2R) = D^2/(8R)$ であり、これを $\lambda/16$ 以下に抑えよという条件を解くだけで $R \ge 2D^2/\lambda$ が出てきます。以降の代数計算は、この 1 枚の絵を式に翻訳しているにすぎません。「開口の端から来る波の道のりが、中心から来る波より $\lambda/16$ だけ長い程度なら、まあ許そう」という取り決めです。$\lambda/16$ という値自体に深い理論的必然性はありませんが、後で数値実験で確かめるように、この基準を満たせば軸上利得の低下は $0.1$ dB を下回り、第一サイドローブレベルの誤差も $0.3$ dB 程度に収まります。実用的にちょうどよい線引きとして 20 世紀半ばから使われ続けている、経験に裏打ちされた基準です。
二次項の最大値を求める
二次位相項による位相誤差は、展開の第 3 項に波数 $k$ を掛けて
$$ \Delta\phi_2(x’, \theta) = k\,\frac{x’^2\cos^2\theta}{2R} $$
です。これを開口全体と全観測角にわたって最大化します。$x’^2$ は $x’ = \pm D/2$ で最大、$\cos^2\theta$ は $\theta = 0$（軸上）で最大 $1$ です。したがって
$$ \Delta\phi_{2,\max} = k\,\frac{(D/2)^2}{2R} = \frac{2\pi}{\lambda}\cdot\frac{D^2/4}{2R} = \frac{\pi D^2}{4\lambda R} $$
ここで $k = 2\pi/\lambda$ を代入し、$(D/2)^2 = D^2/4$ を使って整理しました。分子分母の $2$ を約分して $\pi D^2/(4\lambda R)$ という簡潔な形になっています。
境界距離を解く
この最大位相誤差が $\pi/8$ 以下という条件を課します。
$$ \frac{\pi D^2}{4\lambda R} \le \frac{\pi}{8} $$
両辺を $\pi$ で割ると $D^2/(4\lambda R) \le 1/8$ となります。両辺に $8 \cdot 4\lambda R = 32\lambda R$ を掛けると $8D^2 \le 4\lambda R$、つまり $R \ge 2D^2/\lambda$ が得られます。まとめると
$$ \begin{equation} \boxed{\;R \ge R_2 = \frac{2D^2}{\lambda}\;} \end{equation} $$
これがフラウンホーファー距離あるいは遠方界距離です。$2D^2/\lambda$ の「$2$」は、$\pi/8$ 基準からきっちり出てくる係数だったわけです。もし基準を $\pi/16$（$\lambda/32$）に厳しくすれば $R \ge 4D^2/\lambda$、逆に $\pi/4$ に緩めれば $R \ge D^2/\lambda$ となります。文献によって $D^2/\lambda$ や $4D^2/\lambda$ という値が出てくるのは、この基準の取り方が違うだけです。
数値で実感する
いくつかの実例で $R_2$ を計算してみましょう。
| アンテナ | 口径 $D$ | 周波数 | 波長 $\lambda$ | $R_1 = 0.62\sqrt{D^3/\lambda}$ | $R_2 = 2D^2/\lambda$ |
|---|---|---|---|---|---|
| Wi-Fi 端末アンテナ | 0.3 m | 2.4 GHz | 125 mm | 0.29 m | 1.4 m |
| 衛星放送用パラボラ | 0.6 m | 12 GHz | 25.0 mm | 1.82 m | 28.8 m |
| C バンド地球局 | 1.2 m | 4 GHz | 75.0 mm | 2.98 m | 38.4 m |
| X バンド地球局 | 3 m | 10 GHz | 30.0 mm | 18.6 m | 600 m |
| Ka バンド大口径 | 5 m | 30 GHz | 10.0 mm | 69.3 m | 5.0 km |
表を眺めると、$R_2$ が $D^2/\lambda$ に比例するため、口径と周波数が上がると爆発的に遠くなることがわかります。$0.6$ m のパラボラなら $29$ m の屋外レンジで測れますが、$5$ m の Ka バンドアンテナは $5$ km — 事実上、遠方界での直接測定は不可能です。この壁があるからこそ、近傍界で測って数学的に遠方界へ変換する NF-FF 変換や、放物面鏡で平面波を作るコンパクトレンジといった技術が必要になります。詳しくはアンテナ測定の基礎を参照してください。
外側の境界が決まりました。では内側、つまり「エネルギーが確かに放射として出ていく」と言えるようになる距離はどこでしょうか。次はもう一段深く、三次項に踏み込みます。
内側境界 $0.62\sqrt{D^3/\lambda}$ の導出
なぜ三次項なのか
フレネル領域とは「フレネル近似（二次項まで残した近似）が使える領域」です。逆にいえば、三次項が無視できなくなったらフレネル近似は破綻し、そこから内側はもはやフレネル領域と呼べません。この論理で内側境界を定めるのが、$0.62\sqrt{D^3/\lambda}$ という式の出自です。
同じ $\pi/8$ 基準を三次項に適用します。三次位相項は
$$ \Delta\phi_3(x’,\theta) = k\,\frac{x’^3\sin\theta\cos^2\theta}{2R^2} $$
でした。二次項と違い、こちらは $\theta$ について $\sin\theta\cos^2\theta$ という形をしています。$\theta = 0$ でも $\theta = 90^\circ$ でもゼロになり、途中のどこかで最大値をとります。
角度依存部の最大値
$g(\theta) = \sin\theta\cos^2\theta$ を最大化しましょう。$\sin\theta = s$ と置くと $g = s(1 – s^2) = s – s^3$ です。$s$ で微分してゼロと置くと
$$ \frac{dg}{ds} = 1 – 3s^2 = 0 \quad \Longrightarrow \quad s = \frac{1}{\sqrt{3}} $$
このとき $\theta = \arcsin(1/\sqrt{3}) \approx 35.26^\circ$ です。最大値は $s$ を代入して
$$ g_{\max} = \frac{1}{\sqrt{3}}\left(1 – \frac{1}{3}\right) = \frac{1}{\sqrt{3}}\cdot\frac{2}{3} = \frac{2}{3\sqrt{3}} \approx 0.3849 $$
$1/\sqrt{3}$ を代入して $1 – s^2 = 1 – 1/3 = 2/3$ を使っただけの計算です。この $2/(3\sqrt{3})$ という数が、後で $0.62$ という係数を生みます。
緑の実線が三次項の角度依存 $g(\theta) = \sin\theta\cos^2\theta$ で、$\theta = 0$ と $\theta = 90^\circ$ でゼロ、その中間の $35.26^\circ$ で最大値 $0.3849$ をとる山型になっています。橙の破線は比較のために描いた二次項の角度依存 $\cos^2\theta$ で、こちらは $\theta = 0$（軸上）で最大 $1$ です。二次項の最悪角度は軸上、三次項の最悪角度は斜め $35.26^\circ$ という違いがこの図の要点で、同じ $\pi/8$ 基準を適用しても現れる係数が異なる（$2$ と $0.62$）のは、この角度依存の違いが直接の原因です。
境界距離を解く
$x’ = D/2$、$g = g_{\max}$ を代入すると、三次位相項の最大値は
$$ \Delta\phi_{3,\max} = \frac{2\pi}{\lambda}\cdot\frac{(D/2)^3}{2R^2}\cdot\frac{2}{3\sqrt{3}} $$
$(D/2)^3 = D^3/8$ を代入し、係数を整理します。$2\pi/(\lambda \cdot 8 \cdot 2R^2) = \pi/(8\lambda R^2)$ なので
$$ \Delta\phi_{3,\max} = \frac{\pi D^3}{8\lambda R^2}\cdot\frac{2}{3\sqrt{3}} $$
これを $\pi/8$ 以下と置きます。
$$ \frac{\pi D^3}{8\lambda R^2}\cdot\frac{2}{3\sqrt{3}} \le \frac{\pi}{8} $$
両辺の $\pi/8$ が丸ごと約分されるのが気持ちのよいところです。残るのは
$$ \frac{D^3}{\lambda R^2}\cdot\frac{2}{3\sqrt{3}} \le 1 $$
$R^2$ について解くと
$$ R^2 \ge \frac{2}{3\sqrt{3}}\cdot\frac{D^3}{\lambda} $$
平方根をとって
$$ \begin{equation} \boxed{\;R \ge R_1 = \sqrt{\frac{2}{3\sqrt{3}}}\sqrt{\frac{D^3}{\lambda}} = 0.6204\sqrt{\frac{D^3}{\lambda}}\;} \end{equation} $$
$\sqrt{2/(3\sqrt{3})} = \sqrt{0.38490} = 0.62040$ です。教科書に載っている $0.62$ という半端な係数は、こうして $\sqrt{2/(3\sqrt{3})}$ という閉じた形から出てきます。魔法の数字でも実験式でもなく、$\pi/8$ 基準と三次項の角度依存性だけで決まる純粋に幾何学的な定数です。
二つの境界のスケーリングの違い
$R_1 \propto D^{3/2}\lambda^{-1/2}$、$R_2 \propto D^2\lambda^{-1}$ という違いに注目してください。両者の比は
$$ \frac{R_2}{R_1} = \frac{2D^2/\lambda}{0.62\sqrt{D^3/\lambda}} = \frac{2}{0.62}\sqrt{\frac{D}{\lambda}} \approx 3.23\sqrt{\frac{D}{\lambda}} $$
左の両対数グラフでは、緑の $R_1$ が傾き $3/2$、赤の $R_2$ が傾き $2$ の直線になっており、開口を大きくするほど 2 本が離れていく — つまり橙で塗ったフレネル領域が縦に広がっていくことがわかります。右はその比 $R_2/R_1 = 3.23\sqrt{D/\lambda}$ で、$D = 10\lambda$ では約 $10$ 倍にすぎませんが、$D = 500\lambda$（5 m・30 GHz）では約 $72$ 倍に達します。$D/\lambda$、つまり波長で測った開口の大きさが大きいほど、フレネル領域は相対的に広くなります。$D = 10\lambda$ なら $R_2/R_1 \approx 10$ で、フレネル領域は 1 桁分の距離範囲を占めます。$D = 500\lambda$（Ka バンドの 5 m アンテナ）なら $R_2/R_1 \approx 72$ で、フレネル領域は $69$ m から $5$ km まで、実に 2 桁近い広がりを持ちます。大口径・高周波のアンテナほど、実務で触れる空間のほとんどがフレネル領域なのです。
なお、$D \lesssim \lambda$ の小型アンテナではこの式は使えません。$R_1$ が $\lambda$ より小さくなってしまい、物理的な意味を失うからです。微小ダイポールのような電気的に小さいアンテナでは、リアクティブ近傍界の境界として $R \approx \lambda/(2\pi)$（放射球の半径）を使うのが一般的です。$D^3$ という高次のスケーリングは、あくまで「開口が波長より十分大きい」という前提のもとで成り立ちます。
境界の式が両方揃いました。しかし式が出ただけでは「フレネル領域でパターンがどう崩れるか」はわかりません。次に、二次位相項がパターンに与える具体的な影響を、光学のアナロジーを使って理解しましょう。
フレネル領域でパターンがぼける理由
二次位相＝デフォーカスしたレンズ
フレネル近似での放射積分は、二次項まで残して
$$ E(\theta; R) \simeq \frac{e^{-jkR}}{R}\int_{-D/2}^{D/2} A(x’)\, \exp\!\left(jkx’\sin\theta\right)\exp\!\left(-jk\frac{x’^2\cos^2\theta}{2R}\right) dx’ $$
と書けます。この式の構造を読み解きましょう。第 1 の指数がフーリエ変換の核、第 2 の指数が余計な二次位相です。ここで発想を変えて、第 2 の指数を開口分布のほうに押し込んでしまいます。
$$ A_{\text{eff}}(x’) = A(x’)\exp\!\left(-jk\frac{x’^2\cos^2\theta}{2R}\right) $$
とすれば、$E(\theta;R)$ は「$A_{\text{eff}}$ のフーリエ変換」の形になります（$\theta$ 依存が残るので厳密なフーリエ変換ではありませんが、主ビーム付近では $\cos^2\theta \simeq 1$ とみなせるので実質的にそう扱えます）。つまり、近い距離で測るということは、開口面に二次位相分布を人工的に加えたアンテナを、遠方界で測ることと等価なのです。
開口に $\exp(-jk x’^2/(2R))$ という位相を掛ける操作は、光学ではまさにレンズの作用です。焦点距離 $f$ の薄いレンズの透過位相は $\exp(+jkx’^2/(2f))$ ですから、符号を合わせれば、測定距離 $R$ で測ることは焦点距離 $-R$ のレンズを開口に貼り付けたのと同じ、つまりデフォーカスした状態に相当します。カメラのピントを外したときに点光源がぼやけた円盤になるのと、フレネル領域でビームが太くぼけるのは、まったく同じ物理です。
左は、測定距離ごとに開口へ上乗せされる位相を実際にプロットしたものです。どれも $x’^2$ の放物線で、$R = 10R_2$ では開口端でわずか $2.25^\circ$ しか曲がっていないのに、$R = 0.1R_2$ では $225^\circ$ — 半波長以上 — も曲がってしまいます。右がその物理的な意味で、遠方界では波面が平面（位相が揃う）なのに対し、フレネル領域では波面が球面になり開口端ほど遅れます。開口に $\exp(-jkx’^2/2R)$ を掛けるこの操作は光学でいう焦点距離 $-R$ のレンズを貼る操作そのもので、近づいて測るとはデフォーカスして撮るということだと理解できます。
この見方はいくつもの帰結を一気に説明します。
- ビームが太くなる: デフォーカスすると点像が広がる
- サイドローブが持ち上がる: 位相誤差はメインローブのエネルギーを周辺へ撒き散らす
- ヌルが埋まる: 完全な打ち消しには全開口の位相が正確に揃う必要があり、位相誤差があると打ち消し残りが出る
- 軸上利得が下がる: 開口効率が位相誤差の分だけ落ちる
- パターンが距離ごとに違う: $R$ が変われば「レンズの度」が変わる
コルヌ螺旋との関係
一様励振（$A(x’) = 1$）で軸上（$\theta = 0$）を見ると、積分は
$$ E(0;R) \propto \int_{-D/2}^{D/2}\exp\!\left(-jk\frac{x’^2}{2R}\right)dx’ $$
となります。$x’ = \sqrt{\lambda R}\, t/\sqrt{2}\cdot\sqrt{2}$ のような変数変換で無次元化すると、これはフレネル積分
$$ C(w) = \int_0^{w}\cos\frac{\pi t^2}{2}dt, \qquad S(w) = \int_0^{w}\sin\frac{\pi t^2}{2}dt $$
に帰着します。複素平面上で $C(w) + jS(w)$ を描いた曲線が有名なコルヌ螺旋（クロソイド）です。遠方界極限では積分区間が螺旋の中心付近の直線部分に対応し、寄与がまっすぐ足し合わさって最大値を与えます。近づくにつれて積分区間が螺旋を巻き始め、ベクトル和が「巻いた分だけ短くなる」ため、界の振幅が落ちていきます。ナイフエッジ回折の解析で登場するコルヌ螺旋とまったく同じ道具が、開口アンテナの近傍界解析にも現れるわけです（フレネルゾーンと回折損失を解説 参照）。
左の図で、各色の太線がその測定距離に対応する積分区間、破線がその両端を結ぶ弦で、軸上の界の大きさはこの弦の長さに比例します。$R = 10R_2$（青）では区間が螺旋の中心付近のほぼ直線部分に収まるので弦は区間長そのもの（正規化値 $1.000$）ですが、$R = 0.05R_2$（茶）では区間が螺旋の両端の渦に巻き込まれてしまい、弦の長さは $0.361$ まで縮んでしまいます。右はこれを距離の関数として dB で描いたもので、$R_2$ での低下量は $-0.060$ dB とほぼ無視できる一方、$R_1$ より内側では $-7$ dB 以下まで落ち込み、しかも振動を始めます — この振動が「近づくほど単調に弱くなるとは限らない」という近傍界の厄介さです。
「フレネル領域」と「フレネルゾーン」は別物
用語の整理をしておきます。日本語でも英語でも紛らわしい 2 つの概念があります。
フレネル領域（Fresnel region） は本記事の主題で、1 つのアンテナのまわりの距離 $R_1 < R < R_2$ の空間を指します。アンテナの口径 $D$ と波長 $\lambda$ だけで決まり、送受信の関係は関係ありません。
フレネルゾーン（Fresnel zone） は、送信点と受信点を結ぶ伝搬路のまわりに広がる回転楕円体の領域です。送受信点からの経路差が $n\lambda/2$ となる点の集合として定義され、第 1 フレネルゾーンの半径は、送信点から $d_1$、受信点から $d_2$ の位置で
$$ r_1 = \sqrt{\frac{\lambda d_1 d_2}{d_1 + d_2}} $$
で与えられます。見通し内リンクで「第 1 フレネルゾーンの 60 % を確保せよ」という設計ルールが使われるのはこちらです。障害物の回避に関する概念であって、アンテナの近傍界とは無関係です。
共通しているのは「二次の位相差（経路差）を考える」という数学的な骨格だけで、対象も使いどころもまったく違います。どちらもフレネルの名を冠しているのは、彼が波動光学で導入した二次近似が両方の根底にあるからです。
理屈は揃いました。ここからは Python で実際に放射積分を数値計算し、フレネル領域でパターンがどう崩れるかを自分の目で確かめましょう。
Python で 3 領域を可視化する
距離展開の各項がどれだけ効くか
まず、テイラー展開の各項が距離によってどう変わるかを見ます。$D = 10\lambda$ の線状開口を例に、開口端（$x’ = D/2$）における二次項・三次項の位相誤差を、それぞれが最大になる角度で評価します。
import numpy as np
import matplotlib, matplotlib.pyplot as plt
for cand in ["Hiragino Sans", "Yu Gothic", "Noto Sans CJK JP", "IPAexGothic", "Meiryo"]:
if any(cand == f.name for f in matplotlib.font_manager.fontManager.ttflist):
plt.rcParams["font.family"] = cand
break
plt.rcParams["axes.unicode_minus"] = False
lam = 1.0 # 波長を1に規格化
D = 10 * lam # 開口長
k = 2 * np.pi / lam
R1 = 0.62 * np.sqrt(D**3 / lam) # 内側境界
R2 = 2 * D**2 / lam # 遠方界距離
print(f"R1 = {R1/lam:.1f} λ, R2 = {R2/lam:.1f} λ, R2/R1 = {R2/R1:.2f}")
R = np.logspace(np.log10(2*lam), np.log10(30*R2), 400)
phi2 = k * (D/2)**2 / (2*R) # 二次項（θ=0で最大）
phi3 = k * (D/2)**3 * (2/(3*np.sqrt(3))) / (2*R**2) # 三次項（θ≈35.26°で最大）
plt.figure(figsize=(9, 5.5))
plt.loglog(R/lam, np.degrees(phi2), lw=2, label="二次位相項の最大値 $\\Delta\\phi_2$")
plt.loglog(R/lam, np.degrees(phi3), lw=2, label="三次位相項の最大値 $\\Delta\\phi_3$")
plt.axhline(22.5, color="k", ls="--", lw=1.2, label="許容基準 $\\pi/8$ = 22.5°")
plt.axvline(R1/lam, color="tab:green", ls=":", lw=2, label="$R_1 = 0.62\\sqrt{D^3/\\lambda}$")
plt.axvline(R2/lam, color="tab:red", ls=":", lw=2, label="$R_2 = 2D^2/\\lambda$")
plt.xlabel("距離 $R$ [波長]"); plt.ylabel("最大位相誤差 [度]")
plt.title("開口端の位相誤差と3領域の境界（$D = 10\\lambda$）")
plt.grid(True, which="both", alpha=0.3); plt.legend(fontsize=9)
plt.tight_layout(); plt.show()
このグラフは境界式の意味をそのまま絵にしたものです。二次項の曲線が $22.5^\circ$ の水平線を横切る点がちょうど $R_2 = 200\lambda$、三次項の曲線が横切る点がちょうど $R_1 \approx 19.6\lambda$ になっており、導出が正しいことが確認できます。両対数プロットで二次項は傾き $-1$、三次項は傾き $-2$ の直線になっている点にも注目してください。これは $\Delta\phi_2 \propto 1/R$、$\Delta\phi_3 \propto 1/R^2$ というスケーリングの直接の表れで、遠ざかるにつれて三次項のほうが 2 倍の速さで消えていくことを示しています。だからこそ、外側の境界は二次項が、内側の境界は三次項が支配するのです。
放射積分を厳密に解く関数
次に、近似を一切使わない厳密な距離 $r = \sqrt{R^2 – 2Rx’\sin\theta + x’^2}$ を用いて、放射積分を数値的に評価する関数を用意します。これが本記事の実験の土台になります。
import numpy as np
lam = 1.0
D = 10 * lam
k = 2 * np.pi / lam
R2 = 2 * D**2 / lam
def aperture_field(R, theta, N=40001, taper=None):
"""線状開口の放射積分を厳密な距離で数値評価する。
R : 観測距離
theta : 観測角(rad)のスカラーまたは配列
taper : 開口分布 A(x') を返す関数（Noneなら一様励振）
"""
xp = np.linspace(-D/2, D/2, N)
A = np.ones_like(xp) if taper is None else taper(xp)
theta = np.atleast_1d(np.asarray(theta, dtype=float))
out = np.empty(theta.shape, dtype=complex)
for i, t in enumerate(theta):
# 厳密な距離（近似なし）
r = np.sqrt(R**2 - 2*R*xp*np.sin(t) + xp**2)
out[i] = np.trapz(A * np.exp(-1j*k*r) / r, xp)
return out
# 動作確認：軸上の値が 1/R でスケールすることを見る
for f in [1.0, 10.0, 100.0]:
R = f * R2
print(f"R = {f:6.1f} R2 : R*|E(0)| = {R*abs(aperture_field(R, 0.0))[0]:.5f}")
出力される R*|E(0)| は、$R = R_2$ で $9.9306$、$10R_2$ で $9.9993$、$100R_2$ で $9.99999$ となります。開口長 $D = 10\lambda$ の一様開口では、完全な遠方界で $R|E| \to D = 10$ に収束するはずで、実際にそうなっています。$R = R_2$ の時点ですでに理論値の $99.3\,\%$（$-0.06$ dB）に達しており、$\pi/8$ 基準が「軸上利得はほぼ正しく測れる」水準の線引きであることがわかります。一方で、後で見るようにサイドローブやヌルはこの距離ではまだ十分に収束していません。
距離を変えて放射パターンを重ねる
いよいよ本題です。$R = 0.5R_2$、$R = R_2$、$R = 10R_2$ の 3 つの距離でパターンを計算し、重ね描きします。比較のためリアクティブ近傍界に近い $R = 0.1R_2$ も加えます。
import numpy as np
import matplotlib.pyplot as plt
theta = np.radians(np.linspace(-25, 25, 1501))
cases = [(0.10, "$0.1R_2$（かなり近い）", "tab:purple"),
(0.50, "$0.5R_2$（フレネル領域）", "tab:orange"),
(1.00, "$R_2 = 2D^2/\\lambda$", "tab:red"),
(10.0, "$10R_2$（十分な遠方界）", "tab:blue")]
plt.figure(figsize=(10, 6))
for f, label, c in cases:
R = f * R2
E = np.abs(aperture_field(R, theta))
EdB = 20*np.log10(E / E.max()) # 各距離でピークを0dBに規格化
plt.plot(np.degrees(theta), EdB, color=c, lw=2, label=label)
plt.axhline(-13.26, color="k", ls="--", lw=1,
label="一様開口の理論SLL $-13.26$ dB")
plt.ylim(-45, 2); plt.xlim(-25, 25)
plt.xlabel("観測角 $\\theta$ [度]"); plt.ylabel("正規化パターン [dB]")
plt.title("測定距離によって放射パターンが変わる（$D = 10\\lambda$ 一様開口）")
plt.grid(alpha=0.3); plt.legend(fontsize=9, loc="upper right")
plt.tight_layout(); plt.show()
このグラフがフレネル領域の本質を示しています。$10R_2$ の青線は理想的な $\mathrm{sinc}$ 型パターンで、$\pm 5.74^\circ$ に深いヌルがあり、第一サイドローブは理論値どおり $-13.26$ dB です。$R_2$（赤）ではヌルが $-22$ dB までしか落ちず、サイドローブは $-12.94$ dB とわずかに持ち上がっています。$0.5R_2$（橙）になるとヌルは $-15.9$ dB まで埋まり、サイドローブは $-12.04$ dB へ上昇して、メインローブとの境目が曖昧になり始めます。$0.1R_2$（紫）に至ってはヌルもサイドローブも完全に消え、メインローブは $\pm 6^\circ$ にわたる平坦な台地状に広がります。しかも最大値は軸上ではなく $\pm 5.9^\circ$ にあり、軸上は $-0.59$ dB、$\pm 2.5^\circ$ には $-0.86$ dB の浅い落ち込みができています。開口端の位相誤差が $225^\circ$ にも達し、開口中心からの寄与と端からの寄与が互いに打ち消し始めるためです。この距離ではもはや「放射パターンを測っている」とは言えません。
重要なのは、ピーク方向とヌルの角度位置はほとんど動いていないという点です。動くのは主に「深さ」と「太さ」です。二次位相誤差は開口の左右対称な誤差なので、ビームを傾けたりせず、コントラストだけを落とします。これはデフォーカスしたカメラの写真で、被写体の位置は変わらずコントラストだけが失われるのと同じです。
サイドローブ・ヌル・利得の距離依存を定量化する
目で見た変化を数値で押さえましょう。距離を連続的に振って、第一サイドローブレベル、第一ヌルの深さ、軸上利得の低下量を測ります。
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import minimize_scalar
def metrics(R):
"""第一ヌル深さ・第一サイドローブレベル・軸上利得損失を測る"""
p0 = abs(aperture_field(R, 0.0))[0]
f = lambda t: 20*np.log10(abs(aperture_field(R, t, N=20001))[0] / p0)
null = minimize_scalar(f, bounds=np.radians([4.5, 7.0]),
method="bounded", options={"xatol": 1e-10})
slb = minimize_scalar(lambda t: -f(t), bounds=np.radians([7.0, 9.5]),
method="bounded", options={"xatol": 1e-10})
return null.fun, -slb.fun, p0 * R
ratios = np.array([0.25, 0.5, 1.0, 2.0, 5.0, 10.0])
ref = abs(aperture_field(1000*R2, 0.0))[0] * 1000 * R2 # 遠方界の基準値
print(" R/R2 | 端の位相誤差 | ヌル深さ | 第一SLL | 軸上利得損失")
rows = []
for f in ratios:
R = f * R2
nd, sll, axial = metrics(R)
loss = 20*np.log10(axial / ref)
rows.append((f, np.degrees(k*(D/2)**2/(2*R)), nd, sll, loss))
print(f"{f:5.2f} | {rows[-1][1]:9.2f}° | {nd:8.1f} dB | {sll:6.2f} dB | {loss:+7.3f} dB")
実行結果は次の表のとおりです。
| $R/R_2$ | 開口端の二次位相誤差 | 第一ヌル深さ | 第一サイドローブレベル | 軸上利得損失 |
|---|---|---|---|---|
| 0.25 | 90.0° | $-9.6$ dB | $-9.06$ dB | $-0.977$ dB |
| 0.50 | 45.0° | $-15.9$ dB | $-12.04$ dB | $-0.242$ dB |
| 1.00 | 22.5° | $-22.0$ dB | $-12.94$ dB | $-0.060$ dB |
| 2.00 | 11.25° | $-28.0$ dB | $-13.18$ dB | $-0.015$ dB |
| 5.00 | 4.50° | $-36.0$ dB | $-13.25$ dB | $-0.002$ dB |
| 10.0 | 2.25° | $-42.0$ dB | $-13.26$ dB | $-0.001$ dB |
この表を距離について連続的に描き直すと、次のようになります。
3 つのパネルは「同じ距離でも、測りたい量によって合格・不合格が変わる」ことを示しています。右の軸上利得は $R/R_2 = 0.9$ あたりで早々に $-0.1$ dB の線を上回って収束してしまうのに対し、中央の第一サイドローブは $R_2$ を過ぎてもまだ理論値 $-13.26$ dB の黒破線に届かず、$3R_2$ 付近でようやく重なります。左のヌル深さに至っては片対数上で完全な直線のまま改善し続け、頭打ちになる気配がありません。利得は $R_2$ で足りる／サイドローブは数倍必要／ヌルは 10 倍必要という指針は、この 3 枚の傾きの違いをそのまま言葉にしたものです。
この表から 3 つの重要な事実が読み取れます。第 1 に、$R = R_2$ での軸上利得損失は $0.060$ dB で、$\pi/8$ 基準が「利得測定なら $0.1$ dB 以内」という実用的な保証を与えていることが確認できます。第 2 に、第一サイドローブレベルは $R_2$ で $-12.94$ dB、理論値 $-13.26$ dB に対して誤差 $0.32$ dB です。サイドローブを $0.1$ dB の精度で測りたいなら $R_2$ では足りず、$2R_2$ 以上が必要になります。第 3 に、ヌル深さは距離の 1 桁ごとに正確に $20$ dB 改善しています（$R_2$ で $-22.0$、$10R_2$ で $-42.0$、比較計算では $100R_2$ で $-62.0$ dB）。これは、ヌルの埋まり具合が二次位相誤差の 1 次に比例する、つまり $1/R$ に比例することを意味します。ヌルの深さを $-40$ dB 以上取りたいレーダーや干渉抑圧の用途では、$2D^2/\lambda$ の 10 倍の距離が必要になる、という設計指針がここから出てきます。
3 領域の境界を口径と周波数の平面に描く
最後に、実務で使いやすい形として、口径 $D$ と周波数を軸にとって境界距離を等高線で描いてみます。
import numpy as np
import matplotlib.pyplot as plt
c = 2.99792458e8
D_m = np.logspace(np.log10(0.05), np.log10(10.0), 300) # 口径 0.05〜10 m
f_hz = np.logspace(np.log10(0.5e9), np.log10(100e9), 300) # 0.5〜100 GHz
DD, FF = np.meshgrid(D_m, f_hz)
LL = c / FF
R2_map = 2 * DD**2 / LL
plt.figure(figsize=(9, 6))
levels = [0.1, 1, 10, 100, 1e3, 1e4, 1e5]
cs = plt.contour(DD, FF/1e9, R2_map, levels=levels,
colors="tab:blue", linewidths=1.5)
plt.clabel(cs, fmt=lambda v: f"{v:g} m", fontsize=9)
plt.contourf(DD, FF/1e9, np.log10(R2_map), levels=30, cmap="viridis", alpha=0.5)
plt.colorbar(label="$\\log_{10}(2D^2/\\lambda)$ [m]")
for (d, fg, name) in [(0.3, 2.4, "Wi-Fi端末"), (0.6, 12, "BSパラボラ"),
(3.0, 10, "Xバンド地球局"), (5.0, 30, "Kaバンド大口径")]:
plt.plot(d, fg, "r*", ms=14)
plt.annotate(name, (d, fg), textcoords="offset points",
xytext=(8, 6), color="white", fontsize=10)
plt.xscale("log"); plt.yscale("log")
plt.xlabel("開口径 $D$ [m]"); plt.ylabel("周波数 [GHz]")
plt.title("遠方界距離 $2D^2/\\lambda$ のマップ（★は代表的なアンテナ）")
plt.tight_layout(); plt.show()
このマップは、アンテナ測定の現実的な難しさを一目で示します。等高線が右上へ行くほど急激に大きな値になり、口径 3 m・10 GHz の X バンド地球局はすでに $600$ m の等高線付近、口径 5 m・30 GHz の Ka バンドアンテナは $5000$ m の等高線に乗っています。逆に左下の Wi-Fi 端末アンテナは $1.4$ m で、机の上でも遠方界測定が成立します。等高線の間隔が対数軸で等間隔になっているのは $R_2 \propto D^2 f$ という単純なべき則の帰結で、口径を 2 倍にすれば距離は 4 倍、周波数を 2 倍にすれば距離も 2 倍になります。
フレネル領域とフレネルゾーンを並べて描く
用語の混同を防ぐため、2 つの「フレネル」を同じ図に並べます。
import numpy as np
import matplotlib.pyplot as plt
fig, ax = plt.subplots(1, 2, figsize=(13, 5))
# 左：1つのアンテナのまわりの3領域（半径は対数スケールで描く）
Da, lam_m, rmin = 1.0, 0.03, 0.1
r1, r2 = 0.62*np.sqrt(Da**3/lam_m), 2*Da**2/lam_m
lg = lambda r: np.log10(r/rmin) # 実距離 → 描画半径
for rad, col, lab in [(lg(600), "#cfe8ff", f"フラウンホーファー領域\n遠方界（$R>${r2:.1f} m）"),
(lg(r2), "#ffe0b3", f"フレネル領域＝放射近傍界\n({r1:.2f} m 〜 {r2:.1f} m)"),
(lg(r1), "#ffb3b3", f"リアクティブ近傍界\n（$R<${r1:.2f} m）")]:
ax[0].add_patch(plt.Circle((0, 0), rad, color=col, zorder=1))
for rad, lab in [(lg(600), "遠方界"), (lg(r2), "フレネル領域"), (lg(r1), "リアクティブ近傍界")]:
ax[0].annotate(lab, xy=(0, rad), xytext=(0, rad+0.12), ha="center", fontsize=10, zorder=5)
ax[0].plot([0, 0], [-0.25, 0.25], "k-", lw=6, zorder=4)
ax[0].text(0, -lg(600)-0.35, f"開口 $D$=1 m, 10 GHz 境界: {r1:.2f} m と {r2:.1f} m\n（半径は対数スケール）",
ha="center", fontsize=10)
ax[0].set_xlim(-lg(600)*1.15, lg(600)*1.15); ax[0].set_ylim(-lg(600)*1.35, lg(600)*1.2)
ax[0].set_aspect("equal"); ax[0].axis("off")
ax[0].set_title("フレネル「領域」: 1本のアンテナのまわりの距離区分")
# 右：伝搬路のまわりの第1フレネルゾーン
Ltot, lam2 = 10e3, 0.06
d1 = np.linspace(1, Ltot-1, 800)
rF = np.sqrt(lam2 * d1 * (Ltot-d1) / Ltot)
ax[1].fill_between(d1/1e3, rF, -rF, color="#cfe8ff", label="第1フレネルゾーン")
ax[1].plot([0, Ltot/1e3], [0, 0], "k--", lw=1.2, label="見通し線")
ax[1].plot(0, 0, "ro", ms=9); ax[1].plot(Ltot/1e3, 0, "ro", ms=9)
ax[1].text(0, 12, "送信点", ha="center"); ax[1].text(10, 12, "受信点", ha="center")
ax[1].set_xlabel("距離 [km]"); ax[1].set_ylabel("見通し線からの半径 [m]")
ax[1].set_title("フレネル「ゾーン」: 伝搬路のまわりの回転楕円体（5 GHz, 10 km）")
ax[1].legend(loc="lower left", fontsize=9); ax[1].grid(alpha=0.3)
plt.tight_layout(); plt.show()
左右で対象がまったく違うことが視覚的にはっきりします。左のフレネル領域はアンテナを中心とした同心球で、境界はアンテナ自身の口径と波長だけで決まります（$D=1$ m・10 GHz なら $3.6$ m から $66.7$ m）。右のフレネルゾーンは送受信 2 点を焦点とする回転楕円体で、経路の中央で最も太くなり、$5$ GHz・$10$ km のリンクでは中央の半径が約 $12.2$ m になります（横軸が km、縦軸が m なので図では大きく膨らんで見えますが、実際は長さ $10\,$km に対して太さ $12$ m の、極端に細長い葉巻形です）。前者は「アンテナ単体の測定距離をどう取るか」、後者は「伝搬路にどれだけ空間を空けるか」という、答えるべき問いが根本的に異なります。名前が似ているだけで混同しないよう注意してください。
理論の導出と数値実験が揃いました。最後に、これらの結果を現場でどう使うか — 測定距離をどう決め、距離が足りないときに何をするか、そしてフレネル領域を逆手に取る使い方 — をまとめておきましょう。
実務での使いどころ
測定距離をどう決めるか
これまでの数値実験から、測りたい量ごとに必要な距離が違うことがわかりました。実務的な目安をまとめると次のようになります。
利得だけを測りたいなら $R = 2D^2/\lambda$ で十分です。軸上利得損失は $0.06$ dB で、通常の利得測定の不確かさ（$\pm 0.2$ dB 程度）に埋もれます。ただし、この距離では自由空間伝搬損失の式（フリスの伝達公式）に $1/R^2$ 則を適用する前提が成り立っていることも確認しておく必要があります。
メインローブの形（ビーム幅）を測りたいなら、$R_2$ でも $-3$ dB ビーム幅の誤差は $5.087^\circ$ 対 $5.070^\circ$、率にして $0.3\,\%$ です。十分実用になります。
サイドローブレベルを $\pm 0.1$ dB で測りたいなら $2R_2$ 以上が必要です。表のとおり $R_2$ では $0.32$ dB の誤差が出ます。低サイドローブアンテナ（$-30$ dB 級）を評価する場合は、要求精度に応じてさらに遠くする必要があります。
深いヌルや相殺特性を評価したいなら $10R_2$ 級が必要です。$-40$ dB のヌルを正しく測るには、$R_2$ の 10 倍の距離を確保しなければなりません。アダプティブアレーの干渉抑圧性能を実測する場面などがこれに該当します。
距離が足りないときの対処
$R_2$ が現実的でないほど大きいときは、次のいずれかを選びます。
近傍界測定＋NF-FF 変換では、フレネル領域どころかもっと近い距離（典型的には数波長）でプローブを走査して振幅と位相を測り、それをフーリエ変換で遠方界パターンに変換します。平面走査・円筒走査・球面走査の 3 方式があり、大口径アンテナの標準的な測定法になっています。近傍界の測定には位相情報が必須で、これが振幅だけを測ればよい遠方界測定との最大の違いです。
コンパクトアンテナテストレンジ（CATR）では、放物面反射鏡を使って点源からの球面波を平面波に変換し、数 m の暗室内に「疑似遠方界」を作ります。反射鏡の焦点に給電ホーンを置くと、鏡の開口から出る波はほぼ平面波になり、被測定アンテナは実質的に無限遠を見ることになります。
フレネル領域測定＋距離補正という手も使えます。測定距離が既知なら、二次位相項 $\exp(-jkx'^2/(2R))$ の効果は計算で除去できます。実際、フレネル領域で測った複数距離のパターンから遠方界パターンを推定する手法（外挿法や、二次位相を逆補正する手法）が使われています。デフォーカスした画像をデコンボリューションで復元するのと同じ発想です。
近傍界を積極的に使う
フレネル領域は「不便な中途半端な領域」ではなく、積極的に利用される領域でもあります。開口に意図的に二次位相を与えて $\exp(+jkx'^2/(2R_f))$ とすれば、距離 $R_f$ に電磁的な焦点を作れます。これがビームフォーカシングで、マイクロ波無線電力伝送（WPT）、近距離ミリ波イメージング、医療用ハイパーサーミアなどで使われています。遠方界では「ビームは広がるだけ」ですが、フレネル領域なら「一点に絞れる」というのは大きな違いです。焦点でのスポットサイズはおおよそ $\lambda R_f / D$ 程度になり、開口が大きく波長が短いほど鋭く絞れます。
また、大型フェーズドアレーが至近距離の目標を扱う場面（近距離レーダー、車載レーダーの校正など）では、素子ごとの位相を「平面波の到来」ではなく「球面波の到来」を前提に設定しないとビームが目標に当たりません。これも二次位相項を明示的に扱う設計で、フレネル領域の理解がそのまま実装になります。
こうして見ると、フレネル領域は「遠方界の劣化版」ではなく、二次位相という自由度を持つ独立した設計空間であることがわかります。ここまでの議論を最後に整理しておきましょう。
まとめ
本記事では、アンテナ周辺の空間を距離で 3 つに切り分ける考え方と、その境界式の導出を扱いました。
- アンテナのまわりは、リアクティブ近傍界（$R < 0.62\sqrt{D^3/\lambda}$）、フレネル領域＝放射近傍界（$0.62\sqrt{D^3/\lambda} < R < 2D^2/\lambda$）、フラウンホーファー領域＝遠方界（$R > 2D^2/\lambda$）の 3 つに分かれる
- 境界を決めるのは、開口面上の点から観測点までの距離 $r = \sqrt{R^2 - 2Rx'\sin\theta + x'^2}$ のテイラー展開である。$r \simeq R - x'\sin\theta + x'^2\cos^2\theta/(2R) + x'^3\sin\theta\cos^2\theta/(2R^2)$ の各項が、共通位相・フーリエ核・二次位相誤差・三次位相誤差に対応する
- 二次位相項の最大値（$\theta = 0$、$x' = D/2$）が $\pi/8$（$= \lambda/16$ の経路差）を超えない条件から、遠方界距離 $R \ge 2D^2/\lambda$ が出る
- 三次位相項の最大値（$\sin\theta\cos^2\theta$ の最大 $2/(3\sqrt{3})$、$\theta \approx 35.26^\circ$）に同じ $\pi/8$ 基準を課すと、内側境界 $R \ge \sqrt{2/(3\sqrt{3})}\sqrt{D^3/\lambda} = 0.6204\sqrt{D^3/\lambda}$ が出る。$0.62$ は経験値ではなく閉じた形を持つ幾何学的定数である
- フレネル領域で二次位相が残ることは、開口にデフォーカスしたレンズを貼るのと等価である。だからビームが太くなり、サイドローブが持ち上がり、ヌルが埋まり、パターンが距離ごとに変わる
- 数値実験（$D = 10\lambda$ 一様開口）では、$R = 2D^2/\lambda$ で軸上利得損失 $0.060$ dB、第一サイドローブ $-12.94$ dB（理論値 $-13.26$ dB）、第一ヌル $-22.0$ dB。ヌル深さは距離が 10 倍になるごとに $20$ dB 改善する
- したがって、利得なら $2D^2/\lambda$ で十分、サイドローブ精度が要るなら $2R_2$、深いヌルを測るなら $10R_2$ という測定距離の指針が立つ
- 「フレネル領域」（1 本のアンテナのまわりの距離区分）と「フレネルゾーン」（伝搬路のまわりの回転楕円体）は別概念である。共通するのは二次の経路差を扱うという数学的骨格だけ
3 領域の切り分けは、アンテナ測定の実務だけでなく、開口面アンテナの設計、近距離レーダー、無線電力伝送のビームフォーカシングまで広く効いてきます。特に「開口分布のフーリエ変換が遠方界パターンになる」という関係は、開口アンテナ理論の中心にある事実です。
次のステップとして、以下の記事も参考にしてください。
- アンテナの近傍界と遠方界とは？フラウンホーファー距離を図解で理解する — 3 領域のうち近傍界と遠方界の物理（蓄積エネルギー・波動インピーダンス）を掘り下げます
- アンテナ測定の基礎 — 放射パターン・利得・インピーダンスの測定法 — 実際の測定系（電波暗室・NF-FF 変換・利得測定法）を解説します
- パラボラアンテナと開口面アンテナの理論を解説 — 大口径開口アンテナの代表格。開口効率と本記事の位相誤差の関係が見えてきます
- アレイアンテナとビームフォーミングの理論を解説 — 離散素子で同じ議論を行うと、素子ごとの位相制御と二次位相補正の話につながります
- アレイファクタの理論 — アレイアンテナの指向性パターンを数式から導出する — 本記事の連続開口積分の離散版にあたる考え方です
- フレネルゾーンと回折損失を解説 — 混同しやすい「もう一つのフレネル」を確認しておきましょう