JPWO2012063532A1 - Arrival angle calculation device - Google Patents
Arrival angle calculation device Download PDFInfo
- Publication number
- JPWO2012063532A1 JPWO2012063532A1 JP2012542831A JP2012542831A JPWO2012063532A1 JP WO2012063532 A1 JPWO2012063532 A1 JP WO2012063532A1 JP 2012542831 A JP2012542831 A JP 2012542831A JP 2012542831 A JP2012542831 A JP 2012542831A JP WO2012063532 A1 JPWO2012063532 A1 JP WO2012063532A1
- Authority
- JP
- Japan
- Prior art keywords
- unit
- output
- arrival angle
- calculation
- component
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/46—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
- G01S3/48—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems the waves arriving at the antennas being continuous or intermittent and the phase difference of signals derived therefrom being measured
Abstract
高い精度で到来角度を算出できる到来角度算出装置を提供することを目的とする。本発明の到来角度算出装置(1)は、複数のアンテナと、複数の受信信号処理部と、到来角度算出部と、を備え、各受信信号処理部は、受信部(12a、12b)と、相関処理部(21a、21b)と、ピーク検出部(22a、22a)と、タイミング制御部(23a、23b)と、を備え、タイミング制御部(23a、23b)は、情報単位に相当する期間におけるピーク期間の電力と当該ピーク期間を除いた期間の電力との比が、しきい値より大きければ、相関処理部(21a、21b)からの信号を到来角度算出部(24)に出力することを特徴とする。An object of the present invention is to provide an arrival angle calculation device that can calculate an arrival angle with high accuracy. The arrival angle calculation device (1) of the present invention includes a plurality of antennas, a plurality of reception signal processing units, and an arrival angle calculation unit, and each reception signal processing unit includes a reception unit (12a, 12b), A correlation processing unit (21a, 21b), a peak detection unit (22a, 22a), and a timing control unit (23a, 23b); and the timing control unit (23a, 23b) in a period corresponding to an information unit If the ratio between the power in the peak period and the power in the period excluding the peak period is larger than the threshold, the signal from the correlation processing unit (21a, 21b) is output to the arrival angle calculating unit (24). Features.
Description
本発明は、到来する電波の位相を検出して電波到来角度の算出に用いる到来角度算出装置に関する。 The present invention relates to an arrival angle calculation device that detects the phase of an incoming radio wave and uses it to calculate the radio wave arrival angle.
従来の到来方向推定装置においては、相互相関係数の算出や逆行列演算等の演算量の大きい演算が用いられており、数百シンボル分もの演算が必要であった。このため、簡便な演算で到来方向を推定できる到来方向推定装置が望まれていた。 In the conventional direction-of-arrival estimation apparatus, calculation with a large calculation amount such as calculation of cross-correlation coefficient and inverse matrix calculation is used, and calculation for several hundred symbols is required. For this reason, an arrival direction estimation device capable of estimating the arrival direction by simple calculation has been desired.
特許文献1において、演算規模を縮小した到来方向推定装置が提案されている。特許文献1に記載の到来方向推定装置では、二つのアンテナで受信した受信信号に対して、複素共役回路と乗算回路によって到来方向の係数を算出し、到来方向検出回路において逆正接演算と逆余弦演算を行うことにより、受信波の到来方向を推定している。
しかしながら、特許文献1では、1スロット間の到来方向ベクトルのパワーとしきい値とを比較して、しきい値より大きい場合に到来方向を更新する構成を採用しているため、受信波中の希望波を正確に検出して到来方向を更新することができない可能性がある。例えば、受信波のバックグラウンドレベルが高い場合には、希望波の信号レベルに関わらず、受信波のパワーがしきい値より大きくなることがある。このような状況では、バックグラウンドレベルにおいても到来方向が算出、更新されるため、到来方向を正しく推定することができない。
However,
本発明はかかる点に鑑みてなされたものであり、受信波のバックグラウンドレベルの影響を抑えることができ、高い精度で到来角度を算出できる到来角度算出装置を提供することを目的とする。 The present invention has been made in view of this point, and an object of the present invention is to provide an arrival angle calculation apparatus that can suppress the influence of the background level of the received wave and can calculate the arrival angle with high accuracy.
本発明の到来角度算出装置は、ある位置から送信された電波を受信する複数のアンテナと、前記各アンテナに対応して設けられた複数の受信信号処理部と、前記複数の受信信号処理部から出力される出力信号から受信信号処理部間で同一情報単位となる信号成分を取り込んで前記電波の到来角度を算出する到来角度算出部と、を備え、前記各受信信号処理部は、対応する前記アンテナで受信した電波を当該電波の位相情報を有する受信信号に変換して出力する受信部と、前記受信部から出力された受信信号を相関処理する相関処理部と、前記相関処理された受信信号のピークを検出するピーク検出部と、前記相関処理部の出力信号から前記受信信号処理部間で同一情報単位となる信号成分が切り出されるように、前記ピーク検出部で検出されたピークのタイミングに合わせて、前記相関処理部から出力される出力信号の取り込みタイミングを制御するタイミング制御部と、を備え、前記タイミング制御部は、前記情報単位に相当する期間におけるピーク期間の電力と当該ピーク期間を除いた期間の電力との比が、しきい値より大きければ、前記相関処理部からの信号を前記到来角度算出部に出力することを特徴とする。 An arrival angle calculation apparatus according to the present invention includes a plurality of antennas that receive radio waves transmitted from a certain position, a plurality of reception signal processing units provided corresponding to the antennas, and a plurality of reception signal processing units. An arrival angle calculation unit that calculates an arrival angle of the radio wave by taking in a signal component that is the same information unit between the reception signal processing units from the output signal that is output, and each reception signal processing unit corresponds to the A reception unit that converts a radio wave received by an antenna into a reception signal having phase information of the radio wave and outputs the signal, a correlation processing unit that performs correlation processing on the reception signal output from the reception unit, and the correlation-processed reception signal Detected by the peak detection unit so that a signal component that is the same unit of information is extracted from the output signal of the correlation processing unit and the received signal processing unit from the output signal of the correlation processing unit. A timing control unit that controls the timing of capturing the output signal output from the correlation processing unit in accordance with the peak timing, and the timing control unit includes power in a peak period in a period corresponding to the information unit, If the ratio with the power in the period excluding the peak period is larger than a threshold value, the signal from the correlation processing unit is output to the arrival angle calculation unit.
この構成によれば、ピーク期間の電力とそれ以外の期間の電力の比をしきい値と比較し、比がしきい値より大きい場合に到来角度を算出するため、受信波の希望波以外の信号レベルが高い場合であっても、希望波のピークを正確に検出し、到来角度を算出することができる。つまり、希望波以外の部分から到来角度を算出することが無いため、到来角度の算出精度を高めることができる。 According to this configuration, the ratio of the power in the peak period and the power in the other periods is compared with a threshold value, and when the ratio is larger than the threshold value, the arrival angle is calculated. Even if the signal level is high, it is possible to accurately detect the peak of the desired wave and calculate the arrival angle. That is, since the arrival angle is not calculated from a portion other than the desired wave, the calculation accuracy of the arrival angle can be improved.
本発明の到来角度算出装置において、前記タイミング制御部は、前記情報単位に相当する期間におけるピーク期間の電力の和ΣP1と前記情報単位に相当する期間において前記ピーク期間を除いた期間における電力の和ΣP2との比ΣP1/ΣP2を、しきい値と比較し、前記比ΣP1/ΣP2が前記しきい値より大きい場合、前記相関処理部からの信号を前記到来角度算出部に出力しても良い。In the arrival angle calculating apparatus of the present invention, the timing control unit, a period in the power of excluding the peak period in the period corresponding to the information unit to the sum .SIGMA.P 1 power peak periods in the period corresponding to the information unit The ratio ΣP 1 / ΣP 2 with respect to the sum ΣP 2 is compared with a threshold value. When the ratio ΣP 1 / ΣP 2 is larger than the threshold value, the signal from the correlation processing unit is sent to the arrival angle calculation unit. It may be output.
本発明の到来角度算出装置において、前記到来角度算出部は、一方のアンテナに対応した一方の受信信号処理部のタイミング制御部からの信号の複素共役をとる複素共役部と、前記複素共役部の出力と、他方のアンテナに対応した他方の受信信号処理部のタイミング制御部からの信号とを乗算する複素乗算部と、前記複素乗算部の出力を用いて逆正接演算を行い、前記アンテナ間での前記受信電波の位相差を算出する逆正接部と、前記逆正接部の出力を平均化する平均化部と、前記平均化部の出力を用いて逆三角関数演算を行い、到来角度に変換する到来角度変換部と、を備えても良い。この構成によれば、相互相関係数の算出や逆行列演算などを用いることなく到来角度を算出できるため、到来角度算出装置の規模を小さくすることができる。 In the arrival angle calculation device according to the present invention, the arrival angle calculation unit includes a complex conjugate unit that takes a complex conjugate of a signal from a timing control unit of one received signal processing unit corresponding to one antenna, and the complex conjugate unit. A complex multiplier that multiplies the output and the signal from the timing control unit of the other received signal processing unit corresponding to the other antenna, and performs an arctangent calculation using the output of the complex multiplier, and An arctangent that calculates the phase difference of the received radio wave, an averaging unit that averages the output of the arctangent, and an inverse trigonometric function calculation using the output of the averaging unit to convert to an arrival angle And an arrival angle conversion unit. According to this configuration, since the arrival angle can be calculated without using cross-correlation coefficient calculation or inverse matrix calculation, the size of the arrival angle calculation device can be reduced.
本発明の到来角度算出装置において、前記到来角度算出部は、前記算出された位相差が、I−Q平面上における+180°及び/又は−180°付近に分布している場合は、各位相差を所定角度回転させた上で平均化し、当該平均値から前記所定角度を減じて逆三角関数演算を行い、到来角度を算出しても良い。この構成によれば、位相差が、到来角度の算出精度が低下する傾向にある位相差領域に分布する場合に、位相差を所定角度だけ回転させて到来角度算出の演算を行うため、到来角度の算出精度が低下せずに済む。その結果、到来角度の算出精度を十分に高めることができる。 In the arrival angle calculation device of the present invention, the arrival angle calculation unit may calculate each phase difference when the calculated phase difference is distributed around + 180 ° and / or −180 ° on the IQ plane. The arrival angle may be calculated by averaging after rotating by a predetermined angle, subtracting the predetermined angle from the average value, and performing an inverse trigonometric function calculation. According to this configuration, when the phase difference is distributed in the phase difference region where the calculation accuracy of the arrival angle tends to decrease, the arrival angle calculation is performed by rotating the phase difference by a predetermined angle. The calculation accuracy does not decrease. As a result, the calculation accuracy of the arrival angle can be sufficiently increased.
本発明の到来角度算出装置において、前記I−Q平面上において、+90°より大きく、または−90°より小さい位相差の数が、+90°より小さくかつ−90°より大きい位相差の数より多い場合、前記I−Q平面上における+180°及び/又は−180°付近に分布していると判断しても良い。 In the arrival angle calculation apparatus of the present invention, the number of phase differences larger than + 90 ° or smaller than −90 ° on the IQ plane is larger than the number of phase differences smaller than + 90 ° and larger than −90 °. In this case, it may be determined that the distribution is in the vicinity of + 180 ° and / or −180 ° on the IQ plane.
本発明の到来角度算出装置において、前記所定角度は、+90°、−90°、+180°、または−180°のいずれかの角度としても良い。 In the arrival angle calculation device of the present invention, the predetermined angle may be any of + 90 °, −90 °, + 180 °, or −180 °.
本発明の到来角度算出装置において、前記複素乗算部の出力のI成分が負であり、前記複素乗算部の出力のI成分の絶対値がQ成分の絶対値と比べて十分に大きい場合、前記Q成分の符号を反転させた上でI成分とQ成分との関係を入れ替えた逆正接演算を行うことにより補正された位相差を算出し、前記補正された位相差を平均化し、当該平均値から90°を減じて逆三角関数演算を行い、到来角度に変換しても良い。この構成によれば、位相差が、到来角度の算出精度が低下する傾向にある位相差領域に分布する場合に、位相差を所定角度だけ回転させて到来角度算出の演算を行うため、到来角度の算出精度が低下せずに済む。その結果、到来角度の算出精度を十分に高めることができる。 In the arrival angle calculation apparatus of the present invention, when the I component of the output of the complex multiplier is negative and the absolute value of the I component of the output of the complex multiplier is sufficiently larger than the absolute value of the Q component, A phase difference corrected by performing an arctangent calculation in which the relationship between the I component and the Q component is reversed after inverting the sign of the Q component, averaging the corrected phase difference, and calculating the average value Alternatively, 90 ° may be subtracted from the above and inverse trigonometric function calculation may be performed to convert the angle of arrival. According to this configuration, when the phase difference is distributed in the phase difference region where the calculation accuracy of the arrival angle tends to decrease, the arrival angle calculation is performed by rotating the phase difference by a predetermined angle. The calculation accuracy does not decrease. As a result, the calculation accuracy of the arrival angle can be sufficiently increased.
本発明の到来角度算出装置において、前記複素乗算部の出力のI成分が負であり、前記複素乗算部の出力のI成分の絶対値がQ成分の絶対値と比べて十分に大きい場合、前記I成分の符号を反転させた上でI成分とQ成分との関係を入れ替えた逆正接演算を行うことにより補正された位相差を算出し、前記補正された位相差を平均化し、当該平均値に90°を加えて逆三角関数演算を行い、到来角度に変換しても良い。この構成によれば、位相差が、到来角度の算出精度が低下する傾向にある位相差領域に分布する場合に、位相差を所定角度だけ回転させて到来角度算出の演算を行うため、到来角度の算出精度が低下せずに済む。その結果、到来角度の算出精度を十分に高めることができる。
In the arrival angle calculation apparatus of the present invention, when the I component of the output of the complex multiplier is negative and the absolute value of the I component of the output of the complex multiplier is sufficiently larger than the absolute value of the Q component, A phase difference corrected by performing an arc tangent calculation in which the relationship between the I component and the Q component is exchanged after inverting the sign of the I component, averaging the corrected phase difference, and calculating the
本発明の到来角度算出装置において、前記複素乗算部の出力のI成分が負であり、前記複素乗算部の出力のI成分の絶対値がQ成分の絶対値と比べて十分に大きい場合、前記I成分の符号とQ成分の符号を反転させた上で逆正接演算を行うことにより補正された位相差を算出し、前記補正された位相差を平均化し、当該平均値から180°を減じて逆三角関数演算を行い、到来角度に変換しても良い。この構成によれば、位相差が、到来角度の算出精度が低下する傾向にある位相差領域に分布する場合に、位相差を所定角度だけ回転させて到来角度算出の演算を行うため、到来角度の算出精度が低下せずに済む。その結果、到来角度の算出精度を十分に高めることができる。 In the arrival angle calculation apparatus of the present invention, when the I component of the output of the complex multiplier is negative and the absolute value of the I component of the output of the complex multiplier is sufficiently larger than the absolute value of the Q component, Calculate the phase difference corrected by performing the arctangent calculation after inverting the sign of the I component and the sign of the Q component, averaging the corrected phase difference, and subtracting 180 ° from the average value An inverse trigonometric function calculation may be performed to convert to an arrival angle. According to this configuration, when the phase difference is distributed in the phase difference region where the calculation accuracy of the arrival angle tends to decrease, the arrival angle calculation is performed by rotating the phase difference by a predetermined angle. The calculation accuracy does not decrease. As a result, the calculation accuracy of the arrival angle can be sufficiently increased.
本発明の到来角度算出装置によれば、ピーク期間の電力とピーク期間以外の残りの期間における電力との比を求めて、この求めた比としきい値とを比較し、比がしきい値より大きい場合に到来角度を算出するため、受信波の希望波以外の信号レベルが高い場合であっても、希望波のピークを正確に検出し、到来角度を算出することができる。つまり、希望波以外の部分から到来角度を算出することが無いため、到来角度の算出精度を高めることができる。 According to the arrival angle calculation device of the present invention, the ratio between the power in the peak period and the power in the remaining period other than the peak period is obtained, and the obtained ratio is compared with the threshold value. Since the arrival angle is calculated when it is large, the peak of the desired wave can be accurately detected and the arrival angle can be calculated even when the signal level of the received wave other than the desired wave is high. That is, since the arrival angle is not calculated from a portion other than the desired wave, the calculation accuracy of the arrival angle can be improved.
図1は、本発明の一実施の形態に係る到来角度算出装置の構成例を示すブロック図である。本実施の形態に係る到来角度算出装置1は、所定の発振周波数で基準信号を発振可能な基準信号発生部10と、所定間隔離して配置された受信用アンテナ11a、11bと、受信用アンテナ11a、11bで受けた電波を、基準信号発生部10から出力される基準信号を用いて受信信号に変換し出力する受信部12a、12bと、受信部12a、12bから出力される受信信号から到来角度算出のための各種演算処理を行う演算部13と、を備える。なお、到来角度算出装置1は、電波の伝搬遅延に起因する位相遅れに基づいて到来角度を算出するため、同じ情報を持つ電波を所定間隔離れた二点(または二以上の点)で受信する必要がある。このため、受信電波に対応する二つ(またはそれ以上)のアンテナおよび受信系を備えていることが必要である。なお、同一の到来電波(同じ情報単位)を所定間隔離れた二以上の位置で受信できるのであれば、到来角度算出装置1は、二以上の受信系を備えている構成に限定されない。
FIG. 1 is a block diagram showing a configuration example of an arrival angle calculation apparatus according to an embodiment of the present invention. The arrival
受信部12a、12bは、ローノイズアンプ、ミキサ、バンドパスフィルタなどを含んで構成され、所定周波数の電波を受信できるように構成されている。演算部13は、受信信号の相関処理を行う相関処理部21a、21bと、相関処理された受信信号のピークを検出するピーク検出部22a、22bと、ピーク検出部22a、22bで検出されたピークのタイミングに合わせて相関処理部21a、21bからの信号を出力するタイミング制御部23a、23bと、タイミング制御部23a、23bからの信号に基づいて、到来角度の計算を行う到来角度算出部24と、を含んで構成される。なお、演算部13の構成や機能は、ハードウェアで実現しても良いし、ソフトウェアで実現しても良い。
The receiving
相関処理部21a、21bは、受信部12a、12bからの受信信号と当該受信信号と相関の高い信号とを乗算して出力する。相関処理部21a、21bにおいて乗じられる信号は受信信号との相関が高いため、相関処理部21a、21bから出力される信号は、相関区間でピークとなる。ピーク検出部22a、22bは、相関処理部21a、21bからの出力信号の電力を算出し、出力信号の電力ピークを検出する。タイミング制御部23a、23bは、ピーク検出部22a、22bにおいて検出されたピークタイミングに合わせて、相関処理部21a、21bからの出力信号を到来角度算出部24に出力する。具体的には、検出されたピーク期間の電力から算出される情報を元に、相関処理部21a、21bからの出力信号を到来角度算出部24に出力するか否かを決定する。
The
図2は、変調方式として直接スペクトラム拡散(DSSS)を用いる場合の到来角度算出装置の具体的構成例を示すブロック図である。なお、図2では、図1における演算部13に相当する構成のみを示している。
FIG. 2 is a block diagram illustrating a specific configuration example of an arrival angle calculation apparatus when direct spread spectrum (DSSS) is used as a modulation method. 2 shows only the configuration corresponding to the
図2において、相関処理部21aは、拡散コードを発生する拡散コード発生器31と、受信信号と拡散コードとを乗算する乗算器32aおよび32bと、乗算器32aおよび32bの出力を1ビット期間分だけ足し合わせてピーク検出部22aおよびタイミング制御部23aに出力する加算器33aおよび33bとを備える。ピーク検出部22aは、加算器33aおよび33bから出力された信号の電力を算出する電力算出部34aと、その電力ピークを検出してタイミング制御部23aに出力するピーク電力検出部35aとを備える。タイミング制御部23aは、ピーク電力検出部35aからの信号を元に加算器33aおよび33bからの信号の到来角度算出部24への出力可否を制御するバッファ部36aを備える。同様に、相関処理部21bは、拡散コード発生器31、乗算器32cおよび32d、加算器33cおよび33dを備え、ピーク検出部22bは、電力算出部34b、ピーク電力検出部35bを備え、タイミング制御部23bはバッファ部36bを備える。到来角度算出部24は、バッファ部36aの出力の複素共役をとる複素共役部41と、複素共役部41の出力とバッファ部36bの出力とを複素乗算する複素乗算部42と、複素乗算部42の出力を用いて逆正接演算を行う逆正接部43と、複素乗算部42の出力信号からチップ区間ごとの電力を算出する電力算出部44と、電力算出部44からの情報に基づいて逆正接部43の出力を平均化する平均化部45と、平均化部45の出力を用いて到来角度に変換する到来角度変換部46とを備える。
In FIG. 2, a
拡散コード発生器31は、DSSSによって周波数軸上に拡散された信号を逆拡散するための拡散コードを発生する。当該拡散コードは、送信側でコード変調(拡散)の際に使用された拡散コードに対応するものである。乗算器32aおよび32bは、受信信号に上記拡散コードを乗じて逆拡散を行う。乗算器32aには、受信部12aからの受信信号のうちの同相成分I1が入力される。また、乗算器32bには、受信部12aからの受信信号のうちの直交成分Q1が入力される。加算器33aおよび33bは、乗算器32aおよび32bのチップ区間ごとの出力を1ビットに相当する期間(ビット区間)足し合わせて出力する。図3(a)に加算器33aからの出力波形の例を示す。図3(b)は、図3(a)に示す出力波形の部分拡大図である。また、図3(c)に加算器33bからの出力波形の例を示す。図3(d)は、図3(c)に示す出力波形の部分拡大図である。
The spreading
加算器33aの出力信号および加算器33bの出力信号は、ピーク検出部22aの電力算出部34a、およびタイミング制御部23aのバッファ部36aに入力される。電力算出部34aは、加算器33aおよび33bの出力信号からチップ区間ごとの電力を算出する。具体的には、電力算出部34aは、同相成分に相当する加算器33aの出力信号の絶対値と、直交成分に相当する加算器33bの出力信号の絶対値とを足し合わせ、チップ区間ごとの電力情報としてピーク電力検出部35aに出力する。ピーク電力検出部35aは、チップ区間ごとの電力情報を受け取ると、受信信号中の電力ピークを検出し、電力ピーク情報としてタイミング制御部23aのバッファ部36aに出力する。なお、加算器33aの出力信号の2乗値と、加算器33bの出力信号の2乗値とを足し合わせてピーク電力検出部35aに出力しても良い。
The output signal of the
ピーク検出部22a(ピーク電力検出部35a)から出力される電力ピーク情報は、受信信号のピークの有無を判定する情報である。具体的には、電力ピーク情報は、受信信号のピーク点付近の期間(ピーク期間)における電力の和ΣP1と、DSSSでの情報単位となる1ビット期間からピーク期間を除いた期間における電力の和ΣP2との比R(=ΣP1/ΣP2)がしきい値Rthより大きいか否かを示す情報である。電力ピーク情報において、RがRthより大きい場合には、タイミング制御部23a(バッファ部36a)は、そのタイミングで受信信号がピークを有するものとして、1ビット分の信号Ia1および信号Qa1を到来角度算出部24に出力する。一方で、RがRthより小さい場合には、タイミング制御部23a(バッファ部36a)は、そのタイミングでは受信信号がピークを有しないものとして、到来角度算出部24への出力を停止する。なお、ここでは、ピーク検出部22aが電力ピーク情報に関する演算処理を行っているが、電力ピーク情報に関する演算処理はタイミング制御部23aにおいて行っても良い。The power peak information output from the
相関処理部21b(拡散コード発生器31、乗算器32cおよび32d、加算器33cおよび33d)、ピーク検出部22b(電力算出部34b、ピーク電力検出部35b)、タイミング制御部23b(バッファ部36b)の動作や機能は、上記相関処理部21a(拡散コード発生器31、乗算器32aおよび32b、加算器33aおよび33b)、ピーク検出部22a(電力算出部34a、ピーク電力検出部35a)、タイミング制御部23a(バッファ部36a)の動作や機能と同様である。ただし、相関処理部21bに入力される受信信号と、相関処理部21aに入力される受信信号とは、同一電波を所定間隔離れた二点で受信しているため位相が僅かに異なっている。このため、タイミング制御部23bから出力される信号と、タイミング制御部23aから出力される信号とでは、位相が僅かに相違する。タイミング制御部23aの出力Oa1、およびタイミング制御部23bの出力Oa2を、同相成分に相当する信号を実部、直交成分に相当する信号を虚部として複素数で表現すると、下記式(1)、(2)のようになる。なお、φ1およびφ2は、各信号の位相を表す。
タイミング制御部23aの出力Oa1は、到来角度算出部24の複素共役部41に入力される。複素共役部41は、タイミング制御部23aの出力Oa1の複素共役を複素乗算部42に出力する。つまり、複素共役部41からは、信号Ia1と、信号Qa1の符号が反転した信号が出力される。複素共役部41の出力Oa1´を複素数で表現すると、下記式(3)のようになる。
複素乗算部42は、複素共役部41の出力Oa1´と、タイミング制御部23bの出力Oa2とを複素乗算して、乗算結果である信号Ibおよび信号Qbを逆正接部43および電力算出部44に出力する。複素乗算部42の出力Ob、出力Obの同相成分Ibおよび直交成分Qbは下記式(4)〜(6)のように表される。
逆正接部43は、複素乗算部42の出力を用いて逆正接演算を行う。具体的には、複素乗算部42の出力信号Ibを分母とし、出力信号Qbを分子とした値の逆正接演算を行う。図4(a)に逆正接部43からの出力波形の例を示す。逆正接部43の出力Oarctanは位相差φ2−φ1に相当し、下記式(7)で表される。
電力算出部44は、複素乗算部42の出力信号からチップ区間ごとの電力を算出する。具体的には、電力算出部44は、Ibの絶対値とQbの絶対値とを足し合わせ、チップ区間ごとの電力情報として平均化部45に出力する。なお、Ibの2乗値と、Qbの2乗値とを足し合わせて平均化部45に出力しても良い。図4(b)に電力算出部44からの出力波形の例を示す。平均化部45は、チップ区間ごとの電力情報を受け取ると、その情報に基づいて逆正接部43の出力Oarctanを平均化して到来角度変換部46に出力する。なお、電力算出部44および平均化部45は、適宜省略して良い。The
到来角度変換部46は、平均化部45の出力(平均化部45を有しない場合は逆正接部43の出力)を用いて逆三角関数演算により到来角度に変換する。逆三角関数演算としては、例えば、逆正弦演算を適用することができる。当該演算によって求められる値、すなわち、到来角度変換部46の出力が、到来角度θ(rad)に相当する。到来角度変換部46の出力Oarcsinは下記式(8)で表される。なお、下記式において、λ(m)は受信波の波長であり、d(m)は受信用アンテナ間の距離である。
上記処理により到来角度が得られるのは、図5に示すような幾何学的な関係が成立するためである。所定の方向を基準として間隔d(m)離して配置された二つの受信用アンテナ11a、11bに到来する電波のなす角度をθ(rad)とする。受信用アンテナ11bに到来する電波の伝搬距離は、受信用アンテナ11aに到来する電波の伝搬距離と比べてΔ(m)だけ長くなり、位相遅延(位相差φ2−φ1(rad))が生じる。このモデルにおいて生じる伝搬距離の差分Δと位相差φ2−φ1との関係を受信波の波長λ(m)を用いて表すと、下記式(9)のようになる。なお、下記式において、Δ<λである。
また、上記モデルにおける伝搬距離の差分Δ、アンテナ間隔d、到来角度θの幾何学的な関係から、下記式(10)が成り立つ。
つまり、到来角度θは下記式(11)のように表されることになる。なお、式(11)は、到来角度変換部46における処理に相当する。このように、本実施の形態の到来角度算出装置によって到来角度が算出されることが分かる。
次に、到来角度算出装置を用いた位置検出システムの例について説明する。図6に示される位置検出システム101は、到来角度算出装置1aと、到来角度算出装置1aと所定距離D離して配置される他の到来角度算出装置1bと、アクセスポイント2又はユーザ端末3とを含んで構成される。アクセスポイント2およびユーザ端末3は、それぞれ送信系および受信系を備え(図示せず)、双方向の情報伝送(通信)が可能に構成されている。また、アクセスポイント2およびユーザ端末3は、それぞれが備える送信系によって、到来角度算出装置1aおよび到来角度算出装置1bに到来角度算出用の電波を送信できるように構成されている。位置検出の対象は、アクセスポイント2又はユーザ端末3のいずれでも良い。
Next, an example of a position detection system using the arrival angle calculation device will be described. The
到来角度算出装置1aは、アクセスポイント2の送信用アンテナから送信された電波を受信用アンテナ11aaおよび11abで受信して、到来角度算出装置1aを基準とする到来角度を算出する。また、到来角度算出装置1bは、アクセスポイント2の送信用アンテナから送信された電波を受信用アンテナ11baおよび11bbで受信して、到来角度算出装置1bを基準とする到来角度を算出する。到来角度算出装置1aと到来角度算出装置1bの位置関係が既知であれば、それぞれを基準とする到来角度からアクセスポイント2の位置を決定することができる。
The arrival angle calculation device 1a receives the radio waves transmitted from the transmission antenna of the
また、ユーザ端末3の位置検出の場合は、到来角度算出装置1aおよび到来角度算出装置1bは、ユーザ端末3から送信される電波の到来角度を算出する。
In the case of detecting the position of the user terminal 3, the arrival angle calculation device 1 a and the arrival
図7は本実施の形態に係る到来角度算出装置1における到来角度算出のフロー図である。到来角度算出装置1が到来角度算出対象の電波を受信すると、受信部12a、12bは相関処理部21a、21bに受信信号を出力する。そして、相関処理部21a、21bは、ステップ201において受信信号の相関処理および加算処理を行う。
FIG. 7 is a flowchart of arrival angle calculation in the arrival
その後、ピーク検出部22a、22bは、ステップ202において相関処理部21a、21bの出力信号から電力のピーク値Ppeakを検出する。そして、ピーク点付近の期間(ピーク期間)における電力の和ΣP1と、1ビット期間(情報単位の期間)からピーク期間を除いた期間における電力の和ΣP2とを算出し、これらの比R(=ΣP1/ΣP2)を算出する。図8(a)には、ピーク検出部22a、22bに入力される信号を模式的に示す。ピーク電力Ppeakは、図8(a)におけるピーク点Pにおける電力であり、ΣP1は、ピーク期間t1における電力の和であり、ΣP2は、1ビット期間からピーク期間t1を除いた期間t2における電力の和である。ここで、ピーク期間t1は、ピークの裾野を含む期間である。例えば、図9に示すように、変調方式としてDSSSを用いる場合には拡散コードの周期tcの2倍の裾野が形成される。このため、当該2・tcの期間をピーク期間t1とすることができる。なお、図9において、期間t2は1ビット期間tbを用いてtb−2・tcと表される。Thereafter, in
図10に、AD変換器を用いて受信信号を取り込んだ場合のピーク検出部22a、22bに入力される信号の例を示す。図10の横軸tはサンプル番号を示しており、tは離散的な値をとる。変調方式としてDSSSを用いる場合、例えば拡散コードが11チップであり、1ビット期間が1μsであれば、拡散コードの1チップ期間は0.091μsである。AD変換は、1チップ期間の4倍のオーバーサンプリングであるとすれば、裾野は1チップ分広がり、ips=ip−3、ipe=ip+3となる。この場合、Rは下記式(12)のように表される。
ステップ203において、ピーク検出部22a、22bは、算出された比R(=ΣP1/ΣP2)と所定のしきい値Rthとを比較する。RthよりRが大きい場合、ピーク検出部22a、22bは、その旨の信号をタイミング制御部23a、23bに出力する。タイミング制御部23aは、RthよりRが大きい旨の信号を受けると、受信信号にピークが存在するものとして到来角度の計算に必要な信号を到来角度算出部24に出力する。そして、ステップ204において、到来角度算出部24は到来角度を算出する。一方、RがRth以下の場合、ピーク検出部22a、22bは、その旨の信号をタイミング制御部23a、23bに出力し、タイミング制御部23aは、受信信号にピークが存在しないものとして到来角度算出部24への出力を停止する。そして、到来角度算出装置1は、ステップ201からのフローを再度実行する。Rthは任意の値である。例えば、Rとの比較でピークの有無を判定できる程度の値をRthとして設定することができる。In
このように、検出されたピークに関する指標(R)としきい値(Rth)とを比較してピークの有無を判定することにより、ピークの有無を正確に判定することが可能である。As described above, the presence or absence of the peak can be accurately determined by comparing the index (R) with respect to the detected peak and the threshold (R th ) to determine the presence or absence of the peak.
ここで、単純に、パワー(電力)とパワーのしきい値とを比較してピークの有無を判定する方法について考察する。図8(b)に、バックグラウンドレベルが高い信号(実線)と、バックグラウンドレベルが低い信号(点線)を模式的に示す。図8(b)に点線で示すようにバックグラウンドレベルが低い場合、パワーのピーク値とパワーのしきい値Pthとを比較することでピークを検出することができる。しかしながら、図8(b)に実線で示すようにバックグラウンドレベルがPthを超える程度に高くなると、パワーのピーク値とPthとを比較してもピークを検出することができない。パワーのピーク値とパワーのしきい値との単純な比較では、バックグラウンドレベルを考慮できないためである。そこで、本実施の形態において示すように、ピークの検出にバックグラウンドレベルを考慮した指標を用いることで、ピークの有無を正確に判定することができる。Here, a method of simply determining the presence or absence of a peak by comparing power (power) and a power threshold value will be considered. FIG. 8B schematically shows a signal with a high background level (solid line) and a signal with a low background level (dotted line). When the background level is low as indicated by a dotted line in FIG. 8B, the peak can be detected by comparing the power peak value with the power threshold value Pth . However, as shown by the solid line in FIG. 8B, when the background level becomes high enough to exceed Pth , no peak can be detected even if the power peak value is compared with Pth . This is because the background level cannot be considered in a simple comparison between the power peak value and the power threshold value. Therefore, as shown in this embodiment, the presence or absence of a peak can be accurately determined by using an index that takes the background level into account for peak detection.
以上に示すように、本実施の形態に係る到来角度算出装置は、ピーク期間の電力とピーク期間以外の残りの期間における電力との比を求めて、この求めた比としきい値とを比較してピークの有無を判定することにより、受信波のバックグラウンドレベルが高い場合であっても希望波のピークを正確に検出して到来角度の計算に用いることができる。つまり、希望波以外の信号成分から到来角度を算出することが無いため、到来角度の算出精度を高めることができる。 As described above, the angle-of-arrival calculation device according to the present embodiment obtains a ratio between the power in the peak period and the power in the remaining period other than the peak period, and compares the obtained ratio with a threshold value. By determining whether or not there is a peak, the peak of the desired wave can be accurately detected and used for calculation of the arrival angle even when the background level of the received wave is high. That is, since the arrival angle is not calculated from signal components other than the desired wave, the calculation accuracy of the arrival angle can be improved.
図11は、到来角度算出装置1における到来角度算出部24の別の一態様を説明するブロック図である。図11に示される到来角度算出部24は、タイミング制御部23aの出力Oa1の複素共役をとる複素共役部51と、複素共役部51の出力Oa1´と、タイミング制御部23bの出力Oa2を複素乗算する複素乗算部52と、複素乗算部52の出力を用いて逆正接演算を行う逆正接部53とを備える。複素共役部51、複素乗算部52、逆正接部53の動作や機能は、上述の複素共役部41、複素乗算部42、逆正接部43の動作や機能と同様である。また、逆正接部53の演算結果(位相差)を元に演算結果を補正する位相差補正部54と、位相差補正部54の出力を平均化する平均化部55と、位相差補正部54において補正を行った場合に平均化部55の演算結果(平均値)を補正する位相差再補正部56と、位相差再補正部56の出力を用いて到来角度に変換する到来角度変換部57と、を備える。到来角度変換部57の動作や機能は、上述の到来角度変換部46の動作や機能と同様である。FIG. 11 is a block diagram illustrating another aspect of the arrival
位相差補正部54は、逆正接部53の演算結果である位相差が、+180°(+π)付近や−180°(−π)付近の値になる場合、逆正接部の演算結果に所定の角度(位相差)を加える処理を行う。図12のI−Q平面に示すように、本実施の形態の到来角度算出部24は、位相差を−180°〜+180°(−π〜+π)の位相差範囲の座標上に投影する。このため、例えば、図13(a)に示されるように、逆正接部53によって算出される位相差が+180°および−180°近傍の値にならない場合には、これを平均化することで適切に到来角度を算出することができる。しかし、図13(b)に示されるように、逆正接部53によって算出される位相差が+180および−180近傍の値になる場合、算出される位相差の僅かな誤差が角度算出に大きな影響を与えることになる。ここで、位相差データとして、−178°および+178°の二つの値が得られ、一方の値である+178°は本来の値である−178°から−4°の誤差が生じて+178°になっていると想定する。これらの差は、実際には僅かに4°である。しかし、平均化処理において、−178°と+178°として平均化すると、平均値は0°となる。実際には約180°の位相差が存在するにもかかわらず、平均化処理によって0°として扱われてしまうのである。このように、平均化された位相差が本来の位相差から大幅にずれてしまうと、適切な到来角度算出は困難になる。
When the phase difference that is the calculation result of the
そこで、図11に示される到来角度算出部24は、逆正接部53によって算出される位相差が+180°および−180°付近の値になる場合、位相差補正部54が逆正接部53の演算結果に所定の角度(位相差)を加える補正処理を行って、適切な平均化が行われるようにするのである。逆正接部53の演算結果が+180°または−180°近傍の値であるか否かは、逆正接部53の演算結果として得られる複数の位相差の分布を元に判定することができる。例えば、+90°(+π/2)より大きく、または−90°(−π/2)より小さい位相差の数が、+90°より小さくかつ−90°より大きくなる位相差の数より多い場合には、逆正接部53の演算結果が+180°および−180°付近の値であると判定できる。位相差補正部54が加える角度(位相差)は、例えば+90°とすることができるが、適切な平均化処理が可能な角度であればこれに限られない。好適には、−90°、+180°または−180°のいずれかでも良い。
Therefore, when the phase difference calculated by the
平均化部55は、位相差補正部54の出力を平均化する。本実施の形態の到来角度算出部24は、平均化に適さない位相差が算出される場合に位相差を加える補正を行うため、平均化部55において適切な平均化処理が可能である。位相差再補正部56は、位相差補正部54において位相差の補正を行っている場合に平均化部55の出力を補正する。具体的には、位相差補正部54において補正値として加えた角度(位相差)を減ずる補正を行う。
The averaging
図14に、位相差が+180°および−180°付近となる場合の到来角度算出の概略を模式的に示す。逆正接部53によって算出された位相差が、I−Q平面において+180°および−180°付近の場合、位相差補正部54は位相差に補正値(+90°)を加えて座標軸を回転させ、平均値算出用の座標軸に変換する。平均化部55は、当該データを元に平均値(−92°)を算出する。位相差再補正部56は、位相差補正部54の出力データから補正値(+90°)を減ずる補正を行い、逆正弦部57に補正されたデータ(+178°)を出力する。
FIG. 14 schematically shows an outline of arrival angle calculation when the phase difference is near + 180 ° and −180 °. When the phase difference calculated by the
図15は上記到来角度算出部24における処理フロー図である。到来角度算出部24の複素共役部51は、ステップ301において、タイミング制御部23aの出力Oa1の複素共役を算出する。また、複素乗算部52は、ステップ302において、タイミング制御部23bの出力Oa2と複素共役部51の出力Oa1´とを乗算する。そして、逆正接部53は、ステップ303において、複素乗算部52の出力を用いて逆正接演算を行い、受信信号間の位相差を算出する。FIG. 15 is a processing flowchart in the arrival
ステップ304において、位相差補正部54は、算出された位相差がI−Q平面において+180°および−180°近傍の値であるかを判定する。算出された位相差が+180°および−180°近傍の値でない場合はステップ305に進み、到来角度算出部24は位相差を補正することなく到来角度を算出する。算出された位相差が+180°近傍、または−180°近傍の値の場合はステップ306に進む。当該判定は、上述のように、+90°より大きく、または−90°より小さい位相差の数が、+90°より小さくかつ−90°より大きくなる位相差の数より多いかどうかを基準として行うことができる。
In
ステップ306において、位相差補正部54は、逆正接部53の演算結果である位相差に90°を加える処理を行う。ステップ307において、平均化部55は、位相差補正部54の出力を平均化する。そして、ステップ308において、位相差再補正部56は、平均化部55の演算結果である平均値から90°を減ずる処理を行う。その後、ステップ309において、到来角度変換部57は、位相差再補正部56の出力から到来角度を算出する。このように、図11に示される到来角度算出部24では、所定の位相差を加えて平均化した後に所定の位相差を減ずるという一連の処理によって適切な平均値が算出されるため、到来角度の算出精度が低下せずに済む。その結果、到来角度の算出精度を十分に高めることができる。
In
なお、ここでは、位相差補正部54が、逆正接部53の演算結果に所定の角度を加える処理を行っているが、適切な平均化処理が実現できるのであればこれに限られない。例えば、図16に示すような構成の到来角度算出部24を用いることもできる。図16に示す到来角度算出部24は、タイミング制御部23aの出力Oa1の複素共役をとる複素共役部61と、複素共役部61の出力Oa1´と、タイミング制御部23bの出力Oa2を複素乗算する複素乗算部62とを備える。複素共役部61、複素乗算部62の動作や機能は、上述の複素共役部41、複素乗算部42の動作や機能と同様である。また、複素乗算部62の出力の同相成分(I成分)の絶対値と直交成分(Q成分)の絶対値とを比較するIQ比較部63と、複素乗算部62の出力を用い、IQ比較部63の出力に応じて演算方法を選択、変更して逆正接演算を行う逆正接部64とを備える。また、逆正接部64の演算結果である位相差を平均化する平均化部65と、逆正接部64の演算方法に応じて平均化部65の演算結果である平均値を補正する位相差再補正部66と、位相差再補正部66の出力を用いて到来角度に変換する到来角度変換部67と、を備える。到来角度変換部67の動作や機能は、上述の到来角度変換部46の動作や機能と同様である。Here, the phase
IQ比較部63は、複素乗算部の出力の同相成分(I成分)が負であるか否かを判定すると共に、複素乗算部62の出力の同相成分(I成分)の絶対値と直交成分(Q成分)の絶対値とを比較する。具体的には、IQ比較部63は、同相成分Ibの符号を判定すると共に、同相成分の絶対値|Ib|が直交成分の絶対値|Qb|より十分に大きいか否か(直交成分の絶対値|Qb|が同相成分の絶対値|Ib|より十分に小さいか否か)を判定する。受信信号の位相差がI−Q平面において+180°および−180°近傍の値をとる場合には、同相成分Ibが負になり(Ib<0)、同相成分の絶対値|Ib|が直交成分の絶対値|Qb|より十分に大きくなる。このため、同相成分Ibの符号を判定し、同相成分の絶対値|Ib|が直交成分の絶対値|Qb|より十分に大きいか否かを判定することにより、位相差が+180°および−180°近傍の値をとるか否かを判定することができる。
The IQ comparison unit 63 determines whether or not the in-phase component (I component) of the output of the complex multiplication unit is negative, and the absolute value of the in-phase component (I component) of the output of the
逆正接部64は、複素乗算部62の出力を用い、IQ比較部63の出力に応じて演算方法を選択して逆正接演算を行う。同相成分が正である場合や、同相成分が負であり、かつ同相成分の絶対値|Ib|が直交成分の絶対値|Qb|と同程度であるか、または小さい場合、複素乗算部62の出力Ibを分母とし、出力Qbを分子とした値の逆正接演算を行う。同相成分が負であり、かつ同相成分の絶対値|Ib|が直交成分の絶対値|Qb|より十分に大きい場合、例えば、複素乗算部62の出力Qbの符号を反転させた−Qbを分母とし、出力Ibを分子とした値の逆正接演算を行う。なお、同相成分の絶対値|Ib|が直交成分の絶対値|Qb|より十分に大きい場合の上記処理は、座標軸を+90°回転させて逆正接演算を行う処理に相当する。つまり、当該処理によって得られる位相差は、元来の位相差に+90°が加えられた値である。
The
なお、同相成分の絶対値|Ib|が直交成分の絶対値|Qb|より十分に大きい場合の処理は、上述のものに限られない。例えば、複素乗算部62の出力Qbを分母とし、出力Ibの符号を反転させた−Ibを分子とした値の逆正接演算を行っても良い。当該処理は、座標軸を−90°回転させて逆正接演算を行う処理に相当する。つまり、当該処理によって得られる位相差は、元来の位相差に−90°が加えられた値(+90°が減じられた値)である。また、例えば、複素乗算部62の出力Ibの符号と、出力Qbの符号とを反転させて逆正接演算を行っても良い。当該処理は、座標軸を+180°(または−180°)回転させて逆正接演算を行う処理に相当する。つまり、当該処理によって得られる位相差は、元来の位相差に+180°(または−180°)が加えられた値である。このような処理によっても、適切な平均値を算出することができる。
Note that the processing when the absolute value | Ib | of the in-phase component is sufficiently larger than the absolute value | Qb | of the quadrature component is not limited to the above. For example, an arctangent operation may be performed on a value using the output Qb of the
平均化部65は、逆正接部64の出力を平均化する。本実施の形態の到来角度算出部24は、平均化に適さない位相差が算出される場合に実質的に位相差を加える(または減ずる)補正を行うため、平均化部65において適切な平均化処理が可能である。位相差再補正部66は、逆正接部64が座標軸を+90°回転させる処理を行っている場合に平均化部65の出力を補正する。具体的には、+90°を減ずる補正を行う。なお、逆正接部64が座標軸を−90°回転させる処理を行っている場合には、−90°を減ずる補正(つまり、+90°を加える補正)を行う。同様に、逆正接部64が座標軸を+180°(または−180°)回転させる処理を行っている場合には、+180°(または−180°)を減ずる補正を行う。
The averaging
このように、図16に示す到来角度算出部24も、図11に示される到来角度算出部24と同様に適切な平均値を算出できるため、到来角度の算出精度が低下せずに済む。その結果、到来角度の算出精度を十分に高めることができる。
In this manner, the arrival
図17は、変調方式として直交周波数分割多重(OFDM)を用いる場合の到来角度算出装置の具体的構成例を示すブロック図である。なお、図17では、図1における演算部13に相当する構成のみを示している。
FIG. 17 is a block diagram illustrating a specific configuration example of an arrival angle calculation apparatus when orthogonal frequency division multiplexing (OFDM) is used as a modulation scheme. Note that FIG. 17 shows only the configuration corresponding to the
図17において、相関処理部21aは、受信部12aの出力の複素共役をとる複素共役部71aと、受信部12aの出力を所定期間だけ遅延させて出力する遅延部72aと、複素共役部71aの出力と遅延部72aの出力とを複素乗算する複素乗算部73aと、複素乗算部73aの出力をGI(ガードインターバル)期間だけ足し合わせて出力する加算器74a、74bとを備える。ピーク検出部22aは、加算器74a、74bから出力された信号の電力を算出する電力算出部75aと、その電力ピークを検出してタイミング制御部23aに出力するピーク電力検出部76aとを備える。タイミング制御部23aは、ピーク電力検出部76aからの信号を元に受信部12aからの信号の到来角度算出部24への出力タイミングを制御する遅延部77aを備える。同様に、相関処理部21bは、複素共役部71b、遅延部72b、複素乗算部73b、加算器74c、74dを備え、ピーク検出部22bは、電力算出部75b、ピーク電力検出部76bを備え、タイミング制御部23bは遅延部77bを備える。到来角度算出部24は、遅延部77aの出力の複素共役をとる複素共役部81と、複素共役部81の出力と、遅延部77bの出力を複素乗算する複素乗算部82と、複素乗算部42の出力をGI(ガードインターバル)期間だけ足し合わせて出力する加算部83a、83bと、加算部83a、83bの出力を用いて逆正接演算を行う逆正接部84と、逆正接部84の出力を平均化する平均化部85と、平均化部85の出力を用いて到来角度に変換する到来角度変換部86とを備える。
In FIG. 17, the
遅延部72a、72bは、OFDMシンボル列の自己相関をとるため、受信部12aの出力を所定期間だけ遅延させて出力する。具体的には、遅延部72a、72bは、複素共役部71aが出力するOFDMシンボルの末部と、遅延部72a、72bが出力するGI(ガードインターバル)とが同じタイミングで複素乗算部73aに入力されるように、受信部12aの出力を所定期間だけ遅延させて出力する。複素乗算部73aは、複素共役部71aの出力と遅延部72aの出力とを複素乗算する。加算器74aおよび74bは、複素乗算部73aのチップ区間ごとの出力をGI期間だけ足し合わせて出力する。
The
図18(a)は、OFDMシンボル列の構成を示す模式図である。図18(a)に示すように、OFDMシンボル列は、データ部であるOFDMシンボルと、OFDMシンボルの先頭に配置されるGIとによって構成される。GIはOFDMシンボル末部をコピーしたデータであり、OFDMシンボル間の干渉を防ぐために挿入される。図18(b)は、相関処理部21aにおけるOFDMシンボル列の相関処理(自己相関処理)の様子を示す模式図である。図18(a)に示すように、遅延部72aの出力は、複素共役部71aの出力に対してOFDMシンボル長だけ遅れている。このため、複素乗算部73aにおいて、複素共役部71aの出力と遅延部72aの出力とを乗算することで自己相関をとることができる。自己相関値(GI相関値)は、複素共役部71aの出力と遅延部72aの出力にGIと同じデータが現れたときにピークを示すため、これを用いることで、データ部であるOFDMシンボルの先頭を検出することができる。
FIG. 18A is a schematic diagram showing a configuration of an OFDM symbol sequence. As shown in FIG. 18A, the OFDM symbol string is composed of an OFDM symbol that is a data part and a GI that is arranged at the head of the OFDM symbol. GI is data obtained by copying the end of the OFDM symbol, and is inserted to prevent interference between OFDM symbols. FIG. 18B is a schematic diagram illustrating a state of correlation processing (autocorrelation processing) of the OFDM symbol sequence in the
加算器74aおよび74bの出力信号は、ピーク検出部22aの電力算出部75aに入力される。電力算出部75aは、加算器74aおよび74bの出力信号からチップ区間ごとの電力を算出する。具体的には、電力算出部34aは、同相成分に相当する出力信号の絶対値と、直交成分に相当する出力信号の絶対値とを足し合わせ、チップ区間ごとの電力情報としてピーク電力検出部76aに出力する。なお、同相成分に相当する出力信号の2乗値と、直交成分に相当する出力信号の2乗値とを足し合わせてピーク電力検出部76aに出力しても良い。図19(a)に電力算出部75aからの出力波形の例を示す。図19(b)は、図19(a)に示す出力波形の部分拡大図である。ピーク電力検出部76aは、チップ区間ごとの電力情報を受け取ると、受信信号中の電力ピークを検出し、電力ピーク情報としてタイミング制御部23aの遅延部77aに出力する。
The output signals of the
ピーク検出部22a(ピーク電力検出部35a)から出力される電力ピーク情報は、受信信号のピークの有無を判定する情報である。具体的には、電力ピーク情報は、受信信号のピーク点付近の期間(ピーク期間)における電力の和ΣP1と、OFDMでの情報単位となる1シンボル期間からピーク期間を除いた期間における電力の和ΣP2との比R(=ΣP1/ΣP2)がしきい値Rthより大きいか否かを示す情報である。変調方式としてOFDMを用いる場合、ピーク期間はGI期間に等しくなる。また、1シンボル期間とは、GI期間とデータ期間(OFDMシンボル期間)とを合計した期間に相当する。電力ピーク情報において、RがRthより大きい場合には、タイミング制御部23a(遅延部77a)は、そのタイミングで受信信号がピークを有するものとして、受信部12aからの受信信号を到来角度算出部24に出力する。一方で、RがRthより小さい場合には、タイミング制御部23a(遅延部77a)は、そのタイミングでは受信信号がピークを有しないものとして、到来角度算出部24への出力を停止する。なお、ここでは、ピーク検出部22aが電力ピーク情報に関する演算処理を行っているが、電力ピーク情報に関する演算処理はタイミング制御部23aにおいて行っても良い。The power peak information output from the
相関処理部21b(複素共役部71b、遅延部72b、複素乗算部73b、加算器74c、74d)、ピーク検出部22b(電力算出部75b、ピーク電力検出部76b)、タイミング制御部23b(遅延部77b)の動作や機能は、相関処理部21a(複素共役部71a、遅延部72a、複素乗算部73a、加算器74a、74b)、ピーク検出部22a(電力算出部75a、ピーク電力検出部76a)、タイミング制御部23a(遅延部77a)の動作や機能と同様である。ただし、相関処理部21bに入力される受信信号と、相関処理部21aに入力される受信信号とは、同一電波を所定間隔離れた二点で受信しているため位相が僅かに異なっている。このため、タイミング制御部23bから出力される信号と、タイミング制御部23aから出力される信号とでは、位相が僅かに相違する。
タイミング制御部23aの出力は、到来角度算出部24の複素共役部81に入力される。複素共役部81は、タイミング制御部23aの出力の複素共役を複素乗算部82に出力する。複素乗算部82は、複素共役部81の出力と、タイミング制御部23bの出力とを複素乗算して、演算結果を加算部83aおよび83bに出力する。加算部83aおよび83bは、複素乗算部82のチップ区間ごとの出力をGI期間だけ足し合わせて逆正接部84に出力する。図19(c)に加算部83aおよび83bからの出力波形の例を示す。図中で、加算部83aの出力波形はIで示しており、加算部83bの出力波形はQで示している。
The output of the
逆正接部84は、加算部83aおよび83bの出力を用いて逆正接演算を行い、受信信号の位相差を算出する。図19(d)に逆正接部84からの出力波形の例を示す。平均化部85は、逆正接部84の出力を平均化して到来角度変換部86に出力する。なお、平均化部85は、適宜省略して良い。到来角度変換部86は、平均化部85の出力(平均化部85を有しない場合は逆正接部84の出力)を用いて逆三角関数演算により到来角度に変換する。当該演算によって求められる値、すなわち、到来角度変換部86の出力が、到来角度に相当する。
The
このように、図17の演算部13を有する到来角度算出装置1においても、ピーク期間の電力とピーク期間以外の残りの期間における電力との比を求めて、この求めた比としきい値とを比較してピークの有無を判定することにより、受信波のバックグラウンドレベルが高い場合であっても希望波のピークを正確に検出して、到来角度の計算に用いることができる。つまり、希望波以外の信号成分から到来角度を算出することが無いため、到来角度の算出精度を高めることができる。
As described above, the arrival
図20は、到来角度算出装置1をカプセル内視鏡の位置特定に応用したカプセル内視鏡システムについて示す模式図である。図20に示すカプセル内視鏡システムは、複数のセンサアレイ401と、センサアレイ401からのデータを記録するデータレコーダー402とを備える。センサアレイ401は、到来角度算出装置1の受信用アンテナに相当するアンテナを備えており、患者が飲み込んだカプセル内視鏡からの電波を受信できるように構成されている。データレコーダー402は、センサアレイ401において受信した電波の持つ位相情報から、患者が飲み込んだカプセル内視鏡の位置を特定する。
FIG. 20 is a schematic diagram showing a capsule endoscope system in which the arrival
患者が飲み込んだカプセル内視鏡は、消化管の蠕動運動によって移動する。カプセル内視鏡の位置はモニタされており、診察部位に到達したか否かを確認することができる。カプセル内視鏡が診察部位に到達すると、カプセル内視鏡は診察部位の様子を撮影してデータレコーダー402に送信し、データレコーダー402は画像情報を記録する。このように、カプセル内視鏡の位置をモニタすることで、診察部位を見逃すことなく撮影することができる。また、カプセル内視鏡が診察部位に到達したタイミングでカメラ等の電源を入れ、診察部位をはずれた場合にはカメラ等の電源を切る事が可能になるため、電池容量を小さくする事ができる。また、センサ(アンテナ)の数を削減する事が可能となる。また、電池容量が同じであれば、従来型のカプセル内視鏡と比較して多数の画像を送信でき、鮮明な画像を得ることができる。
The capsule endoscope swallowed by the patient moves by the peristaltic movement of the digestive tract. The position of the capsule endoscope is monitored, and it can be confirmed whether or not the examination site has been reached. When the capsule endoscope reaches the examination site, the capsule endoscope captures the state of the examination site and transmits it to the
このように、到来角度算出装置1をカプセル内視鏡の位置特定に応用することで、優れたカプセル内視鏡システムを構築することができる。
In this way, an excellent capsule endoscope system can be constructed by applying the arrival
以上のように、本発明の到来角度算出装置によれば、ピーク期間の電力とピーク期間以外の残りの期間における電力との比を求めて、この求めた比としきい値とを比較し、比がしきい値より大きい場合に到来角度を算出するため、受信波の希望波以外の信号レベルが高い場合であっても、希望波のピークを正確に検出し、到来角度を算出することができる。つまり、希望波以外の部分から到来角度を算出することが無いため、到来角度の算出精度を高めることができる。 As described above, according to the arrival angle calculation device of the present invention, the ratio between the power in the peak period and the power in the remaining period other than the peak period is obtained, and the obtained ratio is compared with the threshold value. Since the arrival angle is calculated when is greater than the threshold value, the peak of the desired wave can be accurately detected and the arrival angle can be calculated even when the signal level of the received wave other than the desired wave is high. . That is, since the arrival angle is not calculated from a portion other than the desired wave, the calculation accuracy of the arrival angle can be improved.
なお、本発明は上記実施の形態の記載に限定されず、その効果を発揮する態様で適宜変更することができる。例えば、上記実施の形態において、ピーク期間の電力の和と、ピーク期間を除いた期間の電力の和との比をしきい値と比較しているが、希望波以外の信号のレベルを考慮した到来角度算出が可能であればこれに限られない。例えば、ピーク期間のあるタイミングにおける電力と、ピーク期間を除いた期間のあるタイミングにおける電力とをパラメータとして用いても良い。 In addition, this invention is not limited to description of the said embodiment, In the aspect which exhibits the effect, it can change suitably. For example, in the above embodiment, the ratio of the sum of power during the peak period and the sum of power during the period excluding the peak period is compared with the threshold value, but the level of the signal other than the desired wave is considered. The present invention is not limited to this as long as the arrival angle can be calculated. For example, power at a certain timing in the peak period and power at a certain timing other than the peak period may be used as parameters.
また、上記実施の形態において、添付図面に示されている構成などは、これに限定されず、本発明の効果を発揮する範囲内で適宜変更することが可能である。 Moreover, in the said embodiment, the structure etc. which are shown by the attached drawing are not limited to this, It is possible to change suitably in the range which exhibits the effect of this invention.
本発明の到来角度算出装置は、対象の位置を特定するシステム、その他の各種用途に用いることができる。 The arrival angle calculation device of the present invention can be used for a system for specifying a target position and various other uses.
本出願は、2010年11月12日出願の特願2010−254011に基づく。この内容は、全てここに含めておく。 This application is based on Japanese Patent Application No. 2010-254011 of an application on November 12, 2010. All this content is included here.
Claims (9)
前記各受信信号処理部は、対応する前記アンテナで受信した電波を当該電波の位相情報を有する受信信号に変換して出力する受信部と、前記受信部から出力された受信信号を相関処理する相関処理部と、前記相関処理された受信信号のピークを検出するピーク検出部と、前記相関処理部の出力信号から前記受信信号処理部間で同一情報単位となる信号成分が切り出されるように、前記ピーク検出部で検出されたピークのタイミングに合わせて、前記相関処理部から出力される出力信号の取り込みタイミングを制御するタイミング制御部と、を備え、
前記タイミング制御部は、前記情報単位に相当する期間におけるピーク期間の電力と当該ピーク期間を除いた期間の電力との比が、しきい値より大きければ、前記相関処理部からの信号を前記到来角度算出部に出力することを特徴とする到来角度算出装置。A plurality of antennas for receiving radio waves transmitted from a certain position, a plurality of reception signal processing units provided corresponding to the respective antennas, and reception signal processing from output signals output from the plurality of reception signal processing units An angle-of-arrival calculation unit that takes in a signal component that is the same information unit between the units and calculates the angle of arrival of the radio wave,
Each reception signal processing unit performs correlation processing on a reception unit that converts a radio wave received by the corresponding antenna into a reception signal having phase information of the radio wave, and a reception signal output from the reception unit. A processing unit, a peak detection unit that detects a peak of the correlation-processed reception signal, and a signal component that is the same information unit between the reception signal processing units from the output signal of the correlation processing unit, In accordance with the timing of the peak detected by the peak detection unit, a timing control unit that controls the capture timing of the output signal output from the correlation processing unit, and
The timing control unit receives the signal from the correlation processing unit if the ratio between the power in the peak period in the period corresponding to the information unit and the power in the period excluding the peak period is greater than a threshold value. An angle-of-arrival calculation apparatus that outputs to an angle calculation unit.
一方のアンテナに対応した一方の受信信号処理部のタイミング制御部からの信号の複素共役をとる複素共役部と、
前記複素共役部の出力と、他方のアンテナに対応した他方の受信信号処理部のタイミング制御部からの信号とを乗算する複素乗算部と、
前記複素乗算部の出力を用いて逆正接演算を行い、前記アンテナ間での前記受信電波の位相差を算出する逆正接部と、
前記逆正接部の出力を平均化する平均化部と、
前記平均化部の出力を用いて逆三角関数演算を行い、到来角度に変換する到来角度変換部と、を備えることを特徴とする請求項1または請求項2に記載の到来角度算出装置。The arrival angle calculation unit
A complex conjugate unit that takes a complex conjugate of the signal from the timing control unit of one received signal processing unit corresponding to one antenna;
A complex multiplier that multiplies the output of the complex conjugate unit by the signal from the timing control unit of the other received signal processing unit corresponding to the other antenna;
An arc tangent is calculated using the output of the complex multiplier, and an arc tangent for calculating a phase difference of the received radio wave between the antennas;
An averaging unit that averages the output of the arctangent,
The arrival angle calculation device according to claim 1, further comprising: an arrival angle conversion unit that performs an inverse trigonometric function operation using an output of the averaging unit and converts the result into an arrival angle.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012542831A JP5677452B2 (en) | 2010-11-12 | 2011-08-10 | Arrival angle calculation device |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010254011 | 2010-11-12 | ||
JP2010254011 | 2010-11-12 | ||
PCT/JP2011/068281 WO2012063532A1 (en) | 2010-11-12 | 2011-08-10 | Arrival angle calculation device |
JP2012542831A JP5677452B2 (en) | 2010-11-12 | 2011-08-10 | Arrival angle calculation device |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2012063532A1 true JPWO2012063532A1 (en) | 2014-05-12 |
JP5677452B2 JP5677452B2 (en) | 2015-02-25 |
Family
ID=46050690
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012542831A Active JP5677452B2 (en) | 2010-11-12 | 2011-08-10 | Arrival angle calculation device |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP5677452B2 (en) |
CN (1) | CN103180751B (en) |
WO (1) | WO2012063532A1 (en) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2895222B1 (en) * | 2012-09-12 | 2018-12-19 | Silverbow Development LLC | Nasal interface apparatus and systems for use with a respiratory assist device |
US11904096B2 (en) | 2012-09-12 | 2024-02-20 | Inogen, Inc. | Nasal interface apparatus and systems for use with a respiratory assist device |
US9706514B2 (en) | 2014-12-02 | 2017-07-11 | Cisco Technology, Inc. | Wideband angle-of-arrival location determination using bandwidth partitioning |
US10539645B2 (en) * | 2016-01-22 | 2020-01-21 | GM Global Technology Operations LLC | Angle of arrival estimation |
US10509093B2 (en) * | 2016-12-30 | 2019-12-17 | Futurewei Technologies, Inc. | System and method for fast system acquisition and channel estimation |
US10330770B2 (en) | 2017-11-09 | 2019-06-25 | Cisco Technology, Inc. | Channel estimation in OFDMA for switched antenna array based angle-of-arrival location |
JP7114315B2 (en) * | 2018-04-19 | 2022-08-08 | 日本電産サンキョー株式会社 | encoder |
US11095017B2 (en) * | 2018-07-13 | 2021-08-17 | Apple Inc. | Electronic device having angle of arrival detection capabilities |
JP7070243B2 (en) * | 2018-08-27 | 2022-05-18 | 沖電気工業株式会社 | Arrival direction estimation device |
EP3745155A1 (en) | 2019-05-29 | 2020-12-02 | Assa Abloy AB | Determining a position of a mobile key device based on phase difference of samples |
CN111148020B (en) * | 2019-12-30 | 2022-08-12 | 上海美迪索科电子科技有限公司 | Positioning system, method, device and computer readable storage medium |
CN112468962B (en) * | 2020-11-27 | 2023-06-06 | Oppo广东移动通信有限公司 | Angle measuring method, terminal and computer storage medium |
CN114487993B (en) * | 2022-04-18 | 2022-08-02 | 湖南艾科诺维科技有限公司 | Direction finding method and system for uniform circular array correlation interferometer |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6266180A (en) * | 1985-09-19 | 1987-03-25 | Nec Corp | Radar response correlating circuit |
US5406286A (en) * | 1992-11-17 | 1995-04-11 | Honeywell Inc. | Real time passive threat positioning system |
JPH07159521A (en) * | 1993-12-10 | 1995-06-23 | Nec Eng Ltd | Sonar image display processing device |
US6108555A (en) * | 1996-05-17 | 2000-08-22 | Ksi, Inc. | Enchanced time difference localization system |
JP3214481B2 (en) * | 1999-01-28 | 2001-10-02 | 日本電気株式会社 | Target detection device and target detection method |
JP3004015B1 (en) * | 1999-03-18 | 2000-01-31 | 株式会社ワイ・アール・ピー移動通信基盤技術研究所 | Arriving wave measuring method and arriving wave measuring device |
JP2001196835A (en) * | 2000-01-17 | 2001-07-19 | Matsushita Electric Ind Co Ltd | Incoming direction estimating method and radio receiver |
JP2009014688A (en) * | 2007-07-09 | 2009-01-22 | Nippon Telegr & Teleph Corp <Ntt> | Radio wave arrival direction estimating device |
-
2011
- 2011-08-10 JP JP2012542831A patent/JP5677452B2/en active Active
- 2011-08-10 WO PCT/JP2011/068281 patent/WO2012063532A1/en active Application Filing
- 2011-08-10 CN CN201180051904.1A patent/CN103180751B/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP5677452B2 (en) | 2015-02-25 |
CN103180751B (en) | 2015-01-21 |
WO2012063532A1 (en) | 2012-05-18 |
CN103180751A (en) | 2013-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5677452B2 (en) | Arrival angle calculation device | |
JP5860068B2 (en) | Multipath detection method and arrival angle calculation apparatus | |
WO2019038542A2 (en) | Position determining system | |
EP2681582A1 (en) | Method, apparatus, and computer program product for processing received signals for locating | |
US10732273B2 (en) | Radar device for vehicle and method for estimating angle of target using same | |
US20190098507A1 (en) | Secure Channel Estimation Architecture | |
JP5729846B2 (en) | Multipath detection method, multipath detection apparatus, and arrival angle calculation apparatus | |
Liu et al. | Amt: Acoustic multi-target tracking with smartphone mimo system | |
Liu et al. | PAMT: Phase-based acoustic motion tracking in multipath fading environments | |
JP4486950B2 (en) | OQPSK demodulator timing estimator | |
WO2006065084A1 (en) | Synchronization method and apparatus and location awareness method and apparatus in chaotic communication system | |
CN115314986A (en) | Radio receiver synchronization | |
US9191254B2 (en) | Method and apparatus for detecting envelope using difference between sampling signals | |
JP2011053117A (en) | Wireless communication device | |
US20110032974A1 (en) | Carrier frequency synchronization detection circuit and correlation calculator | |
JP2005331466A5 (en) | ||
Shi et al. | Threshold-based resampling for high-speed particle PHD filter | |
KR20170128899A (en) | Apparatus And Method for Calculating Receiving Time Of Wireless Communication Signal | |
CN104246438A (en) | Method and system for energy efficient measurement of sensor signals | |
JP2009250704A (en) | Receiver for distance measurement device and distance measurement device | |
JP2008245253A (en) | Iterative pilot aided frequency offset estimation and c/i measurement for tdma signal with offset larger than nyquist frequency of the reference symbol rate | |
WO2016042606A1 (en) | Distance detection system | |
JP6090036B2 (en) | Spread spectrum receiver | |
US8559483B1 (en) | Signal quality monitor | |
KR102547936B1 (en) | Method of detecting information of target object by using ultrasound and electronic apparatus thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140328 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140526 |
|
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: 20141201 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20141226 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5677452 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |