CN113657138A - Radiation source individual identification method based on equipotential planet chart - Google Patents

Radiation source individual identification method based on equipotential planet chart Download PDF

Info

Publication number
CN113657138A
CN113657138A CN202010401692.6A CN202010401692A CN113657138A CN 113657138 A CN113657138 A CN 113657138A CN 202010401692 A CN202010401692 A CN 202010401692A CN 113657138 A CN113657138 A CN 113657138A
Authority
CN
China
Prior art keywords
radiation source
equipotential
signal
signal data
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010401692.6A
Other languages
Chinese (zh)
Other versions
CN113657138B (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN202010401692.6A priority Critical patent/CN113657138B/en
Publication of CN113657138A publication Critical patent/CN113657138A/en
Application granted granted Critical
Publication of CN113657138B publication Critical patent/CN113657138B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R21/00Arrangements for measuring electric power or power factor
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08Feature extraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Measurement Of Radiation (AREA)

Abstract

A radiation source individual identification method based on an equipotential planet chart comprises the following steps: s1: acquiring signal data of a plurality of radiation source individuals with the same type; s2: carrying out data preprocessing on the signal data to obtain a signal constellation diagram; s3: performing point density calculation on data points in the signal constellation diagram, and coloring the signal constellation diagram according to a point density calculation result to obtain an equipotential sphere diagram; s4: performing feature extraction and classification on the equipotential planet map by using a convolutional neural network; s5: and outputting the recognition result. The method converts IQ two-path signal data from a signal domain to an image domain based on an equipotential planet image, automatically and effectively extracts deep features of radiation source individuals carried by the signals through a convolutional neural network, corresponds each equipotential planet image to the signal data, and corresponds the signal data to the radiation source individuals, thereby more effectively identifying and classifying the radiation source individuals.

Description

Radiation source individual identification method based on equipotential planet chart
Technical Field
The invention relates to a radiation source individual identification method based on an equipotential planet chart, and belongs to the technical field of radiation source signal identification.
Background
Specific Emitter Identification (SEI) refers to a technique of performing characteristic measurement on a received electromagnetic signal and determining an individual radiation source generating the signal according to existing a priori information. The method has extremely important military significance for correctly identifying the detected radiation source in a complex and changeable electromagnetic environment. Meanwhile, with the advent of the age of 5G and the internet of things, the security of wireless networks is also facing a huge challenge. The individual identification technology of the radiation source can be applied to projects such as wireless access authentication, electromagnetic environment supervision and the like, and the safety of a wireless network can be effectively improved. In the civil field, the radiation source individual identification technology is one of key technologies in the fields of cognitive radio, radio positioning, communication equipment fault detection and identification and the like, and has very important function.
With the development of science and technology, the electromagnetic signal environment of the signal is increasingly complex, the density of a radiation source is multiplied, and the pattern of the electromagnetic signal is complex and changeable, which all bring unprecedented challenges to the radiation source identification. The general radiation source identification method is difficult to realize accurate, stable and reliable radiation source identification in a complex electromagnetic environment. Therefore, research and development of new radiation source identification methods and identification devices are urgently needed.
Disclosure of Invention
The invention aims to solve the technical problem that the prior art is not enough, and provides a radiation source individual identification method based on an equipotential planet map.
The technical problem to be solved by the invention is realized by the following technical scheme:
the invention provides a radiation source individual identification method based on an equipotential planet chart, which comprises the following steps of:
s1: acquiring signal data of a plurality of radiation source individuals with the same type;
s2: carrying out data preprocessing on the signal data to obtain a signal constellation diagram;
s3: performing point density calculation on data points in the signal constellation diagram, and coloring the signal constellation diagram according to a point density calculation result to obtain an equipotential sphere diagram;
s4: performing feature extraction and classification on the equipotential planet map by using a convolutional neural network;
s5: and outputting the recognition result.
Preferably, the signal data includes I-path signal data and Q-path signal data.
Preferably, the data preprocessing is performed on the signal data to obtain a signal constellation diagram, and specifically:
and taking the I-path signal data in the signal data of each sampling point of the plurality of radiation source individuals as the abscissa of a signal constellation diagram, and taking the Q-path signal data in the signal data of each sampling point as the ordinate of the signal constellation diagram, thereby mapping the signal data of the plurality of radiation source individuals into the signal constellation diagram.
Preferably, a window function is used for the point density calculation.
Preferably, the window function is a density window function, and the normalized point density ρ at the sampling point is:
Figure BDA0002487787910000021
where ρ (i) represents the normalized point density of the ith sample point, h (i) represents the abscissa of the signal constellation diagram at which the ith sample point is obtained, v (i) represents the ordinate of the signal constellation diagram at which the ith sample point is obtained, r represents half of the side length of the square region in which the normalized point density is calculated, N represents the total number of sample points in the signal constellation diagram, f (x) satisfies the following expression:
Figure BDA0002487787910000022
preferably, the r ═ sqrt ((range (x)/30) ^2+ (range (y)/30) ^2),
wherein x represents the I path signal data, y represents the Q path signal data, range function represents the difference between the maximum value and the minimum value in the solving vector, and sqrt represents the square root function.
Preferably, the coloring the signal constellation according to the result of the point density calculation specifically includes:
and selecting different colors, and coloring the data points with different normalized point density values or different normalized point density ranges in the signal constellation diagram.
Preferably, the normalized dot density from high to low is expressed sequentially using the colors in a yellow to green to blue gradient color bar.
Preferably, the convolutional neural network is AlexNet.
In summary, the IQ two-path signal data are converted from the signal domain to the map domain based on the equipotential star map, the deep features of the radiation source individuals carried by the signals are automatically and effectively extracted through the convolutional neural network, each equipotential star map corresponds to the signal data, and the signal data corresponds to the radiation source individuals, so that the radiation source individuals are more effectively identified and classified.
The technical solution of the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
Drawings
FIG. 1 is a block diagram of the individual identification method of a radiation source based on an equipotential sphere diagram according to the present invention;
fig. 2 is a comparison diagram of the radiation source individual identification method based on the equipotential sphere diagram and the traditional feature extraction method.
Detailed Description
FIG. 1 is a block diagram of the individual identification method of a radiation source based on an equipotential sphere diagram. As shown in fig. 1, the present invention provides an equipotential sphere diagram-based individual radiation source identification method, which includes the following steps:
s1: acquiring signal data of a plurality of radiation source individuals with the same type;
s2: carrying out data preprocessing on the signal data to obtain a signal constellation diagram;
s3: performing point density calculation on data points in the signal constellation diagram, and coloring the signal constellation diagram according to a point density calculation result to obtain an equipotential sphere diagram;
s4: performing feature extraction and classification on the equipotential planet map by using a convolutional neural network;
s5: and outputting the recognition result.
In S1, the signal data includes I-path signal data and Q-path signal data.
In S2, the data preprocessing is performed on the signal data to obtain a signal constellation, specifically, I-path signal data in the signal data of each sampling point of the multiple radiation source individuals is used as an abscissa of the signal constellation, and Q-path signal data in the signal data of each sampling point is used as an ordinate of the signal constellation, so that the signal data of the multiple radiation source individuals is mapped to the signal constellation. It should be added that the above examples are only illustrative, the present invention does not limit the specific way of obtaining the signal constellation through data preprocessing, and those skilled in the art can select a suitable existing method to process the signal data into the signal constellation according to practical situations.
In the signal constellation, different regions have different sample point densities. In the present invention, the point density calculation described in S3 is exemplarily performed using a window function.
Specifically, in the present invention, the window function is preferably a density window function, and when the density window function slides on the signal constellation diagram, the density window function counts the number of points in different region windows, and divides the number of the sampling points of the whole signal constellation diagram by the number of the points to obtain a normalized point density ρ at the sampling points, that is, the normalized point density ρ is obtained
Figure BDA0002487787910000041
Where ρ (i) represents the normalized point density of the ith sample point, h (i) represents the abscissa of the signal constellation diagram at which the ith sample point is obtained, v (i) represents the ordinate of the signal constellation diagram at which the ith sample point is obtained, r represents one half of the side length of the square region in which the normalized point density is calculated, N represents the total number of sample points in the signal constellation diagram, and f (x) satisfies the following expression:
Figure BDA0002487787910000042
preferably, r ═ sqrt ((range (x)/30) ^2+ (range (y)/30) ^2), where x represents the I-path signal data, y represents the Q-path signal data, the range function represents the difference between the maximum value and the minimum value in the solution vector, and sqrt represents the square root function.
By the above method, the point density at each data point in the signal constellation diagram can be obtained, and of course, the invention is not limited to the calculation method of the point density, and those skilled in the art can select other well-known point density calculation methods.
In S3, the coloring the signal constellation according to the result of the point density calculation specifically includes: and selecting different colors, and coloring the data points with different normalized point density values or different normalized point density ranges in the signal constellation diagram. The present invention does not limit the specific rule of coloring, for example, coloring the signal constellation diagram by a preset color bar containing the corresponding relationship between the dot density and the color. Preferably, the normalized dot density from high to low is expressed sequentially using the colors in a yellow to green to blue gradient color bar.
Through the coloring process in S3, an image having RGB information can be formed, that is, an equi-star map can be obtained. The traditional signal constellation diagram has defects when used for extracting signal characteristics, for example, in the actual data acquisition process, the internal noise of a device can seriously pollute radio frequency signals, so that the constellation diagrams of different acquired radio frequency signals have the same graph.
In S4, the present invention does not limit the type of the convolutional neural network as long as it can implement the picture recognition classification. Preferably, the convolutional neural network AlexNet is used for automatically extracting features and classifying the same. Since AlexNet is a public network, the specific process in S4 is not described herein again.
Fig. 2 is a comparison diagram of the radiation source individual identification method based on the equipotential sphere diagram and the traditional feature extraction method. In fig. 2, the method of using a convolutional neural network AlexNet in combination with an equi-potential star map is compared with the conventional Welch method, Yule-Walker method and Burg method for extracting power spectrum features from signal data. The data set used in this comparative experiment is introduced as follows: eight power amplifier output data of the same model and the same batch are collected by a baseband signal receiver to be used as a research object, the operating center frequency of the amplifier module is 433MHz, each individual uses 100 samples to research, and each sample uses 20000 sample points. When the power spectrum feature is extracted, 2048 FFT points are used, PCA is used for dimensionality reduction after a power spectrum is obtained, and finally a K nearest neighbor classifier is used for classification and identification.
As can be seen from FIG. 2, the identification accuracy of the individual identification method of the radiation source of the invention using the convolutional neural network and the equipotential sphere diagram is superior to that of the conventional method.
In summary, the IQ two-path signal data are converted from the signal domain to the map domain based on the equipotential star map, the deep features of the radiation source individuals carried by the signals are automatically and effectively extracted through the convolutional neural network, each equipotential star map corresponds to the signal data, and the signal data corresponds to the radiation source individuals, so that the radiation source individuals are more effectively identified and classified.

Claims (9)

1. A radiation source individual identification method based on an equipotential planet chart is characterized by comprising the following steps:
s1: acquiring signal data of a plurality of radiation source individuals with the same type;
s2: carrying out data preprocessing on the signal data to obtain a signal constellation diagram;
s3: performing point density calculation on data points in the signal constellation diagram, and coloring the signal constellation diagram according to a point density calculation result to obtain an equipotential sphere diagram;
s4: performing feature extraction and classification on the equipotential planet map by using a convolutional neural network;
s5: and outputting the recognition result.
2. The method for identifying individuals as a radiation source based on an equipotential sphere according to claim 1, wherein the signal data includes I-path signal data and Q-path signal data.
3. The method for identifying individuals as a radiation source based on an equi-potential star map as claimed in claim 2, wherein the data preprocessing is performed on the signal data to obtain a signal constellation map specifically as follows:
and taking the I-path signal data in the signal data of each sampling point of the plurality of radiation source individuals as the abscissa of a signal constellation diagram, and taking the Q-path signal data in the signal data of each sampling point as the ordinate of the signal constellation diagram, thereby mapping the signal data of the plurality of radiation source individuals into the signal constellation diagram.
4. The method for identifying individuals as a source of radiation based on an equi-potential star map as claimed in claim 3 wherein said point density calculation is performed using a windowing function.
5. The method for identifying individuals as a radiation source based on an equipotential planet map according to claim 3, wherein the window function is a density window function, and the normalized point density p at the sampling points is:
Figure FDA0002487787880000011
where ρ (i) represents the normalized point density of the ith sample point, h (i) represents the abscissa of the signal constellation diagram at which the ith sample point is obtained, v (i) represents the ordinate of the signal constellation diagram at which the ith sample point is obtained, r represents half of the side length of the square region in which the normalized point density is calculated, N represents the total number of sample points in the signal constellation diagram, f (x) satisfies the following expression:
Figure FDA0002487787880000021
6. the method of claim 5, wherein r ═ sqrt ((range (x)/30) ^2+ (range (y)/30) ^2),
wherein x represents I path signal data, y represents Q path signal data, range function represents the difference between the maximum value and the minimum value in the solving vector, and sqrt represents the root opening number.
7. The method for identifying an individual radiation source based on an equipotential sphere of claim 1, wherein the coloring of the signal constellation according to the result of the point density calculation specifically comprises:
and selecting different colors, and coloring the data points with different normalized point density values or different normalized point density ranges in the signal constellation diagram.
8. The method for identifying individuals as a radiation source based on an equipotential sphere according to claim 7, wherein the normalized dot density from high to low is expressed sequentially by using colors in a color bar with a gradual change from yellow to green to blue.
9. The method for identifying individuals as sources of radiation based on an equinox global map as claimed in claim 1 wherein said convolutional neural network is AlexNet.
CN202010401692.6A 2020-05-12 2020-05-12 Radiation source individual identification method based on equipotential star map Active CN113657138B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010401692.6A CN113657138B (en) 2020-05-12 2020-05-12 Radiation source individual identification method based on equipotential star map

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010401692.6A CN113657138B (en) 2020-05-12 2020-05-12 Radiation source individual identification method based on equipotential star map

Publications (2)

Publication Number Publication Date
CN113657138A true CN113657138A (en) 2021-11-16
CN113657138B CN113657138B (en) 2024-05-21

Family

ID=78476753

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010401692.6A Active CN113657138B (en) 2020-05-12 2020-05-12 Radiation source individual identification method based on equipotential star map

Country Status (1)

Country Link
CN (1) CN113657138B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116628472A (en) * 2023-04-28 2023-08-22 哈尔滨工程大学 Characteristic correlation-based radiation source individual signal diagram structure mapping method

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4516263A (en) * 1982-04-30 1985-05-07 The Charles Stark Draper Laboratory, Inc. Spatially integral, video signal processor
JP2006010418A (en) * 2004-06-24 2006-01-12 Nec Engineering Ltd Fixed star sensor
JP2008040683A (en) * 2006-08-03 2008-02-21 Matsushita Electric Works Ltd Signal identification method and signal identification device
CN101278316A (en) * 2005-08-02 2008-10-01 美国西门子医疗解决公司 System and method for automatic segmentation of vessels in breast MR sequences
CN102231720A (en) * 2011-07-25 2011-11-02 南京信息工程大学 Wavelet blind equalization method for fusing spline function Renyi entropy and time diversity
JP2017150924A (en) * 2016-02-24 2017-08-31 三菱電機株式会社 Angle measurement device
EP3232371A1 (en) * 2016-04-15 2017-10-18 Ricoh Company, Ltd. Object recognition method, object recognition device, and classifier training method
CN108427987A (en) * 2018-03-08 2018-08-21 四川大学 A kind of Modulation Mode Recognition method based on convolutional neural networks
CN110427893A (en) * 2019-08-06 2019-11-08 西安电子科技大学 A kind of specific emitter identification method, apparatus and computer storage medium based on convolutional neural networks

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4516263A (en) * 1982-04-30 1985-05-07 The Charles Stark Draper Laboratory, Inc. Spatially integral, video signal processor
JP2006010418A (en) * 2004-06-24 2006-01-12 Nec Engineering Ltd Fixed star sensor
CN101278316A (en) * 2005-08-02 2008-10-01 美国西门子医疗解决公司 System and method for automatic segmentation of vessels in breast MR sequences
JP2008040683A (en) * 2006-08-03 2008-02-21 Matsushita Electric Works Ltd Signal identification method and signal identification device
CN102231720A (en) * 2011-07-25 2011-11-02 南京信息工程大学 Wavelet blind equalization method for fusing spline function Renyi entropy and time diversity
JP2017150924A (en) * 2016-02-24 2017-08-31 三菱電機株式会社 Angle measurement device
EP3232371A1 (en) * 2016-04-15 2017-10-18 Ricoh Company, Ltd. Object recognition method, object recognition device, and classifier training method
CN108427987A (en) * 2018-03-08 2018-08-21 四川大学 A kind of Modulation Mode Recognition method based on convolutional neural networks
CN110427893A (en) * 2019-08-06 2019-11-08 西安电子科技大学 A kind of specific emitter identification method, apparatus and computer storage medium based on convolutional neural networks

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
J. CHANG, Y. XIAO AND Z. ZHANG: "Wireless Physical-Layer Identification Assisted 5G Network Security", 《IEEE INFOCOM 2019》, pages 1 - 5 *
王威;李诗娴;王新;: "基于星座图的通信辐射源个体识别方法", 湖南城市学院学报(自然科学版), no. 05 *
赵春晖;杜宇;: "基于星座图的混合MPSK信号盲识别算法", 黑龙江大学工程学报, no. 04, 30 November 2012 (2012-11-30) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116628472A (en) * 2023-04-28 2023-08-22 哈尔滨工程大学 Characteristic correlation-based radiation source individual signal diagram structure mapping method
CN116628472B (en) * 2023-04-28 2024-08-30 哈尔滨工程大学 Characteristic correlation-based radiation source individual signal diagram structure mapping method

Also Published As

Publication number Publication date
CN113657138B (en) 2024-05-21

Similar Documents

Publication Publication Date Title
CN109802905B (en) CNN convolutional neural network-based digital signal automatic modulation identification method
CN106845339B (en) Mobile phone individual identification method based on bispectrum and EMD fusion characteristics
CN114564982B (en) Automatic identification method for radar signal modulation type
Sun et al. Automatic modulation classification using techniques from image classification
CN112749633A (en) Separate and reconstructed individual radiation source identification method
CN115294615A (en) Radio frequency fingerprint identification method based on machine learning
CN113657138B (en) Radiation source individual identification method based on equipotential star map
Wang et al. A cooperative spectrum sensing method based on signal decomposition and K-medoids algorithm
Ying et al. Differential complex-valued convolutional neural network-based individual recognition of communication radiation sources
Zhang et al. Data augmentation aided few-shot learning for specific emitter identification
CN109446910B (en) Communication radiation source signal classification and identification method
Geng et al. Spectrum sensing for cognitive radio based on feature extraction and deep learning
Wang et al. A new method of automatic modulation recognition based on dimension reduction
Chang et al. Wireless physical-layer identification assisted 5g network security
CN111310680B (en) Radiation source individual identification method based on deep learning
CN110365434B (en) Multi-antenna cooperative spectrum sensing method based on information geometry and differential evolution clustering algorithm
Zhang et al. A spectrum sensing algorithm for OFDM signal based on deep learning and covariance matrix graph
Wang et al. A radio frequency fingerprinting identification method based on energy entropy and color moments of the bispectrum
CN115809426A (en) Radiation source individual identification method and system
Cai et al. Toward Intelligent Lightweight and Efficient UAV Identification With RF Fingerprinting
Huang et al. Radio frequency fingerprint identification method based on ensemble learning
CN111404852B (en) Modulation mode identification method based on amplitude and spectral amplitude characteristics
CN111191515B (en) High-precision frequency spectrum identification method and system based on deep learning
CN106936744B (en) signal modulation identification method based on dynamic ideal solution
CN112083448B (en) Satellite navigation system-oriented interference signal classification recognition feature extraction method and system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant