JP7638460B2 - Radio wave arrival angle estimation device, radio wave arrival angle estimation method, control circuit, and storage medium - Google Patents
Radio wave arrival angle estimation device, radio wave arrival angle estimation method, control circuit, and storage medium Download PDFInfo
- Publication number
- JP7638460B2 JP7638460B2 JP2024568159A JP2024568159A JP7638460B2 JP 7638460 B2 JP7638460 B2 JP 7638460B2 JP 2024568159 A JP2024568159 A JP 2024568159A JP 2024568159 A JP2024568159 A JP 2024568159A JP 7638460 B2 JP7638460 B2 JP 7638460B2
- Authority
- JP
- Japan
- Prior art keywords
- interference
- frequency channel
- frequency
- unit
- arrival angle
- 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.)
- Active
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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Noise Elimination (AREA)
- Radar Systems Or Details Thereof (AREA)
Description
本開示は、複数の周波数チャネルを使用する電波到来角推定装置、電波到来角推定方法、制御回路および記憶媒体に関する。 The present disclosure relates to a radio wave arrival angle estimation device, a radio wave arrival angle estimation method, a control circuit, and a storage medium that use multiple frequency channels.
近年、電波の到来角を受信側で推定する技術が利用されている。電波の到来角の推定結果は、例えば、移動局の位置推定のためや、送信アンテナと受信アンテナとの間の遅延プロファイルを取得するためなどに使用される。Bluetooth(登録商標)では、電波の到来角から位置推定を行うAoA(Angle of Arrival)、AoD(Angle of Departure)といった技術が使用されることがある。電波の到来角を推定する際、屋内では、電波が壁面などに反射して受信側では複数の経路で到達するマルチパスの影響を受けやすく、電波の到来角の推定精度が低下する場合がある。また、他システムの干渉を受ける環境では、電波の到来角の推定精度が低下してしまう。In recent years, a technology has been used to estimate the angle of arrival of radio waves at the receiving side. The estimated result of the angle of arrival of radio waves is used, for example, to estimate the position of a mobile station or to obtain a delay profile between a transmitting antenna and a receiving antenna. In Bluetooth (registered trademark), technologies such as AoA (Angle of Arrival) and AoD (Angle of Departure) that estimate the position from the angle of arrival of radio waves are sometimes used. When estimating the angle of arrival of radio waves, indoors, the radio waves are easily affected by multipath, where the radio waves are reflected off walls and reach the receiving side via multiple paths, and the accuracy of the estimation of the angle of arrival of radio waves may decrease. In addition, in an environment where interference from other systems occurs, the accuracy of the estimation of the angle of arrival of radio waves decreases.
マルチパスの対策として、非特許文献1には、複数の周波数チャネルを用いて既知信号を送信し、受信した既知信号からチャネルインパルス応答を求め、チャネルインパルス応答から直接波成分と遅延波成分とを分離する技術が開示されている。このとき受信側で既知信号を周波数結合することで、1つの周波数チャネルの周波数帯域幅よりも広い周波数帯域幅でチャネルインパルス応答を求めることができる。チャネルインパルス応答の時間分解能は周波数帯域幅に反比例することから、周波数帯域幅が広くなるとチャネルインパルス応答の時間分解能が高くなり、直接波成分と遅延波成分との分離精度が向上するため、遅延波成分を抑圧することが可能になる。As a countermeasure against multipath, Non-Patent
干渉の対策として、特許文献1には、干渉の到来方向を推定して、複数の受信アンテナを用いたビームパターン制御を利用して、推定した到来方向にビームパターンでヌルを向けることにより、干渉の影響を軽減することができる。As a countermeasure against interference,
しかしながら、マルチパスおよび他システムからの干渉の両方の影響を受ける環境では、電波の到来角を高精度に推定することが困難な場合があるという問題があった。例えば、周波数チャネル毎に異なる干渉源からの干渉を受ける可能性がある場合、周波数チャネル毎に干渉抑圧ウェイトを適用して干渉抑圧処理を行うことが考えられる。しかしながら、異なる干渉抑圧ウェイトで干渉抑圧処理された周波数チャネルのチャネルインパルス応答は、干渉抑圧ウェイトに依存してしまう。このため、複数の周波数チャネルで周波数結合をすることができなくなってしまう。 However, in an environment affected by both multipath and interference from other systems, there is a problem in that it can be difficult to estimate the angle of arrival of radio waves with high accuracy. For example, when there is a possibility that each frequency channel may be subject to interference from a different interference source, it is possible to perform interference suppression processing by applying an interference suppression weight to each frequency channel. However, the channel impulse response of frequency channels that have been subjected to interference suppression processing with different interference suppression weights depends on the interference suppression weight. As a result, frequency coupling cannot be performed across multiple frequency channels.
本開示は、上記に鑑みてなされたものであって、マルチパスおよび他システムからの干渉の両方の影響を受ける環境であっても、電波の到来角を高精度に推定することが可能な電波到来角推定装置を得ることを目的とする。The present disclosure has been made in consideration of the above, and aims to provide a radio wave arrival angle estimation device that can estimate the arrival angle of radio waves with high accuracy even in an environment affected by both multipath and interference from other systems.
上述した課題を解決し、目的を達成するために、本開示にかかる電波到来角推定装置は、複数の周波数チャネルを使用して送信された信号を受信した受信信号から電波の到来方向を推定する電波到来角推定装置であって、周波数チャネル毎に干渉の状態を測定して干渉情報を生成する干渉測定部と、周波数チャネル毎に生成された干渉情報に基づいて、複数の周波数チャネルをグループ分けして周波数チャネル群を生成する周波数チャネル群生成部と、周波数チャネル群単位で干渉抑圧ウェイトを生成して干渉抑圧処理を行う干渉抑圧処理部と、干渉抑圧処理を行った後の信号に対して、周波数チャネル群単位で周波数チャネルを結合する周波数チャネル結合部と、周波数チャネルを結合した後の信号に対して、直接波成分を抽出する直接波成分抽出部と、直接波成分から電波の到来角を推定する電波到来角推定部と、を備えることを特徴とする。In order to solve the above-mentioned problems and achieve the objective, the radio wave arrival angle estimation device disclosed herein is a radio wave arrival angle estimation device that estimates the arrival direction of radio waves from a received signal that receives a signal transmitted using multiple frequency channels, and is characterized by comprising an interference measurement unit that measures the interference state for each frequency channel and generates interference information, a frequency channel group generation unit that groups the multiple frequency channels and generates frequency channel groups based on the interference information generated for each frequency channel, an interference suppression processing unit that generates interference suppression weights on a frequency channel group basis and performs interference suppression processing, a frequency channel combining unit that combines frequency channels on a frequency channel group basis for the signal after the interference suppression processing has been performed, a direct wave component extraction unit that extracts a direct wave component from the signal after the frequency channels have been combined, and a radio wave arrival angle estimation unit that estimates the arrival angle of radio waves from the direct wave component.
本開示にかかる電波到来角推定装置は、マルチパスおよび他システムからの干渉の両方の影響を受ける環境であっても、電波の到来角を高精度に推定することが可能であるという効果を奏する。The radio wave arrival angle estimation device disclosed herein has the advantage of being able to estimate the arrival angle of radio waves with high accuracy even in an environment affected by both multipath and interference from other systems.
以下に、本開示の実施の形態にかかる電波到来角推定装置、電波到来角推定方法、制御回路および記憶媒体を図面に基づいて詳細に説明する。 Below, the radio wave arrival angle estimation device, radio wave arrival angle estimation method, control circuit, and storage medium relating to the embodiments of the present disclosure are described in detail with reference to the drawings.
実施の形態1.
図1は、実施の形態1にかかる電波到来角推定装置1の構成を示す図である。電波到来角推定装置1は、送信装置10と、受信装置20とを有する。電波到来角推定装置1は、電波到来角推定システムと呼ぶこともできる。電波到来角推定装置1は、複数の周波数チャネルを使用して通信を行う。電波到来角推定装置1において、送信装置10は、複数の周波数チャネルを使用して信号を電波として送信し、受信装置20は、送信装置10が送信した信号を受信し、受信信号から電波の到来角を推定する。
1 is a diagram showing a configuration of a radio wave arrival
図2は、実施の形態1にかかる送信装置10の構成を示す図である。送信装置10は、既知信号生成部100と、変調部101と、送信フィルタ部102と、周波数変換部103と、送信アンテナ部104とを有する。2 is a diagram showing the configuration of the transmitting
既知信号生成部100は、受信装置20がチャネルインパルス応答を推定するために使用する既知信号を生成し、生成した既知信号を変調部101へ出力する。変調部101は、既知信号生成部100が出力する既知信号を変調し、変調された既知信号を送信フィルタ部102へ出力する。送信フィルタ部102は、変調部101が出力する変調後の既知信号を送信用のフィルタでフィルタリングし、フィルタリングされた信号を周波数変換部103へ出力する。周波数変換部103は、送信フィルタ部102が出力するフィルタリング後の信号の周波数を、複数の周波数チャネルのうち送信する周波数チャネルに合わせて変換し、周波数変換後の信号を送信アンテナ部104へ出力する。送信アンテナ部104は、周波数変換部103が出力する周波数変換後の信号を電波として送信する。The known signal generating unit 100 generates a known signal used by the
図3は、実施の形態1にかかる受信装置20の構成を示す図である。受信装置20は、受信アンテナ部200と、周波数変換部201と、受信フィルタ部202と、受信信号蓄積メモリ部203と、干渉測定部204と、周波数チャネル群生成部205と、干渉抑圧処理部206と、周波数チャネル結合部207と、直接波成分抽出部208と、電波到来角推定部209と、到来角合成部210とを有する。3 is a diagram showing the configuration of the
受信アンテナ部200は、送信装置10が電波として送信した信号を1つ以上の受信アンテナで受信し、受信信号を周波数変換部201へ出力する。周波数変換部201は、受信アンテナ部200が出力する受信信号をベースバンド信号に周波数変換し、周波数変換後の受信信号を受信フィルタ部202へ出力する。受信フィルタ部202は、周波数変換部201が出力する周波数変換後の受信信号を受信用のフィルタで通信波帯域にフィルタリングし、フィルタリング後の受信信号を受信信号蓄積メモリ部203および干渉測定部204のそれぞれへ出力する。The receiving
受信信号蓄積メモリ部203は、受信フィルタ部202が出力するフィルタリング後の受信信号を一時記憶する。干渉測定部204は、受信フィルタ部202が出力するフィルタリング後の受信信号から観測される干渉波を測定して干渉の状態を示す干渉情報を生成し、生成した干渉情報を周波数チャネル群生成部205および到来角合成部210のそれぞれに出力する。The received signal
周波数チャネル群生成部205は、干渉測定部204が生成した干渉情報に基づいて、電波到来角推定装置1が使用する複数の周波数チャネルをグループ分けして、1または複数の周波数チャネルから成る周波数チャネル群を生成する。周波数チャネル群生成部205は、生成した周波数チャネル群を示す情報を干渉抑圧処理部206に出力する。The frequency channel
干渉抑圧処理部206は、周波数チャネル群生成部205で生成した周波数チャネル群を示す情報に基づいて、受信信号蓄積メモリ部203に蓄積された受信信号を読み出し、周波数チャネル群毎に周波数チャネル群に属する受信信号から干渉抑圧ウェイトを生成して干渉抑圧処理を行い、干渉抑圧処理後の信号を周波数チャネル結合部207へ出力する。The interference
周波数チャネル結合部207は、干渉抑圧処理部206が出力する干渉抑圧処理後の信号に対して、周波数チャネル群毎に、周波数チャネルの結合を行い、結合後の信号を直接波成分抽出部208へ出力する。The frequency
直接波成分抽出部208は、周波数チャネル結合部207が出力する結合後の信号に対して周波数チャネル群毎にチャネルインパルス応答を求めて、チャネルインパルス応答に基づいて直接波成分を抽出し、抽出した直接波成分に相当するチャネルインパルス応答を電波到来角推定部209へ出力する。The direct wave
電波到来角推定部209は、直接波成分抽出部208が出力する直接波成分に相当するチャネルインパルス応答から電波の到来角を推定し、周波数チャネル群毎の電波の到来角を到来角合成部210へ出力する。The radio wave arrival
到来角合成部210は、電波到来角推定部209が出力する周波数チャネル群毎の電波の到来角の推定結果を、干渉測定部204が出力する干渉情報に基づいて合成して到来角推定結果を出力する。The arrival angle synthesis unit 210 synthesizes the estimated results of the arrival angle of radio waves for each frequency channel group output by the radio wave arrival
ここで、干渉測定部204の詳細な処理について説明する。
Here, we will explain the detailed processing of the
図4は、図3に示す干渉測定部204の詳細な構成例を示す図である。干渉測定部204は、SINR(Signal Interference Noise Ratio)推定部211と、干渉電力推定部212と、干渉到来角推定部213と、干渉複素信号推定部214とを有する。干渉測定部204の各部には、受信フィルタ部202の出力する受信信号が入力される。
Figure 4 is a diagram showing a detailed configuration example of the
SINR推定部211は、受信信号からSINRを推定し、SINR推定値を干渉情報として出力する。SINRは、信号対干渉雑音電力比とも呼ばれ、信号電力に対する干渉電力および雑音電力の合計の比率を表す。以下、干渉電力と雑音電力との合計のことを干渉雑音電力と称する。例えば、SINR推定部211は、既知信号を参照系列として、受信信号に対して参照系列による相関電力を算出することで、信号電力Si,mを推定することができる。ここでiは、周波数チャネルインデックスであり、Nfを周波数チャネル数とした場合、0以上Nf未満の整数である。またmは、受信アンテナインデックスであり、Mを受信アンテナ本数とした場合、0以上M未満の整数である。干渉雑音電力INi,mに関しては、SINR推定部211は、例えば、既知信号に用いる変調方式がFSK(Frequency Shift Keying)変調であれば、FSK変調における周波数ビンのうち、信号を送信しない周波数ビンである無信号領域を使用することで、干渉雑音電力INi,mを推定することができる。SINR推定部211は、信号電力Si,mおよび干渉雑音電力INi,mのそれぞれを推定すると、以下の数式(1)を用いてi番目の周波数チャネルについてのSINRであるSINRiを推定することができる。 The SINR estimation unit 211 estimates the SINR from the received signal and outputs the SINR estimation value as interference information. The SINR is also called the signal-to-interference-plus-noise power ratio, and represents the ratio of the total of the interference power and the noise power to the signal power. Hereinafter, the total of the interference power and the noise power is referred to as the interference plus noise power. For example, the SINR estimation unit 211 can estimate the signal power S i,m by using a known signal as a reference sequence and calculating the correlation power of the reference sequence for the received signal. Here, i is a frequency channel index, and is an integer equal to or greater than 0 and less than N f when N f is the number of frequency channels. Also, m is a receiving antenna index, and is an integer equal to or greater than 0 and less than M when M is the number of receiving antennas. Regarding the interference plus noise power IN i,m , if the modulation method used for the known signal is, for example, FSK (Frequency Shift Keying) modulation, the SINR estimation unit 211 can estimate the interference plus noise power IN i,m by using a non-signal region, which is a frequency bin in the FSK modulation that does not transmit a signal. When the SINR estimator 211 estimates each of the signal power S i,m and the interference noise power IN i,m , it can estimate SINR i , which is the SINR for the i-th frequency channel, using the following equation (1).
なお、SINR推定部211がSINRを推定する方法については、上述の方法に限定されるものではない。 Note that the method by which the SINR estimation unit 211 estimates the SINR is not limited to the above-mentioned method.
干渉電力推定部212は、受信信号に含まれる干渉電力を推定し、干渉電力推定値を干渉情報として出力する。干渉電力推定部212は、SINR推定部211と同様の方法で推定した干渉雑音電力INi,mから干渉電力Iiを推定することができる。既にSINR推定部211が推定した干渉雑音電力INi,mの値がある場合、干渉電力推定部212は、SINR推定部211から干渉雑音電力INi,mを取得して用いてもよい。干渉電力推定部212は、以下の数式(2)を用いて、M本のアンテナの干渉雑音電力INi,mの平均を、干渉電力Iiとすることができる。
The interference
干渉到来角推定部213は、受信信号に含まれる干渉波の到来角である干渉到来角を推定し、干渉到来角推定値を干渉情報として出力する。干渉到来角推定部213は、受信アンテナ毎の無信号領域を使用することで、受信アンテナ毎の干渉複素信号を測定し、測定した干渉複素信号を用いて、干渉到来角を推定することができる。受信アンテナが2本以上ある場合、受信アンテナmと受信アンテナlとの間隔dm,lと、受信アンテナmと受信アンテナlとの間の干渉における位相差φi,m,lと、i番目の周波数チャネルに対応する波長λiとを用いて、干渉到来角θi,m,lは、以下の数式(3)を用いて表される。位相差φi,m,lは、受信アンテナmの干渉複素信号Ci,m(n)と、受信アンテナlの干渉複素信号Ci,l(n)とから、以下の数式(4)を用いて求めることができる。なお、X*は、Xの複素共役を表す。なお、nは、無信号領域におけるシンボルインデックスであり、無信号領域におけるシンボル数をNsとすると、0以上Ns未満の整数である。なお、干渉到来角推定部213が干渉到来角θi,m,lを推定する方法は、上記の方法に限定されるものではない。 The interference arrival angle estimator 213 estimates an interference arrival angle, which is an arrival angle of an interference wave included in a received signal, and outputs an interference arrival angle estimate value as interference information. The interference arrival angle estimator 213 measures an interference complex signal for each receiving antenna by using a signal-free region for each receiving antenna, and can estimate an interference arrival angle using the measured interference complex signal. When there are two or more receiving antennas, the interference arrival angle θ i,m ,l is expressed using the following formula (3) using the interval d m,l between receiving antenna m and receiving antenna l, the phase difference φ i, m,l in the interference between receiving antenna m and receiving antenna l, and the wavelength λ i corresponding to the i-th frequency channel. The phase difference φ i,m,l can be obtained from the interference complex signal C i,m (n) of receiving antenna m and the interference complex signal C i,l (n) of receiving antenna l using the following formula (4). Note that X * represents the complex conjugate of X. Here, n is a symbol index in a no-signal region, and is an integer equal to or greater than 0 and less than Ns , where Ns is the number of symbols in a no-signal region. Note that the method by which the interference arrival angle estimator 213 estimates the interference arrival angle θ i,m,l is not limited to the above method.
干渉複素信号推定部214は、受信アンテナ毎の干渉複素信号を推定し、干渉複素信号推定値を干渉情報として出力する。干渉複素信号推定部214は、干渉到来角推定部213と同様の方法で干渉複素信号を測定することができる。既に干渉到来角推定部213が干渉複素信号を測定している場合、干渉複素信号推定部214は、干渉到来角推定部213から干渉複素信号を取得して、干渉情報として出力することができる。The interference complex signal estimation unit 214 estimates the interference complex signal for each receiving antenna and outputs the interference complex signal estimate value as interference information. The interference complex signal estimation unit 214 can measure the interference complex signal in a manner similar to that of the interference arrival angle estimation unit 213. If the interference arrival angle estimation unit 213 has already measured the interference complex signal, the interference complex signal estimation unit 214 can acquire the interference complex signal from the interference arrival angle estimation unit 213 and output it as interference information.
干渉測定部204は、SINR推定値、干渉電力推定値、干渉到来角推定値、および干渉複素信号推定値のうちの少なくとも1つを含む干渉情報を生成して、周波数チャネル群生成部205および到来角合成部210のそれぞれに出力する。The
続いて周波数チャネル群生成部205の詳細な構成について説明する。
Next, we will explain the detailed configuration of the frequency channel
図5は、図3に示す周波数チャネル群生成部205の詳細な構成例を示す図である。周波数チャネル群生成部205は、干渉クラスタリング部215と、周波数チャネル群決定部216とを有する。周波数チャネル群生成部205には、干渉測定部204が出力する干渉情報として、SINR推定値、干渉電力推定値、干渉到来角推定値および干渉複素信号推定値が入力される。これらの干渉情報は、干渉クラスタリング部215に入力される。また、干渉情報のうちSINR推定値は、周波数チャネル群決定部216にも入力される。周波数チャネル群生成部205は、SINR推定値、干渉電力推定値、干渉到来角推定値および干渉複素信号推定値のうちいずれか1つ以上を使用して、周波数チャネル群を生成する。
Figure 5 is a diagram showing a detailed configuration example of the frequency channel
ここで、送信装置10が送信する既知信号について説明する。
Here, we will explain the known signal transmitted by the transmitting
図6は、16個の周波数チャネルを用いて送信される既知信号を示す図である。既知信号は、周波数チャネルインデックス#0~#15の16個の周波数チャネルを用いて、通信スロット300として送信される。図6の横軸は周波数であり、縦軸は時間である。周波数ホッピング方式であれば、使用する周波数チャネルインデックスの順番は、周波数ホッピング系列として定義され、パケット単位で使用する周波数チャネルを変更しながら送信する。図6に示す例では、周波数ホッピング系列は、#0,#4,#8,#12,#1,#5,#9,#13,#2,#6,#10,#14,#3,#7,#11,#15である。なお、ここでは16個の周波数チャネルの例を示しているが、例えば、Bluetoothの場合、周波数チャネルは#0~#39の40チャネルであり、1チャネルあたり2MHzの周波数帯域幅を有している。
Figure 6 is a diagram showing a known signal transmitted using 16 frequency channels. The known signal is transmitted as a communication slot 300 using 16 frequency channels with frequency
図7は、図6に示す通信スロット300を周波数スペクトラムとして表した図である。信号301は、図6の通信スロット300の部分だけを周波数変換したスペクトラムであり、周波数チャネル毎の通信スロット300に対してスペクトラムが存在する。図7の横軸は周波数であり、縦軸は電力である。なお、図6に示す各通信スロット300は異なる時間で送信されているため、図7に示す周波数スペクトラムは、異なる時間の信号を周波数スペクトラムとして並べたものである。 Figure 7 shows the communication slot 300 shown in Figure 6 as a frequency spectrum. Signal 301 is a spectrum obtained by frequency converting only the communication slot 300 portion of Figure 6, and a spectrum exists for the communication slot 300 for each frequency channel. The horizontal axis of Figure 7 is frequency, and the vertical axis is power. Note that since each communication slot 300 shown in Figure 6 is transmitted at a different time, the frequency spectrum shown in Figure 7 is a frequency spectrum of signals at different times.
図8は、図6に示す既知信号が干渉を受けた場合の周波数スペクトラムの一例を示す図である。図8の横軸は周波数であり、縦軸は電力である。干渉波302Aは、周波数チャネルインデックス#1~#7の通信スロットに干渉しており、干渉波302Bは、周波数チャネルインデックス#11~#15の通信スロットに干渉している。それ以外の周波数チャネルインデックス#0,#8~#10の通信スロットは、干渉の影響を受けていない状態である。ここで、干渉波302Aと干渉波302Bとは異なる干渉源の干渉波であり、干渉波302Aと干渉波302Bとは、干渉電力、干渉到来角などが共通するものではない。
Fig. 8 is a diagram showing an example of a frequency spectrum when the known signal shown in Fig. 6 is subjected to interference. The horizontal axis of Fig. 8 is frequency, and the vertical axis is power. An interference wave 302A interferes with communication slots of frequency
周波数チャネル群生成部205における周波数チャネル群の決定方法について、図8に示す干渉を受けた場合を例に説明する。The method of determining the frequency channel group in the frequency channel
干渉クラスタリング部215は、SINR推定値、干渉電力推定値、干渉到来角推定値、および干渉複素信号推定値のいずれか1つ以上を用いて、干渉の有無と、干渉源とに基づいて、周波数チャネルをグループ分けする。具体的には、干渉クラスタリング部215は、干渉を受けた周波数チャネルと、干渉を受けていない周波数チャネルとをグループ分けする。さらに、干渉クラスタリング部215は、干渉を受けた周波数チャネルに関しては、干渉源毎のグループに分ける。干渉クラスタリング部215は、例えば、クラスタリングを使用してグループ分けを行うことができる。クラスタリングは、データの類似度に基づいてそれぞれのデータをグループ分けする手法であり、機械学習における教師なし学習に相当する。クラスタリングには、階層的クラスタリング、非階層的クラスタリングがある。階層的クラスタリングの手法の一例としては、群平均法、ウォード法などが挙げられる。非階層的クラスタリングの手法の一例としては、k-means法が挙げられる。ここでは、干渉クラスタリング部215が使用する手法については、特に限定しない。The
図8に示す干渉波302Aおよび干渉波302Bの影響を受けているとき、正しくグループ分けが行われると、周波数チャネルは、干渉波302Aのグループである周波数チャネルインデックス#1~#7のグループ#1と、干渉波302Bのグループである周波数チャネルインデックス#11~#15のグループ#2と、干渉なしのグループである周波数チャネルインデックス#0,#8~#10のグループ#3とに分けられる。干渉クラスタリング部215は、グループ分けの結果を周波数チャネル群決定部216に出力する。8, if grouping is performed correctly, the frequency channels are divided into
周波数チャネル群決定部216は、干渉クラスタリング部215によるグループ分けの結果から周波数チャネル群を決定する。具体的には、周波数チャネル群決定部216は、干渉を受けているグループと干渉を受けていないグループとを、周波数チャネルが隣接するように1つのグループにまとめる。グループのまとめ方のパターンが複数存在する場合、周波数チャネル群決定部216は、周波数チャネル数が最大になるパターンを選択する。具体的には、周波数チャネル群決定部216は、干渉クラスタリング部215によりグループ分けされた各グループのSINR推定値の平均値を算出し、算出した平均値に基づいて、各グループの干渉の有無を判定する。例えば、周波数チャネル群決定部216は、対象のグループのSINR推定値の平均値が所定の閾値σ以下であれば、「干渉なし」のグループであると判断し、SINR推定値の平均値が閾値σを超える場合には、「干渉あり」のグループであると判断する。図8に示す例では、グループ#1、グループ#2は「干渉あり」のグループであると判断され、グループ#3は「干渉なし」のグループであると判断される。The frequency channel group determination unit 216 determines a frequency channel group based on the grouping result by the
次に、周波数チャネル群決定部216は、「干渉あり」のグループと、「干渉なし」のグループに属する周波数チャネルとを、周波数チャネルが隣接するように1つのグループにまとめたときに、周波数チャネル数が最大となるようなグループ化を行う。なお、″「干渉あり」のグループと、「干渉なし」のグループに属する周波数チャネルとをまとめる″と表現する場合、まとめた後のグループには、「干渉あり」のグループに属する全ての周波数チャネルと、「干渉なし」のグループに属する一部または全部の周波数チャネルとが含まれることを意味する。図8の例では、干渉クラスタリング部215により生成されたグループは、周波数チャネルインデックス#1~#7の「干渉あり」のグループ#1、周波数チャネルインデックス#11~#15の「干渉あり」のグループ#2、周波数チャネルインデックス#0,#8~#10の「干渉なし」のグループ#3である。この場合、まず、「干渉あり」のグループ#1と、「干渉なし」のグループ#3に属する周波数チャネルとを、周波数チャネルが隣接するようにグループ化する。この場合、周波数チャネルインデックス#0~#10を1つのグループにまとめることができ、周波数チャネル数は11である。次に、「干渉あり」のグループ#2と、「干渉なし」のグループ#3に属する周波数チャネルとを、周波数チャネルが隣接するようにグループ化する。この場合、周波数チャネルインデックス#8~#15を1つのグループにまとめることができ、周波数チャネル数は8である。図8に示す例では、グループ#1と干渉なしの周波数チャネルとをまとめた場合の方が、グループ#2と干渉なしの周波数チャネルとをまとめた場合よりも周波数チャネル数が多くなる。この場合、周波数チャネル群決定部216は、周波数チャネルインデックス#0~#10を周波数チャネル群#1とし、周波数チャネルインデックス#11~#15を周波数チャネル群#2として周波数チャネル群を決定し、決定した周波数チャネル群を示す情報を干渉抑圧処理部206に出力する。Next, the frequency channel group determination unit 216 performs grouping such that the number of frequency channels is maximized when the frequency channels belonging to the "interference" group and the "no interference" group are grouped together so that the frequency channels are adjacent to each other. When it is expressed as "grouping the "interference" group and the "no interference" group together", it means that the group after grouping includes all the frequency channels belonging to the "interference" group and some or all the frequency channels belonging to the "no interference" group. In the example of FIG. 8, the groups generated by the
なお、周波数チャネルが「隣接するように」グループをまとめるとは、まとめた後のグループに属する周波数チャネルの周波数チャネルインデックスが連続数となることを指す。ここで、前提として周波数チャネルインデックスが周波数の大きさの順に付与されているものとする。例えば、上記のグループ#2とグループ#3に属する周波数チャネルとをまとめる際、周波数チャネルインデックス#0をまとめると、周波数チャネルインデックスが連続数とならないため、周波数チャネルインデックス#0は除外して、周波数チャネルインデックス#8~#10および#11~#15をまとめて、周波数チャネルインデックス#8~#15のグループにまとめている。
Note that grouping so that the frequency channels are "adjacent" means that the frequency channel indexes of the frequency channels belonging to the group after grouping are consecutive. Here, it is assumed that the frequency channel indexes are assigned in order of frequency magnitude. For example, when grouping the frequency channels belonging to the
次に、干渉抑圧処理部206の詳細について説明する。Next, we will explain the details of the interference
干渉抑圧処理部206は、周波数チャネル群生成部205が生成した周波数チャネル群単位で、干渉抑圧ウェイトを生成して、干渉抑圧処理を行う。具体的には、干渉抑圧処理部206は、生成された複数の周波数チャネル群のそれぞれに対して、干渉抑圧ウェイトを生成して、干渉抑圧処理を行う。干渉抑圧処理部206は、対象の周波数チャネル群に属する受信信号を受信信号蓄積メモリ部203から読み出し、対象の周波数チャネル群に属する受信信号を用いて、干渉抑圧処理を行う。干渉抑圧処理は、様々な方法が考えられ、干渉抑圧処理部206が用いる方法は、特定の方法に限定されない。例えば、受信装置20が複数の受信アンテナを備える場合、干渉波は受信アンテナ間で空間相関を有することから、干渉抑圧処理部206は、受信アンテナ毎の受信信号に含まれる干渉波を無相関にするような重みベクトルを乗算することにより、干渉波を抑圧することができる。まず、干渉抑圧処理部206は、干渉抑圧の対象となる受信信号r(n)を、周波数領域の信号に変換する。ここで受信アンテナ番号mおよびサブキャリア番号kに対する周波数領域表現の受信信号をRm(k)とする。受信アンテナ番号mは0以上アンテナ本数M未満の整数である。サブキャリア番号kは、周波数チャネル群のそれぞれに付された番号であり、0以上サブキャリア数Nb未満の整数である。干渉抑圧処理部206は、以下の数式(5)を用いて、受信信号Rm(k)に対して干渉抑圧処理を行うことができる。
The interference
なお、Qkは、サブキャリア番号kに対する干渉抑圧後の周波数領域表現の受信信号ベクトルであり、以下の数式(6)に示すようにM行1列の行列で表される。数式(6)のQm(k)は、受信アンテナ番号mおよびサブキャリア番号kに対する干渉抑圧後の周波数領域表現の受信信号である。 Q k is a received signal vector in the frequency domain after interference suppression for subcarrier number k, and is expressed as a matrix with M rows and 1 column as shown in the following formula (6): Q m (k) in formula (6) is a received signal in the frequency domain after interference suppression for receiving antenna number m and subcarrier number k.
また、Rkは、サブキャリア番号kに対する干渉抑圧前の周波数領域表現の受信信号ベクトルであり、以下の数式(7)に示すように、M行1列の行列で表される。 Furthermore, R k is a received signal vector in the frequency domain before interference suppression for subcarrier number k, and is expressed as a matrix with M rows and 1 column, as shown in the following equation (7).
Wkは、サブキャリア番号kに対する重みベクトルであり、以下の数式(8)に示すように、M行M列の行列で表される。Wm1,m2(k)は受信アンテナ番号m1およびサブキャリア番号kで受信した通信波に介在する干渉を抑圧するための受信アンテナm2に対応する重み信号である。m1およびm2のそれぞれは、0以上M未満の整数である。また、WH kはWkの複素共役転置を表す。 W k is a weight vector for subcarrier number k, and is expressed as a matrix with M rows and M columns as shown in the following formula (8). W m1,m2 (k) is a weight signal corresponding to receiving antenna number m1 and receiving antenna m2 for suppressing interference present in a communication wave received at subcarrier number k. Each of m1 and m2 is an integer equal to or greater than 0 and less than M. Furthermore, W H k represents the complex conjugate transpose of W k .
干渉抑圧処理部206は、以下の数式(9)を満たすように重みベクトルWkを生成する。重みベクトルWkは、干渉抑圧ウェイトの一例である。
The interference
ここで、Gkは、サブキャリア番号kに対する干渉波の周波数領域信号Bm(k)の共分散行列であり、IはM行M列の単位行列である。Bm(k)は、受信アンテナmおよびサブキャリア番号kの干渉波の周波数領域信号である。数式(9)を満たす重みベクトルWkを用いて、数式(5)~(8)により干渉抑圧を行うことで、受信アンテナ間で無相関な干渉波に変換することができる。 Here, G k is the covariance matrix of the frequency domain signal B m (k) of the interference wave for subcarrier number k, and I is a unit matrix with M rows and M columns. B m (k) is the frequency domain signal of the interference wave of receiving antenna m and subcarrier number k. By using the weight vector W k that satisfies equation (9) and performing interference suppression using equations (5) to (8), it is possible to convert to an interference wave that is uncorrelated between receiving antennas.
続いて、周波数チャネル結合部207の詳細について説明する。
Next, we will explain the details of the frequency
周波数チャネル結合部207は、干渉抑圧処理された受信信号に対して、周波数チャネル群単位で周波数チャネルの結合処理を行う。ここでは周波数チャネル結合部207は、複数の周波数チャネル群のそれぞれについて、周波数チャネル群毎に結合処理を行う。周波数チャネル結合部207は、例えば、干渉抑圧処理された周波数領域の信号である受信信号に対して、既知信号による逆変調処理を行うことで、チャネル周波数応答を周波数チャネル毎に推定する。続いて周波数チャネル結合部207は、対象の周波数チャネル群に属する周波数チャネルのチャネル周波数応答を結合する。このとき、チャネル周波数応答の初期位相は、周波数チャネルごとに異なるため、周波数チャネル結合部207は、周波数チャネルごとに初期位相を合わせる必要がある。周波数チャネルごとに初期位相を合わせる方法としては、例えば、非特許文献1に記載されている、送信装置10および受信装置20以外の無線機を用意して初期位相を推定する方法を使用してもよい。また、上記では、同じ周波数チャネル群に属する周波数チャネルは、周波数チャネルインデックスが連続するように、周波数チャネル群を生成しているが、周波数チャネルインデックスが連続しない場合、周波数チャネル結合部207は、チャネル周波数応答の補間処理を行ってもよい。例えば、ある周波数チャネル群の周波数チャネルインデックスが#1,#3,#4である場合、周波数チャネルインデックス#2のチャネル周波数応答がないため、周波数チャネル結合部207は、周波数チャネルインデックス#1,#3のチャネル周波数応答を用いて、周波数チャネルインデックス#2のチャネル周波数応答を補間する。チャネル周波数応答を結合することで、周波数チャネル群の周波数帯域幅に相当するチャネル周波数応答を得ることができる。周波数チャネル結合部207は、結合したチャネル周波数応答にIDFT(Inverse Division Fourier Transform)処理を適用し、時間領域信号に変換することで、周波数チャネル群のチャネルインパルス応答を得ることができる。
The frequency
図9は、チャネルインパルス応答の一例を示す図である。図9の横軸は時間(サンプル)であり、縦軸はチャネルインパルス応答の電力値である。ここで1サンプルの時間は、周波数チャネル群で結合したときの周波数帯域幅[Hz]の逆数に相当する。このため、結合したときの周波数帯域幅が大きければ、チャネルインパルス応答の1サンプル当たりの時間が短くなる。つまり、チャネルインパルス応答の時間分解能が高くなる。周波数チャネル結合部207は、周波数チャネル群毎にチャネルインパルス応答を求めて出力する。
Figure 9 is a diagram showing an example of a channel impulse response. The horizontal axis of Figure 9 is time (samples), and the vertical axis is the power value of the channel impulse response. Here, the time for one sample corresponds to the inverse of the frequency bandwidth [Hz] when combined in a frequency channel group. Therefore, if the frequency bandwidth when combined is large, the time per sample of the channel impulse response becomes shorter. In other words, the time resolution of the channel impulse response becomes higher. The frequency
ここで直接波成分抽出部208の詳細について説明する。
Here we will explain the details of the direct wave
直接波成分抽出部208は、周波数チャネル結合部207で求めた周波数チャネル群毎のチャネルインパルス応答から直接波成分を抽出する。例えば、直接波成分抽出部208は、図9に示すチャネルインパルス応答に対して予め定められた閾値αを基準として、閾値αを超えるチャネルインパルス応答のサンプルから最も小さい時間サンプルインデックスを直接波成分とすることができる。直接波成分抽出部208は、決定した直接波成分となる時間インデックス「0」のチャネルインパルス応答のみを抽出する。直接波成分抽出部208は、上述の処理を周波数チャネル群毎に行い、周波数チャネル群の数分の直接波成分を出力することになる。The direct wave
続いて電波到来角推定部209の詳細について説明する。
Next, we will explain the details of the radio wave arrival
電波到来角推定部209は、直接波成分抽出部208で求めた周波数チャネル群毎の直接波成分を用いて、直接波成分の到来角を電波到来角として推定する。例えば、電波到来角推定部209は、以下の数式(10)を用いて、電波到来角βj,a,bを推定することができる。
The radio wave arrival
なお、電波到来角βj,a,bは、周波数チャネル群インデックスjにおける受信アンテナの組のインデックスa,b間のチャネルインパルス応答から求めた電波到来角推定値である。なお、a,bは、干渉抑圧に用いた受信アンテナの組のインデックスであり、aはbと異なる値である。周波数チャネル群インデックスjは、0以上周波数チャネル群の数Ng未満の整数である。ρjは周波数チャネル群インデックスjにおける結合後の周波数帯域幅の逆数に相当する波長であり、γj,a,bは周波数チャネル群インデックスjにおける受信アンテナの組のインデックスa,b間のチャネルインパルス応答に対する位相差であり、Da,bは受信アンテナの組aとbとの間隔である。なお、受信アンテナの組aとbとの間隔とは、受信アンテナの組aが受信アンテナ「1」,「2」の組であって、受信アンテナの組bが受信アンテナ「2」,「3」の組である場合、受信アンテナ「1」,「2」の中間点と、受信アンテナ「2」,「3」の中間点との間の距離を意味する。 The radio wave arrival angle β j,a,b is an estimated value of the radio wave arrival angle obtained from the channel impulse response between indexes a and b of the pair of receiving antennas in frequency channel group index j. Note that a and b are indexes of the pair of receiving antennas used for interference suppression, and a is a value different from b. The frequency channel group index j is an integer equal to or greater than 0 and less than the number N g of frequency channel groups. ρ j is a wavelength corresponding to the inverse of the frequency bandwidth after combination in frequency channel group index j, γ j,a,b is a phase difference with respect to the channel impulse response between indexes a and b of the pair of receiving antennas in frequency channel group index j, and D a,b is the interval between the pair of receiving antennas a and b. Note that the interval between the pair of receiving antennas a and b means the distance between the midpoint of the receiving antennas "1" and "2" and the midpoint of the receiving antennas "2" and "3" when the pair of receiving antennas a is the pair of receiving antennas "1" and "2" and the pair of receiving antennas b is the pair of receiving antennas "2" and "3".
なお、既知信号のシンボルインデックスをp、既知信号のシンボル数Np、周波数チャネル群インデックスjであって受信アンテナの組のインデックスaのシンボルインデックスpにおけるチャネルインパルス応答をSj,a(p)とした場合、位相差γj,a,bは、以下の数式(11)で求めることができる。電波到来角推定部209は、周波数チャネル群毎の電波到来角を推定し、推定した電波到来角を出力する。
If the symbol index of the known signal is p, the number of symbols of the known signal is N p , the frequency channel group index is j, and the channel impulse response at symbol index p of index a of the receiving antenna set is S j,a (p), the phase difference γ j,a,b can be obtained by the following formula (11). The radio wave arrival
次に、到来角合成部210の詳細について説明する。Next, we will explain the details of the arrival angle synthesis unit 210.
到来角合成部210は、電波到来角推定部209で求めた周波数チャネル群毎の電波到来角の推定結果を合成する。到来角合成部210は、干渉測定部204で算出したSINR推定値を用いて、電波到来角の推定結果を合成することができる。SINR推定値は、周波数チャネル群ごとに推定されているため、到来角合成部210は、SINR推定値に従って、電波到来角推定値を重み付けして合成することができる。例えば、xから求められる合成ウェイト関数G(x)、周波数チャネル群インデックスjのSINR推定値の平均値SINRj、周波数チャネル群インデックスjの電波到来角推定値の平均βjとすると、合成後の電波到来角βは、以下の数式(12)で表される。なお、周波数チャネル群インデックスjのSINR推定値の平均値SINRjは数式(13)を用いて求めることができ、周波数チャネル群インデックスjの電波到来角推定値の平均βjは数式(14)を用いて求めることができる。
The arrival angle synthesis unit 210 synthesizes the estimation results of the radio wave arrival angle for each frequency channel group obtained by the radio wave arrival
なお、Fjは周波数チャネル群インデックスjにおける周波数チャネルインデックスの集合、Yは干渉抑圧処理の受信アンテナの組に含まれる受信アンテナインデックスの集合、Z(a)は干渉抑圧処理における受信アンテナの組に含まれる受信アンテナインデックスaの対となる受信アンテナインデックスの集合、Nhは干渉抑圧処理の受信アンテナの組の数を表す。 Here, Fj represents a set of frequency channel indexes in frequency channel group index j, Y represents a set of receiving antenna indexes included in the receiving antenna pair in the interference suppression process, Z(a) represents a set of receiving antenna indexes that are pairs with receiving antenna index a included in the receiving antenna pair in the interference suppression process, and Nh represents the number of receiving antenna pairs in the interference suppression process.
また、合成ウェイト関数G(x)は、例えば以下の数式(15)で表される。 Furthermore, the composite weight function G(x) is expressed, for example, by the following equation (15).
図10は、実施の形態1にかかる受信装置20の動作を説明するためのフローチャートである。なお、各処理の詳細については上述しているので、ここでは処理の流れについて主に説明する。
Figure 10 is a flowchart for explaining the operation of the receiving
受信装置20の受信アンテナ部200が送信装置10から信号を受信すると、周波数変換部201は、周波数変換処理を行う(ステップS101)。受信フィルタ部202は、周波数変換後の信号に対して受信フィルタリング処理を行う(ステップS102)。受信フィルタリング処理後の受信信号は、受信信号蓄積メモリ部203に蓄積されると共に、干渉測定部204へ出力される。干渉測定部204は、受信信号に含まれる干渉の状態を測定して干渉情報を生成する(ステップS103)。干渉情報は、周波数チャネル群生成部205および到来角合成部210のそれぞれへ出力される。When the receiving
周波数チャネル群生成部205は、周波数チャネル群を生成し(ステップS104)、生成した周波数チャネル群を干渉抑圧処理部206に通知する。干渉抑圧処理部206は、周波数チャネル群毎に干渉抑圧処理を行い(ステップS105)、干渉抑圧処理後の信号を周波数チャネル結合部207へ出力する。The frequency channel
周波数チャネル結合部207は、周波数チャネル群毎に、周波数チャネルを結合する(ステップS106)。直接波成分抽出部208は、周波数チャネル群毎に直接波成分を抽出する(ステップS107)。The frequency
電波到来角推定部209は、周波数チャネル群毎に電波の到来角を推定し(ステップS108)、周波数チャネル群の数分の電波到来角推定値を出力する。到来角合成部210は、電波到来角推定部209が出力する周波数チャネル群毎の電波到来角推定値を、干渉測定部204が生成する干渉情報に基づいて合成し(ステップS109)、到来角推定結果を出力する。The radio wave arrival
以上説明したように、実施の形態1にかかる電波到来角推定装置1は、複数の周波数チャネルを使用して送信された信号を受信した受信信号から電波の到来方向を推定する電波到来角推定装置であって、周波数チャネル毎に干渉の状態を測定して干渉情報を生成する干渉測定部204と、周波数チャネル毎に生成された干渉情報に基づいて、複数の周波数チャネルをグループ分けして周波数チャネル群を生成する周波数チャネル群生成部205と、周波数チャネル群単位で干渉抑圧ウェイトを生成して干渉抑圧処理を行う干渉抑圧処理部206と、干渉抑圧処理を行った後の信号に対して、周波数チャネル群単位で周波数チャネルを結合する周波数チャネル結合部207と、周波数チャネルを結合した後の信号に対して、直接波成分を抽出する直接波成分抽出部208と、直接波成分から電波の到来角を推定する電波到来角推定部209と、を備える。このような構成を有することにより、電波到来角推定装置1は、周波数チャネル毎の干渉情報に基づいて、周波数チャネルをグループ分けした周波数チャネル群単位で、干渉抑圧処理および周波数チャネルの結合が行われる。したがって、複数の干渉源が存在した場合など周波数チャネル毎に干渉の状態が異なる場合であっても、干渉の状態に合わせて周波数チャネル群単位で干渉抑圧ウェイトが生成され、生成された干渉抑圧ウェイトを用いて干渉抑圧処理が行われる。そして、干渉抑圧処理後の信号に対して、周波数チャネル群単位で周波数チャネルが結合される。したがって、周波数チャネルを結合する際には、共通の干渉抑圧ウェイトを用いて干渉抑圧処理が行われた周波数チャネル同士が結合されるため、干渉抑圧処理を行った後の信号であっても、周波数チャネルの結合を行うことができる。したがって、干渉抑圧処理によって他システムからの干渉の影響を低減し、さらに、周波数チャネルの結合によりマルチパスの影響を低減することができるため、マルチパスおよび他システムからの干渉の両方の影響を受ける環境であっても、電波の到来角を高精度に推定することが可能になる。As described above, the radio wave arrival
実施の形態1にかかる電波到来角推定装置1の具体的な構成としては、干渉抑圧処理部206は、複数の周波数チャネル群のそれぞれに対して干渉抑圧ウェイトを生成して干渉抑圧処理を行い、周波数チャネル結合部207は、周波数チャネル群毎に周波数チャネルを結合し、直接波成分抽出部208は、周波数チャネル群毎に直接波成分を抽出し、電波到来角推定部209は、周波数チャネル群毎に到来角を推定し、周波数チャネル群毎の複数の到来角の推定結果を合成する到来角合成部210、をさらに備える。これにより、複数の周波数チャネル群のそれぞれについて、周波数チャネル群単位で電波の到来角が推定され、複数の到来角の推定結果を合成することにより、最終的な推定結果を得られる。したがって、精度高く電波の到来角を推定することができる。
As a specific configuration of the radio wave arrival
ここで、干渉測定部204は、信号対干渉雑音電力比、干渉電力、干渉到来角および干渉複素信号のうちいずれか1つ以上の推定値を含む干渉情報を生成することができる。周波数チャネル群生成部205は、信号対干渉雑音電力比、干渉電力、干渉到来角および干渉複素信号のうちいずれか1つ以上を使用して、周波数チャネル群を生成することができる。Here, the
周波数チャネル群生成部205は、干渉情報に基づいて、干渉の有無と干渉源とに基づいて周波数チャネルをグループ分けする干渉クラスタリング部215と、グループ分けの結果から周波数チャネル群を決定する周波数チャネル群決定部216と、を有する。より具体的には、干渉クラスタリング部215は、周波数チャネルを、干渉を受けている周波数チャネルのグループと、干渉を受けていない周波数チャネルのグループとに分け、さらに、干渉を受けている周波数チャネルのグループを干渉源毎にグループ分けし、周波数チャネル群決定部216は、干渉を受けているグループと干渉を受けていないグループに属する周波数チャネルとを周波数チャネルが隣接するように1つのグループにまとめたときに周波数チャネル数が最大になるように周波数チャネル群を決定する。これにより、異なる干渉源からの干渉を受けている周波数チャネルを異なる周波数チャネル群にグループ分けすることができると共に、周波数チャネル群単位での周波数チャネル結合を行いやすくなる。The frequency channel
また、実施の形態1によれば、複数の周波数チャネルを使用して送信された信号を受信した受信信号から電波の到来方向を推定する電波到来角推定装置1が実行する電波到来角推定方法を実現することができる。この電波到来角推定方法は、周波数チャネル毎に干渉の状態を測定して干渉情報を生成するステップと、周波数チャネル毎に生成された干渉情報に基づいて、複数の周波数チャネルをグループ分けして周波数チャネル群を生成するステップと、周波数チャネル群単位で干渉抑圧ウェイトを生成して干渉抑圧処理を行うステップと、干渉抑圧処理を行った後の信号に対して、直接波成分を抽出するステップと、直接波成分から電波の到来角を推定するステップと、を含む。
According to the first embodiment, a radio wave arrival angle estimation method can be realized, which is executed by the radio wave arrival
なお、図示していないが、到来角の推定結果は、位置推定のために使用することができる。近年、位置情報を活用した多様なサービスが提供されており、位置情報の需要は高まっている。屋外の位置推定については、GPS(Global Positioning System)を用いた方式が広く普及している。屋内やGPSを利用できない場所などでは、無線LAN(Local Area Network)、Bluetooth、UWB(Ultra Wide Band)などの無線信号を利用した位置推定技術が提案されている。無線信号を利用した位置推定技術には、電波の到来角を利用したAoAおよびAoD、電波の到着時間を利用したToA(Time of Arrival)、TDoA(Time Difference of Arrival)、受信信号の強度を利用したRSSI(Received Signal Strength Indication)などがある。例えばBluetoothでは、バージョン5.1よりも前まではRSSIによる位置推定が採用されていたが、バージョン5.1以降からは電波の到来角を利用したAoA、AoDによる位置推定が採用されており、RSSIよりも高い位置推定精度が期待されている。Bluetoothによる位置推定は、屋内で使用されることが多いと考えられるが、屋内では壁面などで電波が反射して到来するため、受信装置において、直接波と反射波とが到来するマルチパスの影響を受けやすい。マルチパスのある環境では位置推定精度が低下してしまう。また、Bluetoothや無線LANはいずれもアンライセンスバンドを使用した無線通信規格であるため、他システムからの干渉の影響を受ける可能性が高い。このため、屋内の位置推定のために利用される場合、上述のように、複数の干渉源が存在することを想定した、干渉対策とマルチパス対策とを両立することのできる技術は好適である。Although not shown, the estimated result of the angle of arrival can be used for position estimation. In recent years, various services utilizing location information have been provided, and the demand for location information is increasing. For outdoor position estimation, a method using the Global Positioning System (GPS) is widely used. For indoor use or places where GPS cannot be used, position estimation techniques using wireless signals such as wireless LAN (Local Area Network), Bluetooth, and UWB (Ultra Wide Band) have been proposed. Position estimation techniques using wireless signals include AoA and AoD using the angle of arrival of radio waves, ToA (Time of Arrival) and TDoA (Time Difference of Arrival) using the arrival time of radio waves, and RSSI (Received Signal Strength Indication) using the strength of the received signal. For example, in Bluetooth, position estimation using RSSI was adopted until version 5.1, but from version 5.1 onwards, position estimation using AoA and AoD using the angle of arrival of radio waves has been adopted, and higher position estimation accuracy than RSSI is expected. Bluetooth location estimation is likely to be used indoors, but indoors, radio waves arrive after being reflected by walls and the like, so the receiving device is susceptible to multipath effects caused by direct waves and reflected waves. In a multipath environment, the accuracy of location estimation decreases. In addition, since Bluetooth and wireless LAN are both wireless communication standards that use unlicensed bands, they are likely to be affected by interference from other systems. For this reason, when used for indoor location estimation, a technology that can achieve both interference countermeasures and multipath countermeasures assuming the presence of multiple interference sources, as described above, is preferable.
実施の形態2.
図11は、実施の形態2にかかる受信装置20Aの構成を示す図である。なお、図示しないが、実施の形態2にかかる電波到来角推定装置1Aは、実施の形態1の受信装置20の代わりに受信装置20Aを備える。
11 is a diagram showing a configuration of a receiving
受信装置20Aは、受信アンテナ部200と、周波数変換部201と、受信フィルタ部202と、受信信号蓄積メモリ部203と、干渉測定部204と、周波数チャネル群生成部205と、周波数チャネル群選択部217と、干渉抑圧処理部206Aと、周波数チャネル結合部207Aと、直接波成分抽出部208Aと、電波到来角推定部209Aとを有する。つまり、受信装置20Aは、受信装置20の構成に周波数チャネル群選択部217を加え、到来角合成部210を省略し、干渉抑圧処理部206、周波数チャネル結合部207、直接波成分抽出部208および電波到来角推定部209の代わりに、干渉抑圧処理部206A、周波数チャネル結合部207A、直接波成分抽出部208Aおよび電波到来角推定部209Aを備える。以下、実施の形態1と共通の機能については詳細な説明を省略し、実施の形態1と異なる部分について主に説明する。The receiving
実施の形態1では、受信装置20は、周波数チャネル群を生成すると、全ての周波数チャネル群のそれぞれについて、周波数チャネル群毎に、複数の電波到来角推定値を算出し、複数の電波到来角推定値を合成することで到来角推定結果を生成した。この場合、周波数チャネル群の数が多くなると、処理量が多くなる。このため、実施の形態2では、周波数チャネル群の中から先に1つの周波数チャネル群を選択し、選択された周波数チャネル群についてのみ処理を行う。In the first embodiment, when the receiving
周波数チャネル群選択部217は、周波数チャネル群生成部205が生成した周波数チャネル群の中から1つを選択する。例えば、周波数チャネル群選択部217は、周波数チャネル群に属する周波数チャネルの数に基づいて、周波数チャネル群を選択することができる。周波数チャネル群選択部217は、周波数チャネルの数が最も多い周波数チャネル群を選択する。周波数チャネル数が多いと、チャネルインパルス応答の時間分解能が高くなり、マルチパス対策として有効である。ただし、周波数チャネル群の選択方法は、上記の例に限定されない。
The frequency channel
干渉抑圧処理部206A、周波数チャネル結合部207A、直接波成分抽出部208Aおよび電波到来角推定部209Aの処理については、選択された周波数チャネル群についてのみ処理を行う点以外は、干渉抑圧処理部206、周波数チャネル結合部207、直接波成分抽出部208、電波到来角推定部209のそれぞれと同様である。この場合、電波到来角推定部209は、1つの電波到来角の推定値を算出することになるため、この推定値を到来角推定結果として出力し、到来角合成処理は不要となる。The processing of the interference
図12は、実施の形態2にかかる受信装置20Aの動作を説明するためのフローチャートである。ステップS101からステップS104までは、実施の形態1と同様であるため説明を省略する。
Figure 12 is a flowchart for explaining the operation of the receiving
周波数チャネル群生成部205により周波数チャネル群が生成されると、周波数チャネル群選択部217は、周波数チャネル群の1つを選択する(ステップS201)。干渉抑圧処理部206Aは、選択された周波数チャネル群の干渉抑圧処理を行う(ステップS202)。周波数チャネル結合部207Aは、選択された周波数チャネル群の周波数チャネル結合処理を行う(ステップS203)。When the frequency channel groups are generated by the frequency channel
直接波成分抽出部208Aは、選択された周波数チャネル群の直接波成分を抽出する(ステップS204)。電波到来角推定部209Aは、選択された周波数チャネル群の電波の到来角推定処理を行う(ステップS205)。The direct wave
以上説明したように、実施の形態2によれば、実施の形態1と同様の効果を奏することができる。また、実施の形態2によれば、受信装置20Aは、生成された周波数チャネル群の中から処理対象の周波数チャネル群を選択する周波数チャネル群選択部217、をさらに備え、干渉抑圧処理部206A、周波数チャネル結合部207A、直接波成分抽出部208A、および電波到来角推定部209Aのそれぞれは、選択された周波数チャネル群に対して処理を行い、選択されなかった周波数チャネル群に対しては処理を行わない。これにより、実施の形態1と比較して、到来角推定にかかる処理量を削減することができる。特に、周波数チャネル群生成部205が生成する周波数チャネル群の数が多くなった場合、処理量の削減効果は高くなる。As described above, according to the second embodiment, the same effect as that of the first embodiment can be achieved. Also, according to the second embodiment, the receiving
周波数チャネル群選択部217は、周波数チャネル群に属する周波数チャネルの数に基づいて周波数チャネル群を選択することができる。例えば、周波数チャネルの数が最も多い1つの周波数チャネル群を処理対象として選択することができる。周波数チャネルの数が多いほど、チャネルインパルス応答の時間分解能が高くなるため、マルチパス対策として有効である。The frequency channel
実施の形態3.
実施の形態3にかかる電波到来角推定装置1Bは、図2に示す送信装置10と、図11に示す受信装置20Aとを有する。ただし、送信装置10の周波数変換部103が使用する周波数ホッピングシーケンスが実施の形態1,2と異なる。実施の形態1,2では、図6に示したような、周波数ホッピングパターンがランダムな周波数ホッピング系列に従って既知信号が送信されていた。このため、実施の形態2にかかる受信装置20Aで周波数結合する際には、処理対象の周波数チャネル群に含まれる周波数チャネルを対象に結合処理を行うため、対象となる周波数チャネルが揃うまで周波数結合処理を行うことができない。このため、実施の形態3では、処理対象の周波数チャネル群に属する周波数チャネルが優先的に、先行して送信されるような周波数ホッピングシーケンスを用いることで、周波数結合処理を早く開始することができるようにする。
The radio wave arrival angle estimation device 1B according to the third embodiment includes the transmitting
図13は、実施の形態3で送信される既知信号を示す図である。既知信号は、周波数チャネルインデックス#0~#15の16個の周波数チャネルを用いて、通信スロット400として送信される。ここでは、周波数チャネルインデックス#0~#10の周波数チャネルからなる周波数チャネル群#1と、周波数チャネルインデックス#11~#15の周波数チャネルからなる周波数チャネル群#2とが生成されている。このとき、周波数チャネル群#1の方が周波数チャネル群#2よりも周波数チャネル数が多いことから、受信装置20Aの周波数チャネル群選択部217では、周波数チャネル群#1が選択される。この場合、周波数変換部103では、周波数チャネル群#1に属する周波数チャネルインデックス#0~#10の周波数チャネルの中で、ランダムに周波数ホッピングし、周波数チャネル群#1の周波数ホッピングが終わった後に、周波数チャネル群#2に属する周波数チャネルインデックス#11~#15の周波数チャネルの中でランダムに周波数ホッピングするような周波数ホッピングシーケンスに従って、周波数変換処理が行われる。周波数ホッピングシーケンスを示す系列を、周波数ホッピング系列と称する。例えば、周波数変換部103が使用する周波数ホッピング系列は、図13に示す例では、#0,#2,#4,#6,#8,#10,#1,#3,#5,#7,#9,#11,#13,#15,#12,#14である。これにより、周波数チャネル群#1に属する周波数チャネルが周波数チャネル群#2に属する周波数チャネルよりも先に送信され、図13の例では、周波数チャネルインデックス#9の周波数チャネルを受信した時点で、受信装置20Aは、周波数チャネル群#1の周波数結合処理を開始することができる。実施の形態3では、受信装置20Bにおいて、図6に示すような周波数チャネル群に基づかない周波数ホッピング系列を用いるよりも早く、周波数チャネル群#1に属する周波数チャネルが揃う。このため、実施の形態2よりも、周波数チャネル群#1の周波数結合処理を早く開始することができ、結果として、実施の形態2よりも早く到来角推定結果を得ることができる。
Figure 13 is a diagram showing a known signal transmitted in the third embodiment. The known signal is transmitted as a communication slot 400 using 16 frequency channels with frequency
以上説明したように、実施の形態3によれば、電波到来角推定装置1Bは、周波数ホッピング方式により複数の周波数チャネルの中から送信する周波数チャネルに合わせて送信する信号の周波数を変換する周波数変換部103、を備え、周波数変換部103は、受信装置20Aの周波数チャネル群生成部205が生成した周波数チャネル群に基づく切替パターンで、送信する周波数チャネルを切り替える。より具体的には、周波数変換部103は、周波数チャネル群選択部217が選択した周波数チャネル群に属する周波数チャネルを他の周波数チャネル群に属する周波数チャネルよりも先行して送信するように、送信する周波数チャネルを切り替える。これにより、実施の形態1,2と同様の効果を奏するとともに、実施の形態2よりも早く到来角推定結果を得ることができる。As described above, according to the third embodiment, the radio wave arrival angle estimation device 1B includes a frequency conversion unit 103 that converts the frequency of a signal to be transmitted from among a plurality of frequency channels by a frequency hopping method to match the frequency channel to be transmitted, and the frequency conversion unit 103 switches the frequency channel to be transmitted in a switching pattern based on the frequency channel group generated by the frequency channel
続いて、実施の形態1~3にかかる電波到来角推定装置1,1A,1Bのハードウェア構成について説明する。電波到来角推定装置1,1A,1Bの各部は、処理回路により実現される。これらの処理回路は、専用のハードウェアにより実現されてもよいし、CPU(Central Processing Unit)を用いた制御回路であってもよい。Next, the hardware configuration of the radio wave arrival
上記の処理回路が、専用のハードウェアにより実現される場合、これらは、図14に示す処理回路90により実現される。図14は、実施の形態1~3にかかる電波到来角推定装置1,1A,1Bの機能を実現するための専用のハードウェアを示す図である。処理回路90は、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ASIC(Application Specific Integrated Circuit)、FPGA(Field Programmable Gate Array)、またはこれらを組み合わせたものである。When the above processing circuits are realized by dedicated hardware, they are realized by the processing circuit 90 shown in FIG. 14. FIG. 14 is a diagram showing dedicated hardware for realizing the functions of the radio wave arrival
上記の処理回路が、CPUを用いた制御回路で実現される場合、この制御回路は例えば図15に示す構成の制御回路91である。図15は、実施の形態1~3にかかる電波到来角推定装置1,1A,1Bの機能を実現するための制御回路91の構成を示す図である。図15に示すように、制御回路91は、プロセッサ92と、メモリ93とを備える。プロセッサ92は、CPUであり、演算装置、マイクロプロセッサ、マイクロコンピュータ、DSP(Digital Signal Processor)などとも呼ばれる。メモリ93は、例えば、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ、EPROM(Erasable Programmable ROM)、EEPROM(登録商標)(Electrically EPROM)などの不揮発性または揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、DVD(Digital Versatile Disk)などである。
When the above processing circuit is realized by a control circuit using a CPU, the control circuit is, for example, a
上記の処理回路が制御回路91により実現される場合、プロセッサ92がメモリ93に記憶された、各構成要素の処理に対応するプログラムを読み出して実行することにより実現される。プロセッサ92が実行するプログラムは、記憶媒体に記憶された状態で提供されてもよいし、インターネットのような通信路を介して提供されてもよい。また、メモリ93は、プロセッサ92が実行する各処理における一時メモリとしても使用される。なお、実施の形態1~3にかかる電波到来角推定装置1,1A,1Bの機能は、専用のハードウェアである処理回路90とCPUを用いた制御回路91とを組み合わせて実現されてもよい。また、上記の実施の形態で図示した機能部の分け方は一例であり、1つの機能部で示した機能を複数の処理回路で実現してもよいし、複数の機能部で示した機能を複数の処理回路で実現してもよい。When the above processing circuit is realized by the
上記のように、実施の形態1~3によれば、複数の周波数チャネルを使用して送信された信号を受信した受信信号から電波の到来方向を推定する電波到来角推定装置1,1A,1Bを制御する制御回路91を実現することができる。この制御回路91は、周波数チャネル毎に干渉の状態を測定して干渉情報を生成するステップと、周波数チャネル毎に生成された干渉情報に基づいて、複数の周波数チャネルをグループ分けして周波数チャネル群を生成するステップと、周波数チャネル群単位で干渉抑圧ウェイトを生成して干渉抑圧処理を行うステップと、干渉抑圧処理を行った後の信号に対して、直接波成分を抽出するステップと、直接波成分から電波の到来角を推定するステップと、を電波到来角推定装置1,1A,1Bに実行させる。As described above, according to the first to third embodiments, a
また、実施の形態1~3によれば、複数の周波数チャネルを使用して送信された信号を受信した受信信号から電波の到来方向を推定する電波到来角推定装置1,1A,1Bを制御するためのプログラムを記憶した記憶媒体を実現することができる。この記憶媒体が記憶するプログラムは、周波数チャネル毎に干渉の状態を測定して干渉情報を生成するステップと、周波数チャネル毎に生成された干渉情報に基づいて、複数の周波数チャネルをグループ分けして周波数チャネル群を生成するステップと、周波数チャネル群単位で干渉抑圧ウェイトを生成して干渉抑圧処理を行うステップと、干渉抑圧処理を行った後の信号に対して、直接波成分を抽出するステップと、直接波成分から電波の到来角を推定するステップと、を電波到来角推定装置1,1A,1Bに実行させる。
Furthermore, according to the first to third embodiments, a storage medium can be realized that stores a program for controlling the radio wave arrival
以上の実施の形態に示した構成は、一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、実施の形態同士を組み合わせることも可能であるし、要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。 The configurations shown in the above embodiments are merely examples, and may be combined with other known technologies, or the embodiments may be combined with each other. Also, parts of the configurations may be omitted or modified without departing from the spirit of the invention.
例えば、実施の形態2では周波数チャネル群選択部217が1つの周波数チャネル群を選択するものとしたが、生成した周波数チャネル群の数よりも少ない数であれば、複数の周波数チャネル群を選択してもよい。例えば、3つ以上の周波数チャネル群が生成された場合において、各周波数チャネル群に属する周波数チャネルの数が多い方から2つの周波数チャネル群を選択してもよい。この場合、到来角推定装置1Bは、実施の形態1にかかる受信装置20を備え、実施の形態1と同様に、選択された複数の周波数チャネル群のそれぞれについて、干渉抑圧処理部206、周波数チャネル結合部207、直接波成分抽出部208および電波到来角推定部209の処理が行われ、到来角合成部210において、合成処理が行われる。For example, in the second embodiment, the frequency channel
なお、上記の各実施の形態で説明した機能の全てを備えなくても、電波到来角推定装置1,1A,1Bの一部の機能のみを実現することもできる。この場合、例えば、受信装置20,20Aの機能のみを有する装置を「電波到来角推定装置」と称することもある。It is also possible to realize only some of the functions of the radio wave arrival
1,1A,1B 電波到来角推定装置、10 送信装置、20,20A 受信装置、90 処理回路、91 制御回路、92 プロセッサ、93 メモリ、100 既知信号生成部、101 変調部、102 送信フィルタ部、103,201 周波数変換部、104 送信アンテナ部、200 受信アンテナ部、202 受信フィルタ部、203 受信信号蓄積メモリ部、204 干渉測定部、205 周波数チャネル群生成部、206,206A 干渉抑圧処理部、207,207A 周波数チャネル結合部、208,208A 直接波成分抽出部、209,209A 電波到来角推定部、210 到来角合成部、211 SINR推定部、212 干渉電力推定部、213 干渉到来角推定部、214 干渉複素信号推定部、215 干渉クラスタリング部、216 周波数チャネル群決定部、217 周波数チャネル群選択部、300,400 通信スロット、301 信号、302A,302B 干渉波。1, 1A, 1B Radio wave arrival angle estimation device, 10 Transmitting device, 20, 20A Receiving device, 90 Processing circuit, 91 Control circuit, 92 Processor, 93 Memory, 100 Known signal generation unit, 101 Modulation unit, 102 Transmission filter unit, 103, 201 Frequency conversion unit, 104 Transmitting antenna unit, 200 Receiving antenna unit, 202 Receiving filter unit, 203 Received signal storage memory unit, 204 Interference measurement unit, 205 Frequency channel group generation unit, 206, 206A Interference suppression processing unit, 207, 207A Frequency channel combination unit, 208, 208A Direct wave component extraction unit, 209, 209A Radio wave arrival angle estimation unit, 210 Arrival angle synthesis unit, 211 SINR estimation unit, 212 Interference power estimation unit, 213 Interference arrival angle estimation unit, 214 Interference complex signal estimation unit, 215 interference clustering unit, 216 frequency channel group determination unit, 217 frequency channel group selection unit, 300, 400 communication slots, 301 signal, 302A, 302B interference waves.
Claims (14)
前記周波数チャネル毎に干渉の状態を測定して干渉情報を生成する干渉測定部と、
前記周波数チャネル毎に生成された前記干渉情報に基づいて、複数の前記周波数チャネルをグループ分けして周波数チャネル群を生成する周波数チャネル群生成部と、
前記周波数チャネル群単位で干渉抑圧ウェイトを生成して干渉抑圧処理を行う干渉抑圧処理部と、
前記干渉抑圧処理を行った後の信号に対して、前記周波数チャネル群単位で前記周波数チャネルを結合する周波数チャネル結合部と、
前記周波数チャネルを結合した後の信号に対して、直接波成分を抽出する直接波成分抽出部と、
前記直接波成分から電波の到来角を推定する電波到来角推定部と、
を備えることを特徴とする電波到来角推定装置。 A radio wave arrival angle estimation device that estimates an arrival direction of a radio wave from a received signal that is received using a plurality of frequency channels,
an interference measurement unit that measures an interference state for each of the frequency channels and generates interference information;
a frequency channel group generation unit that generates a frequency channel group by grouping the plurality of frequency channels based on the interference information generated for each of the frequency channels;
an interference suppression processing unit that generates an interference suppression weight for each frequency channel group and performs interference suppression processing;
a frequency channel combining unit that combines the frequency channels in units of the frequency channel group for the signal after the interference suppression processing;
a direct wave component extraction unit that extracts a direct wave component from the signal obtained by combining the frequency channels;
a radio wave arrival angle estimating unit that estimates an arrival angle of a radio wave from the direct wave component;
A radio wave arrival angle estimation device comprising:
前記周波数チャネル結合部は、前記周波数チャネル群毎に前記周波数チャネルを結合し、
前記直接波成分抽出部は、前記周波数チャネル群毎に直接波成分を抽出し、
前記電波到来角推定部は、前記周波数チャネル群毎に前記到来角を推定し、
前記周波数チャネル群毎の複数の前記到来角の推定結果を合成する到来角合成部、
をさらに備えることを特徴とする請求項1に記載の電波到来角推定装置。 the interference suppression processing unit generates the interference suppression weight for each of the plurality of frequency channel groups and performs the interference suppression processing;
The frequency channel combining unit combines the frequency channels for each of the frequency channel groups,
The direct wave component extraction unit extracts a direct wave component for each of the frequency channel groups,
the radio wave arrival angle estimation unit estimates the arrival angle for each of the frequency channel groups;
an arrival angle synthesis unit that synthesizes the estimation results of the plurality of arrival angles for each of the frequency channel groups;
The radio wave arrival angle estimation device according to claim 1 , further comprising:
をさらに備え、
前記干渉抑圧処理部、前記周波数チャネル結合部、前記直接波成分抽出部、および前記電波到来角推定部のそれぞれは、選択された前記周波数チャネル群に対して処理を行い、選択されなかった前記周波数チャネル群に対しては処理を行わないことを特徴とする請求項1に記載の電波到来角推定装置。 a frequency channel group selection unit that selects the frequency channel group to be processed from the generated frequency channel group;
Further equipped with
2. The radio wave arrival angle estimation device according to claim 1, wherein each of the interference suppression processing unit, the frequency channel combining unit, the direct wave component extraction unit, and the radio wave arrival angle estimation unit performs processing on the selected frequency channel group and does not perform processing on the unselected frequency channel group.
前記干渉情報に基づいて、干渉の有無と干渉源とに基づいて前記周波数チャネルをグループ分けをする干渉クラスタリング部と、
前記グループ分けの結果から前記周波数チャネル群を決定する周波数チャネル群決定部と、
を有することを特徴とする請求項1から4のいずれか1項に記載の電波到来角推定装置。 The frequency channel group generation unit
an interference clustering unit that groups the frequency channels based on the presence or absence of interference and an interference source based on the interference information;
a frequency channel group determination unit that determines the frequency channel group based on the result of the grouping;
The radio wave arrival angle estimation device according to any one of claims 1 to 4 , further comprising:
前記周波数チャネル群決定部は、干渉を受けている前記グループと干渉を受けていない前記グループに属する周波数チャネルとを周波数チャネルが隣接するように1つのグループにまとめたときに周波数チャネル数が最大になるように前記周波数チャネル群を決定することを特徴とする請求項7に記載の電波到来角推定装置。 The interference clustering unit divides the frequency channels into a group of interfered frequency channels and a group of not interfered frequency channels, and further divides the group of interfered frequency channels into groups for each interference source;
8. The radio wave arrival angle estimation device according to claim 7, wherein the frequency channel group determination unit determines the frequency channel group such that a number of frequency channels is maximized when the group experiencing interference and the frequency channels belonging to the group not experiencing interference are grouped together such that the frequency channels are adjacent to each other.
をさらに備えることを特徴とする請求項1から4のいずれか1項に記載の電波到来角推定装置。 a frequency conversion unit that converts the frequency of the signal to be transmitted according to the frequency channel to be transmitted from among a plurality of the frequency channels by a frequency hopping method;
The radio wave arrival angle estimation device according to any one of claims 1 to 4 , further comprising:
をさらに備え、
前記周波数変換部は、前記周波数チャネル群生成部が生成した前記周波数チャネル群に基づく切替パターンで、送信する前記周波数チャネルを切り替えることを特徴とする請求項3に記載の電波到来角推定装置。 a frequency conversion unit that converts the frequency of the signal to be transmitted according to the frequency channel to be transmitted among the plurality of frequency channels by a frequency hopping method;
Further equipped with
4. The radio wave arrival angle estimation device according to claim 3, wherein the frequency conversion unit switches the frequency channel to be transmitted in a switching pattern based on the frequency channel group generated by the frequency channel group generation unit.
前記周波数チャネル毎に干渉の状態を測定して干渉情報を生成するステップと、
前記周波数チャネル毎に生成された前記干渉情報に基づいて、複数の前記周波数チャネルをグループ分けして周波数チャネル群を生成するステップと、
前記周波数チャネル群単位で干渉抑圧ウェイトを生成して干渉抑圧処理を行うステップと、
前記干渉抑圧処理を行った後の信号に対して、直接波成分を抽出するステップと、
前記直接波成分から電波の到来角を推定するステップと、
を含むことを特徴とする電波到来角推定方法。 A radio wave arrival angle estimation method executed by a radio wave arrival angle estimation device that estimates an arrival direction of a radio wave from a received signal that is received using a plurality of frequency channels, comprising:
measuring an interference state for each of the frequency channels to generate interference information;
generating a frequency channel group by grouping the plurality of frequency channels based on the interference information generated for each of the frequency channels;
generating an interference suppression weight for each frequency channel group and performing an interference suppression process;
extracting a direct wave component from the signal after the interference suppression processing;
estimating an angle of arrival of radio waves from the direct wave component;
A radio wave arrival angle estimation method comprising:
前記周波数チャネル毎に干渉の状態を測定して干渉情報を生成するステップと、
前記周波数チャネル毎に生成された前記干渉情報に基づいて、複数の前記周波数チャネルをグループ分けして周波数チャネル群を生成するステップと、
前記周波数チャネル群単位で干渉抑圧ウェイトを生成して干渉抑圧処理を行うステップと、
前記干渉抑圧処理を行った後の信号に対して、直接波成分を抽出するステップと、
前記直接波成分から電波の到来角を推定するステップと、
を前記電波到来角推定装置に実行させることを特徴とする制御回路。 A control circuit for controlling a radio wave arrival angle estimation device that estimates an arrival direction of a radio wave from a received signal that has been transmitted using a plurality of frequency channels, comprising:
measuring an interference state for each of the frequency channels to generate interference information;
generating a frequency channel group by grouping the plurality of frequency channels based on the interference information generated for each of the frequency channels;
generating an interference suppression weight for each frequency channel group and performing an interference suppression process;
extracting a direct wave component from the signal after the interference suppression processing;
estimating an angle of arrival of radio waves from the direct wave component;
The control circuit causes the radio wave arrival angle estimation device to execute the above.
前記周波数チャネル毎に干渉の状態を測定して干渉情報を生成するステップと、
前記周波数チャネル毎に生成された前記干渉情報に基づいて、複数の前記周波数チャネルをグループ分けして周波数チャネル群を生成するステップと、
前記周波数チャネル群単位で干渉抑圧ウェイトを生成して干渉抑圧処理を行うステップと、
前記干渉抑圧処理を行った後の信号に対して、直接波成分を抽出するステップと、
前記直接波成分から電波の到来角を推定するステップと、
を前記電波到来角推定装置に実行させることを特徴とする記憶媒体。 A storage medium storing a program for controlling a radio wave arrival angle estimation device that estimates an arrival direction of a radio wave from a received signal that has been transmitted using a plurality of frequency channels, the program comprising:
measuring an interference state for each of the frequency channels to generate interference information;
generating a frequency channel group by grouping the plurality of frequency channels based on the interference information generated for each of the frequency channels;
generating an interference suppression weight for each frequency channel group and performing an interference suppression process;
extracting a direct wave component from the signal after the interference suppression processing;
estimating an angle of arrival of radio waves from the direct wave component;
The radio wave arrival angle estimation device according to claim 1,
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/JP2023/011572 WO2024195125A1 (en) | 2023-03-23 | 2023-03-23 | Radio wave arrival angle estimation device, radio wave arrival angle estimation method, control circuit, and storage medium |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JPWO2024195125A1 JPWO2024195125A1 (en) | 2024-09-26 |
| JPWO2024195125A5 JPWO2024195125A5 (en) | 2025-02-26 |
| JP7638460B2 true JP7638460B2 (en) | 2025-03-03 |
Family
ID=92841546
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2024568159A Active JP7638460B2 (en) | 2023-03-23 | 2023-03-23 | Radio wave arrival angle estimation device, radio wave arrival angle estimation method, control circuit, and storage medium |
Country Status (3)
| Country | Link |
|---|---|
| JP (1) | JP7638460B2 (en) |
| TW (1) | TW202438918A (en) |
| WO (1) | WO2024195125A1 (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002502175A (en) | 1998-01-30 | 2002-01-22 | モトローラ・インコーポレイテッド | Method and apparatus for performing target interference suppression |
| JP2016148590A (en) | 2015-02-12 | 2016-08-18 | 株式会社日本自動車部品総合研究所 | Direction of arrival estimation system, measuring device |
| JP2017204662A (en) | 2016-05-09 | 2017-11-16 | 三菱電機株式会社 | Receiving device, receiving method, program, and recording medium |
| US20180131540A1 (en) | 2016-11-04 | 2018-05-10 | Qualcomm Incorporated | Estimating timing and angle information of wireless signals |
-
2023
- 2023-03-23 JP JP2024568159A patent/JP7638460B2/en active Active
- 2023-03-23 WO PCT/JP2023/011572 patent/WO2024195125A1/en not_active Ceased
- 2023-09-15 TW TW112135273A patent/TW202438918A/en unknown
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002502175A (en) | 1998-01-30 | 2002-01-22 | モトローラ・インコーポレイテッド | Method and apparatus for performing target interference suppression |
| JP2016148590A (en) | 2015-02-12 | 2016-08-18 | 株式会社日本自動車部品総合研究所 | Direction of arrival estimation system, measuring device |
| JP2017204662A (en) | 2016-05-09 | 2017-11-16 | 三菱電機株式会社 | Receiving device, receiving method, program, and recording medium |
| US20180131540A1 (en) | 2016-11-04 | 2018-05-10 | Qualcomm Incorporated | Estimating timing and angle information of wireless signals |
Non-Patent Citations (1)
| Title |
|---|
| AYYALASOMAYAJULA, Roshan et al.,"BLoc: CSI-based accurate localization for BLE tags",PROCEEDINGS OF THE 14TH INTERNATIONAL CONFERENCE ON EMERGING NETWORKING EXPERIMENTS AND TECHNOLOGIES,2018年12月04日,Pages 126-138,<URL: https://doi.org/10.1145/3281411.3281428 > |
Also Published As
| Publication number | Publication date |
|---|---|
| TW202438918A (en) | 2024-10-01 |
| WO2024195125A1 (en) | 2024-09-26 |
| JPWO2024195125A1 (en) | 2024-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10505620B2 (en) | Receiving apparatus and receiving method, and program and recording medium | |
| US7376428B2 (en) | Positioning method and radio system | |
| JP4339801B2 (en) | Direction-of-arrival estimation method and reception beam forming apparatus without using eigenvalue decomposition | |
| JP5337111B2 (en) | Signal detection apparatus and signal detection method used in a radio station of a radio communication system | |
| CN106597368B (en) | The localization method and system of indoor interference source | |
| EP0992813A2 (en) | Radio wave arrival direction estimating antenna apparatus | |
| JP3600459B2 (en) | Method and apparatus for estimating direction of arrival of radio wave | |
| RU2232485C2 (en) | Procedure to form directivity pattern of antenna and device for its realization | |
| CA2426650A1 (en) | Method and apparatus for producing a channel estimate with reduced rank | |
| Wang et al. | RIS-enabled integrated sensing and communication for 6G systems | |
| US12040921B1 (en) | Systems and methods for calculating beamforming weights used in wireless network discovery, synchronization, and reference signal waveform identification | |
| JP6211226B2 (en) | Receiving apparatus and receiving method | |
| Kaplan et al. | Detection, identification, and direction of arrival estimation of drone FHSS signals with uniform linear antenna array | |
| JP2013047684A (en) | Method for detecting navigation beacon signals using two antennas or equivalent thereof | |
| JP2001203620A (en) | Wireless base station device and wireless communication method | |
| JP4070074B2 (en) | Method for obtaining transmission gain function and transmission device | |
| JP7638460B2 (en) | Radio wave arrival angle estimation device, radio wave arrival angle estimation method, control circuit, and storage medium | |
| JP6747688B2 (en) | Source estimation method and source estimation apparatus using the same | |
| JP5392994B2 (en) | Radar signal processing apparatus and radar signal processing method | |
| Gentile et al. | Multipath and NLOS mitigation algorithms | |
| JP4070073B2 (en) | Method for obtaining transmission gain function and transmission / reception device | |
| JP2002204193A (en) | Mobile communication system | |
| JP7638466B1 (en) | Receiver and receiving method | |
| US20250310933A1 (en) | Method for multiple antenna joint time of arrival and angle of arrival estimation | |
| JPWO2024195125A5 (en) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20241114 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20241114 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20241114 |
|
| 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: 20250121 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20250218 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7638460 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |