CN104794313B - A kind of method and device for obtaining unidentified system frequency response function - Google Patents
A kind of method and device for obtaining unidentified system frequency response function Download PDFInfo
- Publication number
- CN104794313B CN104794313B CN201510253306.2A CN201510253306A CN104794313B CN 104794313 B CN104794313 B CN 104794313B CN 201510253306 A CN201510253306 A CN 201510253306A CN 104794313 B CN104794313 B CN 104794313B
- Authority
- CN
- China
- Prior art keywords
- signal
- input signal
- response function
- mrow
- frequency response
- 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
Landscapes
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The present invention provides a kind of method and device for obtaining unidentified system frequency response function.This method includes：Linear feedback shift register is triggered by clock signal and generates input signal, reference signal is generated according to the time series cycle of input signal and sample rate triggering software program, input signal and reference signal are inverse repetition M sequence and the cycle is the even-multiple of power frequency period；Input signal, the output signal of collection unidentified system output are inputted to unidentified system；To reference signal, input signal and output signal application Cross-Correlation Theory, the frequency response function of unidentified system is obtained.Because input signal is multicycle inverse repetition M sequence, its frequency spectrum avoids sequentially inputting multiple single frequency signals comprising multiple discrete frequencies, improves the efficiency for obtaining frequency response function.Due to adding reference signal so that obtain that the identification precision of frequency response function is very high, the cycle of input signal and reference signal is the even-multiple of power frequency period, can effectively suppress Hz noise.
Description
Technical field
The present invention relates to Signal acquiring and processing field, System Discrimination field and geophysics electromagnetic prospecting field, tool
For body, it is related to a kind of method and device for obtaining unidentified system frequency response function.
Background technology
At present, in order to understand the system performance of unidentified system, generally require to obtain the frequency response letter of unidentified system
Number, unidentified system can be circuit system, observation instrument system or the earth system to be detected etc..For example, wait to visit to determine
The geology distribution of the earth system is surveyed, and then the ore body in the earth system to be detected is reconnoitred, it is necessary to obtain the earth to be detected
The electromagnetic frequency receptance function of system.
Currently, prior art provides a kind of method for obtaining unidentified system frequency response function, including：With certain frequency
The input signal of rate excites unidentified system, and response output of the collection unidentified system to the input signal is output signal.So
Change the frequency of input signal afterwards, unidentified system, collection unidentified system output are excited with the input signal after change frequency
Output signal.So input the input signal of different frequency to unidentified system successively according to the method described above, and gather respectively
Output signal corresponding to each input signal.Each input signal and its corresponding output signal are analyzed and processed, obtained
The frequency response function of unidentified system.
During the present invention is realized, inventor has found problems with the prior art at least be present：
Multiple input signals are inputted to unidentified system successively and gather output signal, it is necessary to devote a tremendous amount of time, lead
Cause the efficiency for obtaining frequency response function very low.And the industrial ac magnetic field being widely present in modern society can influence input letter
Number and output signal, cause in the serious area of Hz noise, it is difficult to accurately obtain system frequency response function.
The content of the invention
In view of this, the purpose of the embodiment of the present invention is to provide a kind of side for obtaining unidentified system frequency response function
Method and device, the efficiency for obtaining frequency response function is improved, and improve the precision for obtaining frequency response function.
In a first aspect, the embodiments of the invention provide a kind of method for obtaining unidentified system frequency response function, it is described
Method includes：
Linear feedback shift register is triggered by clock signal and generates input signal, according to the time of the input signal
Sequence period and sample rate triggering software program generation reference signal, the input signal and the reference signal are inverse repeat
M sequence and cycle phase are the even-multiple of power frequency period together；
The input signal is inputted to unidentified system, gathers the output signal of the unidentified system output；
To the reference signal, the input signal and the output signal application Cross-Correlation Theory, wait to distinguish described in acquisition
The frequency response function of knowledge system.
With reference in a first aspect, the embodiments of the invention provide the possible implementation of the first of above-mentioned first aspect, its
In, it is described to the reference signal, the input signal and the output signal application Cross-Correlation Theory, obtain described to be identified
The frequency response function of system, including：
The output signal and the input signal are subjected to circulation cross-correlation with the reference signal respectively, obtained described
Cross-correlation time series corresponding to cross-correlation time series corresponding to output signal and the input signal；
To cross-correlation time series corresponding to cross-correlation time series corresponding to the output signal and the input signal
FFT is carried out respectively, obtains the crosspower spectrum of the output signal and the crosspower spectrum of the input signal；
According to the crosspower spectrum of the output signal and the crosspower spectrum of the input signal, the unidentified system is obtained
Frequency response function.
With reference to the first possible implementation of first aspect, the embodiments of the invention provide the of above-mentioned first aspect
Two kinds of possible implementations, wherein, it is described according to the crosspower spectrum of the output signal and the cross-power of the input signal
Spectrum, the frequency response function of the unidentified system is obtained, including：
Obtain output spectral line peak value and its corresponding frequency values from the crosspower spectrum of the output signal, and from described
Input spectral line peak value and its corresponding frequency values are obtained in the crosspower spectrum of input signal；
Frequency values identical is filtered out from the output spectral line peak value of the acquisition and the input spectral line peak value of the acquisition
Spectral line peak value and input spectral line peak value are exported, and the frequency values identical is calculated by equation below (1) and exports spectral line peak value
Ratio between input spectral line peak value, obtains the frequency response function of the unidentified system；
Wherein, in formula (1), ω is frequency values,To export spectral line peak value corresponding to frequency values ω,
To input spectral line peak value corresponding to frequency values ω, H (ω) is the frequency response function of the unidentified system.
With reference to the first possible implementation of first aspect, the embodiments of the invention provide the of above-mentioned first aspect
Three kinds of possible implementations, wherein, it is described to input the input signal to unidentified system, it is defeated to gather the unidentified system
The output signal gone out, including：
The input signal is inputted to the unidentified system, and in a synchronous manner or asynchronous system collection is described to be identified
The output signal of system output；
Accordingly, it is described to the reference signal, the input signal and the output signal application Cross-Correlation Theory, obtain
The frequency response function of the unidentified system is taken, including：
When gathering the output signal with the asynchronous system, according to the crosspower spectrum of the output signal and described defeated
Enter the crosspower spectrum of signal, obtain the frequency response function of the unidentified system, the frequency response function includes amplitude spectrum；
When gathering the output signal with the method for synchronization, according to the crosspower spectrum of the output signal and described defeated
Enter the crosspower spectrum of signal, obtain the frequency response function of the unidentified system, the frequency response function includes phase spectrum
And amplitude spectrum, and by equation below (2), obtain the phase spectrum of the unidentified system frequency response function；
Wherein, in formula (2), mod is modulus operator,For the mutual of the output signal corresponding to frequencies omega
The phase of power spectrum,For the phase of the crosspower spectrum of the input signal corresponding to frequencies omega,Wait to distinguish to be described
Know the phase spectrum of system frequency response function.
With reference in a first aspect, the embodiments of the invention provide the possible implementation of the 4th of above-mentioned first aspect kind, its
In, after the frequency response function for obtaining the unidentified system, in addition to：
When the unidentified system is the earth system to be detected, Inversion Calculation is carried out to the frequency response function, obtained
The geoelectric cross section of the earth system to be detected is taken, and the electricity of the earth system to be detected is calculated according to the frequency response function
Magnetic explores parameter, and the electromagnetic survey parameter includes step response, complex resistivity, percentage frequency effect and late period apparent resistance
Rate.
Second aspect, a kind of device for obtaining unidentified system frequency response function, described device include：
Generation module, input signal is generated for triggering linear feedback shift register by clock signal, according to described
The time series cycle of input signal and sample rate triggering software program generation reference signal, the input signal and the reference
Signal is against repetition M sequence and cycle phase is the same as the even-multiple for being power frequency period；
Acquisition module, for inputting the input signal to unidentified system, gather the defeated of the unidentified system output
Go out signal；
Acquisition module, for the reference signal, the input signal and the output signal application Cross-Correlation Theory,
Obtain the frequency response function of the unidentified system.
With reference to second aspect, the embodiments of the invention provide the possible implementation of the first of above-mentioned second aspect, its
In, the acquisition module includes：
Cross-correlation unit is circulated, for the output signal and the input signal to be carried out with the reference signal respectively
Cross-correlation is circulated, obtains the cross-correlation time corresponding to cross-correlation time series corresponding to the output signal and the input signal
Sequence；
Fast Fourier transform unit, for believing cross-correlation time series corresponding to the output signal and the input
Cross-correlation time series carries out FFT respectively corresponding to number, obtains the crosspower spectrum of the output signal and described
The crosspower spectrum of input signal；
First acquisition unit, for the crosspower spectrum according to the output signal and the crosspower spectrum of the input signal,
Obtain the frequency response function of the unidentified system.
With reference to the first possible implementation of second aspect, the embodiments of the invention provide the of above-mentioned second aspect
Two kinds of possible implementations, wherein, the first acquisition unit includes：
Subelement is obtained, for obtaining output spectral line peak value and its corresponding frequency from the crosspower spectrum of the output signal
Rate value, and input spectral line peak value and its corresponding frequency values are obtained from the crosspower spectrum of the input signal；
Computation subunit, for being screened from the input spectral line peak value of the output spectral line peak value of the acquisition and the acquisition
Go out frequency values identical output spectral line peak value and input spectral line peak value, and it is identical by equation below (1) the calculating frequency values
Output spectral line peak value and input spectral line peak value between ratio, obtain the frequency response function of the unidentified system；
Wherein, in formula (1), ω is frequency values,To export spectral line peak value corresponding to frequency values ω,
To input spectral line peak value corresponding to frequency values ω, H (ω) is the frequency response function of the unidentified system.
With reference to the first possible implementation of second aspect, the embodiments of the invention provide the of above-mentioned second aspect
Three kinds of possible implementations, wherein, the acquisition module, for inputting the input signal to the unidentified system, and
In a synchronous manner or asynchronous system gathers the output signal of the unidentified system output；
Accordingly, the acquisition module includes：
Second acquisition unit, for when with the asynchronous system gather the output signal when, according to the output signal
Crosspower spectrum and the input signal crosspower spectrum, obtain the frequency response function of the unidentified system, the frequency
Receptance function includes amplitude spectrum；
3rd acquiring unit, for when with the method for synchronization gather the output signal when, according to the output signal
Crosspower spectrum and the input signal crosspower spectrum, obtain the frequency response function of the unidentified system, the frequency
Receptance function includes phase spectrum and amplitude spectrum, and by equation below (2), obtains the unidentified system frequency response function
Phase spectrum；
Wherein, in formula (2), mod is modulus operator,For the mutual of the output signal corresponding to frequencies omega
The phase of power spectrum,For the phase of the crosspower spectrum of the input signal corresponding to frequencies omega,Wait to distinguish to be described
Know the phase spectrum of system frequency response function.
With reference to second aspect, the embodiments of the invention provide the possible implementation of the 4th of above-mentioned second aspect kind, its
In, described device also includes：
Computing module, for when the unidentified system is the earth system to be detected, entering to the frequency response function
Row Inversion Calculation, the geoelectric cross section of the earth system to be detected is obtained, and according to calculating the frequency response function
The electromagnetic survey parameter of the earth system to be detected, the electromagnetic survey parameter include step response, complex resistivity, percentage frequency
Effect and late period apparent resistivity.
In embodiments of the present invention, because input signal is inverse repetition M sequence, multiple discrete frequencies are included in its frequency spectrum
Signal, the signal of this multiple discrete frequency amplitude approximately equal in main band, so avoiding defeated successively to unidentified system
Enter multiple single frequency signals, improve the efficiency for obtaining frequency response function.In addition, according to reference signal, output signal and defeated
Entering signal, the precision of frequency response function is very high obtaining by cross-correlation method, and the cycle of input signal and reference signal is work
The even-multiple in frequency cycle, it can effectively suppress Hz noise.
To enable the above objects, features and advantages of the present invention to become apparent, preferred embodiment cited below particularly, and coordinate
Appended accompanying drawing, is described in detail below.
Brief description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below by embodiment it is required use it is attached
Figure is briefly described, it will be appreciated that the following drawings illustrate only certain embodiments of the present invention, therefore be not construed as pair
The restriction of scope, for those of ordinary skill in the art, on the premise of not paying creative work, can also be according to this
A little accompanying drawings obtain other related accompanying drawings.
Fig. 1 shows a kind of method stream for acquisition unidentified system frequency response function that the embodiment of the present invention 1 is provided
Cheng Tu；
Fig. 2A shows a kind of method stream for acquisition unidentified system frequency response function that the embodiment of the present invention 2 is provided
Cheng Tu；
Fig. 2 B show a kind of inverse time-domain waveform diagram for repeating M sequence that the embodiment of the present invention 2 is provided；
Fig. 2 C show a kind of inverse power spectrum schematic diagram for repeating M sequence that the embodiment of the present invention 2 is provided；
Fig. 2 D show a kind of signal for inverse auto-correlation function curve for repeating M sequence that the embodiment of the present invention 2 is provided
Figure；
Fig. 2 E show the signal of the crosspower spectrum of a kind of output signal and reference signal that the embodiment of the present invention 2 is provided
Figure；
Fig. 3 shows the method flow diagram for the acquisition unidentified system frequency response function that the embodiment of the present invention 3 is provided；
Fig. 4 shows a kind of device knot for acquisition unidentified system frequency response function that the embodiment of the present invention 4 is provided
Structure schematic diagram.
Embodiment
Below in conjunction with accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Ground describes, it is clear that described embodiment is only part of the embodiment of the present invention, rather than whole embodiments.Generally exist
The component of the embodiment of the present invention described and illustrated in accompanying drawing can be configured to arrange and design with a variety of herein.Cause
This, the detailed description of the embodiments of the invention to providing in the accompanying drawings is not intended to limit claimed invention below
Scope, but it is merely representative of the selected embodiment of the present invention.Based on embodiments of the invention, those skilled in the art are not doing
The every other embodiment obtained on the premise of going out creative work, belongs to the scope of protection of the invention.
In view of inputting multiple input signals to unidentified system successively in prior art and gathering output signal, it is necessary to flower
Take the substantial amounts of time, cause the efficiency for obtaining frequency response function very low.And industrial ac frequency can influence in correlation technique
The frequency of input signal, cause the precision of the frequency response function of acquisition very low.Based on this, the embodiments of the invention provide one kind
Obtain the method and device of unidentified system frequency response function.It is described below by embodiment.
Embodiment 1
Referring to Fig. 1, the embodiments of the invention provide a kind of method for obtaining unidentified system frequency response function, this method
It can be performed by the device of acquisition unidentified system frequency response function.This method specifically includes following steps：
Step 101：Linear feedback shift register is triggered by clock signal and generates input signal, according to input signal
Time series cycle and sample rate triggering software program generation reference signal, input signal and reference signal are inverse repetition M sequences
Arrange and cycle phase is the same as the even-multiple for being power frequency period；
Step 102：The input signal, the output signal of collection unidentified system output are inputted to unidentified system；
Step 103：To reference signal, input signal and output signal application Cross-Correlation Theory, unidentified system is obtained
Frequency response function.
In embodiments of the present invention, because input signal is inverse repetition M sequence, multiple discrete frequencies are included in its frequency spectrum
Signal, the signal of this multiple discrete frequency amplitude approximately equal in main band, so avoiding defeated successively to unidentified system
Enter multiple single frequency signals, improve the efficiency for obtaining frequency response function.In addition, according to reference signal, output signal and defeated
Entering signal, the precision of frequency response function is very high obtaining by cross-correlation method, and the cycle of input signal and reference signal is work
The even-multiple in frequency cycle, it can effectively suppress Hz noise.
Embodiment 2
The embodiments of the invention provide a kind of method for obtaining unidentified system frequency response function, this method can be by obtaining
The device of unidentified system frequency response function is taken to perform.
Wherein, unidentified system is circuit system, observation instrument system or the earth system to be detected etc..Due to system to be identified
Comprising the medium that electric conductivity is different with magnetic conductivity in system, unidentified system is excited with input signal, due to electromagnetic induction be present, is treated
Identification system can produce response to input signal.Gather unidentified system output in a synchronous manner in embodiments of the present invention
Output signal, and obtain by method provided in an embodiment of the present invention the frequency response function of unidentified system.
Referring to Fig. 2A, this method specifically includes following steps：
Step 201：Linear feedback shift register is triggered by clock signal and generates input signal, the input signal is inverse
Repetition M sequence and cycle are the even-multiple of power frequency period；
Wherein, input signal is periodically pseudorandom inverse repetition M sequence.
The executive agent of the embodiment of the present invention is the device for obtaining unidentified system frequency response function, system to be identified
System can be circuit system, observation instrument system or the earth system to be detected etc..N+1 bit registers are provided with the device, N is
Integer more than or equal to 3.In this N+1 bit register, 1 bit register is used to generate two divided-frequency shift register clock letter
Number, remaining N-bit register forms linear feedback shift register, presets the feedback bit of linear feedback shift register,
This N positions linear feedback shift register constitutes maximum length pseudorandom m-sequence generation circuit.Above-mentioned another 1 bit register two divides
Frequency linear feedback shift register clock signal, by the clock signal and the pseudorandom m of N positions linear feedback shift register generation
Sequence carries out XOR, realizes and pseudorandom m-sequence is negated every position, obtains inverse repetition M sequence.Sample rate is preset, is pressed
The inverse time series for repeating M sequence is sampled according to the sample rate, the time series signal that sampling obtains is defined as inputting
Signal, sample rate are the frequency of collection input signal.Wherein, input signal is the signal of multicycle.
Wherein, the inverse cycle for repeating M sequence and frequency band be adjustable, can be by selecting different linear feedback shifts
Register digit timely clock frequency is adjusted.The digit of linear feedback shift register typically takes 5,6,7,8,9 or 10, clock
Frequency can be selected in frequency separation [0.5Hz, 10kHz].Wherein, as shown in table 1, different linear feedback shifts is posted
Storage digit and clock frequency, the corresponding inverse cycle for repeating M sequence and sequence bit length are different.Primitive polynomial is production in table 1
The raw inverse feedback bit expression formula for repeating M sequence, clock frequency in the inverse corresponding table for repeating the M sequence cycle with simply schematically giving
Go out clock frequency, clock frequency can be selected in frequency separation [0.5Hz, 10kHz] in concrete operations.
Table 1
Wherein, the inverse cycle for repeating M sequence is 2 times of cycle of maximum linear feedback pseudorandom m-sequence.Take inverse repetition M
The cycle of sequence is the even-multiple of power frequency period, i.e. the cycle of input signal is the even-multiple of power frequency period.Commercial power compared with
Under more scenes, input signal can be influenceed by commercial power.And the correlation peak of input signal show as periodically it is positive and negative
It is alternately present, therefore takes the cycle of input signal effectively to suppress Hz noise for the even-multiple of power frequency period, weaken work
Influence of the frequency interference to input signal, and then improve the accuracy of the frequency response function obtained.
Wherein, power frequency is the frequency of AC power industrially used, generally 50Hz.
Wherein, as shown in Figure 2 B, waveform is random to the inverse time domain waveform for repeating M sequence in inverse repetition 1 cycle of M sequence
's.As shown in Figure 2 C, its power spectrum is discrete line spectrum to the inverse power spectrum for repeating M sequence.The inverse auto-correlation letter for repeating M sequence
As shown in Figure 2 D, the curve of the auto-correlation function shows as positive and negative alternate spike to number curve.
Wherein, input signal is the inverse signal for repeating M sequence, multiple discrete frequencies being contained in inverse repetition M sequence signal,
And this multiple signals amplitude in main band is equal, so obtaining frequency response function according to the input signal, can avoid
Multiple single frequency signals one by one are inputted to unidentified system, improve the efficiency for obtaining frequency response function.
Wherein, by after the operation generation input signal of this step, it is necessary to as follows 202 operation obtain
The output signal of unidentified system.
Step 202：The input signal is inputted to unidentified system, and gathers the defeated of unidentified system output in a synchronous manner
Go out signal；
Wherein, the method for synchronization represents, according to gathering input signal identical sample rate and sampling time, to pass through performance one
The sampler of cause gathers the output signal exported after unidentified system responds to input signal with keeping strokes.The output is believed
Number can be electric signal or magnetic signal.
Wherein, the device for obtaining unidentified system frequency response function can be split type apparatus for transmitting signal and letter
Number receiving device or send and receive integral type collecting device.Apparatus for transmitting signal can include high-power inverter
Drive and exported by electric dipole by high-power inverter and make to unidentified system, signal receiver Deng, input signal
The electromagnetic response of unidentified system is received with electric dipole or magnetic dipole, electric dipole can be electrically isolated electrode etc.,
Magnetic dipole can be magnetic probe etc..Wherein, input signal is electric signal, and output signal has electric signal and the species of magnetic signal two
Type.Input signal is sent to unidentified system by apparatus for transmitting signal using inverter by electric dipole.In collection output letter
Number when, the output signal gathered by the electric dipole of signal collecting device is electric signal, passes through the magnetic couple of signal collecting device
The output signal of extremely son collection is magnetic signal.
This step is specifically, control signal transmission equipment is connected on unidentified system by electric dipole, and is controlled
Electric dipole or magnetic probe are connected on unidentified system by signal collecting device.Control signal sends equipment and leads to input signal
Cross electric dipole and input to unidentified system, excite unidentified system, at the same time, control signal collecting device by electrode or
Magnetic probe, the electromagnetic response of unidentified system is gathered to be defeated according to collection input signal identical sample rate and sampling time
Go out signal.
Wherein, in embodiments of the present invention, can be treated by electric dipole or the magnetic dipole collection of signal collecting device
The electric field response or magnetic responsiveness of identification system.In order to improve the precision of the frequency response function of acquisition, in the embodiment of the present invention
In, it is necessary to keep input signal collector consistent with output signal collector systems characteristic.
In addition, being additionally provided with time synchronism equipment in device for obtaining unidentified system frequency response function, control
The time synchronism equipment synchronizing signal sends equipment, and synchronous signal acquisition equipment, excites apparatus for transmitting signal to be identified
Signal collecting device gathers the output signal of unidentified system output simultaneously during system.Further, can also be in time synchronized
Sample rate is preset in equipment, by the time synchronism equipment control signal collecting device according to default sample rate come synchronous acquisition
Input signal and output signal.Sample rate is gathers the frequency of signal, i.e., every preset time period input signal of collection
The amplitude of amplitude and output signal.
Further, signal collecting device can also store the input signal and output signal of collection, and by collection
Input signal and output signal send the processing equipments such as computer to, make the processing equipments such as computer to the input signal of collection and
Output signal handled deeper into ground calculating.
Wherein, obtained by the operation of above-mentioned steps 201 and 202 after input signal and output signal, it is necessary to by as follows
The operation of step 203 generates reference signal.
Step 203：According to the time series cycle of input signal and sample rate triggering software program generation reference signal；
Wherein, the sequence length in 1 cycle of the sequence length in 1 cycle of reference signal and input signal is identical, at the beginning of it
Beginning state can be different from input signal.Reference signal only needs the inverse of 1 or 2 complete cycle to repeat M sequence, and input signal
For the even-multiple that multicycle inverse repeatedly M sequence and cycle are power frequency period.The generating mode of the reference signal can be believed with input
Number generating mode it is identical.Wherein, set the digit for producing the inverse linear feedback shift register for repeating M sequence and feedback bit it
Afterwards, the caused inverse M sequence that repeats is fixed, is simply had between inverse repetition M sequence caused by different dwell vessel original states
Having time postpones, but the spectrum distribution of caused inverse repetition M sequence is identical.
Wherein, in this step, the parameters such as clock frequency and the sample rate needed for reference signal are generated, are inputted with generation
The parameter of signal is identical.In addition, technical staff can develop the software program for generating reference signal, setting clock frequency in advance
Rate and linear feedback shift register digit, the software program is triggered according to the time series cycle of input signal and sample rate
1 or 2 the inverse of complete cycle of generation repeats M sequence, and the inverse M sequence that repeats of generation is defined as into reference signal.
Wherein, the cycle of the reference signal is identical with input signal cycle, is the even-multiple of power frequency period.Follow-up basis
The reference signal carries out circulation cross correlation process to input signal and output signal, can effectively suppress Hz noise, improves
The precision of the frequency response function of acquisition.
Wherein, this step operation can also perform simultaneously with the operation of step 201.
Wherein, can be with as follows 204 operation after getting input signal, output signal and reference signal
To obtain the frequency response function of unidentified system.
Step 204：To reference signal, input signal and output signal application Cross-Correlation Theory, unidentified system is obtained
Frequency response function；
Wherein, this step can be realized by following S1-S4 operation, be specially：
S1：According to the reference signal, input signal and output signal, the cross-correlation time is obtained corresponding to input signal respectively
Cross-correlation time series corresponding to sequence and output signal；
Specifically, according to reference signal and input signal, reference signal and input signal are carried out by equation below (3)
Cross-correlation is circulated, obtains cross-correlation time series corresponding to input signal.And according to reference signal and output signal, pass through
Reference signal and output signal are carried out circulation cross-correlation by equation below (4), obtain cross-correlation time sequence corresponding to output signal
Row.
Cus(t)=cxcorr (u, SS) ... (3)
Cys(t)=cxcorr (yout, SS) ... (4)
Wherein, in formula (3) and (4), u is input signal, and yout is output signal, and SS is reference signal, when t is
Between, Cus(t) it is cross-correlation time series, C corresponding to input signalys(t) it is cross-correlation time series corresponding to output signal,
Cxcorr (u, SS) expressions carry out circulation computing cross-correlation to input signal u and reference signal SS, and cxcorr (yout, SS) is represented
Circulation computing cross-correlation is carried out to output signal yout and reference signal SS.
S2：Cross-correlation time series corresponding to cross-correlation time series and output signal according to corresponding to input signal, point
Huo Qu not the crosspower spectrum of input signal and the crosspower spectrum of output signal；
Specifically, fast Flourier change is carried out to cross-correlation time series corresponding to input signal by equation below (5)
Change, obtain the crosspower spectrum of input signal.And cross-correlation time series corresponding to output signal is entered by equation below (6)
Row FFT, obtain the crosspower spectrum of output signal.
Wherein, in formula (5) and (6), ω is frequency values,For the crosspower spectrum of input signal,For
The crosspower spectrum of output signal, FFT (Cus(t)) represent to cross-correlation time series C corresponding to input signalus(t) carry out quick
Fourier transform, FFT (Cys(t)) represent to cross-correlation time series C corresponding to output signalys(t) fast Flourier change is carried out
Change.
Wherein, the crosspower spectrum of output signal and reference signal as shown in Figure 2 E, with reference signal followed by output signal
Ring cross-correlation obtains cross-correlation time series corresponding to output signal, then the cross-correlation time series of output signal is carried out quickly
Fourier transform, obtain the crosspower spectrum of output signal.Crosspower spectrum shown in Fig. 2 E is the crosspower spectrum of output signal.
Wherein, the crosspower spectrum of input signal and the crosspower spectrum of output signal are got by above-mentioned S1 and S2 operation
Afterwards, according to the crosspower spectrum of input signal and the crosspower spectrum of output signal, can be obtained by following S3 and S4 operation
The frequency response function of unidentified system.
S3：Obtain input spectral line peak value and its corresponding frequency values from the crosspower spectrum of input signal, and from output
Output spectral line peak value and its corresponding frequency values are obtained in the crosspower spectrum of signal；
Wherein, the crosspower spectrum that spectral line peak value is more than predetermined threshold value for the amplitude in the crosspower spectrum of input signal is inputted
Value, output spectral line peak value are more than the cross-power spectrum of above-mentioned predetermined threshold value for the amplitude in the crosspower spectrum of output signal.
This step from the spectral line of the crosspower spectrum of input signal specifically, determine that amplitude is more than all of predetermined threshold value
Spectral line maximum point, cross-power spectrum and frequency values corresponding to each maximum point of determination are obtained, by corresponding to maximum point
Cross-power spectrum is defined as inputting spectral line peak value, and frequency values corresponding to maximum point are defined as into the input spectral line peak value pair
The frequency values answered.Similarly, determine that amplitude is more than all spectrums of predetermined threshold value from the spectral line of the crosspower spectrum of output signal
Line maximum point, cross-power spectrum and frequency values corresponding to each maximum point of determination are obtained, will be mutual corresponding to maximum point
Power spectral value is defined as exporting spectral line peak value, and frequency values corresponding to maximum point are defined as into the output spectral line peak value and corresponded to
Frequency values.
S4：According to input spectral line peak value and its corresponding frequency values and output spectral line peak value and its corresponding frequency values,
Obtain the frequency response function of unidentified system.
Specifically, from input spectral line peak value and output spectral line peak value in filter out frequency values identical input spectral line peak value and
Spectral line peak value is exported, and frequency values identical is calculated by equation below (1) and exported between spectral line peak value and input spectral line peak value
Ratio, obtain the frequency response function of unidentified system.
Wherein, in formula (1), ω is frequency values,To export spectral line peak value corresponding to frequency values ω,
To input spectral line peak value corresponding to frequency values ω, H (ω) is the frequency response function of unidentified system.
Wherein,For the mutual work(of the crosspower spectrum of output signal, herein output signal corresponding to expression frequency values ω
Output spectral line peak value corresponding to rate spectrum, i.e. frequency values ω.For the crosspower spectrum of input signal, frequency values are represented herein
Input spectral line peak value corresponding to the cross-power spectrum of input signal corresponding to ω, i.e. frequency values ω.
Wherein, because the cycle of input signal and reference signal is the even-multiple of power frequency period, so being believed according to input
Number excite unidentified system to obtain output signal, and circulated respectively with reference signal according to the input signal and output signal
Cross-correlation, the frequency response function of acquisition are influenceed very little by Hz noise, and the precision of the frequency response function of acquisition is very high.
Wherein, can be in the frequency response function of acquisition simultaneously when gathering input signal and output signal in a synchronous manner
Reflect the amplitude spectrum and phase spectrum of unidentified system frequency response.Wherein, the amplitude spectrum reflection unidentified system of frequency response
The intensity size of response to input signal, frequency response phase spectrum reflection unidentified system response and input signal when
Between postpone size.So contained in the frequency response function obtained when gathering input signal and output signal in a synchronous manner
Amplitude spectrum and phase spectrum.
Wherein, in embodiments of the present invention, when gathering the output signal of unidentified system in a synchronous manner, can also pass through
The operation of following steps 205 directly to obtain the phase spectrum of the frequency response function of unidentified system.
Step 205：According to the crosspower spectrum of input signal and the crosspower spectrum of output signal, the frequency of unidentified system is obtained
The phase spectrum of rate receptance function；
Wherein, the crosspower spectrum for input signal and the crosspower spectrum of output signal, can directly from step 204 S1 and
S2 is calculated.
This step passes through equation below specifically, according to the crosspower spectrum of input signal and the crosspower spectrum of output signal
(2) phase spectrum of unidentified system frequency response function, is obtained.
Wherein, in formula (2), mod is modulus operator,For the crosspower spectrum of output signal corresponding to frequencies omega
Phase,For the phase of the crosspower spectrum of input signal corresponding to frequencies omega,Rung for the frequency of unidentified system
Answer the phase spectrum of function.
Wherein, the amplitude spectrum of the frequency response function of unidentified system is got by above-mentioned steps 201-205 operation
After phase spectrum, the system performance of unidentified system and medium can be analyzed according to the frequency response function of unidentified system
Situations such as composition.Because the frequency response function amplitude spectrum and the precision of phase spectrum that are obtained in the embodiment of the present invention are very high, so
According to frequency response function come the error very little analyzed unidentified system.
Wherein, in embodiments of the present invention, when unidentified system is the earth system to be detected, the frequency response letter of acquisition
Number is impedance spectrum greatly, carries out Inversion Calculation with the impedance spectrum, obtains the geoelectric cross section of the earth system to be detected, with
And the electromagnetic survey parameter of the earth system to be detected is calculated according to frequency response function, the electromagnetic survey parameter can be that step is rung
Should, complex resistivity, percentage frequency effect and late period apparent resistivity etc..
Specifically, the amplitude spectrum and phase spectrum of the earth system frequency response function to be detected picked out, reflection the earth are
The impedance information of system, Inversion Calculation is carried out with amplitude spectrum and phase spectrum, obtain the geoelectric cross section of the earth system to be detected.It is and logical
The complex resistivity that equation below (7) exports the earth system to be identified is crossed, and the earth system to be identified is exported by equation below (8)
The percentage frequency effect of system.And the step response of the earth system to be detected is calculated according to frequency response function, rung by step
The late period asymptotic value answered calculates the late period apparent resistivity of the earth system to be identified.
ρ (ω)=KZ (ω) ... (7)
Wherein, in formula (7), ρ (ω) is complex resistivity, and K is electrode coefficient, and Z (ω) is the earth system to be identified
Amplitude spectrum.In formula (8), PFE is percentage frequency effect, Z (ωL) it is low bin magnitudes, Z (ωH) it is high frequency points amplitude.
Wherein, also can be to the output signal because signal collecting device is when receiving the output signal of unidentified system
Responded, i.e., there is also frequency response, the frequency response can reduce the frequency of the unidentified system of acquisition for signal collecting device
The accuracy of receptance function.So in embodiments of the present invention, by entering to signal collecting device application the technology of the present invention method
Row demarcation in advance and debugging, the frequency response values for making signal collecting device are always 1, so as to eliminate the frequency of signal collecting device
The influence of response.So the frequency response obtained in embodiments of the present invention is only the frequency response of unidentified system in itself, frequency
The accuracy of rate receptance function is very high, and frequency resolution is very high.
In embodiments of the present invention, because input signal is inverse repetition M sequence, multiple discrete frequencies are included in its frequency spectrum
Signal, avoid and sequentially input multiple single frequency signals, improve the efficiency for obtaining frequency response function.Reference signal and output
Signal and input signal make circulation cross-correlation respectively so that it is very high to obtain the precision of frequency response function, and input signal and ginseng
The cycle for examining signal is the even-multiple of power frequency period, can effectively suppress Hz noise.
Embodiment 3
The embodiments of the invention provide a kind of method for obtaining unidentified system frequency response function, this method can be by obtaining
The device of unidentified system frequency response function is taken to perform.
Wherein, unidentified system is circuit system, observation instrument system or the earth system to be detected etc..Due to system to be identified
Comprising the medium that electric conductivity is different with magnetic conductivity in system, input signal is inputted into unidentified system, due to electromagnetic induction be present,
Unidentified system can produce response to input signal.Gather unidentified system output in an asynchronous manner in embodiments of the present invention
Output signal, and obtain by method provided in an embodiment of the present invention the amplitude spectrum of unidentified system frequency response function.
Referring to Fig. 3, this method specifically includes following steps：
Step 301：It is identical with the operation of step 201, it will not be repeated here；
Step 302：The input signal is inputted to unidentified system, and gathers the defeated of unidentified system output in an asynchronous manner
Go out signal；
Wherein, asynchronous system be start to gather time point of unidentified system input signal and output signal can be different,
But sample rate is identical.
Step 303：It is identical with the operation of step 203, it will not be repeated here；
Step 304：To reference signal, input signal and output signal application Cross-Correlation Theory, unidentified system is obtained
Frequency response function.
Wherein, operation of the mode with S1-S4 in step 204 of frequency response function is specifically obtained.It is pointed out that
When gathering input signal and output signal in an asynchronous manner, the collection of output signal and the collection time delay of input signal be not true
It is fixed, therefore in asynchronous mode, the frequency response function of acquisition can not reflect unidentified system response output and input signal
Between time delay, i.e., the frequency response function obtained in asynchronous mode only includes amplitude spectrum, not comprising phase spectrum.
In embodiments of the present invention, because input signal is inverse repetition M sequence, multiple discrete frequencies are included in its frequency spectrum
Signal, avoid and sequentially input multiple single frequency signals, improve the efficiency for obtaining frequency response function.Reference signal and output
Signal and input signal make circulation cross-correlation respectively so that it is very high to obtain the precision of frequency response function, and input signal and ginseng
The cycle for examining signal is the even-multiple of power frequency period, can effectively suppress Hz noise.
Embodiment 4
Referring to Fig. 4, the embodiments of the invention provide a kind of device for obtaining unidentified system frequency response function, the device
Method for performing above-mentioned acquisition unidentified system frequency response function.The device specifically includes：
Generation module 401, input signal is generated for triggering linear feedback shift register by clock signal, according to defeated
Enter time series cycle and the sample rate triggering software program generation reference signal of signal, input signal and reference signal are inverse
Repeat M sequence and cycle phase is the same as the even-multiple for being power frequency period；
Acquisition module 402, for inputting the input signal, the output letter of collection unidentified system output to unidentified system
Number；
Acquisition module 403, for reference signal, input signal and output signal application Cross-Correlation Theory, obtaining and waiting to distinguish
The frequency response function of knowledge system.
Wherein, acquisition module 403 includes：
Cross-correlation unit is circulated, for output signal and input signal to be carried out into circulation cross-correlation with reference signal respectively,
Obtain cross-correlation time series corresponding to cross-correlation time series corresponding to output signal and input signal；
Fast Fourier transform unit, for corresponding to cross-correlation time series corresponding to output signal and input signal
Cross-correlation time series carries out FFT respectively, obtains the crosspower spectrum of output signal and the cross-power of input signal
Spectrum；
First acquisition unit, for the crosspower spectrum and the crosspower spectrum of input signal according to output signal, obtain and wait to distinguish
The frequency response function of knowledge system.
Wherein, first acquisition unit includes：
Subelement is obtained, for obtaining output spectral line peak value and its corresponding frequency from the crosspower spectrum of output signal
Value, and input spectral line peak value and its corresponding frequency values are obtained from the crosspower spectrum of input signal；
Computation subunit, for filtering out frequency values from the output spectral line peak value of acquisition and the input spectral line peak value obtained
Identical exports spectral line peak value and input spectral line peak value, and calculates frequency values identical by equation below (1) and export spectral line peak
Ratio between value and input spectral line peak value, obtains the frequency response function of unidentified system；
Wherein, in formula (1), ω is frequency values,To export spectral line peak value corresponding to frequency values ω,
To input spectral line peak value corresponding to frequency values ω, H (ω) is the frequency response function of unidentified system.
Wherein, acquisition module 402, for inputting input signal to unidentified system, and in a synchronous manner or asynchronous system
Gather the output signal of unidentified system output；
Accordingly, acquisition module 403 includes：
Second acquisition unit, for when gathering output signal in an asynchronous manner, according to the crosspower spectrum of output signal and
The crosspower spectrum of input signal, obtains the frequency response function of unidentified system, and frequency response function includes amplitude spectrum；
3rd acquiring unit, for when gathering output signal in a synchronous manner, according to the crosspower spectrum of output signal and
The crosspower spectrum of input signal, obtains the frequency response function of unidentified system, and frequency response function includes phase spectrum and amplitude
Spectrum, and pass through equation below (2), the phase spectrum of acquisition unidentified system frequency response function；
Wherein, in formula (2), mod is modulus operator,For the cross-power of output signal corresponding to frequencies omega
The phase of spectrum,For the phase of the crosspower spectrum of input signal corresponding to frequencies omega,Rung for unidentified system frequency
Answer the phase spectrum of function.
Further, the device also includes：
Computing module 404, for when unidentified system is the earth system to be detected, inverting to be carried out to frequency response function
Calculate, obtain the geoelectric cross section of the earth system to be detected, and the electricity of the earth system to be detected is calculated according to frequency response function
Magnetic explores parameter, and electromagnetic survey parameter includes step response, complex resistivity, percentage frequency effect and late period apparent resistivity.
In embodiments of the present invention, because input signal is inverse repetition M sequence, multiple discrete frequencies are included in its frequency spectrum
Signal, avoid and sequentially input multiple single frequency signals, improve the efficiency for obtaining frequency response function.Reference signal and output
Signal and input signal make circulation cross-correlation respectively so that it is very high to obtain the precision of frequency response function, and input signal and ginseng
The cycle for examining signal is the even-multiple of power frequency period, can effectively suppress Hz noise.
The device for the acquisition unidentified system frequency response function that the embodiment of the present invention is provided can be the spy in equipment
Determine hardware or the software being installed in equipment or firmware etc..It is apparent to those skilled in the art that to retouch
The convenience stated and the specific work process of succinct, described above system, device and unit, it may be referred to above method implementation
Corresponding process in example.
In several embodiments provided herein, it should be understood that disclosed apparatus and method, can be by other
Mode realize.Device embodiment described above is only schematical, for example, the division of the unit, only one
Kind of division of logic function, can there is other dividing mode when actually realizing, in another example, multiple units or component can combine or
Person is desirably integrated into another system, or some features can be ignored, or does not perform.Another, shown or discussed is mutual
Between coupling or direct-coupling or communication connection can be the INDIRECT COUPLING or logical of device or unit by some communication interfaces
Letter connection, can be electrical, mechanical or other forms.
The unit illustrated as separating component can be or may not be physically separate, show as unit
The part shown can be or may not be physical location, you can with positioned at a place, or can also be distributed to multiple
On NE.Some or all of unit therein can be selected to realize the mesh of this embodiment scheme according to the actual needs
's.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, can also
That unit is individually physically present, can also two or more units it is integrated in a unit.
If the function is realized in the form of SFU software functional unit and is used as independent production marketing or in use, can be with
It is stored in a computer read/write memory medium.Based on such understanding, technical scheme is substantially in other words
The part to be contributed to prior art or the part of the technical scheme can be embodied in the form of software product, the meter
Calculation machine software product is stored in a storage medium, including some instructions are causing a computer equipment (can be
People's computer, server, or network equipment etc.) perform all or part of step of each embodiment methods described of the present invention.
And foregoing storage medium includes：USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only Memory), arbitrary access are deposited
Reservoir (RAM, Random Access Memory), magnetic disc or CD etc. are various can be with the medium of store program codes.
The foregoing is only a specific embodiment of the invention, but protection scope of the present invention is not limited thereto, any
Those familiar with the art the invention discloses technical scope in, change or replacement can be readily occurred in, should all be contained
Cover within protection scope of the present invention.Therefore, protection scope of the present invention described should be defined by scope of the claims.
Claims (8)
- A kind of 1. method for obtaining unidentified system frequency response function, it is characterised in that methods described includes：Linear feedback shift register is triggered by clock signal and generates input signal, according to the time series of the input signal Cycle and sample rate triggering software program generation reference signal, the input signal and the reference signal are inverse repetition M sequences Arrange and cycle phase is the same as the even-multiple for being power frequency period；The input signal is inputted to unidentified system, gathers the output signal of the unidentified system output；The output signal and the input signal are subjected to circulation cross-correlation with the reference signal respectively, obtain the output Cross-correlation time series corresponding to cross-correlation time series corresponding to signal and the input signal；Cross-correlation time series corresponding to cross-correlation time series corresponding to the output signal and the input signal is distinguished FFT is carried out, obtains the crosspower spectrum of the output signal and the crosspower spectrum of the input signal；According to the crosspower spectrum of the output signal and the crosspower spectrum of the input signal, the frequency of the unidentified system is obtained Rate receptance function.
- 2. according to the method for claim 1, it is characterised in that according to the crosspower spectrum of the output signal and the input The crosspower spectrum of signal, the frequency response function of the unidentified system is obtained, including：Obtain output spectral line peak value and its corresponding frequency values from the crosspower spectrum of the output signal, and from the input Input spectral line peak value and its corresponding frequency values are obtained in the crosspower spectrum of signal；The output of frequency values identical is filtered out from the output spectral line peak value of the acquisition and the input spectral line peak value of the acquisition Spectral line peak value and input spectral line peak value, and by equation below (1) calculate the frequency values identical export spectral line peak value with it is defeated Enter the ratio between spectral line peak value, obtain the frequency response function of the unidentified system；<mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>P</mi> <msub> <mi>C</mi> <mrow> <mi>y</mi> <mi>s</mi> </mrow> </msub> </msub> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>P</mi> <msub> <mi>C</mi> <mrow> <mi>u</mi> <mi>s</mi> </mrow> </msub> </msub> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>...</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>Wherein, in formula (1), ω is frequency values,To export spectral line peak value corresponding to frequency values ω,For frequency Input spectral line peak value corresponding to rate value ω, H (ω) are the frequency response function of the unidentified system.
- 3. according to the method for claim 1, it is characterised in that it is described to input the input signal to unidentified system, adopt Collect the output signal of the unidentified system output, including：The input signal is inputted to the unidentified system, and in a synchronous manner or asynchronous system gathers the unidentified system The output signal of output；Accordingly, when gathering the output signal with the asynchronous system, the frequency response function includes amplitude spectrum；When with When the method for synchronization gathers the output signal, the frequency response function includes phase spectrum and amplitude spectrum, and by such as Lower formula (2), obtain the phase spectrum of the unidentified system frequency response function；Wherein, in formula (2), mod is modulus operator,For the cross-power of the output signal corresponding to frequencies omega The phase of spectrum,For the phase of the crosspower spectrum of the input signal corresponding to frequencies omega,For the system to be identified The phase spectrum of system frequency response function.
- 4. according to the method described in any one of claim 1-3 claim, it is characterised in that described to obtain the system to be identified After the frequency response function of system, in addition to：When the unidentified system is the earth system to be detected, Inversion Calculation is carried out to the frequency response function, obtains institute The geoelectric cross section of the earth system to be detected is stated, and the electricity of the earth system to be detected is calculated according to the frequency response function Magnetic explores parameter, and the electromagnetic survey parameter includes step response, complex resistivity, percentage frequency effect and late period apparent resistance Rate.
- 5. a kind of device for obtaining unidentified system frequency response function, it is characterised in that described device includes：Generation module, input signal is generated for triggering linear feedback shift register by clock signal, according to the input The time series cycle of signal and sample rate triggering software program generation reference signal, the input signal and the reference signal It is inverse to repeat M sequence and cycle phase is the same as the even-multiple for being power frequency period；Acquisition module, for inputting the input signal to unidentified system, gather the output letter of the unidentified system output Number；Acquisition module, for carrying out circulating mutually the output signal and the input signal with the reference signal respectively Close, obtain cross-correlation time series corresponding to cross-correlation time series corresponding to the output signal and the input signal；It is right Cross-correlation time series corresponding to cross-correlation time series corresponding to the output signal and the input signal is carried out soon respectively Fast Fourier transform, obtain the crosspower spectrum of the output signal and the crosspower spectrum of the input signal；According to the output The crosspower spectrum of the crosspower spectrum of signal and the input signal, obtain the frequency response function of the unidentified system.
- 6. device according to claim 5, it is characterised in that the acquisition module includes：Subelement is obtained, for obtaining output spectral line peak value and its corresponding frequency from the crosspower spectrum of the output signal Value, and input spectral line peak value and its corresponding frequency values are obtained from the crosspower spectrum of the input signal；Computation subunit, for filtering out frequency from the input spectral line peak value of the output spectral line peak value of the acquisition and the acquisition Rate value identical exports spectral line peak value and input spectral line peak value, and defeated by equation below (1) the calculating frequency values identical Go out spectral line peak value and the ratio inputted between spectral line peak value, obtain the frequency response function of the unidentified system；<mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msub> <mi>P</mi> <msub> <mi>C</mi> <mrow> <mi>y</mi> <mi>s</mi> </mrow> </msub> </msub> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>P</mi> <msub> <mi>C</mi> <mrow> <mi>u</mi> <mi>s</mi> </mrow> </msub> </msub> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>...</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>Wherein, in formula (1), ω is frequency values,To export spectral line peak value corresponding to frequency values ω,For frequency Input spectral line peak value corresponding to rate value ω, H (ω) are the frequency response function of the unidentified system.
- 7. device according to claim 5, it is characterised in that the acquisition module, for defeated to the unidentified system Enter the input signal, and in a synchronous manner or asynchronous system gathers the output signal of the unidentified system output；Accordingly, the acquisition module includes：Second acquisition unit, for when with the asynchronous system gather the output signal when, according to the mutual of the output signal The crosspower spectrum of power spectrum and the input signal, obtain the frequency response function of the unidentified system, the frequency response Function includes amplitude spectrum；3rd acquiring unit, for when with the method for synchronization gather the output signal when, according to the mutual of the output signal The crosspower spectrum of power spectrum and the input signal, obtain the frequency response function of the unidentified system, the frequency response Function includes phase spectrum and amplitude spectrum, and by equation below (2), obtains the phase of the unidentified system frequency response function Position spectrum；Wherein, in formula (2), mod is modulus operator,For the crosspower spectrum of the output signal corresponding to frequencies omega Phase,For the phase of the crosspower spectrum of the input signal corresponding to frequencies omega,For the unidentified system The phase spectrum of frequency response function.
- 8. according to the device described in any one of claim 5-7 claim, it is characterised in that described device also includes：Computing module, for when the unidentified system is the earth system to be detected, being carried out to the frequency response function anti- Calculating is drilled, obtains the geoelectric cross section of the earth system to be detected, and wait to visit according to calculating the frequency response function The electromagnetic survey parameter of the earth system is surveyed, the electromagnetic survey parameter includes step response, complex resistivity, percentage frequency effect With late period apparent resistivity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510253306.2A CN104794313B (en) | 2015-05-18 | 2015-05-18 | A kind of method and device for obtaining unidentified system frequency response function |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510253306.2A CN104794313B (en) | 2015-05-18 | 2015-05-18 | A kind of method and device for obtaining unidentified system frequency response function |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104794313A CN104794313A (en) | 2015-07-22 |
| CN104794313B true CN104794313B (en) | 2018-01-23 |
Family
ID=53559104
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510253306.2A Expired - Fee Related CN104794313B (en) | 2015-05-18 | 2015-05-18 | A kind of method and device for obtaining unidentified system frequency response function |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104794313B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105607131B (en) * | 2016-01-11 | 2018-03-23 | 甘肃省有色地质调查院 | Coding source electromagnetic depth method obtains the method and device of the earth APPARENT RESISTIVITY spectrum |
| CN106054269B (en) * | 2016-07-14 | 2019-01-08 | 中国地质大学（北京） | A kind of Correlation Identification system and method for frequency segmentation |
| CN108801296B (en) * | 2018-06-13 | 2020-06-05 | 安徽大学 | Sensor frequency response function calculation method based on error model iterative compensation |
| CN110187394B (en) * | 2019-06-20 | 2020-11-27 | 甘肃省地震局（中国地震局兰州地震研究所） | Method and device for acquiring formation resistivity anisotropy by double-field source electromagnetic sounding method |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102426393A (en) * | 2011-11-16 | 2012-04-25 | 中国地质大学（北京） | Electric prospecting method and device |
-
2015
- 2015-05-18 CN CN201510253306.2A patent/CN104794313B/en not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102426393A (en) * | 2011-11-16 | 2012-04-25 | 中国地质大学（北京） | Electric prospecting method and device |
Non-Patent Citations (4)
| Title |
|---|
| 《Circular correlation based identification of switching power converter with uncertainty analysis using fuzzy density approach》;Tomi Roinila等;《Simulation Modelling Practice & Theory》;20090731;第19卷(第6期);第1043-1058页 * |
| 《基于伪随机***辨识的电磁法仪器标定》;刘立超等;《仪器仪表学报》;20140831;第35卷(第8期);第1721-1728页 * |
| 《基于相关辨识的逆重复m序列伪随机电磁法》;汤井田等;《地球物理学报》;20080731;第51卷(第4期);第1226-1233页 * |
| 《编码电磁测深》;罗维斌等;《地球物理学报》;20120131;第55卷(第1期);第341页摘要，第343-344页第2节，图2-3 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104794313A (en) | 2015-07-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104794313B (en) | A kind of method and device for obtaining unidentified system frequency response function | |
| CN105607131B (en) | Coding source electromagnetic depth method obtains the method and device of the earth APPARENT RESISTIVITY spectrum | |
| CN110529087B (en) | Method and device for evaluating hydraulic fracturing effect of stratum | |
| CN104391330B (en) | Relevant identification technology based electrical prospecting system | |
| CN103823244B (en) | Magnetic resonance three-component noise removing device and noise removing method | |
| CN108267784A (en) | A kind of seismic signal random noise compression process method | |
| CN104597499B (en) | Adjacent shot interference suppressing method and device for independent simultaneous sweeping seismic data of seismic vibrators | |
| CN105640545A (en) | Fetal electrocardiosignal extraction method and device | |
| CN102053278A (en) | Electrical prospecting method and measuring device | |
| CN203054093U (en) | Apparatus for extracting DC (direct current) signal from mixed signal by using synchronous sampling | |
| CN102809687B (en) | Digital measurement method for alternating-current frequency | |
| CN112285786B (en) | Method for improving exploration depth based on ultra-high density electrical equipment | |
| CN102520451B (en) | High-precision double-way ground voice frequency electromagnetic physical geography prospecting instrument | |
| CN103713310A (en) | Nuclear pulse signal high speed digital separation method and system and digital spectrometer | |
| CN103944655A (en) | Noise identification method in power line carrier communication fault detection system | |
| CN105467460A (en) | Method and device for electromagnetic prospecting | |
| CN102916683A (en) | Parameter-adjustable nuclear pulse simulation method | |
| CN102520449A (en) | Receiver device of high density electric method device | |
| CN109001828A (en) | A kind of recognition methods of buried target body and device, storage medium | |
| CN106772193B (en) | Measuring method using current transformer frequency characteristic measuring device | |
| CN109948223A (en) | A kind of impulse amplitude acquisition methods based on Lagrange's interpolation | |
| CN106125142A (en) | A kind of Underground medium system based on Correlation Identification frequency domain method and method thereof | |
| CN104459485B (en) | The treating method and apparatus of shelf depreciation ultrasonic signal | |
| CN104502986B (en) | Physical prospecting induced electrical sounding data chromatography processing method | |
| CN112764109B (en) | Separation and extraction method and device for dipole transverse wave reflected waves |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate 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: 20180123 Termination date: 20200518 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |