衛星放送のパラボラアンテナを設置するとき、アンテナの「向き」だけでなく「角度の微調整」が求められます。この微調整は、電波の 偏波（polarization） を正しく受信するために必要な操作です。偏波がずれていると、アンテナの向きが完璧でも受信電力が大幅に低下してしまいます。
偏波とは、電磁波の電界ベクトルが空間中でどのような軌跡を描くかを表す性質です。「電波の振動する向き」と言い換えることもできます。偏波を理解すると、以下のような幅広い応用が見えてきます。
- 衛星通信の周波数有効利用: 同じ周波数帯で垂直偏波と水平偏波を同時に使い、伝送容量を2倍にする（偏波多重）
- レーダーによる目標識別: 偏波特性の違いから、航空機の形状や地表の植生を分類する（偏波SAR）
- 光通信・光学: 偏光フィルタ、液晶ディスプレイ、偏光顕微鏡など、光の偏波を制御する技術の基盤
本記事の内容
- 偏波とは何か — 電界ベクトルの軌跡としての理解
- 直線偏波、円偏波（右旋・左旋）、楕円偏波の分類
- ジョーンズベクトルとジョーンズ行列による偏波の記述
- ストークスパラメータとポアンカレ球
- 偏波不整合損失（偏波効率）の計算
- アンテナへの応用: 円偏波アンテナ、衛星通信での偏波利用
- Pythonによる偏波楕円の可視化とポアンカレ球の3Dプロット
前提知識
この記事を読む前に、以下の知識があると理解が深まります。
- 電磁気学の基礎（マクスウェル方程式、平面波の概念）
- ベクトルと複素数の基本的な演算
偏波とは — 電界ベクトルの軌跡
直感的な理解
ロープの一端を持って上下に振ると、波がロープに沿って進みます。このとき、波の振動方向は「上下」です。横に振れば「左右」に振動する波が進みます。電磁波も同じように、電界の振動方向を持っています。ただし、ロープの波は1次元的（上下のみ）ですが、電磁波の電界は進行方向に垂直な平面内の 任意の方向 に振動できます。この「電界がどの方向にどう振動するか」が偏波です。
より正確に言えば、偏波とは、電磁波の進行方向の一定点において、時間とともに電界ベクトルの先端が描く軌跡のことです。この軌跡が直線なら「直線偏波」、円なら「円偏波」、楕円なら「楕円偏波」と呼ばれます。
一般的な平面波の電界
$z$ 方向に伝搬する単色平面波の電界は、$x$ 成分と $y$ 成分の重ね合わせとして
$$ \bm{E}(z, t) = E_x(z, t)\, \hat{\bm{x}} + E_y(z, t)\, \hat{\bm{y}} $$
と書けます。各成分は
$$ E_x(z, t) = E_{0x} \cos(\omega t – kz + \delta_x) $$
$$ E_y(z, t) = E_{0y} \cos(\omega t – kz + \delta_y) $$
です。ここで $E_{0x}$, $E_{0y}$ は各成分の振幅（非負の実数）、$\omega$ は角周波数、$k$ は波数、$\delta_x$, $\delta_y$ はそれぞれの初期位相です。
偏波を決定するのは、2つの成分の 振幅比 $E_{0y}/E_{0x}$ と 位相差 $\delta = \delta_y – \delta_x$ の2つのパラメータです。位相差 $\delta$ が偏波の形状を劇的に変えることが、これから見る分類の鍵となります。
偏波楕円の方程式
$z = 0$ の位置で電界ベクトルの軌跡を追うと、一般に楕円を描きます。これを示すために、$E_x$ と $E_y$ の式から時間 $t$ を消去します。
$E_x = E_{0x}\cos(\omega t + \delta_x)$, $E_y = E_{0y}\cos(\omega t + \delta_y)$ において、$\phi = \omega t + \delta_x$ とおくと
$$ \frac{E_x}{E_{0x}} = \cos\phi, \quad \frac{E_y}{E_{0y}} = \cos(\phi + \delta) $$
ここで $\delta = \delta_y – \delta_x$ です。右辺を加法定理で展開すると
$$ \frac{E_y}{E_{0y}} = \cos\phi\cos\delta – \sin\phi\sin\delta $$
$\cos\phi = E_x/E_{0x}$ を代入し、$\sin\phi = \pm\sqrt{1 – E_x^2/E_{0x}^2}$ を使って $\sin\phi$ を消去します。両辺を整理すると、次の偏波楕円の方程式が得られます。
$$ \left(\frac{E_x}{E_{0x}}\right)^2 + \left(\frac{E_y}{E_{0y}}\right)^2 – 2\frac{E_x}{E_{0x}}\frac{E_y}{E_{0y}}\cos\delta = \sin^2\delta $$
この式は、$E_x$-$E_y$ 平面上の楕円を表しています。$\delta$ と振幅比によって楕円の形状と傾きが決まります。特殊な場合として、$\delta = 0$ や $\delta = \pi$ のとき楕円は直線に退化し、$\delta = \pm\pi/2$ かつ $E_{0x} = E_{0y}$ のとき円になります。
偏波楕円の一般的な形がわかったところで、具体的な偏波の種類を分類していきましょう。
偏波の分類 — 直線・円・楕円
直線偏波
最もシンプルな偏波は、電界ベクトルが常に一定の方向に振動する 直線偏波（linear polarization） です。これは2つの成分の位相差 $\delta$ が $0$ または $\pi$ の整数倍のときに生じます。
水平偏波（$E_{0y} = 0$）: 電界が $x$ 軸方向のみに振動します。
$$ \bm{E} = E_{0x}\cos(\omega t – kz)\, \hat{\bm{x}} $$
垂直偏波（$E_{0x} = 0$）: 電界が $y$ 軸方向のみに振動します。
斜め直線偏波（$\delta = 0$）: $E_{0x} \neq 0$, $E_{0y} \neq 0$ で位相差がゼロの場合、電界ベクトルは
$$ \frac{E_y}{E_x} = \frac{E_{0y}}{E_{0x}} = \text{const.} $$
という一定の傾きの直線上を往復します。傾き角 $\tau = \arctan(E_{0y}/E_{0x})$ です。$\delta = \pi$ の場合は $E_y/E_x = -E_{0y}/E_{0x}$ で、傾きが反転します。
直線偏波は地上テレビ放送や多くの無線LANで使われています。アンテナ（ダイポールアンテナなど）を水平に設置すれば水平偏波、垂直に設置すれば垂直偏波を送受信します。
円偏波
2つの成分の振幅が等しく（$E_{0x} = E_{0y} = E_0$）、位相差が $\delta = \pm\pi/2$ のとき、電界ベクトルの先端は円を描きます。これが 円偏波（circular polarization） です。
右旋円偏波（RHCP: Right-Hand Circular Polarization）: $\delta = -\pi/2$ のとき
$$ E_x = E_0 \cos(\omega t – kz), \quad E_y = E_0 \sin(\omega t – kz) $$
波の進行方向（$+z$ 方向）に向かって見ると、電界ベクトルは 時計回り に回転します。この回転方向は、右手の親指を進行方向に向けたとき、指の巻く方向に一致するので「右旋」と呼ばれます。
左旋円偏波（LHCP: Left-Hand Circular Polarization）: $\delta = +\pi/2$ のとき
$$ E_x = E_0 \cos(\omega t – kz), \quad E_y = -E_0 \sin(\omega t – kz) $$
進行方向に向かって 反時計回り に回転します。
円偏波は、衛星通信やGPSで広く利用されています。その理由は、円偏波は送受信アンテナの相対的な回転角度に影響を受けにくいためです。衛星と地上局の間では、電波が電離層を通過する際にファラデー回転（偏波面の回転）が発生しますが、円偏波はこの影響を受けません。
楕円偏波
直線偏波と円偏波は、楕円偏波の特殊な場合です。一般の振幅比と位相差の組み合わせでは、電界ベクトルは楕円を描きます。
楕円偏波の形状は、以下の2つのパラメータで特徴づけられます。
傾斜角（tilt angle） $\tau$: 偏波楕円の長軸が $x$ 軸となす角度。$-\pi/2 < \tau \leq \pi/2$ の範囲で定義されます。
$$ \tan 2\tau = \frac{2E_{0x}E_{0y}\cos\delta}{E_{0x}^2 – E_{0y}^2} $$
楕円率角（ellipticity angle） $\varepsilon$: 楕円のアスペクト比を表す角度。$-\pi/4 \leq \varepsilon \leq \pi/4$ の範囲で定義されます。
$$ \sin 2\varepsilon = \frac{2E_{0x}E_{0y}\sin\delta}{E_{0x}^2 + E_{0y}^2} $$
$\varepsilon = 0$ のとき直線偏波、$\varepsilon = \pm\pi/4$ のとき円偏波です。$\varepsilon > 0$ が右旋、$\varepsilon < 0$ が左旋に対応します。
楕円の長軸と短軸の比は $|\tan\varepsilon|$ で与えられます。軸比（Axial Ratio: AR） は
$$ \text{AR} = \frac{E_{\text{major}}}{E_{\text{minor}}} = \frac{1}{|\tan\varepsilon|} $$
と定義され、$1 \leq \text{AR} \leq \infty$ の範囲を取ります。$\text{AR} = 1$ が円偏波、$\text{AR} = \infty$ が直線偏波です。デシベルで表すと $\text{AR}_{\text{dB}} = 20\log_{10}(\text{AR})$ であり、実用的な円偏波アンテナでは $\text{AR}_{\text{dB}} < 3$ dB が目安です。
偏波の分類ができたところで、これらをコンパクトに数学的に記述する方法が必要です。ジョーンズベクトルは、偏波状態を2次元の複素ベクトルで表現する強力なツールです。
ジョーンズベクトルとジョーンズ行列
ジョーンズベクトルの定義
完全偏波の電界を複素表現で書くと
$$ \tilde{\bm{E}} = \begin{pmatrix} E_{0x} e^{j\delta_x} \\ E_{0y} e^{j\delta_y} \end{pmatrix} e^{j(\omega t – kz)} $$
となります。伝搬因子 $e^{j(\omega t – kz)}$ は偏波状態に影響しないため、これを除いた部分
$$ \bm{J} = \begin{pmatrix} E_{0x} e^{j\delta_x} \\ E_{0y} e^{j\delta_y} \end{pmatrix} $$
を ジョーンズベクトル（Jones vector） と呼びます。ジョーンズベクトルは偏波状態の全情報を2つの複素数で表現します。
通常、エネルギーを正規化して $|\bm{J}|^2 = |J_x|^2 + |J_y|^2 = 1$ とします。正規化後の主要な偏波のジョーンズベクトルは以下の通りです。
水平直線偏波（H）:
$$ \bm{J}_H = \begin{pmatrix} 1 \\ 0 \end{pmatrix} $$
垂直直線偏波（V）:
$$ \bm{J}_V = \begin{pmatrix} 0 \\ 1 \end{pmatrix} $$
45度直線偏波:
$$ \bm{J}_{45} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 1 \end{pmatrix} $$
右旋円偏波（RHCP）:
$$ \bm{J}_R = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ -j \end{pmatrix} $$
左旋円偏波（LHCP）:
$$ \bm{J}_L = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ j \end{pmatrix} $$
ジョーンズベクトルの便利な点は、偏波間の内積が直感的に計算できることです。2つの偏波 $\bm{J}_1$ と $\bm{J}_2$ の内積 $\bm{J}_1^H \bm{J}_2$（$H$ はエルミート転置）がゼロのとき、2つの偏波は 直交 しています。例えば、$\bm{J}_H$ と $\bm{J}_V$ は直交し、$\bm{J}_R$ と $\bm{J}_L$ も直交します。
ジョーンズ行列
光学素子や伝搬媒質が偏波状態を変化させる場合、その変換は $2 \times 2$ の複素行列（ジョーンズ行列）で表されます。入力偏波 $\bm{J}_{\text{in}}$ が素子を通過すると、出力偏波は
$$ \bm{J}_{\text{out}} = \bm{T} \bm{J}_{\text{in}} $$
となります。ここで $\bm{T}$ がジョーンズ行列です。
代表的なジョーンズ行列をいくつか示します。
水平偏光子（水平成分のみ通過）:
$$ \bm{T}_H = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix} $$
角度 $\theta$ の回転（ファラデー回転など）:
$$ \bm{T}_{\text{rot}}(\theta) = \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix} $$
1/4波長板（高速軸が水平）:
$$ \bm{T}_{\lambda/4} = e^{-j\pi/4}\begin{pmatrix} 1 & 0 \\ 0 & j \end{pmatrix} $$
1/4波長板は $y$ 成分に $\pi/2$ の位相遅れを導入します。45度直線偏波に1/4波長板を通すと
$$ \bm{T}_{\lambda/4} \bm{J}_{45} = e^{-j\pi/4} \begin{pmatrix} 1 & 0 \\ 0 & j \end{pmatrix} \frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ 1 \end{pmatrix} = \frac{e^{-j\pi/4}}{\sqrt{2}} \begin{pmatrix} 1 \\ j \end{pmatrix} $$
となり、左旋円偏波が得られます。この原理は円偏波アンテナの設計で実際に利用されます。
複数の素子を直列に配置する場合、全体のジョーンズ行列は各素子の行列の積で表されます。
$$ \bm{T}_{\text{total}} = \bm{T}_n \bm{T}_{n-1} \cdots \bm{T}_2 \bm{T}_1 $$
ジョーンズベクトルは完全偏波を簡潔に記述できますが、部分偏波（偏波度が100%未満の波）は扱えません。部分偏波やランダム偏波を含む一般的な状況を記述するには、ストークスパラメータが必要です。
ストークスパラメータとポアンカレ球
ストークスパラメータの定義
ジョーンズベクトルが「偏波の複素振幅」を直接記述するのに対し、ストークスパラメータ は「偏波の電力（強度）」で記述する手法です。電力量なので直接測定が容易であり、部分偏波も自然に扱えるという大きなメリットがあります。
4つのストークスパラメータ $S_0$, $S_1$, $S_2$, $S_3$ は次のように定義されます。
$$ S_0 = E_{0x}^2 + E_{0y}^2 $$
$$ S_1 = E_{0x}^2 – E_{0y}^2 $$
$$ S_2 = 2E_{0x} E_{0y} \cos\delta $$
$$ S_3 = 2E_{0x} E_{0y} \sin\delta $$
各パラメータの物理的意味は以下の通りです。
$S_0$ は波の 全電力 です。偏波状態に依存しない総エネルギーを表します。
$S_1$ は 水平/垂直偏波の優位性 を表します。$S_1 > 0$ なら水平偏波成分が強く、$S_1 < 0$ なら垂直偏波成分が強いことを意味します。
$S_2$ は 45度/-45度直線偏波の優位性 を表します。$S_2 > 0$ なら45度方向が、$S_2 < 0$ なら-45度方向が支配的です。
$S_3$ は 右旋/左旋円偏波の優位性 を表します。$S_3 > 0$ なら右旋、$S_3 < 0$ なら左旋が支配的です。
ストークスベクトル
4つのパラメータをまとめて ストークスベクトル と呼びます。
$$ \bm{S} = \begin{pmatrix} S_0 \\ S_1 \\ S_2 \\ S_3 \end{pmatrix} $$
完全偏波の場合、次の関係が成り立ちます。
$$ S_0^2 = S_1^2 + S_2^2 + S_3^2 $$
これは、4つのパラメータが独立ではなく、3つの自由度しかないことを意味します（偏波状態は実質的に $E_{0y}/E_{0x}$ と $\delta$ の2パラメータで決まるため、全電力 $S_0$ を加えて3つ）。
部分偏波の場合は
$$ S_0^2 \geq S_1^2 + S_2^2 + S_3^2 $$
となり、偏波度（Degree of Polarization: DOP） は
$$ \text{DOP} = \frac{\sqrt{S_1^2 + S_2^2 + S_3^2}}{S_0} $$
で定義されます。$\text{DOP} = 1$ が完全偏波、$\text{DOP} = 0$ が完全無偏波（自然光）です。
ジョーンズベクトルとの関係
ジョーンズベクトル $\bm{J} = (J_x, J_y)^T$ からストークスパラメータへの変換は
$$ S_0 = |J_x|^2 + |J_y|^2 $$
$$ S_1 = |J_x|^2 – |J_y|^2 $$
$$ S_2 = 2\,\text{Re}(J_x J_y^*) $$
$$ S_3 = 2\,\text{Im}(J_x J_y^*) $$
で与えられます。ここで $J_y^*$ は $J_y$ の複素共役です。
主要な偏波のストークスベクトルを確認しておきましょう（正規化: $S_0 = 1$）。
| 偏波 | ジョーンズベクトル | ストークスベクトル $(S_1, S_2, S_3)$ |
|---|---|---|
| 水平直線 | $(1, 0)^T$ | $(1, 0, 0)$ |
| 垂直直線 | $(0, 1)^T$ | $(-1, 0, 0)$ |
| 45度直線 | $\frac{1}{\sqrt{2}}(1, 1)^T$ | $(0, 1, 0)$ |
| -45度直線 | $\frac{1}{\sqrt{2}}(1, -1)^T$ | $(0, -1, 0)$ |
| 右旋円 | $\frac{1}{\sqrt{2}}(1, -j)^T$ | $(0, 0, 1)$ |
| 左旋円 | $\frac{1}{\sqrt{2}}(1, j)^T$ | $(0, 0, -1)$ |
ポアンカレ球
正規化されたストークスパラメータ $(S_1/S_0, S_2/S_0, S_3/S_0)$ を3次元空間の座標と見なすと、完全偏波は半径1の球面上の点に対応します。この球を ポアンカレ球（Poincare sphere） と呼びます。
ポアンカレ球の幾何学的な構造は非常に直感的です。
赤道: 全ての直線偏波が赤道上に並びます。水平偏波が $(1, 0, 0)$、垂直偏波が $(-1, 0, 0)$、45度偏波が $(0, 1, 0)$ です。
北極と南極: 北極 $(0, 0, 1)$ が右旋円偏波、南極 $(0, 0, -1)$ が左旋円偏波です。
その他の点: 赤道と極以外の全ての点が楕円偏波に対応します。北半球は右旋、南半球は左旋の楕円偏波です。
直交偏波: ポアンカレ球上で原点を挟んで対称な2点（正反対の点）は、互いに直交する偏波を表します。
偏波楕円のパラメータ $(\tau, \varepsilon)$ とポアンカレ球の座標の関係は
$$ S_1 = S_0 \cos 2\varepsilon \cos 2\tau $$
$$ S_2 = S_0 \cos 2\varepsilon \sin 2\tau $$
$$ S_3 = S_0 \sin 2\varepsilon $$
です。これは球面座標（経度 $2\tau$、緯度 $2\varepsilon$）そのものです。角度が2倍になっている点に注意してください。傾斜角 $\tau$ が $\pi/2$ 変化すると、ポアンカレ球上では $\pi$ 回転（反対側に移動）します。これは $\tau$ が $\pi/2$ 異なる2つの直線偏波（水平と垂直）が直交偏波であることと整合しています。
ストークスパラメータとポアンカレ球という記述法が揃ったところで、実用上極めて重要な「偏波不整合損失」の概念に進みましょう。送受信アンテナの偏波が一致しないとき、どれだけの電力が失われるかを定量化します。
偏波不整合損失（偏波効率）
偏波不整合の物理
アンテナが受信できる電力は、到来波の偏波とアンテナの偏波の「一致度」に依存します。完全に一致していれば電力損失はゼロですが、偏波がずれていると受信電力が低下します。極端な例として、水平偏波の波を垂直偏波のアンテナで受信しようとすると、理論上は受信電力がゼロになります。
偏波効率の定義
到来波のジョーンズベクトルを $\bm{J}_w$、受信アンテナの偏波（アンテナが最も感度の高い偏波）のジョーンズベクトルを $\bm{J}_a$ とするとき、偏波効率（偏波整合係数） $\eta_p$ は
$$ \eta_p = |\bm{J}_a^H \bm{J}_w|^2 $$
と定義されます。$\eta_p$ は $0 \leq \eta_p \leq 1$ の範囲を取り、$\eta_p = 1$ が完全整合、$\eta_p = 0$ が完全不整合（直交偏波）です。
デシベルで表した偏波不整合損失は
$$ L_{\text{pol}} = -10\log_{10}(\eta_p) \quad [\text{dB}] $$
です。$\eta_p = 0.5$（-3 dB）は、受信電力が半分に落ちることを意味します。
具体例の計算
いくつかの重要なケースで偏波効率を計算してみましょう。
ケース1: 直線偏波の角度ずれ
到来波が水平偏波 $\bm{J}_w = (1, 0)^T$、アンテナが角度 $\theta$ だけ傾いた直線偏波 $\bm{J}_a = (\cos\theta, \sin\theta)^T$ の場合
$$ \eta_p = |(\cos\theta)(1) + (\sin\theta)(0)|^2 = \cos^2\theta $$
これはマルス（Malus）の法則そのものです。$\theta = 0$ で $\eta_p = 1$（完全整合）、$\theta = 45°$ で $\eta_p = 0.5$（-3 dB）、$\theta = 90°$ で $\eta_p = 0$（完全不整合）です。
ケース2: 円偏波を直線偏波で受信
到来波が右旋円偏波 $\bm{J}_w = \frac{1}{\sqrt{2}}(1, -j)^T$、アンテナが水平直線偏波 $\bm{J}_a = (1, 0)^T$ の場合
$$ \eta_p = \left|\frac{1}{\sqrt{2}}(1)(1) + (0)(-j)\right|^2 = \left|\frac{1}{\sqrt{2}}\right|^2 = \frac{1}{2} $$
常に 3 dB の損失が発生します。これは直感的にも納得できます。円偏波は水平成分と垂直成分を均等に含むので、直線偏波のアンテナは半分のエネルギーしか捕捉できません。
ケース3: 右旋円偏波を左旋円偏波で受信
$\bm{J}_w = \frac{1}{\sqrt{2}}(1, -j)^T$（RHCP）、$\bm{J}_a = \frac{1}{\sqrt{2}}(1, j)^T$（LHCP）の場合
$$ \eta_p = \left|\frac{1}{2}\left[(1)(1) + (j)(-j)\right]\right|^2 = \left|\frac{1}{2}(1 + 1)\right|^2 $$
ここで内積を計算しましょう。$\bm{J}_a^H = \frac{1}{\sqrt{2}}(1, -j)$ なので
$$ \bm{J}_a^H \bm{J}_w = \frac{1}{2}\left[(1)(1) + (-j)(-j)\right] = \frac{1}{2}(1 + j^2) = \frac{1}{2}(1 – 1) = 0 $$
$\eta_p = 0$ です。右旋と左旋の円偏波は完全に直交しており、受信電力はゼロになります。衛星通信ではこの直交性を利用して、同じ周波数帯で右旋と左旋の2つの信号を同時に伝送する 偏波多重 を行っています。
ストークスパラメータによる偏波効率
ストークスパラメータを用いた偏波効率の表現も有用です。到来波と受信アンテナのストークスベクトルをそれぞれ $\bm{S}_w$, $\bm{S}_a$（正規化: $S_{0w} = S_{0a} = 1$）とすると
$$ \eta_p = \frac{1}{2}(1 + \hat{\bm{s}}_w \cdot \hat{\bm{s}}_a) $$
ここで $\hat{\bm{s}} = (S_1, S_2, S_3)^T / S_0$ は正規化ストークスベクトル（ポアンカレ球上の位置ベクトル）です。
この式の幾何学的な意味は明快です。ポアンカレ球上で2つの偏波が同じ点にあれば $\hat{\bm{s}}_w \cdot \hat{\bm{s}}_a = 1$ で $\eta_p = 1$、正反対の点（直交偏波）にあれば $\hat{\bm{s}}_w \cdot \hat{\bm{s}}_a = -1$ で $\eta_p = 0$ です。つまり、ポアンカレ球上での2点間の角度距離が偏波効率を直接表している のです。
偏波不整合損失の定量的な計算方法がわかったところで、これらの理論が実際のアンテナ設計でどう活用されているかを見ていきましょう。
アンテナへの応用
円偏波アンテナの設計原理
円偏波を生成するには、振幅が等しく位相が 90 度異なる2つの直交する直線偏波成分を合成する必要があります。これを実現するアンテナ設計の手法はいくつかあります。
直交給電パッチアンテナ: 正方形パッチアンテナの2つの直交する辺から、90度位相差をつけて同時に給電します。一方の辺から $E_x$ 成分を、もう一方の辺から $E_y$ 成分を、$\pi/2$ の位相差で励振することで円偏波を生成します。
シーケンシャル回転アレイ: 直線偏波の素子を $N$ 個配置し、隣接素子を $360°/N$ ずつ物理的に回転させ、給電位相も同じ量だけシフトします。たとえば4素子の場合、各素子を 0°, 90°, 180°, 270° 回転し、給電位相も 0°, 90°, 180°, 270° とします。このアレイ全体で見ると、電界ベクトルが回転する円偏波が合成されます。
摂動型パッチアンテナ: 正方形パッチの対角を切り落とす（truncated corner）か、対角にスロットを入れることで、2つの直交モードの共振周波数をわずかにずらし、結果として 90 度の位相差を実現します。この方法は単一の給電点で済むため、給電回路がシンプルになる利点があります。
衛星通信での偏波利用
衛星通信では偏波の選択が通信リンクの設計に直接影響します。
円偏波の利点: 衛星通信では、電離層のファラデー回転により直線偏波の偏波面が回転します。回転量は周波数と電離層の状態に依存し、予測が困難です。円偏波を使えば、ファラデー回転の影響を受けません。これは、円偏波が偏波面の回転に対して本質的に不変だからです。ジョーンズ行列で言えば、回転行列 $\bm{T}_{\text{rot}}(\theta)$ を右旋円偏波のジョーンズベクトルに作用させると
$$ \bm{T}_{\text{rot}}(\theta) \bm{J}_R = \frac{1}{\sqrt{2}} \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix} \begin{pmatrix} 1 \\ -j \end{pmatrix} $$
を計算すると、展開して
$$ = \frac{1}{\sqrt{2}} \begin{pmatrix} \cos\theta + j\sin\theta \\ \sin\theta – j\cos\theta \end{pmatrix} = \frac{e^{j\theta}}{\sqrt{2}} \begin{pmatrix} 1 \\ -j \end{pmatrix} = e^{j\theta} \bm{J}_R $$
となります。全体の位相 $e^{j\theta}$ がかかるだけで、偏波状態は変化しません。
偏波多重による容量倍増: 右旋円偏波と左旋円偏波は直交しているため、同じ周波数帯で2つの独立な信号を同時に伝送できます。BS/CS 衛星放送では、この偏波多重を使ってチャンネル数を2倍にしています。典型的な偏波間のアイソレーション（交差偏波識別度: XPD）は 25-30 dB 程度であり、これにより2つの偏波チャンネル間の干渉は十分に小さく保たれます。
交差偏波識別度（XPD）: 実際のアンテナや伝搬路では、偏波の純度が完全ではありません。所望偏波の受信電力 $P_{\text{co}}$ と交差偏波の受信電力 $P_{\text{cross}}$ の比
$$ \text{XPD} = 10\log_{10}\frac{P_{\text{co}}}{P_{\text{cross}}} \quad [\text{dB}] $$
が偏波多重の品質を決定します。降雨時には雨滴の非球形性により偏波回転が生じ、XPD が劣化します。これを 降雨偏波減衰 と呼び、周波数が高いほど影響が大きくなります。
GPSにおける円偏波
GPS衛星は右旋円偏波でL1信号（1575.42 MHz）とL2信号（1227.60 MHz）を送信しています。受信機も右旋円偏波のアンテナを使用します。円偏波を採用する理由は、電離層のファラデー回転の影響を回避するためだけでなく、地上の建物などからの反射波（マルチパス）の影響を軽減するためでもあります。
電磁波が金属面で反射すると、円偏波の旋回方向が反転します。つまり、RHCP の直接波が反射すると LHCP になります。RHCP アンテナは LHCP の反射波に対して偏波不整合（$\eta_p = 0$）となるため、反射波の影響が大幅に抑制されます。この特性は、精密測位において非常に有利です。
理論の解説が一通り完了しました。ここからは Python を使って、偏波楕円の可視化とポアンカレ球の3Dプロットを実装していきましょう。
Pythonでの実装 — 偏波の可視化
偏波楕円の描画
まず、振幅比と位相差を変えたときの偏波楕円がどう変化するかを可視化します。
import numpy as np
import matplotlib.pyplot as plt
def polarization_ellipse(E0x, E0y, delta, n_points=500):
"""偏波楕円の座標を計算"""
t = np.linspace(0, 2 * np.pi, n_points)
Ex = E0x * np.cos(t)
Ey = E0y * np.cos(t + delta)
return Ex, Ey
fig, axes = plt.subplots(2, 3, figsize=(14, 9))
# 6つの代表的な偏波状態
cases = [
(1.0, 0.0, 0, 'Horizontal LP\n($E_{0y}=0$)'),
(1.0, 1.0, 0, '45-deg LP\n($\\delta=0$)'),
(1.0, 0.5, 0, 'Tilted LP\n($E_{0y}/E_{0x}=0.5, \\delta=0$)'),
(1.0, 1.0, -np.pi/2, 'RHCP\n($\\delta=-\\pi/2$)'),
(1.0, 1.0, np.pi/2, 'LHCP\n($\\delta=+\\pi/2$)'),
(1.0, 0.6, np.pi/4, 'Elliptical\n($E_{0y}/E_{0x}=0.6, \\delta=\\pi/4$)'),
]
for idx, (E0x, E0y, delta, title) in enumerate(cases):
ax = axes[idx // 3][idx % 3]
Ex, Ey = polarization_ellipse(E0x, E0y, delta)
ax.plot(Ex, Ey, 'b-', linewidth=2)
# 回転方向を矢印で示す（円偏波・楕円偏波の場合）
if abs(delta) > 0.01 and E0y > 0.01:
n_arr = len(Ex) // 4
ax.annotate('', xy=(Ex[n_arr+1], Ey[n_arr+1]),
xytext=(Ex[n_arr], Ey[n_arr]),
arrowprops=dict(arrowstyle='->', color='red', lw=2))
ax.set_xlim(-1.3, 1.3)
ax.set_ylim(-1.3, 1.3)
ax.set_aspect('equal')
ax.set_xlabel('$E_x$', fontsize=11)
ax.set_ylabel('$E_y$', fontsize=11)
ax.set_title(title, fontsize=11)
ax.grid(True, alpha=0.3)
ax.axhline(y=0, color='k', linewidth=0.5)
ax.axvline(x=0, color='k', linewidth=0.5)
plt.suptitle('Polarization Ellipses for Various States', fontsize=14, y=1.01)
plt.tight_layout()
plt.savefig('polarization_ellipses.png', dpi=150, bbox_inches='tight')
plt.show()
6つのパネルに、代表的な偏波状態の偏波楕円が描かれています。左上の水平偏波は $x$ 軸上の直線です。中央上の45度直線偏波は対角線方向の直線です。右上の斜め直線偏波は、振幅比 $E_{0y}/E_{0x} = 0.5$ に対応した傾きを持つ直線です。左下の RHCP と中央下の LHCP はいずれも正円ですが、赤矢印で示された回転方向が逆です。右下の楕円偏波は、振幅比が1でなく位相差も $\pi/2$ でない場合の一般的な楕円です。位相差 $\delta$ がゼロから変化するに従い、直線から楕円、そして $\delta = \pm\pi/2$ かつ等振幅で円に至るという偏波の連続的な変化が視覚的に理解できます。
位相差による偏波変化のアニメーション的可視化
位相差 $\delta$ を連続的に変化させたときの偏波楕円の遷移を1枚のグラフにまとめます。
import numpy as np
import matplotlib.pyplot as plt
fig, ax = plt.subplots(figsize=(8, 8))
E0x = 1.0
E0y = 0.8
delta_values = np.linspace(0, 2 * np.pi, 13)[:-1] # 0からほぼ2piまで12点
colors = plt.cm.hsv(np.linspace(0, 1, len(delta_values)))
t = np.linspace(0, 2 * np.pi, 500)
for i, delta in enumerate(delta_values):
Ex = E0x * np.cos(t)
Ey = E0y * np.cos(t + delta)
label = f'$\\delta = {np.degrees(delta):.0f}°$'
ax.plot(Ex, Ey, color=colors[i], linewidth=1.5, alpha=0.8, label=label)
ax.set_xlim(-1.3, 1.3)
ax.set_ylim(-1.3, 1.3)
ax.set_aspect('equal')
ax.set_xlabel('$E_x$', fontsize=12)
ax.set_ylabel('$E_y$', fontsize=12)
ax.set_title(f'Polarization Ellipses: $E_{{0x}}={E0x}$, $E_{{0y}}={E0y}$, varying $\\delta$',
fontsize=13)
ax.legend(fontsize=8, loc='upper right', ncol=2)
ax.grid(True, alpha=0.3)
ax.axhline(y=0, color='k', linewidth=0.5)
ax.axvline(x=0, color='k', linewidth=0.5)
plt.tight_layout()
plt.savefig('polarization_delta_sweep.png', dpi=150, bbox_inches='tight')
plt.show()
位相差 $\delta$ を 0° から 330° まで 30° 刻みで変化させたときの偏波楕円が、虹色のグラデーションで重ねて描かれています。$\delta = 0°$ と $\delta = 180°$ では楕円が直線に退化し（ただし傾きが異なる）、$\delta = 90°$ と $\delta = 270°$ 付近では楕円が最も膨らんで円に近づきます（$E_{0y}/E_{0x} = 0.8$ なので完全な円にはなりません）。位相差だけでこれほど多様な偏波状態が生まれることが、この図から直感的にわかります。
ポアンカレ球の3Dプロット
偏波状態をポアンカレ球上に3Dプロットし、各偏波の幾何学的な配置を確認します。
import numpy as np
import matplotlib.pyplot as plt
fig = plt.figure(figsize=(10, 10))
ax = fig.add_subplot(111, projection='3d')
# ポアンカレ球（ワイヤフレーム）
u = np.linspace(0, 2 * np.pi, 50)
v = np.linspace(0, np.pi, 30)
x_sphere = np.outer(np.cos(u), np.sin(v))
y_sphere = np.outer(np.sin(u), np.sin(v))
z_sphere = np.outer(np.ones_like(u), np.cos(v))
ax.plot_surface(x_sphere, y_sphere, z_sphere, alpha=0.05, color='lightblue')
ax.plot_wireframe(x_sphere, y_sphere, z_sphere, alpha=0.1, color='gray',
linewidth=0.3)
# 赤道（直線偏波の軌跡）
theta_eq = np.linspace(0, 2 * np.pi, 100)
ax.plot(np.cos(theta_eq), np.sin(theta_eq), np.zeros_like(theta_eq),
'k-', linewidth=1, alpha=0.5)
# 座標軸
ax.plot([-1.3, 1.3], [0, 0], [0, 0], 'k-', linewidth=0.5, alpha=0.3)
ax.plot([0, 0], [-1.3, 1.3], [0, 0], 'k-', linewidth=0.5, alpha=0.3)
ax.plot([0, 0], [0, 0], [-1.3, 1.3], 'k-', linewidth=0.5, alpha=0.3)
# 代表的な偏波状態をプロット
def jones_to_stokes(Jx, Jy):
"""ジョーンズベクトルから正規化ストークスパラメータを計算"""
S0 = abs(Jx)**2 + abs(Jy)**2
S1 = (abs(Jx)**2 - abs(Jy)**2) / S0
S2 = 2 * np.real(Jx * np.conj(Jy)) / S0
S3 = 2 * np.imag(Jx * np.conj(Jy)) / S0
return S1, S2, S3
# 各偏波状態
polarizations = {
'H (horizontal)': (1, 0),
'V (vertical)': (0, 1),
'45-deg': (1/np.sqrt(2), 1/np.sqrt(2)),
'-45-deg': (1/np.sqrt(2), -1/np.sqrt(2)),
'RHCP': (1/np.sqrt(2), -1j/np.sqrt(2)),
'LHCP': (1/np.sqrt(2), 1j/np.sqrt(2)),
}
colors_pol = ['red', 'blue', 'green', 'orange', 'purple', 'cyan']
markers_pol = ['o', 's', '^', 'v', '*', 'D']
for (name, (Jx, Jy)), color, marker in zip(polarizations.items(), colors_pol, markers_pol):
S1, S2, S3 = jones_to_stokes(Jx, Jy)
ax.scatter([S1], [S2], [S3], color=color, marker=marker, s=150,
edgecolors='black', linewidth=0.5, zorder=10, label=name)
ax.text(S1*1.15, S2*1.15, S3*1.15, name.split('(')[0].strip(),
fontsize=8, ha='center')
# 楕円偏波の軌跡（E0y/E0xを固定、deltaを変化）
delta_range = np.linspace(0, 2*np.pi, 200)
E0x_trace, E0y_trace = 1.0, 0.6
S1_trace = (E0x_trace**2 - E0y_trace**2) / (E0x_trace**2 + E0y_trace**2)
S2_trace = 2*E0x_trace*E0y_trace*np.cos(delta_range) / (E0x_trace**2 + E0y_trace**2)
S3_trace = 2*E0x_trace*E0y_trace*np.sin(delta_range) / (E0x_trace**2 + E0y_trace**2)
ax.plot(S1_trace, S2_trace, S3_trace, 'm-', linewidth=2, alpha=0.6,
label=f'$E_{{0y}}/E_{{0x}}$={E0y_trace}')
ax.set_xlabel('$S_1$', fontsize=12)
ax.set_ylabel('$S_2$', fontsize=12)
ax.set_zlabel('$S_3$', fontsize=12)
ax.set_title('Poincare Sphere', fontsize=14)
ax.legend(fontsize=8, loc='upper left')
ax.view_init(elev=20, azim=45)
plt.tight_layout()
plt.savefig('poincare_sphere.png', dpi=150, bbox_inches='tight')
plt.show()
ポアンカレ球上に6つの基本偏波状態がプロットされています。水平偏波（H, 赤丸）と垂直偏波（V, 青四角）は $S_1$ 軸の正と負に、45度偏波と-45度偏波は $S_2$ 軸の正と負に、RHCP（紫星）と LHCP（シアンダイヤ）は $S_3$ 軸の正と負に配置されています。これら6点はちょうど座標軸上の点であり、3つの直交偏波ペア（H/V、45/-45、R/L）がポアンカレ球上で正反対に位置していることが視覚的に確認できます。マゼンタの曲線は、振幅比 $E_{0y}/E_{0x} = 0.6$ を固定して位相差を0から $2\pi$ まで変化させた軌跡です。この軌跡は球面上の円（緯度一定ではないが $S_1$ 一定の円）を描き、偏波楕円の形状変化をポアンカレ球上で追跡できます。
偏波不整合損失の可視化
最後に、直線偏波のアンテナ角度ずれと円偏波の軸比による偏波損失を可視化します。
import numpy as np
import matplotlib.pyplot as plt
fig, axes = plt.subplots(1, 2, figsize=(14, 6))
# 左: 直線偏波の角度ずれ vs 偏波効率
theta_deg = np.linspace(0, 90, 200)
theta_rad = np.radians(theta_deg)
eta_linear = np.cos(theta_rad) ** 2
loss_linear_db = -10 * np.log10(np.maximum(eta_linear, 1e-10))
ax = axes[0]
ax.plot(theta_deg, eta_linear, 'b-', linewidth=2, label='$\\eta_p = \\cos^2\\theta$')
ax.set_xlabel('Angle mismatch $\\theta$ [deg]', fontsize=12)
ax.set_ylabel('Polarization efficiency $\\eta_p$', fontsize=12)
ax.set_title('Linear Polarization Mismatch', fontsize=13)
ax2 = ax.twinx()
ax2.plot(theta_deg, loss_linear_db, 'r--', linewidth=2, label='Loss [dB]')
ax2.set_ylabel('Polarization loss [dB]', fontsize=12, color='red')
ax2.tick_params(axis='y', labelcolor='red')
ax2.set_ylim(0, 30)
# 重要な点にマーカー
for theta_mark in [15, 30, 45, 60]:
eta_mark = np.cos(np.radians(theta_mark))**2
ax.plot(theta_mark, eta_mark, 'ko', markersize=6)
ax.annotate(f'{theta_mark}deg: {eta_mark:.3f} ({-10*np.log10(eta_mark):.1f}dB)',
xy=(theta_mark, eta_mark),
xytext=(theta_mark + 5, eta_mark + 0.05),
fontsize=8, arrowprops=dict(arrowstyle='->', color='gray'))
ax.grid(True, alpha=0.3)
ax.set_xlim(0, 90)
ax.set_ylim(0, 1.05)
# 右: 軸比(AR) vs 円偏波の偏波効率
ar_db = np.linspace(0, 10, 200)
ar_linear = 10 ** (ar_db / 20)
# 楕円偏波の到来波を理想的な円偏波アンテナで受信した場合の偏波効率
# AR -> 楕円率角 epsilon = arctan(1/AR)
epsilon = np.arctan(1 / ar_linear)
# RHCP受信の場合、到来波が同旋のとき
eta_cp = np.cos(np.pi/4 - epsilon) ** 2
loss_cp_db = -10 * np.log10(np.maximum(eta_cp, 1e-10))
ax = axes[1]
ax.plot(ar_db, eta_cp, 'b-', linewidth=2, label='$\\eta_p$ (co-polar)')
ax.set_xlabel('Axial Ratio (AR) [dB]', fontsize=12)
ax.set_ylabel('Polarization efficiency $\\eta_p$', fontsize=12)
ax.set_title('Circular Polarization: Effect of Axial Ratio', fontsize=13)
ax3 = ax.twinx()
ax3.plot(ar_db, loss_cp_db, 'r--', linewidth=2, label='Loss [dB]')
ax3.set_ylabel('Polarization loss [dB]', fontsize=12, color='red')
ax3.tick_params(axis='y', labelcolor='red')
ax3.set_ylim(0, 5)
ax.grid(True, alpha=0.3)
ax.set_xlim(0, 10)
ax.set_ylim(0.5, 1.05)
# 3dB ARラインを強調
ax.axvline(x=3, color='green', linestyle=':', linewidth=1.5, label='AR = 3 dB')
plt.tight_layout()
plt.savefig('polarization_mismatch_loss.png', dpi=150, bbox_inches='tight')
plt.show()
左図では、直線偏波の角度ずれに対する偏波効率（青実線）と損失（赤破線）を示しています。角度ずれ15度で損失は約0.3 dB と小さいですが、30度で1.3 dB、45度で3.0 dB と急速に増加します。60度では6.0 dBの損失となり、受信電力が1/4にまで低下します。実用上、直線偏波のアンテナ設置では15度以内の精度が望まれることがこの図から理解できます。右図では、円偏波の軸比（AR）が偏波効率に与える影響を示しています。AR = 0 dB（完全な円偏波）では損失はゼロですが、AR が増加すると損失が発生します。AR = 3 dB の緑破線は実用的な円偏波アンテナの性能目安であり、この時点での損失は約 0.5 dB です。AR = 6 dB では約 1.0 dB、AR = 10 dB では約 1.7 dB の損失が生じます。
偏波多重の可視化
衛星通信での偏波多重の原理を視覚的に示します。
import numpy as np
import matplotlib.pyplot as plt
fig, axes = plt.subplots(1, 2, figsize=(14, 6))
t = np.linspace(0, 4 * np.pi, 500)
# 左: 時間波形（RHCPとLHCPの電界成分）
ax = axes[0]
# RHCP信号（データ: +1）
Ex_R = np.cos(t)
Ey_R = np.sin(t) # delta = -pi/2 -> sin
# LHCP信号（データ: -1）
Ex_L = 0.8 * np.cos(t)
Ey_L = -0.8 * np.sin(t) # delta = +pi/2 -> -sin
# 合成波
Ex_total = Ex_R + Ex_L
Ey_total = Ey_R + Ey_L
ax.plot(t / np.pi, Ex_total, 'b-', linewidth=1.5, label='$E_x$ (combined)', alpha=0.8)
ax.plot(t / np.pi, Ey_total, 'r-', linewidth=1.5, label='$E_y$ (combined)', alpha=0.8)
ax.plot(t / np.pi, Ex_R, 'b--', linewidth=1, label='RHCP $E_x$', alpha=0.4)
ax.plot(t / np.pi, Ey_R, 'r--', linewidth=1, label='RHCP $E_y$', alpha=0.4)
ax.set_xlabel('Time [$\\pi$]', fontsize=12)
ax.set_ylabel('Electric field', fontsize=12)
ax.set_title('Dual Circular Polarization Multiplexing', fontsize=13)
ax.legend(fontsize=9)
ax.grid(True, alpha=0.3)
ax.set_xlim(0, 4)
# 右: 偏波分離後の信号空間
ax = axes[1]
# 受信アンテナで偏波分離した後の出力を計算
# RHCP受信機の出力: RHCP成分を取り出す
n_samples = 50
t_samples = np.linspace(0, 4*np.pi, n_samples)
# RHCP信号 + LHCP干渉 + 雑音
np.random.seed(42)
noise_std = 0.15
# RHCP受信機出力（RHCPはそのまま、LHCPは抑圧）
rhcp_signal = 1.0 + noise_std * np.random.randn(n_samples) # 信号+雑音
lhcp_leakage = 0.05 * np.random.randn(n_samples) # XPD=26dBの漏洩
rhcp_output = rhcp_signal + lhcp_leakage
# LHCP受信機出力
lhcp_signal = -0.8 + noise_std * np.random.randn(n_samples)
rhcp_leakage = 0.05 * np.random.randn(n_samples)
lhcp_output = lhcp_signal + rhcp_leakage
ax.scatter(rhcp_output, np.zeros_like(rhcp_output) + 0.3,
c='blue', alpha=0.6, s=30, label='RHCP channel output')
ax.scatter(lhcp_output, np.zeros_like(lhcp_output) - 0.3,
c='red', alpha=0.6, s=30, label='LHCP channel output')
ax.axvline(x=0, color='k', linestyle=':', linewidth=1)
ax.axvline(x=1.0, color='blue', linestyle='--', linewidth=1, alpha=0.5)
ax.axvline(x=-0.8, color='red', linestyle='--', linewidth=1, alpha=0.5)
ax.set_xlabel('Demodulated signal amplitude', fontsize=12)
ax.set_ylabel('Channel', fontsize=12)
ax.set_yticks([-0.3, 0.3])
ax.set_yticklabels(['LHCP ch.', 'RHCP ch.'])
ax.set_title('Polarization Demultiplexing (XPD $\\approx$ 26 dB)', fontsize=13)
ax.legend(fontsize=10)
ax.grid(True, alpha=0.3)
ax.set_xlim(-2, 2)
ax.set_ylim(-1, 1)
plt.tight_layout()
plt.savefig('polarization_multiplexing.png', dpi=150, bbox_inches='tight')
plt.show()
左図は、RHCP 信号と LHCP 信号が重畳された合成電界の時間波形です。$E_x$（青）と $E_y$（赤）の両成分に2つの偏波信号が混在していますが、偏波分離によって各チャンネルの信号を独立に取り出せます。右図は、偏波分離後の2つのチャンネルの出力を散布図で示しています。RHCP チャンネル（青）の出力は振幅 +1.0 付近に、LHCP チャンネル（赤）の出力は振幅 -0.8 付近にそれぞれ集中しており、XPD が約 26 dB のため交差偏波からの漏洩は非常に小さいことがわかります。偏波の直交性を利用して、同一周波数帯で2つの独立な通信チャンネルを確保できることが視覚的に確認できます。
まとめ
本記事では、電磁波の偏波について理論からアンテナへの応用まで体系的に解説しました。
- 偏波の定義: 電界ベクトルの先端が描く軌跡。振幅比 $E_{0y}/E_{0x}$ と位相差 $\delta = \delta_y – \delta_x$ の2パラメータで偏波状態が決まる
- 偏波の分類: $\delta = 0, \pi$ で直線偏波、$\delta = \pm\pi/2$ かつ $E_{0x} = E_{0y}$ で円偏波（右旋/左旋）、それ以外で楕円偏波。軸比 AR は楕円の扁平率を表す
- ジョーンズベクトル: 完全偏波を2次元複素ベクトルで表現。ジョーンズ行列による偏波変換（回転、位相板など）が行列演算で計算できる
- ストークスパラメータとポアンカレ球: 電力量（$S_0, S_1, S_2, S_3$）で偏波を記述。部分偏波も扱え、ポアンカレ球上の幾何学で偏波間の関係が直感的に理解できる
- 偏波不整合損失: $\eta_p = |\bm{J}_a^H \bm{J}_w|^2$ で定量化。直線偏波の45度ずれで 3 dB 損失、右旋-左旋間は完全不整合（$\eta_p = 0$）
- アンテナへの応用: 円偏波アンテナ（パッチ、シーケンシャル回転）、衛星通信での偏波多重（容量2倍）、GPS での反射波抑制
偏波は電磁波の基本的な性質でありながら、アンテナ設計、レーダー信号処理、光通信など極めて多くの工学分野に直結する重要な概念です。特に衛星通信では、偏波の選択と管理がリンクバジェットに直接影響するため、偏波不整合損失の正確な見積りが不可欠です。
次のステップとして、以下の記事も参考にしてください。