CN106778549A - A kind of Harmonic Detection object detection algorithms and device based on multidimensional information fusion - Google Patents
A kind of Harmonic Detection object detection algorithms and device based on multidimensional information fusion Download PDFInfo
- Publication number
- CN106778549A CN106778549A CN201611078917.9A CN201611078917A CN106778549A CN 106778549 A CN106778549 A CN 106778549A CN 201611078917 A CN201611078917 A CN 201611078917A CN 106778549 A CN106778549 A CN 106778549A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- detection
- energy
- sense channel
- triple
- 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
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
The present invention relates to Harmonic Detection technology, for the problem that prior art is present, there is provided a kind of Harmonic Detection object detection algorithms and device based on multidimensional information fusion.Multidimensional information is obtained, from different dimension cascading judgements, the object resolving power of Harmonic Detection is improved, false-alarm probability is reduced, antijamming capability is improved, and then improve the actual handlability of Harmonic Detection.The present invention first crosses whether threshold point gross energy is 0 to discriminate whether with the presence of target according to each passage, if, show with the presence of target, threshold point number, second harmonic are then crossed according to second harmonic cross threshold point total energy value, triple-frequency harmonics and cross threshold point number, triple-frequency harmonics and cross threshold point total energy value, second harmonic energy to the contribution degree of total harmonic energy and the relation of second harmonic coefficient of kurtosis and triple-frequency harmonics coefficient of kurtosis, detect irregular many metal knot objects(Such as key), regular object based on metal junction characteristic(Such as metal-back electronic equipment), regular object based on semiconductor junction characteristic;Otherwise, then it is judged as without detection target, and exports without detection objective result.
Description
Technical field
Calculated the present invention relates to Harmonic Detection technology, especially a kind of Harmonic Detection object detection based on multidimensional information fusion
Method and device.
Background technology
Research shows that artificial metal's target and semiconductor targets have nonlinear characteristic, when the electromagnetism using CF
During ripple irradiation artificial metal's target and semiconductor targets in addition to incidence wave is scattered, can also be externally unexistent in incident ripple
Frequency component --- harmonic component.Its V-I characteristic can be expressed as exponential form for semiconductor, can produce abundant harmonic wave frequently
Rate, and the transition loss of lowest frequency harmonic wave is minimum, therefore its second-harmonic power is maximum.For artificial metal's target, its V-I
Characteristic is an antisymmetry curve, so it can only produce odd harmonic, and triple-frequency harmonics therein is most strong.
Because the harmonic generation efficiency of artificial metal's knot and semiconductor junction is relatively low, so Harmonic Detection is visited suitable for short range
Survey.Using the electromagnetic wave irradiation detected object of special frequency channel, target is found by receiving second harmonic and triple-frequency harmonics.Mesh
Before, conventional harmonic detecting method mainly has the harmonic detecting side based on instantaneous reactive power theory and based on Fourier transformation
Method.Harmonic detecting method principle based on instantaneous reactive power theory is simple, rapid dynamic response speed, but its accuracy of detection is filtered
The influence of ripple device is very big, while only entering line translation in time domain, is unfavorable for spectrum analysis;Harmonic detecting side based on Fourier transformation
Method, by docking collect mail number energy for carrying out Fourier transformation acquisition second harmonic and triple-frequency harmonics, the detection side based on energy
Formula principle is simple, and engineering realizability is high, but the method simply can only adjudicate target whether there is, and object resolving power is relatively low, and this is serious
Have impact on the practice of Harmonic Detection technology.
The content of the invention
The technical problems to be solved by the invention are:For the problem that prior art is present, there is provided a kind of based on multidimensional letter
Cease the Harmonic Detection object detection algorithms and device of fusion.Solution is based purely on the harmonic detecting method object resolving power of energy
Low, false alarm rate problem high, by taking multipath reception signal various processing modes, obtains multidimensional information, from different dimensions
Degree cascading judgement, improves Harmonic Detection object resolving power, reduces false-alarm probability, improves antijamming capability, and then improves harmonic wave and visits
The actual handlability surveyed.
The technical solution adopted by the present invention is as follows:
A kind of Harmonic Detection object detection algorithms based on multidimensional information fusion include:
Step 1:Detection zone is divided into N number of subregion, each region one group of transceiver channel of correspondence, N roads transceiver channel
It is separate;N number of transceiver channel produces the harmonic signal of N number of sense channel;
Step 2:First the harmonic signal to N number of sense channel is pre-processed, and then extracts N number of sense channel detection frequency range
Whether the spectrum information of harmonic signal, judged threshold point harmonic wave gross energy and be not 0 port number more than 0 according to spectrum information;
If threshold point harmonic wave gross energy is for 0 sense channel quantity is not 0;Then it is judged as without detection target, and exports without detection target
As a result;Otherwise, N number of sense channel is crossed into threshold point harmonic wave gross energy ESiDescending is arranged, i=1, flag=0, performs step 3:
Wherein i is less than or equal to N;
Step 3:Whether the harmonic wave gross energy of sense channel corresponding to harmonic wave gross energy row i-th bit is not 0 after rearrangement;If
It is then to perform step 4;Otherwise, then it is judged as without detection target, and exports without detection objective result:
Step 4:The second harmonic calculated according to spectrum information crosses threshold point total energy value, that triple-frequency harmonics crosses threshold point is total
The contribution degree and second harmonic coefficient of kurtosis and triple-frequency harmonics coefficient of kurtosis of energy value, second harmonic energy to total harmonic energy
Relation, detect irregular many metal knot objects (such as key), regular object (such as metal-back electronics based on metal junction characteristic
Equipment), the regular object based on semiconductor junction characteristic.
Further, the N roads transceiver channel course of work is:Object to be detected in N number of parallel subregion is subject to electricity
Magnetic wave irradiate when produce harmonic signal, harmonic signal to enter receiving channel by reception antenna, in receiving channel complete frequency dividing,
Intermediate frequency is changed to from radio frequency after filtering and down-converted, then is data signal by AD sample conversions;Data signal after filtering,
FFT is done after down coversion and extraction treatment, spectrum information is obtained, second-harmonic detection frequency range spectrum information and three are then extracted
Subharmonic detects frequency range spectrum information.
Further, each sense channel of statistics is crossed threshold point number and is counted by test of many times and determines a fixed threshold
Points Threshold, counts each passage second harmonic and triple-frequency harmonics detection frequency range crosses the points N 2 of thresholdingiAnd N3i, i represents logical
Taoist monastic name, i=1,2 ... N;
Calculate the gross energy ES2 that each passage second harmonic crosses threshold pointiWith the gross energy ES3 that triple-frequency harmonics crosses threshold pointi;
Calculate each passage second harmonic signal and harmonic signal crosses the energy and ES of threshold pointi, and each passage two
Contribution degree P2 of the subharmonic energy to harmonic wave gross energyi;
Calculate coefficient of kurtosis:Coefficient of kurtosis computing formula isWherein YiIt is sample measures value,For equal
Value, S is standard variance, and N is sample number.Second harmonic coefficient of kurtosis Pk2 is calculated respectivelyiWith triple-frequency harmonics coefficient of kurtosis Pk3i。
Further, step 4 detailed process is:
Step 41:Judge that the second harmonic of a certain sense channel crosses threshold point gross energy and triple-frequency harmonics crosses thresholding total energy
Whether amount is not 0;If not being 0, step 43 is performed;Otherwise, judge that the sense channel second harmonic crosses thresholding gross energy and is
It is no to be equal to 0, perform step 42;
Step 42:When the sense channel second harmonic crosses thresholding gross energy equal to 0, represent that object to be detected is metal knot
Then regular object based on characteristic, i=i+1 performs step 3;Otherwise, the sense channel second harmonic crosses thresholding gross energy not
Equal to 0, represent that object to be detected is irregular many metal objects, then i=i+1 performs step 3;
Step 43:Judge to cross threshold point number by the sense channel second harmonic whether cross threshold point number sum with triple-frequency harmonics big
In P;When more than P, then judge that second harmonic crosses whether threshold point gross energy crosses thresholding gross energy less than triple-frequency harmonics, performs step
Rapid 44;Otherwise, i=i+1, then performs step 3;P is more than or equal to 4;
Step 44:When the sense channel second harmonic crosses threshold point gross energy crosses thresholding gross energy less than triple-frequency harmonics,
Then judge the sense channel second harmonic to total harmonic energy contribution degree whether less than M and the sense channel and adjacency channel three times
Harmonic wave and second harmonic coefficient of kurtosis are all higher than H, perform step 47;Otherwise, step 45 is performed;M ∈ 0.3 to 0.5, H ∈ 80 is arrived
100
Step 45:Judge whether flag is 1;When flag is equal to 1, judge the sense channel second harmonic to total harmonic wave energy
Whether amount contribution degree is more than S, performs step 46;Otherwise, object to be detected be semiconductor junction characteristic based on regular object, i=i
+ 1, then perform step 3;S belongs to 0.6 to 0.8
Step 46:When the sense channel second harmonic is more than S to total harmonic energy contribution degree, then the object to be detected is
Then regular object based on semiconductor junction characteristic, i=i+1 performs step 3;Otherwise, i=i+1, then performs step 3;
Step 47:When the sense channel second harmonic energy to gross energy contribution degree less than M and the sense channel and adjacent
When passage triple-frequency harmonics and second harmonic coefficient of kurtosis are all higher than H, then the sense channel testing result is irregular many metal knots
Object, flag=1;Otherwise, judge whether flag is 1, perform step 48;
Step 48:When flag is 1, then i=i+1, then performs step 3;Otherwise, the sense channel object to be detected is
Then regular object based on metal junction characteristic, i=i+1 performs step 3.
A kind of Harmonic Detection article detection device based on multidimensional information fusion includes:
Signal pre-processing module, for detection zone to be divided into N number of subregion, the one group of transmitting-receiving of correspondence of each region is led to
Road, N roads transceiver channel is separate;N number of transceiver channel produces the harmonic signal of N number of sense channel;N number of transceiver channel produces N
The harmonic signal of individual sense channel;First the harmonic signal to N number of sense channel is pre-processed, and then extracts N number of sense channel
The spectrum information of frequency range harmonic signal is detected, judged threshold point harmonic wave gross energy not as 0 port number is according to spectrum information
It is no to be more than 0;If threshold point harmonic wave gross energy is for 0 sense channel quantity is not 0;Then it is judged as without detection target, and exports nothing
Detection objective result;Otherwise, N number of sense channel is crossed into threshold point harmonic wave gross energy ESiDescending is arranged, i=1, flag=0, so
Processed by Harmonic Detection object detection module afterwards, wherein i is less than or equal to N;
Harmonic Detection object detection module, for the harmonic wave of sense channel corresponding to harmonic wave gross energy row i-th bit after rearrangement
Whether gross energy is not 0;If so, the second harmonic for then being calculated according to spectrum information crosses threshold point total energy value, triple-frequency harmonics
Cross threshold point total energy value, second harmonic energy humorous to the contribution degree and second harmonic coefficient of kurtosis of total harmonic energy and three times
The relation of ripple coefficient of kurtosis, detects irregular many metal knot objects (such as key), the regular object based on metal junction characteristic (such as
Metal-back electronic equipment), the regular object based on semiconductor junction characteristic;Otherwise, then it is judged as without detection target, and exports without inspection
Survey objective result.
Further, the N roads transceiver channel course of work is:Object to be detected in N number of parallel subregion is subject to electricity
Magnetic wave irradiate when produce harmonic signal, harmonic signal to enter receiving channel by reception antenna, in receiving channel complete frequency dividing,
Intermediate frequency is changed to from radio frequency after filtering and down-converted, then is data signal by AD sample conversions;Data signal after filtering,
FFT is done after down coversion and extraction treatment, spectrum information is obtained, second-harmonic detection frequency range spectrum information and three are then extracted
Subharmonic detects frequency range spectrum information.
Further, each sense channel of statistics is crossed threshold point number and is counted by test of many times and determines a fixed threshold
Points Threshold, counts each passage second harmonic and triple-frequency harmonics detection frequency range crosses the points N 2 of thresholdingiAnd N3i, i represents logical
Taoist monastic name, i=1,2 ... N;
Calculate the gross energy ES2 that each passage second harmonic crosses threshold pointiWith the gross energy ES3 that triple-frequency harmonics crosses threshold pointi;
Calculate each passage second harmonic signal and harmonic signal crosses the energy and ES of threshold pointi, and each passage two
Contribution degree P2 of the subharmonic energy to harmonic wave gross energyi;
Calculate coefficient of kurtosis:Coefficient of kurtosis computing formula isWherein YiIt is sample measures value,For equal
Value, S is standard variance, and N is sample number.Second harmonic coefficient of kurtosis Pk2 is calculated respectivelyiWith triple-frequency harmonics coefficient of kurtosis Pk3i。
Further, the second harmonic calculated according to spectrum information crosses threshold point total energy value, triple-frequency harmonics mistake
The contribution degree and second harmonic coefficient of kurtosis and triple-frequency harmonics of threshold point total energy value, second harmonic energy to total harmonic energy
The relation of coefficient of kurtosis, detects irregular many metal knot objects (such as key), the regular object based on metal junction characteristic (such as gold
Category shell electronic equipment), the regular object detailed process based on semiconductor junction characteristic be:
Step 41:Judge that the second harmonic of a certain sense channel crosses threshold point gross energy and triple-frequency harmonics crosses thresholding total energy
Whether amount is not 0;If not being 0, step 43 is performed;Otherwise, judge that the sense channel second harmonic crosses thresholding gross energy and is
It is no to be equal to 0, perform step 42;
Step 42:When the sense channel second harmonic crosses thresholding gross energy equal to 0, represent that object to be detected is metal knot
Then regular object based on characteristic, i=i+1 performs step 3;Otherwise, the sense channel second harmonic crosses thresholding gross energy not
Equal to 0, represent that object to be detected is irregular many metal objects, then i=i+1 performs step 3;
Step 43:Judge to cross threshold point number by the sense channel second harmonic whether cross threshold point number sum with triple-frequency harmonics big
In P;When more than P, then judge that second harmonic crosses whether threshold point gross energy crosses thresholding gross energy less than triple-frequency harmonics, performs step
Rapid 44;Otherwise, i=i+1, then performs step 3;P is more than or equal to 4;
Step 44:When the sense channel second harmonic crosses threshold point gross energy crosses thresholding gross energy less than triple-frequency harmonics,
Then judge the sense channel second harmonic to total harmonic energy contribution degree whether less than M and the sense channel and adjacency channel three times
Harmonic wave and second harmonic coefficient of kurtosis are all higher than H, perform step 47;Otherwise, step 45 is performed;M ∈ 0.3 to 0.5, H ∈ 80 is arrived
100
Step 45:Judge whether flag is 1;When flag is equal to 1, judge the sense channel second harmonic to total harmonic wave energy
Whether amount contribution degree is more than S, performs step 46;Otherwise, object to be detected be semiconductor junction characteristic based on regular object, i=i
+ 1, then perform step 3;S belongs to 0.6 to 0.8
Step 46:When the sense channel second harmonic is more than S to total harmonic energy contribution degree, then the object to be detected is
Then regular object based on semiconductor junction characteristic, i=i+1 performs step 3;Otherwise, i=i+1, then performs step 3;
Step 47:When the sense channel second harmonic energy to gross energy contribution degree less than M and the sense channel and adjacent
When passage triple-frequency harmonics and second harmonic coefficient of kurtosis are all higher than H, then the sense channel testing result is irregular many metal knots
Object, flag=1;Otherwise, judge whether flag is 1, perform step 48;
Step 48:When flag is 1, then i=i+1, then performs step 3;Otherwise, the sense channel object to be detected is
Then regular object based on metal junction characteristic, i=i+1 performs step 3.
In sum, by adopting the above-described technical solution, the beneficial effects of the invention are as follows:
The present invention realizes Harmonic Detection using the method for multidimensional information fusion, according to the nonlinear characteristic of different type object
Existing difference, the detection and identification of object are realized from different dimensions.With the Harmonic Detection method that tradition is based purely on energy
Compare, the present invention has object resolution capability higher, can effectively exclude irregular many metal knot objects such as key such as and bring
Interference, and the classification of detected object is can recognize that, facilitate subsequent treatment, it is adapted to apply to the fortune higher to detection false alarm rate requirement
Use place.
Brief description of the drawings
Examples of the present invention will be described by way of reference to the accompanying drawings, wherein:
Fig. 1 is flow chart of the present invention.
Specific embodiment
All features disclosed in this specification, or disclosed all methods or during the step of, except mutually exclusive
Feature and/or step beyond, can combine by any way.
Any feature disclosed in this specification, unless specifically stated otherwise, can be equivalent or with similar purpose by other
Alternative features are replaced.I.e., unless specifically stated otherwise, each feature is an example in a series of equivalent or similar characteristics
.
Artificial many metal knot objects in irregular shape in the electromagnetic wave irradiation by certain frequency and energy, scattering
Harmonic signal occurs the phenomenon of video stretching, so whole spectral change can than shallower, and not only comprising odd harmonic but also
Comprising even-order harmonic;And general artificial metal's knot object and semiconductor junction object shine in the electromagnetic wave by same frequency and energy
Penetrate that time-frequency spectrum widening phenomenon is not obvious, spectral change is more precipitous.Coefficient of kurtosis is to characterize probability density distribution curve average
The characteristic of peak value height at value.Intuitively apparently, coefficient of kurtosis reflects the thickness of afterbody.Coefficient of kurtosis is bigger, and distribution just has
More extremums, then its residual value must be more concentrated on around mode, its Energy distribution is inevitable just more precipitous.Instead
It, Energy distribution is gentler.The present invention crosses threshold point number, mistake by calculating the coefficient of kurtosis combination harmonic wave that harmonic spectrum is distributed
Threshold point energy and and second harmonic the contribution degree of harmonic wave gross energy is detected and differentiates object.
1) spectrum information of detection frequency range is obtained
Fourier transformation is done to echo-signal, the spectrum distribution of second harmonic signal and harmonic signal is obtained, according to
The frequency and bandwidth selection special frequency channel of second harmonic signal and harmonic signal are used as detection frequency range;
2) multidimensional information is calculated
Second harmonic is calculated respectively and triple-frequency harmonics detects the coefficient of kurtosis of frequency range.Determine a fixed threshold by testing,
Statistics second-harmonic detection frequency range and triple-frequency harmonics detection frequency range cross threshold point number, cross threshold point gross energy respectively, calculate secondary
Contribution degree of the harmonic signal to total harmonic signal energy.
3) multidimensional information fusion judgement
Threshold point gross energy is crossed by harmonic wave to judge whether to detect target;Multidimensional information is combined if target is detected
Target type is differentiated;Selectively alarmed according to concrete application scene.
Fig. 1 is multidimensional information fusion harmonic detecting and object resolution procedure figure, and by taking the sense channel of N roads as an example, each passage is humorous
Ripple signal is first pre-processed and FFT, then extracts detection frequency range, counts multidimensional information.Each Air conduct measurement result is initial
Value assigns 0, indicates without detection target.During judgement, first determine whether whether there is detection target, terminate if being adjudicated without if, otherwise by each passage
Threshold point harmonic wave gross energy descending is crossed successively to crossing threshold point harmonic wave gross energy not for 0 passage makes decisions.Tu Zhong branches 1 are sentenced
Certainly whether detection target is irregular many metal knot objects, and testing result assigns 1.In branch 2, when second harmonic crosses thresholding gross energy
Represent that detection object is the regular object based on metal junction characteristic during equal to 0, testing result assigns 3, is otherwise semiconductor junction characteristic
Based on regular object, testing result assign 2.Branch 3 is the passage judgement that characteristic of semiconductor is better than metal junction characteristic, if flag=
0 to indicate and go out irregular many metal knot objects without Air conduct measurement, and testing result assigns 2;Otherwise indicate that Air conduct measurement goes out irregularly
Many metal knot objects, represent that this channel semiconductor characteristic is strong, detection knot if second harmonic is more than 0.7 to the contribution of total harmonic energy
Fruit assigns 2, otherwise it is assumed that the passage is without detection target, the harmonic signal for detecting is by the irregular many metal knot objects of other passages
Influence.Branch 4 is that metal junction characteristic is better than semiconductor junction characteristic but is not belonging to sentencing for irregular many metal knot object channels
Certainly, if flag=1, then it is assumed that the harmonic signal of the passage is influenceed by the irregular many metal knot objects of other passages, is detected
Result assigns 0;Otherwise it is assumed that being the regular object based on metal junction characteristic, testing result assigns 3.
Embodiment one:
Step 1:Spectrum information is extracted;
Whole detection zone is divided into N number of subregion by the system, each region one group of transceiver channel of correspondence, the transmitting-receiving of N roads
Passage is separate, is taken turns to operate by certain scanning sequency.When system works, metal knot or semiconductor junction are subject to CF
Electromagnetic wave irradiation, so as to scatter second harmonic signal and harmonic signal, the received antenna of harmonic signal enters radio frequency
End receiving channel, after completion frequency dividing, filtering and down-converted in receiving channel from rf conversion to intermediate frequency, then adopts by AD
Sample is changed into data signal.Data signal enters numerical control plate, and filtering, down coversion and extraction treatment are realized in numerical control plate, so that
Data volume and throughput are reduced, FFT is then done, spectrum information is obtained.Second harmonic and triple-frequency harmonics detection frequency range are extracted,
Retain standby.
Step 2:Multidimensional information is calculated
1) count each sense channel and cross threshold point number.Counted by test of many times and determine a fixed threshold Threshold,
Count each passage second harmonic and triple-frequency harmonics detection frequency range crosses the points N 2 of thresholdingiAnd N3i, i represents channel number, i=1,2 ...
N;
2) the gross energy ES2 that each passage second harmonic crosses threshold point is calculatediWith the gross energy that triple-frequency harmonics crosses threshold point
ES3i;
3) calculate each passage second harmonic signal and harmonic signal crosses the energy and ES of threshold pointi, and each passage
Contribution degree P2 of the second harmonic energy to harmonic wave gross energyi;
4) coefficient of kurtosis is calculated.Coefficient of kurtosis computing formula isWherein YiIt is sample measures value,For
Average, S is standard variance, and N is sample number.Second harmonic coefficient of kurtosis Pk2 is calculated respectivelyiWith triple-frequency harmonics coefficient of kurtosis Pk3i。
Step 3:Multidimensional information fusion, object detection and identification.
Each passage is crossed into threshold point harmonic wave gross energy ESiCarry out descending arrangement.Judgement flow according to Fig. 1, according to ESiDrop
Channel sequence after sequence arrangement makes decisions to each passage successively.Object to be detected is divided into three major types by court verdict, respectively
For the regular object (such as metal-back electronic equipment) based on irregular many metal knot objects (such as key), metal junction characteristic, partly lead
Regular object (such as plastic housing electronic equipment) based on body junction characteristic.Phase is made to testing result finally according to practice scene
Should process.
The invention is not limited in foregoing specific embodiment.The present invention is expanded to and any in this manual disclosed
New feature or any new combination, and disclose any new method or process the step of or any new combination.
Claims (8)
1. a kind of Harmonic Detection object detection algorithms based on multidimensional information fusion, it is characterised in that including:
Step 1:Detection zone is divided into N number of subregion, each region one group of transceiver channel of correspondence, N roads transceiver channel is mutual
It is independent;N number of transceiver channel produces the harmonic signal of N number of sense channel;
Step 2:First the harmonic signal to N number of sense channel is pre-processed, and then extracts N number of sense channel detection frequency range harmonic wave
Whether the spectrum information of signal, judged threshold point harmonic wave gross energy and be not 0 port number more than 0 according to spectrum information;If door
Point of accumulation harmonic wave gross energy is for 0 sense channel quantity is not 0;Then it is judged as without detection target, and exports without detection objective result;
Otherwise, N number of sense channel is crossed into threshold point harmonic wave gross energy ESiDescending is arranged, i=1, flag=0, performs step 3:Wherein i
Less than or equal to N;
Step 3:Whether the harmonic wave gross energy of sense channel corresponding to harmonic wave gross energy row i-th bit is not 0 after rearrangement;If so, then
Perform step 4;Otherwise, then it is judged as without detection target, and exports without detection objective result:
Step 4:The second harmonic calculated according to spectrum information crosses threshold point total energy value, triple-frequency harmonics and crosses threshold point gross energy
Value, second harmonic energy are to the contribution degree of total harmonic energy and the pass of second harmonic coefficient of kurtosis and triple-frequency harmonics coefficient of kurtosis
System, detects regular object, the regular thing based on semiconductor junction characteristic based on irregular many metal knot objects, metal junction characteristic
Body.
2. a kind of Harmonic Detection object detection algorithms based on multidimensional information fusion according to claim 1, its feature exists
It is in the N roads transceiver channel course of work:Object to be detected in N number of parallel subregion is subject to be produced during electromagnetic wave irradiation humorous
Ripple signal, harmonic signal enters receiving channel by reception antenna, and frequency dividing, filtering and down-converted are completed in receiving channel
Intermediate frequency is changed to from radio frequency afterwards, then is data signal by AD sample conversions;Data signal after filtering, down coversion and extraction treatment
After do FFT, obtain spectrum information, then extract second-harmonic detection frequency range spectrum information and triple-frequency harmonics detection frequency range frequently
Spectrum information.
3. a kind of Harmonic Detection object detection algorithms based on multidimensional information fusion according to claim 2, its feature exists
Cross threshold point number and counted by test of many times in each sense channel of statistics and determine a fixed threshold points Threshold,
Count each passage second harmonic and triple-frequency harmonics detection frequency range crosses the points N 2 of thresholdingiAnd N3i, i represents channel number, i=1,2 ...
N;
Calculate the gross energy ES2 that each passage second harmonic crosses threshold pointiWith the gross energy ES3 that triple-frequency harmonics crosses threshold pointi;
Calculate each passage second harmonic signal and harmonic signal crosses the energy and ES of threshold pointi, and each passage second harmonic
Contribution degree P2 of the energy to harmonic wave gross energyi;
Calculate coefficient of kurtosis:Coefficient of kurtosis computing formula isWherein YiIt is sample measures value,It is average, S
It is standard variance, N is sample number;
Second harmonic coefficient of kurtosis Pk2 is calculated respectivelyiWith triple-frequency harmonics coefficient of kurtosis Pk3i。
4. a kind of Harmonic Detection object detection algorithms based on multidimensional information fusion according to claim 2, its feature exists
It is in step 4 detailed process:
Step 41:Judge that the second harmonic of a certain sense channel crosses threshold point gross energy and triple-frequency harmonics crosses thresholding gross energy and is
No is not 0;If not being 0, step 43 is performed;Otherwise, judge the sense channel second harmonic cross thresholding gross energy whether etc.
In 0, step 42 is performed;
Step 42:When the sense channel second harmonic crosses thresholding gross energy equal to 0, represent that object to be detected is metal junction characteristic
Based on regular object, i=i+1, then perform step 3;Otherwise, the sense channel second harmonic is crossed thresholding gross energy and is not equal to
0, represent that object to be detected is irregular many metal objects, then i=i+1 performs step 3;
Step 43:Judge that whether the sense channel second harmonic crosses threshold point number and crosses threshold point number sum more than P with triple-frequency harmonics;
When more than P, then judge that second harmonic crosses whether threshold point gross energy crosses thresholding gross energy less than triple-frequency harmonics, performs step
44;Otherwise, i=i+1, then performs step 3;P is more than or equal to 4;
Step 44:When the sense channel second harmonic crosses threshold point gross energy crosses thresholding gross energy less than triple-frequency harmonics, then sentence
Whether the disconnected sense channel second harmonic is less than M and the sense channel and adjacency channel triple-frequency harmonics to total harmonic energy contribution degree
H is all higher than with second harmonic coefficient of kurtosis, step 47 is performed;Otherwise, step 45 is performed;M ∈ 0.3 to 0.5, H ∈ 80 to 100;
Step 45:Judge whether flag is 1;When flag is equal to 1, judge the sense channel second harmonic to total harmonic energy tribute
Whether degree of offering is more than S, performs step 46;Otherwise, object to be detected be semiconductor junction characteristic based on regular object, i=i+1,
Then step 3 is performed;S belongs to 0.6 to 0.8
Step 46:When the sense channel second harmonic is more than S to total harmonic energy contribution degree, then the object to be detected is partly to lead
Then regular object based on body junction characteristic, i=i+1 performs step 3;Otherwise, i=i+1, then performs step 3;
Step 47:When the sense channel second harmonic energy is less than M and the sense channel and adjacency channel to gross energy contribution degree
When triple-frequency harmonics and second harmonic coefficient of kurtosis are all higher than H, then the sense channel testing result is irregular many metal knot objects,
Flag=1;Otherwise, judge whether flag is 1, perform step 48;
Step 48:When flag is 1, then i=i+1, then performs step 3;Otherwise, the sense channel object to be detected is metal
Then regular object based on junction characteristic, i=i+1 performs step 3.
5. a kind of Harmonic Detection article detection device based on multidimensional information fusion, it is characterised in that including:
Signal pre-processing module, for detection zone to be divided into N number of subregion, each region one group of transceiver channel of correspondence, N roads
Transceiver channel is separate;N number of transceiver channel produces the harmonic signal of N number of sense channel;N number of transceiver channel produces N number of detection
The harmonic signal of passage;First the harmonic signal to N number of sense channel is pre-processed, and then extracts N number of sense channel detection frequency
The spectrum information of section harmonic signal, threshold point harmonic wave gross energy was judged not as whether 0 port number is more than according to spectrum information
0;If threshold point harmonic wave gross energy is for 0 sense channel quantity is not 0;Then it is judged as without detection target, and exports without detection mesh
Mark result;Otherwise, N number of sense channel is crossed into threshold point harmonic wave gross energy ESiDescending is arranged, and then i=1, flag=0 pass through
Harmonic Detection object detection module is processed, and wherein i is less than or equal to N;
Harmonic Detection object detection module, for the harmonic wave total energy of sense channel corresponding to harmonic wave gross energy row i-th bit after rearrangement
Whether amount is not 0;If so, the second harmonic for then being calculated according to spectrum information crosses threshold point total energy value, triple-frequency harmonics moving into one's husband's household upon marriage
The contribution degree and second harmonic coefficient of kurtosis and triple-frequency harmonics peak of point of accumulation total energy value, second harmonic energy to total harmonic energy
The relation of state coefficient, detects that the regular object based on irregular many metal knot objects, metal junction characteristic, semiconductor junction characteristic are
Main regular object;Otherwise, then it is judged as without detection target, and exports without detection objective result.
6. a kind of Harmonic Detection article detection device based on multidimensional information fusion according to claim 5, its feature exists
It is in the N roads transceiver channel course of work:Object to be detected in N number of parallel subregion is subject to be produced during electromagnetic wave irradiation humorous
Ripple signal, harmonic signal enters receiving channel by reception antenna, and frequency dividing, filtering and down-converted are completed in receiving channel
Intermediate frequency is changed to from radio frequency afterwards, then is data signal by AD sample conversions;Data signal after filtering, down coversion and extraction treatment
After do FFT, obtain spectrum information, then extract second-harmonic detection frequency range spectrum information and triple-frequency harmonics detection frequency range frequently
Spectrum information.
7. a kind of Harmonic Detection article detection device based on multidimensional information fusion according to claim 6, its feature exists
Cross threshold point number and counted by test of many times in each sense channel of statistics and determine a fixed threshold points Threshold,
Count each passage second harmonic and triple-frequency harmonics detection frequency range crosses the points N 2 of thresholdingiAnd N3i, i represents channel number, i=1,2 ...
N;
Calculate the gross energy ES2 that each passage second harmonic crosses threshold pointiWith the gross energy ES3 that triple-frequency harmonics crosses threshold pointi;
Calculate each passage second harmonic signal and harmonic signal crosses the energy and ES of threshold pointi, and each passage second harmonic
Contribution degree P2 of the energy to harmonic wave gross energyi;
Calculate coefficient of kurtosis:Coefficient of kurtosis computing formula isWherein YiIt is sample measures value,It is average, S
It is standard variance, N is sample number;Second harmonic coefficient of kurtosis Pk2 is calculated respectivelyiWith triple-frequency harmonics coefficient of kurtosis Pk3i。
8. a kind of Harmonic Detection article detection device based on multidimensional information fusion according to claim 6, its feature exists
In the second harmonic calculated according to spectrum information cross threshold point total energy value, triple-frequency harmonics cross threshold point total energy value,
Second harmonic energy to the contribution degree of total harmonic energy and the relation of second harmonic coefficient of kurtosis and triple-frequency harmonics coefficient of kurtosis,
Detect irregular many metal knot objects (such as key), the regular object (such as metal-back electronic equipment) based on metal junction characteristic,
Regular object detailed process based on semiconductor junction characteristic is:
Step 41:Judge that the second harmonic of a certain sense channel crosses threshold point gross energy and triple-frequency harmonics crosses thresholding gross energy and is
No is not 0;If not being 0, step 43 is performed;Otherwise, judge the sense channel second harmonic cross thresholding gross energy whether etc.
In 0, step 42 is performed;
Step 42:When the sense channel second harmonic crosses thresholding gross energy equal to 0, represent that object to be detected is metal junction characteristic
Based on regular object, i=i+1, then perform step 3;Otherwise, the sense channel second harmonic is crossed thresholding gross energy and is not equal to
0, represent that object to be detected is irregular many metal objects, then i=i+1 performs step 3;
Step 43:Judge that whether the sense channel second harmonic crosses threshold point number and crosses threshold point number sum more than P with triple-frequency harmonics;
When more than P, then judge that second harmonic crosses whether threshold point gross energy crosses thresholding gross energy less than triple-frequency harmonics, performs step
44;Otherwise, i=i+1, then performs step 3;P is more than or equal to 4;
Step 44:When the sense channel second harmonic crosses threshold point gross energy crosses thresholding gross energy less than triple-frequency harmonics, then sentence
Whether the disconnected sense channel second harmonic is less than M and the sense channel and adjacency channel triple-frequency harmonics to total harmonic energy contribution degree
H is all higher than with second harmonic coefficient of kurtosis, step 47 is performed;Otherwise, step 45 is performed;M ∈ 0.3 to 0.5, H ∈ 80 to 100
Step 45:Judge whether flag is 1;When flag is equal to 1, judge the sense channel second harmonic to total harmonic energy tribute
Whether degree of offering is more than S, performs step 46;Otherwise, object to be detected be semiconductor junction characteristic based on regular object, i=i+1,
Then step 3 is performed;S belongs to 0.6 to 0.8
Step 46:When the sense channel second harmonic is more than S to total harmonic energy contribution degree, then the object to be detected is partly to lead
Then regular object based on body junction characteristic, i=i+1 performs step 3;Otherwise, i=i+1, then performs step 3;
Step 47:When the sense channel second harmonic energy is less than M and the sense channel and adjacency channel to gross energy contribution degree
When triple-frequency harmonics and second harmonic coefficient of kurtosis are all higher than H, then the sense channel testing result is irregular many metal knot objects,
Flag=1;Otherwise, judge whether flag is 1, perform step 48;
Step 48:When flag is 1, then i=i+1, then performs step 3;Otherwise, the sense channel object to be detected is metal
Then regular object based on junction characteristic, i=i+1 performs step 3.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611078917.9A CN106778549B (en) | 2016-11-30 | 2016-11-30 | Harmonic detection object detection algorithm and device based on multi-dimensional information fusion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611078917.9A CN106778549B (en) | 2016-11-30 | 2016-11-30 | Harmonic detection object detection algorithm and device based on multi-dimensional information fusion |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106778549A true CN106778549A (en) | 2017-05-31 |
CN106778549B CN106778549B (en) | 2020-04-21 |
Family
ID=58898917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611078917.9A Active CN106778549B (en) | 2016-11-30 | 2016-11-30 | Harmonic detection object detection algorithm and device based on multi-dimensional information fusion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106778549B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113498482A (en) * | 2020-10-15 | 2021-10-12 | 深圳市安卫普科技有限公司 | Scanning type nonlinear node detection method and device |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013138640A1 (en) * | 2012-03-16 | 2013-09-19 | Newport Corporation | Continuous wave ultraviolet laser based on stimulated raman scattering |
CN104849764A (en) * | 2015-05-26 | 2015-08-19 | 上海海潮新技术研究所 | Array non-linear object detection system and method |
CN105116455A (en) * | 2015-09-24 | 2015-12-02 | 长春市泽安科技有限公司 | Detector for composite metals and non-linear nodes |
CN106054177A (en) * | 2016-05-24 | 2016-10-26 | 广东电网有限责任公司信息中心 | Electronic equipment detection device based on semiconductor component detection |
-
2016
- 2016-11-30 CN CN201611078917.9A patent/CN106778549B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013138640A1 (en) * | 2012-03-16 | 2013-09-19 | Newport Corporation | Continuous wave ultraviolet laser based on stimulated raman scattering |
CN104849764A (en) * | 2015-05-26 | 2015-08-19 | 上海海潮新技术研究所 | Array non-linear object detection system and method |
CN105116455A (en) * | 2015-09-24 | 2015-12-02 | 长春市泽安科技有限公司 | Detector for composite metals and non-linear nodes |
CN106054177A (en) * | 2016-05-24 | 2016-10-26 | 广东电网有限责任公司信息中心 | Electronic equipment detection device based on semiconductor component detection |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113498482A (en) * | 2020-10-15 | 2021-10-12 | 深圳市安卫普科技有限公司 | Scanning type nonlinear node detection method and device |
WO2022077328A1 (en) * | 2020-10-15 | 2022-04-21 | 深圳市安卫普科技有限公司 | Scanning-type non-linear node detection method and apparatus |
Also Published As
Publication number | Publication date |
---|---|
CN106778549B (en) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107358250B (en) | Body gait recognition methods and system based on the fusion of two waveband radar micro-doppler | |
CN104158611B (en) | Wireless signal Interference Detection system and method based on spectrum analysis | |
CN101975940B (en) | Segmentation combination-based adaptive constant false alarm rate target detection method for SAR image | |
CN106153340B (en) | A kind of Fault Diagnosis of Roller Bearings | |
CN104135327B (en) | frequency spectrum sensing method based on support vector machine | |
CN108564005A (en) | A kind of human body tumble discrimination method based on convolutional neural networks | |
CN101299237B (en) | High spectroscopic data supervision classifying method based on information quantity dimensionality sequence | |
CN103399297A (en) | Machine learning based ultra-wideband NLOS (non line of sight) identification method | |
Jadhav et al. | A fruit quality management system based on image processing | |
CN103134765A (en) | Chinese medicine sample authenticity preliminary screening method based on terahertz time-domain spectrum | |
CN110187378B (en) | Neutron detector for identifying neutrons and gammas | |
CN106291275A (en) | A kind of local discharge superhigh frequency single waveform frequency domain character extracts and recognition methods | |
CN110514975A (en) | A kind of shelf depreciation genre classification methods of GIS device | |
CN106772666B (en) | A kind of new violation electronic equipment detection system and detection method | |
Gao et al. | Fast detecting and locating groups of targets in high-resolution SAR images | |
CN103675610A (en) | Method for extracting characteristic factors in online local discharge detection | |
Conning et al. | Analysis of measured radar data for specific emitter identification | |
CN102129685B (en) | Method for detecting irregular circle based on Gauss pyramid decomposition | |
CN102208017A (en) | Ship detection method based on high-resolution synthetic aperture radar images | |
CN102279146A (en) | Blood cell five classification method based on laser sheath flow technology | |
CN107067039A (en) | SAR image Ship Target quick determination method based on super-pixel | |
CN103020592A (en) | Remote sensing image ship detection false-alarm rejecting method based on local context | |
CN109886319A (en) | A kind of concealment object detecting method and system based on millimeter wave safety check image | |
CN106778549A (en) | A kind of Harmonic Detection object detection algorithms and device based on multidimensional information fusion | |
CN105022041A (en) | A towing interference sub-type identification method based on a empirical mode decomposition feature library |
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 |