JP6597505B2 - Radio source identification device and identification system - Google Patents

Radio source identification device and identification system Download PDF

Info

Publication number
JP6597505B2
JP6597505B2 JP2016139890A JP2016139890A JP6597505B2 JP 6597505 B2 JP6597505 B2 JP 6597505B2 JP 2016139890 A JP2016139890 A JP 2016139890A JP 2016139890 A JP2016139890 A JP 2016139890A JP 6597505 B2 JP6597505 B2 JP 6597505B2
Authority
JP
Japan
Prior art keywords
model
inspected
feature
signal
radio wave
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
Application number
JP2016139890A
Other languages
Japanese (ja)
Other versions
JP2018011241A (en
Inventor
哲也 高尾
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2016139890A priority Critical patent/JP6597505B2/en
Publication of JP2018011241A publication Critical patent/JP2018011241A/en
Application granted granted Critical
Publication of JP6597505B2 publication Critical patent/JP6597505B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Circuits Of Receivers In General (AREA)

Description

この発明は、電波源が発する電波を受信して生成した受信信号の特徴量から電波源の機種、可能な場合には個体を識別する装置およびシステムに関する。   The present invention relates to an apparatus and a system for identifying a radio wave source model, and, if possible, an individual from a feature quantity of a received signal generated by receiving a radio wave emitted from the radio wave source.

無線RF信号を発する電波源からの電波をアンテナで受信して生成した受信信号には、無線機などで設定する周波数や帯域幅などの電波諸元の違いの他に、機種ごとまたは個体ごとに、意図しないで伴う特徴が発生することが知られている。個体とは、電波源となるそれぞれの通信装置やレーダ装置などである。この意図しない特徴を解析し無線機の個体を識別するシステムとして、電波源が発する電波を受信して生成した受信信号を収集し、受信信号の信号立ち上がりにおける特徴ベクトル波形を用いて識別する手法が提案されている。この手法では、周波数ごとに予め保存したデータベースにおける既知特徴ベクトル波形との照合により電波源の個体識別を行う(例えば、特許文献1参照)。   The received signal generated by receiving the radio wave from the radio source that emits the radio RF signal with the antenna is different for each model or individual in addition to the difference in radio wave specifications such as the frequency and bandwidth set by the radio etc. It is known that unintentional features occur. An individual is a communication device or a radar device that serves as a radio wave source. As a system for analyzing unintended features and identifying individual wireless devices, there is a method of collecting received signals generated by receiving radio waves emitted from radio sources and identifying them using feature vector waveforms at the rising edge of the received signals. Proposed. In this method, individual identification of radio wave sources is performed by collating with known feature vector waveforms in a database stored in advance for each frequency (see, for example, Patent Document 1).

特許第4760967号Japanese Patent No. 4760967

特許文献1では、受信信号の立ち上がりを用いて個体識別を行うため、収集対象の電波を通信開始時点から補足している必要があるという課題がある。   In patent document 1, since individual identification is performed using a rising edge of a received signal, there is a problem that it is necessary to supplement radio waves to be collected from the communication start time.

この発明は、上記のような問題点を解決するためになされたものであり、通信開始時点から電波が収集できていない場合でも電波源の機種を識別できる電波源の識別装置及び電波源の識別システムを得ることを目的とする。 The present invention has been made to solve the above-described problems, and a radio wave source identification device and radio wave source identification that can identify the radio wave source model even when radio waves have not been collected since the start of communication. The purpose is to obtain a system .

この発明に係る電波源の識別装置及び電波源の識別システムは、アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、ある機種の電波源が送信する電波が受信されて得られる受信信号の複数の特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を電波源の機種ごとに保存する機種別特徴量保存部とを備える。さらに、被検査信号から抽出した特徴量である被検査特徴量と機種別特徴量保存部に保存された機種別特徴量範囲とを比較して、複数の被検査特徴量のそれぞれが機種別特徴量範囲に入る機種を、被検査信号を送信した電波源の機種として識別する機種識別部とを備える。 The radio wave source identification apparatus and radio wave source identification system according to the present invention are configured to transmit at least one feature quantity of a signal to be inspected detected from a discretized reception signal generated from radio waves received by an antenna during at least one communication. A feature value extraction unit that extracts and includes a feature value, and a model-specific feature value range that is a range that each of a plurality of feature values of a reception signal obtained by receiving a radio wave transmitted by a certain type of radio wave source can receive. And a model-specific feature amount storage unit that stores each source model. Further, the feature quantity extracted from the signal to be inspected is compared with the machine-specific feature quantity range stored in the model-specific feature quantity storage unit, and each of the plurality of feature quantities to be inspected is a model-specific feature. A model identifying unit that identifies a model that falls within the quantity range as a model of the radio wave source that transmitted the signal to be inspected;

この発明に係る電波源の識別装置及び電波源の識別システムによれば、通信開始時点から電波が収集できていない場合でも電波源の機種を識別できる。
According to the radio wave source identification device and radio wave source identification system according to the present invention, the radio wave source model can be identified even when radio waves have not been collected since the start of communication.

この発明の実施の形態1に係る電波源の識別システムの構成図である。It is a block diagram of the identification system of the radio wave source which concerns on Embodiment 1 of this invention. 実施の形態1に係る電波源の識別システムの機種別特徴量範囲の例を説明する図である。It is a figure explaining the example of the model-specific feature-value range of the identification system of the radio wave source which concerns on Embodiment 1. FIG. 実施の形態1に係る電波源の識別システムの概略処理を説明するフロー図である。FIG. 3 is a flowchart illustrating a schematic process of the radio wave source identification system according to the first embodiment. 実施の形態1に係る電波源の識別システムの機種識別処理を説明するフロー図である。3 is a flowchart for explaining model identification processing of the radio wave source identification system according to Embodiment 1. FIG. 実施の形態1に係る電波源の識別システムの機種識別処理で通信中での特徴量ごとに候補機種を絞る処理を説明するフロー図である。FIG. 10 is a flowchart for explaining processing for narrowing down candidate models for each feature amount during communication in the model identification process of the radio wave source identification system according to the first embodiment. 実施の形態1に係る電波源の識別システムの機種識別処理で被検査信号の立ち上がりまたは立ち下がりでの特徴量ごとに候補機種を絞る処理を説明するフロー図である。FIG. 6 is a flowchart for explaining processing for narrowing down candidate models for each feature amount at the rise or fall of an inspected signal in the model identification process of the radio wave source identification system according to the first embodiment. 実施の形態1に係る電波源の識別システムの機種識別処理でIQ軌跡を用いて候補機種を絞る処理を説明するフロー図である。6 is a flowchart for explaining processing for narrowing down candidate models using an IQ trajectory in model identification processing of the radio wave source identification system according to Embodiment 1. FIG. 実施の形態1に係る電波源の識別システムの機種識別処理で機種を識別する例を説明する図である。It is a figure explaining the example which identifies a model by the model identification process of the identification system of the radio wave source which concerns on Embodiment 1. FIG. 実施の形態1に係る電波源の識別システムの個体識別処理を説明するフロー図である。3 is a flowchart for explaining individual identification processing of the radio wave source identification system according to Embodiment 1. FIG. 実施の形態1に係る電波源の識別システムの個体識別処理で通信中での特徴量ごとに候補個体を絞る処理を説明するフロー図である。6 is a flowchart for explaining processing for narrowing down candidate individuals for each feature amount during communication in the individual identification processing of the radio wave source identification system according to Embodiment 1. FIG. 実施の形態1に係る電波源の識別システムの個体識別処理で被検査信号の立ち上がりまたは立ち下がりでの特徴量ごとに候補個体を絞る処理を説明するフロー図である。FIG. 6 is a flowchart for explaining processing for narrowing down candidate individuals for each feature amount at the rise or fall of a signal under test in the individual identification processing of the radio wave source identification system according to the first embodiment. 実施の形態1に係る電波源の識別システムの個体識別処理でIQ軌跡を用いて候補個体を絞る処理を説明するフロー図である。FIG. 6 is a flowchart illustrating processing for narrowing down candidate individuals using IQ trajectories in the individual identification processing of the radio wave source identification system according to the first embodiment. この発明の実施の形態2に係る電波源の識別システムの構成図である。It is a block diagram of the identification system of the radio wave source which concerns on Embodiment 2 of this invention. 実施の形態2に係る電波源の識別システムの概略処理を説明するフロー図である。FIG. 10 is a flowchart for explaining schematic processing of the radio wave source identification system according to the second embodiment. 実施の形態2に係る電波源の識別システムの機種識別処理を説明するフロー図である。FIG. 10 is a flowchart illustrating model identification processing of the radio wave source identification system according to the second embodiment. 実施の形態2に係る電波源の識別システムの個体識別処理を説明するフロー図である。FIG. 10 is a flowchart for explaining individual identification processing of the radio wave source identification system according to the second embodiment. この発明の実施の形態3に係る電波源の識別システムの構成図である。It is a block diagram of the identification system of the radio wave source which concerns on Embodiment 3 of this invention. 実施の形態3に係る電波源の識別システムの概略処理を説明するフロー図である。FIG. 10 is a flowchart for explaining schematic processing of the radio wave source identification system according to the third embodiment. 実施の形態3に係る電波源の識別システムの機種識別処理を説明するフロー図である。FIG. 10 is a flowchart illustrating model identification processing of the radio wave source identification system according to the third embodiment. 実施の形態3に係る電波源の識別システムの個体識別処理を説明するフロー図である。FIG. 10 is a flowchart for explaining individual identification processing of the radio wave source identification system according to the third embodiment. この発明の実施の形態4に係る電波源の識別システムの構成図である。It is a block diagram of the identification system of the radio wave source which concerns on Embodiment 4 of this invention. 実施の形態4に係る電波源の識別システムの概略処理を説明するフロー図である。FIG. 10 is a flowchart for explaining schematic processing of the radio wave source identification system according to the fourth embodiment. 実施の形態4に係る電波源の識別システムの機種識別処理を説明するフロー図である。FIG. 10 is a flowchart illustrating model identification processing of the radio wave source identification system according to the fourth embodiment. 実施の形態4に係る電波源の識別システムの個体識別処理を説明するフロー図である。FIG. 10 is a flowchart for explaining individual identification processing of the radio wave source identification system according to the fourth embodiment.

実施の形態1.
この発明の実施の形態1に係る電波源の識別装置について、図1から図12までを参照して説明する。図1は、この発明の実施の形態1に係る電波源の識別システムの構成図である。
Embodiment 1 FIG.
A radio wave source identification device according to Embodiment 1 of the present invention will be described with reference to FIGS. 1 is a configuration diagram of a radio wave source identification system according to Embodiment 1 of the present invention.

図1において、電波源の識別システム50は、電波源からの電波を受信するアンテナ1、アンテナ1が受信した電波から生成したアナログの受信信号が入力されてデジタル化した被検査信号を出力する受信装置2、被検査信号が入力されて被検査信号の電波源を識別する電波源の識別装置3で構成される。披検査信号が電波源を識別する対象の信号である。   In FIG. 1, a radio wave source identification system 50 includes an antenna 1 that receives radio waves from a radio wave source, and an analog reception signal generated from radio waves received by the antenna 1 and receives a digital signal to be inspected. The apparatus 2 includes a radio wave source identification device 3 that receives a signal to be inspected and identifies the radio wave source of the signal to be inspected. The demonstration signal is a signal for identifying a radio wave source.

受信装置2は、アンテナ1で収集したRF信号である受信信号をIF信号に周波数変換する周波数変換器4、IF信号を決められたサンプリングレートで決められたビット数のデジタル信号に変換するA/D変換器5、A/D変換器5が出力する受信信号の中からノイズレベルを超えた信号を直交検波してIQ信号を出力する検波器6を有する。検波器6が出力するIQ信号が、アンテナ1が受信した電波から生成した離散化された受信信号から検出された被検査信号である。   The receiving apparatus 2 has a frequency converter 4 that converts the received signal, which is an RF signal collected by the antenna 1, into an IF signal, and an A / A that converts the IF signal into a digital signal having a predetermined number of bits at a predetermined sampling rate. A D converter 5 and a detector 6 that outputs an IQ signal by quadrature detection of a signal that exceeds the noise level from the received signals output by the A / D converter 5. The IQ signal output from the detector 6 is an inspected signal detected from the discretized received signal generated from the radio wave received by the antenna 1.

電波源の識別装置3は、被検査信号から特徴量を抽出する特徴量抽出部としての周期性検出部7、振幅特性検出部8、位相特性検出部9、周波数特性検出部10、IQ軌跡特性検出部11および帯域幅特性検出部12を有する。さらに、機種ごとに特徴量に関する情報を保存する機種別特徴量保存部13、個体ごとに特徴量に関する情報を保存する個体別特徴量保存部14、被検査信号の電波源の機種を識別する機種識別部15および個体を識別する個体識別部16を有する。   The radio wave source identification device 3 includes a periodicity detection unit 7, an amplitude characteristic detection unit 8, a phase characteristic detection unit 9, a frequency characteristic detection unit 10, an IQ trajectory characteristic as a feature quantity extraction unit that extracts a feature quantity from an inspected signal. A detection unit 11 and a bandwidth characteristic detection unit 12 are included. Further, a model-specific feature amount storage unit 13 that stores information about feature amounts for each model, an individual feature amount storage unit 14 that stores information about feature amounts for each individual, and a model that identifies the radio wave source model of the signal to be inspected It has the identification part 15 and the individual identification part 16 which identifies an individual.

周期性検出部7は、被検査信号と被検査信号を遅延させた信号との積の決められた時間での時間積分である自己相関値を、遅延時間を変えて計算する。自己相関値が極大値をとる遅延時間を特徴量として検出する。特に振幅変調系の変調方式では、通信の情報を載せる際のビットの入れ替わり時に機種または個体ごとの特徴が現れると考えられる。機種または個体ごとの特徴により、被検査信号に周期性が表れる場合がある。その周期性を検出するのが周期性検出部7である。周期性は、通信中の特性を評価する特徴量である。被検査被検査信号の立ち上がりまたは立ち下がりを収集できない場合でも、通信中の被検査信号を収集できれば電波源を識別できる利点がある。   The periodicity detection unit 7 calculates an autocorrelation value, which is a time integration over a predetermined time of a product of the signal to be inspected and a signal obtained by delaying the signal to be inspected, by changing the delay time. A delay time at which the autocorrelation value takes a maximum value is detected as a feature amount. In particular, in the modulation method of the amplitude modulation system, it is considered that characteristics of each model or individual appear when bits are exchanged when carrying communication information. Depending on the characteristics of each model or individual, periodicity may appear in the signal under test. The periodicity detection unit 7 detects the periodicity. Periodicity is a feature value for evaluating characteristics during communication. Even when the rising or falling edge of the signal under inspection cannot be collected, there is an advantage that the radio wave source can be identified if the signal under inspection during communication can be collected.

振幅特性検出部8、位相特性検出部9、周波数特性検出部10およびIQ軌跡特性検出部11は、検出した立ち上がりまたは立ち下がりの被検査信号のそれぞれに対して処理する。ここで、信号の立ち上がりとは、電波源が電波の送信を開始してからの決められた期間である。被検査信号の中である時刻よりも前には検出されていない被検査信号に関して、検出し始める時刻から決められた期間すなわち送信開始時の決められた期間の被検査信号が、立ち上がりの被検査信号である。被検査信号の立ち下がりとは、電波源が電波の送信を終了する前の決められた期間である。被検査信号の中で検出している被検査信号がある時刻よりも後には検出できなくなる場合に、検出できなくなる時刻よりも前の決められた期間すなわち送信終了時の決められた期間の被検査信号が、立ち下がりの被検査信号である。   The amplitude characteristic detection unit 8, the phase characteristic detection unit 9, the frequency characteristic detection unit 10, and the IQ trajectory characteristic detection unit 11 perform processing on each detected signal to be inspected rising or falling. Here, the rise of the signal is a predetermined period after the radio wave source starts transmitting radio waves. With respect to the signal to be inspected that is not detected before a certain time in the signal to be inspected, the signal to be inspected in the period determined from the time to start detection, that is, the period determined at the start of transmission, Signal. The fall of the signal to be inspected is a predetermined period before the radio wave source ends transmission of radio waves. If the signal under test detected in the signal under test cannot be detected after a certain time, the device is inspected for a predetermined period before the time when it cannot be detected, that is, for a predetermined period at the end of transmission. The signal is a falling signal under test.

各機種または個体で用いるRF部品の性能のバラつきなどから、立ち上がりおよび立ち下がりの信号には、その振幅、位相、周波数それぞれの標準偏差、歪度、尖度の統計量に、機種または個体差が現れると考えられる。   Due to variations in the performance of the RF components used in each model or individual, the rise and fall signals have standard deviation, skewness, and kurtosis statistics for their amplitude, phase, and frequency, and there is a model or individual difference. It is thought to appear.

振幅特性検出部8は、立ち上がりの被検査信号が得られている場合には、立ち上がりの被検査信号の振幅の値xiの集合に対して平均μ、標準偏差σ、歪度γ、尖度γを以下のように計算する。振幅の値xiの個数をNとする。ΣはN個の要素の和をとることを意味する。
μ=(Σxi)/N (1)
μ=(Σ(xi−μ))/N (2)
σ=√(μ) (3)
μ=(Σ(xi−μ))/N (4)
γ=μ/(σ) (5)
μ=(Σ(xi−μ))/N (6)
γ=μ/(σ)−3 (7)
立ち下がりの被検査信号が得られている場合には、立ち下がりの被検査信号の振幅の値xiの集合に対して、同様に平均μ、標準偏差σ、歪度γ、尖度γを計算する。
When the rising signal under test is obtained, the amplitude characteristic detector 8 calculates the average μ, standard deviation σ, skewness γ 1 , kurtosis with respect to the set of amplitude values xi of the rising signal under test. γ 2 is calculated as follows: Let N be the number of amplitude values xi. Σ means taking the sum of N elements.
μ = (Σxi) / N (1)
μ 2 = (Σ (xi−μ) 2 ) / N (2)
σ = √ (μ 2 ) (3)
μ 3 = (Σ (xi−μ) 3 ) / N (4)
γ 1 = μ 3 / (σ 3 ) (5)
μ 4 = (Σ (xi−μ) 4 ) / N (6)
γ 2 = μ 4 / (σ 4 ) -3 (7)
When a falling signal under test is obtained, the average μ, standard deviation σ, skewness γ 1 , skewness γ 1 , kurtosis γ 2 are similarly applied to the set of amplitude values xi of the falling signal under test. Calculate

位相特性検出部9は、立ち上がりの被検査信号または立ち下がりの被検査信号のそれぞれ得られている位相の値δiの集合に対して、平均μ、標準偏差σ、歪度γ、尖度γを計算する。 The phase characteristic detector 9 calculates the mean μ, standard deviation σ, skewness γ 1 , skewness γ 1 , and kurtosis γ with respect to the set of phase values δi obtained for the rising or falling inspection signals. 2 is calculated.

周波数特性検出部10は、立ち上がりの被検査信号または立ち下がりの被検査信号のそれぞれ得られている周波数の値fiの集合に対して、平均μ、標準偏差σ、歪度γ、尖度γを計算する。 The frequency characteristic detection unit 10 calculates an average μ, a standard deviation σ, a skewness γ 1 , and a kurtosis γ with respect to a set of frequency values fi obtained from the rising signal or the falling signal to be tested. 2 is calculated.

被検査信号のIQ軌跡は、決められた期間の被検査信号(I信号およびQ信号)を複素平面で表した軌跡である。IQ軌跡特性検出部11は、立ち上がりの被検査信号または立ち下がりの被検査信号のそれぞれ得られているものについて、IQ軌跡を生成する。被検査信号から生成したIQ軌跡を被検査軌跡と呼ぶ。IQ軌跡は、信号の振幅および位相の時間変化から得られる情報である。振幅特性または位相特性だけで機種または個体を判定できない場合にも、IQ軌跡により振幅特性および位相特性を同時考慮することで機種または個体により表れる機種または個体による特徴が検出できる場合がある。   The IQ trajectory of the signal to be inspected is a trajectory representing the signal to be inspected (I signal and Q signal) in a predetermined period on a complex plane. The IQ trajectory characteristic detecting unit 11 generates an IQ trajectory for each of the rising signal to be inspected or the falling signal to be inspected. The IQ trajectory generated from the signal to be inspected is called an inspected trajectory. The IQ trajectory is information obtained from the time change of the amplitude and phase of the signal. Even when the model or individual cannot be determined only by the amplitude characteristic or the phase characteristic, there may be a case where the characteristic by the model or the individual represented by the model or the individual can be detected by simultaneously considering the amplitude characteristic and the phase characteristic from the IQ locus.

帯域幅特性検出部12は、被検査信号の帯域幅を算出する。帯域幅は、通信中の信号から得られる特徴量である。無線機によっては帯域幅の設定を変更可能なものもあるが、帯域幅の変更を行えない機種も存在する。他の特徴量で機種を1個に絞れない場合に、帯域幅の違いで機種を識別できる場合も考えられるので、帯域幅も特徴量とする。なお、帯域幅の設定を変更可能な機種では、取りうる可能性があるすべての帯域幅を、その機種の帯域幅とする。   The bandwidth characteristic detector 12 calculates the bandwidth of the signal to be inspected. The bandwidth is a feature amount obtained from a signal being communicated. Some wireless devices can change the bandwidth setting, but there are models that cannot change the bandwidth. Since it is possible that the model can be identified by the difference in bandwidth when the number of models cannot be reduced to one with other feature values, the bandwidth is also used as the feature value. Note that in a model whose bandwidth setting can be changed, all the possible bandwidths are set as the bandwidth of the model.

帯域幅は、例えば以下のようにして求める。検波器6が検波に使用する周波数は、被検査信号の中で電力レベルがピークとなる周波数である。ピークの電力レベルを基準とした雑音レベルと比較し、雑音レベル以上と判断できる周波数の幅を帯域幅とする。   The bandwidth is obtained as follows, for example. The frequency used by the detector 6 for detection is a frequency at which the power level peaks in the signal to be inspected. Compared with the noise level based on the peak power level, the frequency width that can be determined to be equal to or higher than the noise level is defined as the bandwidth.

機種別特徴量保存部13は、機種ごとに、周期性、振幅特性、位相特性、周波数特性および帯域幅特性である複数の特徴量のそれぞれについて、その機種での特徴量の取りうる範囲である機種別特徴量範囲を保存する。機種別特徴量保存部13には、機種ごとの標準的なIQ軌跡である機種別軌跡も保存される。機種別特徴量範囲は、その機種に属する個体で得られる特徴量の最大値、最小値、平均および標準偏差などから決める。機種別特徴量範囲に特徴量が入らない場合に、決められた確率以上でその機種では無いと判断できるように、機種別特徴量範囲を決める。機種別特徴量範囲に特徴量が入る場合に、決められた確率以上でその機種であると判断できるような、第2の機種別特徴量範囲を決めてもよい。   The model-specific feature amount storage unit 13 is a range that the feature amount of each model can take for each of a plurality of feature amounts that are periodicity, amplitude characteristic, phase characteristic, frequency characteristic, and bandwidth characteristic for each model. Save the model-specific feature range. The model-specific feature amount storage unit 13 also stores a model-specific trajectory that is a standard IQ trajectory for each model. The model-specific feature amount range is determined from the maximum value, minimum value, average, standard deviation, and the like of the feature amounts obtained by individuals belonging to the model. If the feature quantity does not enter the model-specific feature quantity range, the model-specific feature quantity range is determined so that it can be determined that the model is not the model with a predetermined probability or more. The second model-specific feature amount range may be determined so that when the feature amount enters the model-specific feature amount range, it can be determined that the model is more than a predetermined probability.

機種別特徴量範囲の例を図2に示す。図2は、振幅特性に関して、機種Cの機種別特徴量範囲X1、機種Dの機種別特徴量範囲X2、機種Eの機種別特徴量範囲X3を示す。図2では、標準偏差σ、歪度γ、尖度γの3次元空間で、機種別特徴量範囲X1、X2、X3をそれぞれ分離した楕円球として表現している。図2に示すような場合であれば、被検査信号の振幅特性の被検査特徴量(標準偏差σ、歪度γ、尖度γ)と機種別特徴量範囲X1、X2、X3を比較することで、機種C、機種D、機種Eどれか、または他の機種かを識別できる。複数の機種の機種別特徴量範囲が重なる部分を持つ場合もある。 An example of the model-specific feature amount range is shown in FIG. FIG. 2 shows a model-specific feature amount range X1 of model C, a model-specific feature amount range X2 of model D, and a model-specific feature amount range X3 of model E regarding amplitude characteristics. In FIG. 2, in the three-dimensional space of the standard deviation σ, the skewness γ 1 , and the kurtosis γ 2 , the model-specific feature amount ranges X1, X2, and X3 are represented as separated elliptic spheres. In the case as shown in FIG. 2, the feature quantity to be inspected (standard deviation σ, skewness γ 1 , kurtosis γ 2 ) of the amplitude characteristic of the signal to be inspected is compared with the model-specific feature quantity ranges X1, X2, and X3. By doing so, it can be identified whether it is model C, model D, model E, or another model. In some cases, the model-specific feature ranges of a plurality of models have overlapping portions.

個体別特徴量保存部14は、機種に属する個体ごとに、周期性、振幅特性、位相特性、周波数特性および帯域幅特性である複数の特徴量のそれぞれについて、その特徴量の取りうる範囲である個体別特徴量範囲を保存する。個体別特徴量保存部14には、個体ごとのIQ軌跡である個体軌跡も保存される。ここでは、特徴量ごとに検出誤差を決めておき、各個体で検出された特徴量である個体特徴量に検出誤差を加算および減算して個体別特徴量範囲を決めている。個体別特徴量保存部14は個体特徴量も保存している。なお、特徴量を複数回、検出した場合は、それらの平均値を個体特徴量とする。検出された特徴量の範囲から個体特徴範囲を決めてもよい。個体別特徴量保存部14は、個体特徴量を保存しなくてもよい。   The individual feature amount storage unit 14 is a range that the feature amount can take for each of a plurality of feature amounts that are periodicity, amplitude characteristic, phase characteristic, frequency characteristic, and bandwidth characteristic for each individual belonging to the model. Save individual feature range. The individual feature amount storage unit 14 also stores an individual locus that is an IQ locus for each individual. Here, a detection error is determined for each feature amount, and the individual feature amount range is determined by adding and subtracting the detection error to the individual feature amount that is a feature amount detected by each individual. The individual feature amount storage unit 14 also stores individual feature amounts. In addition, when the feature amount is detected a plurality of times, the average value thereof is set as the individual feature amount. The individual feature range may be determined from the detected feature amount range. The individual feature amount storage unit 14 may not store the individual feature amount.

振幅特性、位相特性および周波数特性のそれぞれは、標準偏差、歪度、尖度のすべてについて機種別特徴量範囲を決めなくてもよく、少なくとも1つについて決めればよい。個体別特徴量範囲についても同様である。   Each of the amplitude characteristic, the phase characteristic, and the frequency characteristic may not be determined for each standard deviation, skewness, and kurtosis, but may be determined for at least one. The same applies to the individual feature amount range.

機種識別部15は、周期性検出部7などの特徴量検出部により被検査信号から抽出される特徴量である被検査特徴量と、機種別特徴量保存部13を参照して得られる機種別特徴量範囲とを比較して機種を識別する。特徴量は1個ずつ順番に処理する。最初はすべての機種の可能性がある。ある特徴量に関して、被検査特徴量がその特徴量の機種別特徴量範囲に入らない機種を除いていく。被検査特徴量がその特徴量の機種別特徴量範囲に入る機種が1個だけになれば、その機種を被検査信号の機種と識別する。   The model identification unit 15 is a feature type extracted from a signal to be inspected by a feature amount detection unit such as the periodicity detection unit 7 and a model type obtained by referring to the model-specific feature amount storage unit 13. The model is identified by comparing with the feature value range. The feature values are processed one by one in order. Initially all models are possible. With respect to a certain feature value, models whose inspected feature value does not fall within the model-specific feature value range of the feature value are excluded. If there is only one model whose feature quantity to be inspected falls within the model-specific feature quantity range of the feature quantity, that model is identified as the model of the signal to be inspected.

個体識別部16は、被検査特徴量と、個体別特徴量保存部14を参照して得られる、機種識別部が識別した機種に属する個体の個体別特徴量範囲を比較して個体を識別する。特徴量は1個ずつ順番に処理する。最初は、その機種に属するすべての個体に可能性がある。ある特徴量に関して、被検査特徴量がその特徴量の個体別特徴量範囲に入らない個体を除いていく。被検査特徴量が個体別特徴量範囲に入る個体が1個だけになれば、その個体を被検査信号の個体とする。すべての特徴量で検査しても個体が決められない場合は、被検査特徴量を持つ新たな個体として個体別特徴量保存部14に登録する。   The individual identifying unit 16 identifies the individual by comparing the individual feature amount ranges of the individuals belonging to the model identified by the model identifying unit, which are obtained by referring to the inspected feature amount and the individual feature storing unit 14. . The feature values are processed one by one in order. Initially, all individuals belonging to that model have the potential. With respect to a certain feature amount, individuals whose inspection feature amount does not fall within the individual feature amount range of the feature amount are excluded. If there is only one individual whose inspected feature quantity falls within the individual feature quantity range, that individual is set as the individual of the inspected signal. If an individual cannot be determined even after inspection with all the feature quantities, it is registered in the individual-specific feature quantity storage unit 14 as a new individual having the inspected feature quantity.

動作を説明する。図3は、実施の形態1に係る電波源の識別システムの概略処理を説明するフロー図である。まず、ステップS01で、受信信号から被検査信号を抽出する。ステップS02で、機種判定処理を実施する。ステップS03で、個体判定処理を実施する。なお、ステップS02以降は、検出した被検査信号ごとに処理する。   The operation will be described. FIG. 3 is a flowchart for explaining schematic processing of the radio wave source identification system according to the first embodiment. First, in step S01, a signal to be inspected is extracted from the received signal. In step S02, model determination processing is performed. In step S03, an individual determination process is performed. In addition, after step S02, it processes for every detected to-be-inspected signal.

S01では、受信装置2が以下のようにして被検査信号を検出する。周波数変換器4で、アンテナ1から入力されるRF信号である受信信号をIF信号に周波数変換する。次に、A/D変換器5で、IF信号を決められたサンプリングレートで決められたビット数のデジタル信号の受信信号に変換する。デジタル化した受信信号は記憶装置に保存する。検波器6は、デジタル化した受信信号を直交検波してIQ信号を被検査信号として出力する。被検査信号は記憶装置に保存する。検波器6は、デジタル化した受信信号を異なる周波数で検波するごとに1個の被検査信号を出力する。   In S01, the receiving apparatus 2 detects the signal to be inspected as follows. The frequency converter 4 converts the received signal, which is an RF signal input from the antenna 1, into an IF signal. Next, the A / D converter 5 converts the IF signal into a received signal of a digital signal having a determined number of bits at a determined sampling rate. The digitized received signal is stored in a storage device. The detector 6 performs quadrature detection on the digitized reception signal and outputs an IQ signal as a signal to be inspected. The signal under test is stored in a storage device. The detector 6 outputs one inspected signal every time the digitized received signal is detected at a different frequency.

検出した被検査信号ごとに、S02で、電波源の識別装置3が機種を識別する。図4は、実施の形態1に係る電波源の識別システムの機種識別処理を説明するフロー図である。図4に示す処理フローでは、周期性、振幅特性、位相特性、周波数特性、IQ特性および帯域幅特性を、この順番に考慮する。ステップS11で、すべての機種を被検査信号の機種である可能性がある機種(候補機種と呼ぶ)とする。ステップS12で、周期性についての特徴量を使用して機種を識別する。その後、ステップS13で候補機種が2個以上かどうかチェックする。候補機種が2個以上である場合(S13でYES)は、ステップS14で、振幅特性についての特徴量を使用して機種を識別する。その後、ステップS15で候補機種が2個以上かどうかチェックする。以降、同様にして、ステップS16、S18、S20およびS22で、それぞれ振幅特性、位相特性、周波数特性、IQ特性および帯域幅特性に関して機種識別を実施する。ステップSS17、S19、S21およびS23で、候補機種が2個以上かどうかチェックする。   For each detected signal to be inspected, in S02, the radio wave source identifying device 3 identifies the model. FIG. 4 is a flowchart for explaining model identification processing of the radio wave source identification system according to the first embodiment. In the processing flow shown in FIG. 4, periodicity, amplitude characteristics, phase characteristics, frequency characteristics, IQ characteristics, and bandwidth characteristics are considered in this order. In step S11, all models are set as models that are likely to be models of signals to be inspected (referred to as candidate models). In step S12, the model is identified using the feature value for periodicity. Thereafter, in step S13, it is checked whether there are two or more candidate models. If there are two or more candidate models (YES in S13), the model is identified using the feature quantity for the amplitude characteristic in step S14. Thereafter, in step S15, it is checked whether there are two or more candidate models. Thereafter, in the same manner, in steps S16, S18, S20, and S22, model identification is performed with respect to amplitude characteristics, phase characteristics, frequency characteristics, IQ characteristics, and bandwidth characteristics, respectively. In steps SS17, S19, S21 and S23, it is checked whether there are two or more candidate models.

すべての特徴量を考慮しても候補機種が2個以上である場合(S23がYES)は、S24で、機種識別不能とする。すべての特徴量を考慮した後の候補機種が2個以上でも、被検査信号の機種として識別し、候補機種のどれかに属する個体を個体識別の対象とするようにしてもよい。   If there are two or more candidate models even if all feature quantities are taken into consideration (YES in S23), it is determined that the model cannot be identified in S24. Even if there are two or more candidate models after considering all the feature quantities, the model of the signal to be inspected may be identified, and an individual belonging to any of the candidate models may be an individual identification target.

S13、S15、S17、S19、S21およびS23のいずれかで、候補機種が2個未満である場合は、ステップS25で候補機種が1個かチェックする。候補機種が1個の場合(S25でYES)は、S26で、その1個の機種を、被検査信号を送信した電波源の機種として識別する。候補機種が0個の場合(S25でNO)は、ステップS27で機種識別不能とする。   If there are less than two candidate models in any of S13, S15, S17, S19, S21 and S23, it is checked in step S25 whether there is one candidate model. If there is one candidate model (YES in S25), that one model is identified as the model of the radio wave source that transmitted the signal to be inspected in S26. If the number of candidate models is zero (NO in S25), model identification is not possible in step S27.

各特徴量(IQ軌跡は除く)では、機種別特徴量範囲と被検査特徴量とを比較して、候補機種を絞っていく。図5は、実施の形態1に係る電波源の識別システムの機種識別処理で通信中での特徴量ごとに候補機種を絞る処理を説明するフロー図である。まず、ステップS41で、周期性または帯域幅特性である被検査特徴量を取得する。ステップS42で、候補機種から1個の機種(処理中機種と呼ぶ)を取り出す。ステップS43で、処理中機種の機種別特徴量範囲を機種別特徴量保存部13から取得する。ステップS44で、被検査特徴量が処理中機種の機種別特徴量範囲に含まれるか否かチェックする。含まれない場合(S44でNO)は、ステップS45で、処理中機種を候補機種から除く。ステップS46で、未処理の候補機種があるかチェックする。未処理の候補機種がある場合(S46でYES)は、S42に戻る。未処理の候補機種が無い場合(S46でNO)は、終了する。   For each feature value (excluding IQ trajectory), the model-specific feature value range and the feature value to be inspected are compared to narrow down candidate models. FIG. 5 is a flowchart for explaining processing for narrowing down candidate models for each feature amount during communication in the model identification process of the radio wave source identification system according to the first embodiment. First, in step S41, an inspected feature value that is periodicity or bandwidth characteristics is acquired. In step S42, one model (referred to as a model being processed) is extracted from the candidate models. In step S 43, the model-specific feature amount range of the model being processed is acquired from the model-specific feature amount storage unit 13. In step S44, it is checked whether the inspected feature value is included in the model-specific feature value range of the model being processed. If not included (NO in S44), the model being processed is excluded from the candidate models in step S45. In step S46, it is checked whether there is an unprocessed candidate model. If there is an unprocessed candidate model (YES in S46), the process returns to S42. If there is no unprocessed candidate model (NO in S46), the process ends.

図6は、実施の形態1に係る電波源の識別システムの機種識別処理で被検査信号の立ち上がりまたは立ち下がりでの特徴量ごとに候補機種を絞る処理を説明するフロー図である。ステップS51で、被検査信号の立ち上がりが得られているかチェックする。被検査信号の立ち上がりが得られている場合(S51でYES)は、ステップS52で立ち上がり時の振幅特性、位相特性、周波数特性のいずれかである被検査特徴量を取得する。ステップS53で、候補機種から1個の機種を取り出す。ステップS54で、処理中機種の機種別特徴量範囲を機種別特徴量保存部13から取得する。ステップS55で、被検査特徴量が処理中機種の機種別特徴量範囲に含まれるか否かチェックする。含まれない場合(S55でNO)は、ステップS56で、処理中機種を候補機種から除く。ステップS57で、未処理の候補機種があるかチェックする。未処理の候補機種がある場合(S57でYES)は、S53に戻る。   FIG. 6 is a flowchart for explaining processing for narrowing down candidate models for each feature amount at the rise or fall of the signal under test in the model identification process of the radio wave source identification system according to the first embodiment. In step S51, it is checked whether the rising edge of the signal under test is obtained. If the rising edge of the inspected signal is obtained (YES in S51), in step S52, the inspected characteristic amount that is one of the amplitude characteristic, phase characteristic, and frequency characteristic at the time of rising is acquired. In step S53, one model is extracted from the candidate models. In step S54, the model-specific feature value range of the model being processed is acquired from the model-specific feature value storage unit 13. In step S55, it is checked whether the inspected feature value is included in the model-specific feature value range of the model being processed. If not included (NO in S55), the model being processed is excluded from the candidate models in step S56. In step S57, it is checked whether there is an unprocessed candidate model. If there is an unprocessed candidate model (YES in S57), the process returns to S53.

未処理の候補機種が無い場合(S57でNO)または被検査信号の立ち上がりが得られていない場合(S51でNO)は、ステップS58で、被検査信号の立ち下がりが得られているかチェックする。被検査信号の立ち下がりが得られている場合(S58でYES)は、ステップS59で立ち下がり時の振幅特性、位相特性、周波数特性のいずれかである被検査特徴量を取得する。ステップS60で、候補機種から1個の機種を取り出す。ステップS61で、処理中機種の機種別特徴量範囲を機種別特徴量保存部13から取得する。ステップS62で、被検査特徴量が処理中機種の機種別特徴量範囲に含まれるか否かチェックする。含まれない場合(S62でNO)は、ステップS63で、処理中機種を候補機種から除く。ステップS64で、未処理の候補機種があるかチェックする。未処理の候補機種がある場合(S64でYES)は、S59に戻る。未処理の候補機種が無い場合(S64でNO)は、終了する。   If there is no unprocessed candidate model (NO in S57) or the rising edge of the signal to be inspected is not obtained (NO in S51), it is checked in step S58 whether the falling edge of the signal to be inspected is obtained. If the fall of the signal to be inspected is obtained (YES in S58), the feature quantity to be inspected that is one of the amplitude characteristic, phase characteristic, and frequency characteristic at the time of falling is acquired in step S59. In step S60, one model is extracted from the candidate models. In step S 61, the model-specific feature amount range of the model being processed is acquired from the model-specific feature amount storage unit 13. In step S62, it is checked whether the inspected feature value is included in the model-specific feature value range of the model being processed. If not included (NO in S62), the model being processed is excluded from the candidate models in step S63. In step S64, it is checked whether there is an unprocessed candidate model. If there is an unprocessed candidate model (YES in S64), the process returns to S59. If there is no unprocessed candidate model (NO in S64), the process ends.

IQ軌跡を用いて候補種別を絞る処理について説明する。図7は、実施の形態1に係る電波源の識別システムの機種識別処理でIQ軌跡を用いて候補機種を絞る処理を説明するフロー図である。ステップS71で、被検査信号の立ち上がりが得られているかチェックする。被検査信号の立ち上がりが得られている場合(S71でYES)は、ステップS72で、IQ軌跡特性検出部11により立ち上がりの被検査信号からIQ軌跡を取得する。被検査信号から取得したIQ軌跡を被検査軌跡と呼ぶ。ステップS73で、候補機種から1個の機種を取り出す。ステップS74で、処理中機種のIQ軌跡を機種別特徴量保存部13から取得する。処理中機種のIQ軌跡を機種別軌跡と呼ぶ。ステップS75で、被検査軌跡と機種別軌跡の相関値を求める。ここで、相関が最大になるように被検査軌跡を拡大または縮小と回転させた時の相関値を、被検査軌跡と機種別軌跡の相関値とする。ステップS76で、相関値が閾値以上か否かをチェックする。相関値が閾値未満の場合(S76でNO)は、ステップS77で、処理中機種を候補機種から除く。ステップS78で、未処理の候補機種があるかチェックする。未処理の候補機種がある場合(S78でYES)は、S73に戻る。   Processing for narrowing candidate types using IQ trajectories will be described. FIG. 7 is a flowchart for explaining processing for narrowing down candidate models using the IQ trajectory in the model identification process of the radio wave source identification system according to the first embodiment. In step S71, it is checked whether the rising edge of the signal under test is obtained. If the rising edge of the signal to be inspected is obtained (YES in S71), the IQ locus is detected from the rising inspection signal by the IQ locus characteristic detecting unit 11 in step S72. The IQ trajectory acquired from the signal to be inspected is called an inspected trajectory. In step S73, one model is extracted from the candidate models. In step S74, the IQ trajectory of the model being processed is acquired from the model-specific feature amount storage unit 13. The IQ trajectory of the model being processed is called a model-specific trajectory. In step S75, a correlation value between the inspection trajectory and the model-specific trajectory is obtained. Here, the correlation value when the inspection trajectory is enlarged or reduced and rotated so that the correlation is maximized is set as the correlation value between the inspection trajectory and the model-specific trajectory. In step S76, it is checked whether or not the correlation value is greater than or equal to a threshold value. If the correlation value is less than the threshold (NO in S76), the model being processed is excluded from the candidate models in step S77. In step S78, it is checked whether there is an unprocessed candidate model. If there is an unprocessed candidate model (YES in S78), the process returns to S73.

未処理の候補機種が無い場合(S78でNO)または被検査信号の立ち上がりが得られていない場合(S71でNO)は、S79で、被検査信号の立ち下がりが得られているかチェックする。被検査信号の立ち下がりが得られている場合(S79でYES)は、ステップS80で、IQ軌跡特性検出部11により立ち下がりの被検査信号からIQ軌跡を取得する。ステップS81で、候補機種から1個の機種を取り出す。ステップS82で、処理中機種のIQ軌跡を機種別特徴量保存部13から取得する。処理中機種のIQ軌跡を機種別軌跡と呼ぶ。ステップS83で、被検査軌跡と機種別軌跡の相関値を求める。ここで、相関が最大になるように被検査軌跡を拡大または縮小と回転させた時の相関値を、被検査軌跡と機種別軌跡の相関値とする。ステップS84で、相関値が閾値以上か否かをチェックする。相関値が閾値未満の場合(S84でNO)は、ステップS85で、処理中機種を候補機種から除く。ステップS86で、未処理の候補機種があるかチェックする。未処理の候補機種がある場合(S86でYES)は、S81に戻る。未処理の候補機種が無い場合(S86でNO)は、終了する。   If there is no unprocessed candidate model (NO in S78) or the rising edge of the signal to be inspected is not obtained (NO in S71), it is checked in S79 whether the falling edge of the signal to be inspected is obtained. If the falling edge of the signal to be inspected is obtained (YES in S79), the IQ trajectory is acquired from the falling signal to be inspected by the IQ trajectory characteristic detecting unit 11 in step S80. In step S81, one model is extracted from the candidate models. In step S82, the IQ trajectory of the model being processed is acquired from the model-specific feature amount storage unit 13. The IQ trajectory of the model being processed is called a model-specific trajectory. In step S83, a correlation value between the inspection trajectory and the model-specific trajectory is obtained. Here, the correlation value when the inspection trajectory is enlarged or reduced and rotated so that the correlation is maximized is set as the correlation value between the inspection trajectory and the model-specific trajectory. In step S84, it is checked whether or not the correlation value is greater than or equal to a threshold value. If the correlation value is less than the threshold (NO in S84), the model being processed is excluded from the candidate models in step S85. In step S86, it is checked whether there is an unprocessed candidate model. If there is an unprocessed candidate model (YES in S86), the process returns to S81. If there is no unprocessed candidate model (NO in S86), the process ends.

図8を使用して、機種識別処理で機種を識別する例を説明する。機種としてAからHが存在する場合で説明する。周期性特性判定では機種A、機種Bおよび他の機種とでその機種別特徴量範囲が分離している。そのため、被検査信号の周期性の特徴量により、周期性特性判定で被検査信号が機種Aまたは機種Bと識別できる場合がある。振幅特性判定では、機種Cと他の機種とで機種別特徴量範囲が分離している。そのため、被検査信号の振幅特性の特徴量が機種Cの機種別特徴量範囲に入る場合は、機種Cと識別できる。   An example of identifying a model by the model identification process will be described with reference to FIG. A case where A to H exist as models will be described. In the periodicity characteristic determination, the model-specific feature range is separated for model A, model B, and other models. Therefore, the signal to be inspected may be identified from the model A or the model B in the periodicity characteristic determination based on the periodic characteristic amount of the signal to be inspected. In the amplitude characteristic determination, the model-specific feature range is separated between the model C and the other models. Therefore, when the feature value of the amplitude characteristic of the signal to be inspected falls within the model-specific feature value range of model C, it can be identified as model C.

位相特性判定では、機種Eと機種Fの機種別特徴量範囲が、それぞれだけの部分と重複する部分を有して、他の機種の機種別特徴量範囲と分離している。そのため、機種E、機種Fと識別できる場合、機種Eまたは機種Fと識別できる場合、機種D、機種Gおよび機種Hのどれかと識別できる場合がある。周波数特性判定により、位相特性判定では、機種Eまたは機種Fだったものが、機種Eまたは機種Fのどちらであるかが分る。IQ軌跡特性判定では、機種Gが識別できる場合がある。帯域幅特性判定により、機種Dと機種Hとを区別(識別)できる。   In the phase characteristic determination, the model-specific feature amount ranges of the model E and the model F have portions that overlap with the respective portions, and are separated from the model-specific feature amount ranges of other models. Therefore, when it can be identified as model E or model F, when it can be identified as model E or model F, it may be identified as any of model D, model G, or model H. According to the frequency characteristic determination, in the phase characteristic determination, it can be determined whether the model E or the model F is the model E or the model F. In the IQ trajectory characteristic determination, the model G may be identified. The model D and the model H can be distinguished (identified) by the bandwidth characteristic determination.

被検査信号を送信した電波源の機種が識別できると、識別した機種に属する個体を識別する。個体識別は機種識別と同様に実施する。図9は、実施の形態1に係る電波源の識別システムの個体識別処理を説明するフロー図である。機種を識別する場合の図4と同様な処理である。異なる点を説明する。すべての特徴量を考慮しても候補個体が1個以下にならない場合(S109でNO)は、ステップS104で、候補個体に残っている各個体の個体特徴量と被検査信号の特徴量との距離を求める。距離は、例えば、各特徴量の差の重み付け二乗和とする。各特徴量の差は、個体別特徴量範囲の幅で割ったものでもよい。ステップS105で、距離が最小である個体として識別する。すべての特徴量を考慮しても候補個体に残る個体は、各特徴量が個体別特徴範囲に入るので、距離が決められた上限以下であるといえる。   When the model of the radio wave source that transmitted the signal to be inspected can be identified, the individual belonging to the identified model is identified. Individual identification is performed in the same way as model identification. FIG. 9 is a flowchart for explaining individual identification processing of the radio wave source identification system according to the first embodiment. The process is the same as that in FIG. 4 for identifying the model. Explain the differences. If the number of candidate individuals does not become one or less even if all feature quantities are considered (NO in S109), in step S104, the individual feature quantity of each individual remaining in the candidate individual and the feature quantity of the signal to be inspected Find the distance. The distance is, for example, the weighted square sum of the difference between the feature amounts. The difference between the feature amounts may be divided by the width of the individual feature amount range. In step S105, the individual with the smallest distance is identified. Individuals that remain in the candidate individuals even if all feature quantities are taken into account can be said to be below the upper limit for which the distance is determined because each feature quantity falls within the individual feature range.

以下のように、被検査信号の特徴量との距離を求めることで、個体を識別してもよい。すべての個体との間でその個体の個体特徴量と被検査信号の特徴量との距離を求める。最小の距離を有する個体を最小距離個体と呼ぶ。求められた最小の距離が決められた上限以下である場合に、最小距離個体を被検査信号の個体と識別する。最小の距離が上限以下かどうかチェックする替わりに、各特徴量で最小距離個体の個体別特徴量範囲に被検査特徴量が入るかチェックしてもよい。   As described below, an individual may be identified by obtaining a distance from the feature amount of the signal to be inspected. The distance between the individual feature quantity of the individual and the feature quantity of the signal to be inspected is obtained for all individuals. An individual having the minimum distance is called a minimum distance individual. When the determined minimum distance is less than or equal to the determined upper limit, the minimum distance individual is identified as the individual of the signal under test. Instead of checking whether the minimum distance is less than or equal to the upper limit, it may be checked whether the feature quantity to be inspected falls within the individual feature quantity range of the individual with the minimum distance for each feature quantity.

なお、すべての特徴量を考慮しても候補個体が1個以下にならない場合は、個体識別不能としてもよい。候補個体が0個の場合(S106でNO)は、ステップS108で、被検査信号を新たな個体として個体別特徴量保存部14に保存する。候補個体が0個の場合(S106でNO)に、個体を識別できないとしてもよい。   In addition, when all the feature quantities are taken into consideration and the number of candidate individuals does not become one or less, individual identification may not be possible. If the number of candidate individuals is zero (NO in S106), the inspected signal is stored in the individual-specific feature storage unit 14 as a new individual in Step S108. If there are no candidate individuals (NO in S106), the individuals may not be identified.

図10は、実施の形態1に係る電波源の識別システムの個体識別処理で通信中での特徴量ごとに候補個体を絞る処理を説明するフロー図である。図11は、実施の形態1に係る電波源の識別システムの個体識別処理で被検査信号の立ち上がりまたは立ち下がりでの特徴量ごとに候補個体を絞る処理を説明するフロー図である。図12は、実施の形態1に係る電波源の識別システムの個体識別処理でIQ軌跡を用いて候補個体を絞る処理を説明するフロー図である。機種を識別する場合の図5から図7と同様な処理である。   FIG. 10 is a flowchart for explaining processing for narrowing down candidate individuals for each feature amount during communication in the individual identification processing of the radio wave source identification system according to the first embodiment. FIG. 11 is a flowchart for explaining processing for narrowing down candidate individuals for each feature amount at the rise or fall of the signal under test in the individual identification processing of the radio wave source identification system according to the first embodiment. FIG. 12 is a flowchart illustrating processing for narrowing down candidate individuals using IQ trajectories in the individual identification processing of the radio wave source identification system according to the first embodiment. The process is the same as that shown in FIGS. 5 to 7 for identifying the model.

周期性および帯域幅特性という通信中に得られる特徴量で機種および個体を識別するので、被検査信号の立ち上がりまたは立ち下がりが取得できない場合でも、被検査信号の機種および個体を識別できる。また、複数の特徴量を考慮して、機種および個体を精度よく識別できる。機種を先に識別するので、個体の識別まではできない場合も、機種を識別できる。機種を先に識別しているので、識別された機種に属する個体だけを識別対象にするので、効率的に個体を識別できる。
以上のことは、他の実施の形態にもあてはまる。
Since the model and the individual are identified by the feature values obtained during communication such as periodicity and bandwidth characteristics, the model and the individual of the signal to be inspected can be identified even when the rising or falling edge of the signal to be inspected cannot be acquired. In addition, the model and the individual can be accurately identified in consideration of a plurality of feature amounts. Since the model is identified first, the model can be identified even when the individual cannot be identified. Since the model is identified first, only individuals belonging to the identified model are targeted for identification, so that the individual can be identified efficiently.
The above also applies to other embodiments.

実施の形態2.
実施の形態2は、被検査信号の到来方位および変調方式を検査して、到来方位が合致する機種および個体を対象として、変調方式を考慮して選択した特徴量を使用して機種識別および個体識別を実施する場合である。
Embodiment 2. FIG.
In the second embodiment, the arrival direction and the modulation method of the signal to be inspected are inspected, and the model and the individual using the feature amount selected in consideration of the modulation method are targeted for the model and the individual whose arrival direction matches. This is a case where identification is performed.

図13は、この発明の実施の形態2に係る電波源の識別システムの構成図である。実施の形態1の場合の図1と異なる点を説明する。アンテナ1Aは、鉛直な方位軸の周りに回転可能なアンテナである。受信装置2Aは、被検査信号の到来方位を決める方位測定部17を有する。方位測定部17は、アンテナ1Aが回転するように制御する。また、アンテナ1Aが指向する方向において検波器6が検波する被検査信号の信号強度を調べ、信号強度が最大になる方位を被検査信号の到来方位として決める。到来方位は、幅を持つ方位の範囲として決められる。アンテナがアレイアンテナであり、電子的に指向する方位角を変化させて到来方位を決めてもよい。   FIG. 13 is a configuration diagram of a radio wave source identification system according to Embodiment 2 of the present invention. A different point from FIG. 1 in the case of Embodiment 1 is demonstrated. The antenna 1A is an antenna that can rotate around a vertical azimuth axis. The receiving apparatus 2A has an azimuth measuring unit 17 that determines the arrival azimuth of the signal to be inspected. The direction measuring unit 17 controls the antenna 1A to rotate. Further, the signal intensity of the signal under inspection detected by the detector 6 in the direction in which the antenna 1A is directed is examined, and the direction in which the signal intensity is maximized is determined as the arrival direction of the signal under inspection. The arrival direction is determined as a range of directions having a width. The antenna may be an array antenna, and the direction of arrival may be determined by changing the azimuth angle directed electronically.

電波源の識別装置3Aは、変調方式判定部18と方位限定検索部19も有する。変調方式判定部18は、被検査信号の変調方式を判定する。方位限定検索部19は、到来方位を考慮して機種を機種別特徴量保存部13Aから検索し、到来方位を考慮して個体を個体別特徴量保存部14Aから検索する。   The radio wave source identification device 3 </ b> A also includes a modulation method determination unit 18 and an orientation limitation search unit 19. The modulation method determination unit 18 determines the modulation method of the signal under inspection. The azimuth limited search unit 19 searches for the model from the model-specific feature storage unit 13A in consideration of the arrival azimuth, and searches the individual feature storage unit 14A in consideration of the arrival azimuth.

変調方式判定部18は、被検査信号に対して時間波形分析、周波数スペクトラム分析およびコンスタレーション分析などを実施して、変調方式を判定する。例えば、時間波形分析を適用する場合には、PSK(Phase Shift Keying)系のデジタル変調では位相変化時に振幅がゼロとなる点(ゼロクロス)が存在し、時間波形ではゼロクロス通過により振幅の落ち込み(ヌル)が発生する特徴が現れる。時間波形にヌルが存在する場合に、PSK変調と判断する。また、例えば、周波数スペクトラム分析を適用する場合には、周波数スペクトラムでのピークの数によりFSK(Frequency Shift Keying)系のデジタル変調方式で何値のFSK化を判断する。これは、FSK系のデジタル変調は伝送するビットに合わせて搬送波の周波数を変化させるため、周波数スペクトラム上で複数のピークが生じることに基づく。コンスタレーション分析を適用する場合には、IQ平面におけるシンボル点の集まり具合で、変調方式を判断する。例えば、QAM(Quadrature Amplitude Modulation)系の変調方式では、コンスタレーション分析によりシンボル点の数の推定が行える。   The modulation method determination unit 18 performs time waveform analysis, frequency spectrum analysis, constellation analysis, and the like on the signal to be inspected to determine the modulation method. For example, when applying time waveform analysis, there is a point where the amplitude becomes zero (zero cross) at the time of phase change in digital modulation of PSK (Phase Shift Keying) system, and in the time waveform, the amplitude drops due to zero crossing (null) ) Appears. When null exists in the time waveform, it is determined as PSK modulation. For example, when frequency spectrum analysis is applied, the FSK (Frequency Shift Keying) type digital modulation scheme is used to determine the value of FSK based on the number of peaks in the frequency spectrum. This is based on the fact that a plurality of peaks occur on the frequency spectrum because the FSK digital modulation changes the frequency of the carrier wave in accordance with the bit to be transmitted. When constellation analysis is applied, the modulation method is determined based on the collection of symbol points on the IQ plane. For example, in a QAM (Quadrature Amplitude Modulation) modulation method, the number of symbol points can be estimated by constellation analysis.

変調方式判定部18は、振幅変調、位相変調、周波数変調、不明のどれかを出力する。振幅変調、位相変調、周波数変調のそれぞれについて、その変調方式である可能性があるか否かを出力するようにしてもよい。   The modulation method determination unit 18 outputs one of amplitude modulation, phase modulation, frequency modulation, and unknown. For each of amplitude modulation, phase modulation, and frequency modulation, whether or not there is a possibility of the modulation method may be output.

個体別特徴量保存部14Aには、各個体の位置が固定かどうかと、固定である場合の位置も保存される。機種別特徴量保存部13Aには、位置が固定である個体の有無と、位置が固定な個体が属する機種である場合には、それらの個体の位置を参照できるデータが保存される。   Whether or not the position of each individual is fixed and the position when it is fixed are also stored in the individual feature amount storage unit 14A. In the model-specific feature amount storage unit 13A, data indicating whether or not there is an individual whose position is fixed and a model to which the individual whose position is fixed belong can be stored.

機種識別部15Aは、被検査信号の到来方位および変調方式を考慮して、被検査信号の機種を識別する。個体識別部16Aは、被検査信号の到来方位および変調方式を考慮して、被検査信号の機種を識別する。   The model identifying unit 15A identifies the model of the signal to be inspected in consideration of the arrival direction of the signal to be inspected and the modulation method. The individual identifying unit 16A identifies the model of the signal to be inspected in consideration of the arrival direction of the signal to be inspected and the modulation method.

動作を説明する。図14は、実施の形態2に係る電波源の識別システムの概略処理を説明するフロー図である。実施の形態1の場合の図3と異なる点を説明する。S01で抽出された被検査信号に対して、ステップS04で、方位測定部17により被検査信号の到来方位を検出する。ステップS05で、変調方式判定部18により被検査信号の変調方式を検出する。   The operation will be described. FIG. 14 is a flowchart illustrating a schematic process of the radio wave source identification system according to the second embodiment. A different point from FIG. 3 in the case of Embodiment 1 is demonstrated. In step S04, the direction measuring unit 17 detects the arrival direction of the signal to be inspected with respect to the signal to be inspected extracted in S01. In step S05, the modulation system determination unit 18 detects the modulation system of the signal to be inspected.

ステップS02Aで、被検査信号の到来方位および変調方式を考慮して、機種識別部15Aが被検査信号の機種を識別する。図15は、実施の形態2に係る電波源の識別システムの機種識別処理を説明するフロー図である。実施の形態1の場合の図4と異なる点を説明する。   In step S02A, the model identifying unit 15A identifies the model of the signal to be inspected in consideration of the arrival direction of the signal to be inspected and the modulation method. FIG. 15 is a flowchart for explaining model identification processing of the radio wave source identification system according to the second embodiment. A different point from FIG. 4 in the case of Embodiment 1 is demonstrated.

ステップS11Aで、被検査信号の到来方位を考慮して被検査信号の候補機種を決める。具体的には、個体の位置が固定であってアンテナ1Aから見て到来方位の範囲内に存在する個体、または位置が固定でない個体が属する機種が候補機種である。アンテナ1Aから見て到来方位にある位置を到来方位と整合する位置と呼ぶ。到来方位と整合する位置に存在する個体を方位整合個体と呼ぶ。方位限定検索部19が、到来方位に整合する機種を検索する。到来方位に整合する機種は、方位整合個体を有する機種である。   In step S11A, a candidate model of the signal to be inspected is determined in consideration of the arrival direction of the signal to be inspected. Specifically, a candidate model is an individual whose position is fixed and exists within the range of the arrival direction as viewed from the antenna 1A, or a model to which an individual whose position is not fixed belongs. A position in the arrival direction as viewed from the antenna 1A is referred to as a position that matches the arrival direction. An individual that exists at a position that matches the arrival direction is called an orientation-matched individual. The azimuth limited search unit 19 searches for models that match the incoming azimuth. A model that matches the direction of arrival is a model that has an orientation-matched individual.

S13でNOの場合に、ステップS28で、変調方式が振幅変調または不明かチェックする。変調方式が振幅変調または不明である場合(S28でYES)は、S14で振幅特性についての特徴量を使用して機種を識別する。S28で、変調方式が振幅変調である可能性があるか否かをチェックしてもよい。   If NO in S13, it is checked in step S28 whether the modulation method is amplitude modulation or unknown. If the modulation method is amplitude modulation or unknown (YES in S28), the model is identified using the feature value for the amplitude characteristic in S14. In S28, it may be checked whether or not the modulation method may be amplitude modulation.

S28またはS15でNOの場合に、ステップS29で、変調方式が位相変調または不明かチェックする。変調方式が位相変調または不明である場合(S29でYES)は、S16で位相特性についての特徴量を使用して機種を識別する。S29で、変調方式が位相変調である可能性があるか否かをチェックしてもよい。   If NO in S28 or S15, it is checked in step S29 whether the modulation method is phase modulation or unknown. If the modulation method is phase modulation or unknown (YES in S29), the model is identified using the feature amount for the phase characteristic in S16. In S29, it may be checked whether or not the modulation method may be phase modulation.

S29またはS17でNOの場合に、ステップS30で、変調方式が周波数変調または不明かチェックする。変調方式が周波数変調または不明である場合(S30でYES)は、S18で周波数特性についての特徴量を使用して機種を識別する。S30で、変調方式が周波数変調である可能性があるか否かをチェックしてもよい。   If NO in S29 or S17, it is checked in step S30 whether the modulation method is frequency modulation or unknown. If the modulation method is frequency modulation or unknown (YES in S30), the model is identified using the feature amount for the frequency characteristic in S18. In S30, it may be checked whether or not the modulation method may be frequency modulation.

ここでは、変調方式により振幅変調、位相変調、周波数変調のどれかであることが明確な場合は、それぞれ振幅特性、位相特性、周波数特性だけを使用して機種を識別した。検出した変調方式によらず振幅特性、位相特性、周波数特性で機種を識別するようにし、変調方式により振幅特性、位相特性、周波数特性という3個の特徴量を処理する順番を変更するようにしてもよい。   Here, when it is clear that the modulation method is one of amplitude modulation, phase modulation, and frequency modulation, the model is identified using only the amplitude characteristics, phase characteristics, and frequency characteristics, respectively. Regardless of the detected modulation method, the model is identified by the amplitude characteristic, phase characteristic, and frequency characteristic, and the order of processing the three feature quantities of the amplitude characteristic, the phase characteristic, and the frequency characteristic is changed by the modulation method. Also good.

ステップS03Aで、被検査信号の到来方位および変調方式を考慮して、個体識別部16Aが被検査信号の個体を識別する。図16は、実施の形態2に係る電波源の識別システムの個体識別処理を説明するフロー図である。実施の形態1の場合の図9と異なる点を説明する。   In step S03A, the individual identifying unit 16A identifies the individual of the signal under test in consideration of the arrival direction and modulation scheme of the signal under inspection. FIG. 16 is a flowchart for explaining individual identification processing of the radio wave source identification system according to the second embodiment. A different point from FIG. 9 in the case of Embodiment 1 is demonstrated.

ステップS91Aで、被検査信号の到来方位を考慮して被検査信号の候補個体を決める。具体的には、個体の位置が固定であってアンテナ1Aから見て到来方位の範囲内に存在する個体、または位置が固定でない個体である。方位限定検索部19が、到来方位に整合する個体を検索する。   In step S91A, a candidate individual of the signal to be inspected is determined in consideration of the arrival direction of the signal to be inspected. Specifically, it is an individual whose position is fixed and exists within the range of the arrival direction as viewed from the antenna 1A, or an individual whose position is not fixed. The azimuth limited search unit 19 searches for an individual that matches the arrival azimuth.

S93でYESの場合に、ステップS109で、変調方式が振幅変調または不明かチェックする。変調方式が振幅変調または不明である場合(S109でYES)は、S94で振幅特性についての特徴量を使用して個体を識別する。   If YES in step S93, it is checked in step S109 if the modulation method is amplitude modulation or unknown. If the modulation method is amplitude modulation or unknown (YES in S109), an individual is identified using the feature amount for the amplitude characteristic in S94.

S109でNOまたはS95でYESの場合に、ステップS110で、変調方式が位相変調または不明かチェックする。変調方式が位相変調または不明である場合(S110でYES)は、S96で位相特性についての特徴量を使用して個体を識別する。   If NO in S109 or YES in S95, it is checked in step S110 whether the modulation method is phase modulation or unknown. If the modulation method is phase modulation or unknown (YES in S110), the individual is identified using the feature amount for the phase characteristic in S96.

S110でNOまたはS97でYESの場合に、ステップS111で、変調方式が周波数変調または不明かチェックする。変調方式が周波数変調または不明である場合(S111でYES)は、S98で周波数特性についての特徴量を使用して個体を識別する。   If NO in S110 or YES in S97, it is checked in step S111 whether the modulation method is frequency modulation or unknown. If the modulation method is frequency modulation or unknown (YES in S111), an individual is identified using the feature amount for the frequency characteristic in S98.

この実施の形態では、電波源の到来方位を考慮して、到来方位が整合する位置を有する機種および個体だけを対象とするので、処理を高速にかつより高精度に識別できる。変調方式を考慮して被検査信号の変調方式に応じて検査する特徴量や検査する順番を変更できるので、処理を高速にまたはより高精度に識別できる。   In this embodiment, taking into account the arrival direction of the radio wave source, only models and individuals having a position where the arrival direction matches are targeted, so that the processing can be identified at high speed and with higher accuracy. Considering the modulation method, the feature quantity to be inspected and the inspection order can be changed according to the modulation method of the signal to be inspected, so that the processing can be identified at high speed or with higher accuracy.

実施の形態3.
実施の形態3は、被検査信号の到来方位を検査して、到来方位が合致する機種および個体を対象として機種識別および個体識別を実施する場合である。
Embodiment 3 FIG.
The third embodiment is a case in which the arrival direction of a signal to be inspected is inspected, and model identification and individual identification are performed for models and individuals whose arrival directions match.

図17は、この発明の実施の形態3に係る電波源の識別システムの構成図である。実施の形態2の場合と同様な、アンテナ1A、方位測定部17、個体別特徴量保存部14Aおよび機種別特徴量保存部13Aを有する。機種識別部15Bは、被検査信号の到来方位を考慮して、被検査信号の機種を識別する。個体識別部15Bは、被検査信号の到来方位を考慮して、被検査信号の機種を識別する。   FIG. 17 is a configuration diagram of a radio wave source identification system according to Embodiment 3 of the present invention. As in the case of the second embodiment, the antenna 1A, the azimuth measuring unit 17, the individual feature amount storing unit 14A, and the model-specific feature amount storing unit 13A are provided. The model identifying unit 15B identifies the model of the signal to be inspected in consideration of the arrival direction of the signal to be inspected. The individual identification unit 15B identifies the model of the signal to be inspected in consideration of the arrival direction of the signal to be inspected.

動作を説明する。図18は、実施の形態3に係る電波源の識別システムの概略処理を説明するフロー図である。実施の形態2の場合の図14と同様に、ステップS04で、方位測定部17により被検査信号の到来方位を判定する。   The operation will be described. FIG. 18 is a flowchart for explaining schematic processing of the radio wave source identification system according to the third embodiment. Similarly to FIG. 14 in the case of the second embodiment, the azimuth measuring unit 17 determines the arrival azimuth of the signal to be inspected in step S04.

ステップS02Bで、被検査信号の到来方位を考慮して、機種識別部15Bが被検査信号の機種を識別する。図19は、実施の形態3に係る電波源の識別システムの機種識別処理を説明するフロー図である。実施の形態2の場合と同様に、ステップS11Aで、被検査信号の到来方位を考慮して被検査信号の候補機種を決める。他は実施の形態1の場合の図4と同様である。   In step S02B, the model identifying unit 15B identifies the model of the signal under inspection in consideration of the arrival direction of the signal under inspection. FIG. 19 is a flowchart for explaining model identification processing of the radio wave source identification system according to the third embodiment. As in the case of the second embodiment, in step S11A, a candidate model of the signal to be inspected is determined in consideration of the arrival direction of the signal to be inspected. Others are the same as those in FIG. 4 in the first embodiment.

ステップS03Bで、被検査信号の到来方位を考慮して、個体識別部16Bが被検査信号の個体を識別する。図20は、実施の形態2に係る電波源の識別システムの個体識別処理を説明するフロー図である。ステップS91Aで、被検査信号の到来方位を考慮して被検査信号の候補個体を決める。他は実施の形態1の場合の図9と同様である。   In step S03B, the individual identifying unit 16B identifies the individual of the signal under test in consideration of the arrival direction of the signal under inspection. FIG. 20 is a flowchart for explaining individual identification processing of the radio wave source identification system according to the second embodiment. In step S91A, a candidate individual of the signal to be inspected is determined in consideration of the arrival direction of the signal to be inspected. Others are the same as those in FIG. 9 in the first embodiment.

電波源の到来方位を考慮して、到来方位が整合する位置を有する機種および個体だけを対象とするので、処理を高速にかつより高精度に識別できる。   Considering the arrival direction of the radio wave source, only the model and the individual having the position where the arrival direction matches are targeted, so that the processing can be identified at high speed and with higher accuracy.

実施の形態4.
実施の形態4は、被検査信号の変調方式を検査して、変調方式を考慮して選択した特徴量を使用して機種識別および個体識別を実施する場合である。
Embodiment 4 FIG.
The fourth embodiment is a case where the modulation scheme of a signal to be inspected is inspected, and model identification and individual identification are performed using feature quantities selected in consideration of the modulation scheme.

図21は、この発明の実施の形態4に係る電波源の識別システムの構成図である。実施の形態1の場合の図1と異なる点を説明する。実施の形態2と同様な変調方式判定部18を有する。機種識別部15Cは、被検査信号の変調方式を考慮して被検査信号の機種を識別する。個体識別部16Cは、被検査信号の変調方式を考慮して被検査信号の機種を識別する。   FIG. 21 is a configuration diagram of a radio wave source identification system according to Embodiment 4 of the present invention. A different point from FIG. 1 in the case of Embodiment 1 is demonstrated. A modulation scheme determination unit 18 similar to that of the second embodiment is included. The model identifying unit 15C identifies the model of the signal to be inspected in consideration of the modulation method of the signal to be inspected. The individual identification unit 16C identifies the model of the signal to be inspected in consideration of the modulation method of the signal to be inspected.

動作を説明する。図22は、実施の形態4に係る電波源の識別システムの概略処理を説明するフロー図である。実施の形態2の場合と同様に、ステップS05で、変調方式判定部18により被検査信号の変調方式を判定する。   The operation will be described. FIG. 22 is a flowchart illustrating a schematic process of the radio wave source identification system according to the fourth embodiment. As in the case of the second embodiment, in step S05, the modulation method determination unit 18 determines the modulation method of the signal to be inspected.

ステップS02Cで、被検査信号の変調方式を考慮して、機種識別部15Cが被検査信号の機種を識別する。図23は、実施の形態4に係る電波源の識別システムの機種識別処理を説明するフロー図である。実施の形態2の場合の図15と異なる点を説明する。被検査信号の到来方位を求めないので、ステップS11で、すべての機種を候補機種とする。   In step S02C, the model identifying unit 15C identifies the model of the signal to be inspected in consideration of the modulation method of the signal to be inspected. FIG. 23 is a flowchart for explaining model identification processing of the radio wave source identification system according to the fourth embodiment. A different point from FIG. 15 in the case of Embodiment 2 is demonstrated. Since the arrival direction of the signal to be inspected is not obtained, all models are set as candidate models in step S11.

ステップS03Cで、被検査信号の変調方式を考慮して、個体識別部16Cが被検査信号の個体を識別する。図24は、実施の形態4に係る電波源の識別システムの個体識別処理を説明するフロー図である。実施の形態2の場合の図16と異なる点を説明する。被検査信号の到来方位を求めないので、ステップS91で、すべての個体を候補個体とする。   In step S03C, the individual identifying unit 16C identifies the individual of the signal under test in consideration of the modulation scheme of the signal under test. FIG. 24 is a flowchart for explaining individual identification processing of the radio wave source identification system according to the fourth embodiment. A different point from FIG. 16 in the case of Embodiment 2 is demonstrated. Since the arrival direction of the signal to be inspected is not obtained, all individuals are set as candidate individuals in step S91.

変調方式を考慮して被検査信号の変調方式に応じて検査する特徴量や検査する順番を変更できるので、処理をより高速またはより高精度に識別できる。   Considering the modulation method, the feature quantity to be inspected and the order of inspection can be changed according to the modulation method of the signal to be inspected, so that the processing can be identified at higher speed or higher accuracy.

本発明はその発明の精神の範囲内において各実施の形態の自由な組み合わせ、あるいは各実施の形態の変形や省略が可能である。   The present invention can be freely combined with each other, or can be modified or omitted within the spirit of the invention.

50、50A、50B、50C 電波源の識別システム
1、1A アンテナ
2、2A 受信装置
3、3A、3B、3C 電波源の識別装置
4 周波数変換器
5 A/D変換器
6 検波器
7 周期性検出部
8 振幅特性検出部
9 位相特性検出部
10 周波数特性検出部
11 IQ軌跡特性検出部
12 帯域幅特性検出部
13、13A 機種別特徴量保存部
14、14A 個体別特徴量保存部
15、15A、15B、15C 機種識別部
16、16A、16B、16C 個体識別部
17 方位測定部
18 変調方式判定部
19 方位限定検索部
50, 50A, 50B, 50C Radio wave source identification system 1, 1A antenna 2, 2A receiver 3, 3A, 3B, 3C radio wave source identification device 4 Frequency converter 5 A / D converter 6 Detector 7 Periodic detection Unit 8 amplitude characteristic detection unit 9 phase characteristic detection unit 10 frequency characteristic detection unit 11 IQ trajectory characteristic detection unit 12 bandwidth characteristic detection unit 13, 13A model-specific feature amount storage unit 14, 14A individual feature amount storage unit 15, 15A, 15B, 15C Model identification unit 16, 16A, 16B, 16C Individual identification unit 17 Orientation measurement unit 18 Modulation method determination unit 19 Orientation limited search unit

Claims (14)

アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、
ある機種の電波源が送信する電波が受信されて得られる前記受信信号の複数の前記特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を前記電波源の前記機種ごとに保存する機種別特徴量保存部と、
前記被検査信号から抽出した前記特徴量である被検査特徴量と前記機種別特徴量保存部に保存された前記機種別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記機種別特徴量範囲に入る前記機種を、前記被検査信号を送信した前記電波源の前記機種として識別する機種識別部とを備え
前記特徴量として、前記被検査信号と遅延時間だけ遅延させた前記被検査信号との相関値が極大値をとる前記遅延時間を使用する、電波源の識別装置。
A feature quantity extraction unit that extracts a plurality of feature quantities of a signal under test detected from a discretized received signal generated from a radio wave received by an antenna, including at least one feature quantity during communication;
A model-specific feature range that is a range that can be taken by each of the plurality of feature quantities of the received signal obtained by receiving radio waves transmitted from a radio wave source of a certain model. A feature storage unit;
By comparing the inspected feature quantity, which is the feature quantity extracted from the inspected signal, with the model-specific feature quantity range stored in the model-specific feature quantity storage unit, each of the plurality of inspected feature quantities is A model identification unit that identifies the model that falls within the model-specific feature amount range as the model of the radio wave source that has transmitted the signal to be inspected ;
A radio wave source identification device using, as the feature amount, the delay time at which a correlation value between the signal to be inspected and the signal to be inspected delayed by a delay time takes a maximum value .
アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、
ある機種の電波源が送信する電波が受信されて得られる前記受信信号の複数の前記特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を前記電波源の前記機種ごとに保存する機種別特徴量保存部と、
前記被検査信号から抽出した前記特徴量である被検査特徴量と前記機種別特徴量保存部に保存された前記機種別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記機種別特徴量範囲に入る前記機種を、前記被検査信号を送信した前記電波源の前記機種として識別する機種識別部とを備え、
前記特徴量として、前記被検査信号の帯域幅を使用する、電波源の識別装置。
A feature quantity extraction unit that extracts a plurality of feature quantities of a signal under test detected from a discretized received signal generated from a radio wave received by an antenna, including at least one feature quantity during communication;
A model-specific feature range that is a range that can be taken by each of the plurality of feature quantities of the received signal obtained by receiving radio waves transmitted from a radio wave source of a certain model. A feature storage unit;
By comparing the inspected feature quantity, which is the feature quantity extracted from the inspected signal, with the model-specific feature quantity range stored in the model-specific feature quantity storage unit, each of the plurality of inspected feature quantities is A model identification unit that identifies the model that falls within the model-specific feature amount range as the model of the radio wave source that has transmitted the signal to be inspected;
A radio wave source identification device that uses the bandwidth of the signal to be inspected as the feature amount .
アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、
ある機種の電波源が送信する電波が受信されて得られる前記受信信号の複数の前記特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を前記電波源の前記機種ごとに保存する機種別特徴量保存部と、
前記被検査信号から抽出した前記特徴量である被検査特徴量と前記機種別特徴量保存部に保存された前記機種別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記機種別特徴量範囲に入る前記機種を、前記被検査信号を送信した前記電波源の前記機種として識別する機種識別部と、
それぞれの前記機種に属する前記電波源であるそれぞれの個体が取りうる範囲である個体別特徴量範囲を前記個体ごとに保存する個体別特徴量保存部と、
前記被検査特徴量と前記個体別特徴量保存部に保存された前記個体別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記個体別特徴量範囲に入る前記個体を、前記被検査信号を送信した前記電波源であるとして識別する個体識別部とを備えた電波源の識別装置。
A feature quantity extraction unit that extracts a plurality of feature quantities of a signal under test detected from a discretized received signal generated from a radio wave received by an antenna, including at least one feature quantity during communication;
A model-specific feature range that is a range that can be taken by each of the plurality of feature quantities of the received signal obtained by receiving radio waves transmitted from a radio wave source of a certain model. A feature storage unit;
By comparing the inspected feature quantity, which is the feature quantity extracted from the inspected signal, with the model-specific feature quantity range stored in the model-specific feature quantity storage unit, each of the plurality of inspected feature quantities is A model identifying unit that identifies the model that falls within the model-specific feature amount range as the model of the radio wave source that has transmitted the signal to be inspected;
An individual feature storage unit that stores, for each individual, an individual feature amount range that can be taken by each individual that is the radio wave source belonging to each of the models;
The individual feature amount range is compared with the individual feature amount range stored in the individual feature amount storage unit, and each of the plurality of inspected feature amounts falls within the individual feature amount range. A radio wave source identification device comprising: an individual identification unit that identifies the radio wave source that has transmitted the signal to be inspected .
前記個体別特徴量保存部が、前記個体の前記特徴量である個体特徴量を前記個体ごとに保存し、
前記個体識別部が、前記機種識別部が識別した前記機種に属する前記個体の中で前記個体特徴量と前記被検査特徴量との距離が最も小さい前記個体を、前記被検査信号を送信した前記電波源であるとして識別する、請求項3に記載の電波源の識別装置。
The individual feature amount storage unit stores the individual feature amount that is the feature amount of the individual for each individual,
The individual identification unit has transmitted the signal to be inspected, the individual having the smallest distance between the individual feature quantity and the feature quantity to be inspected among the individuals belonging to the model identified by the model identification unit. 4. The radio wave source identification device according to claim 3, wherein the radio wave source is identified as being a radio wave source.
アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、
ある機種の電波源が送信する電波が受信されて得られる前記受信信号の複数の前記特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を前記電波源の前記機種ごとに保存する機種別特徴量保存部と、
前記被検査信号から抽出した前記特徴量である被検査特徴量と前記機種別特徴量保存部に保存された前記機種別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記機種別特徴量範囲に入る前記機種を、前記被検査信号を送信した前記電波源の前記機種として識別する機種識別部と、
それぞれの前記機種に属する前記電波源であるそれぞれの個体の前記特徴量である個体特徴量を前記個体ごとに保存する個体別特徴量保存部と、
前記機種識別部が識別した前記機種に属する前記個体の前記個体特徴量の中で前記被検査特徴量との距離が決められた上限以下で最も小さい前記個体を、前記被検査信号を送信した前記電波源であるとして識別する個体識別部とを備えた電波源の識別装置。
A feature quantity extraction unit that extracts a plurality of feature quantities of a signal under test detected from a discretized received signal generated from a radio wave received by an antenna, including at least one feature quantity during communication;
A model-specific feature range that is a range that can be taken by each of the plurality of feature quantities of the received signal obtained by receiving radio waves transmitted from a radio wave source of a certain model. A feature storage unit;
By comparing the inspected feature quantity, which is the feature quantity extracted from the inspected signal, with the model-specific feature quantity range stored in the model-specific feature quantity storage unit, each of the plurality of inspected feature quantities is A model identifying unit that identifies the model that falls within the model-specific feature amount range as the model of the radio wave source that has transmitted the signal to be inspected;
An individual feature storage unit that stores, for each individual, an individual feature that is the individual feature of each individual that is the radio wave source belonging to each of the models;
The individual that is the smallest of the individual feature amounts of the individual belonging to the model identified by the model identification unit is less than or equal to a predetermined upper limit of the distance to the inspected feature amount is transmitted to the inspected signal. 1. A radio wave source identification device comprising an individual identification unit that identifies a radio wave source.
前記個体別特徴量保存部が前記個体の位置および位置が固定か否かを保存し
前記機種識別部が、前記位置が固定であって前記被検査信号が到来する方位である到来方位と整合する前記個体または前記位置が固定でない前記個体である方位整合個体が属する前記機種を対象にして処理し、
前記個体識別部は、前記機種識別部が識別した前記機種に属する前記方位整合個体を対象にして処理する、請求項3から請求項5のいずれか1項に記載の電波源の識別装置。
The individual feature amount storage unit stores whether or not the position and position of the individual are fixed.
The model identification unit targets the model to which the individual that matches the arrival azimuth that is the azimuth where the signal to be inspected is fixed or the azimuth matching individual that is the individual whose position is not fixed. Process
The radio wave source identification device according to any one of claims 3 to 5, wherein the individual identification unit processes the orientation matching individual belonging to the model identified by the model identification unit.
アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、
ある機種の電波源が送信する電波が受信されて得られる前記受信信号の複数の前記特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を前記電波源の前記機種ごとに保存する機種別特徴量保存部と、
前記被検査信号から抽出した前記特徴量である被検査特徴量と前記機種別特徴量保存部に保存された前記機種別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記機種別特徴量範囲に入る前記機種を、前記被検査信号を送信した前記電波源の前記機種として識別する機種識別部とを備え、
前記特徴量として、送信開始時の決められた期間の前記被検査信号の振幅、位相、周波数のいずれか少なくとも1つについて標準偏差、歪度および尖度の少なくとも1つを使用する、電波源の識別装置。
A feature quantity extraction unit that extracts a plurality of feature quantities of a signal under test detected from a discretized received signal generated from a radio wave received by an antenna, including at least one feature quantity during communication;
A model-specific feature range that is a range that can be taken by each of the plurality of feature quantities of the received signal obtained by receiving radio waves transmitted from a radio wave source of a certain model. A feature storage unit;
By comparing the inspected feature quantity, which is the feature quantity extracted from the inspected signal, with the model-specific feature quantity range stored in the model-specific feature quantity storage unit, each of the plurality of inspected feature quantities is A model identification unit that identifies the model that falls within the model-specific feature amount range as the model of the radio wave source that has transmitted the signal to be inspected;
A radio wave source that uses at least one of standard deviation, skewness, and kurtosis for at least one of amplitude, phase, and frequency of the signal to be inspected for a predetermined period at the start of transmission as the feature amount Identification device.
アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、
ある機種の電波源が送信する電波が受信されて得られる前記受信信号の複数の前記特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を前記電波源の前記機種ごとに保存する機種別特徴量保存部と、
前記被検査信号から抽出した前記特徴量である被検査特徴量と前記機種別特徴量保存部に保存された前記機種別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記機種別特徴量範囲に入る前記機種を、前記被検査信号を送信した前記電波源の前記機種として識別する機種識別部とを備え、
前記特徴量として、送信終了時の決められた期間の前記被検査信号の振幅、位相、周波数のいずれか少なくとも1つについて標準偏差、歪度および尖度の少なくとも1つを使用する、電波源の識別装置。
A feature quantity extraction unit that extracts a plurality of feature quantities of a signal under test detected from a discretized received signal generated from a radio wave received by an antenna, including at least one feature quantity during communication;
A model-specific feature range that is a range that can be taken by each of the plurality of feature quantities of the received signal obtained by receiving radio waves transmitted from a radio wave source of a certain model. A feature storage unit;
By comparing the inspected feature quantity, which is the feature quantity extracted from the inspected signal, with the model-specific feature quantity range stored in the model-specific feature quantity storage unit, each of the plurality of inspected feature quantities is A model identification unit that identifies the model that falls within the model-specific feature amount range as the model of the radio wave source that has transmitted the signal to be inspected;
A radio wave source that uses at least one of standard deviation, skewness, and kurtosis for at least one of amplitude, phase, and frequency of the signal to be inspected for a predetermined period at the end of transmission as the feature amount Identification device.
アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、
ある機種の電波源が送信する電波が受信されて得られる前記受信信号の複数の前記特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を前記電波源の前記機種ごとに保存する機種別特徴量保存部と、
前記被検査信号から抽出した前記特徴量である被検査特徴量と前記機種別特徴量保存部に保存された前記機種別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記機種別特徴量範囲に入る前記機種を、前記被検査信号を送信した前記電波源の前記機種として識別する機種識別部とを備え、
前記被検査信号は直交検波され、
前記特徴量として、送信開始時の決められた期間の前記被検査信号の複素平面での軌跡である被検査軌跡と、前記機種別特徴量保存部に保存された前記機種ごとの複素平面での軌跡である機種別軌跡との類似度を使用する、電波源の識別装置。
A feature quantity extraction unit that extracts a plurality of feature quantities of a signal under test detected from a discretized received signal generated from a radio wave received by an antenna, including at least one feature quantity during communication;
A model-specific feature range that is a range that can be taken by each of the plurality of feature quantities of the received signal obtained by receiving radio waves transmitted from a radio wave source of a certain model. A feature storage unit;
By comparing the inspected feature quantity, which is the feature quantity extracted from the inspected signal, with the model-specific feature quantity range stored in the model-specific feature quantity storage unit, each of the plurality of inspected feature quantities is A model identification unit that identifies the model that falls within the model-specific feature amount range as the model of the radio wave source that has transmitted the signal to be inspected;
The inspected signal is quadrature detected,
As the feature amount, an inspected trajectory that is a trajectory in the complex plane of the signal to be inspected for a predetermined period at the start of transmission, and a complex plane for each model stored in the model-specific feature amount storage unit. A radio wave source identification device that uses a similarity with a trajectory of a model, which is a trajectory .
前記被検査信号は直交検波され、
前記特徴量として、送信開始時の決められた期間の前記被検査信号の複素平面での軌跡である被検査軌跡と、前記個体別特徴量保存部に保存された前記個体ごとの複素平面での軌跡である個体別軌跡との類似度を使用する、請求項3から請求項6のいずれか1項に記載の電波源の識別装置。
The inspected signal is quadrature detected,
As the feature amount, an inspected trajectory that is a trajectory in the complex plane of the signal to be inspected for a predetermined period at the start of transmission, and a complex plane for each individual stored in the individual feature amount storage unit. The radio wave source identification device according to any one of claims 3 to 6, wherein a similarity with an individual-specific trajectory that is a trajectory is used .
アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、
ある機種の電波源が送信する電波が受信されて得られる前記受信信号の複数の前記特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を前記電波源の前記機種ごとに保存する機種別特徴量保存部と、
前記被検査信号から抽出した前記特徴量である被検査特徴量と前記機種別特徴量保存部に保存された前記機種別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記機種別特徴量範囲に入る前記機種を、前記被検査信号を送信した前記電波源の前記機種として識別する機種識別部と、
前記被検査信号の変調方式を判定する変調方式判定部とを備え、
前記変調方式判定部が判定した変調方式に応じて、前記機種識別部が使用する前記特徴量を決める、電波源の識別装置。
A feature quantity extraction unit that extracts a plurality of feature quantities of a signal under test detected from a discretized received signal generated from a radio wave received by an antenna, including at least one feature quantity during communication;
A model-specific feature range that is a range that can be taken by each of the plurality of feature quantities of the received signal obtained by receiving radio waves transmitted from a radio wave source of a certain model. A feature storage unit;
By comparing the inspected feature quantity, which is the feature quantity extracted from the inspected signal, with the model-specific feature quantity range stored in the model-specific feature quantity storage unit, each of the plurality of inspected feature quantities is A model identifying unit that identifies the model that falls within the model-specific feature amount range as the model of the radio wave source that has transmitted the signal to be inspected;
A modulation method determination unit for determining a modulation method of the signal to be inspected,
A radio wave source identifying apparatus that determines the feature quantity used by the model identifying unit in accordance with a modulation scheme determined by the modulation scheme determining unit.
前記被検査信号の変調方式を判定する変調方式判定部をさらに備え、
前記変調方式判定部が判定した変調方式に応じて、前記機種識別部および前記個体識別部のいずれか少なくとも1つが使用する前記特徴量を決める、請求項から請求項のいずれか1項に記載の電波源の識別装置。
A modulation method determination unit for determining a modulation method of the signal to be inspected;
Depending on the modulation scheme the modulation scheme determining section determines the one of the model identification section and the individual identification unit at least one of determining the feature quantity to be used, in any one of claims 6 claims 3 The radio wave source identification device described.
アンテナが受信した電波から生成した離散化された受信信号から検出された被検査信号の複数の特徴量を少なくとも一つの通信中の特徴量を含んで抽出する特徴量抽出部と、
ある機種の電波源が送信する電波が受信されて得られる前記受信信号の複数の前記特徴量のそれぞれが取りうる範囲である機種別特徴量範囲を前記電波源の前記機種ごとに保存する機種別特徴量保存部と、
前記被検査信号から抽出した前記特徴量である被検査特徴量と前記機種別特徴量保存部に保存された前記機種別特徴量範囲とを比較して、複数の前記被検査特徴量のそれぞれが前記機種別特徴量範囲に入る前記機種を、前記被検査信号を送信した前記電波源の前記機種として識別する機種識別部と、
前記受信信号から前記被検査信号を検出する被検査信号検出部とを備えた電波源の識別装置。
A feature quantity extraction unit that extracts a plurality of feature quantities of a signal under test detected from a discretized received signal generated from a radio wave received by an antenna, including at least one feature quantity during communication;
A model-specific feature range that is a range that can be taken by each of the plurality of feature quantities of the received signal obtained by receiving radio waves transmitted from a radio wave source of a certain model. A feature storage unit;
By comparing the inspected feature quantity, which is the feature quantity extracted from the inspected signal, with the model-specific feature quantity range stored in the model-specific feature quantity storage unit, each of the plurality of inspected feature quantities is A model identifying unit that identifies the model that falls within the model-specific feature amount range as the model of the radio wave source that has transmitted the signal to be inspected;
Discrimination device electric wave source and a inspection signal detector for detecting the inspection signal from the received signal.
電波を受信する前記アンテナと、
前記アンテナが前記電波を受信して生成する信号を離散化した前記受信信号を生成し、前記受信信号から前記被検査信号を検出する被検査信号検出部と、
請求項1から請求項13のいずれか1項に記載の電波源の識別装置とを備えた電波源の識別システム。
It said antenna for receiving radio waves,
An inspected signal detector for generating the received signal obtained by discretizing a signal generated by the antenna receiving the radio wave, and detecting the inspected signal from the received signal;
A radio wave source identification system comprising: the radio wave source identification device according to any one of claims 1 to 13 .
JP2016139890A 2016-07-15 2016-07-15 Radio source identification device and identification system Active JP6597505B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016139890A JP6597505B2 (en) 2016-07-15 2016-07-15 Radio source identification device and identification system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016139890A JP6597505B2 (en) 2016-07-15 2016-07-15 Radio source identification device and identification system

Publications (2)

Publication Number Publication Date
JP2018011241A JP2018011241A (en) 2018-01-18
JP6597505B2 true JP6597505B2 (en) 2019-10-30

Family

ID=60995779

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016139890A Active JP6597505B2 (en) 2016-07-15 2016-07-15 Radio source identification device and identification system

Country Status (1)

Country Link
JP (1) JP6597505B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7447669B2 (en) 2020-05-13 2024-03-12 日本電気株式会社 Transmitting device verification device, transmitting device verification method, and program

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0545390A (en) * 1991-08-09 1993-02-23 Mitsubishi Heavy Ind Ltd Classifying and discriminating method of radio wave source
JP2005033312A (en) * 2003-07-08 2005-02-03 Interenergy Co Ltd Radio wave identification apparatus, identification switching terminal, switch apparatus, and identification method
JP2007248110A (en) * 2006-03-14 2007-09-27 Nec Corp Device and method for monitoring electric wave
JP4367479B2 (en) * 2006-11-20 2009-11-18 日本電気株式会社 Digital radio station identification method and digital radio station identification apparatus
JP6143545B2 (en) * 2013-05-15 2017-06-07 三菱電機株式会社 Radio frequency signal identification apparatus and radio frequency signal identification method
JP6503794B2 (en) * 2015-03-06 2019-04-24 日本電気株式会社 Identification device, identification method and program for identification device

Also Published As

Publication number Publication date
JP2018011241A (en) 2018-01-18

Similar Documents

Publication Publication Date Title
AU2007345621B2 (en) Method and apparatus for target discrimination within return signals
US9807549B2 (en) Systems and methods for adaptive multi-feature semantic location sensing
JP2006521756A (en) System and method for waveform classification and characterization using multidimensional higher order statistics
WO2016035439A1 (en) Interference identification device, radio communication device, and interference identification method
EP3791315A1 (en) Radio frequency (rf) object detection using radar and machine learning
US9954714B2 (en) System and method for blindly acquiring frequency hopped spread spectrum signals
JP6597505B2 (en) Radio source identification device and identification system
Sardar et al. Vehicle detection and classification using LTE‐CommSense
WO2021002049A1 (en) Sensor
US11817977B2 (en) Method of determining single sampling frequency of classification target signal in order to predict modulation type of classification target signal, and method and apparatus for predicting modulation type by using classification target signal sampled with single sampling frequency
CN111401185B (en) Electromagnetic signal and equipment association analysis method and system based on radio frequency fingerprint SEV
KR101300650B1 (en) A method of recognizing pri modulation type based on support vector machine, and an apparatus of recognizing pri modulation type
JP6503794B2 (en) Identification device, identification method and program for identification device
KR101317887B1 (en) Radar modulation identification method using gini&#39;s coefficient
Alammar et al. Automatic extraction of signal areas from radio spectrograms based on the Hough transform
KR20170054168A (en) Method and apparatus for processing signal based CFAR in radar system
Kumar et al. A machine learning approach to modulation detection in wireless communications
JP5979942B2 (en) Intra-pulse modulation analysis apparatus and intra-pulse modulation analysis method
Uppal et al. Rich feature deep learning classifier for multiple simultaneous radio signals
CN109495201B (en) Station identification method, device and equipment
JP5757295B2 (en) Database automatic creation device, radio station identification device, radio wave monitoring device, database automatic creation method, and database automatic creation program
US20230246889A1 (en) Modulation scheme identification apparatus, control circuit, storage medium, and modulation scheme identification method
CN116719061A (en) GNSS induced deception detection method based on RMS sliding envelope and SVM
CN118033652A (en) Echo data processing method, device, medium and equipment
JP2005337854A (en) Target categorizing method and device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180724

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190513

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190528

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190708

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190916

R151 Written notification of patent or utility model registration

Ref document number: 6597505

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250