パルス圧縮のソースを表示

{{otheruses|信号処理技術|パルスレーザー増幅技術|チャープパルス増幅}}
{{翻訳直後|[[:en:Special:Redirect/revision/695804642|en:Pulse compression]]|date=2016年10月}}

'''パルス圧縮'''（パルスあっしゅく{{lang-en-short|Pulse compression}}）とは、[[レーダー]]・[[ソナー]]・[[超音波検査]]などに広く応用される、距離[[分解能]]および[[SN比|S/N比]]向上のための[[信号処理]]技術である。送信パルスを[[変調方式|変調]]し、受信信号と送信パルスとの[[相互相関関数|相関]]をとることにより達成される<ref>J. R. Klauder, A. C, Price, S. Darlington and W. J. Albersheim, ‘The Theory and Design of Chirp Radars,” Bell System Technical Journal 39, 745 (1960).</ref>。

== 単純パルス ==

=== 信号 ===
パルスレーダーが送信できる最も単純な信号は振幅 {{Mvar|A}} で{{仮リンク|搬送波周波数|en|Carrier frequency}}が {{Math|''f''<sub>0</sub>}} の正弦波を幅 {{Mvar|T}} の[[矩形関数]]で切り取った正弦波パルスである。このパルスは周期的に送信されるが、そのことはこの項における主題に関係がないので、単パルス {{Mvar|s}} のみを考えるものとする。パルスが時刻 {{Math|1=''t'' = 0}} から開始するものと仮定すれば、この信号は[[複素数]]表記を用いて次のように書ける。
: <math>
s(t) =  \begin{cases}
	A \exp(2 i \pi f_0 t) &\text{if} \; 0 \leq t < T \\
	0 &\text{otherwise}
\end{cases}</math>

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

ここで整合フィルタとは、受信信号と送信信号との[[相互相関関数]]を計算することである。これは入力信号に送信信号の[[複素共役|共役]]をとって時間反転を施したものと[[畳み込み|畳み込む]]ことにより達成できる。この演算はソフトウェアによって行うことも、ハードウェアによって行うこともできる。この相互相関関数を {{Math|&lang;''s'', ''r''&rang;}} と書くことにすると、以下を得る。
: <math>\langle s,r \rangle (t) = \int_{t'\,=\,0}^{+\infty} s^\star(t')r(t+t') \mathrm dt'</math>
反射された信号が時刻 {{Math|''t''<sub>''r''</sub>}} に受信機に戻ってきたとし、減衰因子を {{Mvar|K}} とすると次を得る。
: <math>r(t)=
\begin{cases}
	K A \exp\{2 i \pi f_0 (t\,-\,t_r)\} +B(t) &\mbox{if} \; t_r \leq t < t_r+T \\
	B(t) &\mbox{otherwise}
\end{cases}</math>
送信信号は既知であるから、次を得る。
: <math>\langle s,r \rangle (t) = KA^2\Lambda\left (\frac{t-t_r}{T} \right)\exp\{2 i \pi f_0 (t\,-\,t_r)\} + B'(t)</math>
ここで、 {{Math|''B''&prime;(''t'')}} は送信信号とノイズとの相互相関関数である。 関数 {{Math|&Lambda;}} は[[三角形関数]]で、 {{Math|[&minus;&infin;, &minus;{{sfrac|1|2}}]∪[{{sfrac|1|2}}, &infin;]}} における値は 0 で {{Math|[&minus;{{sfrac|1|2}}, 0]}} において線形に増加して最大値 1 に達し、 {{Math|[0, {{sfrac|1|2}}]}} においては線形に減少して 0 にもどる。この節の末尾の図にサンプル信号（赤線）の相互相関関数の形を示す。サンプル信号は長さ {{Math|1=''T'' = 1}} 秒に切り取られた、振幅 {{Math|1}}、周波数 {{Math|1=''f''<sub>0</sub> = {{val|10|u=Hz}}}} の正弦波である。二つのエコー（青線）が3秒および5秒遅れで、それぞれ振幅が送信信号の 0.5 倍および 0.3 倍になって戻ってきている。 これらは例示のために選んだランダムな値である。信号は実数であるから、相互相関関数には追加で {{Frac|1|2}} 倍がかかっている。

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

時間 {{Mvar|T}} の間に波が伝播する距離は {{Mvar|cT}} （{{Mvar|c}} はその媒質中における波の速さ）であり、これが往復時間に対応することから次を得る。
{| class="wikitable" style="margin: auto auto 10px;"
! 結果 1
|-
| 正弦波パルスの距離分解能は {{Math|{{sfrac|1|2}}''cT''}} である。ただし {{Mvar|T}} はパルス長、{{Mvar|c}} は波の速さとする。

結論: 分解能を挙げるためには、パルス長を減らさなければならない。
|}

{| style="margin: 1em auto; margin-bottom: 10px;" border="0"
|+ '''例 （単純パルス）: 送信信号は赤（搬送波 10 ヘルツ、振幅 1、パルス長 1 秒）で二つのエコーは青。'''
!整合フィルタリング前
! 整合フィルタリング後
|-
|[[ファイル:Two_targets_before.jpg|サムネイル|300x300ピクセル|標的が十分離れていれば…]] 
| [[ファイル:Two_targets_resoluted.jpg|サムネイル|300x300ピクセル|…エコーを弁別することが可能である。]]
|-
| [[ファイル:Two_targets_below_resolution_before.jpg|サムネイル|300x300ピクセル|標的が互いに近すぎると…]]
| [[ファイル:Two_targets_below_resolution.jpg|サムネイル|300x300ピクセル|…エコーは混りあってしまう。]]
|}

=== 信号の送信に必要なエネルギー ===
送信パルスの瞬時電力は {{Math|1=''P''(''t'') = {{mabs|''s''}}<sup>2</sup>(''t'')}} で与えられる。このパルスの送信に必要なエネルギーは以下のように書ける。
: <math>E = \int_0^T P(t)\mathrm dt = A^2 T</math>
同様に、 受信パルスのエネルギーは {{Math|1=''E''<sub>''r''</sub> = ''K''<sup>2</sup>''A''<sup>2</sup>T}} である。{{Mvar|&sigma;}} をノイズの標準偏差とすれば、受信信号の[[SN比|S/N比]]（SNR）は以下の通りとなる。
: <math>SNR = \frac{E_r}{\sigma^{2}} = \frac{K^2 A^2 T}{\sigma^{2}}</math>
SNR は、他のパラメータを一定に保てばパルス長 {{Mvar|T}} に比例する。このことから、次のトレードオフが帰結する。 {{Mvar|T}} が長くなれば SNR は改善するが、分解能が犠牲になる。また、逆も同様である。

== 線形周波数変調（チャープ）によるパルス圧縮 ==
=== 基本的原理 ===
解像度を損うことなく十分な SNR を得られるだけの大きなパルスを得るためにはどのようにすればよいだろうか？ここにパルス圧縮が登場する。基本的原理は次のようなものである。

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

[[レーダー]]や[[ソナー]]への応用においては、パルス圧縮には線形[[チャープ信号|チャープ]]が最も一般的に用いられる。パルスは有限の長さを持ち、振幅は[[矩形関数]]とする。 もし送信信号の長さが {{Mvar|T}} で時刻 {{Math|1=''t'' = 0}} に開始し、周波数帯 {{Math|&Delta;''f''}} だけ搬送波周波数 {{Math|''f''<sub>0</sub>}} を中心として掃引するものとすると、次のように書ける。
: <math>
s_c(t) = 
\begin{cases}
	A \exp\left[ i 2 \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}</math>
上の定義から、チャープ信号の位相（つまり、複素べき指数の偏角）は二次式となる。
: <math>\phi(t) = 2\pi \left( \left( f_0 \,-\, \frac{\Delta f}{2}\right) t \, + \, \frac{\Delta f}{2T}t^2 \, \right) </math>
よって、瞬時周波数は（定義より）次のように書ける。
: <math>f(t) = \frac{1}{2\pi}\left[\frac{\mathrm d\phi}{\mathrm dt}\right ]_t = f_0-\frac{\Delta f}{2}+\frac{\Delta f}{T}t</math>
これは、意図どおり {{Math|1=''t'' = 0}} における {{Math|''f''<sub>0</sub> &minus; {{sfrac|&Delta;''f''|2}}}} から {{Math|1=''t'' = ''T''}} における {{Math|''f''<sub>0</sub> + {{sfrac|&Delta;''f''|2}}}} へと線形な坂を登っている。

位相と周波数の関係式を逆に用いることも多い。所望のチャープ周波数を {{Math|''f''(''t'')}} とすると、チャープ位相は積分を用いて次のように書ける。
: <math>\phi(t) = 2 \pi \int_0^t f(u)\,\mathrm du </math>

=== 送信信号と受信信号の相互相関関数 ===
単純パルスの場合と同様、送信信号と受信信号の相互相関関数を計算する。簡略化のため、チャープを上述のものではなく次の別の形とすることにしよう（最終的な結果は同じである）。
: <math>
s_{c'}(t) = \begin{cases}
	A \exp\left[ 2 i \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}</math>
相互相関関数は、（減衰因子 {{Mvar|K}} を除けば）{{Math|''s''<sub>''c''&prime;</sub>}} の自己相関関数と等しいから、以降は次の自己相関関数について考える。
: <math>\langle s_{c'},s_{c'}\rangle (t) = \int_{-\infty}^{+\infty}s_{c'}^\star(-t')s_{c'}(t-t')\mathrm dt'</math>
{{Math|''s''<sub>''c''&prime;</sub>}} の自己相関関数は次のように書けることが知られている<ref>Achim Hein, ''Processing of SAR Data: Fundamentals, Signal Processing, Interferometry'', Springer, 2004, ISBN 3-540-05043-4, pages 38 to 44. </ref>。
: <math>\langle s_{c'}, s_{c'}\rangle(t) = A^2T \Lambda \left(\frac{t}{T} \right) \mathrm{sinc} \left[ \pi \Delta f t \Lambda \left( \frac{t}{T}\right) \right] \exp\left( 2 i \pi f_0 t \right) </math>
{{Math|''s''<sub>''c''&prime;</sub>}} の自己相関関数は {{Math|1=''t'' = 0}} において最大値をとる。{{Math|1=''t'' = 0}} の近傍では、この関数は [[Sinc関数|{{Math|sinc}} 関数]]項のように振る舞う。この {{Math|sinc}} 関数の {{Val|−3|u=dB}} 時間幅はおおよそ {{Math|1=''T&prime;'' = {{sfrac|1|&Delta;''f''}}}} に等しい。整合フィルタリング後には全てがパルス長 {{Mvar|T&prime;}} の単純パルスであった場合と同様になる。通常用いられる {{Math|&Delta;''f''}} の値の場合、 {{Mvar|T&prime;}} は {{Mvar|T}} よりも短く、したがってこれを「パルス圧縮」と呼ぶ。

{{Math|sinc}} 関数は厄介な{{仮リンク|サイドローブ|en|Sidelobe}}を持っているので、通常の運用では適当な[[窓関数]]（[[窓関数#ハミング窓|ハミング窓]]や[[窓関数#ハン窓|ハン窓]]など）を用いてフィルタリングを行う。実用上、これは適応フィルタリングと同時にフィルタに基準チャープを積算することで行うことができる。結果として最大値は若干減るが、サイドローブがフィルタリングされることのほうが重要である。
{| class="wikitable" style="margin: auto auto 10px;"
! 結果 2
|-
| あるバンド幅 {{Math|&Delta;''f''}} の線形周波数変調により達成できる距離分解能は {{Math|{{sfrac|''c''|2&Delta;''f''}}}} である。ただし、 {{Mvar|c}} は波の速さ。
|}
{| class="wikitable" style="margin: auto auto 10px;"
!定義
|-
| 比 {{Math|1={{sfrac|''T''|''T''&prime;}} = ''T''&Delta;''f''}} をパルス圧縮比と呼ぶ。一般にこの値は {{Math|1}} よりも大きい（通常、 {{Math|20}} から {{Math|30}} 程度）。
|}

{| style="margin: 1em auto; margin-bottom: 10px;" border="0"
|+ 例（チャープパルス）: 送信信号は赤線（搬送波周波数 10 ヘルツ、変調周波数 16 ヘルツ、振幅 1、 パルス長 1 秒) 、エコーは青線
|[[ファイル:Chirp_before.jpg|サムネイル|300x300ピクセル|整合フィルタリング前]] 
| [[ファイル:Chirp_compr.jpg|サムネイル|300x300ピクセル|整合フィルタリング後: エコーが時間的に短縮されている。]]
|}

=== パルス圧縮による SNR の改善 ===
パルス圧縮により、信号のエネルギー総量が変更されることはない。しかし、信号のエネルギーは幅およそ {{Math|''T&prime;'' &cong; {{sfrac|1|&Delta;''f''}}}} の {{Math|sinc}} 関数のメインローブに詰め込まれることになる。{{Mvar|P}} を圧縮前の信号パワー、{{Mvar|P&prime;}} を圧縮後の信号パワーとすると、次の関係式を得る。
: <math>P\times T = P' \times T' </math>
したがって、次のように変形できる。
: <math>P'= P\times \frac{T}{T'} </math>
一方で、ノイズは送信周波数と全く相関しない（完全ランダム）のでノイズのパワーは変化しない。結果として、次のように言える。 <br>
{| class="wikitable" style="margin: auto auto 10px;"
! 結果 3
|-
| パルス圧縮後、受信信号のパワーは {{Math|''T''&Delta;''f''}} だけ増幅されたものと考えてよい。この追加ゲインは[[レーダー|レーダー方程式]]に代入できる。
|}

{| style="margin: 1em auto; margin-bottom: 10px;" border="0"
|+ 例: 上図と同じ信号に加法的ガウシアンホワイトノイズ ({{Math|1=&sigma; = 0.5}}) を重畳した信号波形
|[[ファイル:Chirp_noise.jpg|サムネイル|300x300ピクセル|整合フィルタリング前: 信号はノイズに埋もれて見えない。]] 
| [[ファイル:Chirp_compr_noise.jpg|サムネイル|300x300ピクセル|整合フィルタリング後: エコーが見えるようになった。]]
|}

== 位相符号化によるパルス圧縮 ==
他にも信号を変調する方法は存在する。[[位相変調]]はよく用いられる手法で、この場合はパルスを {{Mvar|N}} 個の長さ {{Math|{{sfrac|''T''|''N''}}}} のタイムスロットに分割して、それぞれの開始時の位相をあらかじめ定めた規約に基いて決定する。たとえば、位相を変えないタイムスロット（そのスロットでは信号はそのまま送信される）と、信号を {{Mvar|&pi;}} だけシフト（信号の符号反転と等価）するタイムスロットを並べる。位相シフト {{Math|{0, ''&pi;''} }} の配列は{{仮リンク|バーカー符号|en|Barker code}}として知られる手法に則って選ぶ。二種類の位相だけでなくより多くの位相を用いて符号化することもある（多相符号化）。線形チャープの場合と同様、パルス圧縮は相互相関を通じて達成される。

バーカー符号の利点<ref>J.-P. Hardange, P. Lacomme, J.-C. Marchais, ''Radars aéroportés et spatiaux'', Masson, Paris, 1995, ISBN 2-225-84802-5, page 104. </ref>はその単純さ（上述のとおり、 {{Mvar|&pi;}} 位相シフトは符号反転で済む）にあるが、パルス圧縮比はチャープの場合よりも低くなり、また[[ドップラー効果]]による {{Math|{{sfrac|1|''T''}}}} を超える周波数変化に対して非常に鋭敏になる。

== 出典 ==
{{脚注ヘルプ}}
{{reflist}}

== 関連文献 ==
* Nadav Levanon, and Eli Mozeson. Radar signals. Wiley. com, 2004.
* Hao He, Jian Li, and Petre Stoica. [http://www.sal.ufl.edu/book/ Waveform design for active sensing systems: a computational approach]. Cambridge University Press, 2012.
* M. Soltanalian. [http://theses.eurasip.org/theses/573/signal-design-for-active-sensing-and/download/ Signal Design for Active Sensing and Communications]. Uppsala Dissertations from the Faculty of Science and Technology (printed by Elanders Sverige AB), 2014.
* Solomon W. Golomb, and Guang Gong. [http://www.cambridge.org/us/academic/subjects/computer-science/cryptography-cryptology-and-coding/signal-design-good-correlation-wireless-communication-cryptography-and-radar Signal design for good correlation: for wireless communication, cryptography, and radar]. Cambridge University Press, 2005.
* Fulvio Gini, Antonio De Maio, and Lee Patton, eds. Waveform design and diversity for advanced radar systems. Institution of engineering and technology, 2012.
* John J. Benedetto, Ioannis Konstantinidis, and Muralidhar Rangaswamy. "[http://ieeexplore.ieee.org/xpl/login.jsp?tp=&arnumber=4775877 Phase-coded waveforms and their design]." IEEE Signal Processing Magazine, 26.1 (2009): 22-31.
* Ducoff, Michael R., and Byron W. Tietjen. "Pulse compression radar." Radar Handbook (2008): 8-3.

== 関連項目 ==
* [[スペクトラム拡散]]
* {{仮リンク|チャープ圧縮|en|Chirp compression}}

{{DEFAULTSORT:はるすあつしゆく}}
[[Category:信号処理]]