JP3682398B2 - Wireless direction measuring device - Google Patents

Wireless direction measuring device Download PDF

Info

Publication number
JP3682398B2
JP3682398B2 JP2000087082A JP2000087082A JP3682398B2 JP 3682398 B2 JP3682398 B2 JP 3682398B2 JP 2000087082 A JP2000087082 A JP 2000087082A JP 2000087082 A JP2000087082 A JP 2000087082A JP 3682398 B2 JP3682398 B2 JP 3682398B2
Authority
JP
Japan
Prior art keywords
phase difference
evaluation function
antenna element
antenna
receivers
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 - Lifetime
Application number
JP2000087082A
Other languages
Japanese (ja)
Other versions
JP2001272447A (en
Inventor
慎太郎 荒田
Original Assignee
株式会社光電製作所
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 株式会社光電製作所 filed Critical 株式会社光電製作所
Priority to JP2000087082A priority Critical patent/JP3682398B2/en
Publication of JP2001272447A publication Critical patent/JP2001272447A/en
Application granted granted Critical
Publication of JP3682398B2 publication Critical patent/JP3682398B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Description

【0001】
【発明の属する技術分野】
本発明は、電波監視などの分野で利用される無線方位測定装置に関するものである。
【0002】
【従来の技術】
電波監視などの分野では、電波の到来方向を探知する無線方位測定装置が使用されている。このような無線方位測定装置の一つの典型例であるインタ−フェロメ−タ方式の無線方位測定装置の一例は、3個以上の複数のアンテナ素子を所定の配列で配置し、各アンテナ素子が受信した電波の位相差を検出することにより、この電波の到来方向を探知するように構成されている。
【0003】
すなわち、このようなインタ−フェロメ−タ方式の無線方位測定装置の典型的な一例は、図2に示すように、n個のアンテナ素子A1〜Anと、n個の受信機R1〜Rnと、n個のサンプリング部SP1〜SPnと、位相差検出部PDと、方位算定部DDと、較正信号発生部CSと、n個の信号切替部S1〜Snとから構成されている。アンテナ素子A1〜Anから出力される受信信号x1〜xnが対応の受信機R1〜Rnとサンプリング部SP1〜SPnで処理され、複素受信電圧ベクトルX=〔X1,X2,X3・・・・Xn〕t が生成される。
【0004】
得られた複素受信電圧ベクトルX=〔X1,X2,X3・・・・Xn〕t から任意のアンテナ素子対Ai , Aj に対応する複素受信電圧Xi , Xj を選択すれば、このアンテナ素子対の位相差Yijは(1)式のように表現される。
ij=<Xi ・Xj * > /ABS[<Xi ・Xj * >] ・・・・(1)
ただし、<z>はzの算定回数にわたる平均値を表し、ABS[z] はzの絶対値を表わす。
【0005】
また、アレーアンテナのモードベクトルをA(θ)=〔A1(θ),A2(θ)・・・・An(θ)〕t と表すと、任意のアンテナ素子対Ai , Aj に対応する位相差Bij (θ) は次式で与えられる。
ij (θ) =Ai (θ) ・Aj (θ) * /ABS[Ai (θ) ・Aj (θ) * ] ・・・(2)
【0006】
ただし、上記アレーアンテナのモードベクトルの各要素は、電波の到来方向を既知の値θとした場合に、実測や計算によって得られる各アンテナ素子から出力される複素受信電圧である。
【0007】
上記(1)式と(2)式とから次の(3)式で与えられる評価関数PIF(θ)が定義される。
IF(θ)=Σ i Σj ABS[Yij−Bij (θ)] ・・・・(3)
ただし、記号Σ i Σj z は(i,j) の一部又は全ての組合せについてのzの総和を意味する。そして、(3)式の評価関数PIF(θ)を最小にするθの値が電波の到来方位として検出される。
【0008】
図2の無線方位測定装置では、受信機R1〜Rnから複素受信電圧X1〜Xnが出力され、これらについて、サンプリング部SP1〜SPnにおいてサンプリングとディジタル化が行われ、この複素受信電圧について、位相差検出部PDと方位算定部DDにおいて、(1)式〜(3)式の演算が行われる。
【0009】
実際には、各受信機や、各サンプリング部の特性や接続ケーブルの長さのばらつきなどによって位相差が発生し、これが検出方位に測定誤差を発生させる。このような位相差を較正するために、図2に示すように、各アンテナ素子の後に切替部S1〜Snを設置し、対応のアンテナ素子A1〜Anから出力される受信信号の代わりに、較正信号発生部CSから出力される較正信号を信号経路の切替えによって各受信機R1〜Rnに入力させ、この場合の位相差(Yij0 を算定してメモリに保存しておく。そして、引き続いて行う実測によって得た実測値から、保存中の位相差(Yij0 を減算することにより、実測値に対する較正が行われる。
【0010】
図2に示すように、実測値の較正のために必要となる構成要素なども含めると無線方位測定装置全体のハードウエア量が増大し、測定装置全体が大型かつ高価になる。そこで、ハードウエア量の圧縮を図るために、図3に示すような構成が採用さている。
【0011】
すなわち、受信機とサンプリング部を2系統だけ設置し、アンテナ素子対選択部PSによって各アンテナ素子の中からアンテナ素子対を選択して2系統の受信機とサンプリング部とに接続することにより、限られた台数の受信機とサンプリング部とを時分割式に共用し、これによって、高価な受信機とサンプリング部の設置台数を節減している。ただし、この時分割共用化の構成では、全アンテナ素子対の選択にある程度の時間がかかるので、その間到来電波が存続していることが必要である。
【0012】
【発明が解決しようとする課題】
従来の無線方位測定装置では、図2の構成から図3の構成のように変更を加えることにより、受信機とサンプリング部のハードウエア量はかなり低減される。しかしながら、較正信号発生部や、切替部など較正用のハードウエアは依然として残り、装置全体の製造費用のかなりの部分を占めるという問題がある。
【0013】
特に、ビルの鉄塔上などに設置される切替部及びアンテナ素子選択部と、屋上や屋内に設置される受信機との間は、場合によっては百メートルにも及ぶ長大なケーブルで接続される。また、受信機と同一箇所に設置される較正信号発生部からの較正信号を切替部に供給するために、同一の長さのケーブルを敷設しなければならない。ところが、このようなケーブルは、一般に広帯域特性が要求されるために高価なものとなり、装置全体の費用を増大させるという問題がある。従って、本発明の一つの目的は、較正用のハードウエア量を低減することにより、装置全体の製造費用を低減することにある。
【0014】
また、較正対象の位相差(Yij0 は信号の周波数や、受信信号のレベルや、温度変化などの環境に依存して変化する。このため、較正信号について予め算定し保存しておくデータ量を低減したり、較正に要する時間を短縮するために、受信機の周波数特性を均一化したり、広い入力レベルにわたって入出力特性の直線化を図ったり、温度変化に対する安定化などを図る必要があり、この結果、受信機の製造費用が高価になるという問題がある。従って、本発明の他の目的は、製造・維持費用の高価な特別な受信機を必要としない無線方位測定装置を提供することにある。
【0015】
【課題を解決するための手段】
上記従来技術の課題を解決する本発明の無線方位測定装置は、複数のアンテナ素子から成るアレーアンテナと、前記各アンテナ素子から出力される受信信号が選択的に入力される2個の受信機と、前記各アンテナ素子の中から1対(Ai ,Aj )を順次選択して、それぞれの受信信号を前記2個の受信機に選択的に入力させるアンテナ素子対選択部と、前記2個の受信機の出力から、前記選択されたアンテナ素子対(Ai ,Aj )の受信信号の位相差Yijを算定すると共に、前記アレーアンテナのモードベクトルから得られる前記選択されたアンテナ素子対(Ai ,Aj )の出力の位相差Bij(θ)との差に基づく評価関数P IF (θ)を算定し、この評価関数P IF (θ)から入射角θを測定方位として得る演算部とを備えたインタ−フェロメ−タ方式の無線方位測定装置である。
【0016】
そして、上記演算部は、前記順次選択され、算定されたアンテナ素子対の受信信号の位相差 ij のうち任意の一つYpqを選択し、前記位相差Yijの代わりにこの位相差Yijと上記位相差Ypqとの差 ij (=Y ij ・Y pq * を用いると共に、上記位相差Bij(θ)の代わりにこの位相差Bij(θ)と上記位相差Ypqに対応するBpq(θ)との差 ij (θ)(=B ij (θ)・B pq * (θ))を用いて評価関数P IF (θ)の代わりに評価関数PP IF (θ)=Σ i Σ j ABS[ ij −C ij (θ) ] から入射角θを算定する手段を備えることにより、受信系の較正を不要にするように構成されている。
【0017】
【発明の実施の形態】
本発明の好適な実施の形態によれば、上記演算部は、上記評価関数の算定に際し、上記位相差Ypqを複数個選択して、それぞれの場合の評価関数を算定し、これらを平均化したものを最終的な評価関数とすることにより、雑音の影響を一層軽減して検出精度を高めるように構成されている。
【0018】
【実施例】
図1は、本発明の一実施例の無線方位測定装置の構成を示す機能ブロック図である。この無線方位測定装置は、n個のアンテナ素子A1〜Anと、アンテナ素子対選択部PSと、2個の受信機R1,R2と、同じく2個のサンプリング部SP1,SP2と、位相差検出部PDと、方位算定部DDとを備えている。
【0019】
n個のアンテナ素子A1〜Anからは、受信信号x1〜xnが出力され、これらの受信信号のうちの1対xi ,xj がアンテナ素子対選択部PSで選択され、受信信号xi は受信機R1に、受信信号xj は受信機R2に供給される。さらに受信機R1,R2の出力がそれぞれサンプリング部SP1,SP2でディジタル信号に変換され、得られた複素電圧X1iとX2jが後段の位相差検出部PDに供給される。
【0020】
上記1対のアンテナ素子は、アレーアンテナにおけるアンテナ素子の配置と受信周波数においてアンテナ素子間に発生する位相差を考慮し、隣接するものどうしがペアリングされるだけでなく、最遠方のものどうしがペアリングされたりする。また、アンテナ素子対選択部の製造の容易さなども考慮され、全ての組合せが実現される場合だけでなく、一部の組合せだけが実現される場合もある。
【0021】
位相差検出部PDは、前段のサンプリング部SP1,SP2のそれぞれから供給された複素受信電圧X1i,X2jから、次の(4)式で表される位相差Yijを算定する。
ij=<X1i・X2j * >/ABS[<X1i・X2j * >] ・・・・(4)
さらに、位相差検出部PDは、上記アンテナ素子対の位相差のうち任意のアンテナ素子対間の位相差Ypqを選択し、以下の(5)式で与えられるYijとYpq との位相差Zijを算定し、この算定結果を後段の方位算定部DDに渡す。
ij=Yij・Ypq * ・・・・(5)
【0022】
さらに、位相差検出部PDは、上記(2)式に従って算定済みのモードベクトル間の位相差Bij (θ) と、上記任意のアンテナ素子対Ypqに対応する位相差Bpq(θ)との位相差Cij(θ)を算定し、この算定結果を後段の方位算定部DDに渡す。
ij (θ) =Bij (θ) ・Bpq *(θ) ・・・・(6)
【0023】
方位算定部DDは、前述した(1)式で与えられる位相差Yijと、(2)式で与えられる位相差Bij (θ) の代わりに、上記(5)式で与えられるZijと(6)式で与えられるCij (θ) とから評価関数PPIF(θ)を算定する。
PPIF(θ)=Σ i Σj ABS[Zij−Cij (θ)] ・・・・(7)
そして、方位算定部DDは、(7)式の評価関数PPIF(θ)を最小にするθの値を電波の到来方位として検出する。
【0024】
上記(7)式では、受信系に含まれる位相遅延量が相殺されて除去され、このため、自動的に較正が行われたと同一の結果となる。すなわち、受信機R1、サンプリング部SP1及び信号線とから成る受信系1内の位相遅延量をαとし、受信機R2、サンプリング部SP2及び信号線から成る受信系2内の位相遅延量をβとする。この場合、上記実測された複素受信電圧X1iとX2jは、次式で与えられる。
1i=V1i exp(j( φi +α)) ・・・・(8)
2j=V2j exp(j( φj +β)) ・・・・(9)
【0025】
(8)式と(9) 式とを(4) 式、に代入すると、

Figure 0003682398
【0026】
(10) 式と(11)式を(5) 式に代入すると、
Figure 0003682398
【0027】
(12)式から明らかなように、受信系1,2の各位相遅延量α,βは、互いに相殺されて除去されることにより、Zij中には含まれなくなる。このように、受信系1,2の各位相遅延量α,βが互いに相殺されて除去されるため、従来の無線方位測定装置においてこのα、βを予め測定しておくために必要であった較正用のハードウエアが一切不要になる。
【0028】
(7)式では、アンテナ素子対 Ap , A qを任意の対に固定した。しかしながら、このアンテナ素子対を固定することなく、測定に使用するアンテナ素子対からアンテナ素子対 Ap , A qを複数選択しながら評価関数を算定し、これらの評価関数を平均したものを最終的な評価関数として算定することにより、雑音の影響を更に軽減して測定精度を高めることができる。
【0029】
すなわち、測定に使用するアンテナ素子対から選択されたアンテナ素子対 Ap , A qごとに(7)式の評価関数PPIF(θ)が算定され、これらの総和が新たな評価関数
PPPIF(θ)= pΣq iΣj ABS[Zij−Cij (θ)] ・・・・(13)
として算定される。そして、この新たな評価関数PPPIF(θ)を最小とするθが到来方位として検出される。この新たな評価関数においては、特定のアンテナ素子に発生する雑音の影響がアンテナ素子間の平均化によって軽減され、方位検出精度が一層向上する。
【0030】
以下では、図1に示した実施例について行った計算機シミュレーションによる確認の結果を、図4乃至図7を参照しながら説明する。
【0031】
図4は、この計算機シミュレーションで想定した各アンテナ素子の配置を示す平面図である。5個のアンテナ素子1〜5が水平面内に想定した円周上に等間隔で配置されている。到来電波の周波数は、対応の半波長(λ/2)が、最隣接のアンテナ素子間の直線距離に等しくなる値とする。また、アンテナ素子の対は、素子1と3、2と4、3と5、4と1、5と2という一つ跳びに組合された5対とする。また、受信系1と2との間に250°の位相差が設定されている。さらに、受信電波は方位150 °から到来し、信号のS/Nは10dBで、データの平均回数は20回とする。
【0032】
図5は、(3)式のPIF(θ)について、LOG10(1/PIF(θ))の計算結果をプロットしたグラフである。図6は、(7)式のPPIF(θ)について、LOG10(1/PPIF(θ))の計算結果をプロットしたグラフである。図7は(13)式のPPPIF(θ)について、LOG10(1/PPPIF(θ))の計算結果をプロットしたグラフである。
【0033】
図5のグラフは、図2や図3に示すように、相当のハードウエアを使用して理想的な較正を行った場合の結果である。これに対して、この実施例の無線方位測定装置について行ったシミュレーション結果を示す図6と図7のグラフは、較正を全く行わない場合の結果であるが、これらは図5の場合と遜色ない結果となっている。
【0034】
このように、実測開始前の較正が不要となったため、較正用のハードウエアが不要になる。同時に、従来装置のように、メモリに保存しておく較正用データの量を低減したり、較正に要する時間を短縮するために、受信機R1〜Rnの周波数特性を均一化したり、入出力特性の直線領域を拡大したり、温度特性を安定にしたりすることの必要性が軽減される。この結果、受信機R1〜Rnとして、特別な受信機を必要とせず、普通の周波数特性や温度安定性を有する簡易・安価なものを使用できる。
【0035】
以上、説明の便宜上、位相差検出部と、方位算定部が個別に構成される場合を例示した。しかしながら、これら各部をマイクロコンピュータなどによって実現される一体の処理部として構成することもできる。
【0036】
【発明の効果】
以上詳細に説明したように、本発明の無線方位測定装置は、二つの受信系から得られる受信電圧の位相差Yijと、モードベクトルの位相差Bij(θ)の代わりに、Yijpq * とBij(θ)Bpq(θ)* を用いて評価関数P(θ)を算定する構成であるから、この評価関数の中で、各受信系に含まれる位相差が互いに相殺され、除去される。この結果、各受信系に含まれる位相差を較正するための、較正信号発生部や、高性能で高価なケーブルや、切替部などのハードウエアが一切不要となり、装置全体の製造費用が大幅低減される。
【0037】
さらに、本発明の無線方位測定装置では、受信機として普通の周波数特性、入出力の直線性や温度安定性を有する簡易・安価なものを使用でき、この点においても、測定装置全体としての製造費用が低減されるという効果も奏される。
【図面の簡単な説明】
【図1】本発明の一実施例の無線方位測定装置の構成を示す機能ブロック図である。
【図2】従来の無線方位測定装置の構成の一例を示す機能ブロック図である。
【図3】従来の無線方位測定装置の構成の他の一例を示す機能ブロック図である。
【図4】本発明の効果を確認するための電子計算機シミュレーションに使用したアンテナ素子の配列を示す平面図である。
【図5】上記電子計算機シミュレーションによって得られた従来の評価関数である。
【図6】上記電子計算機シミュレーションによって得られた本発明の一実施例による評価関数である。
【図7】上記電子計算機シミュレーションによって得られた本発明の他の実施例による評価関数である。
【符号の説明】
A1〜An アンテナ素子
R1〜Rn 受信機
SP1 〜SPn サンプリング部
PS アンテナ素子対選択部
PD 位相差検出部
DD 方位算定部
CS 較正信号発生部
S1〜Sn 切替部[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a radio direction measuring device used in fields such as radio wave monitoring.
[0002]
[Prior art]
In fields such as radio wave monitoring, wireless azimuth measuring devices that detect the direction of arrival of radio waves are used. One example of such a radio direction measuring apparatus is an interferometer type radio direction measuring apparatus, in which three or more antenna elements are arranged in a predetermined arrangement, and each antenna element receives a signal. By detecting the phase difference between the radio waves, the direction of arrival of the radio waves is detected.
[0003]
That is, as shown in FIG. 2, a typical example of such an interferometer-type radio direction measuring device includes n antenna elements A1 to An, n receivers R1 to Rn, It consists of n sampling units SP1 to SPn, a phase difference detection unit PD, an orientation calculation unit DD, a calibration signal generation unit CS, and n signal switching units S1 to Sn. Received signals x1 to xn output from antenna elements A1 to An are processed by corresponding receivers R1 to Rn and sampling units SP1 to SPn, and complex received voltage vector X = [X1, X2, X3... Xn]. t is generated.
[0004]
The resulting complex reception voltage vector X = [X1, X2, X3 ···· Xn] Any antenna element pair from t A i, the received complex voltage X i corresponding to A j, by selecting the X j, the The phase difference Y ij of the antenna element pair is expressed as in equation (1).
Y ij = <X i · X j *> / ABS [<X i · X j *>] ···· (1)
Here, <z> represents an average value over the number of times z is calculated, and ABS [z] represents an absolute value of z.
[0005]
Further, when the mode vector of the array antenna is expressed as A (θ) = [A1 (θ), A2 (θ)... An (θ)] t , it corresponds to an arbitrary antenna element pair A i , A j . The phase difference B ij (θ) is given by the following equation.
B ij (θ) = A i (θ) A j (θ) * / ABS [A i (θ) A j (θ) * ] (2)
[0006]
However, each element of the mode vector of the array antenna is a complex reception voltage output from each antenna element obtained by actual measurement or calculation when the arrival direction of the radio wave is a known value θ.
[0007]
The evaluation function P IF (θ) given by the following equation (3) is defined from the above equations (1) and (2).
P IF (θ) = Σ i Σ j ABS [Y ij −B ij (θ)] (3)
Here, the symbol Σ i Σ j z means the total sum of z for some or all combinations of (i, j) . Then, the value of θ that minimizes the evaluation function P IF (θ) of the expression (3) is detected as the arrival direction of the radio wave.
[0008]
In the radio direction measuring apparatus of FIG. 2, complex reception voltages X1 to Xn are output from the receivers R1 to Rn, and these are sampled and digitized by the sampling units SP1 to SPn. In the detection unit PD and the direction calculation unit DD, calculations of formulas (1) to (3) are performed.
[0009]
Actually, a phase difference occurs due to variations in the characteristics of each receiver, each sampling unit, and the length of the connection cable, and this causes a measurement error in the detection direction. In order to calibrate such a phase difference, as shown in FIG. 2, switching units S1 to Sn are installed after each antenna element, and calibration is performed instead of the reception signals output from the corresponding antenna elements A1 to An. The calibration signal output from the signal generator CS is input to each of the receivers R1 to Rn by switching the signal path, and the phase difference (Y ij ) 0 in this case is calculated and stored in the memory. Then, the measured value is calibrated by subtracting the phase difference (Y ij ) 0 during storage from the actually measured value obtained by the subsequent actual measurement.
[0010]
As shown in FIG. 2, when the components necessary for calibration of the actual measurement values are included, the amount of hardware of the entire wireless azimuth measuring apparatus increases, and the entire measuring apparatus becomes large and expensive. Therefore, in order to compress the amount of hardware, a configuration as shown in FIG. 3 is adopted.
[0011]
That is, only two systems of the receiver and the sampling unit are installed, and the antenna element pair selection unit PS selects the antenna element pair from each antenna element and connects it to the two systems of the receiver and the sampling unit. The same number of receivers and sampling units are shared in a time-sharing manner, thereby reducing the number of expensive receivers and sampling units installed. However, in this time-sharing shared configuration, it takes a certain amount of time to select all antenna element pairs, and thus it is necessary that the incoming radio waves continue during that time.
[0012]
[Problems to be solved by the invention]
In the conventional radio direction measuring apparatus, the hardware amount of the receiver and the sampling unit is considerably reduced by changing the configuration of FIG. 2 to the configuration of FIG. However, the calibration hardware such as the calibration signal generator and the switching unit still remains, and there is a problem that it accounts for a considerable part of the manufacturing cost of the entire apparatus.
[0013]
In particular, a switching unit and an antenna element selection unit installed on a steel tower of a building and a receiver installed on a rooftop or indoors are connected with a long cable of up to 100 meters in some cases. Moreover, in order to supply the calibration signal from the calibration signal generation unit installed at the same location as the receiver to the switching unit, a cable having the same length must be laid. However, such a cable is generally expensive because it requires a broadband characteristic, and there is a problem that the cost of the entire apparatus is increased. Accordingly, one object of the present invention is to reduce the manufacturing cost of the entire apparatus by reducing the amount of hardware for calibration.
[0014]
The phase difference (Y ij ) 0 to be calibrated changes depending on the environment such as the frequency of the signal, the level of the received signal, and the temperature change. For this reason, to reduce the amount of data that is calculated and stored in advance for the calibration signal, to reduce the time required for calibration, the receiver frequency characteristics are made uniform, and the input / output characteristics are linearized over a wide input level. There is a problem that the manufacturing cost of the receiver becomes high as a result. Accordingly, another object of the present invention is to provide a radio direction measuring apparatus that does not require a special receiver that is expensive to manufacture and maintain.
[0015]
[Means for Solving the Problems]
A radio direction measuring apparatus of the present invention that solves the above-described problems of the prior art includes an array antenna including a plurality of antenna elements, and two receivers to which reception signals output from the antenna elements are selectively input. A pair (A i , A j ) are sequentially selected from the antenna elements, and an antenna element pair selector for selectively inputting each received signal to the two receivers; The phase difference Y ij of the received signal of the selected antenna element pair (A i , A j ) is calculated from the output of the receiver of the selected antenna element pair, and the selected antenna element pair obtained from the mode vector of the array antenna An evaluation function P IF (θ) is calculated based on the difference between the output of (A i , A j ) and the phase difference B ij (θ), and the incident angle θ is obtained as a measurement direction from the evaluation function P IF (θ). Interferometric wireless method with computing unit It is a measuring device.
[0016]
Then, the calculation unit selects any one Y pq of the phase differences Y ij of the received signals of the antenna element pairs that are sequentially selected and calculated, and replaces the phase difference Y ij with this phase difference Y ij. ij and with use of the difference Z ij (= Y ij · Y pq *) between the phase difference Y pq, the phase difference between the phase difference B ij (theta) in place of the retardation B ij (θ) Y pq the difference C ij of the corresponding B pq (theta) to (θ) (= B ij ( θ) · B pq * (θ)) in place of the evaluation function P IF (theta) using the evaluation function PP IF (theta ) = by comprising means for calculating the incident angle theta from Σ i Σ j ABS [Z ij -C ij (θ)], and is configured to eliminate the need for calibration of the receiving system.
[0017]
DETAILED DESCRIPTION OF THE INVENTION
According to a preferred embodiment of the present invention, the calculation unit selects a plurality of the phase differences Y pq when calculating the evaluation function, calculates the evaluation function in each case, and averages these By using this as the final evaluation function, the influence of noise is further reduced and the detection accuracy is increased.
[0018]
【Example】
FIG. 1 is a functional block diagram showing a configuration of a radio direction measuring apparatus according to an embodiment of the present invention. This radio direction measuring device includes n antenna elements A1 to An, an antenna element pair selection unit PS, two receivers R1 and R2, two sampling units SP1 and SP2, and a phase difference detection unit. PD and azimuth | direction calculation part DD are provided.
[0019]
From the n antenna elements Al-An, the received signal x1~xn is outputted, a pair x i of the received signals, x j is selected by the antenna element pair selection unit PS, the received signal x i is The received signal x j is supplied to the receiver R1 to the receiver R1. Further, the outputs of the receivers R1 and R2 are converted into digital signals by the sampling units SP1 and SP2, respectively, and the obtained complex voltages X 1i and X 2j are supplied to the subsequent phase difference detection unit PD.
[0020]
In consideration of the arrangement of the antenna elements in the array antenna and the phase difference generated between the antenna elements at the reception frequency, the pair of antenna elements is not only paired between adjacent elements but also between the farthest ones. It is paired. Considering the ease of manufacturing the antenna element pair selection unit, etc., not only all combinations may be realized, but also some combinations may be realized.
[0021]
The phase difference detection unit PD calculates a phase difference Y ij represented by the following equation (4) from the complex reception voltages X 1i and X 2j supplied from the sampling units SP1 and SP2 in the previous stage.
Y ij = <X 1i · X 2j * > / ABS [<X 1i · X 2j * >] (4)
Further, the phase difference detection unit PD selects a phase difference Y pq between any antenna element pair among the phase differences of the antenna element pair, and Y ij and Y pq given by the following equation (5) The phase difference Z ij is calculated and the calculation result is passed to the azimuth calculation unit DD in the subsequent stage.
Z ij = Y ij · Y pq * (5)
[0022]
Further, the phase difference detector PD calculates the phase difference B ij (θ) between the mode vectors calculated according to the above equation (2) and the phase difference B pq (θ) corresponding to the arbitrary antenna element pair Y pq. The phase difference C ij (θ) is calculated, and the calculation result is passed to the azimuth calculation unit DD in the subsequent stage.
C ij (θ) = B ij (θ) • B pq * (θ) (6)
[0023]
The azimuth calculating unit DD replaces the phase difference Y ij given by the above-described formula (1) and the phase difference B ij (θ) given by the formula (2) with Z ij given by the formula (5) above. The evaluation function PP IF (θ) is calculated from C ij (θ) given by equation (6).
PP IF (θ) = Σ i Σ j ABS [Z ij −C ij (θ)] (7)
Then, the azimuth calculation unit DD detects the value of θ that minimizes the evaluation function PP IF (θ) of the equation (7) as the arrival direction of the radio wave.
[0024]
In the above equation (7), the amount of phase delay included in the receiving system is canceled out and removed, so that the result is the same as when calibration was automatically performed. That is, the phase delay amount in the reception system 1 composed of the receiver R1, the sampling unit SP1 and the signal line is α, and the phase delay amount in the reception system 2 composed of the receiver R2, the sampling unit SP2 and the signal line is β. To do. In this case, the actually measured complex reception voltages X 1i and X 2j are given by the following equations.
X 1i = V 1i exp (j (φ i + α)) (8)
X 2j = V 2j exp (j (φ j + β)) (9)
[0025]
Substituting Equation (8) and Equation (9) into Equation (4),
Figure 0003682398
[0026]
Substituting (10) and (11) into (5),
Figure 0003682398
[0027]
As is clear from the equation (12), the phase delay amounts α and β of the receiving systems 1 and 2 are canceled out and removed, so that they are not included in Z ij . As described above, since the phase delay amounts α and β of the receiving systems 1 and 2 are canceled out and removed, it is necessary to previously measure α and β in the conventional radio direction measuring device. No calibration hardware is required.
[0028]
In the equation (7), the antenna element pair A p , A q is fixed to an arbitrary pair. However, without fixing the antenna element pair, the evaluation function is calculated while selecting a plurality of antenna element pairs A p and A q from the antenna element pairs used for measurement, and the average of these evaluation functions is finally obtained. By calculating as a simple evaluation function, it is possible to further reduce the influence of noise and increase the measurement accuracy.
[0029]
That is, the evaluation function PP IF (θ) of the equation (7) is calculated for each antenna element pair A p , A q selected from the antenna element pair used for measurement, and the sum of these is calculated as a new evaluation function PPP IF ( θ) = p Σ qi Σ j ABS [Z ij −C ij (θ)] (13)
Calculated as Then, θ that minimizes the new evaluation function PPP IF (θ) is detected as the arrival direction. In this new evaluation function, the influence of noise generated in a specific antenna element is reduced by averaging between the antenna elements, and the direction detection accuracy is further improved.
[0030]
In the following, the results of confirmation by computer simulation performed on the embodiment shown in FIG. 1 will be described with reference to FIGS.
[0031]
FIG. 4 is a plan view showing the arrangement of each antenna element assumed in this computer simulation. Five antenna elements 1 to 5 are arranged at equal intervals on the assumed circumference in a horizontal plane. The frequency of the incoming radio wave is set so that the corresponding half wavelength (λ / 2) is equal to the linear distance between the adjacent antenna elements. In addition, the antenna element pairs are five pairs that are combined into one jump of the elements 1, 3, 2, 4, 3, 3, 5, 4, 1, 5, and 2. Further, a phase difference of 250 ° is set between the receiving systems 1 and 2. Further, the received radio wave comes from an azimuth 150 °, the signal S / N is 10 dB, and the average number of data is 20 times.
[0032]
FIG. 5 is a graph in which the calculation result of LOG 10 (1 / P IF (θ)) is plotted for P IF (θ) in equation (3). Figure 6 is a graph plotting the calculation results for (7) of the PP IF (θ), LOG 10 (1 / PP IF (θ)). FIG. 7 is a graph in which the calculation result of LOG 10 (1 / PPP IF (θ)) is plotted for PPP IF (θ) in equation (13).
[0033]
The graph of FIG. 5 shows the result when ideal calibration is performed using a considerable amount of hardware as shown in FIGS. On the other hand, the graphs of FIG. 6 and FIG. 7 showing the results of the simulation performed on the radio direction measuring apparatus of this example are the results in the case where no calibration is performed, but these are the same as in the case of FIG. It is the result.
[0034]
As described above, since calibration before the start of actual measurement is unnecessary, calibration hardware is not required. At the same time, as in the conventional device, in order to reduce the amount of calibration data stored in the memory, to shorten the time required for calibration, the frequency characteristics of the receivers R1 to Rn are made uniform, and the input / output characteristics The necessity of enlarging the straight line region and stabilizing the temperature characteristic is reduced. As a result, no special receiver is required as the receivers R1 to Rn, and simple and inexpensive ones having normal frequency characteristics and temperature stability can be used.
[0035]
Heretofore, for convenience of explanation, the case where the phase difference detection unit and the direction calculation unit are individually configured has been illustrated. However, these units can be configured as an integrated processing unit realized by a microcomputer or the like.
[0036]
【The invention's effect】
As described above in detail, the radio direction measuring device of the present invention, a phase difference Y ij reception voltage obtained from the two reception systems, instead of the phase difference B ij mode vector (θ), Y ij Y Since the evaluation function P (θ) is calculated using pq * and B ij (θ) B pq (θ) * , the phase difference included in each receiving system is canceled out in the evaluation function. Removed. As a result, no calibration signal generator, high-performance and expensive cable, or hardware such as a switching unit is required to calibrate the phase difference included in each receiving system, greatly reducing the manufacturing cost of the entire device. Is done.
[0037]
Furthermore, in the radio direction measuring device of the present invention, a simple and inexpensive device having normal frequency characteristics, input / output linearity and temperature stability can be used as a receiver. There is also an effect that the cost is reduced.
[Brief description of the drawings]
FIG. 1 is a functional block diagram showing a configuration of a radio direction measuring apparatus according to an embodiment of the present invention.
FIG. 2 is a functional block diagram showing an example of the configuration of a conventional radio direction measuring apparatus.
FIG. 3 is a functional block diagram showing another example of the configuration of a conventional radio direction measuring apparatus.
FIG. 4 is a plan view showing an arrangement of antenna elements used in an electronic computer simulation for confirming the effect of the present invention.
FIG. 5 is a conventional evaluation function obtained by the electronic computer simulation.
FIG. 6 is an evaluation function according to an embodiment of the present invention obtained by the computer simulation.
FIG. 7 is an evaluation function according to another embodiment of the present invention obtained by the computer simulation.
[Explanation of symbols]
A1 ~ An antenna element
R1 ~ Rn receiver
SP1 to SPn sampling section
PS antenna element pair selector
PD phase difference detector
DD bearing calculator
CS calibration signal generator
S1 ~ Sn switching part

Claims (2)

複数のアンテナ素子から成るアレーアンテナと、
前記各アンテナ素子から出力される受信信号が選択的に入力される2個の受信機と、
前記各アンテナ素子の中から1対(Ai ,Aj )を順次選択して、それぞれの受信信号を前記2個の受信機に選択的に入力させるアンテナ素子対選択部と、
前記2個の受信機の出力から、前記選択されたアンテナ素子対(Ai ,Aj )の受信信号の位相差Yijを算定すると共に、前記アレーアンテナのモードベクトルから得られる前記選択されたアンテナ素子対(Ai ,Aj )の出力の位相差Bij(θ)との差に基づく評価関数P IF (θ)=Σ i Σ j ABS[ ij −B ij (θ) ] (ただし、 ABS[ z ] はzの絶対値、記号Σ i Σ j zは (i,j) の一部又は全ての組合せについてのzの総和を意味する)を算定し、この評価関数 P IF (θ)から入射角θを測定方位として得る演算部とを備えたインタ−フェロメ−タ方式の無線方位測定装置に変更を加えることにより受信系の較正を不要にした無線方位測定装置であって、
前記変更は、前記演算部が、
前記順次選択され、算定されたアンテナ素子対の受信信号の位相差 ij のうち任意の一つYpqを選択し、前記位相差Yijの代わりにこの位相差Yijと前記位相差Ypqとの差 ij (=Y ij ・Y pq * を用いると共に、前記位相差Bij(θ)の代わりにこの位相差Bij(θ)と前記位相差Ypqに対応するBpq(θ)との差 ij (θ)(=B ij (θ)・B pq * (θ))を用いて前記評価関数P IF (θ)の代わりに評価関数PP IF (θ)=Σ i Σ j ABS[ ij −C ij (θ) ] から入射角θを算定する手段を備えたことであることを特徴とする無線方位測定装置
An array antenna comprising a plurality of antenna elements;
Two receivers to which a reception signal output from each antenna element is selectively input;
An antenna element pair selector for sequentially selecting a pair (A i , A j ) from the antenna elements and selectively inputting the received signals to the two receivers;
The phase difference Y ij of the received signal of the selected antenna element pair (A i , A j ) is calculated from the outputs of the two receivers, and the selected from the mode vector of the array antenna Evaluation function P IF (θ) = Σ i Σ j ABS [ Y ij −B ij (θ) ] (wherein the difference from the output phase difference B ij (θ) of the antenna element pair (A i , A j )) , ABS [z] is the absolute value of z, and the symbol Σ i Σ j z is the sum of z for some or all combinations of (i, j)) , and this evaluation function P IF (θ A wireless azimuth measuring apparatus that eliminates the need for calibration of the receiving system by adding a change to the interferometer type wireless azimuth measuring apparatus provided with a calculation unit that obtains the incident angle θ as a measurement azimuth from
The change is performed by the calculation unit ,
Wherein are sequentially selected, calculated by selecting an arbitrary one Y pq of the phase difference Y ij antenna element pairs of the received signal, the phase difference Y the phase difference between the phase difference Y ij instead of ij Y pq the difference with using Z ij (= Y ij · Y pq *), the phase difference B ij (θ) and corresponding to the phase difference Y pq B pq instead of the phase difference B ij (θ) (θ with ) Using the difference C ij (θ) (= B ij (θ) · B pq * (θ) ) , instead of the evaluation function P IF (θ) , the evaluation function PP IF (θ) = Σ i Σ j A wireless azimuth measuring apparatus comprising means for calculating an incident angle θ from ABS [ Z ij −C ij (θ) ] .
請求項1において、
前記演算部が、
前記評価関数の算定に際し、前記位相差Ypqを複数個選択して、それぞれの場合の評価関数を算定し、これらを平均化したものを最終的な評価関数とする手段を更に備えたことを特徴とする無線方位測定装置。
In claim 1,
The computing unit is
The calculation of the evaluation function further includes means for selecting a plurality of the phase differences Y pq , calculating the evaluation function in each case, and averaging these to obtain a final evaluation function A wireless azimuth measuring device.
JP2000087082A 2000-03-23 2000-03-23 Wireless direction measuring device Expired - Lifetime JP3682398B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2000087082A JP3682398B2 (en) 2000-03-23 2000-03-23 Wireless direction measuring device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2000087082A JP3682398B2 (en) 2000-03-23 2000-03-23 Wireless direction measuring device

Publications (2)

Publication Number Publication Date
JP2001272447A JP2001272447A (en) 2001-10-05
JP3682398B2 true JP3682398B2 (en) 2005-08-10

Family

ID=18603151

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2000087082A Expired - Lifetime JP3682398B2 (en) 2000-03-23 2000-03-23 Wireless direction measuring device

Country Status (1)

Country Link
JP (1) JP3682398B2 (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4003229B2 (en) 2002-06-20 2007-11-07 日本電気株式会社 Array antenna receiver
JP2006121525A (en) * 2004-10-22 2006-05-11 Koden Electronics Co Ltd All direction spatial diversity receiving apparatus
JP4907910B2 (en) * 2005-07-04 2012-04-04 株式会社東芝 Radio wave emission source detection device
JP4980094B2 (en) * 2007-02-16 2012-07-18 三菱電機株式会社 Displacement measuring device
JP5047661B2 (en) * 2007-03-27 2012-10-10 三菱電機株式会社 Angle measuring device and angle measuring method
JP2012078163A (en) * 2010-09-30 2012-04-19 Toshiba Corp Composite guidance system and composite guidance method
JP5955140B2 (en) * 2012-07-10 2016-07-20 三菱電機株式会社 Angle measuring device and angle measuring method
JP2015087132A (en) * 2013-10-28 2015-05-07 株式会社東芝 Signal detection device and signal detection method

Also Published As

Publication number Publication date
JP2001272447A (en) 2001-10-05

Similar Documents

Publication Publication Date Title
CN111044971B (en) Two-dimensional interferometer direction finding system
CN108776330B (en) High-precision calibration method and device for multiple receiving channels of FMCW radar
JP7356510B2 (en) Formaldehyde concentration detection method, device and air purifier
US10732253B2 (en) Antenna arrangements for measurement of angle of arrival
US9551777B2 (en) Direction finding using antenna array rotation
US6169519B1 (en) TCAS bearing measurement receiver apparatus with phase error compensation method
CN105738851B (en) A kind of synthetic aperture microwave radiometer joint bearing calibration
CN111903218B (en) Interferometer direction finding method based on Chinese remainder theorem
JP3682398B2 (en) Wireless direction measuring device
JP2006284591A (en) System and method for improving precision of electric equipment by correcting environmental condition
US10545182B2 (en) Crosstalk calibration for multi-channel systems
CN113258288A (en) Phased array antenna beam control device and control method
KR101640787B1 (en) The system for finding direction of a wireless signal
US5122808A (en) Phase only bearing mesurement with amiguity correction in a collision avoidance system
Will et al. Segmental polynomial approximation based phase error correction for precise near field displacement measurements using Six-Port microwave interferometers
JP3808431B2 (en) Direction finding device
CN110896334B (en) High-precision channel calibration method based on channel phase noise statistics
CN113567916B (en) Correlation interferometer direction finding method, system, equipment and storage medium
CN106918826B (en) A kind of calibration range-measurement system of interferometric radiometer, calibration and H_2O maser method
CN102508060B (en) External pulse interference resistant array calibration pretreatment method
US6542832B1 (en) Error detection and correction system for use with dual-pulse output metering devices
JP2016197053A (en) Bearing detection device and bearing detection method
KR101022358B1 (en) Least square digital frequency discriminator in phase calibration and method thereof
CN107144343A (en) Low-frequency vibration displacement transducer network-building method, system and device
US3303503A (en) Single-site monopulse ranging system

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20041202

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050208

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050325

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: 20050517

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050523

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 3682398

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090527

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090527

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100527

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110527

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120527

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130527

Year of fee payment: 8

EXPY Cancellation because of completion of term