CN111563227B - Fingerprint characteristic parameter extraction method for radiation source signal and radiation source identification - Google Patents

Fingerprint characteristic parameter extraction method for radiation source signal and radiation source identification Download PDF

Info

Publication number
CN111563227B
CN111563227B CN202010365329.3A CN202010365329A CN111563227B CN 111563227 B CN111563227 B CN 111563227B CN 202010365329 A CN202010365329 A CN 202010365329A CN 111563227 B CN111563227 B CN 111563227B
Authority
CN
China
Prior art keywords
signal
radiation source
instantaneous phase
amplitude
dimension
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
CN202010365329.3A
Other languages
Chinese (zh)
Other versions
CN111563227A (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.)
Evaluation Argument Research Center Academy Of Military Sciences Pla China
Original Assignee
Evaluation Argument Research Center Academy Of Military Sciences Pla China
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 Evaluation Argument Research Center Academy Of Military Sciences Pla China filed Critical Evaluation Argument Research Center Academy Of Military Sciences Pla China
Priority to CN202010365329.3A priority Critical patent/CN111563227B/en
Publication of CN111563227A publication Critical patent/CN111563227A/en
Application granted granted Critical
Publication of CN111563227B publication Critical patent/CN111563227B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Computing Systems (AREA)
  • Complex Calculations (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

The application discloses a fingerprint characteristic parameter extraction method of a radiation source signal, which comprises the following steps: (1) carrying out ITD decomposition on the radio frequency digital signal of the radiation source to obtain a series of rotation components with gradually reduced instantaneous frequency, namely: any radio frequency digital signal data sn=(s1,...,sn) Decomposition by ITD into:
Figure DDA0002476551740000011
wherein the content of the first and second substances,
Figure DDA0002476551740000012
refers to the i-th order baseline signal,
Figure DDA0002476551740000013
is the ith order rotational component;
Figure DDA0002476551740000014
is the residual component; here, the radio frequency digital signal data snIs the real part I of any complex digital signalnImaginary part QnOr amplitude En(ii) a (2) Calculating by utilizing the different-order rotation components in the step (1) to obtain corresponding energy distribution, box dimension, information dimension, instantaneous phase mean value, instantaneous phase standard deviation and envelope R characteristics; (3) directly comparing the (1) step signal data snReal part InImaginary part QnOr amplitude EnAnd calculating the corresponding box dimension and information dimension.By using the extracted fingerprint characteristic parameters and combining a machine learning algorithm, the judgment and identification of the radiation signal source can be accurately realized.

Description

Fingerprint characteristic parameter extraction method for radiation source signal and radiation source identification
Technical Field
The technical scheme of the application belongs to the technical field of signal analysis, signal processing and signal source distinguishing and identifying, and particularly aims at obtaining a series of fingerprint characteristic parameters of energy distribution, box dimension, information dimension, instantaneous phase mean value, instantaneous phase standard deviation and envelope R characteristic of different electromagnetic radiation signals after digital processing, and establishing a learning machine by utilizing a machine learning method to determine the type of a signal source according to the signal fingerprint characteristic parameters through autonomous comparison and learning, namely a fingerprint characteristic parameter extraction method and radiation source identification of a radiation source signal.
Background
In modern communication and data transmission technology, the realization of information transmission by electromagnetic field transmission has become a standard information transmission mode. In the transmission of electromagnetic field information, it is necessary to use a device or apparatus that transmits a different electromagnetic field radiation signal, i.e., a signal radiation source. The electromagnetic radiation source emits electromagnetic signals which necessarily carry fine characteristic information of hardware characteristics of equipment due to differences of circuits and devices of the electromagnetic radiation source, and the individual characteristic information of the radiation source is the characteristic which is different from other individuals and exists intrinsically. Therefore, the purpose of identifying the individual identity of the radiation source can be realized in principle by analyzing the fine characteristics of the radiation source, and the core of the method is that the characteristics of the received electromagnetic signals are measured, and the individual process of the radiation source generating the signals is determined according to the existing prior information. In order to realize the process, the characteristics of the fine characteristics of the radiation signals of the radiation source are firstly analyzed: (1) independence: the fine characteristic of the signal of the radiation source is irrelevant to the form of the emission signal, and when the waveform of the emission signal changes, the characteristic is still unchanged; (2) stability: the fine characteristics are stable and do not change significantly due to environmental changes such as temperature, vibration and the like; (3) the testability: means that the characteristic is measurable, and the measurement precision can meet the requirement of individual classification. The subtle features of the radiation source signal are more represented by irregular non-stationary, non-linear, non-gaussian signals, and the nature of the irregular non-stationary, non-linear, non-gaussian signals is difficult to be revealed more deeply by the traditional first-order and second-order analysis methods. The high-order statistical theory is used as an analysis tool of non-stationary signals and has great advantages for analyzing the fine characteristics of the radiation source. The high-order statistical features can be mainly divided into integral bispectrum, selective bispectrum, signal kurtosis and the like. The fine feature set of commonly used radiation sources is characterized by: a carrier frequency offset of the signal; modulation parameters of the signal (frequency modulation index of the FM signal, code rate of the FSK signal); the high order J characteristics of the signal envelope; hilbert spectral symmetry parameters based on EMD; box dimension, information dimension, LZC complexity of the signal envelope; kurtosis of the signal envelope; bispectral features of the signal; the n-th order high order spectral features of the signal. When the radiation source is distinguished and identified, firstly, the characteristic extraction is needed to be carried out on the radiation source fine signals, namely, the most effective characteristic subset is extracted from the original characteristic signal set according to a certain calculation rule, which is equivalent to the nonlinear transformation from the original measurement space to the characteristic space, the transformation processing from the high-dimensional measurement space to the low-dimensional characteristic space is realized, and the information with classification significance is obtained. After the extracted feature information is obtained, the source resolution is realized through classification decision by combining with the artificial intelligence of machine learning, which is the key of the final purpose. The classification decision is the last link of the radiation source identification task and is also a crucial link. The method has the function of establishing a decision rule according to the feature vectors obtained in the feature extraction link, and realizing the discrimination and classification of the tested objects. The main task of classifier design is to complete the nonlinear relation approximation of input-target pairs quickly according to limited training samples, the approximation is a training and learning process of a high-dimensional space, a certain generalization performance needs to be ensured in the approximation process, and the phenomena of over-learning and over-fitting cannot occur.
The existing radiation signal identification technology has the defect of inaccuracy and completeness in the selection of signal fine characteristic information, so that the identification precision and accuracy are low; aiming at the problems, according to the technical scheme, firstly, parameters of characteristic fine information are obtained through mathematical principle analysis, a calculation method of the parameters is provided, and on the basis, machine learning algorithm is combined to obtain more accurate judgment and identification of the radiation signal source.
Disclosure of Invention
The invention aims to further improve the identification accuracy of the radiation source, extract fine characteristic information (radiation source fingerprint information) contained in a radiation signal of the radiation source and realize the identification of a signal source by combining the information with a machine learning algorithm. The scheme for realizing the aim of the invention comprises a fingerprint information extraction method and a signal source identification process, and the specific technical scheme is as follows:
the fingerprint characteristic parameter extraction method of the radiation source signal comprises the following steps of energy distribution, box dimension, information dimension, instantaneous phase mean value, instantaneous phase standard deviation and envelope R characteristic:
(1) carrying out ITD decomposition on the radio frequency digital signal of the radiation source to obtain a series of rotation components with gradually reduced instantaneous frequency, namely: any radio frequency digital signal data sn=(s1,...,sn) Decomposition by ITD into:
Figure GDA0003202620540000021
wherein,
Figure GDA0003202620540000022
refers to the i-th order baseline signal,
Figure GDA0003202620540000023
is the ith order rotational component;
Figure GDA0003202620540000031
is the residual component; here, the radio frequency digital signal data snIs any complex signal;
(2) calculating to obtain energy distribution, box dimension, information dimension, instantaneous phase mean value, instantaneous phase standard deviation and envelope R characteristics corresponding to different rotation components by using the rotation components of different orders obtained in the step (1);
(3) radio frequency digital signal data s directly aiming at step (1)nReal part of (I)nImaginary part QnOr amplitude EnAnd calculating the corresponding box dimension and information dimension.
According to the technical scheme, energy distribution, box dimension, information dimension, instantaneous phase mean value, instantaneous phase standard deviation and envelope R characteristics of different signal data are selected as characteristic parameters in the fingerprint characteristics of the signal transmitted by the signal source, the parameters not only aim at real part and imaginary part data of an original time domain signal, but also can obtain a series of rotating component parts with gradually reduced instantaneous frequency through ITD decomposition of the data, the completeness of fingerprint characteristic parameters is guaranteed, and the parameter data quantity of subsequent judgment is improved.
Radio frequency digital signal s of the above-mentioned radiation sourcen=[s(1),s(2),...s(n)]Expressed in complex numbers: sn=In+jQnAmplitude of which
Figure GDA0003202620540000032
The corresponding energy distribution, box dimension, information dimension, instantaneous phase mean, instantaneous phase standard deviation and envelope R characteristics can also be obtained for the real part, imaginary part and amplitude of these signals.
The energy distribution corresponding to the rotation component in the step (2) is characterized by energy entropy:
Figure GDA0003202620540000033
herein, the
Figure GDA0003202620540000034
Representing the sum of energies, entropy of energies I, of the rotational componentsEERepresenting the energy distribution of the respective order rotation components of the signal by ITD decomposition, IEESmaller, represents a more concentrated energy distribution.
The method for calculating the instantaneous phase mean value and the instantaneous phase standard deviation corresponding to the rotation component in the step (2) comprises the following steps: (S1) extracting the rotation component signal segment between the zero-crossing points, the instantaneous phase is calculated as:
Figure GDA0003202620540000041
where A is1>0,A2> 0 represents the maximum amplitude in the positive and negative half-waves of the signal segment, respectively, [ t [1,t2,…,t5]Being a point of time on a signal segment, t1,t5At the beginning and end of the signal segment, t2,t4At the time point of maximum positive and negative signal amplitude, t3For time points at which the signal amplitude is zero, xtSignal amplitude at time t; (S2) calculating the mean and variance of the instantaneous phase theta to obtain the instantaneous phase mean and the instantaneous phase standard deviation. The calculation method can accurately obtain the phase information of different areas through the phase differential calculation of a plurality of signal segments, and can obtain the instantaneous mean value and the instantaneous standard deviation of the phase through the average and variance calculation, so that the calculation precision is greatly improved.
The box dimension calculation method of different signal data in the step (2) and the step (3) comprises the following steps:
d (Delta) and d (2 Delta) are defined,
Figure GDA0003202620540000042
Figure GDA0003202620540000043
box dimension
Figure GDA0003202620540000044
S herenIs the real part I of the signal or signal data of the rotating componentnImaginary part QnOr amplitude EnFloor (·) denotes the lower rounding operation, and Δ represents the sampling interval.
The information dimension calculation method of different signal data in the step (2) and the step (3) comprises the following steps: time domain continuous signal snConversion to the frequency domain fsThen to the frequency domain signal fsSampling is carried out to obtain a discrete sequence { fs(i) I 1,2,. T }, where T is the length of the pre-processed signal sequence, Δ f (i) fs(i+1)-fs(i) Definition of
Figure GDA0003202620540000045
To calculate the information dimension
Figure GDA0003202620540000046
S herenIs the real part I of the rotating component signal or signal data of the time domainnImaginary part QnOr amplitude En
The envelope R characteristic calculation method in the step (2) comprises
Figure GDA0003202620540000051
That is, it is the ratio of the variance D of the envelope of the signal s to the square of the envelope mean E, where s is any of the different rotation component signals, if the signal is represented as s (t) ═ a (t) × cos (wt) + δ (t), δ (t) is noise, and a (t) is the amplitude of the useful signal, then:
Figure GDA0003202620540000052
where E is the signal amplitude calculation.
In the method, the traditional method for extracting the characteristics from the real signals is utilized, the characteristics of the complex signals obtained by the orthogonal demodulation receiver are extracted to obtain more abundant information for individual identification of the radiation source, and the fine signals of the digital communication signals can be extracted with high precision due to good demodulation performance of ITD orthogonal demodulation.
By combining the extracted fingerprint characteristic parameters with a machine learning algorithm, the application also provides a method for identifying the radiation source, and the specific technical scheme is as follows: (1) constructing a radiation source signal fingerprint characteristic parameter classifier, and inputting known radiation source signal fingerprint characteristic parameters for classifier training; (2) and (3) receiving the radiation source signal by using the classifier completed in the step (1) to perform signal fingerprint characteristic parameter comparison and identification, and determining the radiation source. The radiation source signal fingerprint characteristic parameter classifier is constructed as a Kernel Extreme Learning Machine (KELM).
The kernel-extreme learning machine (KELM) projects the network input into a high-dimensional space by introducing a kernel function to replace random mapping in the ELM, so that in some occasions, a better classification effect is achieved compared with the traditional ELM, and the result is more stable. In the identification method, after a classifier is obtained by giving a fully trained sample for training, the radiation source individual to which the new signal received by the receiver belongs can be judged by analyzing the new signal.
Drawings
FIG. 1 is a schematic diagram of a real part of a signal to obtain a rotation component through ITD decomposition and to obtain a fingerprint characteristic parameter through direct calculation;
FIG. 2 is a schematic diagram of obtaining rotation components through ITD decomposition of the imaginary signal part and obtaining fingerprint characteristic parameters through direct calculation;
FIG. 3 is a schematic diagram of obtaining rotation components through ITD decomposition after calculating amplitudes of real and imaginary parts of signals and obtaining fingerprint characteristic parameters through direct calculation;
fig. 4 is a schematic diagram of instantaneous phase mean and instantaneous phase standard deviation calculation.
Detailed Description
The present invention will be further described in detail with reference to the drawings and examples, as shown in fig. 1,2 and 3. Any one of the digital signals rn=[r(1),r(2),...r(n)]Can be represented in complex form: r isn=In+jQnIn which In=[I(1),I(2),...,I(n)]Represented as a co-directional output signal; and Qn=[Q(1),Q(2),...,Q(n)]Representing quadrature output signals (in-and-out with respect to the received signal). I isn、QnAnd its envelope signal EnFeatures can be extracted by respective rotational components. Wherein the envelope signal En=[E(1),E(2),...,E(n)]The calculation formula of (2) is as follows:
Figure GDA0003202620540000061
due to the envelope signal EnR is discardednThus in feature extraction, EnAnd its rotation components of each order no longer extract the phase information of the signal. Before extracting the fingerprint characteristic parameters, firstly, the signal is subjected to inherent time scale decomposition (ITD) and any radio frequency digital signal data sn=(s1,...,sn) Decomposition by ITD into:
Figure GDA0003202620540000062
wherein,
Figure GDA0003202620540000063
refers to the i-th order baseline signal,
Figure GDA0003202620540000064
is the ith order rotational component;
Figure GDA0003202620540000065
is the residual component; here, the radio frequency digital signal data snIs a complex signal; and then calculating the different-order rotation components obtained by ITD decomposition to obtain the energy distribution, box dimension, information dimension, instantaneous phase mean, instantaneous phase standard deviation and envelope R characteristics corresponding to the components. The same calculation method is directly used for the radio frequency digital signal data snReal part of (I)nImaginary part QnOr amplitude EnAnd calculating the corresponding box dimension and information dimension. The specific calculation method of the parameters comprises the following steps: (1) calculation of energy distribution, characterized by energy entropy:
Figure GDA0003202620540000066
herein, the
Figure GDA0003202620540000067
Representing the sum of energies, entropy of energies I, of the rotational componentsEERepresenting the energy distribution of the respective order rotation components of the signal by ITD decomposition, IEEThe smaller the size, the more concentrated the energy distribution; (2) as shown in fig. 4, (S1) extracting the rotation component signal segment between the zero-crossing points, the instantaneous phase mean and the instantaneous phase standard deviation, and the instantaneous phase calculation formula is:
Figure GDA0003202620540000071
where A is1>0,A2> 0 represents the maximum amplitude in the positive and negative half-waves of the signal segment, respectively, [ t [1,t2,…,t5]Being a point of time on a signal segment, t1,t5At the beginning and end of the signal segment, t2,t4At the time point of maximum positive and negative signal amplitude, t3For time points at which the signal amplitude is zero, xtSignal amplitude at time t; (S2) calculating the mean value and the variance of the instantaneous phase theta to obtain the instantaneous phase mean value and the instantaneous phase standard deviation; (3) the box dimension, as a parameter of the fractal dimension, the size and the change of which can reflect the irregularity and the complexity of the signal, is mainly used for describing the geometrical scale condition of a graph or a structure, and the specific calculation method comprises the following steps: d (Delta) and d (2 Delta) are defined,
Figure GDA0003202620540000072
Figure GDA0003202620540000073
box dimension
Figure GDA0003202620540000074
S herenIs the real part I of the signal or signal data of the rotating componentnImaginary part QnOr amplitude EnFloor (·) denotes the lower rounding operation, Δ represents the sampling interval; (4) the information dimension, which is also a parameter of the fractal dimension, reflects the density degree of spatial distribution, and the influence of carrier frequency change can be eliminated by extracting the information dimension in a frequency domain, and the calculation method comprises the following steps: time domain continuous signal snConversion to the frequency domain fsThen to the frequency domain signal fsSampling is carried out to obtain a discrete sequence { fs(i) I 1,2,. T }, where T is the length of the pre-processed signal sequence, Δ f (i) fs(i+1)-fs(i) Definition of
Figure GDA0003202620540000075
To calculate the information dimension
Figure GDA0003202620540000076
S herenIs a rotating component signal of the time domain, containing the real part InImaginary part QnOr amplitude En(ii) a (5) In general, parameters such as a rising edge, a falling edge, a top fall, and a pulse width are extracted in extracting a waveform feature of a signal, but these parameters have a large variation in an unsteady signal and cannot be used as a signalAnd fingerprint feature extraction, wherein the signal waveform feature can be extracted by extracting high-order statistic features of the envelope. The method has the advantages that the stray output of different radiation source individuals is different, high-order statistic characteristics are different due to the fact that the stray output is reflected on a signal envelope, common signal envelope high-order statistic refers to a J value and an R value, and the envelope R value which is relatively simple to calculate is used as the extracted fingerprint characteristics. The calculation method comprises the following steps:
Figure GDA0003202620540000081
that is, it is the ratio of the variance D of the envelope of the signal s to the square of the envelope mean E, where s is any of the different rotation component signals, if the signal is represented as s (t) ═ a (t) × cos (wt) + δ (t), δ (t) is noise, and a (t) is the amplitude of the useful signal, then:
Figure GDA0003202620540000082
where E is the signal amplitude. The envelope R value can better reflect the envelope characteristic of the radiation source signal and reduce the influence of additive noise, so that the envelope R value has a certain inhibiting effect on the additive noise. For different radiation source signals, due to different spurious output noises, different extra modulations are added to each radiation source signal even under the same working mode, the extra modulations affect the envelope R value, and the envelope R value can realize the identification and classification of different radiation source signals.
After the process is carried out, the fine fingerprint characteristic parameters of the radiation source contained in the signal are extracted, and the parameters are used for radiation source identification: (1) constructing a radiation source signal fingerprint characteristic parameter classifier, and inputting known radiation source signal fingerprint characteristic parameters for classifier training; (2) and (3) receiving the radiation source signal by using the classifier completed in the step (1) to perform signal fingerprint characteristic parameter comparison and identification, and determining the radiation source. The radiation source signal fingerprint characteristic parameter classifier is constructed as a Kernel Extreme Learning Machine (KELM). The kernel-extreme learning machine (KELM) projects the network input into a high-dimensional space by introducing a kernel function to replace random mapping in the ELM, so that in some occasions, a better classification effect is achieved compared with the traditional ELM, and the result is more stable. In the identification method, after a classifier is obtained by giving a fully trained sample for training, the radiation source individual to which the new signal received by the receiver belongs can be judged by analyzing the new signal.
The technical scheme of the application adopts a layered feature extraction method for individual identification of the radiation source. A series of rotating components with gradually reduced instantaneous frequency are obtained by carrying out ITD decomposition on a received radiation source radio frequency signal, then combination characteristics are respectively extracted from an original signal and each rotating component, and the signal is depicted in a multi-level, multi-dimensional and multi-angle mode. The selected characteristics comprise fractal characteristics, statistical characteristics, energy distribution characteristics and the like so as to retain the fine fingerprint characteristics of the signal, thereby realizing individual identification of the radiation source.
The embodiments of the present invention are merely illustrative and not restrictive, and those skilled in the art can modify the embodiments without inventive contribution as required after reading the present specification, but the present invention is protected by patent law within the scope of the appended claims.

Claims (6)

1. The fingerprint characteristic parameter extraction method of the radiation source signal comprises the following steps of energy distribution, box dimension, information dimension, instantaneous phase mean value, instantaneous phase standard deviation and envelope R characteristic:
(1) carrying out ITD decomposition on the radio frequency digital signal of the radiation source to obtain a series of rotation components with gradually reduced instantaneous frequency, namely: any radio frequency digital signal data sn=[s(1),s(2),...s(n)]Decomposition by ITD into:
Figure FDA0003202620530000011
wherein,
Figure FDA0003202620530000012
refers to the i-th order baseline signal,
Figure FDA0003202620530000013
is the ith order rotational component;
Figure FDA0003202620530000014
is the residual component; here, the radio frequency digital signal data snIs a complex signal;
(2) calculating to obtain energy distribution, box dimension, information dimension, instantaneous phase mean value, instantaneous phase standard deviation and envelope R characteristics corresponding to different rotation components by using the rotation components of different orders obtained in the step (1);
(3) directly aiming at the (1) step radio frequency digital signal data snReal part of (I)nImaginary part QnOr amplitude EnCalculating the box dimension and the information dimension corresponding to the box dimension;
the box dimension calculation method of different signal data in the step (2) and the step (3) comprises the following steps: d (Delta) and d (2 Delta) are defined,
Figure FDA0003202620530000015
Figure FDA0003202620530000016
box dimension
Figure FDA0003202620530000017
floor (·) denotes the lower rounding operation, Δ represents the sampling interval.
2. The method of claim 1, wherein the radiation source has a digital RF signal data sn=[s(1),s(2),...s(n)]Expressed in complex numbers: sn=In+jQnAmplitude of which
Figure FDA0003202620530000018
3. The method for extracting fingerprint characteristic parameters of radiation source signals according to claim 1, wherein the energy distribution corresponding to the rotation component in the step (2) is characterized by energy entropy:
Figure FDA0003202620530000019
herein, the
Figure FDA00032026205300000110
Representing the sum of energies, entropy of energies I, of the rotational componentsEERepresenting the energy distribution of the respective order rotation components of the signal by ITD decomposition, IEESmaller, represents a more concentrated energy distribution.
4. The method for extracting fingerprint characteristic parameters of radiation source signals according to claim 1, wherein the method for calculating the instantaneous phase mean and the instantaneous phase standard deviation corresponding to the rotation component in the step (2) comprises: (S1) extracting the rotation component signal segment between the zero-crossing points, the instantaneous phase is calculated as:
Figure FDA0003202620530000021
where A is1>0、A2> 0 represents the maximum amplitude in the positive and negative half-waves of the signal segment, respectively, [ t [1,t2,…,t5]Being a point of time on a signal segment, t1、t5At the beginning and end of the signal segment, t2、t4At the time point of maximum positive and negative signal amplitude, t3For time points at which the signal amplitude is zero, xtSignal amplitude at time t; (S2) calculating the mean and variance of the instantaneous phase theta to obtain the instantaneous phase mean and the instantaneous phase standard deviation.
5. The method for extracting fingerprint characteristic parameters of radiation source signals according to claim 1, wherein the information dimension calculation method of different signal data in the step (2) and the step (3) is: time domain continuous signal snConversion to the frequency domain fsThen to the frequency domain signal fsSampling is carried out to obtain a discrete sequence { fs(i) I 1,2,. T }, where T is the length of the pre-processed signal sequence, Δ f (i) fs(i+1)-fs(i) Definition of
Figure FDA0003202620530000022
To calculate the information dimension
Figure FDA0003202620530000023
6. The method for extracting fingerprint characteristic parameters of radiation source signal according to claim 1, wherein the envelope R characteristic calculation method in step (2) is
Figure FDA0003202620530000024
I.e. it is the ratio of the variance D of the envelope of the signal s to the square of the envelope mean E, if the signal is represented as s (t) ═ a (t) · cos (wt) + δ (t), δ (t) is noise, a (t) is the amplitude of the useful signal, then:
Figure FDA0003202620530000031
where E is the signal amplitude calculation.
CN202010365329.3A 2020-04-30 2020-04-30 Fingerprint characteristic parameter extraction method for radiation source signal and radiation source identification Active CN111563227B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010365329.3A CN111563227B (en) 2020-04-30 2020-04-30 Fingerprint characteristic parameter extraction method for radiation source signal and radiation source identification

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010365329.3A CN111563227B (en) 2020-04-30 2020-04-30 Fingerprint characteristic parameter extraction method for radiation source signal and radiation source identification

Publications (2)

Publication Number Publication Date
CN111563227A CN111563227A (en) 2020-08-21
CN111563227B true CN111563227B (en) 2021-10-15

Family

ID=72070778

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010365329.3A Active CN111563227B (en) 2020-04-30 2020-04-30 Fingerprint characteristic parameter extraction method for radiation source signal and radiation source identification

Country Status (1)

Country Link
CN (1) CN111563227B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112733613A (en) * 2020-12-18 2021-04-30 国网浙江省电力有限公司武义县供电公司 Radiation source identification method based on Hilbert transform and Helbert coefficient characteristics
CN113011346A (en) * 2021-03-19 2021-06-22 电子科技大学 Radiation source unknown signal identification method based on metric learning
CN112911597B (en) * 2021-03-24 2022-11-01 上海电机学院 Internet of things physical layer multilevel feature extraction method based on radio frequency signal fine portrait
CN113219415A (en) * 2021-04-30 2021-08-06 北京理工大学 Interference source individual identification method based on envelope fingerprint characteristics
CN113205076A (en) * 2021-05-31 2021-08-03 西南交通大学 Method and system for extracting behavior characteristics of wireless frame based on MEMD decomposition
CN113343868A (en) * 2021-06-15 2021-09-03 四川九洲电器集团有限责任公司 Radiation source individual identification method and device, terminal and storage medium
CN114531328B (en) * 2022-04-24 2022-07-01 四川轻化工大学 Modulation format identification method based on signal envelope flatness
CN114866162B (en) * 2022-07-11 2023-09-26 中国人民解放军国防科技大学 Signal data enhancement method and system and communication radiation source identification method and system
CN115169406B (en) * 2022-07-15 2023-04-07 中国人民解放军国防科技大学 Instantaneous phase fingerprint feature enhancement method based on empirical mode decomposition

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090462A (en) * 2017-12-29 2018-05-29 哈尔滨工业大学 A kind of Emitter Fingerprint feature extracting method based on box counting dimension

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090462A (en) * 2017-12-29 2018-05-29 哈尔滨工业大学 A kind of Emitter Fingerprint feature extracting method based on box counting dimension

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
基于信号指纹的通信辐射源个体识别技术研究;徐书华;《中国博士学位论文全文数据库(电子期刊)信息科技辑》;20091215(第12期);I138-66 *
基于固有时间尺度分解模型的通信辐射源特征提取算法;桂云川等;《计算机应用研究》;20170430;第34卷(第4期);1172-1175 *
基于复杂度特征的通信辐射源个体识别;李靖超;《中国博士学位论文全文数据库(电子期刊)信息科技辑》;20171115(第11期);I136-24 *
基于小波包变换的爆破振动信号能量熵特征分析;王伟等;《爆破器材》;20191231;第48卷(第6期);19-23 *
桂云川等.基于固有时间尺度分解模型的通信辐射源特征提取算法.《计算机应用研究》.2017,第34卷(第4期),1172-1175. *

Also Published As

Publication number Publication date
CN111563227A (en) 2020-08-21

Similar Documents

Publication Publication Date Title
CN111563227B (en) Fingerprint characteristic parameter extraction method for radiation source signal and radiation source identification
CN109254274A (en) A kind of Radar emitter discrimination method based on Fusion Features
CN109559758B (en) Method for converting texture image into tactile signal based on deep learning
CN110244271A (en) Radar emitter sorting recognition methods and device based on multiple simultaneous compressed transform
CN110443293A (en) Based on double zero sample image classification methods for differentiating and generating confrontation network text and reconstructing
CN105572499B (en) A kind of eye pattern generation method based on vector network analyzer
CN103941254A (en) Soil physical property classification recognition method and device based on geological radar
CN108919203B (en) Radar active interference identification method and system
Li A new robust signal recognition approach based on holder cloud features under varying SNR environment
CN107561497B (en) FSK and multiple non-linear frequency modulation signal identification and parameter estimation method
CN109409525B (en) Quantum bit quantum state determines that method and quantum bit read signal resolution method
CN108171119B (en) SAR image change detection method based on residual error network
CN104796365A (en) Modulating signal recognition method based on complexity feature under low signal to noise ratio
CN114692665A (en) Radiation source open set individual identification method based on metric learning
CN106940688B (en) A kind of mathematical model modeling method of EMD chaff component identification
CN104655583A (en) Fourier-infrared-spectrum-based rapid coal quality recognition method
CN104143115A (en) Technological method for achieving soil water content classified identification through geological radar technology
CN110988511A (en) Power electronic converter nonlinear identification method based on multiple entropy characteristic extraction
CN110160778A (en) Gearbox fault state identification method based on sequential hypothesis testing
CN111191502B (en) Stick-slip and jump drill abnormal working condition identification method based on drill string vibration signal
CN117331047A (en) Human behavior data analysis method and system based on millimeter wave radar
CN116340533B (en) Satellite-borne electromagnetic spectrum big data intelligent processing system based on knowledge graph
Xu et al. Individual radio transmitter identification based on spurious modulation characteristics of signal envelop
CN107944469A (en) A kind of low signal-to-noise ratio microseismic event discrimination method based on rearrangement ST
CN115563465A (en) Characteristic coding method for individual identification of radiation source

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