CN112924964A - Data fusion method and device for stepping frequency radar signal and storage medium - Google Patents
Data fusion method and device for stepping frequency radar signal and storage medium Download PDFInfo
- Publication number
- CN112924964A CN112924964A CN202110286813.1A CN202110286813A CN112924964A CN 112924964 A CN112924964 A CN 112924964A CN 202110286813 A CN202110286813 A CN 202110286813A CN 112924964 A CN112924964 A CN 112924964A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- component
- signals
- data fusion
- 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.)
- Pending
Links
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 17
- 230000005540 biological transmission Effects 0.000 claims abstract description 31
- 230000004927 fusion Effects 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 15
- 238000001514 detection method Methods 0.000 claims abstract description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000000149 penetrating effect Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000002592 echocardiography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a data fusion method, a data fusion device and a storage medium for a stepping frequency radar signal, wherein the method comprises the following steps of: acquiring sinusoidal continuous wave transmitting signals with different frequencies; acquiring a reflection signal corresponding to the emission signal; carrying out orthogonal decomposition on each frequency reflection signal to obtain corresponding in-phase component InAnd quadrature component Qn(ii) a For each frequency fnCorresponding to InAnd QnRespectively with a cut-off frequency of fnIs processed to obtain InAnd QnDifference frequency component I 'in the result'nAnd Q'n(ii) a Passing the reflected signal of each frequency through I 'in the frequency domain'n+jQ′nIs expressed in terms of form; and performing signal fusion on the conjugate symmetric sequence composed of the reflection signals of all the frequencies in the frequency domain by inverse Fourier transform. The scheme makes the frequency of the signal obtained by fusionThe method is flexible and controllable, and the synthesized radar image has higher signal-to-noise ratio compared with the conventional pulse signal; since the continuous wave signal is transmitted, the transmission power of the signal is increased in use, and the detection depth is larger.
Description
Technical Field
The invention relates to the field of exploration geophysical data processing, in particular to a multi-frequency data fusion method and device for a stepping frequency radar signal and a storage medium.
Background
The ground penetrating radar is used as a high-efficiency nondestructive detection technology, and is widely applied to the aspects of engineering investigation, concrete member defect, pavement detection, roadbed disease exploration and the like, and a good detection effect is obtained. Generally, if a high-resolution detection result is required, a radar system with a high center frequency is required, but the detection depth is shallow. If a larger exploration depth is required, a radar system with a lower center frequency is selected, but the resolution is lower.
However, the performance indexes of the ground penetrating radar, namely the detection resolution and the exploration depth, cannot be considered together with the requirement of practical application. The traditional pulse system radar has extremely short pulse transmitting time, small average power of transmitted signals and low signal to noise ratio of obtained echoes, and influences the exploration depth and the resolution capability of the ground penetrating radar. The step frequency ground penetrating radar gradually becomes a research hotspot in the field in recent years as one of frequency domain ground penetrating radars, and the advantages of the step frequency ground penetrating radar in detection depth and resolution compared with an impulse system ground penetrating radar are also reflected in research application. However, for data fusion of step frequency radar signals containing multiple frequencies, no mature method exists at present, so that the interpretation of radar signals of multiple frequencies received by the step frequency radar is difficult at present.
Disclosure of Invention
The invention provides a data fusion method, a data fusion device and a storage medium for a stepping frequency radar signal, which are used for solving the problem that the radar signal with various frequencies received by a stepping frequency radar is difficult to interpret at present.
In a first aspect, a data fusion method for stepped frequency radar signals is provided, including:
sinusoidal continuous wave transmit signals of different frequencies are acquired and denoted as Tx (f)1)，Tx(f2)……Tx(fn) Wherein f isnAn nth frequency representing a stepped frequency;
reflected signals corresponding to sinusoidal continuous wave transmission signals of different frequencies are obtained and are denoted as Rx (f)1)，Rx(f2)……Rx(fn)；
Orthogonal decomposition is carried out on the reflected signal of each frequency respectively to obtain corresponding in-phase component InAnd quadrature component Qn；
For each frequency fnThe corresponding in-phase component and orthogonal component respectively adopt a cut-off frequency fnIs processed to obtain the corresponding difference frequency component I 'in the results of the in-phase component and the quadrature component'nAnd Q'n；
Passing the reflected signal of each frequency through I 'in the frequency domain'n+jQ′nIs expressed in terms of form;
and performing signal fusion on the conjugate symmetric sequence composed of the reflection signals of all the frequencies in the frequency domain by inverse Fourier transform.
Further, the transmission signal is obtained by the following method:
selecting the bandwidth range of the stepping frequency continuous wave signal;
sinusoidal continuous wave signals of different frequencies are generated as transmission signals based on the bandwidth range.
Further, the performing quadrature decomposition on the reflected signal of each frequency to obtain a corresponding in-phase component and a corresponding quadrature component specifically includes:
reflection signal Rx (f) for each frequencyn) With the same frequency of the transmitted signal Tx (f)n) Multiplying to obtain an in-phase component;
reflection signal Rx (f) for each frequencyn) Multiplying the quadrature component by the cosine signal with the same frequency; wherein the cosine signal is the sum transmission signal Tx (f)n) Same frequency but phase difference
Of the signal of (1).
Further, the transmission signal is represented as:
wherein t is the duration, and k is a constant less than 1; through kfnt to control the amplitude of the signal at the previous stage of transmission, when kfnt<When 1, the amplitude of the transmitted signal is increased at any time for a long time t, and the smaller k is, the slower the amplitude is increased; when kfnt>At 1, the transmit signal amplitude takes 1.
Further, the reflection signals are obtained by forward modeling based on the constructed two-dimensional geological model and by adopting a forward modeling program based on a finite difference method, and the reflection signals are expressed as follows:
wherein A isnIs the echo amplitude; f. ofn＝f0+(n-1)Δf，f0Step starting frequency, n represents the nth frequency point of the step frequency, and delta f is the frequency step length;
the phase change of the reflected signal relative to the transmitted signal is shown as v, the propagation speed of the transmitted signal or the reflected signal in the medium to be detected is shown as d, and the position of the detection point and the target to be detected is shown as d.
Further, the in-phase component InAnd quadrature component QnAre respectively obtained by the following formulas:
for each frequency fnCorresponding in-phase component InAnd quadrature component QnSeparately miningWith a cut-off frequency fnIs processed to obtain the corresponding difference frequency component I 'in the results of the in-phase component and the quadrature component'nAnd Q'nExpressed as:
then, the reflected signal of each frequency is represented in the frequency domain as:
Rx(n)＝I′n+jQ′n。
in a second aspect, there is provided a data fusion apparatus for stepped frequency radar signals, comprising:
a transmission signal acquisition module for acquiring sinusoidal continuous wave transmission signals of different frequencies, denoted as Tx (f)1)，Tx(f2)……Tx(fn) Wherein f isnAn nth frequency representing a stepped frequency;
a reflected signal acquiring module, configured to acquire reflected signals corresponding to sinusoidal continuous wave transmitting signals with different frequencies, and denoted as Rx (f1), Rx (f2) … … Rx (fn);
the orthogonal decomposition module is used for respectively carrying out orthogonal decomposition on the reflected signal of each frequency to obtain a corresponding in-phase component InAnd quadrature component Qn；
A filtering module for each frequency fnThe corresponding in-phase component and orthogonal component respectively adopt a cut-off frequency fnIs processed to obtain the corresponding difference frequency component I 'in the results of the in-phase component and the quadrature component'nAnd Q'n；
A multi-frequency data fusion module for passing the reflection signal of each frequency through I 'in the frequency domain'n+jQ′nIs expressed in terms of form; the conjugate symmetrical sequence composed of the reflection signals of all frequencies in the frequency domain is transformed into the time domain by inverse Fourier transformAnd (4) signal fusion.
In a second aspect, a computer-readable storage medium is provided, having stored thereon a computer program adapted to be loaded by a processor and to execute the method for data fusion of stepped frequency radar signals as described above.
Advantageous effects
The invention provides a data fusion method, a data fusion device and a storage medium for a stepping frequency radar signal, which are used for obtaining a corresponding in-phase component I by performing orthogonal decomposition on a reflection signal of each frequencynAnd quadrature component QnAnd for it respectively a cut-off frequency fnIs processed to obtain the corresponding difference frequency component I 'in the results of the in-phase component and the quadrature component'nAnd Q'n(ii) a The reflected signal of each frequency is then passed through I 'in the frequency domain'n+jQ′nIs expressed in terms of form; and finally, performing signal fusion on a conjugate symmetric sequence formed by the reflection signals of all frequencies in a frequency domain through inverse Fourier transform to a time domain, and obtaining a time domain radar image. The proposal ensures that the frequency components of the signals obtained by fusion are flexible and controllable, the anti-noise performance is good, and the synthesized radar image has higher signal-to-noise ratio compared with the conventional pulse signals under the noise condition. And because the continuous wave signal is transmitted, the transmission power of the signal is increased in practical use, and the detection depth is also larger. The scheme realizes the consideration of two performance indexes of detectable resolution and exploration depth.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of a data fusion method for stepped frequency radar signals according to an embodiment of the present invention;
FIG. 2 is a two-dimensional geological model map for forward modeling provided by an embodiment of the present invention;
FIG. 3 is a 100MHz transmit signal waveform provided by embodiments of the present invention;
FIG. 4 is a 100MHz first pass reflected signal waveform provided by embodiments of the present invention;
FIG. 5 is a time domain signal fused with 201 frequencies at a first pass location according to an embodiment of the present invention;
FIG. 6 is a Bscan image generated by step frequency fusion data under noiseless conditions according to an embodiment of the present invention;
FIG. 7 is a Bscan image generated by step frequency fusion data under noisy conditions according to an embodiment of the present invention;
FIG. 8 is a Bscan image of a Ricker wavelet signal under noiseless conditions according to the present invention;
fig. 9 is a Bscan image obtained by forward modeling of the Ricker wavelet signal under the noise condition provided by the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be described in detail below. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the examples given herein without any inventive step, are within the scope of the present invention.
Example 1
As shown in fig. 1, the present embodiment provides a data fusion method for a stepped frequency radar signal, including:
(1) sinusoidal continuous wave transmit signals of different frequencies are acquired and denoted as Tx (f)1)，Tx(f2)……Tx(fn) Wherein f isnAn nth frequency representing a stepped frequency;
(2) reflected signals corresponding to sinusoidal continuous wave transmission signals of different frequencies are obtained and are denoted as Rx (f)1)，Rx(f2)……Rx(fn)；
(3) Are respectively pairedThe reflected signal of each frequency is subjected to orthogonal decomposition to obtain a corresponding in-phase component InAnd quadrature component Qn；
(4) For each frequency fnThe corresponding in-phase component and orthogonal component respectively adopt a cut-off frequency fnIs processed to obtain the corresponding difference frequency component I 'in the results of the in-phase component and the quadrature component'nAnd Q'n；
(5) Passing the reflected signal of each frequency through I 'in the frequency domain'n+jQ′nIs expressed in terms of form;
(6) and performing signal fusion on the conjugate symmetric sequence composed of the reflection signals of all the frequencies in the frequency domain by inverse Fourier transform.
Specifically, the transmission signal is obtained by the following method:
selecting the bandwidth range of the stepping frequency continuous wave signal;
sinusoidal continuous wave signals of different frequencies are generated as transmission signals based on the bandwidth range.
The quadrature decomposition of the reflected signal of each frequency to obtain a corresponding in-phase component and a corresponding quadrature component specifically includes:
reflection signal Rx (f) for each frequencyn) With the same frequency of the transmitted signal Tx (f)n) Multiplying to obtain an in-phase component;
reflection signal Rx (f) for each frequencyn) Multiplying the quadrature component by the cosine signal with the same frequency; wherein the cosine signal is the sum transmission signal Tx (f)n) Same frequency but phase difference
Of the signal of (1).
To further understand the technical solution of the present invention, a forward simulation example for a step frequency ground penetrating radar is described below.
Firstly, constructing a two-dimensional geological model for forward modeling, wherein the set model is a layered model, the upper layer is a medium 1, the thickness of the layer is 2m, the relative dielectric constant is 7, and the conductivity is 0.005S/m as shown in figure 2; the lower layer is medium 2, the thickness of the layer is 2m, the relative dielectric constant is 12, and the conductivity is 0.002S/m; a circular cavity with the radius of 0.3m is arranged in the middle of the model, the relative dielectric constant of the cavity is 1, and the conductivity is 0.
A ground penetrating radar forward modeling program based on a finite difference algorithm is adopted to carry out grid subdivision on a model, the selected stepping frequency range is 100 MHz-300 MHz, the step length is 1MHz, 201 frequency points are shared, the size of the model is 10m multiplied by 4.2m, the position of a transmitting point is 0.1m x, 4m y, the position of a receiving point is 0.2m x, 4m y, and the moving step length of the transmitting point and the receiving point is 0.05 m. 194 channels of data are collected, and the number of sampling points of each channel of data is 8192 (a time window is 96 ns).
The form of the transmission signal is a continuous sine wave, taking 100MHz as an example, the waveform of the transmission signal is shown in fig. 3, in order to avoid introducing high-frequency noise in the forward modeling result, the amplitude of the transmission signal needs to be controlled to slowly increase from zero to 1, and the transmission signal is expressed as:
where t is the duration, k is a constant of less than 1 (0.25 in this example), fnFor the nth frequency of the stepped frequency, using kfnt to control the amplitude of the signal at the previous stage of transmission, when kfnt<When 1, the amplitude of the transmitted signal is slowly increased along with the time t, and the smaller k is, the slower the amplitude is increased; when kfnt>At 1, the transmit signal amplitude takes 1.
With v representing the velocity of the signal propagating in a homogeneous medium, a target reflection signal at distance d can be represented as:
wherein A isnAs echo amplitude, f0Is the step starting frequency, n represents the nth frequency point of the step frequency, deltaf is the frequency step length,
is the signal two-way travel time.
The above-mentioned reflected signal can be simply expressed as:
wherein f isn＝f0+ (n-1) Δ f, which is the nth frequency of the step frequency;
is the phase change of the reflected signal relative to the transmitted signal.
The reflected signal waveform corresponding to the transmitted signal at 100MHz at the first pass position is shown in figure 4.
After obtaining the reflected signal of each frequency, it needs to be subjected to quadrature decomposition, wherein the in-phase component InAnd quadrature component QnAre respectively obtained by the following formulas:
that is, the obtained reflection signal of each frequency is multiplied by the transmission signal of the same frequency to obtain the in-phase component, and the reflection signal of each frequency is multiplied by the cosine signal of the same frequency (the cosine signal is the same frequency as the transmission signal but has a phase difference with the transmission signal)
Cosine signal of (c) to obtain the quadrature component. Then for the in-phase component I of each frequencynAnd quadrature component QnRespectively with a cut-off frequency of fnIs processed by the low-pass filter to obtain a difference frequency component I'nAnd Q'n：
Then, the reflected signal of each frequency is represented in the frequency domain as:
Rx(n)＝I′n+jQ′n
where j is the complex imaginary sign.
For the sequence Rx (n) generated by 201 frequency points, a conjugate symmetric sequence is formed in a frequency domain, the conjugate symmetric sequence is converted into a time domain through inverse Fourier transform for data fusion, and the synthesized first time domain data is shown in figure 5. From the resultant data it can be seen that the horizon boundary reflection signal occurs at a position of approximately 36ns, corresponding to a depth of 2.0394m, which is consistent with the model setting of the boundary at a depth of 2 m. A Bscan image synthesized from 194 scan lane position fused data is shown in FIG. 6. For comparison with the forward modeling result of pulse signals commonly used in ground penetrating radar, a Ricker wavelet with a center frequency of 200MHz is used as a source signal for forward modeling, and 194 channels of data are collected to synthesize a Bscan image as shown in FIG. 8.
In order to compare the stepping frequency data fusion result with the pulse signal forward result under the noise condition, noise with the same degree is added in the forward processes of the stepping frequency data fusion result and the pulse signal forward result respectively. For the Ricker wavelet signal, noise is added to each reflection signal, and then a Bscan image is synthesized (see FIG. 9); after the same level of noise is added to the continuous wave reflection signal of each frequency point of the SFCW (stepped frequency continuous wave) signal, the SFCW signal is processed according to the flow (3) → (4) → (5) → (6) in fig. 1, and then converted into single channel data synthesized in the time domain, and the 194 channel data obtained by fusion is synthesized into a Bscan image (see fig. 7). The comparison of the Bscan results obtained by the two methods before and after the noise is added shows that the signal-to-noise ratio of the result corresponding to the Ricker wavelet signal is obviously reduced after the noise is added, and radar echoes reflected by a horizon interface and an abnormal body are difficult to identify (as shown in FIG. 9); and the multi-frequency data fusion result of the stepping frequency signal still keeps higher signal-to-noise ratio under the condition of adding noise, and the echo characteristics are clear (as shown in figure 7). The above results show that the proposed multi-frequency data fusion method for stepped frequency radar signals has better noise immunity.
Example 2
The embodiment provides a data fusion device for a step frequency radar signal, which comprises:
a transmission signal acquisition module for acquiring sinusoidal continuous wave transmission signals of different frequencies, denoted as Tx (f)1)，Tx(f2)……Tx(fn) Wherein f isnAn nth frequency representing a stepped frequency;
a reflected signal acquiring module, configured to acquire reflected signals corresponding to sinusoidal continuous wave transmitting signals with different frequencies, and denoted as Rx (f1), Rx (f2) … … Rx (fn);
the orthogonal decomposition module is used for respectively carrying out orthogonal decomposition on the reflected signal of each frequency to obtain a corresponding in-phase component InAnd quadrature component Qn；
A filtering module for each frequency fnThe corresponding in-phase component and orthogonal component respectively adopt a cut-off frequency fnIs processed to obtain the corresponding difference frequency component I 'in the results of the in-phase component and the quadrature component'nAnd Q'n；
A multi-frequency data fusion module for passing the reflection signal of each frequency through I 'in the frequency domain'n+jQ′nIs expressed in terms of form; and performing signal fusion on the conjugate symmetric sequence composed of the reflection signals of all the frequencies in the frequency domain by inverse Fourier transform.
Example 3
The present embodiment provides a computer-readable storage medium storing a computer program adapted to be loaded by a processor and to perform the data fusion method for stepped frequency radar signals as described above.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It is understood that the same or similar parts in the above embodiments may be mutually referred to, and the same or similar parts in other embodiments may be referred to for the content which is not described in detail in some embodiments.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and alternate implementations are included within the scope of the preferred embodiment of the present invention in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the present invention.
Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention, and that variations, modifications, substitutions and alterations can be made to the above embodiments by those of ordinary skill in the art within the scope of the present invention.
Claims (8)
1. A data fusion method for stepped frequency radar signals, comprising:
sinusoidal continuous wave transmit signals of different frequencies are acquired and denoted as Tx (f)1)，Tx(f2)……Tx(fn) Wherein f isnAn nth frequency representing a stepped frequency;
acquiring reflection signals corresponding to sinusoidal continuous wave emission signals with different frequencies, and expressing the reflection signals as Rx (f1), Rx (f2) … … Rx (fn);
orthogonal decomposition is carried out on the reflected signal of each frequency respectively to obtain corresponding in-phase component InAnd quadrature component Qn；
For each frequency fnThe corresponding in-phase component and orthogonal component respectively adopt a cut-off frequency fnIs processed to obtain the corresponding difference frequency component I 'in the results of the in-phase component and the quadrature component'nAnd Q'n；
Passing the reflected signal of each frequency through I 'in the frequency domain'n+jQ′nIs expressed in terms of form;
and performing signal fusion on the conjugate symmetric sequence composed of the reflection signals of all the frequencies in the frequency domain by inverse Fourier transform.
2. The data fusion method for stepped frequency radar signals according to claim 1, wherein the transmission signal is obtained by:
selecting the bandwidth range of the stepping frequency continuous wave signal;
sinusoidal continuous wave signals of different frequencies are generated as transmission signals based on the bandwidth range.
3. The data fusion method for stepped frequency radar signals according to claim 1, wherein the performing quadrature decomposition on the reflected signal of each frequency to obtain corresponding in-phase component and quadrature component specifically comprises:
reflection signal Rx (f) for each frequencyn) With the same frequency of the transmitted signal Tx (f)n) Multiplying to obtain an in-phase component;
reflection signal Rx (f) for each frequencyn) Multiplying the quadrature component by the cosine signal with the same frequency; wherein the cosine signal is the sum transmission signal Tx (f)n) Same frequency but phase difference
Of the signal of (1).
4. The data fusion method for stepped-frequency radar signals according to claim 1, wherein the transmission signal is represented as:
wherein t is the duration, and k is a constant less than 1; through kfnt to control the amplitude of the signal at the previous stage of transmission, when kfnt<When 1, the amplitude of the transmitted signal is increased at any time for a long time t, and the smaller k is, the slower the amplitude is increased; when kfnt>At 1, the transmit signal amplitude takes 1.
5. The data fusion method for the stepped-frequency radar signal according to claim 4, wherein the reflection signal is obtained by forward modeling based on a two-dimensional geological model and using a forward modeling procedure based on a finite difference method, and the reflection signal is represented as follows:
wherein A isnIs the echo amplitude; f. ofn＝f0+(n-1)△f，f0Step starting frequency, n represents the nth frequency point of the step frequency, and delta f is the frequency step length;
the phase change of the reflected signal relative to the transmitted signal is shown as v, the propagation speed of the transmitted signal or the reflected signal in the medium to be detected is shown as d, and the position of the detection point and the target to be detected is shown as d.
6. The data fusion method for stepped-frequency radar signals according to claim 5, characterized in that the in-phase component InAnd quadrature component QnAre respectively obtained by the following formulas:
for each frequency fnCorresponding in-phase component InAnd quadrature component QnRespectively with a cut-off frequency of fnIs processed to obtain the corresponding difference frequency component I 'in the results of the in-phase component and the quadrature component'nAnd Q'nExpressed as:
then, the reflected signal of each frequency is represented in the frequency domain as:
Rx(n)＝I′n+jQ′n。
7. a data fusion apparatus for stepped frequency radar signals, comprising:
a transmission signal acquisition module for acquiring sinusoidal continuous wave transmission signals of different frequencies, denoted as Tx (f)1)，Tx(f2)……Tx(fn) Wherein f isnAn nth frequency representing a stepped frequency;
a reflected signal acquiring module, configured to acquire reflected signals corresponding to sinusoidal continuous wave transmitting signals with different frequencies, and denoted as Rx (f1), Rx (f2) … … Rx (fn);
the orthogonal decomposition module is used for respectively carrying out orthogonal decomposition on the reflected signal of each frequency to obtain a corresponding in-phase component InAnd quadrature component Qn；
A filtering module for each frequency fnThe corresponding in-phase component and orthogonal component respectively adopt a cut-off frequency fnIs processed to obtain the corresponding difference frequency component I 'in the results of the in-phase component and the quadrature component'nAnd Q'n；
A multi-frequency data fusion module for passing the reflection signal of each frequency through I 'in the frequency domain'n+jQ′nIs expressed in terms of form; and performing signal fusion on the conjugate symmetric sequence composed of the reflection signals of all the frequencies in the frequency domain by inverse Fourier transform.
8. A computer-readable storage medium, in which a computer program is stored which is adapted to be loaded by a processor and to carry out a method for data fusion of stepped frequency radar signals according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110286813.1A CN112924964A (en) | 2021-03-17 | 2021-03-17 | Data fusion method and device for stepping frequency radar signal and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110286813.1A CN112924964A (en) | 2021-03-17 | 2021-03-17 | Data fusion method and device for stepping frequency radar signal and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112924964A true CN112924964A (en) | 2021-06-08 |
Family
ID=76175874
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110286813.1A Pending CN112924964A (en) | 2021-03-17 | 2021-03-17 | Data fusion method and device for stepping frequency radar signal and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112924964A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113390371A (en) * | 2021-06-02 | 2021-09-14 | 郑州大学 | Microwave sensing-based target multipoint deformation rapid detection and analysis method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5499029A (en) * | 1992-07-14 | 1996-03-12 | Eg&G Energy Measurements, Inc. | Wide band stepped frequency ground penetrating radar |
| CN104181521A (en) * | 2014-06-24 | 2014-12-03 | 合肥工业大学 | High range resolution radar capable of emitting multifrequency carrier wave |
| CN109738878A (en) * | 2019-02-27 | 2019-05-10 | 电子科技大学 | Radar range profile's recognition methods based on compressed sensing and frequency stepping waveform |
| CN109782233A (en) * | 2019-01-25 | 2019-05-21 | 北京电子工程总体研究所 | A kind of radar working method and system based on Fourier transformation |
| CN111273250A (en) * | 2020-03-04 | 2020-06-12 | 中南大学 | Nonlinear frequency stepping method and system for stepping frequency radar |
| CN111427018A (en) * | 2020-04-22 | 2020-07-17 | 安徽华可智能科技有限公司 | Radar interference equipment interference effect evaluation method |
-
2021
- 2021-03-17 CN CN202110286813.1A patent/CN112924964A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5499029A (en) * | 1992-07-14 | 1996-03-12 | Eg&G Energy Measurements, Inc. | Wide band stepped frequency ground penetrating radar |
| CN104181521A (en) * | 2014-06-24 | 2014-12-03 | 合肥工业大学 | High range resolution radar capable of emitting multifrequency carrier wave |
| CN109782233A (en) * | 2019-01-25 | 2019-05-21 | 北京电子工程总体研究所 | A kind of radar working method and system based on Fourier transformation |
| CN109738878A (en) * | 2019-02-27 | 2019-05-10 | 电子科技大学 | Radar range profile's recognition methods based on compressed sensing and frequency stepping waveform |
| CN111273250A (en) * | 2020-03-04 | 2020-06-12 | 中南大学 | Nonlinear frequency stepping method and system for stepping frequency radar |
| CN111427018A (en) * | 2020-04-22 | 2020-07-17 | 安徽华可智能科技有限公司 | Radar interference equipment interference effect evaluation method |
Non-Patent Citations (2)
| Title |
|---|
| 罗兴华: "超宽带步进频雷达前端的实现及其相关问题研究", 《中国优秀硕士学位论文全文数据库信息科技辑》, 15 December 2014 (2014-12-15), pages 136 - 140 * |
| 陈波: "步进频率连续波探地雷达数字中频接收机设计与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》, 15 December 2011 (2011-12-15), pages 140 - 853 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113390371A (en) * | 2021-06-02 | 2021-09-14 | 郑州大学 | Microwave sensing-based target multipoint deformation rapid detection and analysis method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108549067B (en) | Coherent accumulation detection method applied to third-order maneuvering target | |
| Zheng et al. | Radar high-speed target detection based on the scaled inverse Fourier transform | |
| CN106371097A (en) | Radar system | |
| US5760732A (en) | Method and apparatus for enhanced resolution of range estimates in echo location for detection and imaging systems | |
| CN109738894B (en) | High squint multi-angle imaging method for large-field-of-view synthetic aperture radar | |
| Zheng et al. | ISAR imaging of targets with complex motions based on a noise-resistant parameter estimation algorithm without nonuniform axis | |
| WO2009122193A1 (en) | Multi-range object location estimation | |
| US9024815B2 (en) | Direct-to-digital software-defined radar | |
| CN113805166B (en) | Target tracking and ranging method and system of radar level gauge | |
| Xenaki et al. | Compressive synthetic aperture sonar imaging with distributed optimization | |
| JP6179940B2 (en) | Doppler imaging signal transmitter, Doppler imaging signal receiver, Doppler imaging system and method | |
| CN104730498A (en) | Target detection method based on Keystone and weighting rotating FFT | |
| CN108169746B (en) | Linear frequency modulation pulse semi-active radar seeker signal processing method | |
| CN109633622A (en) | Across the unit accumulation method of radar maneuvering target based on time reversal second order Keystone transformation | |
| KR20200075858A (en) | System and method for composite aperture radar image formation | |
| CN115508820A (en) | Target detection method of linear frequency modulation pulse radar | |
| CN110879391B (en) | Radar image data set manufacturing method based on electromagnetic simulation and missile-borne echo simulation | |
| CN106950600B (en) | A kind of minimizing technology of near surface scattering surface wave | |
| CN112924964A (en) | Data fusion method and device for stepping frequency radar signal and storage medium | |
| CN113567983B (en) | Radar synthetic aperture sparse imaging method and system using slide rail | |
| CN108845318B (en) | Satellite-borne high-resolution wide-range imaging method based on Relax algorithm | |
| CN108594229B (en) | Satellite-borne SAR intra-pulse Doppler effect two-dimensional compensation method and device and storage medium | |
| US5091890A (en) | Method of extracting target range and Doppler information from a Doppler-spread signal | |
| Li et al. | A fast non-searching method for maneuvering target detection based on ACCF | |
| CN110308448B (en) | Method for enhancing two-dimensional image of inverse synthetic aperture radar |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |