(57)【特許請求の範囲】
【請求項１】
２つの信号間の時間シフトを求めるための装置であって、
各チャネルが対応する周波数に関連付けられる複数のチャネルであって、
チャネルに対応する周波数において第１の変調信号によって前記２つの信号の一方を変調すると共に、チャネルに対応する周波数において第２の変調信号によって前記２つの信号の前記一方を変調して、被変調信号の対を生成し、前記第１の変調信号及び前記第２の変調信号は、同位相でもなく逆位相でもない（out of phase and out of anti-phase）固定された位相関係を有する、第１の変調手段、
チャネルに対応する周波数において前記２つの信号の他方を変調して、被変調信号を生成する、第２の変調手段、並びに
さまざまな時間シフトにおいて、前記第２の変調手段によって出力された前記被変調信号と、前記第１の変調手段によって出力された前記被変調信号の対とを比較して、比較関数の対を生成する、比較手段、
を有する、複数のチャネルと、
前記複数のチャネルからの前記複数の比較関数を結合して合成比較関数を生成する、結合手段と、
前記合成比較関数から時間シフトを抽出する、抽出手段と、
を備える、２つの信号間の時間シフトを求めるための装置。
【請求項２】
前記結合手段によって、前記合成比較関数のエネルギーが前記複数のチャネルにおいて生成された前記複数の比較関数のエネルギーと比較して時間領域において集中化されるように、前記複数のチャネルに対応する各周波数が配置される、請求項１に記載の装置。
【請求項３】
前記チャネルのそれぞれは、すべてのチャネルにおいて実質的に同じタイミングの主ピークの大きさと、すべてのチャネルにおいてそれぞれ異なるタイミングの副ピークとを有する複素関数に対応する比較関数の対を生成するように動作可能であり、
それによって、前記結合手段が前記複数の比較関数の対を結合すると、前記主ピークの大きさは、前記副ピークの大きさと比較して増強される、請求項１又は２に記載の装置。
【請求項４】
前記比較手段は、少なくとも１つの相互相関器を備える、請求項１〜３のいずれか一項に記載の装置。
【請求項５】
前記第１の変調信号及び前記第２の変調信号は、互いの位相が直交（in phase quadrature with each other）している、請求項１〜４のいずれか一項に記載の装置。
【請求項６】
前記チャネルのそれぞれは、該チャネルに対応する周波数のローカル発振器を備える、請求項１〜５のいずれか一項に記載の装置。
【請求項７】
すべての変調動作がコヒーレントである、請求項１〜６のいずれか一項に記載の装置。
【請求項８】
前記結合手段は、
前記チャネルのそれぞれから前記比較関数の対の一方を受け取って、該受け取った比較関数に対して加算演算を実行する、第１の加算手段と、
前記チャネルのそれぞれから前記比較関数の対の他方を受け取って、該受け取った比較関数に対して加算演算を実行する、第２の加算手段と、
前記第１の加算手段の出力及び前記第２の加算手段の出力を２乗して、その結果に対して加算演算を実行して前記合成比較関数を生成するように動作可能である、関数結合器と、
を備える、請求項１〜７のいずれか一項に記載の装置。
【請求項９】
前記第１の加算手段及び前記第２の加算手段の少なくとも一方は、重み係数のセットに従って重み付け加算を実行するように動作可能である、請求項８に記載の装置。
【請求項１０】
前記重み係数のセットは、ｐ乗された２項係数を含み、ｐは０から１までである、請求項９に記載の装置。
【請求項１１】
ベースバンドに関連付けられる前記重み係数は、対応する２項係数の２分の１である、請求項１０に記載の装置。
【請求項１２】
ｐの値を変化させるように動作可能な制御手段をさらに備える、請求項１０又は１１に記載の装置。
【請求項１３】
前記抽出手段は、前記合成比較関数の最大ピークのタイミングを求めるように動作可能である、請求項１〜１２のいずれか一項に記載の装置。
【請求項１４】
前記２つの信号は帯域通過信号である、請求項１〜１３のいずれか一項に記載の装置。
【請求項１５】
２つの信号間のシフトを求める方法であって、
複数の周波数のそれぞれについて、前記２つの信号が該周波数に変換されたバージョン間の複素相互相関を実行することと、
そのようにして導出された複数の複素相互相関関数を結合することによって、前記シフトを示す位置においてピークを有する合成相互相関関数を導出することと
含む、２つの信号間のシフトを求める方法。
【請求項１６】
２つの信号間のシフトを求める方法であって、
複数の周波数のそれぞれについて、
前記２つの信号が該周波数に変換されたバージョンの第１の成分を比較することによって、第１の比較関数を導出することと、
前記２つの信号の前記バージョンの第２の成分を比較することによって、第２の比較関数を導出することと、
前記複数の比較関数を結合して、前記シフトを示す位置においてピークを有する合成比較関数を導出することと、
を含み、
前記第２の成分は、前記第１の成分と同位相でもなく逆位相でもない（out of phase or antiphase）、２つの信号間のシフトを求める方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、相互相関関数の圧縮バージョン（compacted version）を利用することによって、帯域通過信号とその時間シフトされたレプリカとの間の時間シフトを求めるための方法及び装置に関する。この方法は、他を除外するものではないが特に、物体の検出及び測距のタスクを実行するシステムにおける時間遅延の測定に適用可能である。
【背景技術】
【０００２】
多くの計測システムは、帯域通過信号という特定の形態の信号を利用する。帯域通過信号とは、低域周波数ｆL及び高域周波数ｆHという２つの特定の周波数によって定められる一定の周波数帯域の外側において、周波数スペクトルが実質的にゼロになる信号である。帯域通過信号は、その中央帯域（又は中心）周波数ｆC＝（ｆH＋ｆL）／２に対する信号帯域幅（ｆH−ｆL）の比として定義されるその比帯域幅の値に応じて、狭帯域、広帯域、又は超広帯域として分類することができる。狭帯域信号は小さな（例えば、１パーセントの）比帯域幅を有するが、超広帯域信号は０．２５を超える比帯域幅によって特徴付けられる。通例、帯域通過信号が狭帯域でも超広帯域でもないとき、それらの帯域通過信号は広帯域信号と呼ばれる。
【０００３】
帯域通過信号を活用する１つの重要な用途は、自動車におけるマイクロ波センサを用いた障害物の検出及び測距の用途である。図１は、マイクロ波自動車レーダの簡略化したブロック図である。
【０００４】
波形ジェネレータＷＧＴによって供給された適切な波形ＭＷは、マイクロ波発振器ＯＳＣによって生成された搬送波Ｆ０を変調するために、変調器ＭＯＤにおいて使用される。変調器ＭＯＤの出力信号ＳＴは、増幅器／ドライバＤＲＶで増幅され、適切に変調された形態ＳＳで送信素子ＴＥＬに印加される。その結果、マイクロ波インターロゲート信号が障害物ＯＢＪに向けて送信される。この送信信号は、連続的なものとすることもできるし、追加のパルス内変調が施されたパルス列若しくはパルス内変調が施されないパルス列の形態とすることもできる。
【０００５】
障害物ＯＢＪから反射された信号は、受信素子ＲＥＬによって捕捉され、変調された形態ＲＲで増幅され、信号調節ユニットＳＣＵにおいて適切に成形される。ユニットＳＣＵの出力信号ＳＲは、発振器ＯＳＣによって生成された搬送波Ｆ０によって同様に駆動されるダウンコンバータＤＣＲにおいてベースバンドへ変換される。
【０００６】
受信信号のベースバンドバージョンＢＲは、適切な相関器ＣＯＲにおいて、変調波形ＭＷと併せて処理される。当業者には既知であるように、信号ＢＲとＭＷとの間の相互相関関数ＲＥを使用して、送信マイクロ波信号と反射マイクロ波信号との間の時間遅延を求めることができる。したがって、障害物までの距離を求めることができる。
【０００７】
より具体的には、時間遅延は、得られた相互相関関数ＲＥの時間領域におけるピークの位置から求められる。反射信号が送信信号の形態を保持していれば、相互相関関数はその信号の自己相関関数と同じ形状を有する。当業者には既知であるように、インターロゲート信号が占有する帯域幅が広くなると、相関ピークの幅は狭くなる。したがって、時間遅延の測定の精度を改善するためには、広帯域信号が用いられる。
【０００８】
図１における信号ＳＴとＳＲのような送信帯域通過信号と反射帯域通過信号との間の相互相関関数の形状を解析すると、相互相関関数の形状は、搬送波周波数におけるマイクロ波パルスの形態をとり、このパルスの包絡線は、ベースバンドにおいて求められる相互相関関数を複製することが示される。
【０００９】
図２ａは、送信帯域通過信号の電力スペクトル密度を示す。このスペクトルは、搬送波周波数ｆCにおいて帯域幅２Ｂを占有する。図２ｂは、ベースバンドへ変換されたこの信号のスペクトルを示す。図２ｃは、この帯域通過信号の相関関数を示し、図２ｄは、この帯域通過信号のベースバンドバージョンの相関関数を示す。相関包絡線内で観測される搬送波周波数ｆCのサイクル数は、比ｆC／Ｂに比例する。
【００１０】
帯域通過信号とそのベースバンドバージョンのそれぞれの相関関数の形状の詳細な解析は、S. Sugasawa著：「Time Difference Measurement of Ultrasonic Pulses Using Cross-Correlation Function between Analytic Signals. Jpn. J. Appl. Phys. Vol. 41 (2002) pp. 3299-3307」に示されている。
【００１１】
この解析によって、ベースバンド相関関数のピークは放物線（１−γＢ2τ2）によって近似できるのに対して、帯域通過相関関数のピークは別の放物線［１−γ（ｆC2＋Ｂ2）τ2］によって近似できることが示される。ここで、γは定数である。図から分かるように、帯域通過相関は、ベースバンド相関関数が示すよりもはるかに鋭いピークを常に示す。実際、比ｆC／Ｂの大きな値については、ピークの鋭さに対する信号帯域幅Ｂの寄与は無視することができ、ピークは主として搬送波周波数ｆCの値によって決まる。したがって、この潜在的に役立つ特性を何とかして活用することに関心が持たれる。しかしながら、帯域通過相関関数を時間遅延の測定のタスクに直接適用する試みは、以下の２つの主な理由から失敗するであろう。
【００１２】
−  相関包絡線内に搬送波の多くのサイクルがあるとき、最も大きな値を有するものを見つけることは、事実上不可能である。
−  搬送波サイクル数が比較的少なく、さらに個々のサイクルがはっきりと認識できるとき、搬送波の未知の位相が、包絡線の最大値（正確な遅延を示す）と搬送波の最大値との間に未知のオフセットをもたらす。この現象によって、時間遅延の測定は非常に不正確なものとなる。
【００１３】
もちろん、第１の場合でも搬送波の位相は未知であるが、搬送波のピークの「密度」が高いためにそれほど問題にはならない。
【発明の概要】
【発明が解決しようとする課題】
【００１４】
したがって、帯域通過相関関数の特徴的な構造を十分に活用することによって、帯域通過信号とその時間シフトされたレプリカとの間の時間シフトを求めるための方法及び装置を提供することが望ましい。
【００１５】
本発明の態様は、添付の特許請求の範囲に明確に記載されている。
【図面の簡単な説明】
【００１６】
ここで、添付図面を参照して本発明のさまざまな実施形態を述べることにする。
【００１７】
【図１】マイクロ波自動車レーダの簡略化したブロック図である。
【図２】ａ）は、帯域通過信号の電力スペクトル密度を示す図である。ｂ）は、ダウンコンバートの結果として帯域通過信号から得られたベースバンド信号の電力スペクトル密度を示す図である。ｃ）は、帯域通過信号の自己相関関数を示す図である。ｄ）は、ダウンコンバートの結果として帯域通過信号から得られたベースバンド信号の自己相関関数を示す図である。
【図３】搬送波周波数ｆCにおける帯域通過信号を使用して、ベースバンドバージョンと、中間周波数ｆkに位置するＫ個の周波数シフトされたレプリカとを生成する概念を模式的に示す図である。
【図４ａ】ベースバンド（ｋ＝０）及び９つの異なる中間周波数において得られる同相相互相関関数及び直交相互相関関数の例を示す図である。
【図４ｂ】ベースバンド（ｋ＝０）及び９つの異なる中間周波数において得られる同相相互相関関数及び直交相互相関関数の例を示す図の続きである。
【図５】すべての重みａkが等しいときに、９つの異なる中間周波数について得られる周期的櫛型関数ＤK（τ）の形状を示す図である。
【図６】重みａkが本発明に従って構成されたときに、９つの異なる中間周波数について得られる周期的櫛型関数ＤK（τ）の形状を示す図である。
【図７】位相角θの任意に選択された６つの値について、２つの成分ＲR（τ）及びＲI（τ）、並びに、圧縮相互相関関数ＲC（τ）の大きさ｜ＲC（τ）｜を示す図である。
【図８】相互相関包絡線ＲE（τ）と、等しい重みについて得られる圧縮相互相関関数ＲC（τ）の大きさ｜ＲC（τ）｜とを示す図である。
【図９】制御パラメータｐの３つの値０．０５、０．２、及び０．４について、Ｋ＝９のパスカル窓の形状を示す図である。
【図１０】種々のパスカル窓において得られる圧縮相互相関関数ＲC（τ）の大きさ｜ＲC（τ）｜を示す図である。
【図１１】相互相関包絡線ＲE（τ）と、制御パラメータｐ＝０．０５及びｐ＝０．４におけるパスカル窓において得られるその２つの圧縮バージョンとを示す図である。
【図１２】図１１のすべての相互相関関数を２乗したものを示す図である。
【図１３】本発明に従って構成された相関圧縮システムＣＣＲの機能ブロック図である。
【図１４】本発明に従って構成された相関圧縮システムＣＣＲを組み込むように修正されたマイクロ波自動車レーダのブロック図である。
【発明を実施するための形態】
【００１８】
搬送波周波数（すなわち、中心周波数）ｆCを有する帯域通過信号ｓ（ｔ）が、未知の量Δだけ遅延し、時間遅延したレプリカｓ（ｔ−Δ）が生成されるものと仮定する。信号ｓ（ｔ）は、対象となる或る障害物に向けて送信される無線周波数（ＲＦ）帯域通過信号又はマイクロ波帯域通過信号とすることができ、信号ｓ（ｔ−Δ）は、障害物によって反射されて戻ってくるものとみなすことができる。よく知られているように、障害物までの距離は、送信信号が受ける遅延Δから求めることができる。
【００１９】
例えば、送信信号ｓ（ｔ）は、変調波形によって適切に振幅変調及び／又は周波数変調された搬送波周波数ｆCの正弦搬送波の形態を有することができる。当業者にはよく知られているように、このような変調を搬送波に与えると、搬送波の電力が拡散される。このとき、当初は周波数ｆCに集中していた電力は、周波数ｆCを中心とした周波数帯域Ｂ0内に分散される。その結果、変調プロセスは、コヒーレント搬送波の１次余弦自己相関関数ｃｏｓ２πｆCτを、時間の経過により減衰する包絡線Ｒ（τ）を有する別の余弦関数Ｒ（τ）ｃｏｓ２πｆCτに変換する。包絡線Ｒ（τ）の形状は電力スペクトル分布の形状に依存し、包絡線の時間幅及び範囲は周波数帯域Ｂ0の増加と共に減少する。
【００２０】
元の送信信号ｓ（ｔ）とその遅延レプリカｓ（ｔ−Δ）の両方が、未知の遅延Δの値を求めるその後の処理に利用可能であると仮定する。本発明の別の態様によれば、遅延信号ｓ（ｔ−Δ）は、ベースバンド及び所定のＫ個の中間周波数ｆ1，ｆ2，…，ｆk，…，ｆKへダウンコンバートされる。その結果、信号ｓ（ｔ−Δ）についての以下に示す（Ｋ＋１）個の表現が得られる。
【００２１】
【数１】
0005383683.tif 000002
【００２２】
送信信号ｓ（ｔ）もベースバンド及び同じＫ個の中間周波数ｆ1，ｆ2，…，ｆk，…，ｆKへダウンコンバートされる。しかしながら、この場合には、信号ｓ（ｔ）についての２つの直交成分が生成される。したがって、送信信号をダウンコンバートした結果、以下に示す２（Ｋ＋１）個の表現が得られる。
【００２３】
【数２】
0005383683.tif 000003
【００２４】
ここで、ｘk（ｔ）及びｙk（ｔ）のそれぞれは、いわゆる中間周波数ｆkにおいて得られた同相成分（in-phase component）及び直交成分（quadrature component）である（明らかに、ｘ0（ｔ）及びｙ0（ｔ）はベースバンド成分である）。或いはまた、信号ｓ（ｔ）についてのこれら２（Ｋ＋１）個の表現は、（Ｋ＋１）個の複素表現とみなすこともできる。適切に選ばれた中間周波数において信号の直交成分を生成する技法としては、多くの種々の技法が知られている。
【００２５】
必須ではないが、或る適切に選ばれた基本周波数ｆ1の連続した高調波である中間周波数、すなわち、
【００２６】
【数３】
0005383683.tif 000004
【００２７】
を利用することが便利である。
【００２８】
本発明のさらなる態様によれば、ベースバンド及びＫ個の中間周波数において、２（Ｋ＋１）個の相関器を含むバンクを利用して、１次信号ｓ（ｔ）についての（Ｋ＋１）個の各複素表現と、遅延信号ｓ（ｔ−Δ）についての対応する各表現との間の相互相関関数が求められる。
【００２９】
図３は、開示される本概念を概略的に示す。搬送波周波数ｆCにおける帯域通過信号ｓ（ｔ）を使用して、そのベースバンドバージョンを生成すると共に、中間周波数ｆkに位置するＫ個の周波数シフトされたレプリカを生成する。このような操作は、搬送波周波数ｆCを中心として有限の周波数帯域を占有する帯域通過信号を、信号の帯域幅が使用される中間周波数のうちで最大のもの、すなわちｆKにまで（ほぼ）拡張された概念上のベースバンド「信号」への変換とみなすことができる。したがって、時間遅延を求める目的で帯域通過信号の帯域幅を拡張するには、この手法を使用することが有利である。１次送信信号ｓ（ｔ）は、それ自体は有用な情報を全く含まず、その特徴は効率的な電力送信を得ると共に、時間遅延の測定を容易にするように選択されることに留意されたい。
【００３０】
図３は、遅延信号ｓ（ｔ−Δ）が受ける同様の帯域幅拡張も示す。結果として生じる２つの信号ｓ（ｔ）及びｓ（ｔ−Δ）のスペクトルにおける特定の「タイル」構造は、対応する「タイル」（スペクトルレプリカ）によって表される信号成分を相互相関させるのに活用することができる。そして、得られた複数の相互相関関数を適切に結合して、明白な鋭いピークを有する「圧縮された（compacted）」関数を得ることができ、時間遅延の測定を容易にすることができる。
【００３１】
ダウンコンバート及び相互相関の操作を組み合わせて実行した結果として、相互相関関数の以下のセットが得られる。
−  ベースバンドにおいて、
【００３２】
【数４】
0005383683.tif 000005
【００３３】
−  Ｋ個の中間周波数ｆk、ｋ＝１，２，…，Ｋのそれぞれにおいて、
【００３４】
【数５】
0005383683.tif 000006
【００３５】
上記数式におけるθは、中心周波数ｆCにおける位相遅延に起因すると共に、対象となる物体による送信信号ｓ（ｔ）の反射の位相にも起因して生じる未知の位相角である。
【００３６】
関数ＲE（τ−Δ）は、上記すべての相互相関関数：Ｃk（τ）、Ｓk（τ）、ｋ＝０，１，…，Ｋの包絡線である。以下に示すように、この包絡線は２つのベースバンド相互相関関数から求めることができる。
【００３７】
【数６】
0005383683.tif 000007
【００３８】
ＲE（τ−Δ）の時間領域におけるピークの位置は、信号ｓ（ｔ）のいわゆる群遅延Δに対応し、従来から時間遅延の測定に使用されている。
【００３９】
以下では、余弦項を含む相互相関関数Ｃk（τ）、ｋ＝０，１，…，Ｋを同相相互相関関数（in phase cross-correlation function）と呼ぶと都合が良い。同様に、正弦項を含む相互相関関数Ｓk（τ）は、直交相互相関関数（quadrature cross-correlation function）と呼ばれる。
【００４０】
例示のために、かつ、本発明の理解を容易にするために、ベースバンド（ｋ＝０）及び９つの異なる中間周波数ｋｆ1、ｋ＝１，２，…，９において得られる同相相互相関関数及び直交相互相関関数の例を図４に示す。未知の位相角θの値はπ／３に設定されている。ここでは、遅延Δは０に設定されている。また、相対時間τは、包絡線ＲE（τ）が時間区間（−１，１）に及ぶように正規化され、基本周波数ｆ1の周期は２に等しい。図から分かるように、いずれの中間周波数ｆkにおいても、それぞれの相互相関関数Ｃk（τ）、Ｓk（τ）の最大値の位置と、その包絡線ＲE（τ）の最大値の位置とは一致しない。これは、観測不可能で且つ一般に未知の位相角θの存在によるものである。
【００４１】
本発明のさらなる態様によれば、すべての同相相互相関関数Ｃk（τ）は、適切な重みａkで平均化されて、圧縮相互相関関数（compacted cross-correlation function）ＲC（τ）の実数成分ＲR（τ）が生成される。
【００４２】
【数７】
0005383683.tif 000008
【００４３】
ここで、重みａkは、それらの重みの合計が１に等しくなるように正規化されている。したがって、
【００４４】
【数８】
0005383683.tif 000009
【００４５】
となる。重みａkの最適な選択は、この後で詳細に論じられる。
【００４６】
同様に、すべての直交相互相関関数Ｓk（τ）も、同じ所定の重みａkで平均化されて、圧縮相互相関関数ＲC（τ）の虚数成分ＲI（τ）が生成される。
【００４７】
【数９】
0005383683.tif 000010
【００４８】
本発明に従って構成された圧縮相互相関関数ＲC（τ）は、成分ＲR（τ）とＲI（τ）を含む複素関数であり、したがって、
【００４９】
【数１０】
0005383683.tif 000011
【００５０】
となる。
【００５１】
したがって、圧縮相互相関関数ＲC（τ）の大きさ｜ＲC（τ）｜は
【００５２】
【数１１】
0005383683.tif 000012
【００５３】
によって求めることができる。
【００５４】
大きさ｜ＲC（τ）｜が
【００５５】
【数１２】
0005383683.tif 000013
【００５６】
の形になることは容易に示すことができる。
【００５７】
上記数式は、次のように解釈することができる。圧縮相互相関関数ＲC（τ）の観測される大きさ｜ＲC（τ）｜は、相互相関包絡線ＲE（τ）による周期的櫛型関数ＤK（τ）
【００５８】
【数１３】
0005383683.tif 000014
【００５９】
の時間「ゲーティング（gating）」の結果である。
【００６０】
したがって、圧縮相互相関関数ＲC（τ）の大きさ｜ＲC（τ）｜のピークは、相互相関包絡線ＲE（τ）のピークと一致し、未知の遅延Δを求めるのに使用することができる。圧縮相互相関関数｜ＲC（τ）｜のピークは、相互相関包絡線ＲE（τ）のピークよりも大幅に狭くすることができるため、時間遅延測定の精度は、従来の時間遅延測定の技法よりも向上する。
【００６１】
さらに、圧縮相互相関関数ＲC（τ）の観測される大きさ｜ＲC（τ）｜は、未知の位相角θに依存しない。したがって、例えば送信パスの微細な変動、ドップラー効果、或いは周波数変動に起因する位相θの変化は、時間遅延の測定に悪影響を与えない。圧縮相互相関関数ＲC（τ）におけるこの特有の特性によって、未知の位相角θを取り扱う問題は解決される。
【００６２】
本発明に従って構成された圧縮相互相関関数ＲC（τ）の上記２つの特性（明白な鋭いピーク及び位相からの独立性）によって、開示される本技法は、従来技術の技法よりも優れたものとなる。
【００６３】
すべての重みａk、ｋ＝０，１，…，Ｋが同じ値であり、１／（Ｋ＋１）に等しい場合には、周期的櫛型関数ＤK（τ）は、以下のような閉じた式で表すことができる。
【００６４】
【数１４】
0005383683.tif 000015
【００６５】
上記数式は、相互相関包絡線ＲE（τ）において指定された圧縮を達成するのに必要な中間周波数の個数Ｋを求めるのに使用することができる。
【００６６】
図５は、基本周波数ｆ1の周期が２に設定され、すべての重みａkが等しいときに、９つの異なる中間周波数ｋｆ1、ｋ＝１，２，…，９について得られる周期的櫛型関数ＤK（τ）の形状を示す。この形状は、複数の特徴的なサイドローブを呈する。櫛型関数ＤK（τ）の形状は、未知の位相角θの値に関係なく不変であることを指摘しておく。
【００６７】
比較のために、以降でより詳細に説明する本発明のさらに別の態様に従って重みａkが構成されたときの周期的櫛型関数ＤK（τ）の形状を図６に示す。このとき、図から分かるように、櫛型ＤK（τ）は、目に見えるサイドローブをもたない単峰性の「歯」（unimodal ‘teeth’）を備える。
【００６８】
本発明に従って構成された圧縮相互相関関数ＲC（τ）の大きさが未知の位相角θに依存しないことを示すために、位相角θの任意に選択された６つの値について、２つの成分ＲR（τ）とＲI（τ）、及び、関数ＲC（τ）の大きさ｜ＲC（τ）｜を図７に示す。図から分かるように、成分ＲR（τ）とＲI（τ）の形状は大きく変化するが、大きさの形状は位相角θによる影響を受けていない。図７に示す圧縮相互相関関数は、図４に示した形状を有する１０対の相互相関関数を処理することによって得られたものである。平均化に使用された重みａkは、図６に示した櫛型ＤK（τ）を生成するのに使用されたものと同じである。
【００６９】
図８は、等しい重みについて得られた圧縮相互相関関数ＲC（τ）の大きさ｜ＲC（τ）｜を示す。比較のために、相互相関包絡線ＲE（τ）も示されている。図から分かるように、圧縮相互相関関数の大きさ｜ＲC（τ）｜は、図５に示した櫛型関数ＤK（τ）のサイドローブに由来する特徴的なサイドローブを呈する。
【００７０】
いくつかの実際の用途では、圧縮相互相関関数の大きさ｜ＲC（τ）｜に現れるサイドローブは、相互相関ピークを位置特定する際に不確定性をもたらさない限り、許容することができる。しかしながら、対象となる遅延信号ｓ（ｔ−Δ）が雑音や他の干渉によって破損されるときには、相互相関ピークを広げるという不可避な犠牲を払ってでも、サイドローブを抑制することが好ましい。
【００７１】
したがって、単一の「制御」パラメータを変化させることによって、重みの配分を一様なもの（狭い応答を生成するすべて等しい重み）から一様でないものに次第に変換し、その結果、応答がサイドローブを有しない単峰性となるように重みａkを構成するための方法を開発することが役に立つ。
【００７２】
本発明の別の態様によれば、重みａkは、以下の方法で２項係数（すなわち、パスカルの三角形の行の要素）から求められる。初めに、重みａkの個数（Ｋ＋１）が指定されると、１次係数Ａkを以下のように求める。
【００７３】
【数１５】
0005383683.tif 000016
【００７４】
ここで、ｐは区間（０，１）から選択される値を取る制御パラメータである。ｐ＝０のとき、すべての重み（Ａ0を除く）は等しい。（正確な演算のために、係数Ａ0は常に、それぞれの２項係数がｐ乗されたものの２分の１に等しい）。次に、以下のように、１次係数Ａkを正規化することによって、重みａkが得られる。
【００７５】
【数１６】
0005383683.tif 000017
【００７６】
信号処理において、解析される波形を修正するのに使用される係数のセットは、窓（window）と呼ばれる。したがって、本発明に従って構成された重みのセット｛ａk｝は、重みがパスカルの三角形の適切な要素から得られたものであるため、（１パラメータ）パスカル窓と呼ばれる。
【００７７】
図９は、制御パラメータｐの３つの異なる値０．０５、０．２、及び０．４について、ｋ＝０，１，…，９におけるパスカル窓の形状を示したものである。図１０は、それぞれのパスカル窓に対応する圧縮相互相関関数ＲC（τ）の大きさ｜ＲC（τ）｜を示したものである。図から分かるように、ｐ＝０．４のとき、関数｜ＲC（τ）｜は、明らかに滑らかな傾斜を有する単峰性である。ｐ＝０．２のとき、この関数は依然として単峰性であるが、傾斜はいくつかの小さな起伏効果を呈する。最後に、ｐ＝０．０５のとき、小さなサイドローブが現れているため、この関数はもはや単峰性ではない。しかしながら、予想されたように、大きさ｜ＲC（τ）｜の幅は制御パラメータｐの値が大きいほど広くなる。したがって、パラメータｐの値は、干渉の予想レベル及び意図した用途の他の態様を考慮に入れて、賢明な方法で選択されるべきである。
【００７８】
当業者に知られている他の窓も利用できることを指摘しておく。さらに、数値最適化を行うことによって、一組の特定の基準を満たす窓を設計することもできることを指摘しておく。しかしながら、本発明に従って構成されたパスカル窓は、設計が容易であり、開示された方法を使用する時間遅延の測定によく適しているように思われる。
【００７９】
周期的櫛型関数ＤK（τ）は、その櫛型関数ＤK（τ）の唯一の主ピークが「捕捉」されるように、相互相関包絡線ＲE（τ）によって「ゲーティング」されるべきである。したがって、基本周波数ｆ1等の使用される最小周波数の周期は、好ましくは、相互相関包絡線ＲE（τ）の時間範囲、すなわち台（support）に等しくあるべきである。一般に、圧縮相互相関関数｜ＲC（τ）｜の幅は、使用されるすべての中間周波数における最大のもの、すなわちｆKに反比例する。
【００８０】
時間遅延の測定の改善を示すために、相互相関包絡線ＲE（τ）と、その２つの圧縮バージョンとを図１１に示す。ここで、２つの圧縮バージョンは、ベースバンドに９つの中間周波数を加えたものを利用して、制御パラメータの２つの異なる値、ｐ＝０．０５及びｐ＝０．４を有するパスカル窓を適用することによって得られたものである。基礎となる関数ＲE（τ）の幅は、（２分の１の高さの）４．５倍（ｐ＝０．０５）及び３倍（ｐ＝０．４）に削減されている。明らかに、このような削減によってピーク位置の局所性は大幅に改善され、したがって、時間遅延の測定はより正確になる。
【００８１】
圧縮相互相関関数ＲC（τ）の大きさ｜ＲC（τ）｜の代わりに、大きさの２乗である
【００８２】
【数１７】
0005383683.tif 000018
【００８３】
を使用することによって、演算が簡単になると共に、関数のピークがより明白になる。したがって、時間遅延の測定がさらに改善される。比較のために、図１１のすべての相互相関関数を２乗した結果を図１２に示す。
【００８４】
開示される技法は、送信信号の帯域幅を増加させることなく、改善された時間遅延計測を提供する。逆に言えば、より狭帯域の信号を送信して、必要とされる精度を得ることができる。この特性は、マルチユーザ環境或いは高密度信号環境を対象としたシステムを構成するときに特に有利である。
【００８５】
図１３は、本発明に従って構成された相関圧縮システムＣＣＲの機能ブロック図である。このシステムは、
−  多重周波数ジェネレータＢＬＦ、
−  （Ｋ＋１）個の直交ダウンコンバータＤＴ０、ＤＴ１、…、ＤＴＫのバンク、
−  （Ｋ＋１）個のダウンコンバータＤＲ０、ＤＲ１、…、ＤＲＫのバンク、
−  相関器の２つのバンク、第１のバンクは（Ｋ＋１）個の相関器ＸＣ０、ＸＣ１、…、ＸＣＫを備え、第２のバンクは相関器ＸＳ０、ＸＳ１、…、ＸＳＫを備える、
−  ２つの平均化回路ＡＶＣ及びＡＶＳ、
−  （ピタゴラス）関数結合器ＰＴＨ、
−  制御ブロックＣＴＢ、
を備える。
【００８６】
システムの主要部は、（Ｋ＋１）個の同一のチャネルを含むものとみなすことができる。各チャネルは、異なる中間周波数ｆkにおいてではあるが、同じ機能及び動作を実行する。したがって、一例として、チャネル１の動作をより詳細に説明することにする。このチャネルは、回路ＤＴ１、ＸＣ１、ＸＳ１、及びＤＲ１から成る。
【００８７】
送信信号ｓ（ｔ）等の基準信号ＳＴは、直交ダウンコンバータＤＴ１の一方の入力に印加される。直交ダウンコンバータＤＴ１の他方の入力は、多重周波数ジェネレータＢＬＦによって供給される正弦波信号Ｆ１によって駆動される。正弦波信号Ｆ１の周波数は、必要とされる所定の中間周波数ｆ1を得るように選択される。ダウンコンバータＤＴ１は、入力信号ＳＴの周波数ｆ1における２つの直交バージョンＩ１及びＱ１を生成する。
【００８８】
受信信号ｓ（ｔ−Δ）等の遅延信号ＳＲは、ダウンコンバータＤＲ１の一方の入力に印加される。ダウンコンバータＤＲ１の他方の入力は、多重周波数ジェネレータＢＬＦによって供給される正弦波信号Ｆ１によって駆動される。正弦波信号Ｆ１の周波数は、必要とされる所定の中間周波数ｆ1を得るように選択される。ダウンコンバータＤＲ１は、入力信号ＳＲの周波数ｆ1におけるバージョンＲ１を生成する。
【００８９】
直交ダウンコンバータＤＴ１によって生成された信号Ｉ１及びＱ１と、ダウンコンバータＤＲ１によって生成された信号Ｒ１は、すべて同じ中心周波数ｆ1を有する。信号Ｉ１及びＲ１が相関器ＸＣ１で処理されて、相互相関関数Ｃ１が得られる。同様に、信号Ｑ１及びＲ１が相関器ＸＳ１で処理されて、相互相関関数Ｓ１が得られる。
【００９０】
（Ｋ＋１）個のすべてのチャネルは並列に動作しているため、これらのチャネルは、（Ｋ＋１）個の同相相互相関関数Ｃ０、Ｃ１、…、ＣＫ、及び（Ｋ＋１）個の直交相互相関関数Ｓ０、Ｓ１、…、ＳＫを併せて生成する。得られたすべての同相相互相関関数は、平均化回路ＡＶＣにおいて適切な重みで平均化される。同様に、得られたすべての直交相互相関関数も、平均化回路ＡＶＳにおいて適切な重みで平均化される。
【００９１】
関数結合器ＰＴＨは、平均化回路ＡＶＣ及びＡＶＳによってそれぞれ供給された平均ＲＲ及びＲＩを利用して、圧縮相互相関関数を表す出力信号ＲＣを生成する。圧縮相互相関関数のピークの位置は、２つの入力信号ＳＴとＳＲとの間の未知の遅延Δを示す。そして、適切な外部プロセッサがこの遅延の値を求めるために、出力ＲＣを使用することができる。
【００９２】
また、相関圧縮システムＣＣＲは、システムの正確な動作に必要な適切なタイミング及び同期信号を供給する制御ブロックＣＴＢも備える。制御信号ＴＣは、（Ｋ＋１）個の相関器ＸＣ０、ＸＣ１、…、ＸＣＫのバンクに印加される。それに対して、制御信号ＴＳは、（Ｋ＋１）個の相関器ＸＳ０、ＸＳ１、…、ＸＳＫのバンクに印加される。また、２つの平均化回路ＡＶＣ及びＡＶＳは、ブロックＣＴＢから制御信号ＴＡを受け取る。制御ブロックＣＴＢは、自律的に動作してもよいし、或いは入力ＣＣに現れる外部信号によって制御されてもよい。
【００９３】
当業者には既知であるように、上記の機能及び動作のうちのいくつか、さらにはすべてを、汎用プロセッサ又は特注プロセッサによってディジタルに実装することができる。しかしながら、ギガヘルツ領域では、依然として、必要な動作をアナログハードウェアで少なくとも部分的に実装するのがより便利である。
【００９４】
当業者には既知であるように、相互相関処理は適切に構成されたマッチドフィルタによって実行することもできる。しかしながら、一般に、マッチドフィルタリングは、継続時間が有限であり形状が既知の信号に使用される。したがって、マッチドフィルタリングは連続信号には適しておらず、特に、非決定論的信号又は非決定論的変調がなされた信号には適していない。
【００９５】
図１４は、図１に示した従来技術のシステムを修正した機能ブロック図である。図から分かるように、ダウンコンバータＤＣＲと相関器ＣＯＲは、本発明に従って構成された相関圧縮システムＣＣＲに置き換えられている。従来技術のシステムと比較すると、修正されたシステムは、求める時間遅延の改善された推定値を提供する。
【００９６】
本発明の好ましい実施形態の上記説明は、例示及び説明の目的で提示されたものである。この説明は、網羅的であることも、開示される正確な形態に本発明を限定することも意図するものではない。上記説明に鑑みると、予想される特定の使用に適したさまざまな実施形態で本発明を利用するために、多くの変更、修正、及び変形が当業者に可能であることは明らかである。
【図１】
0005383683.tif 000019
【図２】
0005383683.tif 000020
【図３】
0005383683.tif 000021
【図４ａ】
0005383683.tif 000022
【図４ｂ】
0005383683.tif 000023
【図５】
0005383683.tif 000024
【図６】
0005383683.tif 000025
【図７】
0005383683.tif 000026
【図８】
0005383683.tif 000027
【図９】
0005383683.tif 000028
【図１０】
0005383683.tif 000029
【図１１】
0005383683.tif 000030
【図１２】
0005383683.tif 000031
【図１３】
0005383683.tif 000032
【図１４】
0005383683.tif 000033