パルス圧縮

パルス圧縮（パルスあっしゅく英: Pulse compression）とは、レーダー・ソナー・超音波検査などに広く応用される、距離分解能およびS/N比向上のための信号処理技術である。送信パルスを変調し、受信信号と送信パルスとの相関をとることにより達成される^[1]。

単純パルス[編集]

信号[編集]

パルスレーダーが送信できる最も単純な信号は振幅 $A$ で搬送波周波数（英語版）が $f 0$ の正弦波を幅 $T$ の矩形関数で切り取った正弦波パルスである。このパルスは周期的に送信されるが、そのことはこの項における主題に関係がないので、単パルス $s$ のみを考えるものとする。パルスが時刻 $t = 0$ から開始するものと仮定すれば、この信号は複素数表記を用いて次のように書ける。

s(t)={\begin{cases}A\exp(2i\pi f_{0}t)&{\text{if}}\;0\leq t<T\\0&{\text{otherwise}}\end{cases}}

距離分解能[編集]

このような信号で得ることのできる距離分解能を決定しよう。帰還信号 $r (t)$ は、送信信号のコピーに減衰と時間遅れを加えたものである（現実にはドップラー効果も作用するが、ここでは重要ではない）。入力信号には実部にも虚部にもノイズが乗っていると考え、これを白色かつガウス的なノイズとして（現実でも一般的にはこれが成り立つ）、 $B (t)$ と書くことにする。入力信号を検波するには、整合フィルタ（英語版）が用いられることが一般的である。この方法は加法的ホワイトガウシアンノイズ（英語版）中に埋もれた既知信号を検波するのに最適である。

ここで整合フィルタとは、受信信号と送信信号との相互相関関数を計算することである。これは入力信号に送信信号の共役をとって時間反転を施したものと畳み込むことにより達成できる。この演算はソフトウェアによって行うことも、ハードウェアによって行うこともできる。この相互相関関数を $⟨ s, r ⟩$ と書くことにすると、以下を得る。

\langle s,r\rangle (t)=\int _{t'\,=\,0}^{+\infty }s^{\star }(t')r(t+t')\mathrm {d} t'

反射された信号が時刻 $t r$ に受信機に戻ってきたとし、減衰因子を $K$ とすると次を得る。

r(t)={\begin{cases}KA\exp\{2i\pi f_{0}(t\,-\,t_{r})\}+B(t)&{\mbox{if}}\;t_{r}\leq t<t_{r}+T\\B(t)&{\mbox{otherwise}}\end{cases}}

送信信号は既知であるから、次を得る。

\langle s,r\rangle (t)=KA^{2}\Lambda \left({\frac {t-t_{r}}{T}}\right)\exp\{2i\pi f_{0}(t\,-\,t_{r})\}+B'(t)

ここで、 $B'(t)$ は送信信号とノイズとの相互相関関数である。関数 $Λ$ は三角形関数で、 $[−∞, −.mw-parser-output .sfrac{white-space:nowrap}.mw-parser-output .sfrac.tion,.mw-parser-output .sfrac .tion{display:inline-block;vertical-align:-0.5em;font-size:85%;text-align:center}.mw-parser-output .sfrac .num,.mw-parser-output .sfrac .den{display:block;line-height:1em;margin:0 0.1em}.mw-parser-output .sfrac .den{border-top:1px solid}.mw-parser-output .sr-only{border:0;clip:rect(0,0,0,0);height:1px;margin:-1px;overflow:hidden;padding:0;position:absolute;width:1px}1/2]∪[1/2, ∞]$ における値は 0 で $[- 1 / 2, 0]$ において線形に増加して最大値 1 に達し、 $[0, 1 / 2]$ においては線形に減少して 0 にもどる。この節の末尾の図にサンプル信号（赤線）の相互相関関数の形を示す。サンプル信号は長さ $T = 1$ 秒に切り取られた、振幅 $1$ 、周波数 $f 0 = 10 Hz$ の正弦波である。二つのエコー（青線）が3秒および5秒遅れで、それぞれ振幅が送信信号の 0.5 倍および 0.3 倍になって戻ってきている。これらは例示のために選んだランダムな値である。信号は実数であるから、相互相関関数には追加で 1⁄2 倍がかかっている。

もしふたつのパルスが（ほぼ）同時に戻ってきた場合、相互相関関数はふたつの素相互相関関数の和となる。この中から、ひとつのパルスの「三角型」の包絡線を別のパルスのものと弁別するためには、二つのパルスの最大値が分離できなければならないので、パルスの到着時刻が最低でも $T$ だけ離れていなければならないことは明らかである。この条件が満たされない場合、二つの三角形は互いに混りあい分離不可能となる。

時間 $T$ の間に波が伝播する距離は $cT$ （ $c$ はその媒質中における波の速さ）であり、これが往復時間に対応することから次を得る。

結果 1
正弦波パルスの距離分解能は $1 / 2 cT$ である。ただし $T$ はパルス長、 $c$ は波の速さとする。結論: 分解能を挙げるためには、パルス長を減らさなければならない。

**例（単純パルス）: 送信信号は赤（搬送波 10 ヘルツ、振幅 1、パルス長 1 秒）で二つのエコーは青。**
整合フィルタリング前	整合フィルタリング後
標的が十分離れていれば…	…エコーを弁別することが可能である。
標的が互いに近すぎると…	…エコーは混りあってしまう。

信号の送信に必要なエネルギー[編集]

送信パルスの瞬時電力は $P (t) = | s | 2 (t)$ で与えられる。このパルスの送信に必要なエネルギーは以下のように書ける。

E=\int _{0}^{T}P(t)\mathrm {d} t=A^{2}T

同様に、受信パルスのエネルギーは $E r = K 2 A 2 T$ である。 $σ$ をノイズの標準偏差とすれば、受信信号のS/N比（SNR）は以下の通りとなる。

SNR={\frac {E_{r}}{\sigma ^{2}}}={\frac {K^{2}A^{2}T}{\sigma ^{2}}}

SNR は、他のパラメータを一定に保てばパルス長 $T$ に比例する。このことから、次のトレードオフが帰結する。 $T$ が長くなれば SNR は改善するが、分解能が犠牲になる。また、逆も同様である。

線形周波数変調（チャープ）によるパルス圧縮[編集]

基本的原理[編集]

解像度を損うことなく十分な SNR を得られるだけの大きなパルスを得るためにはどのようにすればよいだろうか？ここにパルス圧縮が登場する。基本的原理は次のようなものである。

送信信号はエネルギー予算が正しくなよう十分に長く送信する。
この信号は整合フィルタリング後の相互相関信号の幅が標準的な上述した正弦波パルスよりも小さくなるように設計する（このことにパルス圧縮の名前は由来する）。

レーダーやソナーへの応用においては、パルス圧縮には線形チャープが最も一般的に用いられる。パルスは有限の長さを持ち、振幅は矩形関数とする。もし送信信号の長さが $T$ で時刻 $t = 0$ に開始し、周波数帯 $Δ f$ だけ搬送波周波数 $f 0$ を中心として掃引するものとすると、次のように書ける。

s_{c}(t)={\begin{cases}A\exp \left[i2\pi \left\{\left(f_{0}\,-\,{\frac {\Delta f}{2}}\right)t\,+\,{\frac {\Delta f}{2T}}t^{2}\,\right\}\right]&{\mbox{if}}\;0\leq t<T\\0&{\mbox{otherwise}}\end{cases}}

上の定義から、チャープ信号の位相（つまり、複素べき指数の偏角）は二次式となる。

\phi (t)=2\pi \left(\left(f_{0}\,-\,{\frac {\Delta f}{2}}\right)t\,+\,{\frac {\Delta f}{2T}}t^{2}\,\right)

よって、瞬時周波数は（定義より）次のように書ける。

f(t)={\frac {1}{2\pi }}\left[{\frac {\mathrm {d} \phi }{\mathrm {d} t}}\right]_{t}=f_{0}-{\frac {\Delta f}{2}}+{\frac {\Delta f}{T}}t

これは、意図どおり $t = 0$ における $f 0 - Δ f / 2$ から $t = T$ における $f 0 + Δ f / 2$ へと線形な坂を登っている。

位相と周波数の関係式を逆に用いることも多い。所望のチャープ周波数を $f (t)$ とすると、チャープ位相は積分を用いて次のように書ける。

\phi (t)=2\pi \int _{0}^{t}f(u)\,\mathrm {d} u

送信信号と受信信号の相互相関関数[編集]

単純パルスの場合と同様、送信信号と受信信号の相互相関関数を計算する。簡略化のため、チャープを上述のものではなく次の別の形とすることにしよう（最終的な結果は同じである）。

s_{c'}(t)={\begin{cases}A\exp \left[2i\pi \left(f_{0}\,+\,{\frac {\Delta f}{2T}}t\right)t\right]&{\mbox{if}}\;-{\frac {T}{2}}\leq t<{\frac {T}{2}}\\0&{\mbox{otherwise}}\end{cases}}

相互相関関数は、（減衰因子 $K$ を除けば） $s c'$ の自己相関関数と等しいから、以降は次の自己相関関数について考える。

\langle s_{c'},s_{c'}\rangle (t)=\int _{-\infty }^{+\infty }s_{c'}^{\star }(-t')s_{c'}(t-t')\mathrm {d} t'

$s c'$ の自己相関関数は次のように書けることが知られている^[2]。

\langle s_{c'},s_{c'}\rangle (t)=A^{2}T\Lambda \left({\frac {t}{T}}\right)\mathrm {sinc} \left[\pi \Delta ft\Lambda \left({\frac {t}{T}}\right)\right]\exp \left(2i\pi f_{0}t\right)

$s c'$ の自己相関関数は $t = 0$ において最大値をとる。 $t = 0$ の近傍では、この関数は $sinc$ 関数項のように振る舞う。この $sinc$ 関数の −3 dB 時間幅はおおよそ $T' = 1 / Δ f$ に等しい。整合フィルタリング後には全てがパルス長 $T'$ の単純パルスであった場合と同様になる。通常用いられる $Δ f$ の値の場合、 $T'$ は $T$ よりも短く、したがってこれを「パルス圧縮」と呼ぶ。

$sinc$ 関数は厄介なサイドローブ（英語版）を持っているので、通常の運用では適当な窓関数（ハミング窓やハン窓など）を用いてフィルタリングを行う。実用上、これは適応フィルタリングと同時にフィルタに基準チャープを積算することで行うことができる。結果として最大値は若干減るが、サイドローブがフィルタリングされることのほうが重要である。

結果 2
あるバンド幅 $Δ f$ の線形周波数変調により達成できる距離分解能は $c / 2Δ f$ である。ただし、 $c$ は波の速さ。

定義
比 $T / T' = T Δ f$ をパルス圧縮比と呼ぶ。一般にこの値は $1$ よりも大きい（通常、 $20$ から $30$ 程度）。

例（チャープパルス）: 送信信号は赤線（搬送波周波数 10 ヘルツ、変調周波数 16 ヘルツ、振幅 1、パルス長 1 秒) 、エコーは青線
整合フィルタリング前	整合フィルタリング後: エコーが時間的に短縮されている。