JP3973902B2 - センサ配置決定方法及びその装置並びに受信装置 - Google Patents
センサ配置決定方法及びその装置並びに受信装置 Download PDFInfo
- Publication number
- JP3973902B2 JP3973902B2 JP2001400964A JP2001400964A JP3973902B2 JP 3973902 B2 JP3973902 B2 JP 3973902B2 JP 2001400964 A JP2001400964 A JP 2001400964A JP 2001400964 A JP2001400964 A JP 2001400964A JP 3973902 B2 JP3973902 B2 JP 3973902B2
- Authority
- JP
- Japan
- Prior art keywords
- steering vector
- correlation amount
- incident
- sensor
- vector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Description
【発明の属する技術分野】
本発明は、複数のセンサの各々に入射される複数の入射信号の到来方位を推定する際に発生する偽像を最小限に押えるように複数のセンサの配置を決定するセンサ配置決定方法及びその装置並びに受信装置に関する。
【0002】
【従来の技術】
従来、スーパーレゾリューション(超分解能)による入射信号の入射方位の測定や、干渉波抑圧により複数の入射信号から1つの入射信号を分離するヌルステアリングという技術が知られている。特開2000−216620号公報は、このような技術を採用した「受信装置」を開示している。
【0003】
この従来の受信装置は、スーパーレゾリューションの1つであるMUSIC(Multiple Signal Classification)法と、ヌルステアリングの1つであるアダプティブアレイ技術といった各々独立した技術を組合せることにより、複数のアンテナで受信された信号を用いて、同一周波数の複数の到来波の方位を測定し、且つ複数の到来波の中から不要波を抑圧して希望波のみを分離するといった2つの機能を同時に実現している。
【0004】
【発明が解決しようとする課題】
また、複数のセンサを所定間隔で配置したセンサアレイシステムにより、広帯域の周波数範囲の入射信号を受信しようとした場合には、センサ間隔を受信信号の半波長間隔以内に配置することが不可能となる。この場合、入射信号の到来方位を求めるアルゴリズムを用いたMUSIC法により、到来角に対するスペクトラムを求め、該スペクトラムに基づいて複数の入射信号の到来方位を推定すると、スペクトラム上には本来の到来方位に相当する部分にピークを持つ真像が発生するとともに、本来の到来方位とは異なる方位に相当する部分にピークを持つ偽像が発生する。このため、この偽像を真像と誤って認識してしまい、信号の到来していない方位を推定する場合がある。
【0005】
本発明は、入射信号の到来方位を推定する際に発生する偽像を最小限に押えることができ、これによって到来方位を正確に推定できるセンサ配置決定方法及びその装置並びに受信装置を提供することにある。
【0006】
【課題を解決するための手段】
本発明の第1の態様に係るセンサ配置決定方法は、上記課題を達成するために、同一帯域内に複数の入射信号が混信している場合に、同一帯域内の複数の入射信号が同時に入射される複数のセンサを配置し、前記各センサにおける前記複数の入射信号の位相応答を行列で表したステアリングベクトルを生成し、生成されたステアリングベクトルに基づいて前記ステアリングベクトルの自己相関量と相互相関量とからなる自己内相関量を算出し、算出された前記相互相関量が所定値か否かを判定し、前記相互相関量が所定値以下になっていない場合には前記複数のセンサの配置を変更することを特徴とする。
【0007】
この本発明の第1の態様に係るセンサ配置決定方法によれば、ステアリングベクトルに基づいて算出された自己内相関量の内の相互相関量が所定値以下か否かを判定し、相互相関量が所定値以下になっていない場合には複数のセンサの配置を変更するので、入射信号の到来方位を推定する際に発生する偽像を最小限に押えることができ、これによって到来方位を正確に推定できる。
【0008】
また、第1の態様に係るセンサ配置決定方法においては、前記ステアリングベクトルは、仰角を固定し且つ試験電波を各到来方向から前記各センサに入射したときに前記各センサにおいて所定の角度毎に測定したデータから構成できる。また、生成されたステアリングベクトルに基づいて前記ステアリングベクトルの複素共役転置行列を求め、前記複素共役転置行列と前記ステアリングベクトルとを演算することにより前記ステアリングベクトルの前記自己内相関量を求めるように構成できる。
【0009】
本発明の第2の態様に係るセンサ配置決定装置は、上記課題を達成するために、同一帯域内に複数の入射信号が混信している場合に、同一帯域内の複数の入射信号が同時に入射される複数のセンサと、前記各センサにおける前記複数の入射信号の位相応答を行列で表したステアリングベクトルを生成するベクトル生成部と、前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの自己相関量と相互相関量とからなる自己内相関量を算出する相関量算出部と、前記相関量算出部で算出された前記相互相関量が所定値か否かを判定し、前記相互相関量が所定値以下になっていない場合には前記複数のセンサの配置を変更するセンサ配置決定部とを備えたことを特徴とする。
【0010】
また、第2の態様に係るセンサ配置決定装置においては、前記ベクトル生成部は、仰角を固定し且つ試験電波を各到来方向から前記各センサに入射したときに前記各センサにおいて所定の角度毎に測定したデータを前記ステアリングベクトルとするように構成できる。また、前記相関量算出部は、前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの複素共役転置行列を求め、前記複素共役転置行列と前記ステアリングベクトルとを演算することにより前記ステアリングベクトルの前記自己内相関量を求めるように構成できる。
【0011】
本発明の第3の態様に係る受信装置は、同一帯域内に複数の入射信号が混信している場合に、同一帯域内の複数の入射信号が同時に入射される複数のセンサと、前記各センサにおける前記複数の入射信号の位相応答を行列で表したステアリングベクトルを生成するベクトル生成部と、前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの自己相関量と相互相関量とからなる自己内相関量を算出する相関量算出部と、前記相関量算出部で算出された前記相互相関量が所定値以下か否かを判定し、前記相互相関量が所定値以下になっていない場合には前記複数のセンサの配置を変更するセンサ配置決定部と、前記センサ配置決定部で前記複数のセンサの配置が決定されたときにおける、前記相関量算出部で算出された前記自己内相関量に基づいて到来角に対するスペクトラムを求め、該スペクトラムに基づいて前記複数の入射信号の到来方位を求める到来方位演算部とを備えたことを特徴とする。
【0012】
この構成によれば、到来方位演算部は、センサ配置決定部で複数のセンサの配置が決定されたときにおける、相関量算出部で算出された自己内相関量に基づいて到来角に対するスペクトラムを求め、スペクトラムに基づいて複数の入射信号の到来方位を求めるので、スペクトラム上における偽像は最小限に押えられるから、複数の入射信号の到来方位を正確に推定できる。
【0013】
また、本発明の第3の態様に係る受信装置は、同一帯域内に複数の入射信号が混信している場合に、同一帯域内の複数の入射信号が同時に入射される複数のセンサと、前記各センサにおける前記複数の入射信号の位相応答を行列で表したステアリングベクトルを生成するベクトル生成部と、前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの自己相関量と相互相関量とからなる自己内相関量を算出する相関量算出部と、前記相関量算出部で算出された前記相互相関量が所定値以下か否かを判定し、前記相互相関量が所定値以下になっていない場合には前記複数のセンサの配置を変更するセンサ配置決定部と、前記複数のセンサからの前記複数の入射信号が混合された複数の混合信号を所定の狭帯域にそれぞれ制限して出力する複数の帯域制限ろ波器と、前記複数の帯域制限ろ波器からの狭帯域に制限された複数の信号に基づいて、独立成分分析の手法により前記複数の入射信号を分離・抽出する独立成分分析部とを備えたことを特徴とする。
【0014】
この構成によれば、第1の態様に係るセンサ配置決定方法の作用及び効果と同様な作用及び効果が得られるとともに、受信した各センサからの混合信号を狭帯域に帯域制限し、その狭帯域に帯域制限された各アンテナの信号を用いて独立成分分析を行うので、センサ毎の入射信号の変調成分で見たときに、事実上時間差が無視できるようになり、センサ毎の信号到達時間の差が生じる条件下でも一般的なブラインド信号分離の手法によりアンテナからの観測信号である混合信号中に含まれる原信号を分離・抽出できる。
【0015】
また、本発明の第3の態様に係る受信装置においては、前記ベクトル生成部は、仰角を固定し且つ試験電波を各到来方向から前記各センサに入射したときに前記各センサにおいて所定の角度毎に測定したデータを前記ステアリングベクトルとするように構成できる。また、前記相関量算出部は、前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの複素共役転置行列を求め、前記複素共役転置行列と前記ステアリングベクトルとを演算することにより前記ステアリングベクトルの前記自己内相関量を求めるように構成できる。
【0016】
【発明の実施の形態】
以下、本発明の実施の形態を図面を参照しながら詳細に説明する。
【0017】
(第1の実施の形態)
本発明の第1の実施の形態に係るセンサ配置決定装置は、センサとしてアンテナを使用し、複数のアンテナの各々に複数の入射信号を入射し、広帯域の周波数範囲の信号を受信する上で、複数のアンテナの配置と入射信号の条件から求められるステアリングベクトルの自己内相関量の内の相互相関量に基づいて、偽像を最小とするように複数のアンテナの配置を決定する。
【0018】
なお、以下では、説明を簡単にするために、アンテナの数を「4」とし、入射信号の数を「n=4」として説明するが、アンテナの数及び入射信号の数はこれらに限定されず任意である。
【0019】
図1は、本発明の第1の実施の形態に係るセンサ配置決定装置の構成を示すブロック図である。このセンサ配置決定装置は、第1〜第4アンテナ11〜14、第1〜第4帯域制限ろ波器21〜24、ベクトル生成部27、相関量算出部28、センサ配置決定部29から構成されている。
【0020】
第1〜第4アンテナ11〜14としては、バーチカルアンテナ、ダイポールアンテナといった無指向性のアンテナ、及び任意の指向性を持ったアンテナ等が用いられ、種々の方位からの電波を受信する。これら第1〜第4アンテナ11〜14を設置する間隔や高さは任意である。第1〜第4アンテナ11〜14は、空中からの複数の入射信号(電波)S1〜S4を受信し、これらが混合された混合信号を第1〜第4帯域制限ろ波器21〜24に送る。
【0021】
第1〜第4帯域制限ろ波器21〜24は、第1〜第4アンテナ11〜14からの混合信号に含まれる所定帯域の周波数成分のみを通過させてベクトル生成部27に送る。なお、各第1〜第4帯域制限ろ波器21〜24が通過させる周波数帯域は同じである。
【0022】
帯域制限制御器25は、通過させる周波数帯域、つまり通過させる周波数成分の範囲を指定するための制御信号を生成する。この帯域制限制御器25で生成された制御信号は、第1〜第4帯域制限ろ波器21〜24に送られる。第1〜第4帯域制限ろ波器21〜24は、帯域制限制御器25からの制御信号に従って、入力された混合信号に含まれる所定帯域の周波数成分のみを通過させる。
【0023】
この帯域制限制御器25は、任意の帯域の周波数成分を通過させるような制御信号を生成できるように構成されている。したがって、帯域制限制御器25からの制御信号を適宜変更することにより、第1〜第4帯域制限ろ波器21〜24を通過する周波数帯域を任意に変化させることができる。
【0024】
ベクトル生成部27は、第1〜第4帯域制限ろ波器21〜24からの信号に基づいて、各アンテナ11〜14における複数の入射信号の位相応答を行列で表したステアリングベクトルを生成する。相関量算出部28は、ベクトル生成部27で生成されたステアリングベクトルに基づいてステアリングベクトルの複素共役転置行列を求め、複素共役転置行列とステアリングベクトルとを演算することによりステアリングベクトルの自己内相関量を求める。アンテナ配置決定部29は、相関量算出部28で算出された相互相関量が所定値以下になるように第1〜第4アンテナ11〜14の配置を決定する。
【0025】
次に、図3のフローチャートを参照しながらセンサ配置決定装置により実現されるセンサ配置決定方法を説明する。
【0026】
今、n個(nは2以上の整数)のアンテナのアレイに同時にm個(mは2以上の整数)の信号が入射する場合を仮定する。この状態は、下記式(1)で表すことができる。
【0027】
【数1】
ここで、X1〜nは各アンテナで観測される時系列データ、A=[a(θ1),…a(θm)]はアンテナの配置と特性とで決まるn行×m列の信号混合の行列、S1〜m(t)はm個の入射信号、N1〜n(t)は各アンテナにおける雑音である。
【0028】
一般に、スーパーレゾリューションによる方位測定やヌルステアリングにおいては、上記式(1)における「A」に相当する情報が直接的又は間接的に推定される。
【0029】
ここで、[A]をステアリングベクトルaとする。このステアリングベクトルaは、ステアリングベクトル生成部27で生成される。図2に示すように、仰角を固定し、試験アンテナ10を回転させながら試験アンテナ10から試験電波を第1〜第4アンテナ11〜14に向けて放射し、第1〜第4アンテナ11〜14が試験アンテナ10から試験電波を各到来方向から入射したときに所定の角度毎(例えば1degステップ)にデータを測定し、測定されたたデータをステアリングベクトルaとする(ステップS11)。このステアリングベクトルaは、各アンテナ11〜14における複数の入射信号の位相応答を行列で表したものであり、[4×360]の行列、即ち、式(2)で表される。
【0030】
【数2】
なお、このステアリングベクトルaは、例えばシュミレーションにより測定しても良い。このステアリングベクトルaは、アンテナのメインローブの方向を決めるものである。
【0031】
次に、相関量算出部28は、ベクトル生成部27で生成されたステアリングベクトルに基づいてステアリングベクトルの複素共役転置行列を求め、複素共役転置行列とステアリングベクトルとを演算することによりステアリングベクトルの自己内相関量Pを求める(ステップS13)。
【0032】
この自己内相関量Pは、式(3)で表される。
【0033】
【数3】
ここで、上添字Hは複素共役転置を表す。そして、自己内相関量Pは、式(4)で示すように[360×360]の行列となる。
【0034】
【数4】
この行列の対角成分が自己相関部分(本発明の自己相関量に対応)であり、その他の成分が相互相関部分(本発明の相互相関量に対応)となる。そして、自己内相関量から相互相関量が抽出される(ステップS15)。この相互相関量からステアリングベクトルaの相関関係が求まる。例えば、相互相関量p40・2は、ステアリングベクトルaの40degと2degとの成分の相関値を表す。相関関係が強い場合には、相互相関量p40・2の値は大きくなり、相関関係が弱い場合には、相互相関量p40・2の値は小さくなる。
【0035】
相互相関量p40・2が大きい値を示した場合には、到来波が40degから到来した場合、40deg以外の2degからも到来したという、本来の到来方位以外にも偽像として方位を示すアンビギュティ(曖昧性)が発生する。また、相互相関量が大きい場合には、第2の実施の形態において後述するように、MUSICスペクトラム上の偽像の値も大きくなる。逆に、相互相関量が小さいほど、MUSICスペクトラム上の偽像が小さい。
【0036】
このため、アンテナ配置決定部29は、相関量算出部28で算出された相互相関量が所定値以下か否かを判定し(ステップS17)、相互相関量が所定値以下になっていない場合には、第1〜第4アンテナ11〜14の配置を変更し(ステップS19)、ステップS11に戻って、再度、ステップS11〜ステップS17までの処理を繰り返し行う。そして、相互相関量が所定値以下になった場合には、そのときの第1〜第4アンテナ11〜14の配置を最適な配置として決定する(ステップS21)。
【0037】
このように第1の実施の形態に係るセンサ配置決定方法及びその装置によれば、生成されたステアリングベクトルに基づいてステアリングベクトルの自己相関量と相互相関量とからなる自己内相関量を算出し、算出された相互相関量が所定値以下になるように複数のアンテナの配置を決定するので、所望の周波数範囲において、入射信号の到来方位を推定する際に発生する偽像を最小限に押えることができる。また、複数のアンテナの最適な配置としては、一直線上に所定間隔で配置される場合、円周上に所定間隔で配置される場合、あるいは、円弧上に所定間隔で配置される場合等が挙げられる。
【0038】
なお、一例として、MUSIC法の場合を説明したが、MUSIC法以外の他のアルゴリズムでも同様な関係を求めることができる。また、アンテナ配置決定部29は、例えば、相互相関量が自己相関量の1/n(n>1で例えばn=7)以下になるように第1〜第4アンテナ11〜14の配置を決定しても良く、同様な効果が得られる。
【0039】
(第2の実施の形態)
図4に示す第2の実施の形態に係る受信装置は、図1に示す第1の実施の形態に係るセンサ配置決定装置を備えた受信装置であり、第1の実施の形態のセンサ配置決定装置に対して、さらにサンプリング部30及び到来方位演算部47を追加した点が異なる。なお、図4に示す部分において、図1に示す部分と同一部分には同一符号を付し、その説明は省略する。ここでは、異なる部分についてのみ説明する。
【0040】
サンプリング部30は、第1〜第4中間周波数変換器31〜34、局部発振器35、第1〜第4A/D変換器41〜44及び発振器45から構成されている。
【0041】
局部発振器35は、受信電波を中間周波数に変換するために必要とする発信周波数を有する信号を生成する。この局部発振器35で生成された信号は、第1〜第4中間周波数変換器31〜34に送られる。
【0042】
第1〜第4中間周波数変換器31〜34の各々は、何れも図示を省略するが、高周波増幅器、周波数混合器及び中間周波数増幅器から構成されている。
【0043】
高周波増幅器は、受信周波数帯の高周波を、次段の周波数混合器の入力電圧として適当な大きさになるように増幅する。周波数混合器は、高周波増幅器で増幅された信号と局部発振器35の出力信号とを混合し、それらの和又は差の周波数を作ることにより中間周波数の信号に変換する。中間周波数増幅器は、受信電波の周波数を、より低い中間周波数に変換して増幅する。これにより、安定で高利得の増幅を行うことができ、感度を向上させることができる。
【0044】
上記のように構成される第1〜第4中間周波数変換器31〜34から出力される信号は、第1〜第4A/D変換器41〜44にそれぞれ送られる。
【0045】
発振器45は、第1〜第4中間周波数変換器31〜34からの信号をサンプリングするためのサンプリングクロックを生成する。この発振器45で生成されたサンプリングクロックは第1〜第4A/D変換器41〜44に送られる。
【0046】
第1〜第4A/D変換器41〜44は、発振器45からの信号をサンプリングクロックとして、第1〜第4中間周波数変換器31〜34からのアナログ信号をサンプリングすることにより、デジタル信号にそれぞれ変換する。第1〜第4A/D変換器41〜44の各々から出力されるデジタル信号は、到来方位演算部47及びベクトル生成部27に供給される。
【0047】
ベクトル生成部27は、第1〜第4A/D変換器41〜44からの信号に基づいて、各アンテナ11〜14における複数の入射信号の位相応答を行列で表らしたステアリングベクトルを生成する。
【0048】
このような構成の受信装置によれば、ベクトル生成部27、相関量算出部28及びアンテナ配置決定部29を設けたので、第1の実施の形態に係るセンサ配置決定装置の効果と同様な効果が得られる。即ち、入射信号の到来方位を推定する際に発生する偽像を最小限に押えることができる。
【0049】
また、到来方位演算部47は、第1〜第4A/D変換器41〜44からの信号と相関量算出部28で算出された自己内相関量とに基づいて、MUSIC等のアルゴリズムを用いて入射信号の到来方位を測定する。この到来方位演算部47の処理を図5を参照して詳細に説明する。
【0050】
MUSIC法は相関行列の固有値と固有ベクトルとを用いた推定法である。図5のようにアンテナ間隔dのM素子等間隔リニアアレーに平面波がK波到来していて、各到来波の信号波形と到来角がFk(t),θk(k=1,2…K)と表されるとき、各アンテナにおける各到来波の位相応答を表す方向ベクトルa(θk)は、式(5)で与えられる。
【0051】
【数5】
ここで、上添字Tは転置を表す。よって、入力ベクトルは式(6)〜式(10)で表される。
【0052】
【数6】
上式においてN(t)は熱雑音ベクトルであり、その成分は平均が0で分散(電力)がσ2の独立な複素ガウス過程である。このとき、アンテナ間の相関特性を表す相関行列は式(11)〜式(12)で与えられる。
【0053】
【数7】
ここで、上添字Hは複素共役転置を表す。到来波が互いに無関係であれば信号相関行列SのランクはKとなる。また、方向行列AもランクはKである。従って、この場合の相関行列RxxはランクKの非負定値エルミート行列となる。この行列の固有値λi(i=1,2…,M)は実数となり、下記の式(13)の関係を有する。
【0054】
【数8】
従って、相関行列の固有値を求め、熱雑音電力σ2より大きい固有値の数から到来波数Kを推定することができる。また、固有値λi(i=1,2…,M)に対応する固有ベクトルをei(i=1,2…,M)とすると、M次元のエルミート空間の正規直交基底ベクトルとして扱われる。この空間は信号空間span{e1,…eK}と雑音空間span{eK+ 1,…eM}との二つの部分空間にわけることができ、信号空間と雑音空間とは互いに直交補空間の関係にある。
【0055】
span{e1,…eK}はベクトルei(i=1,2…,M)で張られる空間とする。また、信号空間は方向ベクトルを用いて、span{a(θ1),…,a(θk)}と表すことができる。従って、熱雑音電力に等しい固有値に対応する固有ベクトルは全て到来波の方向ベクトルと直交することになる。そこで、式(14)のような評価関数を定義する。
【0056】
【数9】
これはMUSICスペクトラムと呼ばれ、到来角θに対するスペクトラムのK個のピークが到来方位θk(k=1,2…K)となる。なお、式(13)からもわかるように、熱雑音電力に等しい最小固有値が少なくとも一つ必要なので、アレーのセンサ数はM≧K+1が必要条件となる。
【0057】
ここで、式(14)の分子の部分が、式(3)の自己内相関量Pとなる。このため、到来方位演算部47は、アンテナ配置決定部27で第1〜第4アンテナ11〜14の配置が決定されたときにおける、相関量算出部28で算出された自己内相関量Pを用いて式(14)に示すMUSICスペクトラムを求め、該MUSICスペクトラムに基づいて複数の入射信号の到来方位を求める。
【0058】
即ち、第1〜第4アンテナ11〜14の配置が決定されたときにおける自己内相関量Pの内の相互相関量は、所定値以下に設定されているので、図6に示すように、MUSICスペクトラム上における偽像は最小限に押えらて、該スペクトラム上には所望信号(複数の入射信号)のみのピークP1のみが現れる。このため、ピークP1における方位が複数の入射信号の到来方位となり、複数の入射信号の到来方位を正確に推定できる。
【0059】
(第3の実施の形態)
本発明の第3の実施の形態に係る受信装置は、センサとしてアンテナを使用し、複数のアンテナに入射された複数の入射信号の混合信号から原信号を分離して出力することを特徴とし、入射信号S1,・・・,nをブラインド信号分離の手法により求めるものである。ブラインド信号分離とは、センサの応答性や信号の性質、入射信号の予備知識なしでセンサ数と同数までの入射信号を分離するアルゴリズムである。このアルゴリズムは、ブラインド信号分離や独立成分分析(ICA:Independent Component Analysis)と呼ばれ、多くの公知文献が発表されている。
【0060】
図7に示す本発明の第3の実施の形態に係る受信装置は、図1に示す第1の実施の形態に係るセンサ配置決定装置を備えた受信装置であり、図4に示す第2の実施の形態に係る受信装置に対して、到来方位演算部47に代えて、独立成分分析部50及び出力処理部60を設けた点が異なる。図7に示す部分において、図4に示す部分と同一部分には同一符号を付し、その説明は省略する。ここでは、異なる部分についてのみ説明する。
【0061】
独立成分分析部50は、ブラインド信号分離の手法により、第1〜第4A/D変換器41〜44からの原信号(入射信号S1〜S4)が混合されてなる観測信号から原信号を分離・抽出する。この独立成分分析部50の詳細は後述する。この独立成分分析部50で抽出された信号は出力処理部60に送られる。
【0062】
出力処理部60は、第1〜第4D/A変換器61〜64から構成されている。第1〜第4D/A変換器61〜64は、独立成分分析部50からのデジタル信号をアナログ信号にそれぞれ変換し、分離信号O1〜O4としてそれぞれ出力する。
【0063】
次に、このように構成された第3の実施の形態に係る受信装置の動作を説明する。
【0064】
各アンテナ11〜14で観測される入射信号S1,・・・,nの混合信号x1,・・・,n(t)は、第1〜第4帯域制限ろ波器21〜24にそれぞれ送られる。第1〜第4帯域制限ろ波器21〜24は、混合信号x1,・・・,n(t)の所定帯域の周波数成分のみをそれぞれ通過させ、サンプリング部30に送る。
【0065】
サンプリング部30の第1中間周波数変換器31は、第1帯域制限ろ波器21からの混合信号x1(t)を入力して中間周波数の信号に変換する。詳しくは、第1中間周波数変換器31内の高周波増幅器は混合信号x1(t)を高周波増幅する。この高周波増幅された信号は、局部発振器35からの信号と混合されてそれらの和又は差の周波数が作られることにより中間周波数の信号に変換される。中間周波数の信号は中間周波数増幅器により増幅されて第1A/D変換器41に送られる。
【0066】
第1A/D変換器41は、第1中間周波数変換器31からのアナログ信号を、発振器45からのサンプリングクロックを用いてサンプリングすることによりデジタル信号に変換し、観測信号X1(t)として出力する。
【0067】
同様に、第2〜第4中間周波数変換器32〜34は、第2〜第4帯域制限ろ波器22〜24からの混合信号x2,3,4(t)をそれぞれ入力して中間周波数の信号に変換し、第2〜第4A/D変換器41にそれぞれ送る。
【0068】
第2〜第4A/D変換器42〜44は、第2〜第4中間周波数変換器32〜34からのアナログ信号を、発振器45からのサンプリングクロックを用いてサンプリングすることによりデジタル信号に変換し、観測信号X2,3,4(t)として出力する。この際、サンプリング間隔(サンプリングクロックの周波数)は、各アンテナ11〜14に到達する入射信号S1,・・・,nの時間差の影響が出ない値に調整される。このようにしてサンプリング部30で生成された観測信号X1,・・・,n(t)は、独立成分分析部50に送られる。
【0069】
独立成分分析部50は、観測信号X1,・・・,n(t)として送られてくるn個のサンプリング時系列を用いて、ブラインド信号分離のアルゴリズムにより、選択された帯域内の信号y(t)を分離・抽出する。以下、分離・抽出の手順を詳細に説明する。
【0070】
この第3の実施の形態におけるブラインド信号分離は、時間遅れなく混合された信号の分離に対して用いられる。先ず、その基本となるブラインド信号分離の問題についてここで定義する。信号源が下記式(15)のベクトルで与えられるとする。
【0071】
【数10】
但し、s(t)は、n個の入射信号であり、平均「0」であって、互いに独立であるとする。また、Tは転置を表す。
【0072】
観測は、各アンテナ11〜14で観測され帯域制限された時系列データを意味しており、
【数11】
で表すものとする。これは、1,・・・,nの各アンテナ11〜14で観測された信号であると考えることができる。一般には、アンテナの数と信号源の数とは必ずしも一致しないが、ここでは一致しているものとする。
【0073】
単純なICAの問題では、s(t)とx(t)との間に、
【数12】
なる単純な線形関係を仮定する。Aは、各アンテナ11〜14の配置と特性で決まる信号混合行列(n行×n列)の実数行列である。s(t)やAに関する知識を持たずx(t)を独立な信号成分に分離する。
【0074】
即ち、あるn×nの実数行列を求めることにより、
【数13】
で求まる互いに独立なy(t)を再構成することがICAの目的である。Bは理想的にはA−1となればよいわけだが、そうはならず順番の入れ違い(permutation)と大きさ(amplitude)の任意性は残ってしまう。
【0075】
この問題の解法の1つとして、確率分布の独立性に基づく分離法がある。各si(t)が(強)定常でガウシアン(Gaussian)でないという仮定のもとで、yi(t)が互いに独立になるようにBを求める手法がさまざまに提案されているが、それらの多くは次のようにまとめることができる。y(t)を強定常過程として、その同時分布の密度関数を、
【数14】
とすると、独立性の定義はp(yi)をp(y)のyiについての周辺分布として、
【数15】
とかける。
【0076】
同時分布と周辺分布の積との間Kullback-Leibler 発散(divergence)は、
【数16】
となる。但し、H(Y;B)は同時分布p(y)のエントロピー、H(Yi;B)は周辺分布p(yi)のエントロピーである。
【0077】
これは{Yi}(i=1,・・・,n)の相互情報量である。信号源が正規分布でないという仮定からKL(B)はp(yi)が互いに独立な場合に限り「0」となる。これらはp(x)とBによって定まるものである。
【0078】
ここで、p(y)dy=p(x)dx、p(y)=p(x)/|B|(|B|はBの行列式)であることに注意すると、
【数17】
となる。
【0079】
一方、周辺分布のエントロピーは、
【数18】
である。よって、
【数19】
となり、
【数20】
のようにすれば最急降下法として正しいBを求めることができる。上記式(25)の中で問題となるのは逆行列(BT)−1を計算している点である。
【0080】
収束性に関しては、これにいかなる正定値行列を掛けても構わないことから、BTBを掛ければ(これは正則な行列の多様体上でのリーマン(Rieman)計量に対応している)、
【数21】
が新たな学習則となる。定常性の仮定よりp(s)、p(x)とp(y)は時間的に独立である。この仮定のもと、アンサンブル平均を時間平均に置き換えることができる。
【0081】
【数22】
したがって、ηを正の定数とし、データが観測される毎に下記式(28)に従ってパラメータを更新すればBtが得られる。
【0082】
【数23】
ここで、当然問題になるのは、上記式(24)のp(yi)或いはψ(y)をいかに定義するかである。通常、これはパラメトリックな非線形関数や統計的な展開法が用いられる。
【0083】
大雑把な考え方を示すと、もしp(yi)が正規分布ならはψ(y)は線形関数となる。一方、正規分布より裾が”重い”場合(sub-Gaussian)多項式などで近似するのがよく、正規分布より裾が”軽い”場合(super-Gaussian)シグモイド(Sigmoid)関数などで近似するのがよいとされている。音声信号などは裾が”軽い”ので、シグモイド関数などがうまく働く。
【0084】
以上のようにして独立成分分析部50で分離・抽出された信号y1,・・・,n(t)は出力処理部60に送られる。出力処理部60では、独立成分分析部50からのデジタル信号として送られてくる信号y1,・・・,n(t)をアナログ信号に変換する。より詳しくは、出力処理部60から出力される信号y1(t)は、第1D/A変換器61でデジタル信号に変換され、分離信号O1として外部に送出される。同様に、信号y2,3,4(t)は、第2〜第4D/A変換器62〜64でデジタル信号にそれぞれ変換され、分離信号O2〜O4として外部にそれぞれ送出される。
【0085】
このように、第3の実施の形態に係る受信装置によれば、ベクトル生成部27、相関量算出部28及びアンテナ配置決定部29を設けたので、第1の実施の形態に係るセンサ配置決定装置の効果と同様な効果が得られる。即ち、入射信号の到来方位を推定する際に発生する偽像を最小限に押えることができる。
【0086】
また、受信した各アンテナの受信信号を各アンテナの信号の到達時間に差がなくなる狭帯域信号にまで帯域制限を行い、その狭帯域に帯域制限された各アンテナの信号を用いて独立成分分析を行うようにしたので、各アンテナの信号到達時間の差が生じる条件下でも一般的なブラインド信号分離の手法によりアンテナからの観測信号である混合信号中に含まれる原信号を分離・抽出できる。
【0087】
なお、上述した実施の形態では、センサとして、アンテナを用いた例を説明したが本発明はこれに限定されるものではない。例えば、センサとして、音声を検知するマイクロフォン、生体の種々の状態を検知するセンサ等を用いることができる。また、本発明は、第2の実施の形態に係る受信装置と第3の実施の形態に係る受信装置とを組み合わせた受信装置にも適用可能であるのは勿論である。
【0088】
【発明の効果】
以上詳述したように、本発明によれば、ステアリングベクトルに基づいて算出された自己内相関量の内の相互相関量が所定値以下になるように複数のセンサの配置を決定するので、入射信号の到来方位を推定する際に発生する偽像を最小限に押えることができ、これによって到来方位を正確に推定できるセンサ配置決定方法及びその装置並びに受信装置を提供できる。
【図面の簡単な説明】
【図1】本発明の第1の実施の形態に係るセンサ配置決定装置の構成を示すブロック図である。
【図2】試験電波を用いてステアリングベクトルを測定する例を示す図である。
【図3】本発明の第1の実施の形態に係るセンサ配置決定方法を説明するフローチャートである。
【図4】本発明の第2の実施の形態に係る受信装置の構成を示すブロック図である。
【図5】MUSICによる到来方位推定を説明するための図である。
【図6】到来角に対するMUSICスペクトラムを示す図である。
【図7】本発明の第3の実施の形態に係る受信装置の構成を示すブロック図である。
【符号の説明】
11〜14 第1〜第4アンテナ
21〜24 第1〜第4帯域制限ろ波器
25 帯域制限制御器
27 ベクトル生成部
28 相関量算出部
29 アンテナ配置決定部
31〜34 第1〜第4中間周波数変換器
35 局部発振器
41〜44 第1〜第4A/D変換器
45 発振器
47 到来方位演算部
50 独立成分分析部
60 出力処理部
61〜64 第1〜第4D/A変換器
Claims (10)
- 同一帯域内に複数の入射信号が混信している場合に、同一帯域内の複数の入射信号が同時に入射される複数のセンサを配置し、
前記各センサにおける前記複数の入射信号の位相応答を行列で表したステアリングベクトルを生成し、
生成されたステアリングベクトルに基づいて前記ステアリングベクトルの自己相関量と相互相関量とからなる自己内相関量を算出し、
算出された前記相互相関量が所定値以下か否かを判定し、前記相互相関量が所定値以下になっていない場合には前記複数のセンサの配置を変更することを特徴とするセンサ配置決定方法。 - 前記ステアリングベクトルは、仰角を固定し且つ試験電波を各到来方向から前記各センサに入射したときに前記各センサにおいて所定の角度毎に測定したデータであることを特徴とする請求項1に記載のセンサ配置決定方法。
- 前記生成されたステアリングベクトルに基づいて前記ステアリングベクトルの複素共役転置行列を求め、前記複素共役転置行列と前記ステアリングベクトルとを演算することにより前記ステアリングベクトルの前記自己内相関量を求めることを特徴とする請求項1又は請求項2に記載のセンサ配置決定方法。
- 同一帯域内に複数の入射信号が混信している場合に、同一帯域内の複数の入射信号が同時に入射される複数のセンサと、
前記各センサにおける前記複数の入射信号の位相応答を行列で表したステアリングベクトルを生成するベクトル生成部と、
前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの自己相関量と相互相関量とからなる自己内相関量を算出する相関量算出部と、
前記相関量算出部で算出された前記相互相関量が所定値以下か否かを判定し、前記相互相関量が所定値以下になっていない場合には前記複数のセンサの配置を変更するセンサ配置決定部と、
を備えたことを特徴とするセンサ配置決定装置。 - 前記ベクトル生成部は、仰角を固定し且つ試験電波を各到来方向から前記各センサに入射したときに前記各センサにおいて所定の角度毎に測定したデータを前記ステアリングベクトルとすることを特徴とする請求項4に記載のセンサ配置決定装置。
- 前記相関量算出部は、前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの複素共役転置行列を求め、前記複素共役転置行列と前記ステアリングベクトルとを演算することにより前記ステアリングベクトルの前記自己内相関量を求めることを特徴とする請求項4又は請求項5に記載のセンサ配置決定装置。
- 同一帯域内に複数の入射信号が混信している場合に、同一帯域内の複数の入射信号が同時に入射される複数のセンサと、
前記各センサにおける前記複数の入射信号の位相応答を行列で表したステアリングベクトルを生成するベクトル生成部と、
前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの自己相関量と相互相関量とからなる自己内相関量を算出する相関量算出部と、
前記相関量算出部で算出された前記相互相関量が所定値以下か否かを判定し、前記相互相関量が所定値以下になっていない場合には前記複数のセンサの配置を変更するセンサ配置決定部と、
前記センサ配置決定部で前記複数のセンサの配置が決定されたときにおける、前記相関量算出部で算出された前記自己内相関量に基づいて到来角に対するスペクトラムを求め、該スペクトラムに基づいて前記複数の入射信号の到来方位を求める到来方位演算部と、
を備えたことを特徴とする受信装置。 - 同一帯域内に複数の入射信号が混信している場合に、同一帯域内の複数の入射信号が同時に入射される複数のセンサと、
前記各センサにおける前記複数の入射信号の位相応答を行列で表したステアリングベクトルを生成するベクトル生成部と、
前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの自己相関量と相互相関量とからなる自己内相関量を算出する相関量算出部と、
前記相関量算出部で算出された前記相互相関量が所定値以下か否かを判定し、前記相互相関量が所定値以下になっていない場合には前記複数のセンサの配置を変更するセンサ配置決定部と、
前記複数のセンサからの前記複数の入射信号が混合された複数の混合信号を所定の狭帯域にそれぞれ制限して出力する複数の帯域制限ろ波器と、
前記複数の帯域制限ろ波器からの狭帯域に制限された複数の信号に基づいて、独立成分分析の手法により前記複数の入射信号を分離・抽出する独立成分分析部と、
を備えたことを特徴とする受信装置。 - 前記ベクトル生成部は、仰角を固定し且つ試験電波を各到来方向から前記各センサに入射したときに前記各センサにおいて所定の角度毎に測定したデータを前記ステアリングベクトルとすることを特徴とする請求項7又は請求項8に記載の受信装置。
- 前記相関量算出部は、前記ベクトル生成部で生成されたステアリングベクトルに基づいて前記ステアリングベクトルの複素共役転置行列を求め、前記複素共役転置行列と前記ステアリングベクトルとを演算することにより前記ステアリングベクトルの前記自己内相関量を求めることを特徴とする請求項7乃至請求項9のいずれか1項に記載の受信装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001400964A JP3973902B2 (ja) | 2001-12-28 | 2001-12-28 | センサ配置決定方法及びその装置並びに受信装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001400964A JP3973902B2 (ja) | 2001-12-28 | 2001-12-28 | センサ配置決定方法及びその装置並びに受信装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2003202368A JP2003202368A (ja) | 2003-07-18 |
JP3973902B2 true JP3973902B2 (ja) | 2007-09-12 |
Family
ID=27640009
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001400964A Expired - Fee Related JP3973902B2 (ja) | 2001-12-28 | 2001-12-28 | センサ配置決定方法及びその装置並びに受信装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3973902B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4920270B2 (ja) * | 2006-03-06 | 2012-04-18 | Kddi株式会社 | 信号到来方向推定装置及び方法、並びに信号分離装置及び方法、コンピュータプログラム |
JP5519132B2 (ja) * | 2008-07-28 | 2014-06-11 | 株式会社デンソー | レーダ装置 |
-
2001
- 2001-12-28 JP JP2001400964A patent/JP3973902B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2003202368A (ja) | 2003-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chargé et al. | A non-circular sources direction finding method using polynomial rooting | |
JP6395677B2 (ja) | 到来方向推定装置 | |
WO2006067869A1 (ja) | 到来方向推定装置及びプログラム | |
JP2009210410A (ja) | 探知測距装置および探知測距方法 | |
KR102183439B1 (ko) | Beamspace MUSIC과 TMA를 결합한 도래각 추정 방법 및 장치 | |
JP2009162689A (ja) | 電子走査型レーダ装置、受信波方向推定方法及び受信波方向推定プログラム | |
CN111308424A (zh) | 一种基于相加求和与music联合算法的变电站设备可听声源定位方法 | |
JP3766006B2 (ja) | 受信装置 | |
US6446025B1 (en) | Multiple propagation wave parameter measuring method and apparatus and machine-readable recording medium recording multiple propagation wave parameter measuring program | |
JP5022943B2 (ja) | 方向測定装置 | |
JP5553980B2 (ja) | 電波方向探知装置およびビーム形成装置 | |
JP3973902B2 (ja) | センサ配置決定方法及びその装置並びに受信装置 | |
JP3867069B2 (ja) | 電波到来方向探知方法及び装置 | |
Kulhandjian et al. | 2-D DOA estimation of coherent wideband signals with auxiliary-vector basis | |
JP3959357B2 (ja) | 受信システム | |
Pal et al. | Gridless methods for underdetermined source estimation | |
KR20190134893A (ko) | 안테나 어레이 외삽을 이용한 레이더 수신신호의 도착방향 추정 방법 및 장치 | |
JP3977751B2 (ja) | センサ配置最適化方法及びその装置並びに到来方位測定装置 | |
JP2003185726A (ja) | 受信装置及び検出装置 | |
Bloemendal et al. | Blind source extraction for a combined fixed and wireless sensor network | |
Jiang et al. | A new source number estimation method based on the beam eigenvalue | |
RU2410707C2 (ru) | Способ поляризационно-независимого обнаружения и локализации широкополосных радиосигналов | |
JP2005049123A (ja) | 信号処理装置 | |
JP2008160448A (ja) | 受信装置及び混信分離方法 | |
JP4576742B2 (ja) | 送受信周波数分割多重無線装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040513 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20060220 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060307 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070313 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070508 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070605 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070613 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100622 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100622 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110622 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120622 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120622 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130622 Year of fee payment: 6 |
|
LAPS | Cancellation because of no payment of annual fees |