CN106155983B - A kind of Correlation Identification detection system and its method based on FFT - Google Patents
A kind of Correlation Identification detection system and its method based on FFT Download PDFInfo
- Publication number
- CN106155983B CN106155983B CN201610556801.5A CN201610556801A CN106155983B CN 106155983 B CN106155983 B CN 106155983B CN 201610556801 A CN201610556801 A CN 201610556801A CN 106155983 B CN106155983 B CN 106155983B
- Authority
- CN
- China
- Prior art keywords
- fft
- frequency domain
- correlation
- module
- unit
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The Correlation Identification detection system based on FFT that the present invention relates to a kind of, including：First m-sequence generation module, the first the earth system, first time synchronization module, the first signal receiving module, the first FFT auto-correlations module, the first frequency domain conversion module, the first FFT cross-correlation modules, the second frequency domain conversion module and the first system function solve module；The first time synchronization module is connected with the first m-sequence generation module and the first signal receiving module, and the first m-sequence generation module is connected with the first FFT auto-correlation modules, and the first FFT auto-correlations module is connected with the first frequency domain conversion module；First signal receiving module is connected with the first FFT cross-correlation modules, and the first FFT cross-correlation modules are connected with the second frequency domain conversion module；First frequency domain conversion module and the second frequency domain conversion module solve module with the first system function respectively and are connected.The present invention, which has, recognizes the high characteristic of efficient, high-precision, noiseproof feature.
Description
Technical field
The present invention relates to a kind of Correlation Identification detection system and its method based on FFT.
Background technology
During the system function of detection system, especially when recognizing larger system (such as big ground system), often
It needs to acquire a large amount of data progress operation.Different operands is just had using different schemes, is usually distinguished in calculating correlation
It is calculated using the definition of auto-correlation and cross-correlation function when the correlation of knowledge, summation is used to be calculated.It is cumulative
During calculating, due to needing constantly to keep in, be multiplied and be added, computational efficiency substantially reduces.When calculating related,
There is a kind of quickly and effectively method, is exactly to calculate correlation using FFT by correlation theorem.First input and output signal
Fast Fourier Transform (FFT) is done, input signal is then divided into two groups, the input signal after one of which transformation takes conjugation, another group
Retain transformation results.It finds out that cycle is relevant according to correlation theorem as a result, last another mistake switches back to time domain, obtains relevant knot
Fruit.It is seen on surface, is transformed into frequency domain and has converted back more than time domain a step again, but due to the high efficiency of FFT calculating so that meter
The efficiency of calculation greatly improves.
Currently, the existing system function that Correlation Identification technology is applied to the earth recognizes field.Including being pressed down with Correlation Identification
Noise processed, the method to achieve the purpose that identification system function.But applied to the earth recognize when, often because
It is excessive so that the system seems time-consuming and laborious for data volume.One of them important reason is exactly because calculating relevant mistake
Cheng Zhong is calculated using cumulative definition.And inventor has found that the effect of calculating can be greatly improved when calculating related with FFT
Rate, thus get down to this method within the system.Practice and emulation prove that this method can ensure really and former result phase
With under the premise of, the efficiency of calculating is improved.
Invention content
The technology of the present invention solves the problems, such as：It solves the problems, such as that the Correlation Identification efficiency under big data quantity is low, remains phase
Close the characteristics of recognizing original effective inhibition noise.
Technical solution of the invention is：A kind of Correlation Identification detection system based on FFT, including：First m-sequence is given birth to
At module, the first the earth system, first time synchronization module, the first signal receiving module, the first FFT- auto-correlations module, first
Frequency domain conversion module, the first FFT- cross-correlation modules, the second frequency domain conversion module and the first system function solve module, wherein the
One FFT- auto-correlations module, the first frequency domain conversion module, the first FFT- cross-correlation modules, the first frequency domain conversion module and the first system
System function solves module and belongs to data processing section；The first time synchronization module and the first m-sequence generation module and first
Signal receiving module is connected, and the first m-sequence generation module is connected with the first FFT- auto-correlation modules again, the first FFT- auto-correlation moulds
Block is connected with the first frequency domain conversion module；First signal receiving module is connected with the first cross-correlation module, the first FFT- cross-correlation
Module is connected with the second frequency domain conversion module；First frequency domain conversion module and the second frequency domain conversion module solve mould with system function
Block is connected；
The function of the first m-sequence generation module is to be sent to the earth system after generating m-sequence according to setup parameter；It is first
First first time synchronization module setting predetermined instant makes the m of the first m-sequence generation module generation relevant parameter at the time of setting
Sequence, after generating m-sequence, m-sequence is just sent to the earth system by electrode.In this way, just completing generation and the hair of m-sequence
It send；The voltage data packet x of transmission can be collected simultaneously；
The function of first signal receiving module is to complete the reception through current signal too greatly.According at the first time
The reservation moment of synchronization module setting, the first signal receiving module are received through voltage signal too greatly, to form data packet
y；
The first FFT- auto-correlation modules include the first FFT frequency domain converting units, the first multiplication unit, the first IFFT changes
Change unit and the first data cutout unit.Wherein the first FFT frequency domain converting units are connected with the first multiplication unit, the first multiplication list
Member is connected with the first IFFT converter units, and the first IFFT units and the first data cutout unit are connected.First m-sequence generation module
After being connected with the first FFT frequency domain converting units, by the first FFT- auto-correlation modules, auto-correlation has been calculated with the method for frequency domain
Autocorrelative frequency-domain result, intercept out it is therein be lap, obtained accurate time domain autocorrelation result；
The first frequency domain conversion module includes the first FFT operators unit and frequency domain conversion results.First FFT operator units
It is connected respectively with the first data cutout unit and the first auto-correlation frequency domain unit.First data cutout unit sends first to
The first auto-correlation frequency domain unit is obtained after FFT operator units.In this way, just having obtained autocorrelative frequency-domain result.
The first FFT- cross-correlation modules include the 2nd FFT frequency domain converting units, the second multiplication unit, the 2nd IFFT changes
Change unit and the second data cutout unit.2nd FFT frequency domain converting units are connected with the second multiplication unit, the second multiplication unit with
2nd IFFT converter units, the 2nd IFFT converter units are connected with the second data cutout unit.Second frequency-domain transform unit simultaneously
It is connected with the second multiplication unit, the first signal receiving module is connected with the 2nd FFT frequency-domain transform units.Second signal receiving module
Obtained received data packet y by the 2nd FFT frequency domain converting units obtain frequency domain as a result, being then transferred to the second multiplication list
Member, while result is also transferred to the second multiplication unit by the second frequency-domain transform unit, is exactly the frequency domain of cross-correlation after the two multiplication
As a result, then obtaining the transformation results of time domain by the 2nd IFFT converter units, finally intercepted out with the second data cutout unit
The underlapped accurate cross correlation results of time domain.
The second frequency domain conversion module includes the 2nd FFT operators unit and the first cross-correlation frequency domain unit.2nd FFT is calculated
Subelement is connected with the 2nd FFT frequency domain converting units and the first cross-correlation frequency domain unit respectively.2nd FFT frequency domain converting units pass through
It crosses the 2nd FFT operator units and obtains the first cross-correlation frequency domain unit, the frequency-domain result of cross-correlation has just been calculated in this.
It includes the first divider and the first the earth frequency domain identification unit that system function, which solves module,.First divider point
It is not connected with the first the earth frequency domain identification unit with the first cross-correlation frequency domain unit, the first auto-correlation frequency domain unit.First mutually
It closes frequency domain unit and the first auto-correlation frequency domain unit and passes through the first divider, the first cross-correlation frequency domain unit divided by first is from phase
After closing frequency domain unit, the first the earth frequency domain identification unit is obtained, the first the earth frequency domain identification unit is as a result ultimately stored on.
The relevant treatment work of the Correlation Identification of detection system is completed in frequency domain by FFT operators, operation efficiency pole
Height is particularly suitable for using when the data volume of correlation computations is larger.The first FFT- auto-correlations module, the first FFT- are mutual
Close module solve detection system system function when the specific implementation process is as follows：First, transmitting terminal collects data packet x
After receiving terminal all collects data packet y simultaneously, data are given respectively the first FFT- auto-correlations module and the first FFT- cross-correlation
Module, after the first FT- auto-correlation moulds the first FFT frequency-domain transform units in the block of F obtain the result X (w) of frequency domain, transmitting terminal
Data packet x consign to the first multiplication unit, at this moment, multiplication unit can carry out conjugation processing to one of which data, then borrow
Correlation theorem is helped to calculate the autocorrelative frequency-domain result of cycle.Calculation formula is： To follow
The autocorrelative result of ring.By correlation theorem it is recognised that cycle auto-correlation and the equal condition of cross-correlation areWherein M refers to that M dot cycles are related, LxRefer to the points of x；Similarly, it recycles
The equal condition of cross-correlation and cross-correlation isWherein M is M dot cycle cross-correlation, Ly
For the points of data packet y.According to this condition, non-overlapping correlated results is intercepted in data cutout unit, that is,
Linearly related result.Auto-correlation and cross-correlation thus are calculated with FFT, then according to the frequency domain identification of Correlation Identification public affairs
Formula calculates the system function of the earth.
The default setting of first time synchronization module in the case of no modification is：12:01, i.e. predetermined instant 12:01 hair
Sending end and receiving terminal work at the same time.
The advantages of the present invention over the prior art are that：
(1) this computational efficiency that relevant Correlation Identification is calculated based on FFT is apparently higher than the cumulative method of traditional use
Correlation is calculated, when big data quantity, operational efficiency is very high, can greatly save the identification time.
(2) when identification system, noise can be inhibited by Correlation Identification combination m-sequence.
(3) present invention can show the system function of examining system in final display module.
Description of the drawings
Fig. 1 is a kind of composition frame chart of the Correlation Identification detection system based on FFT of the present invention；
Fig. 2 is the exemplary method flow diagram of specific implementation of the present invention.
Specific implementation mode
As shown in Figure 1, a kind of Correlation Identification detection system based on FFT of the present invention includes：First m-sequence generation module 1,
First the earth system 2, first time synchronization module 3, the first signal receiving module 4, the first the 5, first frequency of FFT- auto-correlations module
Domain conversion module 6, the first FFT- cross-correlation modules 7, the second frequency domain conversion module 8 and the first system function solve module 9, wherein
First FFT- auto-correlations module 5, the first frequency domain conversion module 6, the first FFT- cross-correlation modules 7,8 and of the second frequency domain conversion module
The first system function solves module 9 and belongs to data processing section；3 and first m-sequence of the first time synchronization module generates mould
Block 1 is connected with the first signal receiving module 2, and the first signal generator module 3 is connected with the first FFT- auto-correlations module 5 again, and first
FFT- auto-correlations module 5 is connected with the first frequency domain conversion module 6；First signal receiving module 4 is connected with the first cross-correlation module
7, the first FFT- cross-correlation modules 7 are connected with the second frequency domain conversion module 8；First frequency domain conversion module 6 and the conversion of the second frequency domain
Module 8 solves module 9 with the first system function and is connected；
The function of the first m-sequence generation module 1 is to be sent to the first the earth system after generating m-sequence according to setup parameter
System 2；Predetermined instant, which is arranged, in first time synchronization module 3 first makes the first m-sequence generation unit 13 generate phase at the time of setting
The m-sequence of parameter is answered, after generating m-sequence, m-sequence is just sent to the first the earth system 2 by electrode.In this way, just completing m sequences
The generation and transmission of row；The voltage data packet x of transmission can be collected simultaneously；
The function of the signal receiving module 4 is to complete the reception through current signal too greatly.According to time synchronization mould
The predetermined instant that block 3 is arranged, signal receiving module 4 are received through voltage signal too greatly, to form data packet y；
The first FFT- auto-correlations module 5 includes the first FFT frequency domain converting units 51, the first multiplication unit 52, first
IFFT converter units 53 and the first data cutout unit 54.Wherein the first FFT frequency domain converting units 51 and the first multiplication unit 52
It is connected, the first multiplication unit 52 is connected with the first IFFT converter units 53, the first IFFT units 53 and the first data cutout unit
54 are connected.After first m-sequence generation module 1 is connected with the first FFT frequency domain converting units 51, by the first FFT- auto-correlation modules
5, autocorrelative frequency-domain result has been calculated with the method for frequency domain, underlapped part therein has been intercepted out, obtains
Accurate time domain autocorrelation result；
The first frequency domain conversion module 6 includes FFT operators unit the 1st and frequency domain conversion results the 1st.First
FFT operators unit 61 is connected with the first data cutout unit 54 and the first auto-correlation frequency domain unit 62 respectively.First data are cut
The first auto-correlation frequency domain unit 62 is obtained after taking unit 54 to send the first FFT operators unit 61 in this way, just having obtained auto-correlation
Frequency-domain result.
The first FFT- cross-correlation modules 7 include the 2nd FFT frequency domain converting units 71, the second multiplication unit 72, second
IFFT converter units 73 and the second data cutout unit 74.2nd FFT frequency domain converting units 71 are connected with the second multiplication unit 72,
Second multiplication unit 72 and the 2nd IFFT converter units 73, the 2nd IFFT converter units 73 and the second data cutout unit 74.Together
When the first frequency-domain transform unit 51 be connected with the second multiplication unit 72, the first signal receiving module 4 and the 2nd FFT frequency-domain transform lists
Member 71 is connected.The received data packet y that first signal receiving module 4 obtains obtains frequency domain by the 2nd FFT frequency domain converting units 71
As a result, be then transferred to the second multiplication unit 72, while result is also transferred to square law list by the first frequency-domain transform unit 51
Member 72 is exactly the frequency-domain result of cross-correlation after the two multiplication, then obtains the transformation of time domain by the 2nd IFFT converter units 73
As a result, the last underlapped accurate cross correlation results for intercepting out time domain with the second data cutout unit 74.
The second frequency domain conversion module 8 includes the 2nd FFT operators unit 81 and the first cross-correlation frequency domain unit 82.Second
FFT operators unit 81 is connected with the 2nd FFT frequency domain converting units 71 and the first cross-correlation frequency domain unit 82 respectively.2nd FFT frequencies
Domain converting unit 71 obtains the first cross-correlation frequency domain unit 82 by the 2nd FFT operators unit 81, this has just been calculated mutually
The frequency-domain result of pass.
It includes the first divider 91 and the first the earth frequency domain identification unit 92 that the first system function, which solves module 9,.First
Divider 91 respectively with the first cross-correlation frequency domain unit 82, the first auto-correlation frequency domain unit 62 and the first the earth frequency domain identification list
Member 92 is connected.First cross-correlation frequency domain unit 82 and the first auto-correlation frequency domain unit 62 pass through the first divider 91, and first mutually
After related frequency domain unit 82 divided by the first auto-correlation frequency domain unit 62, the first the earth frequency domain identification unit 92 is obtained, as a result finally
It is stored in the first the earth frequency domain identification unit 92.
The relevant treatment work of the Correlation Identification of detection system is completed in frequency domain by FFT operators, operation efficiency pole
Height is particularly suitable for using when the data volume of correlation computations is larger.The first FFT- auto-correlations module 5, the first FFT- are mutual
Correlation module 7 solve detection system system function when the specific implementation process is as follows：First, transmitting terminal collects data
After packet x while receiving terminal all collect data packet y, data is given respectively the first FFT- auto-correlations module 5 and the first FFT- is mutual
Correlation module 7, after the first FFT frequency-domain transform units 51 in the first FFT- auto-correlations module 5 obtain the result X (w) of frequency domain,
The data packet x of transmitting terminal is consigned to the first multiplication unit 52, at this moment, the first multiplication unit 52 can carry out one of which data
Then conjugation processing calculates the autocorrelative frequency-domain result of cycle by correlation theorem.Calculation formula is： To recycle autocorrelative result.By correlation theorem it is recognised that cycle auto-correlation
The condition equal with cross-correlation beWherein M refers to that M dot cycles are related, LxRefer to
It is the points of x；Similarly, it recycles cross-correlation and the equal condition of cross-correlation isWherein
M is M dot cycle cross-correlation, LyFor the points of data packet y.According to this condition, intercepted out in data cutout unit non-overlapping
Correlated results come.Auto-correlation and cross-correlation thus are calculated with FFT, then according to the frequency domain identification formula of Correlation Identification
Calculate the system function of the earth.
The default setting of first time synchronization module 3 in the case of no modification is：12:01, i.e. predetermined instant 12:01 hair
Sending end and receiving terminal work at the same time.
According to another embodiment of the invention, the first time synchronization module 3 is for being arranged predetermined instant so that the
One m-sequence generation module 1 generates m-sequence at the time of setting；And make the first signal receiving module 4 is received in predetermined instant to pass through
Cross the electric signal of the first the earth system 2；
The first m-sequence generation module 1 is used to generate m-sequence in predetermined instant and is sent to the first the earth via electrode
System 2；And acquire the voltage data packet x of transmission；
The predetermined instant that first signal receiving module 4 is used to be arranged according to first time synchronization module 3 receives and passes through
The voltage signal of first the earth system 2, to obtain voltage data packet y；
The first FFT- auto-correlations module 5 includes the first FFT frequency-domain transform units 51, the first multiplication unit 52, first
IFFT converter units 53 and the first data cutout unit 54；
Wherein the first FFT frequency domain converting units 51 are connected with the first multiplication unit 52, the first multiplication unit 52 and first
IFFT converter units 53 are connected, and the first IFFT units 53 are connected with the first data cutout unit 54；
The first FFT frequency domain converting units 51 carry out FFT transform for voltage data packet x, and time-domain signal is converted to
Frequency-region signal；
First multiplication unit 52, for the frequency-region signal of voltage data packet x to be carried out conjugate operation；And conjugation is transported
Result after calculation is multiplied with the frequency-region signal of voltage data packet x, and calculation formula is
Wherein X (w) is the frequency domain of input signal, X*(w) it is the conjugation of the frequency-domain result of input signal；
The first IFFT converter units 53, the signal for exporting the first multiplication unit 52 carry out IFFT transformation, obtain
To time-domain signal, transformation for mula is
To recycle autocorrelative result；
The first data cutout unit 54, for intercepting the first IFFT converter units 53 cycle auto-correlation and linear
The equal part of auto-correlation obtains the autocorrelative result of time domain linear；
The first frequency domain conversion module 6 includes the first FFT operators unit 61 and the first auto-correlation frequency domain unit 62；Wherein
The first FFT operators unit 61 is connected with the first data cutout unit 54 and the first auto-correlation frequency domain unit 62 respectively；
The time domain autocorrelation result that the first FFT operators unit 61 is used to export the first data cutout unit 54 carries out
FFT operations, to obtain frequency domain autocorrelation result；And frequency domain autocorrelation result is sent to the first auto-correlation frequency domain unit and is protected
It deposits；
The first FFT- cross-correlation modules 7 include the 2nd FFT frequency-domain transform units 71, the second multiplication unit 72, second
IFFT converter units 73 and the second data cutout unit 74；
The 2nd FFT frequency-domain transform units 71 are connected with the second multiplication unit 72, second multiplication unit 72 and
Two IFFT converter units 73 are connected；The 2nd IFFT converter units 73 are connected with the second data cutout unit 74；
First frequency-domain transform unit 51 is connected with the second multiplication unit 72, the first signal receiving module 4 and the 2nd FFT
Frequency-domain transform unit 71 is connected；
The 2nd FFT frequency domain converting units 71 are used to carry out the voltage data packet y that the first signal receiving module 4 obtains
Time-domain signal is converted to frequency-region signal by FFT transform；
Second multiplication unit 72, the frequency domain letter of the voltage data packet x for exporting the first frequency-domain transform unit 51
It number is multiplied with the frequency-region signal of voltage data packet y of the 2nd FFT frequency domain converting units 71 output, obtains the frequency domain knot of cross-correlation
Fruit；
The 2nd IFFT converter units 73, the signal for exporting the second multiplication unit 72 carry out IFFT transformation, obtain
To time-domain signal；
The second data cutout unit 74, for intercepting in the time-domain signal that the 2nd IFFT converter units 73 export
Underlapped part obtains time domain cross correlation results；
The second frequency domain conversion module 8 includes the 2nd FFT operators unit 81 and the first cross-correlation frequency domain unit 82；
The 2nd FFT operators unit 81 respectively with 82 phase of the second data cutout unit 74 and the first cross-correlation frequency domain unit
Even；
The time domain cross correlation results that the 2nd FFT operators unit 81 is used to export the second data cutout unit 74 carry out
FFT operations obtain frequency domain cross correlation results；And frequency domain cross correlation results are sent to the first cross-correlation frequency domain unit 82 and are preserved；
It includes the first divider 91 and the first the earth frequency domain identification unit 92 that the first system function, which solves module 9,；It is described
First divider 91 respectively with the first cross-correlation frequency domain unit 82, the first auto-correlation frequency domain unit 62 and the first the earth frequency domain
Identification unit 92 is connected；
First divider 91 be used for the frequency domain cross correlation results that export the first cross-correlation frequency domain unit 82 divided by
The frequency domain cross correlation results of first auto-correlation frequency domain unit 62 output obtain the first the earth frequency domain identification as a result, and being stored in the
In one the earth frequency domain identification unit 92.
Below by embodiment to a kind of detection method of the Correlation Identification detection system based on FFT of the present invention further
It is described in detail；
As shown in Figure 2：
Step 201：First, the battery of instrument is connected, big voltage electrical generation equipment and sender electrode outside transmitter are connected,
After inspection connects, the power on button of transmitter on (X-on)；
Step 202：After connecting modules booting, transmitter will do initial work to peripheral apparatus.Peripheral apparatus
Including temperature sensor, AD1274, SD card, buzzer, PCF8563,24C02 and TFLCD display screen etc.；Initial work includes
The pretreatment operated is done to these peripheral hardwares, sets the communication protocol of peripheral hardware and equipment, checks whether peripheral hardware is normal etc.；
Step J01：When transmitter initializes, the actual conditions according to peripheral hardware are needed to judge whether modules initialize
Success, if so, thening follow the steps 204；If not, thening follow the steps 203；
Step 203：Peripheral hardware failure is initialized, then shows the corresponding error message of the peripheral hardware；
Step 204：After initializing peripheral hardware success, then show that the main interface of transmitter, main interface are broadly divided into three pieces, setting
Parameter (sample rate, chip width, m-sequence period, m-sequence exponent number, the voltage value of transmission, delivery time, transmission time interval),
Transmitter icon shows, operates brief description, and step 205 is then executed after instrument stabilizer；
Step 205：The parameter of m-sequence sent, sample rate, delivery time, transmission time interval (how long send) are set,
After setting completed according to prompt message, step 206 is then executed；
Step 206：GPS module obtains the time of GPS in real time, the present moment of not short more new instrument；
Step J03：Judged whether to reach delivery time according to the time of real-time update.If so, 216 are thened follow the steps,
If not, thening follow the steps 218；
Step 216：M-sequence is sent according to the parameter of the m-sequence of setting；
Step 218：Execution etc. is to be operated, and the stand-by period reaches delivery time；
Step 217：According to the big voltage signal of m-sequence that the sample rate acquisition of setting is sent, while as required acquisition
Data be transferred to receiver；
Step J05：Judge whether the time is more than or equal to time interval set by user, if not, thening follow the steps 216；
If so, thening follow the steps 220；
Step 208：First, the battery and receiving electrode of instrument are connected, after inspection connects, the booting of transmitter on (X-on) is pressed
Button；
Step 209：After connecting modules booting, receiver will do initial work to peripheral apparatus.Peripheral apparatus
Including temperature sensor, AD1274, SD card, buzzer, PCF8563,24C02 and TFLCD display screen etc.；Initial work includes
The pretreatment operated is done to these peripheral hardwares, sets the communication protocol of peripheral hardware and equipment, checks whether peripheral hardware is normal etc.；
Step J02：Judge whether modules initialization succeeds, if so, thening follow the steps 211；If it is not, then executing
Step 210；
Step 210：The initialization failure of some peripheral hardware, then show the corresponding error message of peripheral hardware；
Step 211：After initializing peripheral hardware success, then show that the main interface of receiver, main interface are broadly divided into three pieces, setting
Parameter (sample rate, the time of reception, receiving time interval), receiver icon show, operates brief description, after instrument stabilizer then
Execute step 212；
Step 212：The sample rate of receiver, the time of reception and receiving time interval are set.Note that these numerical value should be with
The numerical value of transmitter is identical；
Step 213：GPS module obtains the time of GPS in real time, constantly updates the present moment of receiver；
Step J04：Judge whether to reach the time of reception, if so, thening follow the steps 215；If not, thening follow the steps
214；
Step 214：According to the present moment that GPS is obtained, waits for and reach the time of reception；
Step 215：According to the sample rate collection voltages signal of user setting；
Step J06：Judge whether acquisition time is more than or equal to receiving time interval.If not, thening follow the steps 215；Such as
Fruit is to then follow the steps 219；
Step 219：The voltage signal for receiving voltage signal and receiver acquisition that transmitter transmits, is calculated big
The system function of ground system；
It is described calculate the earth system function mode be specially：
According to the relationship of input and output signal：
(1) u (t)=m (t) * h (t)+n (t)
Wherein, m (t) is the pseudorandom m signals for inputting the earth, has rhythmic pseudo-randomness；H (t) is to need to calculate
The ssystem transfer function come, n (t) are random noise signal, and y (t) is the collected signal of instrument；
(2) after being then FFT to equation both ends, multiplied by with Y*(w) it obtains
Wherein,For the frequency-domain result of the cyclic convolution of the cross-correlation function of input signal and output signal, X (w)
For the frequency domain of input signal, X*(w) it is the conjugation of the frequency-domain result of input signal, N (w) is the frequency-domain result of noise, due to making an uproar
Sound and input signal are uncorrelated, last N (w) X*(w) product is 0；
(3) after, Fast Fourier Transform (FFT) IFFT is carried out to equation both sides, obtains the Wiener Hopf equation of time domain
WhereinFor the cycle cross-correlation of input signal and output signal as a result,It is recycled from phase for input signal
The result of pass；
Then it is linear correlation, conversion stripes cycle dependent conversion according to the related and linearly related equal condition of cycle
Part isWherein M is that M dot cycles are related, LyFor the points of data packet y,Wherein M is that M dot cycles are related, LxFor the points of data packet y, then to equation two
FFT transform is done at end, is then calculated
System function and the display of examining system can be calculated by above formula.
Step 221：Terminate identification work；
When according to one embodiment of present invention, using detection system of the present invention, examining system selected first.Then it is arranged
The time of transmission and the parameter for sending signal.Examining system is a second-order system, and system function is：Sys_fun=
11.5e-5tSin (8.6t), it is 1V that the parameter of m-sequence, which is amplitude, and sample rate 1000Hz, exponent number is 12~14, and the period is 3, nothing
Interpolation, noise are 30V random noises.The hardware realization of core calculations FFT is mainly 32 single essences by STM32F407ZGT6
Spend hardware FPU units (FPU Float Point Unit), traditional fixed-point CPU calculate floating point unit when, it is relatively time consuming, and for by
For the CPU of FPU, then only need several instructions, budget speed quite fast.Used FPU is relative to traditional calculating sides CPU
Formula is higher by the decades of times even operational performance of hundreds of times.Meanwhile in written in code control, traditional CPU is compared also simply very
It is more.
The prior art well known to the technical staff of the detailed professional domain is not done in the specification of the present invention.
Schemes described above is a kind of application scenario of the system, which is mainly the correlation being used in big data quantity
The invention is used in other field or does corresponding modification and several improvement should also be as being considered as the invention by the System Discrimination of identification
Protection domain.
Claims (4)
1. a kind of Correlation Identification detection system based on FFT, it is characterised in that including：First m-sequence generation module (1), first
The earth system (2), first time synchronization module (3), the first signal receiving module (4), the first FFT- auto-correlations module (5),
One frequency domain conversion module (6), the first FFT- cross-correlation modules (7), the second frequency domain conversion module (8) and the first system function solve
Module (9)；The first time synchronization module (3) and the first m-sequence generation module (1) and the first signal receiving module (4) phase
Even, the first m-sequence generation module (1) is connected with the first FFT- auto-correlations module (5), the first FFT- auto-correlations module (5) and the
One frequency domain conversion module (6) is connected；First signal receiving module (4) is connected (7) with the first FFT- cross-correlation modules, the first FFT-
Cross-correlation module (7) is connected with the second frequency domain conversion module (8)；First frequency domain conversion module (6) and the second frequency domain conversion module
(8) module (9) is solved with the first system function respectively to be connected；
The first time synchronization module (3) is for being arranged predetermined instant so that the first m-sequence generation module (1) is in setting
Moment generates m-sequence；And the first signal receiving module (4) is made to receive the electricity Jing Guo the first the earth system (2) in predetermined instant
Signal；
The first m-sequence generation module (1) is used to generate m-sequence in predetermined instant and be sent to the first the earth via electrode
It unites (2)；And acquire the voltage data packet x of transmission；
The predetermined instant that first signal receiving module (4) is used to be arranged according to first time synchronization module (3) receives and passes through
The voltage signal of first the earth system (2), to obtain voltage data packet y；
The first FFT- auto-correlations module (5) includes the first FFT frequency-domain transform units (51), the first multiplication unit (52), the
One IFFT converter units (53) and the first data cutout unit (54)；
Wherein the first FFT frequency-domain transform units (51) are connected with the first multiplication unit (52), the first multiplication unit (52) and first
IFFT converter units (53) are connected, and the first IFFT converter units (53) are connected with the first data cutout unit (54)；
The first FFT frequency-domain transform units (51) carry out FFT transform for voltage data packet x, and time-domain signal is converted to frequency
Domain signal；
First multiplication unit (52), for the frequency-region signal of voltage data packet x to be carried out conjugate operation；And by conjugate operation
Result afterwards is multiplied with the frequency-region signal of voltage data packet x, and calculation formula is
Wherein X (w) is the frequency domain of input signal, X*(w) it is the conjugation of the frequency-domain result of input signal；
The first IFFT converter units (53), the signal for exporting the first multiplication unit (52) carry out IFFT transformation, obtain
To time-domain signal, transformation for mula is：
To recycle autocorrelative result；
The first data cutout unit (54), for intercepting the first IFFT converter units (53) cycle auto-correlation and linear
The equal part of auto-correlation obtains the autocorrelative result of time domain linear；
The first frequency domain conversion module (6) includes the first FFT operators unit (61) and the first auto-correlation frequency domain unit (62)；Its
Described in the first FFT operators unit (61) respectively with the first data cutout unit (54) and the first auto-correlation frequency domain unit (62)
It is connected；
The time domain autocorrelation result that the first FFT operators unit (61) is used to export the first data cutout unit (54) carries out
FFT operations, to obtain frequency domain autocorrelation result；And frequency domain autocorrelation result is sent to the first auto-correlation frequency domain unit (62)
It preserves；
The first FFT- cross-correlation modules (7) include the 2nd FFT frequency-domain transform units (71), the second multiplication unit (72), the
Two IFFT converter units (73) and the second data cutout unit (74)；
The 2nd FFT frequency-domain transform units (71) are connected with the second multiplication unit (72), second multiplication unit (72) with
2nd IFFT converter units (73) are connected；The 2nd IFFT converter units (73) are connected with the second data cutout unit (74)；
The first FFT frequency-domain transform units (51) are connected with the second multiplication unit (72), the first signal receiving module (4) and the
Two FFT frequency-domain transform units (71) are connected；
The 2nd FFT frequency domain converting units (71) are used to carry out the voltage data packet y that the first signal receiving module (4) obtains
Time-domain signal is converted to frequency-region signal by FFT transform；
Second multiplication unit (72), the frequency domain of the voltage data packet x for exporting the first FFT frequency-domain transform units (51)
Signal is multiplied with the frequency-region signal for the voltage data packet y that the 2nd FFT frequency domain converting units (71) export, and obtains the frequency domain of cross-correlation
As a result；
The 2nd IFFT converter units (73), the signal for exporting the second multiplication unit (72) carry out IFFT transformation, obtain
To time-domain signal；
The second data cutout unit (74), for intercepting in the time-domain signal that the 2nd IFFT converter units (73) export
Underlapped part obtains time domain cross correlation results；
The second frequency domain conversion module (8) includes the 2nd FFT operators unit (81) and the first cross-correlation frequency domain unit (82)；
The 2nd FFT operators unit (81) respectively with the second data cutout unit (74) and the first cross-correlation frequency domain unit (82)
It is connected；
The time domain cross correlation results that the 2nd FFT operators unit (81) is used to export the second data cutout unit (74) carry out
FFT operations obtain frequency domain cross correlation results；And frequency domain cross correlation results are sent to the first cross-correlation frequency domain unit (82) and are protected
It deposits；
It includes the first divider (91) and the first the earth frequency domain identification unit (92) that the first system function, which solves module (9),；Institute
State the first divider (91) respectively with the first cross-correlation frequency domain unit (82), the first auto-correlation frequency domain unit (62) and first
The earth frequency domain identification unit (92) is connected；
First divider (91) be used for the frequency domain cross correlation results that export the first cross-correlation frequency domain unit (82) divided by
The frequency domain cross correlation results of first auto-correlation frequency domain unit (62) output, obtain the first the earth frequency domain identification as a result, and being stored in
In first the earth frequency domain identification unit (92).
2. a kind of Correlation Identification detection system based on FFT according to claim 1, it is characterised in that：
First time synchronization module (3) generates m-sequence according to parameter preset.
3. a kind of Correlation Identification detection system based on FFT according to claim 2, it is characterised in that：
The default setting of first time synchronization module (3) in the case of no modification is：12:01, i.e., in predetermined instant 12:01 hair
Sending end and receiving terminal work at the same time.
4. according to the detection method of any Correlation Identification detection systems based on FFT of claim 1-3, feature exists
In：
Step 201：It powers on；
Step 202：Peripheral apparatus initial work；The peripheral apparatus include temperature sensor, AD1274, SD card, buzzer,
PCF8563,24C02 and TFLCD display screen；The initial work includes doing the pretreatment operated to these peripheral hardwares, is set
The communication protocol of peripheral hardware and equipment checks whether peripheral hardware is normal etc.；
Step J01：Judge whether modules initialize success, if so, thening follow the steps 204；If not, thening follow the steps
203；
Step 203：Peripheral apparatus failure is initialized, then shows the corresponding error message of the peripheral apparatus；
Step 204：After initializing peripheral hardware success, then main interface is shown, then execute step 205；
Step 205：The parameter of m-sequence sent, sample rate, delivery time, transmission time interval are set, step 206 is executed；
Step 206:GPS time is obtained, time synchronization is carried out；
Step J03：Judged whether to reach delivery time according to the time obtained in real time；If so, 216 are thened follow the steps, if
It is no, then follow the steps 218；
Step 216：M-sequence is sent according to the parameter of the m-sequence of setting, executes step 217；
Step 218：Execution etc. is to be operated, and the stand-by period reaches delivery time；
Step 217：According to the signal for the m-sequence that the sample rate acquisition of setting is sent, while the data of acquisition are passed as required
It is defeated by receiving terminal；
Step J05：Judge whether the time is more than or equal to time interval set by user, if not, thening follow the steps 216；If
It is to then follow the steps 220；
Step 220：Stop sending m-sequence；
Step 208：It powers on；
Step 209：Peripheral apparatus initial work；The peripheral apparatus include temperature sensor, AD1274, SD card, buzzer,
PCF8563,24C02 and TFLCD display screen；Initial work includes doing the pretreatment operated to these peripheral hardwares, sets peripheral hardware
With the communication protocol of equipment, check whether peripheral hardware is normal etc.；
Step J02：Judge whether modules initialization succeeds, if so, thening follow the steps 211；If not, thening follow the steps
210；
Step 210：The initialization failure of some peripheral hardware, then show the corresponding error message of peripheral hardware；
Step 211：After initializing peripheral hardware success, then shows the main interface of receiver, step 212 is then executed after instrument stabilizer；
Step 212：Setting sample rate identical with the first signal generator module, the time of reception and receiving time interval；
Step 213：GPS module obtains the time of GPS in real time；
Step J04：Judged whether to reach the time of reception according to the time obtained in real time, if so, thening follow the steps 215；If
It is no, then follow the steps 214；
Step 214：According to the present moment that GPS is obtained, waits for and reach the time of reception；
Step 215：According to the sample rate collection voltages signal of setting；
Step J06：Judge whether acquisition time is more than or equal to receiving time interval；If not, thening follow the steps 215；If so,
Then follow the steps 219；
Step 219：According to the electric signal for sending voltage data packet x and the receiving voltage data packet y of acquisition of acquisition, it is calculated
The system function of the earth system and display；
The mode of system function for calculating the earth system is specially：
According to the relationship of input and output signal：
(1) y (t)=x (t) * h (t)+n (t)
Wherein, x (t) is the pseudorandom m signals for inputting the earth, has rhythmic pseudo-randomness；H (t) needs calculate
Ssystem transfer function, n (t) are random noise signal, and y (t) is the collected signal of instrument；
(2) after being then FFT to equation both ends, multiplied by with Y*(w) it obtains
Wherein,For the frequency-domain result of the cyclic convolution of the cross-correlation function of input signal and output signal, X (w) is input
The frequency domain of signal, X*(w) be input signal frequency-domain result conjugation, N (w) is the frequency-domain result of noise, due to noise and defeated
It is uncorrelated to enter signal, last N (w) X*(w) product is 0；
(3) after, Fast Fourier Transform (FFT) IFFT is carried out to equation both sides, obtains the Wiener Hopf equation of time domain
WhereinFor the result of input signal and the cycle cross-correlation of output signal；
Then it is linear correlation cycle dependent conversion, switch condition is according to the related and linearly related equal condition of cycle
Wherein M is that M dot cycles are related, LyFor the points of data packet y；
Wherein M is that M dot cycles are related, LxFor the points of data packet y；
Then FFT transform is done to equation both ends, is then calculated
(4)
The system function of examining system can be calculated by above formula；
Step 221：Terminate identification.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610556801.5A CN106155983B (en) | 2016-07-14 | 2016-07-14 | A kind of Correlation Identification detection system and its method based on FFT |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610556801.5A CN106155983B (en) | 2016-07-14 | 2016-07-14 | A kind of Correlation Identification detection system and its method based on FFT |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106155983A CN106155983A (en) | 2016-11-23 |
| CN106155983B true CN106155983B (en) | 2018-09-14 |
Family
ID=58059141
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610556801.5A Expired - Fee Related CN106155983B (en) | 2016-07-14 | 2016-07-14 | A kind of Correlation Identification detection system and its method based on FFT |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106155983B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106649200B (en) * | 2016-12-27 | 2019-06-04 | 南京大学 | One kind being based on time-multiplexed auto-correlation computation VLSI design method |
| CN107219300B (en) * | 2017-05-23 | 2019-09-03 | 徐工集团工程机械股份有限公司 | A kind of compactness detection on locomotive system based on the degree of correlation |
| CN109444961A (en) * | 2018-09-25 | 2019-03-08 | 中国地质大学（北京） | A kind of multithreading Correlation Identification electrical prospecting reception system and method based on all phase |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202330736U (en) * | 2011-11-16 | 2012-07-11 | 中国地质大学（北京） | Electric-method exploration device |
| CN104391330A (en) * | 2014-12-01 | 2015-03-04 | 中国地质大学（北京） | Relevant identification technology based electrical prospecting system |
| CN104391329A (en) * | 2014-12-01 | 2015-03-04 | 中国地质大学（北京） | Receiving device for electric exploration system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008081469A2 (en) * | 2007-01-03 | 2008-07-10 | Council Of Scientific & Industrial Research | A process and device for measurement of spectral induced polarization response using pseudo random binary sequence (prbs) current source |
-
2016
- 2016-07-14 CN CN201610556801.5A patent/CN106155983B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202330736U (en) * | 2011-11-16 | 2012-07-11 | 中国地质大学（北京） | Electric-method exploration device |
| CN104391330A (en) * | 2014-12-01 | 2015-03-04 | 中国地质大学（北京） | Relevant identification technology based electrical prospecting system |
| CN104391329A (en) * | 2014-12-01 | 2015-03-04 | 中国地质大学（北京） | Receiving device for electric exploration system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106155983A (en) | 2016-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106155983B (en) | A kind of Correlation Identification detection system and its method based on FFT | |
| CN111083080B (en) | Method and device for synchronizing PRACH (physical random Access channel) of broadband satellite communication system | |
| CN104391330B (en) | Relevant identification technology based electrical prospecting system | |
| CN110988602B (en) | S-transformation-based traveling wave protection method for hybrid direct current transmission line | |
| CN106534002A (en) | Power line channel estimation method based on compressed sensing | |
| CN104635109A (en) | High-accuracy method for directly calculating and extracting amplitude value and time constant of transient aperiodic component of power grid from fault recording data | |
| CN104796250A (en) | Side channel attack method for implementation of RSA (Rivest, Shamir and Adleman) cipher algorithms M-ary | |
| CN104979829A (en) | Harmonic current estimation method based on improved FastICA algorithm | |
| CN104483539A (en) | Active power rapid measuring method based on Taylor expansion | |
| CN102944802B (en) | Estimation method for dynamic frequency response function of voltage transformer and realizing device thereof | |
| CN104811297A (en) | Method for modular multiplication remainder input side channel attacks aiming at M-ary implementation of RSA | |
| CN106125142B (en) | A kind of Underground medium system and its method based on Correlation Identification frequency domain method | |
| CN104535989B (en) | The acoustical signal decomposition method of way more than shallow sea | |
| CN104393947B (en) | The implementation method and device of frequency domain ZC sequence | |
| CN107703358B (en) | A kind of phase locked algorithm based on improvement Second Order Generalized Integrator | |
| CN106772193B (en) | Measuring method using current transformer frequency characteristic measuring device | |
| CN107770104A (en) | A kind of channel estimation pilot optimization method and device based on compressed sensing | |
| CN106021183A (en) | Mixed base FFT-combined relevant identification processing system and method | |
| CN104506271A (en) | Method and device for generating ZC (Zadoff-Chu) sequence | |
| CN112615648B (en) | Accurate narrowband noise resistance method used in power line carrier communication Rx process | |
| CN114035766A (en) | Disposable random number generating device and application method | |
| CN108362940A (en) | A kind of dynamic harmonic frequency extraction method of half cycles distorted signal | |
| CN103105532A (en) | Harmonic electric energy measuring system of improved Fourier transform and control method thereof | |
| CN203012015U (en) | Harmonic energy metering system improving Fourier transformation | |
| CN106019395A (en) | Electrical method receiving and processing system and method thereof based on correlation identification |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180914 Termination date: 20200714 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |