CN108108712A - A kind of Emitter Fingerprint feature extracting method based on variance dimension - Google Patents
A kind of Emitter Fingerprint feature extracting method based on variance dimension Download PDFInfo
- Publication number
- CN108108712A CN108108712A CN201711489288.3A CN201711489288A CN108108712A CN 108108712 A CN108108712 A CN 108108712A CN 201711489288 A CN201711489288 A CN 201711489288A CN 108108712 A CN108108712 A CN 108108712A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- subinterval
- emitter
- variance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/211—Selection of the most significant subset of features
- G06F18/2113—Selection of the most significant subset of features by ranking or filtering the set of features, e.g. using a measure of variance or of feature cross-correlation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/16—Classification; Matching by matching signal segments
- G06F2218/18—Classification; Matching by matching signal segments by plotting the signal segments against each other, e.g. analysing scattergrams
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Signal Processing (AREA)
- Complex Calculations (AREA)
- Image Analysis (AREA)
Abstract
A kind of Emitter Fingerprint feature extracting method based on variance dimension, the present invention relates to Emitter Fingerprint feature extracting methods.The problem of being difficult to meet specific emitter identification validity and reliability demand, cause specific emitter identification accuracy low the purpose of the present invention is to solve traditional characteristic parameter.A kind of Emitter Fingerprint feature extracting method detailed process based on variance dimension is:First, segment processing is carried out to the one-dimensional emitter Signals received, obtains one-dimensional emitter Signals section;2nd, the one-dimensional emitter Signals section obtained to one carries out variance Dimension Characteristics extraction, obtains variance Dimension Characteristics vector.The present invention is used for specific emitter identification field.
Description
Technical field
The present invention relates to Emitter Fingerprint feature extracting methods.
Background technology
Specific emitter identification, also known as " Emitter Fingerprint identification " or " Specific Emitter Identification ", it is main using passive
The emitter Signals of reception analyze its fine feature parameter, and the difference of same type radiation source is told using Finger print characteristic abstract
Individual.Due to the difference of radiation source physical component so that each radiation source waveform suffers from specific characteristic, is known as fingerprint characteristic,
And the key technology of specific emitter identification is the selection and extraction of fingerprint characteristic.Specific emitter identification is dual-use
Aspect has many applications.Military aspect carries out individual identification mainly for Radar emitter;Civilian aspect essentially consists in radio frequency
Identification, network node is recognized each other, information is collected evidence and the fields such as communication security.In addition, fractal dimension is as the important interior of fractal theory
Hold, compared and be widely applied to the fields such as seismic wave detection, image procossing and material science.Radar emitter with noise
Signal can regard a kind of time series with certain self-similarity nature as, it be able to effectively be portrayed with point shape.
Existing radiation source Feature Extraction Technology is mainly special using the temporal signatures of emitter Signals, frequency domain character, statistics
The traditional parameters such as sign, with increasingly complicated and hardware technique the raising of radiation source waveform, traditional characteristic parameter is difficult to meet
Specific emitter identification validity and reliability demand causes specific emitter identification accuracy low.
The content of the invention
The purpose of the present invention is to solve traditional characteristic parameters to be difficult to meet specific emitter identification validity and reliable
Property demand, the problem of causing specific emitter identification accuracy low, and propose that a kind of Emitter Fingerprint based on variance dimension is special
Levy extracting method.
A kind of Emitter Fingerprint feature extracting method detailed process based on variance dimension is:
Step 1: carrying out segment processing to the one-dimensional emitter Signals received, one-dimensional emitter Signals section S is obtained1,
S2,…Si…,Sn;
SiFor i-th of one-dimensional emitter Signals section, 1≤i≤n, n values are positive integer;
Process is:
The one-dimensional emitter Signals S that step 1 one, input receive, according to the total of the one-dimensional emitter Signals S received
Long and feature extraction segments demand, setting signal segment length w and segments n, each signal segment length are identical;
Step 1 two determines sliding window step delta, and sliding window step delta is obtained by formula (1):
WhereinRepresent positive bracket function;
Step 1 three, call function G=enframe (S, w, Δ) obtain return value G, the G form of function enframe such as
Under:
G is that a line number is n, and columns is the matrix of w;
Obtain one-dimensional emitter Signals section S1,S2,…Si…,Sn;
Step 2: carrying out variance Dimension Characteristics extraction to the one-dimensional emitter Signals section that step 1 obtains, variance dimension is obtained
Number feature vector.
Beneficial effects of the present invention are:
The present invention starts with from Fractal Domain, by the use of variance dimension as the fingerprint characteristic of emitter Signals, utilizes signal subsection
Processing increases intrinsic dimensionality, completes the extraction of Emitter Fingerprint feature, meets specific emitter identification validity and reliability need
It asks, improves specific emitter identification accuracy.Emulation experiment the results show that in signal-to-noise ratio under conditions of 13dB, to utilize
The variance Dimension Characteristics of extraction are classified, and the accuracy of specific emitter identification is more than 95%.
Specific emitter identification refers to distinguish very high same type of of similarity using fingerprint (individual exclusive) feature
Different radiation source individuals.Specific emitter identification is an important development direction in Radar recognition field, although type identification
Technology is gradually ripe, but the individual information included in type identification is less, and has radiated source information needed for engineer application
Type information is not limited to, often also needs to radiation source individual information.Individual identification precisely in order to obtain radiation source individual information, can
To say that this is the Studies on Emitters ID more higher level than type identification.With military application examples, Radar emitter individual identification can
To distinguish the Different Individual of same kind radar, this is vital for commanding and decision-making, while is also that type identification is difficult
With what is reached.And the existing research to specific emitter identification often utilizes the tradition such as temporal signatures, frequency domain character, statistical nature
Parameter, the present invention improve specific emitter identification validity and reliability using variance Dimension Characteristics
The raising of specific emitter identification validity and reliability, in the equal important in inhibiting of dual-use aspect.Army
Individual identification is carried out mainly for Radar emitter with aspect, helps to improve electronic reconnaissance ability, it is excellent to improve our electronic warfare
Gesture captures electronic warfare battlefield commanding elevation;Civilian aspect contributes to radio frequency identification, network node to recognize each other, information is collected evidence and communication security
Technology being improved and improving.
Description of the drawings
Fig. 1 realizes the general frame for one method of the specific embodiment of the invention;
Fig. 2 is one signal subsection flow chart of the specific embodiment of the invention;
Fig. 3 extracts flow chart for variance Dimension Characteristics of the present invention;
Fig. 4 is the simulation experiment result figure the present invention is based on variance dimension.
Specific embodiment
Specific embodiment one:Illustrate present embodiment with reference to Fig. 1, Fig. 2, Fig. 3, one kind of present embodiment is based on variance
The Emitter Fingerprint feature extracting method detailed process of dimension is:
Step 1: carrying out segment processing to the one-dimensional emitter Signals received, one-dimensional emitter Signals section S is obtained1,
S2,…Si…,Sn, to show the fine feature of signal and increase intrinsic dimensionality;SiFor i-th of one-dimensional emitter Signals section, 1≤i
≤ n, n value are positive integer;
Process is:
In order to describe emitter Signals from subtleer angle, and also between better discriminating between different emitter Signals
Difference, emitter Signals are first divided into several signal segments.The concrete methods of realizing of signal subsection is to utilize to slide window function, again
Framing function (enframe) is, signal adding window and the framing being usually used in Speech processing;
The one-dimensional emitter Signals S that step 1 one, input receive, according to the total of the one-dimensional emitter Signals S received
Long and feature extraction demand (segments) setting signal segment length (slip window width) w and segments n, each signal segment length phase
Together;
Step 1 two determines sliding window step delta, and sliding window step delta is obtained by formula (1):
WhereinRepresent positive bracket function (real number x is mapped to the smallest positive integral more than or equal to x);
Step 1 three, call function G=enframe (S, w, Δ) obtain return value G, the G form of function enframe such as
Under:
G is that a line number is n, and columns is the matrix of w, therefore the calculating of fractal dimension must carry out line by line.
Obtain one-dimensional emitter Signals section S1,S2,…Si…,Sn;
Step 2: carrying out variance Dimension Characteristics extraction to the one-dimensional emitter Signals section that step 1 obtains, variance dimension is obtained
Number feature vector;
Using variance Dimension Characteristics vector as the input of grader, categorized device is exporting the identification of one-dimensional emitter Signals just
True rate.
Radiation source individual can be distinguished by grader, increases feature samples, you can obtain discrimination.
The grader is support vector machines (SVM, support vector machine) or backpropagation (BP, back
Propagation) neural network classifier.
Specific embodiment two:The present embodiment is different from the first embodiment in that:To step in the step 2
One obtained one-dimensional emitter Signals section carries out variance Dimension Characteristics extraction, obtains variance Dimension Characteristics vector;Detailed process is:
Variance dimension is a kind of fractal dimension based on Hurst indexes, can be used to the Fractal Properties of analysis time sequence.Side
Poor dimension DσIt is shown below with the relation of Hurst indexes H:
Dσ=E+1-H
Wherein E represents Euclid dimensions, for One-dimension Time Series, takes E=1;
The classical calculation of Hurst indexes is rescaled range method.R/S analytic approach is initially by Britain hydrologist
The method that H.E.Hurst is proposed when studying the storage capability of Nile reservoir.R/S analytic approach calculates simply, only needs average, pole
Three difference, mean square deviation amounts, have unique advantage, in stock in the fractal characteristic of analysis time sequence and long-term memory process
Have preferable application in the analysis in market, and the present invention uses it for the extraction of radiation source " fingerprint " feature expansionaryly.
The present invention using R/S analysis methods calculate variance dimension basic step be:
Step 2 one, as i=1, by SiThe isometric subinterval that β initial length is A is divided into, with { dtRepresent sub-district
Interior sequence;β values are positive integer;
Step 2 two, calculating add up the very poor of deviation in each subinterval;
For each subinterval, if:
Wherein, EqFor d in q-th of subintervalpAverage value, q=1,2 ... β, dpFor { dtIn p-th of element, Zt,qFor
The accumulative deviation in q subinterval;
The very poor of all accumulative deviation values is obtained by formula (4):
Rq=max (Zt,q)-min(Zt,q) (4)
Step 2 three calculates each subinterval standard deviation:
If SqRepresent sequence { d in q-th of subintervaltStandard deviation, then:
Sequence { dtRepresent sequence in each subinterval;
Step 2 four calculates each subinterval rescaled range, and averages;
The length for making subinterval is A changing with 2 for the exponential law at bottom as shown in formula (6),
Each subinterval rescaled range is calculated, and is averaged:
Step 2 one is repeated to step 2 four, untilPerform step 2 five;
Step 2 five, take the logarithm to the value of corresponding rescaled range R/S and length A after linear fit (it is corresponding be
Cycle-index is identical), obtain Hurst indexes H(i);
Step 2 six utilizes Hurst indexes H(i)Calculate variance Dimension Characteristics vector
It is provided by formula (8):
Wherein E represents Euclid dimensions;
Step 2 seven judges whether i is equal to n, if so,Terminate;If not,
Step 2 one is performed, until i is equal to n.
Other steps and parameter are same as the specific embodiment one.
Specific embodiment three:The present embodiment is different from the first and the second embodiment in that:In the step 2 one
As i=1, by SiIt is divided into the isometric subinterval that β initial length is A;β values are positive integer;Process is:
The one-dimensional emitter Signals section S that signal segment length for interception is wi, by one-dimensional emitter Signals section SiIt is divided into β
A initial length is the isometric subinterval of A, with sequence { dtRepresent sequence in each subinterval.
Other steps and parameter are the same as one or two specific embodiments.
Specific embodiment four:Unlike one of present embodiment and specific embodiment one to three:The step 2
Linear fit (corresponding is that cycle-index is identical) after taking the logarithm in five to the value of corresponding rescaled range R/S and length A,
Obtain Hurst indexes H(i);Process is:
Hurst establishes following relation by prolonged data preparation and practice summary:
Wherein c be constant, H(i)For Hurst indexes;Both members are taken the logarithm simultaneously:
Log (R/S)=H(i)log(A)+log(c) (10)
Found out by formula (10), for dividing for shape (self similarity) sequence, logA and log (R/S) is linear relationship, so as to profit
With least square method linear regression, the absolute value of fitting a straight line slope is required one-dimensional emitter Signals section SiHurst indexes H(i)。
Other steps and parameter are identical with one of specific embodiment one to three.
Specific embodiment five:Unlike one of present embodiment and specific embodiment one to four:The Euclid
Dimension E values are 1.
Other steps and parameter are identical with one of specific embodiment one to four.
Beneficial effects of the present invention are verified using following embodiment:
Embodiment one:Variance dimension emulation experiment compared with traditional characteristic performance parameters
The present embodiment compares variance Dimension Characteristics parameter and traditional characteristic performance parameters, mainly from discrimination angle point
Analysis.Specifically follow the steps below:
Step 1:Emitter Signals under the conditions of the multigroup different signal-to-noise ratio of emulation generation;
Step 2:Using it is according to the present invention based on the Emitter Fingerprint feature extracting method of variance dimension to step 1
The emitter Signals of generation carry out feature extraction, and obtaining feature vector, (segments n is 10, that is, the feature vector dimension generated is
10);
Step 3:Using the feature vector of step 2, categorized device calculates discrimination;
Step 4:The discrimination that step 3 is obtained is compared with the discrimination obtained using traditional characteristic parameter, and analysis is real
Test result.
The simulation experiment result is as shown in figure 4, show the effect and and general characteristics of variance Dimension Characteristics identification in figure
Parameter recognition effect[1](Bihl T J,Bauer K W,Temple M A.Feature Selection for RF
Fingerprinting With Multiple Discriminant Analysis and Using ZigBee Device
Emissions[J].IEEE Transactions on Information Forensics and Security,2016,11
(8):Comparison 1862-1874.).It can be seen that the discrimination curve of variance dimension is under general characteristics discrimination baseline,
This is because general characteristics discrimination baseline is the data in the case where dimension is 729, but this needs to consume a large amount of calculating costs.Therefore
Under conditions of considering recognition performance and calculating cost, the analysis of general characteristics performance often first reduces the dimension of feature vector, dimension
Number 50 can be used as a rational selection.
The experimental results showed that:Under the conditions of all signal-to-noise ratio set by emulation experiment, the variance dimension that dimension is 10 is special
Sign all achieves recognition performance more better than the general characteristics that dimension is 50;And under conditions of sample rate permission, pass through increase
Feature vector dimension can further improve recognition performance.On the whole, with the increase of signal-to-noise ratio, based on variance Dimension Characteristics
Discrimination gradually increase, and be maintained on general characteristics parameter discrimination.When signal-to-noise ratio is more than 13dB, tieed up based on variance
Several recognition correct rates is more than 95%.
The present invention can also have other various embodiments, without deviating from the spirit and substance of the present invention, this field
Technical staff makes various corresponding changes and deformation in accordance with the present invention, but these corresponding changes and deformation should all belong to
The protection domain of appended claims of the invention.
Claims (5)
1. a kind of Emitter Fingerprint feature extracting method based on variance dimension, it is characterised in that:The method detailed process is:
Step 1: carrying out segment processing to the one-dimensional emitter Signals received, one-dimensional emitter Signals section S is obtained1,S2,…
Si…,Sn;
SiFor i-th of one-dimensional emitter Signals section, 1≤i≤n, n values are positive integer;
Process is:
Step 1 one, the one-dimensional emitter Signals S that receives of input, according to the overall length of the one-dimensional emitter Signals S received and
Feature extraction segments demand, setting signal segment length w and segments n, each signal segment length are identical;
Step 1 two determines sliding window step delta, and sliding window step delta is obtained by formula (1):
WhereinRepresent positive bracket function;
Step 1 three, call function G=enframe (S, w, Δ), return value G, the G form for obtaining function enframe are as follows:
G is that a line number is n, and columns is the matrix of w;
Obtain one-dimensional emitter Signals section S1,S2,…Si…,Sn;
Step 2: carrying out variance Dimension Characteristics extraction to the one-dimensional emitter Signals section that step 1 obtains, variance dimension spy is obtained
Sign vector.
2. a kind of Emitter Fingerprint feature extracting method based on variance dimension according to claim 1, it is characterised in that:Institute
It states in step 2 and variance Dimension Characteristics extraction is carried out to the one-dimensional emitter Signals section that step 1 obtains, obtain variance Dimension Characteristics
Vector;Detailed process is:
Step 2 one, as i=1, by SiThe isometric subinterval that β initial length is A is divided into, with { dtRepresent in subinterval
Sequence;β values are positive integer;
Step 2 two, calculating add up the very poor of deviation in each subinterval;
For each subinterval, if:
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>A</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>p</mi>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, EqFor d in q-th of subintervalpAverage value, q=1,2 ... β, dpFor { dtIn p-th of element, Zt,qFor q-th of son
The accumulative deviation in section;
The very poor of all accumulative deviation values is obtained by formula (4):
Rq=max (Zt,q)-min(Zt,q) (4)
Step 2 three calculates each subinterval standard deviation:
If SqRepresent sequence { d in q-th of subintervaltStandard deviation, then:
<mrow>
<msub>
<mi>S</mi>
<mi>q</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>A</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>p</mi>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>/</mo>
<mi>A</mi>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
Sequence { dtRepresent sequence in each subinterval;
Step 2 four calculates each subinterval rescaled range, and averages;
The length for making subinterval is A changing with 2 for the exponential law at bottom as shown in formula (6),
Each subinterval rescaled range is calculated, and is averaged:
<mrow>
<mi>R</mi>
<mo>/</mo>
<mi>S</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>&beta;</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>&beta;</mi>
</munderover>
<msub>
<mi>R</mi>
<mi>q</mi>
</msub>
<mo>/</mo>
<msub>
<mi>S</mi>
<mi>q</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
Step 2 one is repeated to step 2 four, untilPerform step 2 five;
Step 2 five, take the logarithm to the value of corresponding rescaled range R/S and length A after linear fit, obtain Hurst indexes H(i);
Step 2 six utilizes Hurst indexes H(i)Calculate variance Dimension Characteristics vector
It is provided by formula (8):
<mrow>
<msubsup>
<mi>F</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mi>E</mi>
<mo>+</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein E represents Euclid dimensions;
Step 2 seven judges whether i is equal to n, if so,Terminate;If not, it performs
Step 2 one, until i is equal to n.
3. a kind of Emitter Fingerprint feature extracting method based on variance dimension according to claim 2, it is characterised in that:Institute
It states in step 2 one as i=1, by SiIt is divided into the isometric subinterval that β initial length is A;β values are positive integer;Process
For:
The one-dimensional emitter Signals section S that signal segment length for interception is wi, by one-dimensional emitter Signals section SiIt is divided at the beginning of β
Beginning length is the isometric subinterval of A, with sequence { dtRepresent sequence in each subinterval.
4. a kind of Emitter Fingerprint feature extracting method based on variance dimension according to claim 3, it is characterised in that:Institute
Linear fit after taking the logarithm in step 2 five to the value of corresponding rescaled range R/S and length A is stated, obtains Hurst indexes H(i);Process is:
Establish following relation:
<mrow>
<mi>R</mi>
<mo>/</mo>
<mi>S</mi>
<mo>=</mo>
<mi>c</mi>
<mo>&times;</mo>
<msup>
<mi>A</mi>
<msup>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</msup>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein c be constant, H(i)For Hurst indexes;Both members are taken the logarithm simultaneously:
Log (R/S)=H(i)log(A)+log(c) (10)
Found out by formula (10), log A and log (R/S) is linear relationship, and so as to utilize least square method linear regression, fitting is straight
The absolute value of line slope is required one-dimensional emitter Signals section SiHurst indexes H(i)。
5. a kind of Emitter Fingerprint feature extracting method based on variance dimension according to claim 4, it is characterised in that:Institute
Euclid dimension E values are stated as 1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711489288.3A CN108108712B (en) | 2017-12-29 | 2017-12-29 | Variance dimension-based radiation source fingerprint feature extraction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711489288.3A CN108108712B (en) | 2017-12-29 | 2017-12-29 | Variance dimension-based radiation source fingerprint feature extraction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108108712A true CN108108712A (en) | 2018-06-01 |
CN108108712B CN108108712B (en) | 2021-07-06 |
Family
ID=62215246
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711489288.3A Active CN108108712B (en) | 2017-12-29 | 2017-12-29 | Variance dimension-based radiation source fingerprint feature extraction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108108712B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109271855A (en) * | 2018-08-09 | 2019-01-25 | 深圳市永达电子信息股份有限公司 | A kind of feature extracting method of industry control signal |
CN110147848A (en) * | 2019-05-24 | 2019-08-20 | 哈尔滨工业大学 | A kind of radiation source personal feature Enhancement Method based on time-variable filtering theory |
CN111337975A (en) * | 2020-03-24 | 2020-06-26 | 吉林大学 | Automatic microseism event identification method based on variance fractal dimension |
CN111695444A (en) * | 2020-05-21 | 2020-09-22 | 哈尔滨工业大学 | Radiation source individual feature extraction method based on wave atomic transformation |
CN113343897A (en) * | 2021-06-25 | 2021-09-03 | 中国电子科技集团公司第二十九研究所 | Method for accelerating signal processing based on slope change of radiation signal |
CN114139579A (en) * | 2021-11-23 | 2022-03-04 | 吉林大学 | Deep learning microseism event detection method based on variance fractal |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012042393A (en) * | 2010-08-20 | 2012-03-01 | Toshiba Corp | Nondestructive deterioration diagnosis method for solder joint part |
CN103220055A (en) * | 2013-05-10 | 2013-07-24 | 厦门大学 | Multi-fractal gradient characteristic fingerprint identification method of wireless transmitter signal |
CN103226704A (en) * | 2013-05-10 | 2013-07-31 | 厦门大学 | Carrier and clock phase noise fingerprint characteristic combined identification method for wireless transmitter |
CN105426832A (en) * | 2015-11-12 | 2016-03-23 | 上海电机学院 | Communication radar radiation source identification method in presence of unsteady SNR (Signal Noise Ratio) |
CN106991393A (en) * | 2017-03-30 | 2017-07-28 | 上海电机学院 | Emitter Fingerprint feature extracting method based on Multifractal Dimension |
-
2017
- 2017-12-29 CN CN201711489288.3A patent/CN108108712B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012042393A (en) * | 2010-08-20 | 2012-03-01 | Toshiba Corp | Nondestructive deterioration diagnosis method for solder joint part |
CN103220055A (en) * | 2013-05-10 | 2013-07-24 | 厦门大学 | Multi-fractal gradient characteristic fingerprint identification method of wireless transmitter signal |
CN103226704A (en) * | 2013-05-10 | 2013-07-31 | 厦门大学 | Carrier and clock phase noise fingerprint characteristic combined identification method for wireless transmitter |
CN105426832A (en) * | 2015-11-12 | 2016-03-23 | 上海电机学院 | Communication radar radiation source identification method in presence of unsteady SNR (Signal Noise Ratio) |
CN106991393A (en) * | 2017-03-30 | 2017-07-28 | 上海电机学院 | Emitter Fingerprint feature extracting method based on Multifractal Dimension |
Non-Patent Citations (4)
Title |
---|
YIBING LI ET AL.: "The Application of Improved Generalized Fractal Box-counting Dimension Algorithm in Emitter Signals Recognition", 《JOURNAL OF INFORMATION & COMPUTATIONAL SCIENCE》 * |
余 沁 等: "基于排列熵分形维数特征提取的通信辐射源个体识别", 《空军预警学院学报》 * |
周斌: "信号细微特征提取及识别技术研究", 《中国优秀硕士学位论文 信息科技辑》 * |
李楠 等: "基于分形理论的辐射源识别算法", 《航天电子对抗》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109271855A (en) * | 2018-08-09 | 2019-01-25 | 深圳市永达电子信息股份有限公司 | A kind of feature extracting method of industry control signal |
CN109271855B (en) * | 2018-08-09 | 2021-07-09 | 深圳市永达电子信息股份有限公司 | Method for extracting characteristics of industrial control signals |
CN110147848A (en) * | 2019-05-24 | 2019-08-20 | 哈尔滨工业大学 | A kind of radiation source personal feature Enhancement Method based on time-variable filtering theory |
CN111337975A (en) * | 2020-03-24 | 2020-06-26 | 吉林大学 | Automatic microseism event identification method based on variance fractal dimension |
CN111695444A (en) * | 2020-05-21 | 2020-09-22 | 哈尔滨工业大学 | Radiation source individual feature extraction method based on wave atomic transformation |
CN111695444B (en) * | 2020-05-21 | 2023-06-23 | 哈尔滨工业大学 | Wave atom transformation-based radiation source individual feature extraction method |
CN113343897A (en) * | 2021-06-25 | 2021-09-03 | 中国电子科技集团公司第二十九研究所 | Method for accelerating signal processing based on slope change of radiation signal |
CN114139579A (en) * | 2021-11-23 | 2022-03-04 | 吉林大学 | Deep learning microseism event detection method based on variance fractal |
CN114139579B (en) * | 2021-11-23 | 2024-04-30 | 吉林大学 | Deep learning microseism event detection method based on variance fractal |
Also Published As
Publication number | Publication date |
---|---|
CN108108712B (en) | 2021-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108108712A (en) | A kind of Emitter Fingerprint feature extracting method based on variance dimension | |
CN108090462A (en) | A kind of Emitter Fingerprint feature extracting method based on box counting dimension | |
CN102891684B (en) | The construction method of the noisy reconfiguration system of a kind of compressed sensing | |
CN102749616B (en) | Fuzzy-clustering-based Aegis system signal sorting method | |
CN106130942A (en) | A kind of wireless communication signals Modulation Identification based on Cyclic Spectrum and method for parameter estimation | |
Liu | Multi-feature fusion for specific emitter identification via deep ensemble learning | |
CN107122643B (en) | Identity recognition method based on feature fusion of PPG signal and respiratory signal | |
CN110062379B (en) | Identity authentication method based on channel state information under human behavior scene | |
Ghadimi et al. | Deep learning-based approach for low probability of intercept radar signal detection and classification | |
CN107273825A (en) | Personal identification method is merged based on the physiological signal for improving canonical correlation analysis | |
CN104331583B (en) | A kind of Multifractal Modeling method based on Observed sea clutter | |
CN103632160A (en) | Combination-kernel-function RVM (Relevance Vector Machine) hyperspectral classification method integrated with multi-scale morphological characteristics | |
Shi | A Novel Ensemble Learning Algorithm Based on DS Evidence Theory for IoT Security. | |
CN107390259A (en) | A kind of nuclide identification method based on SVD and SVM | |
Yang et al. | Radar emitter signal recognition based on time-frequency analysis | |
CN106443625A (en) | Object identification method based on high-resolution one-dimensional image information fusion | |
CN105718915A (en) | Face identification method and system based on multi-visual-angle typical correlation analysis | |
Lv et al. | An algorithm of Iris feature-extracting based on 2D Log-Gabor | |
Limin et al. | Low probability of intercept radar signal recognition based on the improved AlexNet model | |
CN103076595B (en) | Abnormal type identifying method for multivariate discrete radar emitter signal | |
CN116383719A (en) | MGF radio frequency fingerprint identification method for LFM radar | |
CN114158004B (en) | Indoor passive moving target detection method based on graph convolution neural network | |
CN109541567A (en) | High-speed maneuver object detection method based on deep learning | |
CN114048774A (en) | Se-block-based resnet communication radiation source identification method and system | |
Bian et al. | LDA enhanced one-bit compressive sensing method for high-throughput mass spectrometry data feature selection |
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 |