CN106054269A - Frequency partition correlation identification system and method thereof - Google Patents
Frequency partition correlation identification system and method thereof Download PDFInfo
- Publication number
- CN106054269A CN106054269A CN201610556834.XA CN201610556834A CN106054269A CN 106054269 A CN106054269 A CN 106054269A CN 201610556834 A CN201610556834 A CN 201610556834A CN 106054269 A CN106054269 A CN 106054269A
- Authority
- CN
- China
- Prior art keywords
- unit
- correlation
- frequency
- signal
- module
- 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
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Abstract
The invention relates to a frequency partition correlation identification system comprising a first frequency resolution module, a first signal transmitting system, a first ground system, a first signal receiving module, a first correlation identification module, a second correlation identification module, a first result processing module and a first amplitude spectrum identification module. Amplitude spectrum identification of different detection depths is performed by changing transmitting signal parameters so that the more accurate ground system function identification result is obtained, and thus the system has the performance of high efficiency, high accuracy and simple algorithm.
Description
Technical field
The present invention relates to prospecting instrument field, particularly relate to the Correlation Identification system and method for a kind of frequency segmentation.
Background technology
Baseghi, Behshad etc. in 1992 propose and realize wideband digital channel receiver with FFT, give simultaneously
The structure that system realizes.The principle of frequency domain channel is to utilize FFT to carry out even partition frequency domain so that each channel has identical
Bandwidth and frequency response.In most cases, the splicing of signal is applied in image procossing, mainly utilizes Fourier transformation by signal
Transform to frequency domain to be analyzed.The relative displacement utilizing time domain is equivalent to the phase place of frequency domain crosspower spectrum, can calculate signal
Relative displacement, then utilize weighted mean method method to eliminate splicing gap, simple and preferable result can be reached.
Frequency domain splicing method refers to determining that signal is after the channel occupied by channel signal and signal, effective from these
Signal spectrum is extracted by channel, and is spliced into complete signal spectrum.Frequency domain splicing method advantage be to across
The process of channel signal is all carried out after channelization structure, and handling process is clear, it is easy to hardware realize, and its can process across
The signal of the most multiple channels.Therefore, utilizing frequency domain splicing method that the earth system function carries out identification effectively suppresses noise to do
It is critically important for disturbing and improving identification accuracy on amplitude spectrum.
Summary of the invention
For overcoming the deficiencies in the prior art, the present invention provides the Correlation Identification system and method for a kind of frequency segmentation,
On the basis of can be implemented in single the earth system identification, the earth system of different depth is spliced, effectively on frequency spectrum
Inhibit noise, improve the identification precision of amplitude spectrum.
A kind of Correlation Identification system of frequency segmentation, including: first frequency module resolution, the first signal sending system,
First the earth system, the first signal receiving module, the first Correlation Identification module, the second Correlation Identification module, the first result treatment
Module and the first amplitude spectrum identification module, wherein first frequency module resolution, the first Correlation Identification module, the second Correlation Identification
Module, the first result treatment module and amplitude spectrum identification module belong to data processing section；Described first frequency resolution mould
Block and the first signal transmitting module are connected, the first signal transmitting module and the first the earth system and the first Correlation Identification module and the
Two Correlation Identification modules are connected, and the first the earth system is connected with the first signal receiving module again；First signal receiving module and
One Correlation Identification module is connected with the second Correlation Identification module, the first Correlation Identification module and the second Correlation Identification module and first
Result treatment module is connected；First result treatment module is connected with the first amplitude spectrum identification module；
Described first frequency module resolution includes the first chip width operator unit, the second chip width connected side by side
Unit and the first signal length operator unit being sequentially connected, secondary signal length operator's unit.First chip width operator list
Unit is connected with the first signal length operator unit, and the first signal length operator unit and the first signal transmitting unit are connected；Second
Chip width operator unit is connected with secondary signal length operator's unit, and secondary signal length operator's unit sends with secondary signal
Unit is connected, and it is deep that the first first chip width operator unit and the first chip width operator unit pass through the detection to the earth system
Degree difference arranges different chip width, after determining chip width, owing to needing to send the consistent ability of frequency resolution of signal
Signal is carried out segmentation and splicing, then it is corresponding with the calculating of secondary signal length operator's unit to swim the first signal length operator unit
Send signal length, meanwhile, the result of generation is sent the first signal transmitting module.So, just complete sending signal
The setting of parameter；
Described first signal transmitting module includes the first signal transmitting unit and secondary signal transmitting element.First signal is sent out
Send unit to be connected with the first the earth system with secondary signal transmitting element simultaneously.Arrange not according to first frequency module resolution
Produce different m-sequence with sending signal parameter, by the first signal transmitting unit and secondary signal transmitting element, transmission is believed
Number it is sent to the first the earth system, thus completes the transmission of signal；
Described first signal receiving module includes that the first reception electrode unit and first gathers memory element and second and gathers
Memory element.First reception electrode and the first the earth system and first gather memory element and are connected with the second collection memory element,
First gathers memory element is connected with the first auto-correlation operator unit and the first cross-correlation operator unit, and second gathers memory element
It is connected with the first auto-correlation operator unit and the first cross-correlation operator unit.First receives electrode receives through the first the earth system
Voltage signal, after first gathers memory element and the second collection memory element collection storage, is transferred to the first phase respectively
Close identification module and the second Correlation Identification module, thus complete signal receiving module；
Described first Correlation Identification module include the first auto-correlation operator unit, the first cross-correlation operator unit, first from
Relevant frequency domain unit, the first cross-correlation frequency domain unit and the first system function solve unit.First auto-correlation operator unit and
One auto-correlation frequency domain unit is connected；First cross-correlation operator unit and the first cross-correlation frequency domain unit are connected, the first auto-correlation frequency
Territory unit and the first cross-correlation frequency domain unit all solve unit with the first system function and are connected.First auto-correlation operator unit passes through
The signal that first signal transmitting unit sends obtains autocorrelative frequency-domain result through the first auto-correlation frequency domain unit, and first is mutual
The signal and first that sends that pass operator unit is obtained by the first transmitting element gathers the reception signal process that memory element obtains
First cross-correlation frequency domain unit obtains the frequency-domain result of cross-correlation, then solves unit by the first system function and just obtained correspondence
The frequency domain identification result of the earth system of the degree of depth.
Described second Correlation Identification module include the second auto-correlation operator unit, the second cross-correlation operator unit, second from
Relevant frequency domain unit, the second cross-correlation frequency domain unit and second system function solve unit.Second auto-correlation operator unit and
Two auto-correlation frequency domain unit are connected；Second cross-correlation operator unit and the second cross-correlation frequency domain unit are connected, the second auto-correlation frequency
Territory unit and the second cross-correlation frequency domain unit all solve unit with second system function and are connected.Second auto-correlation operator unit passes through
The signal that secondary signal transmitting element sends obtains autocorrelative frequency-domain result through the second auto-correlation frequency domain unit, and second is mutual
The signal and second that sends that pass operator unit is obtained by the second transmitting element gathers the reception signal process that memory element obtains
Second cross-correlation frequency domain unit obtains the frequency-domain result of cross-correlation, then solves unit by second system function and just obtained correspondence
The frequency domain identification result of the earth system of the degree of depth.
Described first result treatment module include the first system frequency domain memory element, second system frequency domain memory element,
One cut-off frequency operator unit and the second cut-off frequency operator unit.The first system frequency domain memory element is calculated with the first cut-off frequency
Subelement is connected, and second system frequency domain memory element is connected with the second cut-off frequency operator unit.The first system function solves list
The earth system frequency-domain result of the corresponding degree of depth that unit obtains is stored in the first system frequency domain memory element, through the first cutoff frequency
Rate operator unit obtains this degree of depth correspondence cut-off frequency；Second system function solves the earth system of the corresponding degree of depth that unit obtains
Frequency-domain result is stored in second system frequency domain memory element, obtains this degree of depth correspondence through the second cut-off frequency operator unit and cuts
Only frequency, so, has just obtained the system frequency domain identification of different frequency corresponding to different investigation depth.
Described first amplitude spectrum identification module includes the first frequency spectrum concatenation unit and the first amplitude spectrum identification unit.First frequency
Spectrum concatenation unit and the first amplitude spectrum identification unit are connected.First frequency spectrum concatenation unit is by the first cut-off frequency operator unit and
Frequency domain identification result obtained by two cut-off frequency operator unit is spliced, and obtains the frequency domain identification of a new the earth system
As a result, the amplitude spectrum of big system frequency domain identification result is obtained through the first amplitude spectrum identification unit.
Present invention advantage compared with prior art is:
(1) under conditions of the Underground medium system of Correlation Identification frequency domain method, when frequency resolution is certain, chip
Width is different, during signal length difference, and the frequency domain splicing to the system function of different investigation depths, obtain the amplitude spectrum of identification more
For accurately, the identification effect of more preferable useful frequency range can be obtained.
(2) by the Correlation Identification system and method for frequency segmentation, it can be seen that to noise jamming under the conditions of different parameters
Inhibition.
(3) present invention is applicable to the parameter detecting of system function of the earth, solid, fluid media (medium), is mainly used in underground and is situated between
Matter detects.
Accompanying drawing explanation
Fig. 1 is the composition frame chart of the Correlation Identification system of a kind of frequency segmentation of the present invention；
Fig. 2 is the method flow diagram being embodied as example of the present invention.
Detailed description of the invention
As it is shown in figure 1, the Correlation Identification system of a kind of frequency segmentation of the present invention, including: first frequency module resolution 1,
First signal sending system the 2, first the earth system the 3, first signal receiving module the 4, first Correlation Identification module 5, second is correlated with
Identification module the 6, first result treatment module 7 and the first amplitude spectrum identification module 8, wherein first frequency module resolution 1, first
Correlation Identification module the 5, second Correlation Identification module the 6, first result treatment module 7 and amplitude spectrum identification module 8 belong to data
Process part；Described first frequency module resolution 1 is connected with the first signal transmitting module 2, the first signal transmitting module 2 and
One the earth system 3 is connected with the first Correlation Identification module 5 and the second Correlation Identification module 6, and the first the earth system 3 is believed with first again
Number receiver module 4 is connected；First signal receiving module 4 is connected with the first Correlation Identification module 5 and the second Correlation Identification module 6,
First Correlation Identification module 5 is connected with the second Correlation Identification module 6 and the first result treatment module 7；First result treatment module 7
It is connected with the first amplitude spectrum identification module 8；
Described first frequency module resolution 1 includes first chip width operator unit the 11, second chip connected side by side
Width unit (12) and the first signal length operator unit 13 being sequentially connected, secondary signal length operator's unit 14.First yard
Sheet width operator unit 11 is connected with the first signal length operator unit 13, the first signal length operator unit 13 and the first signal
Transmitting element 21 is connected；Second chip width operator unit 12 is connected with secondary signal length operator's unit 14, and secondary signal is long
Degree operator unit 14 is connected with secondary signal transmitting element 22, the first first chip width operator unit 11 and the first chip width
Operator unit 12 is by arranging different chip width to the investigation depth difference of the earth system, after determining chip width, due to
Need the frequency resolution sending signal that signal unanimously could carry out segmentation and splicing, then swim the first signal length operator unit
13 and secondary signal length operator's unit 14 calculate correspondence transmission signal length, meanwhile, will generate result send first letter
Number sending module 2.So, the setting to the parameter sending signal is just completed；
Described first signal transmitting module 2 includes the first signal transmitting unit 21 and secondary signal transmitting element 22.First
Signal transmitting unit 21 is connected with the first the earth system 3 with secondary signal transmitting element 22 simultaneously.According to first frequency resolution
The different signal parameters that send that module 1 is arranged produce different m-sequence, are sent out by the first signal transmitting unit 21 and secondary signal
Send unit 22 that transmission signal is sent to the first the earth system 3, thus complete the transmission of signal；
Described first signal receiving module 4 includes that the first reception electrode unit 41 and first gathers memory element 42 and second
Gather memory element.First receives electrode 41 and the first the earth system 3 and the first collection memory element 42 and the second collection storage
Unit 43 is connected, and first gathers memory element 42 and the first auto-correlation operator unit 51 and the first cross-correlation operator unit 52 phase
Even, the second collection memory element 43 is connected with the first auto-correlation operator unit 61 and the first cross-correlation operator unit 62.First connects
Receive electrode 41 and receive voltage signal through the first the earth system 3, single through the first collection memory element 42 and the second collection storage
After unit 43 gathers storage, it is transferred to the first Correlation Identification module 5 and the second Correlation Identification module 6 respectively, thus completes letter
Number receiver module；
Described first Correlation Identification module 5 include first auto-correlation operator unit the 51, first cross-correlation operator unit 52,
One auto-correlation frequency domain unit the 53, first cross-correlation frequency domain unit 54 and the first system function solve unit 55.First auto-correlation is calculated
Subelement 51 is connected with the first auto-correlation frequency domain unit 53；First cross-correlation operator unit 52 and the first cross-correlation frequency domain unit 54
Being connected, the first auto-correlation frequency domain unit 53 and the first cross-correlation frequency domain unit 54 all solve unit 55 phase with the first system function
Even.First auto-correlation operator unit 51 passes through the signal of the first signal transmitting unit 21 transmission through the first auto-correlation frequency domain unit
53 obtain autocorrelative frequency-domain result, transmission signal that the first cross-correlation operator unit 52 is obtained by the first transmitting element 21 and
The signal that receives that first collection memory element 42 obtains obtains the frequency-domain result of cross-correlation through the first cross-correlation frequency domain unit 54,
The frequency domain identification result that unit 55 has just obtained the earth system of the corresponding degree of depth is solved again by the first system function.
Described second Correlation Identification module 6 include second auto-correlation operator unit the 61, second cross-correlation operator unit 62,
Two auto-correlation frequency domain unit the 63, second cross-correlation frequency domain unit 64 and second system function solve unit 65.Second auto-correlation is calculated
Subelement 61 is connected with the second auto-correlation frequency domain unit 63；Second cross-correlation operator unit 62 and the second cross-correlation frequency domain unit 64
Being connected, the second auto-correlation frequency domain unit 63 and the second cross-correlation frequency domain unit 64 all solve unit 65 phase with second system function
Even.Second auto-correlation operator unit 61 passes through the signal of secondary signal transmitting element 22 transmission through the second auto-correlation frequency domain unit
63 obtain autocorrelative frequency-domain result, transmission signal that the second cross-correlation operator unit 62 is obtained by the second transmitting element 22 and
The signal that receives that second collection memory element 43 obtains obtains the frequency-domain result of cross-correlation through the second cross-correlation frequency domain unit 64,
The frequency domain identification result that unit 65 has just obtained the earth system of the corresponding degree of depth is solved again by second system function.
Described first result treatment module 7 includes the first system frequency domain memory element 71, second system frequency domain memory element
72, the first cut-off frequency operator unit 73 and the second cut-off frequency operator unit 74.The first system frequency domain memory element 71 and
One cut-off frequency operator unit 73 is connected, and second system frequency domain memory element 72 is connected with the second cut-off frequency operator unit 71.
The first system function solves the earth system frequency-domain result of the corresponding degree of depth that unit 55 obtains and is stored in the storage of the first system frequency domain
In unit 71, obtain this degree of depth correspondence cut-off frequency through the first cut-off frequency operator unit 73；Second system function solves list
The earth system frequency-domain result of the corresponding degree of depth that unit 65 obtains is stored in second system frequency domain memory element 72, through second section
Only frequency operator unit 74 obtains this degree of depth correspondence cut-off frequency, so, has just obtained the difference corresponding to different investigation depth
The system frequency domain identification of frequency.
Described first amplitude spectrum identification module 8 includes the first frequency spectrum concatenation unit 81 and the first amplitude spectrum identification unit 82.
First frequency spectrum concatenation unit 81 is connected with the first amplitude spectrum identification unit 82.First frequency spectrum concatenation unit 81 is by the first cut-off frequency
Frequency domain identification result obtained by operator unit 73 and the second cut-off frequency operator unit 74 is spliced, obtain one new big
The frequency domain identification result of ground system, obtains the amplitude spectrum of big system frequency domain identification result through the first amplitude spectrum identification unit 81.
The parameter of the m-sequence of the first signal transmitting unit 21 arranges and is defaulted as: interpolation number is 4, and exponent number is (4095,12 rank
Chip), amplitude is 1V.The parameter of the m-sequence of secondary signal transmitting element 22 arranges and is defaulted as: interpolation number is 2, and exponent number is 13 rank
(8191 chips), amplitude is 1V.
Another embodiment of the present invention, is utilized the present invention to detect the earth system different depth, is spelled by frequency domain
Connection obtains new the earth system function.First, according to the needs of investigation depth, arrange the frequency sending signal m-sequence 0～
40Hz, then represent chip width with interpolation number, if the interpolation number sending signal a is 4, the exponent number of m-sequence is 12 rank, so,
The signal cut-off frequency calculated before frequency spectrum splices is 0～20Hz；If the interpolation number sending signal b is 2, the exponent number of m-sequence is
13 rank, ensure and send signal a identical frequency resolution refused to use, the signal cut-off frequency calculated before frequency spectrum splices is 20
～40Hz.
Below the Correlation Identification method of the frequency segmentation of the present invention is further described in detail；
As shown in Figure 2:
Step 201: first, connects instrument with power supply (accumulator etc.), presses switch, starts instrument, on instrument
Electricity；
Step 202: after instrument energising, instrument will perform system initialization operation on backstage, such as, check the unit of self
Device situation is the most normal, it is judged that ambient parameter is the most normal, checks battery electric quantity etc., and waits that user setup sends letter
Number；
Step 203: according to investigation depth, sends signal frequency, and noise jamming difference calculates chip width；
Step 204: according to formula signal calculated length:
Computing formula isΔ F is frequency resolution, and Δ is chip width, and N is signal length；
Step 205: after completing step 203 and step 204, instrument enters step 205, sends the letter setting parameter
Number；
Step 206: signal is by the earth system:
Step 207: receive the signal by step 207 gained；
Step 208: calculate the auto-correlation function R of the signal received by step 207xx(t):
Formula is
Step 209: calculate the auto-correlation function R of the signal received by step 207xy(t):
Formula is
Step 210: obtain auto-correlation function R by Fourier transformationxx(t) certainly compose Rxx(ω) with cross-correlation function Rxy
The cross-spectrum R of (t)xy(ω)；
Step 211: obtain from composing R completing step 210xx(ω) with cross-spectrum Rxy(ω), after, instrument enters step 210, asks
Solve system function, the system function H (ω) of calculating the earth system:
Computing formula H (ω)=Rxy(ω)/Rxx(ω)；
Step 212: obtained the cutoff frequency of signal by different transmission signals (chip width is different, and signal length figure leads to)
Rate；
Step 213: connect the cut-off frequency of two signals, splices the frequency domain of two signals:
Step 214: the operation that data derive；Data in this step, can be derived with operating instrument, be spelled by user
Connect the amplitude spectrum identification result of rear signal, carry out data process and depth analysis for user's later stage；
Step 215: after completing the operation that data derive, the splicing task of signal completes；During next operation, Ke Yichong
Multiple above step；
The present invention calculates first frequency range and the data result of second frequency range respectively by STM32F407, is stored in
In SD card, in SD card, carry out frequency domain splicing, as a frequency range output.Owing to data volume is relatively big, main control chip low memory,
SD card internal memory is relatively big, and access speed is very fast, can meet the requirement of storage two segment datas.
The realization of the frequency domain segmentation of the present invention is mainly manifested in frequency resolution module and amplitude identification module, divides frequency
Resolution is certain, chip width is different, the signal that signal length is different, and the system of the different investigation depths corresponding to it carries out frequency
Territory is spliced, it can be seen that inhibition to noise jamming under the conditions of different parameters, obtains distinguishing of more preferable useful frequency range simultaneously
Know result.
The content not being described in detail in description of the invention belongs to prior art known to professional and technical personnel in the field.
The above is only the preferred embodiment of the present invention, it is noted that come for those of ordinary skill in the art
Saying, under the premise without departing from the principles of the invention, it is also possible to make some improvements and modifications, these improvements and modifications also should be regarded as
Protection scope of the present invention.
Claims (7)
1. the Correlation Identification system of a frequency segmentation, it is characterised in that including: first frequency module resolution (1), the first letter
Number send system (2), first the earth system (3), the first signal receiving module (4), the first Correlation Identification module (5), the second phase
Close identification module (6), the first result treatment module (7) and the first amplitude spectrum identification module (8)；
Described first frequency module resolution (1) is connected with the first signal transmitting module (2), the first signal transmitting module (2) point
It is not connected with the first the earth system (3), the first Correlation Identification module (5) and the second Correlation Identification module (6)；First the earth system
(3) it is connected with the first signal receiving module (4)；First signal receiving module (4) respectively with the first Correlation Identification module (5) and
Two Correlation Identification modules (6) be connected, the first Correlation Identification module (5) respectively with the second Correlation Identification module (6) and the first result
Processing module (7) is connected；First result treatment module (7) is connected with the first amplitude spectrum identification module (8).
2. the Correlation Identification system of frequency segmentation as claimed in claim 1, it is characterised in that:
Described first frequency module resolution (1) includes the first chip width operator unit (11), the second chip connected side by side
Width unit (12) and the first signal length operator unit (13) being sequentially connected, secondary signal length operator's unit (14)；
Described first signal transmitting module (2) includes the first signal transmitting unit (21) and secondary signal transmitting element (22)；
Described first signal receiving module (4) includes that the first reception electrode unit (41) and first gathers memory element (42) and the
Two gather memory element (43)；
Described first Correlation Identification module (5) include the first auto-correlation operator unit (51), the first cross-correlation operator unit (52),
First auto-correlation frequency domain unit (53), the first cross-correlation frequency domain unit (54) and the first system function solve unit (55)；
Described second Correlation Identification module (6) include the second auto-correlation operator unit (61), the second cross-correlation operator unit (62),
Second auto-correlation frequency domain unit (63), the second cross-correlation frequency domain unit (64) and second system function solve unit (65)；
Described first result treatment module (7) includes the first system frequency domain memory element (71), second system frequency domain memory element
(72), the first cut-off frequency operator unit (73) and the second cut-off frequency operator unit (74)；
Described first amplitude spectrum identification module (8) includes the first frequency spectrum concatenation unit (81) and the first amplitude spectrum identification unit
(82)；
First chip width operator unit (11) is connected with the first signal length operator unit (13), the first signal length operator list
Unit (13) is connected with the first signal transmitting unit (21)；Second chip width operator unit (12) is single with secondary signal length operator
Unit (14) is connected, and secondary signal length operator's unit (14) is connected with secondary signal transmitting element (22)；First chip width is calculated
Subelement (11) and the first chip width operator unit (12) are for arranging different for different the earth system investigation depths
Chip width；First signal length operator unit (13) and secondary signal length operator's unit (14) are for calculating the transmission of correspondence
Signal length, meanwhile, sends the first signal transmitting module (2) to arrange transmission by the chip width of generation and transmission signal length
The parameter of signal；
First signal transmitting unit (21) is connected with the first the earth system (3) respectively with secondary signal transmitting element (22)；According to
The different signal parameters that send that first frequency module resolution (1) is arranged produce different m-sequence, are sent single by the first signal
Transmission signal is sent to the first the earth system (3) by unit (21) and secondary signal transmitting element (22), thus completes signal
Send；
First receives electrode (41) gathers memory element (42) and the second collection storage respectively with the first the earth system (3) and first
Unit (43) be connected, first gather memory element (42) respectively with the first auto-correlation operator unit (51) and the first cross-correlation operator
Unit (52) be connected, second gather memory element (43) respectively with the first auto-correlation operator unit (61) and the first cross-correlation operator
Unit (62) is connected；
First receives electrode (41) receives voltage signal through the first the earth system (3), gathers memory element (42) through first
With second gather memory element (43) gather storage after, be transferred to respectively the first Correlation Identification module (5) and second be correlated with distinguish
Know module (6), thus complete signal receiving module；
First auto-correlation operator unit (51) is connected with the first auto-correlation frequency domain unit (53)；First cross-correlation operator unit (52)
It is connected with the first cross-correlation frequency domain unit (54), the first auto-correlation frequency domain unit (53) and the first cross-correlation frequency domain unit (54) point
Do not solve unit (55) with the first system function to be connected；First auto-correlation operator unit (51) passes through the first signal transmitting unit
(21) signal sent obtains autocorrelative frequency-domain result, the first cross-correlation operator list through the first auto-correlation frequency domain unit (53)
The transmission signal and first that unit (52) is obtained by the first transmitting element (21) gathers the reception signal that memory element (42) obtains
Obtain the frequency-domain result of cross-correlation through the first cross-correlation frequency domain unit (54), then solve unit (55) by the first system function
Just obtain the frequency domain identification result of the earth system of the corresponding degree of depth；
Second auto-correlation operator unit (61) is connected with the second auto-correlation frequency domain unit (63)；Second cross-correlation operator unit (62)
Being connected with the second cross-correlation frequency domain unit (64), the second auto-correlation frequency domain unit (63) and the second cross-correlation frequency domain unit (64) are equal
Solve unit (65) to be connected with second system function；Second auto-correlation operator unit (61) passes through secondary signal transmitting element (22)
The signal sent obtains autocorrelative frequency-domain result, the second cross-correlation operator unit through the second auto-correlation frequency domain unit (63)
(62) the transmission signal obtained by the second transmitting element (22) and second gathers the reception signal warp that memory element (43) obtains
Cross the second cross-correlation frequency domain unit (64) and obtain the frequency-domain result of cross-correlation, more just solve unit (65) by second system function
Obtain the frequency domain identification result of the earth system of the corresponding degree of depth；
The first system frequency domain memory element (71) is connected with the first cut-off frequency operator unit (73), and the storage of second system frequency domain is single
Unit (72) is connected with the second cut-off frequency operator unit (71)；The first system function solves the corresponding degree of depth that unit (55) obtains
During greatly system frequency-domain result is stored in the first system frequency domain memory element (71), through the first cut-off frequency operator unit (73)
Obtain this degree of depth correspondence cut-off frequency；Second system function solves the earth system frequency domain knot of the corresponding degree of depth that unit (65) obtains
Fruit is stored in second system frequency domain memory element (72), obtains this degree of depth through the second cut-off frequency operator unit (74) corresponding
Cut-off frequency, thus obtain the system frequency domain identification of different different frequency corresponding to investigation depth；
First frequency spectrum concatenation unit (81) is connected with the first amplitude spectrum identification unit (82)；First frequency spectrum concatenation unit (81) is by
Frequency domain identification result obtained by one cut-off frequency operator unit (73) and the second cut-off frequency operator unit (74) is spliced,
Obtain the frequency domain identification result of the earth system, obtain big system frequency domain identification result through the first amplitude spectrum identification unit (81)
Amplitude spectrum.
The Correlation Identification system of frequency segmentation the most according to claim 2, it is characterised in that:
When detecting the earth system of different depth, according to skin depth computing formula, the unlike signal that different depth is corresponding is sent out
Send frequency but the identical transmission signal of frequency resolution system, by the system identification to the earth, obtain different investigation depth
The greatly frequency spectrum of system function, and frequency spectrum is spliced, it is thus achieved that identification result.
The Correlation Identification system of frequency segmentation the most according to claim 3, it is characterised in that:
Described first frequency module resolution (1), the first Correlation Identification module (5), the second Correlation Identification module (6), the first knot
Really processing module (7) and amplitude spectrum identification module (8) are splicing frequency obtains new the earth system amplitude spectrum identification when
Implement process as follows: first, by the first frequency domain resolution module (1), transmission signal parameter is set, according to formulaWherein Δ F is frequency resolution, and Δ is chip width, and N is signal length, arranges chip width and signal length
Two of the first signal transmitting module (2) transmission signals are made to have identical frequency resolution, through the first the earth system (3) and the
One signal receiving module (4) calculates two respectively by the first auto-correlation operator unit (51) and the second auto-correlation operator unit (61)
The auto-correlation function R of individual signalxxT (), by the first cross-correlation operator unit (52) and the second cross-correlation operator unit (62) point
Ji Suan the cross-correlation function R of two signalsxyT (), then utilizes auto-correlation frequency domain unit and cross-correlation frequency domain unit to obtain it
Frequency domain from composing Rxx(ω) with cross-spectrum Rxy(ω), solve unit by system function and calculate the system function H of the earth system
(ω), computing formula H (ω)=Rxy(ω)/Rxx(ω)；
The system function obtained is calculated transmission letter through the first result treatment module (7) by the first cut-off frequency operator unit
The cut-off frequency f of number a1, the cut-off frequency f sending signal b is calculated by the second cut-off frequency operator unit2, then spelled by frequency spectrum
The frequency domain of two signals is spliced by order unit (81), obtains amplitude spectrum identification result, and leaves the first amplitude spectrum identification in
In unit (82).
5. according to the Correlation Identification system of the arbitrary described frequency segmentation of claim 1-4, it is characterised in that:
When calculating the earth system of different depth, necessarily change chip width and signal length according to frequency resolution, logical
The frequency-domain result of different cut-off frequencies is processed by the method crossing splicing, obtains new amplitude spectrum identification result.
6. according to the Correlation Identification system of the arbitrary described frequency segmentation of claim 1-5, it is characterised in that:
The parameter of the m-sequence of the first signal transmitting unit (21) arranges and is defaulted as: interpolation number is 4, and exponent number is 12 rank, and amplitude is
1V；The parameter of the m-sequence of secondary signal transmitting element (22) arranges and is defaulted as: interpolation number is 2, and exponent number is 13 rank, and amplitude is 1V.
7., according to the Correlation Identification method of the Correlation Identification system of the arbitrary described frequency segmentation of claim 1-6, its feature exists
In:
Step 201: power on；
Step 202: initialization operation；
Step 203: according to investigation depth, sends signal frequency, and noise jamming difference calculates chip width；
Step 204: according to formula signal calculated length:
Computing formula isΔ F is frequency resolution, and Δ is chip width, and N is signal length；
Step 205: send the signal setting parameter；
Step 206: signal is by the earth system:
Step 207: receive the signal by the earth system gained；
Step 208: the auto-correlation function R of signal calculatedxx(t):
Formula is
Step 209: calculate the auto-correlation function R of the signal received by step 207xy(t):
Formula is
Step 210: obtain auto-correlation function R by Fourier transformationxx(t) certainly compose Rxx(ω) with cross-correlation function Rxy(t)
Cross-spectrum Rxy(ω)；
Step 211: obtain from composing R completing step 210xx(ω) with cross-spectrum Rxy(ω) after, solving system function, calculate the earth system
The system function H (ω) of system:
Computing formula H (ω)=Rxy(ω)/Rxx(ω)；
Step 212: obtain the cut-off frequency of signal by the different signals that sends；
Step 213: connect the cut-off frequency of two signals, splices the frequency domain of two signals:
Step 214: data derive；
Step 215: complete.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610556834.XA CN106054269B (en) | 2016-07-14 | 2016-07-14 | A kind of Correlation Identification system and method for frequency segmentation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610556834.XA CN106054269B (en) | 2016-07-14 | 2016-07-14 | A kind of Correlation Identification system and method for frequency segmentation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106054269A true CN106054269A (en) | 2016-10-26 |
| CN106054269B CN106054269B (en) | 2019-01-08 |
Family
ID=57186913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610556834.XA Expired - Fee Related CN106054269B (en) | 2016-07-14 | 2016-07-14 | A kind of Correlation Identification system and method for frequency segmentation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106054269B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5748471A (en) * | 1996-03-29 | 1998-05-05 | Otatco, Inc. | Well collar identification method |
| CN101881638A (en) * | 2009-05-12 | 2010-11-10 | 石家庄西比克仪表有限公司 | Vortex street flowmeter signal digital processing system based on MCU (Microprogrammed Control Unit) |
| CN103257306A (en) * | 2012-10-25 | 2013-08-21 | 国网电力科学研究院武汉南瑞有限责任公司 | Method for diagnosing direct current partial discharging insulation state of converter transformer and measurement system |
| CN103649780A (en) * | 2011-05-13 | 2014-03-19 | 沙特阿拉伯石油公司 | Couple time-distance dependent swept frequency source acquisition design and data de-noising |
| CN104794313A (en) * | 2015-05-18 | 2015-07-22 | 甘肃省有色地质调查院 | Method and device for obtaining frequency response function of system to be identified |
| CN105607131A (en) * | 2016-01-11 | 2016-05-25 | 甘肃省有色地质调查院 | Method and apparatus for obtaining full-region apparent resistivity spectrum of ground by using coding source electromagnetic sounding method |
-
2016
- 2016-07-14 CN CN201610556834.XA patent/CN106054269B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5748471A (en) * | 1996-03-29 | 1998-05-05 | Otatco, Inc. | Well collar identification method |
| CN101881638A (en) * | 2009-05-12 | 2010-11-10 | 石家庄西比克仪表有限公司 | Vortex street flowmeter signal digital processing system based on MCU (Microprogrammed Control Unit) |
| CN103649780A (en) * | 2011-05-13 | 2014-03-19 | 沙特阿拉伯石油公司 | Couple time-distance dependent swept frequency source acquisition design and data de-noising |
| CN103257306A (en) * | 2012-10-25 | 2013-08-21 | 国网电力科学研究院武汉南瑞有限责任公司 | Method for diagnosing direct current partial discharging insulation state of converter transformer and measurement system |
| CN104794313A (en) * | 2015-05-18 | 2015-07-22 | 甘肃省有色地质调查院 | Method and device for obtaining frequency response function of system to be identified |
| CN105607131A (en) * | 2016-01-11 | 2016-05-25 | 甘肃省有色地质调查院 | Method and apparatus for obtaining full-region apparent resistivity spectrum of ground by using coding source electromagnetic sounding method |
Non-Patent Citations (2)
| Title |
|---|
| 何羚等: "基于定点DSP的高性能FFT谱估计", 《电子科技大学学报》 * |
| 周远田: "测井信号分辨率的频率域匹配", 《测井技术》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109444961A (en) * | 2018-09-25 | 2019-03-08 | 中国地质大学（北京） | A kind of multithreading Correlation Identification electrical prospecting reception system and method based on all phase |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106054269B (en) | 2019-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102565634B (en) | Power cable fault location method based on transfer function method | |
| Garg et al. | Geo-location estimation from electrical network frequency signals | |
| CN102356671B (en) | Method and device of communication system frequency search | |
| CN103454495B (en) | Self-adaptation high-precision fast spectral analytical approach | |
| CN101374131B (en) | Method and apparatus of timing synchronization, leading symbol as well as method and apparatus for generating the same | |
| CN104378764B (en) | A kind of LTE system co-frequency cell detection method and device | |
| CN102353877A (en) | Partial discharge detection method for realizing phase resolved partial discharge (PRPD) pattern and time resolved partial discharge (TRPD) pattern | |
| CN102508031A (en) | Fourier series based measurement method of phase angle of partial discharge pulse | |
| CN104391330A (en) | Relevant identification technology based electrical prospecting system | |
| CN103428819A (en) | Carrier frequency point searching method and device | |
| Spagnol | On distance dependence of pinna spectral patterns in head-related transfer functions | |
| CN106054269A (en) | Frequency partition correlation identification system and method thereof | |
| CN105791194A (en) | Synchronization method and system for anti-narrowband interference | |
| CN106885975A (en) | A kind of high-tension cable ageing testing method and device based on impulse response | |
| CN106033671A (en) | Method and device for determining inter-channel time difference parameter | |
| AU2021289742B2 (en) | Methods, apparatus, and systems for detection and extraction of spatially-identifiable subband audio sources | |
| CN101594168A (en) | A kind of lead code detecting method of Random Access Channel, device and equipment | |
| CN103427915B (en) | A kind of radio-frequency apparatus standing-wave ratio detection in remove interference method and device | |
| CN104199064B (en) | Correlation fluctuation based BOC (binary offset carrier) signal parameter blind estimation method | |
| CN109190605A (en) | A kind of human body continuous action method of counting based on CSI | |
| CN103944845A (en) | Related detection method for full-frequency-spectrum carrier modulation | |
| CN106772193A (en) | A kind of current transformer frequency characteristic measurement devices and methods therefor | |
| CN105680964A (en) | Spectrum sensing method, spectrum sensing system, client and server | |
| CN107329123A (en) | A kind of weak radar pulse envelope signal detecting method and device | |
| CN110190917A (en) | A kind of frequency spectrum cavity-pocket cognitive method, device and the equipment of LTE230MHz electric power wireless private network |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190108 Termination date: 20200714 |