JP2014206440A - Wave source azimuth estimation system, wave source azimuth estimation method and program - Google Patents
Wave source azimuth estimation system, wave source azimuth estimation method and program Download PDFInfo
- Publication number
- JP2014206440A JP2014206440A JP2013083787A JP2013083787A JP2014206440A JP 2014206440 A JP2014206440 A JP 2014206440A JP 2013083787 A JP2013083787 A JP 2013083787A JP 2013083787 A JP2013083787 A JP 2013083787A JP 2014206440 A JP2014206440 A JP 2014206440A
- Authority
- JP
- Japan
- Prior art keywords
- wave
- wave source
- source direction
- source azimuth
- measurement
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005259 measurement Methods 0.000 claims abstract description 44
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 17
- 230000005855 radiation Effects 0.000 claims abstract description 11
- 238000001228 spectrum Methods 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000005315 distribution function Methods 0.000 claims description 9
- 238000011156 evaluation Methods 0.000 claims description 6
- 230000001902 propagating effect Effects 0.000 claims description 5
- 230000002452 interceptive effect Effects 0.000 claims 2
- 230000010287 polarization Effects 0.000 description 11
- 230000005684 electric field Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
本発明は、波源方位推定システム、波源方位推定方法及びプログラムに関する。 The present invention relates to a wave source direction estimation system, a wave source direction estimation method, and a program.
近年は、機密情報漏洩の防止やテレビ、ラジオ等の外来電磁波の遮蔽を目的として、電磁波シールドを施す部屋が増加し、更に大規模化している。 In recent years, the number of rooms for electromagnetic wave shielding has been increased and the scale has been further increased for the purpose of preventing leakage of confidential information and shielding external electromagnetic waves such as televisions and radios.
このような背景に関連する技術としては、様々なものが知られている（例えば、特許文献１参照。）。 Various techniques related to such a background are known (see, for example, Patent Document 1).
例えば、特許文献１には、電磁波シールドの外側から電磁波を照射し、電磁波シールドの内側で電磁波シールドからの漏洩波を受信して電磁波シールドの漏洩波源位置を特定する電磁波シールドにおける漏洩波源位置の特定方法が記載されている。より具体的に説明すると、この電磁波シールドにおける漏洩波源位置の特定方法は、電磁波シールドからの漏洩波を複数の素子を持つアレーアンテナで受信し、各素子で受信された漏洩波の位相差を利用して到来方向推定法により漏洩波の到来方向を推定し、この到来方向から漏洩波源位置を特定する。このようにして、この電磁波シールドにおける漏洩波源位置の特定方法によっては、一度の測定で複数の漏洩箇所を特定することができる。 For example, Patent Document 1 discloses that a leakage wave source position in an electromagnetic wave shield is specified by irradiating an electromagnetic wave from the outside of the electromagnetic wave shield, receiving a leakage wave from the electromagnetic wave shield inside the electromagnetic wave shield, and specifying a leakage wave source position of the electromagnetic wave shield. A method is described. More specifically, this method for identifying the position of the leaky wave source in the electromagnetic wave shield receives the leaky wave from the electromagnetic wave shield with an array antenna having a plurality of elements, and uses the phase difference of the leaky wave received by each element. Then, the arrival direction of the leaky wave is estimated by the arrival direction estimation method, and the leaky wave source position is specified from this arrival direction. In this way, depending on the method of specifying the leaky wave source position in this electromagnetic wave shield, a plurality of leaking locations can be specified by a single measurement.
また、近年は、１〜３ＧＨｚ帯の電波を利用した機器、システムが着実な普及を見せている上、パーソナルコンピュータのクロック周波数の高速化も急速に進み、数ＧＨｚで動作する電子機器が急増している。 In recent years, devices and systems using radio waves in the 1 to 3 GHz band have been steadily spreading, and the speed of the clock frequency of personal computers has rapidly increased, and the number of electronic devices operating at several GHz has increased rapidly. ing.
このような背景に関連する技術としては、様々なものが知られている（例えば、特許文献２参照。）。 Various techniques relating to such a background are known (see, for example, Patent Document 2).
例えば、特許文献２には、電子機器から漏洩する電波の波源位置を推定する電波源位置の推定方法が記載されている。より具体的に説明すると、この電波源位置推定装置は、波源を推定するＳＰＭ法において、計測された電界分布の時空間情報を共分散行列化し、この共分散行列を固有値分解して得られる信号部分空間の最大固有値に伴う固有ベクトルを参照データとして用いる。そして、この電波源位置推定装置は、高周波の電波源を電気ダイポールの集合体として表現する。このようにして、この電波源位置推定装置によっては、波源位相と計測ノイズの影響を除去することができ、計算量が時間サンプリング数に依存しない波源空間分布を推定することができる。 For example, Patent Document 2 describes a radio wave source position estimation method for estimating a radio wave source position leaking from an electronic device. More specifically, this radio wave source position estimation apparatus is a signal obtained by converting the spatio-temporal information of the measured electric field distribution into a covariance matrix and eigenvalue decomposition of the covariance matrix in the SPM method for estimating the wave source. The eigenvector associated with the maximum eigenvalue of the subspace is used as reference data. The radio wave source position estimation apparatus expresses a high frequency radio wave source as an assembly of electric dipoles. In this way, with this radio wave source position estimating apparatus, the influence of the wave source phase and measurement noise can be removed, and the wave source spatial distribution whose calculation amount does not depend on the number of time samplings can be estimated.
また、測定誤差を解消した技術がある。 There is also a technology that eliminates measurement errors.
このような背景に関連する技術としては、様々なものが知られている（例えば、特許文献３参照。）。 Various techniques related to such a background are known (see, for example, Patent Document 3).
例えば、特許文献３には、周波数スペクトラム生成手段と、電波識別手段と、到来波推定手段とを備えた電波到来方向推定装置が記載されている。より具体的に説明すると、周波数スペクトラム生成手段は、複数のアンテナで受信された到来波の受信信号の周波数スペクトラムを生成する。電波識別手段は、信号レベルを比較することによって、周波数スペクトラムに含まれる到来波を識別する。到来波推定手段は、電波識別手段によって到来波が検出された周波数帯域において、到来波の到来数を推定すると共に、ＭＵＳＩＣ（Ｍｕｌｔｉｐｌｅ Ｓｉｇｎａｌ Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）法を用いて到来波の到来方向を推定する。具体的には、到来波推定手段は、到来波の数を指定してＭＵＳＩＣスペクトラムを生成し、ＭＵＳＩＣスペクトラム中のノイズレベルが所定レベル未満となった場合に、指定した到来波の数を到来数と推定する。このようにして、この電波到来方向推定装置によっては、到来波の発信源の位置を高速、且つ、高精度に推定することができる。 For example, Patent Document 3 describes a radio wave arrival direction estimation device including a frequency spectrum generation unit, a radio wave identification unit, and an arrival wave estimation unit. More specifically, the frequency spectrum generation means generates a frequency spectrum of received signals of incoming waves received by a plurality of antennas. The radio wave identification means identifies incoming waves included in the frequency spectrum by comparing the signal levels. The arrival wave estimation means estimates the arrival number of arrival waves in the frequency band in which the arrival wave is detected by the radio wave identification means, and estimates the arrival direction of the arrival waves using a MUSIC (Multiple Signal Classification) method. Specifically, the arrival wave estimating means generates the MUSIC spectrum by designating the number of arrival waves, and when the noise level in the MUSIC spectrum becomes less than a predetermined level, the number of arrival waves designated is the number of arrivals. Estimated. In this way, this radio wave arrival direction estimation apparatus can estimate the position of the transmission source of incoming waves at high speed and with high accuracy.
しかしながら、特許文献１〜３に記載の技術によっては、放射妨害波の波源方位を、短時間に高精度で推定することができない。 However, depending on the techniques described in Patent Documents 1 to 3, the wave source direction of the radiation interference wave cannot be estimated with high accuracy in a short time.
本発明の目的は、上述した課題を解決する波源方位推定システム、波源方位推定方法及びプログラムを提供することにある。 An object of the present invention is to provide a wave source azimuth estimation system, a wave source azimuth estimation method, and a program that solve the above-described problems.
上記課題を解決するために、本発明の第１の形態によると、電子機器から離れた一点に設けられて、電波を受信するアンテナと、妨害波の波源方位を推定する波源方位推定装置とを備え、波源方位推定装置は、アンテナによって受信された電磁界多成分の測定を行う測定部と、測定部による測定値に対してＭＵＳＩＣ法による信号処理を行うことにより、妨害波の波源方位を推定する波源方位推定部とを有する。 In order to solve the above-described problem, according to a first aspect of the present invention, an antenna that receives radio waves and a wave source direction estimation device that estimates a wave source direction of an interference wave are provided at a point apart from an electronic device. The wave source azimuth estimation device estimates the azimuth wave source azimuth by measuring the multi-component electromagnetic field received by the antenna and performing signal processing by the MUSIC method on the measurement value of the measurement unit. And a wave source direction estimating unit.
本発明の第２の形態によると、電子機器から離れた一点に設けられたアンテナによって受信された電磁界多成分の測定を行う測定段階と、得られた測定値に対してＭＵＳＩＣ法による信号処理を行って、電子機器からの放射妨害波の波源方位を推定する信号処理段階とを備える。 According to the second aspect of the present invention, the measurement stage of measuring the electromagnetic field multi-component received by the antenna provided at one point away from the electronic device, and the signal processing by the MUSIC method for the obtained measurement value And a signal processing stage for estimating the source direction of the disturbance wave emitted from the electronic device.
本発明の第３の形態によると、コンピュータを、電子機器から離れた一点に設けられたアンテナによって受信された電磁界多成分の測定を行う測定部、測定部による測定値に対してＭＵＳＩＣ法による信号処理を行うことにより、妨害波の波源方位を推定する波源方位推定部として機能させる。 According to the third aspect of the present invention, a computer uses a MUSIC method to measure a multi-component electromagnetic field received by an antenna provided at one point away from an electronic device, and to measure a value obtained by the measurement unit. By performing signal processing, it functions as a wave source direction estimation unit that estimates the wave source direction of the interference wave.
なおまた、上記の発明の概要は、本発明の必要な特徴の全てを列挙したものではない。また、これらの特徴群のサブコンビネーションもまた、発明となり得る。 The above summary of the invention does not enumerate all necessary features of the present invention. Also, a sub-combination of these feature groups can also be an invention.
以上の説明から明らかなように、この発明によっては、放射妨害波の波源方位を、短時間に高精度で推定することができる。 As is apparent from the above description, according to the present invention, the wave source direction of the radiation interference wave can be estimated with high accuracy in a short time.
以下、発明の実施の形態を通じて本発明を説明するが、以下の実施形態は特許請求の範囲にかかる発明を限定するものではなく、また、実施形態の中で説明されている特徴の組み合わせの全てが発明の解決手段に必須であるとは限らない。 Hereinafter, the present invention will be described through embodiments of the invention. However, the following embodiments do not limit the claimed invention, and all combinations of features described in the embodiments are described below. However, this is not always essential for the solution of the invention.
図１は、一実施形態に係る波源方位推定システム１００のブロック構成の一例を示す。波源方位推定システム１００は、アンテナ１１０及び波源方位推定装置１２０を備える。
FIG. 1 shows an example of a block configuration of a wave source
アンテナ１１０は、空間の電波を高周波エネルギーへ変換する装置である。例えば、アンテナ１１０は、電子機器ＥＤから離れた一点に設けられている。
The
波源方位推定装置１２０は、電子機器ＥＤからの放射妨害波の波源方位を推定する装置である。例えば、波源方位推定装置１２０は、アンテナ１１０と電気的に接続されている。
The wave source
波源方位推定装置１２０は、アンテナ１１０により受信した電磁界多成分の測定を行う測定部１２１と、測定部１２１によって得られた測定値に対してＭＵＳＩＣ法による信号処理を行う波源方位推定部１２２とを備える。
The wave source
図２は、波源方位推定方法を説明するフローチャートである。図２に示すように、まず、波源方位推定装置１２０は、電子機器ＥＤから放射妨害波が発生していることを検出する（Ｓ１０１）。そして、波源方位推定装置１２０は、電子機器ＥＤから発生している放射妨害波をアンテナ１１０で受信する（Ｓ１０２）。次に、測定部１２１は、測定値Ｅを得る（Ｓ１０３）。そして、波源方位推定部１２２は、測定値Ｅを用いて、ＭＵＳＩＣ法を用いた信号処理を行ってスペクトルマトリクスＳを作成する（Ｓ１０４）。波源方位推定部１２２は、作成したスペクトルマトリクスＳから固有値λと固有ベクトルとを計算する（Ｓ１０５）。そして、波源方位推定部１２２は、算出した固有値λと固有ベクトルとを用いて放射妨害波の波源数Ｄを導出する（Ｓ１０６）。波源方位推定部１２２は、最小固有値λｍｉｎに対応する固有ベクトルを用いて評価関数ＰＭＵ（θ）を計算する（Ｓ１０７）。そして、波源方位推定部１２２は、ＰＭＵ（θ）のピークを波源数のＤ個とり、そのときの方位（θ，φ）を求めることにより放射妨害波の波源方位を推定する。
FIG. 2 is a flowchart illustrating the wave source direction estimation method. As shown in FIG. 2, first, the wave source
ここで、測定部１２１は、以下に示す具体的な式を用いて測定を行う。周波数ω、入射角θ、方位角φで伝搬してくる平面波のスペクトルマトリクスの（ｉ，ｊ）成分（ｉ，ｊ＝１，２，…，６）は、波動分布関数を用いて以下の式（１）で表される。
Here, the
ここで、ａ、ｉ、ｊは、積分核と呼ばれ、単位立体角あたりのエネルギーを表す。 Here, a, i, and j are called integral kernels and represent energy per unit solid angle.
次に、波源方位推定部１２２は、以下に示す具体的な式を用いて算出処理を行う。固有値の最も小さい解である最小固有値λｍｉｎ、及び、それに対応する固有ベクトルの数は、スペクトルマトリクスＳで表されるＭ個の連立方程式のうち、条件過多な数Ｎを表している。従って、波源数Ｄは、以下の式（２）で求められる。
Next, the wave source
最小固有値λｍｉｎに対応する固有ベクトルは、放射妨害波源からのノイズベクトルＥＮに相当し、ノイズサブスペースを構成する。また、アンテナ１１０で受信する放射妨害波のモードベクトルａ（θ）で張られる空間は、シグナルサブスペースと呼ばれる。評価関数ＰＭＵ（θ）は、モードベクトルからノイズサブスペースへの距離の２乗の逆数で定義される。このとき、距離は、モードベクトルａ（θ）とノイズベクトルＥＮの内積で表すことができ、以下の式（３）で表される。ここで、＊は、共役転置を表す。
Eigenvector corresponding to the smallest eigenvalue λmin corresponds to the noise vector E N from radiated emissions source, constituting the noise subspace. A space spanned by the mode vector a (θ) of the radiation interference wave received by the
自由空間において、電磁波のＴＥ、ＴＭ波は、互いに独立に伝搬できることから、一点の測定点で測定される電磁界の６成分（Ｅｘ、Ｅｙ、Ｅｚ、Ｂｘ、Ｂｙ、Ｂｚ）は、ＴＥ、ＴＭモード入射波の電界、磁界成分（ＴＥ；Ｂ｜｜、Ｅ⊥：ＴＭ；Ｂ⊥、Ｅ｜｜）により、以下の式（４）で表される。 Since the TE and TM waves of electromagnetic waves can propagate independently from each other in free space, the six components (E x , E y , E z , B x , B y , B z) of the electromagnetic field measured at one measurement point. ) Is expressed by the following formula (4) by the electric field and magnetic field components (TE; B || , E⊥ : TM; B⊥ , E || ) of the TE and TM mode incident waves.
ここで、ｃは、自由空間における光速である。電磁界の６成分は、一般化された電磁界成分、ｅ１、ｅ２、ｅ３、…、ｅ６に書き改めると、以下の式（５）で表される。 Here, c is the speed of light in free space. When the six components of the electromagnetic field are rewritten into generalized electromagnetic field components, e 1 , e 2 , e 3 ,..., E 6 , they are expressed by the following formula (5).
ここで、電磁界の偏波Ｐ０は、ＴＥモード入射波とＴＭモード入射波の比とすることにより、以下の式（６）で定義する。 Here, the polarization P 0 of the electromagnetic field is defined by the following formula (6) by setting the ratio of the TE mode incident wave and the TM mode incident wave.
ここで、偏波Ｐ０は、偏波比ｒ、偏波角δにより複素表示される。ｒは、０〜１の間の値をとる。δは、偏波角と呼ばれ、直線偏波、楕円偏波、円偏波などの偏波状態を表す量であり、−９０°から９０°の間の値をとる。電磁界の６成分は、偏波Ｐ０を用いて式（４）を書き改めると、以下の式（７）で表される。 Here, the polarization P 0 is displayed in a complex manner by the polarization ratio r and the polarization angle δ. r takes a value between 0 and 1. δ is called a polarization angle, and is a quantity representing a polarization state such as linear polarization, elliptical polarization, or circular polarization, and takes a value between −90 ° and 90 °. The six components of the electromagnetic field are expressed by the following equation (7) when the equation (4) is rewritten using the polarization P 0 .
また、電磁界成分の積ｅｉ・ｅｊとエネルギー密度ρとの比は、以下の式（８）のように定義される。 Further, the ratio between the product ei · ej of the electromagnetic field component and the energy density ρ is defined as the following equation (8).
このとき、＊は、複素共役を表す。ここで、各周波数ω、波数ベクトルｋで自由空間を伝搬する平面電磁波のエネルギー密度ρにおいてポインティングベクトル（時間平均）は、以下の式（９）として表される。 At this time, * represents a complex conjugate. Here, the pointing vector (time average) at the energy density ρ of the planar electromagnetic wave propagating in free space at each frequency ω and wave number vector k is expressed as the following equation (9).
また、自由空間での群速度ｖｇは、以下の式（１０）で与えられる。 Further, the group velocity v g of the free space is given by the following equation (10).
ポインティングベクトルは、単位面積を単位時間当たりに通過する電磁波のエネルギーを表し、群速度は、エネルギーがどのくらいの速度で移動するかを表している。従って、ポインティングベクトルは、以下の式（１１）で得られ、群速度は、以下の式（１２）で得られる。 The pointing vector represents the energy of an electromagnetic wave passing through a unit area per unit time, and the group velocity represents how fast the energy moves. Therefore, the pointing vector is obtained by the following equation (11), and the group velocity is obtained by the following equation (12).
よって、式（１１）及び式（１２）より、積分核ａｉｊ（ω、θ、φ；Ｐ０）は、以下の式（１３）のように、偏波Ｐ、入射角θ、方位角φの関数で表される。 Therefore, from Equation (11) and Equation (12), the integral kernel a ij (ω, θ, φ; P 0 ) is polarized P, incident angle θ, azimuth angle φ as shown in Equation (13) below. It is expressed by the function of
これにより、波源方位推定システム１００は、偏波Ｐ、入射角θ、方位角φを取得する。
Thereby, the wave source
以上説明したように、本実施の形態の波源方位推定システム１００によれば、波源方位推定装置１２０は、電子機器ＥＤから離れた一点での電磁界多成分の測定値を信号処理するので、多数点での測定時間を必要とせずに、電子機器ＥＤにおける放射妨害波の波源方位を短時間で推定できる。
As described above, according to the wave source
また、本実施の形態の波源方位推定システム１００によれば、波源方位推定装置１２０は、測定点が一点のため、多数の測定点を設定することによる測定精度の変動がなく、高精度で電子機器ＥＤにおける放射妨害波の波源方位を推定できる。
Further, according to the wave source
そして、本実施の形態の波源方位推定システム１００によれば、波源方位推定装置１２０は、測定部１２１がアンテナ１１０により受信した電磁界多成分の測定を行い、波源方位推定部１２２が測定部１２１による測定値に対してＭＵＳＩＣ法による信号処理を行うので、簡素な回路構成にできる。
And according to the wave source
そしてまた、本実施の形態の波源方位推定システム１００によれば、波源方位推定部１２２は、測定部１２１の測定値に基づいてスペクトルマトリクスを作成して周波数とエネルギーとの関係を求めて、スペクトルマトリクスから固有値と固有ベクトルとを計算する。また、波源方位推定システム１００によれば、波源方位推定部１２２は、計算した固有値と固有ベクトルとを用いて放射妨害波の波源数を導出し、最小固有値に対応する固有ベクトルを用いて評価関数を計算し、評価関数のピークを波源数の数だけとって、そのときの方位を求める。従って、波源方位推定システム１００によれば、波源方位推定装置１２０は、高精度で電子機器ＥＤにおける放射妨害波の波源方位を推定できる。
Moreover, according to the wave source
そしてさらに、本実施の形態の波源方位推定システム１００によれば、波源方位推定部１２２は、自遊空間を伝搬する平面電磁波のエネルギー密度に基づいて、放射妨害波の波源数を算出する。従って、波源方位推定システム１００によれば、波源方位推定装置１２０は、高精度で電子機器ＥＤにおける放射妨害波の波源方位を推定できる。
Furthermore, according to the wave source
そして、本実施の形態の波源方位推定方法によれば、波源方位推定装置１２０は、電子機器ＥＤから離れた一点での電磁界多成分の測定値を信号処理するので、多数点での測定時間を必要とせずに、電子機器ＥＤにおける放射妨害波の波源方位を短時間で推定できる。
And according to the wave source azimuth estimation method of the present embodiment, the wave source
本実施の形態のプログラムによれば、波源方位推定装置１２０は、電子機器ＥＤから離れた一点での電磁界多成分の測定値を信号処理するので、多数点での測定時間を必要とせずに、電子機器ＥＤにおける放射妨害波の波源方位を短時間で推定できる。
According to the program of the present embodiment, the wave source
他の実施の形態の波源方位推定方法は、分布関数として評価する信号処理手法を用いることにより、波源のエネルギー分布を推定する。 The wave source direction estimation method according to another embodiment estimates the energy distribution of the wave source by using a signal processing method that is evaluated as a distribution function.
波源方位推定部１２２は、測定値Ｅを用いてスペクトルマトリクスＳを作成し、周波数とエネルギーとの関係を求める。周波数ω、入射角θ、方位角φで伝搬してくる平面波のスペクトルマトリクスの（ｉ、ｊ）成分（ｉ、ｊ＝１、２、…、６）は、波動分布関数を用いることにより、上記式（１）で表される。ここで、ａｉｊは、積分核と呼ばれ、単位立体角あたりのエネルギーを表す。積分核は、固定された任意の周波数で考えられるので、ωを省略して独立な成分に対する式に書き改めると、以下の式（１４）で表される。
The wave source
Ｐｋは、スペクトルマトリクスＳｉｊ（ω）の、ｑｋは、核ａｉｊ（ω、θ、φ）ｓｉｎθの実部及び虚部にそれぞれ対応しており、以下の式（１５）のような関係にある。 P k corresponds to the spectrum matrix S ij (ω), q k corresponds to the real part and the imaginary part of the nucleus a ij (ω, θ, φ) sin θ, respectively, and There is a relationship.
電界の３成分が測定された場合、積分方程式の個数は、Ｎ＝９となる。一般に、逆問題に対する解は無数に存在するが、その全てが物理的に意味のある解ではない。波動分布関数Ｆ（θ、φ）は、波数空間におけるは導エネルギーの分布を表しているので、滑らかで、至るところ正となる解を求める必要がある。「波動分布関数Ｆ（θ、φ）が未知である。」という状態は、高いエントロピーを持っているといえる。つまり、「波動分布関数Ｆ（θ、φ）を推定する」という作業は、この状態の持つ高いエントロピーを下げる作業であるといえる。従って、エントロピーは、最も効率よく下げられるような分布Ｆ（θ、φ）に推定される。すなわち、推定された分布Ｆ（θ、φ）は、最大となる分布を解として求められる。このようにして得られた解は、予測される最も「ありふれた」分布を表していると解釈できる。波動分布関数Ｆ（θ、φ）のエントロピーＨは、以下の式（１６）で与えられる。 When the three components of the electric field are measured, the number of integral equations is N = 9. In general, there are countless solutions to inverse problems, but not all of them are physically meaningful solutions. Since the wave distribution function F (θ, φ) represents the distribution of the conduction energy in the wave number space, it is necessary to obtain a smooth solution that is positive everywhere. It can be said that the state “the wave distribution function F (θ, φ) is unknown” has high entropy. That is, it can be said that the operation of “estimating the wave distribution function F (θ, φ)” is an operation of reducing the high entropy of this state. Therefore, the entropy is estimated to have a distribution F (θ, φ) that can be lowered most efficiently. That is, the estimated distribution F (θ, φ) is obtained using the maximum distribution as a solution. The solution thus obtained can be interpreted as representing the most “common” distribution to be predicted. The entropy H of the wave distribution function F (θ, φ) is given by the following equation (16).
式（１４）の条件の下に、エントロピーＨを最大とする解は、以下の式（１７）のようにラグランジュ（Ｌａｇｒａｎｇｅ）の未定乗数μｋを用いて得られる。 Under the condition of equation (14), a solution that maximizes entropy H is obtained using Lagrange's undetermined multiplier μk as in equation (17) below.
本実施の形態の波源方位推定方法によれば、波源方位推定装置１２０は、分布関数として評価する信号処理手法を用いて波源のエネルギー分布を推定することにより、電子機器ＥＤにおける放射妨害波の波源方位を短時間で推定できる。
According to the wave source azimuth estimation method of the present embodiment, the wave source
なお、情報端末、アクセスシステム、情報処理方法及びプログラムは、前述した一実施形態に限定するものでなく、適宜な変形や改良等が可能である。 Note that the information terminal, the access system, the information processing method, and the program are not limited to the above-described embodiment, and appropriate modifications and improvements can be made.
１００ 波源方位推定システム
１１０ アンテナ
１２０ 波源方位推定装置
１２１ 測定部
１２２ 波源方位推定部
ＥＤ 電子機器
100 Wave Source
Claims (6)
妨害波の波源方位を推定する波源方位推定装置と
を備え、
前記波源方位推定装置は、
前記アンテナによって受信された電磁界多成分の測定を行う測定部と、
前記測定部による測定値に対してＭＵＳＩＣ法による信号処理を行うことにより、妨害波の波源方位を推定する波源方位推定部と
を有する波源方位推定システム。 An antenna that is provided at one point away from the electronic device and receives radio waves;
A source direction estimation device for estimating the source direction of the jamming wave,
The wave source direction estimating device is:
A measurement unit for measuring an electromagnetic field multi-component received by the antenna;
A wave source azimuth estimation system comprising: a wave source azimuth estimation unit that estimates a wave source azimuth of an interfering wave by performing signal processing using a MUSIC method on a measurement value of the measurement unit.
請求項１に記載の波源方位推定システム。 The wave source azimuth estimation unit creates a spectrum matrix based on the measurement value, obtains a relationship between frequency and energy, calculates an eigenvalue and an eigenvector from the spectrum matrix, and uses the calculated eigenvalue and eigenvector to radiate an interference wave The number of wave sources is calculated, the evaluation function is calculated using the eigenvector corresponding to the minimum eigenvalue, the peak of the evaluation function is taken as many as the number of wave sources, and the direction at that time is obtained. The wave source direction estimation system according to claim 1.
請求項１又は２に記載の波源方位推定システム。 The wave source direction estimation system according to claim 1, wherein the wave source direction estimation unit calculates the number of wave sources of the radiated disturbance wave based on an energy density of a plane electromagnetic wave propagating in a self-playing space.
請求項１又は２に記載の波源方位推定システム。 The wave source direction estimation system according to claim 1, wherein the wave source direction estimation unit evaluates the wave source of the radiation disturbance wave as a wave distribution function.
得られた測定値に対してＭＵＳＩＣ法による信号処理を行って、前記電子機器からの放射妨害波の波源方位を推定する信号処理段階と
を備える波源方位推定方法。 A measurement stage for measuring an electromagnetic multi-component received by an antenna provided at one point away from an electronic device;
A signal source direction estimation method comprising: a signal processing step of performing signal processing on the obtained measurement value by a MUSIC method and estimating a source direction of an interference wave from the electronic device.
電子機器から離れた一点に設けられたアンテナによって受信された電磁界多成分の測定を行う測定部、
前記測定部による測定値に対してＭＵＳＩＣ法による信号処理を行うことにより、妨害波の波源方位を推定する波源方位推定部
として機能させるプログラム。 Computer
A measuring unit for measuring multi-component electromagnetic fields received by an antenna provided at one point away from an electronic device;
A program that functions as a wave source azimuth estimating unit that estimates a wave source azimuth of an interfering wave by performing signal processing using a MUSIC method on a measurement value obtained by the measurement unit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013083787A JP6372874B2 (en) | 2013-04-12 | 2013-04-12 | Source direction estimation system, source direction estimation method, and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013083787A JP6372874B2 (en) | 2013-04-12 | 2013-04-12 | Source direction estimation system, source direction estimation method, and program |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2014206440A true JP2014206440A (en) | 2014-10-30 |
| JP2014206440A5 JP2014206440A5 (en) | 2016-06-09 |
| JP6372874B2 JP6372874B2 (en) | 2018-08-15 |
Family
ID=52120096
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013083787A Active JP6372874B2 (en) | 2013-04-12 | 2013-04-12 | Source direction estimation system, source direction estimation method, and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6372874B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113721201A (en) * | 2021-09-08 | 2021-11-30 | 哈尔滨工程大学 | Estimation method for modulation frequency of linear frequency modulation signal |
| CN115603839A (en) * | 2022-12-12 | 2023-01-13 | 中铁第一勘察设计院集团有限公司（Cn） | Distributed wireless interference source positioning method facing railway communication |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000121716A (en) * | 1998-10-13 | 2000-04-28 | Anritsu Corp | Radio wave propagation estimating equipment |
-
2013
- 2013-04-12 JP JP2013083787A patent/JP6372874B2/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000121716A (en) * | 1998-10-13 | 2000-04-28 | Anritsu Corp | Radio wave propagation estimating equipment |
Non-Patent Citations (2)
| Title |
|---|
| 伊藤聡宏 外２名: ""科学衛星で観測されたプラズマ波動のオンボード到来方向推定システムの検討"", 電子情報通信学会２００６年通信ソサイエティ大会講演論文集１, JPN6017005621, 7 September 2006 (2006-09-07), pages 15, ISSN: 0003503487 * |
| 後藤由貴 外２名: ""ガウス分布モデル波動分布関数法による伝搬ベクトル推定に関する研究"", 電子情報通信学会論文誌Ｂ, vol. Volume J84-B，Number 2, JPN6017005622, 1 February 2001 (2001-02-01), JP, pages 263 - 271, ISSN: 0003503488 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113721201A (en) * | 2021-09-08 | 2021-11-30 | 哈尔滨工程大学 | Estimation method for modulation frequency of linear frequency modulation signal |
| CN113721201B (en) * | 2021-09-08 | 2023-10-13 | 哈尔滨工程大学 | Method for estimating modulation frequency of linear frequency modulation signal |
| CN115603839A (en) * | 2022-12-12 | 2023-01-13 | 中铁第一勘察设计院集团有限公司（Cn） | Distributed wireless interference source positioning method facing railway communication |
| CN115603839B (en) * | 2022-12-12 | 2023-03-07 | 中铁第一勘察设计院集团有限公司 | Distributed wireless interference source positioning method facing railway communication |
Also Published As
| Publication number | Publication date |
|---|---|
| JP6372874B2 (en) | 2018-08-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Russer et al. | Modeling of noisy EM field propagation using correlation information | |
| Russer et al. | An efficient method for computer aided analysis of noisy electromagnetic fields | |
| Tao | Performance analysis for interference and noise canceller based on hypercomplex and spatio‐temporal‐polarisation processes | |
| Malmström et al. | On mutual coupling and coupling paths between antennas using the reaction theorem | |
| Mohammadzadeh et al. | Adaptive beamforming based on theoretical interference‐plus‐noise covariance and direction‐of‐arrival estimation | |
| Zhao et al. | Iteration-free-phase retrieval for directive radiators using field amplitudes on two closely separated observation planes | |
| Arnaut et al. | Empirical emission eigenmodes of printed circuit boards | |
| Mikki et al. | Near-field cross-correlation analysis for MIMO wireless communications | |
| Kitavi et al. | Lower bound of the estimation error of an emitter's direction‐of‐arrival/polarisation, for a collocated triad of orthogonal dipoles/loops that fail randomly | |
| Xian-Dong et al. | A novel in-network caching scheme based on betweenness and replacement rate in content centric networking | |
| Molaei et al. | Two‐dimensional DOA estimation for multi‐path environments by accurate separation of signals using k‐medoids clustering | |
| Ebrahimi et al. | Localisation of mixed near‐field and far‐field sources using the largest aperture sparse linear array | |
| JP6372874B2 (en) | Source direction estimation system, source direction estimation method, and program | |
| Du et al. | Adaptive beamforming algorithm for coprime array based on interference and noise covariance matrix reconstruction | |
| Pan et al. | Source enumeration based on a uniform circular array in a determined case | |
| Zhao et al. | Skeletonization-scheme-based adaptive near field sampling for radio frequency source reconstruction | |
| Goli et al. | Performance analysis of beam scan, MIN-NORM, music and Mvdr DOA estimation algorithms | |
| Yang et al. | Central DOA estimation of incoherently distributed noncircular sources with cross-correlation matrix | |
| Yang et al. | Super‐resolution radar imaging using fast continuous compressed sensing | |
| Li et al. | Robust cyclic MUSIC algorithm for finding directions in impulsive noise environment | |
| Ebrahimi et al. | Generalised two stage cumulants‐based MUSIC algorithm for passive mixed sources localisation | |
| Fontgalland et al. | Combining antenna array elements by using ICA method for remote sensing of sources | |
| Ioannopoulos et al. | Evaluating the effect of small number of snapshots and signal‐to‐noise‐ratio on the efficiency of MUSIC estimations | |
| Hou et al. | Oblique projection for direction‐of‐arrival estimation of hybrid completely polarised and partially polarised signals with arbitrary polarimetric array configuration | |
| Lv et al. | Performance bounds of joint angle‐frequency estimation based on spatio‐temporal co‐prime sampling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160308 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160415 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170210 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20170228 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170404 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20170801 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20171101 |
|
| A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20171109 |
|
| A912 | Re-examination (zenchi) completed and case transferred to appeal board |
Free format text: JAPANESE INTERMEDIATE CODE: A912 Effective date: 20171201 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20180712 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6372874 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |