Background
The modern radar is faced with increasingly complex tasks, and in order to improve the detection performance of the radar in a complex and variable environment, a plurality of advanced theories and methods are proposed, and a clutter suppression technology is one of key technologies. Moving object display technology, which is one of the earliest technologies for suppressing clutter, uses the doppler shift of the echo signal of a moving object to distinguish between a fixed object and a moving object. The low repetition frequency MTI radar increases the transmitted energy by transmitting a long pulse width, theoretically if the duty cycle and the total time are the same, the detection power of different repetition frequency radars is the same, because the average power and energy are the same, and in fact, the actual detection distance of the low repetition frequency radar is farther because there is no accumulated loss. The repetition frequency is low, the unambiguous Doppler frequency of the low repetition frequency radar is low, the velocity ambiguity is serious, the low repetition frequency radar is not suitable for pulse Doppler processing, the clutter resistance is weak, but with the development of the technology, the improvement factor of the current advanced MTI radar can reach 60 decibels, the capability of detecting a moving target in clutter can be effectively improved, and the anti-jamming capability of the radar is improved.
At present, the low-repetition-frequency MTI radar is mainly applied to the fields of remote early warning, space surveying and mapping and the like due to the characteristic of clear ranging, and due to the radar system, the low-repetition-frequency MTI radar has serious speed ambiguity when measuring the speed and is not suitable for measuring the speed. However, if the speed information of the target can be obtained, the point trace filtering can be performed by using the speed difference between the target and the clutter, so that a better clutter suppression effect can be achieved. In addition, the speed information of the target can improve the tracking precision of the target and also can roughly estimate the type of the target.
In recent years, the electromagnetic environment in which radar is located has become increasingly complex. In order to improve the anti-interference capability, the radar mostly adopts a frequency agility working mode. Pulse-to-pulse frequency agility has the important advantage of increasing the detectability of certain targets, and frequency agility also mitigates the deleterious effects of echo flicker in tracking radar, facilitating more accurate target tracking. In military radars, inter-pulse frequency agility will force enemy interference signal energy to be spread out over a wide bandwidth, rather than concentrating all energy within the narrow bandwidth of fixed frequency radars.
At present, the study on MTI radar speed measurement by domestic scholars is less. The MTI radar speed measurement algorithm based on phase unwrapping is provided by the inventor of Zweiweicheng et al in China, the algorithm adopts the phase unwrapping and the repetition frequency dispersion ambiguity solution method to realize the MTI radar speed measurement, but the algorithm needs to send six repetition frequency dispersion periods for ambiguity resolution, compared with a radar without ambiguity, redundant waveforms mean the energy loss with serious difference and can seriously affect the action distance of the radar, and waveforms sending the same working frequency for a long time are easily interfered by narrow-band aiming, so that target detection cannot be effectively realized, the ambiguity resolution effect difference of different repetition periods is large, and in order to improve the effect, special optimization design needs to be carried out on the repetition frequency periods. In recent years, in the field of radar technology, the development of multi-channel transmitting and receiving technology is rapid, but the current research mainly focuses on transmitting and receiving a single working frequency, that is, all transmitting and receiving channels are consistent, and the structural design is simple to implement, but the anti-interference performance is weak. The 2011 li founding nation proposes a linear channel and logarithmic channel dual-channel receiver, which not only completely retains the amplitude signal of the signal, but also expands the dynamic range of the receiver, but does not research the electronic impedance performance.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a speed measuring method of a frequency agile MTI radar, which overcomes the defects of weak anti-interference capability and serious energy loss of transmitted redundant waveform of the MTI radar in the prior art. Because single frequency transmission is adopted and redundant waveforms are transmitted in a repeated frequency staggered mode to solve the real speed of a target, the same frequency signal is transmitted for a long time and is extremely easy to be intercepted by a reconnaissance receiver and implement narrow-band aiming type interference, and the redundant waveforms transmitted in the repeated frequency staggered mode mean serious energy loss and can seriously influence the action power of a radar.
The technical scheme adopted by the invention for solving the technical problem comprises the following specific steps:
S1:in order to reduce the probability of interception of the transmitted signal by the reconnaissance receiver, the waveform generator generates a waveform and divides the waveform into two paths to be respectively sent to two transmitting channels, and the two transmitting channels respectively receive the local oscillation signal f sent by the local oscillation sourceLO1And fLO2Require | fLO1-fLO2If the frequency is more than 1GHz, the two paths of signals respectively radiate signals to a space target through the same antenna at different transmitting frequencies after up-conversion;
s2: in order to resist narrow-band aiming type interference, an antenna collects a target reflection signal and sends the target reflection signal to two receiving channels, and the two receiving channels respectively receive a local oscillation signal f sent by a local oscillation sourceLO1And fLO2Respectively completing the preprocessing of signal amplification, down conversion, AD conversion and digital down conversion, and then sending the signal to a signal processor;
s3: in order to improve the target detection performance, the signals preprocessed by the two receiving channels respectively complete pulse compression and MTI processing, then the signals processed by the MTI processing of the two receiving channels are subjected to non-coherent fusion to improve the signal-to-noise ratio, and target constant false alarm detection is completed, wherein the specific contents of the non-coherent fusion are as follows:
receiving channel-MTI processed signal M11(t)、M12(t) are respectively:
wherein A is
1Representing the amplitude of a signal of the receiving channel, f
d1Representing the Doppler frequency, T
rWhich represents the pulse repetition period, is shown,
representing an initial phase value;
receiving channel two-MTI processed signal M21(t)、M22(t) is:
in the same way, there are
Wherein A is
2Representing the amplitude of the two signals of the receiving channel, f
d2Representing the Doppler frequency, T
rThe pulse weight gain period is shown as,
representing an initial phase value;
receiving channel-MTI processed signal M11(t)、M12(t) and reception channel two MTI processed Signal M21(t)、M22(t) performing modulo value operation to obtain | M11(t)|、|M12(t)|、|M21(t)|、|M22(t) |, and then performing non-coherent fusion, namely:
M(t)=|M11(t)|+|M12(t)|+|M21(t)|+|M22(t)|
completing target detection on the non-coherent fused data M (t) by adopting a unit average constant false alarm algorithm to obtain a detection result D;
s4: in order to obtain the speed information of the target, the signals processed by the two receiving channels MTI in S3 and the detection result of the target constant false alarm are used to complete the resolution of the target speed by using the phase information and the screening method, which includes the following specific contents:
s41: according to the receiving channel-MTI processed signal M11(t) and reception channel two MTI processed Signal M12(t), and let initial time t be 0, calculate the phase value:
Phase11=-2πfd1Tr+π-φ01
Phase12=-4πfd1Tr+π-φ01
the phase difference value is:
phase_dif1=Phase11-Phase12=2πfd1Tr
this gives:
s42: according to the cancellation result M of the three pulses in S321(t) and M22(t), and let initial time t be 0, calculate the phase value:
Phase21=-2πfd2Tr+π-φ01
Phase22=-4πfd2Tr+π-φ01
the phase difference value is as follows:
phase_diff2=Phase21-Phase22=2πfd2Tr
this gives:
s43: will f isd1All possible corresponding target speeds are listed, namely:
N1=floor(2VmaxfR1/(frc) wherein V) ismaxTo target maximum possible speed, fR1For the transmit channel, a transmit frequency, c the speed of light, fr=1/TrIs the pulse repetition frequency;
will f isd2All possible corresponding target speeds are listed, namely:
N2=floor(2VmaxfR2/(frc) wherein V) ismaxTo target maximum possible speed, fR2For the emission channel two emission frequencies, c is the speed of light, fr=1/TrIs the pulse repetition frequency;
will be provided with
Each value of (1) is respectively
The difference value of each value is calculated and the absolute norm is obtained
In that
Search for the minimum value, if
Is the minimum value, then k
1Or k
2The target correct speed is finally obtained for the ambiguity
S5: detecting the target D in S3 and the speed information V of the target in S4rAnd sending the data to a subsequent data processing subsystem together for further completing target tracking.
The invention has the beneficial effects that: first, since the present invention employs simultaneous dual frequency transmission, it is possible
The probability of a reconnaissance receiver intercepting radar transmitted signals is reduced, and meanwhile, narrow-band aiming type interference can be effectively resisted by double-frequency receiving, so that the defect that the electronic countermeasure capability in the existing MTI radar is weak is overcome, and the method has the advantage of strong anti-interference capability; secondly, the invention adopts the dual-channel signal level fusion to realize the target detection and does not need redundant waveforms to realize the speed ambiguity resolution, thereby overcoming the defect of serious energy loss by transmitting the redundant waveforms in the prior art and effectively improving the comprehensive detection performance of the target. Therefore, the invention effectively solves the problems of weak anti-interference capability and low utilization rate of transmitting energy of the MTI radar, and provides powerful technical support for realizing remote target detection in a complex electromagnetic environment.
Detailed Description
The invention is further illustrated with reference to the following figures and examples.
The invention utilizes the technology of simultaneously transmitting and receiving signals with different frequencies by two channels, thereby not only effectively improving the detection and speed measurement performance of the target, but also improving the comprehensive anti-interference capability and providing powerful technical support for the practical use of the method in the MTI radar.
The invention provides a frequency agility MTI radar speed measurement method, which not only can effectively detect the speed of a target, but also can obviously improve the detection performance and the comprehensive anti-interference capability of a radar, and lays a technical foundation for improving the comprehensive detection capability of the radar in a defense area.
Referring to fig. 1, a schematic block diagram of a frequency agile MTI radar system structure and a speed measurement method according to the present invention is shown. The invention relates to a method for measuring the speed of a frequency agility MTI radar, which comprises the following steps:
s1: referring to fig. 1, a waveform generator generates four pulse periodic waveforms and divides the waveforms into two paths of signals, one of which is sent to a first transmission channel and a local oscillation signal f sent from a local oscillation sourceLO1Performing up-conversion, filtering, amplifying, and processing by circulator at RF frequency fR1Radiating a signal to a spatial target; at the same time, the other path of signal is sent to the local oscillation signal f sent by the second transmission channel and the local oscillation sourceLO2Performing up-conversion, filtering, amplifying, and processing by circulator at RF frequency fR2Radiating a signal to a target in space, requiring | fLO1-fLO2| > 1GHz, i.e. | fR1-fR2|>1GHz。
S2: the antenna collects the target reflection signal and sends the target reflection signal to two receiving channels through a circulator: receiving a local oscillation signal f sent by a local oscillation source after a target echo signal is received by a first receiving channel and is subjected to amplification and filtering processing
LO1Finishing down-conversion treatment, and obtaining a four-pulse signal set of a receiving channel I through AD conversion and digital down-conversion filtering treatment after filtering and amplifying
Namely:
T
rdenotes the pulse repetition period, r
1It is indicated that the receiving channel one,
indicates the first pulse signal of the receiving channel, and l indicates the pulse number.
The second receiving channel receives the local oscillation signal f sent by the local oscillation source after the target echo signal is amplified and filtered
LO1Finishing down-conversion treatment, and obtaining a four-pulse signal set of a second receiving channel through AD conversion and digital down-conversion filtering treatment after filtering and amplifying
Namely:
here T
rDenotes the pulse repetition period, r
2It is indicated that the receiving channel one,
indicates the first pulse signal of the receiving channel, and l indicates the pulse number.
S3: the signals preprocessed by the two receiving channels respectively complete pulse compression and MTI processing, and then the signals preprocessed by the two receiving channels MTI are subjected to non-coherent fusion, and the implementation is as follows:
s31: receiving channel one according to the preprocessed signal S in S2
r1(t) respectively completing the pulse compression treatment to obtain the result after the pulse compression
The signal model is represented as:
wherein A is
1Representing the amplitude of a signal of the receiving channel, f
d1Representing the Doppler frequency, T
rWhich represents the pulse repetition period, is shown,
representing the initial phase value.
And (3) three-pulse cancellation MTI treatment:
s32: the second receiving channel is based on the preprocessed signal S in S2
r2(t) respectively completing the pulse compression treatment to obtain the result after the pulse compression
Signal moduleThe type is represented as:
wherein A is
2Representing the amplitude of the two signals of the receiving channel, f
d2Representing the Doppler frequency, T
rThe pulse weight gain period is shown as,
representing the initial phase value.
And (3) three-pulse cancellation MTI treatment:
in the same way, there are
S33: receiving channel-MTI processed signal M11(t)、M12(t) and reception channel two MTI processed Signal M21(t)、M22(t) performing modulo value operation to obtain | M11(t)|、|M12(t)|、|M21(t)|、|M22(t) |, and then subjected to non-coherent fusion, i.e.
M(t)=|M11(t)|+|M12(t)|+|M21(t)|+|M22(t)|
S34: and (4) completing target detection on the data M (t) after non-coherent fusion in the S33 by adopting a traditional unit average constant false alarm algorithm to obtain a detection result D.
S4: and (3) utilizing the signals processed by the two receiving channels MTI in S31 and S32, the target constant false alarm detection result and the phase information to complete the calculation of the target real speed by adopting a screening method, and concretely implementing the following steps:
s41: according to the receiving channel-MTI processed signal M11(t) and reception channel two MTI processed Signal M12(t), and let initial time t be 0, calculate the phase value:
Phase11=-2πfd1Tr+π-φ01
Phase12=-4πfd1Tr+π-φ01
the phase difference value is:
phase_diff1=Phase11-Phase12=2πfd1Tr
this gives:
s42: according to the cancellation result M of the three pulses in S321(t) and M22(t), and let initial time t be 0, calculate the phase value:
Phase21=-2πfd2Tr+π-φ01
Phase22=-4πfd2Tr+π-φ01
the phase difference value is as follows:
phase_diiff2=Phase21-Phase22=2πfd2Tr
this gives:
s43: will f isd1All possible corresponding target speeds are listed, namely:
N1=floor(2VmaxfR1/(frc) wherein V) ismaxTo target maximum possible speed, fR1For the transmit channel, a transmit frequency, c the speed of light, fr=1/TrIs the pulse repetition frequency;
will f isd2All possible corresponding target speeds are listed, namely:
N2=floor(2VmaxfR2/(frc) wherein V) ismaxTo target maximum possible speed, fR2For the emission channel two emission frequencies, c is the speed of light, fr=1/TrIs the pulse repetition frequency;
will be provided with
Each value of (1) is respectively
The difference value of each value is calculated and the absolute norm is obtained
In that
Search for the minimum value, if
Is the minimum value, then k
1Or k
2The target correct speed is finally obtained for the ambiguity
S5: detecting the target D in S3 and the speed information V of the target in S4rConstituent target detection information R ═ { D, VrAnd sending the data to a subsequent data processing subsystem to further complete target tracking.
The effect of the present invention is further explained by simulation experiments.
Simulation experiment contents: setting the transmitting waveform as linear frequency modulation signal with time width of 100 mus, bandwidth of 4MHz and pulse repetition period Tr1000 mus, target speed [100m/s,1500m/s]Randomly generated, radar transmission channel-transmission frequency fR10.9GHz, two transmitting frequencies f of radar transmitting channelR2The number of monte carlo tests was 1000 times at 2 GHz. Experiments are carried out in MATLAB13.0a software, and target speed calculation is carried out according to the method provided by the invention to obtain the real speed of the target. Referring to fig. 2, a schematic diagram of speed correct resolution probability variation under different signal-to-noise ratios is shown, in fig. 2, a horizontal axis represents the signal-to-noise ratio (unit: decibel), and a vertical axis represents the detection probability; referring to fig. 3, it is a detection probability curve of the method of the present invention and the conventional method under the condition of the same false alarm rate and different signal-to-noise ratios, in fig. 3, the horizontal axis represents the signal-to-noise ratio (unit: decibel) and the vertical axis represents the detection probability.
From fig. 2, it can be seen that the method provided by the invention has a high probability of correctly resolving the target speed under a certain signal-to-noise ratio; from fig. 3, it can be seen that the method provided by the present invention has better detection performance compared with the conventional method, thereby proving the effectiveness of the present invention.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.