CN105974389A - Mechanic scan meter-wave radar monopulse angle measuring method based on iterative processing - Google Patents
Mechanic scan meter-wave radar monopulse angle measuring method based on iterative processing Download PDFInfo
- Publication number
- CN105974389A CN105974389A CN201610504677.8A CN201610504677A CN105974389A CN 105974389 A CN105974389 A CN 105974389A CN 201610504677 A CN201610504677 A CN 201610504677A CN 105974389 A CN105974389 A CN 105974389A
- Authority
- CN
- China
- Prior art keywords
- angle
- pulse
- target
- signal
- echo
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- 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
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
Abstract
The invention belongs to the technical field of radars, and discloses a mechanic scan meter-wave radar monopulse angle measuring method based on iterative processing. The method can reduce the influence of noise on signals and reduce the angle measuring error. The method comprises the steps of: uniformly dividing an antenna into two sub-arrays to transmit pulse signals; performing coherent accumulation on the received echo signals to obtain accumulated data; performing phase adjustment on the received echo signals, selecting a maximum amplitude, and finding the angle value theta max of the corresponding pulse; separately selecting (L-1)/2 angle values theta 1 on the left and the right of the angle theta max, and finding corresponding amplitudes, thus obtaining two groups of new data; obtaining a sum beam and a difference beam from the two groups of data; performing sum-difference beam monopulse angle measurement on the sum beam and the difference beam to obtain a target off-axis angle; endowing the off-axis angle to the theta max; and repeating the above steps till the target estimated angle satisfies a preset error condition.
Description
Technical field
The present invention relates to Radar Technology field, particularly relate to oneMetre wave radar Monopulse estimation method swept by machine based on iterative processing, can be used for when center of antenna normal in mechanical scanning radar is not aligned with target, target being carried out pulse method angle measurement.
Background technology
Metre wave radar is little due to signal attenuation, and detection range is remote, has the advantage of uniqueness at the aspect such as over-the-horizon detection, Anti-amyloid-β antibody, but simultaneously because the wavelength of metre wave radar is longer, wave beam is wider so that its angle-resolved rate variance, angle measurement accuracy is low.Mechanical scanning radar is due to low cost, it is achieved the factor such as simple is extensively applied by people always.Mechanic scan radar angle-measuring method typically uses maximum-signal method, but the precision of this method is poor, sweeps metre wave radar especially for machine.Its electricity size is less, and wave beam is wider, and maximum-signal method angle measurement accuracy is the highest.In order to improve angle measurement accuracy, Monopulse estimation technology can be used.
Zhang Guangyi academician gives Monopulse estimation method in " phased-array radar principle [M], Beijing: National Defense Industry Press, 1994 " book.Monopulse estimation method refers to the class angle-measuring method utilizing the echo-signal of individual pulse just can measure target actual position.Width phase and difference beam single pulse method is utilized can quickly to obtain the accurate location of target in metre wave radar swept by machine.Width phase and difference beam single pulse method refer to each pulse echo that radar is received after treatment, obtain each pulse and wave beam and difference beam, by calculating difference and ratio, table look-up and understand the off-axis angle of each pulse target.Adding the reference angle of antenna, just obtain the angle on target that individual pulse records, it is exactly the actual position of target that the angle recording all pulses is averaged.
Owing to width phase and difference beam single pulse method realize conveniently, amount of calculation is little, and angle measurement accuracy is high, and tool has great advantage in actual applications.But swept the shortcomings such as wider, the central beam misalignment target of metre wave radar wave beam to be affected by machine, although width phase and difference beam Monopulse estimation method precision are higher than maximum-signal method, but its angle measurement accuracy is affected by noise the biggest, precision needs to improve further, especially when antenna centre normal is not aligned with target, its angle measurement accuracy is the lowest, thus affects resolving power and detection, the accuracy of tracking target of radar.
Summary of the invention
For the deficiency of above-mentioned prior art, it is an object of the invention to provide oneMetre wave radar Monopulse estimation method swept by machine based on iterative processing, thus reduce the noise impact on signal, reduce angle error, improve angle measurement accuracy.
The technical scheme is that machine is swept antenna equivalence becomes the antenna surface being made up of two submatrixs, during antenna scanning, the echo-signal received is carried out Subarray process, the data obtained after processing are formed and difference beam, with sweeping the array radar echo-signal to receiving opportunity, enter horizontal phasing control, the data that obtain after processing taking maximum amplitude, and finds the angle value of respective pulses, all take (L-1)/2 angle value (L is beam direction in the left and right of this angleFigure 3 dThe pulse number that B bandwidth comprises), and find out the corresponding and amplitude of difference beam, use width phase and difference beam single burst algorithm to carry out angle estimation, obtain angle estimation value θe.Again at θeLeft and right all take (L-1)/2 angle value, and find out the corresponding and amplitude of difference beam, and re-use width phase and difference beam single burst algorithm carries out angle estimation, an angle estimation value of getting back, then proceed to repeat above operation, carry out repeatedly width phase and difference beam pulse iterative processing.
For reaching above-mentioned purpose, embodiments of the invention adopt the following technical scheme that
A kind ofMetre wave radar Monopulse estimation method swept by machine based on iterative processing, described method comprises the steps:
Step 1, the machine of setting sweeps the antenna of metre wave radar as the uniform line-array being made up of N number of array element, described uniform line-array is divided into the left submatrix that array number is N/2 and the right submatrix that array number is N/2, and wherein, N is the even number more than or equal to 4；
Step 2, the measurable angle range-90 °～90 ° that machine is swept metre wave radar is divided into multiple ripple position, is spaced the first preset angle between adjacent wave position；Determine the mirror angular curve of each ripple position, composition mirror angular curve table；
Step 3, machine is swept metre wave radar and is launched K pulse, is spaced the second preset angle, determines the reference angle of i-th pulse between adjacent pulseI=1,2 ..., K；Wherein, reference angle is the angle of center of antenna normal and horizontal reference plane；
Step 4, machine is swept metre wave radar and is received echo-signal X of K pulse by antenna array, determines the target off-axis angle of i-th pulse, i=1,2, ..., K, and echo-signal X is entered horizontal phasing control, obtain echo-signal Y after phase place adjusts, obtain the pulse that in echo-signal Y after phase place adjusts, maximum amplitude is corresponding, and the target off-axis angle that this pulse is corresponding, be designated as first object off-axis angle；Off-axis angle is the angle between target and antenna normal；
Step 5, obtains echo-signal T that left submatrix receives respectively1Echo-signal T received with right submatrix2, on the left of first object off-axis angle, take (L-1)/2 angle value, on the right side of first object off-axis angle, take (L-1)/2 be spaced apart the second preset angle between angle value, and adjacent angular；And echo-signal T received from left submatrix1In obtain corresponding with (L-1)/2 angle value on (L-1)/2 angle value on the left of first object off-axis angle and right side range value composition the first new echo-signal Z1, echo-signal T that receives from right submatrix2In obtain corresponding with (L-1)/2 angle value on (L-1)/2 angle value on the left of first object off-axis angle and right side range value composition the second new echo-signal Z2；L is the pulse number comprised in 3dB beam angle；
Step 6, according to the described first new echo-signal Z1With the second new echo-signal Z2Obtain and wave beam and difference beam, thus differed from and ratio, according to described difference with than target off-axis angle θ obtaining i-th pulse with described mirror angular curve tablei, i=1,2 ..., K；By target off-axis angle θ of i-th pulseiReference angle with i-th pulseIt is added, obtains the i-th pulse measurement angle ψ to targeti, i=1,2 ... L；
Step 7, repeated execution of steps 6, until obtaining L the pulse measurement angle to target, L measurement angle being averaged, obtains the estimation angle of target；
Step 8, the first object off-axis angle that the estimation angle of described target is assigned in step 4；
Step 9, is repeated in performing step 5 and meets default error condition, and the angle estimation value that estimation angle-determining is target that will obtain for the last time to step 8, the estimation angle until target.
The present invention compared with prior art has the advantage that the present invention carries out repeatedly pulse method to receiving data, by the advantage successive ignition of pulse method, thus whether center of antenna is weakened pulse method to the shortcoming that quasi goal is the most sensitive, thus improve the angle measurement accuracy of pulse method.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below by use required in embodiment or description of the prior artAccompanying drawingIt is briefly described, it should be apparent that, in describing belowAccompanying drawingIt is only some embodiments of the present invention, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to according to theseAccompanying drawingObtain otherAccompanying drawing。
Figure 1It it is the one of embodiment of the present invention offerMetre wave radar Monopulse estimation method swept by machine based on iterative processingRealize flow process signalFigure；
Figure 2It is tradition width phase Monopulse estimation method and the signal of difference beam curveFigure；
Figure 3It it is the mirror angular curve signal of tradition width phase Monopulse estimation methodFigure；
Figure 4It it is the root-mean-square error correlation curve signal using the inventive method and tradition width phase single pulse method that one target carries out angle measurementFigure。
Detailed description of the invention
Below in conjunction with in the embodiment of the present inventionAccompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments.Based on the embodiment in the present invention, the every other embodiment that those of ordinary skill in the art are obtained under not making creative work premise, broadly fall into the scope of protection of the invention.
The embodiment of the present invention provides oneMetre wave radar Monopulse estimation method swept by machine based on iterative processing, referenceFigure 1, described method comprises the steps:
Step 1, the machine of setting sweeps the antenna of metre wave radar as the uniform line-array being made up of N number of array element, described uniform line-array is divided into the left submatrix that array number is N/2 and the right submatrix that array number is N/2, and wherein, N is the even number more than or equal to 4.
Step 2, the measurable angle range-90 °～90 ° that machine is swept metre wave radar is divided into multiple ripple position, is spaced the first preset angle between adjacent wave position；Determine the mirror angular curve of each ripple position, composition mirror angular curve table.
Step 2 particularly as follows:
(2a) at interval of the first preset angle △ θ, measurable angle range-90 °～90 ° being divided into a ripple position, wherein, the first preset angle △ θ=50.8/ (N-d/ λ), d is array element interval, and λ is signal wavelength, d=λ/2；
(2b) for any one ripple position, this ripple position synthesis wave beam B in left submatrix is obtainedlWith this ripple position at the synthesis wave beam B of right submatrixr, according to this ripple position at the synthesis wave beam B of left submatrixlWith this ripple position at the synthesis wave beam B of right submatrixrObtain this ripple position and wave beam and difference beam；
(2c) according to this ripple position and wave beam and difference beam obtain difference and the ratio of this ripple positionAnd obtain the mirror angular curve of this ripple position；Wherein imag () expression takes imaginary-part operation；With the angular range of this ripple position as abscissa, with difference and ratio as vertical coordinate, the difference corresponding by this ripple position measurable angle range and compare the mirror angular curve drawing this ripple position；
(2d) obtain the mirror angular curve of all ripple positions, and the mirror angular curve of all ripple positions is formed mirror angular curve table.
Step 3, machine is swept metre wave radar and is launched K pulse, is spaced the second preset angle, determines the reference angle of i-th pulse between adjacent pulseI=1,2 ..., K；Wherein, reference angle is the angle of center of antenna normal and horizontal reference plane.
In step 3, machine sweep metre wave radar launch K pulse, be spaced the second preset angle between adjacent pulse particularly as follows:
Machine sweeps metre wave radar when antenna scanning works, and launches a pulse at interval of the second preset angle △ Φ, and comprises L pulse in 3dB beam angle；Wherein, △ Φ=360*t/T, t is the pulse repetition period, and T is the rotating speed that metre wave radar swept by machine.
The angle of center of antenna normal and horizontal reference plane is as the reference angle of i-th pulseI=1,2 ..., K, owing to being engraved in rotation when metre wave radar swept by machine, so center of antenna normal is change.
Step 4, machine is swept metre wave radar and is received echo-signal X of K pulse by antenna array, determines the target off-axis angle of i-th pulse, i=1,2, ..., K, and echo-signal X is entered horizontal phasing control, obtain echo-signal Y after phase place adjusts, obtain the pulse that in echo-signal Y after phase place adjusts, maximum amplitude is corresponding, and the target off-axis angle that this pulse is corresponding, be designated as first object off-axis angle；Target off-axis angle is the angle between target and antenna normal.
Step 4 particularly as follows:
(4a) machine is swept metre wave radar and is received echo-signal X=AS+n of K pulse, wherein, A=[a (θ by antenna array1),...a(θi)...,a(θK)] represent the target phase information that echo-signal comprises, a (θi)=[1, exp (j2 π d/ λ sin θi),...,exp(j2π(N-1)d/λsinθi)]TFor the target off-axis angle vector of i-th pulse, i=1,2 ..., K；D is array element interval, and λ is signal wavelength, represents the transposition of vector, and exp represents the exponential depth with e as the end, and j represents imaginary unit；
S=[S1,...Si...,SK]TRepresent the complex envelope of echo-signal, SiRepresent the complex envelope information of i-th pulse,Wherein fdRepresent the Doppler frequency of target, fd=2Vf0/ c, V represent that target sweeps the radial velocity of metre wave radar, f relative to machine0Representing the mid frequency of radar emission signal, c represents the light velocity, and t is the pulse repetition period；N represent average be 0, variance be 1 N × K rank white Gaussian noise matrix；
(4b) echo-signal X is entered horizontal phasing control, obtain the echo-signal Y=a (θ after phase place adjusts0) X, it is assumed that the actual angle of target is θ0a(θ0)=[1, exp (j2 π d/ λ sin (θ0)),…,exp(j2π(N-1)d/λsin(θ0))]T；
It should be noted that the angle that the actual angle of the target described in the present embodiment is the target estimated according to prior information in advance, belonging to the category of rough estimate, the present invention is that the angle to target is accurately estimated.
(4c) pulse that in echo-signal Y after position adjusts, maximum amplitude is corresponding is obtained, and first object off-axis angle θ that this pulse is correspondingmax。
Step 5, obtains echo-signal T that left submatrix receives respectively1Echo-signal T received with right submatrix2, on the left of first object off-axis angle, take (L-1)/2 angle value, on the right side of first object off-axis angle, take (L-1)/2 be spaced apart described second preset angle between angle value, and adjacent angular；And echo-signal T received from left submatrix1In obtain corresponding with (L-1)/2 angle value on the left of first object off-axis angle range value composition the first new echo-signal Z1, echo-signal T that receives from right submatrix2In obtain corresponding with (L-1)/2 angle value on the right side of first object off-axis angle range value composition the second new echo-signal Z2；L is the pulse number comprised in 3dB beam angle.
Wherein, T1It is 1 × kth moment battle array, T2It it is 1 × kth moment battle array.
Step 6, according to the described first new echo-signal Z1With the second new echo-signal Z2Obtain and wave beam and difference beam, thus differed from and ratio, according to described difference with than target off-axis angle θ obtaining i-th pulse with described mirror angular curve tablei, i=1,2 ..., K；By target off-axis angle θ of i-th pulseiReference angle with i-th pulseIt is added, obtains the i-th pulse measurement angle ψ to targeti, i=1,2 ... L.
Step 6 particularly as follows:
(6a) according to the described first new echo-signal Z1With the second new echo-signal Z2Obtain and wave beam PΣWith difference beam PΔ, PΣ=Z1+Z2, PΔ=Z1-Z2；
(6b) utilize and wave beam PΣWith difference beam PΔ, calculate difference and than P=imag (PΛ/PΣ), according to described difference and the value of ratio, obtained target off-axis angle θ of i-th pulse by inquiry mirror angular curve tablei, i=1,2 ..., K；
(6c) by target off-axis angle θ of i-th pulseiReference angle with i-th pulseIt is added, obtains the i-th pulse measurement angle to targetI=1,2 ... L.
Step 7, repeated execution of steps 6, until obtaining L the pulse measurement angle to target, L measurement angle being averaged, obtains the estimation angle of target.
Step 8, the first object off-axis angle that the estimation angle of described target is assigned in step 4；
Step 9, is repeated in performing step 5 and meets default error condition, and the angle estimation value that estimation angle-determining is target that will obtain for the last time to step 7, the estimation angle until target.
Until the estimation angle of target meets default error condition particularly as follows: the estimation angle of target and the root-mean-square error of actual angle are less than preset value in step 9, wherein preset value is 10。
The effect of the present invention can be verified by machine calculated below emulation:
One, simulated conditions
Simulated conditions 1: assume the antenna surface that antenna is made up of 10 array elements, it is divided into two submatrixs, 5 array elements of each submatrix, the beam angle of antenna is about 10 °, antenna carries out rotation sweep with the speed of 10s/r, and it is the pulse of 300MHz that radar launches a frequency every 10ms, and antenna receives 21 pulse signals every time, antenna impulse center is not aligned with target, deviation 4 pulses of target of antenna impulse center；Assuming to there is a target, the real angle of target is 10 °；Signal to noise ratio takes-10dB to 10dB, and signal to noise ratio here refers to single array element, individual pulse；Iterations M=4.
Two, emulation content
Emulation 1, utilizes simulated conditions 1 produces traditional single pulse angle-measuring method respectively and difference beam curve and mirror angular curve, distinguishesSuch as figure 2WithFigure 3Shown in, whereinFigure 2Abscissa is angle, and vertical coordinate is amplitude；Figure 3Abscissa is angle, and vertical coordinate is difference and ratio.
ByFigure 2Visible, when antenna centre normal is to quasi goal, echo-signal and beam amplitude at maximum position, difference beam amplitude is at minimum position, if antenna is not to quasi goal, does not haves above-mentioned situation.
ByFigure 3Visible, when antenna centre normal is to quasi goal, and the difference of difference beam and ratio are 0, this withFigure 2Shown in situation match, comparisonFigure 2During it is found that heart normal is to quasi goal in antennas, the difference beam amplitude of echo is 0, then it differs from and is 0 than necessarily；Similarly, if center of antenna finds to be not aligned with target, the poorest and ratio has deviation, it addition,Figure 3Shown mirror angular curve can be as the foundation of tabling look-up of the inventive method.
Utilize simulated conditions 1 that target is respectively adopted traditional and difference beam single pulse method and the inventive method carry out angle measurement, obtain the root-mean-square error curve that two kinds of methods change with signal to noise ratio,Such as figure 4Shown in,Figure 4Middle abscissa is signal to noise ratio, and vertical coordinate is root-mean-square error.
ByFigure 4Visible, the inventive method ratio is traditional and difference beam single pulse method precision is high, particularly when signal to noise ratio is low, effect clearly,Figure 4Fully demonstrating the effectiveness of the inventive method, the situation of center of antenna normal misalignment target is relatively common in actual applications, and the inventive method can reduce effect of noise, improves angle measurement accuracy.
The above; being only the detailed description of the invention of the present invention, but protection scope of the present invention is not limited thereto, any those familiar with the art is in the technical scope that the invention discloses; change can be readily occurred in or replace, all should contain within protection scope of the present invention.Therefore, protection scope of the present invention should be with describedClaimProtection domain be as the criterion.
Claims (6)
1. a metre wave radar Monopulse estimation method swept by machine based on iterative processing, and its feature exists
In, described method comprises the steps:
Step 1, the machine of setting sweeps the antenna of metre wave radar as the uniform line-array being made up of N number of array element,
Described uniform line-array is divided into the left submatrix that array number is N/2 and the right submatrix that array number is N/2,
Wherein, N is the even number more than or equal to 4；
Step 2, the measurable angle range-90 °～90 ° that machine is swept metre wave radar is divided into multiple ripple position, phase
The first preset angle it is spaced between adjacent wave position；Determining the mirror angular curve of each ripple position, composition mirror angle is bent
Line table；
Step 3, machine is swept metre wave radar and is launched K pulse, is spaced second and presets between adjacent pulse
Angle, determines the reference angle of i-th pulseI=1,2 ..., K；Wherein, reference angle is center of antenna method
Line and the angle of horizontal reference plane；
Step 4, machine is swept metre wave radar and is received echo-signal X of K pulse by antenna array, really
Determine the target off-axis angle of i-th pulse, i=1,2 ..., K, and echo-signal X is entered horizontal phasing control,
Obtain echo-signal Y after phase place adjusts, in echo-signal Y after acquisition phase place adjustment the most significantly
The pulse that value is corresponding, and the target off-axis angle that this pulse is corresponding, be designated as first object off-axis angle；
Target off-axis angle is the angle between target and center of antenna normal；
Step 5, obtains echo-signal T that left submatrix receives respectively1The echo letter received with right submatrix
Number T2, on the left of first object off-axis angle, take (L-1)/2 angle value, in first object off-axis angle
Right side takes (L-1)/2 and is spaced apart described second preset angle between angle value, and adjacent angular；
And echo-signal T received from left submatrix1In obtain with first object off-axis angle on the left of
(L-1) the range value composition first echo signal Z that/2 angle values are corresponding1, receive from right submatrix
Echo-signal T2In obtain corresponding with (L-1)/2 angle value on the right side of first object off-axis angle
Range value composition second echo signal Z2；L is the pulse number comprised in 3dB beam angle；
Step 6, according to described first echo signal Z1With second echo signal Z2Obtain with wave beam with
And difference beam, thus differed from and ratio, obtain according to described difference with than with described mirror angular curve table
Target off-axis angle θ of i-th pulsei, i=1,2 ..., K；By target off-axis angle θ of i-th pulseiWith
The reference angle of i-th pulseIt is added, obtains the i-th pulse measurement angle ψ to targeti,
I=1,2 ... L；
Step 7, repeated execution of steps 6, until obtaining L the pulse measurement angle to target,
L measurement angle is averaged, obtains the estimation angle of target；
Step 8, the first object off-axis angle that the estimation angle of described target is assigned in step 4；
Step 9, being repeated in performing step 5 expires to step 8, the estimation angle until target
The error condition that foot is preset, and the angle that estimation angle-determining is target that will obtain for the last time
Estimated value.
Metre wave radar simple venation swept by a kind of machine based on iterative processing the most according to claim 1
Rush angle-measuring method, it is characterised in that step 2 particularly as follows:
(2a) measurable angle range-90 °～90 ° is divided into a ripple position at interval of the first preset angle Δ θ,
Wherein, the first preset angle Δ θ=50.8/ (N-d/ λ), d is array element interval, and λ is signal wavelength,
D=λ/2；
(2b) for any one ripple position, this ripple position synthesis wave beam B in left submatrix is obtainedlWith
This ripple position is at the synthesis wave beam B of right submatrixr, according to this ripple position at the synthesis wave beam B of left submatrixlWith
This ripple position is at the synthesis wave beam B of right submatrixrObtain this ripple position and wave beam and difference beam；
(2c) according to this ripple position and wave beam and difference beam obtain difference and the ratio of this ripple position, and
Mirror angular curve to this ripple position；
(2d) obtain the mirror angular curve of all ripple positions, and the mirror angular curve of all ripple positions is formed
Mirror angular curve table.
Metre wave radar simple venation swept by a kind of machine based on iterative processing the most according to claim 1
Rushing angle-measuring method, it is characterised in that in step 3, machine is swept metre wave radar and is launched K pulse,
Be spaced the second preset angle between adjacent pulse particularly as follows:
Machine sweeps metre wave radar when antenna scanning works, and launches one at interval of the second preset angle △ Φ
Pulse, and in 3dB beam angle, comprise L pulse；Wherein, △ Φ=360*t/T, t is pulse
Repetition period, T is the rotating speed that metre wave radar swept by machine.
Metre wave radar simple venation swept by a kind of machine based on iterative processing the most according to claim 1
Rush angle-measuring method, it is characterised in that step 4 particularly as follows:
(4a) machine is swept metre wave radar and is received echo-signal X=AS+n of K pulse by antenna array,
Wherein, A=[a (θ1) ... a (θi) ..., a (θK)] represent the target phase information that echo-signal comprises,
a(θi)=[1, exp (j2 π d/ λ sin θi) ..., exp (j2 π (N-1) d/ λ sin θi)]TTarget off-axis for i-th pulse
Angle vector, i=1,2 ..., K；D is array element interval, and λ is signal wavelength, represents the transposition of vector；
S=[S1... Si..., SK]TRepresent the complex envelope of echo-signal, SiRepresent the multiple bag of i-th pulse
Network information,Wherein fdRepresent the Doppler frequency of target, fd=2Vf0/ c, V table
Show that target sweeps the radial velocity of metre wave radar, f relative to machine0Represent the center frequency of radar emission signal
Rate, c represents the light velocity, and t is the pulse repetition period；N represent average be 0, variance be the N × K of 1
Rank white Gaussian noise matrix；
(4b) echo-signal X is entered horizontal phasing control, obtain the echo-signal after phase place adjusts
Y=a (θ0) X, wherein, θ0For the actual angle of target,
a(θ0)=[1, exp (j2 π d/ λ sin (θ0)) ..., exp (j2 π (N-1) d/ λ sin (θ0))]T；
(4c) pulse that in echo-signal Y after position adjusts, maximum amplitude is corresponding is obtained, and should
First object off-axis angle θ that pulse is correspondingmax。
Metre wave radar simple venation swept by a kind of machine based on iterative processing the most according to claim 1
Rush angle-measuring method, it is characterised in that step 6 particularly as follows:
(6a) according to described first echo signal Z1With second echo signal Z2Obtain and wave beam PΣWith
Difference beam PΔ, PΣ=Z1+Z2, PΔ=Z1-Z2；
(6b) utilize and wave beam PΣWith difference beam PΔ, calculate difference and than P=imag (PΛ/PΣ), according to
Described difference and the value of ratio, obtain target off-axis angle θ of i-th pulse by inquiry mirror angular curve tablei,
I=1,2 ..., K；
(6c) by target off-axis angle θ of i-th pulseiReference angle with i-th pulseIt is added,
Obtain the i-th pulse measurement angle to targetI=1,2 ... L.
Metre wave radar simple venation swept by a kind of machine based on iterative processing the most according to claim 1
Rush angle-measuring method, it is characterised in that until the estimation angle of target meets default in step 9
Error condition is particularly as follows: the estimation angle of target and the root-mean-square error of actual angle are less than presetting
Value, wherein preset value is 1 °.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610504677.8A CN105974389B (en) | 2016-06-30 | 2016-06-30 | Machine based on iterative processing sweeps metre wave radar Monopulse estimation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610504677.8A CN105974389B (en) | 2016-06-30 | 2016-06-30 | Machine based on iterative processing sweeps metre wave radar Monopulse estimation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105974389A true CN105974389A (en) | 2016-09-28 |
| CN105974389B CN105974389B (en) | 2018-09-25 |
Family
ID=56953605
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610504677.8A Active CN105974389B (en) | 2016-06-30 | 2016-06-30 | Machine based on iterative processing sweeps metre wave radar Monopulse estimation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105974389B (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108594195A (en) * | 2018-04-20 | 2018-09-28 | 西安电子科技大学 | Low repetition frequency modulated continuous wave radar Aircraft Targets sorting technique based on segmentation mixing |
| CN109061594A (en) * | 2018-08-03 | 2018-12-21 | 中国航空工业集团公司雷华电子技术研究所 | A kind of angle-measuring method based on the processing of submatrix space-time adaptive under strong clutter environment |
| CN109521418A (en) * | 2018-12-28 | 2019-03-26 | 西安电子科技大学 | Ground-based radar angle-measuring method based on interference field |
| CN110018439A (en) * | 2019-03-27 | 2019-07-16 | 中国电子科技集团公司第三十八研究所 | The direction-finding method that a kind of and difference beam is combined with wave beam MUSIC |
| CN110082750A (en) * | 2019-03-25 | 2019-08-02 | 西安电子科技大学 | A kind of measuring angle by comparing amplitude method for eliminating inter-channel phase error |
| WO2019179511A1 (en) * | 2018-03-22 | 2019-09-26 | 京东方科技集团股份有限公司 | Signal receiving method and apparatus, and computer-readable storage medium and electronic device |
| CN110596692A (en) * | 2019-08-19 | 2019-12-20 | 电子科技大学 | Self-adaptive monopulse direction finding method based on joint constraint |
| CN110988835A (en) * | 2019-11-27 | 2020-04-10 | 中国船舶重工集团公司第七二四研究所 | Distributed coherent radar angle measurement method |
| CN112083382A (en) * | 2020-09-15 | 2020-12-15 | 四川九洲空管科技有限责任公司 | High-precision small phased array secondary radar azimuth compensation calculation method |
| CN113189579A (en) * | 2021-03-24 | 2021-07-30 | 四川九洲空管科技有限责任公司 | Hybrid monopulse angle measurement system and method |
| CN114675257A (en) * | 2022-05-30 | 2022-06-28 | 中国人民解放军国防科技大学 | Non-rectangular array polarized monopulse radar double-target resolution method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7619556B1 (en) * | 2008-02-28 | 2009-11-17 | Rockwell Collins, Inc. | System and method for synthesizing localizer and glide slope deviations from weather radar |
| CN101825707A (en) * | 2010-03-31 | 2010-09-08 | 北京航空航天大学 | Monopulse angular measurement method based on Keystone transformation and coherent integration |
| CN101887120A (en) * | 2010-06-10 | 2010-11-17 | 哈尔滨工业大学 | Method for measuring target azimuth by single-beam mechanical scanning radar |
| US20110298653A1 (en) * | 2010-06-04 | 2011-12-08 | Denso Corporation | Method and device for detecting azimuth |
| CN103389493A (en) * | 2013-06-25 | 2013-11-13 | 西安电子科技大学 | Multi-beam single-pulse angle measuring method based on beam selection method |
| CN103728614A (en) * | 2014-01-15 | 2014-04-16 | 西安电子科技大学 | Mechanical scanning meter wave radar based method for improving single pulse angle measurement |
-
2016
- 2016-06-30 CN CN201610504677.8A patent/CN105974389B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7619556B1 (en) * | 2008-02-28 | 2009-11-17 | Rockwell Collins, Inc. | System and method for synthesizing localizer and glide slope deviations from weather radar |
| CN101825707A (en) * | 2010-03-31 | 2010-09-08 | 北京航空航天大学 | Monopulse angular measurement method based on Keystone transformation and coherent integration |
| US20110298653A1 (en) * | 2010-06-04 | 2011-12-08 | Denso Corporation | Method and device for detecting azimuth |
| CN101887120A (en) * | 2010-06-10 | 2010-11-17 | 哈尔滨工业大学 | Method for measuring target azimuth by single-beam mechanical scanning radar |
| CN103389493A (en) * | 2013-06-25 | 2013-11-13 | 西安电子科技大学 | Multi-beam single-pulse angle measuring method based on beam selection method |
| CN103728614A (en) * | 2014-01-15 | 2014-04-16 | 西安电子科技大学 | Mechanical scanning meter wave radar based method for improving single pulse angle measurement |
Non-Patent Citations (4)
| Title |
|---|
| ZY JUN ET.AL: "Adaptive Monopulse Angle Measurement for Meter-Wave Radar Based on Differential Constraints Technique", 《ADVANCED MATERIALS RESEARCH》 * |
| 赵永波等: "一种多目标情况下的单脉冲测角方法", 《西安电子科技大学学报》 * |
| 赵永波等: "一种改进的单脉冲测角方法", 《第九届全国雷达学术年会论文集》 * |
| 陈成增等: "基于鉴角曲线拟合的自适应单脉冲测角方法", 《系统工程与电子技术》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019179511A1 (en) * | 2018-03-22 | 2019-09-26 | 京东方科技集团股份有限公司 | Signal receiving method and apparatus, and computer-readable storage medium and electronic device |
| US11522288B2 (en) | 2018-03-22 | 2022-12-06 | Beijing Boe Technology Development Co., Ltd. | Signal receiving method and apparatus, and computer-readable storage medium and electronic device |
| CN108594195A (en) * | 2018-04-20 | 2018-09-28 | 西安电子科技大学 | Low repetition frequency modulated continuous wave radar Aircraft Targets sorting technique based on segmentation mixing |
| CN109061594A (en) * | 2018-08-03 | 2018-12-21 | 中国航空工业集团公司雷华电子技术研究所 | A kind of angle-measuring method based on the processing of submatrix space-time adaptive under strong clutter environment |
| CN109061594B (en) * | 2018-08-03 | 2022-10-28 | 中国航空工业集团公司雷华电子技术研究所 | Angle measuring method based on subarray space-time adaptive processing in strong clutter environment |
| CN109521418A (en) * | 2018-12-28 | 2019-03-26 | 西安电子科技大学 | Ground-based radar angle-measuring method based on interference field |
| CN109521418B (en) * | 2018-12-28 | 2022-12-02 | 西安电子科技大学 | Foundation radar angle measurement method based on interference field |
| CN110082750B (en) * | 2019-03-25 | 2021-05-07 | 西安电子科技大学 | Amplitude comparison angle measurement method capable of eliminating phase errors among channels |
| CN110082750A (en) * | 2019-03-25 | 2019-08-02 | 西安电子科技大学 | A kind of measuring angle by comparing amplitude method for eliminating inter-channel phase error |
| CN110018439A (en) * | 2019-03-27 | 2019-07-16 | 中国电子科技集团公司第三十八研究所 | The direction-finding method that a kind of and difference beam is combined with wave beam MUSIC |
| CN110596692A (en) * | 2019-08-19 | 2019-12-20 | 电子科技大学 | Self-adaptive monopulse direction finding method based on joint constraint |
| CN110596692B (en) * | 2019-08-19 | 2022-10-14 | 电子科技大学 | Self-adaptive monopulse direction finding method based on joint constraint |
| CN110988835A (en) * | 2019-11-27 | 2020-04-10 | 中国船舶重工集团公司第七二四研究所 | Distributed coherent radar angle measurement method |
| CN112083382A (en) * | 2020-09-15 | 2020-12-15 | 四川九洲空管科技有限责任公司 | High-precision small phased array secondary radar azimuth compensation calculation method |
| CN112083382B (en) * | 2020-09-15 | 2022-08-02 | 四川九洲空管科技有限责任公司 | High-precision small phased array secondary radar azimuth compensation calculation method |
| CN113189579A (en) * | 2021-03-24 | 2021-07-30 | 四川九洲空管科技有限责任公司 | Hybrid monopulse angle measurement system and method |
| CN114675257B (en) * | 2022-05-30 | 2022-08-02 | 中国人民解放军国防科技大学 | Non-rectangular array polarized monopulse radar double-target resolution method |
| CN114675257A (en) * | 2022-05-30 | 2022-06-28 | 中国人民解放军国防科技大学 | Non-rectangular array polarized monopulse radar double-target resolution method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105974389B (en) | 2018-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105974389A (en) | Mechanic scan meter-wave radar monopulse angle measuring method based on iterative processing | |
| CN103728614B (en) | The improvement Monopulse estimation method of metre wave radar is swept based on machine | |
| CN111537966B (en) | Array antenna error correction method suitable for millimeter wave vehicle-mounted radar field | |
| CN103744077B (en) | The angle-measuring method of metre wave radar under multi-target condition swept by machine | |
| CN105974390B (en) | Ji Sao metre wave radars mass centre angle-measuring method based on doppler information | |
| CN102156279B (en) | Method for detecting moving target on ground by utilizing bistatic radar based on MIMO (Multiple Input Multiple Output) | |
| CN107561508B (en) | Coherent accumulation detection method for uniformly accelerated moving target | |
| CN103823217B (en) | Based on the bistatic MIMO radar high-speed moving object method for parameter estimation of double frequency transmitting | |
| CN104155647B (en) | The method sweeping the estimation azimuth of target of meter wave array radar based on machine | |
| CN107976660B (en) | Missile-borne multi-channel radar ultra-low-altitude target analysis and multi-path echo modeling method | |
| CN103412301B (en) | Metrewave radar angle measurement method | |
| CN105445701B (en) | The pulse angle estimating method of DDMA MIMO radar targets | |
| CN112444811A (en) | Target detection and imaging method integrating MIMO radar and ISAR | |
| CN107329138B (en) | Distance walking correction and coherent accumulation detection method for PD radar | |
| CN109212489B (en) | FDA-MIMO radar fuzzy clutter suppression method based on auxiliary pulse | |
| CN106872969B (en) | Radar target angle estimation method based on MTD pulse accumulation and sliding processing | |
| CN112834980B (en) | Vortex electromagnetic wave super-resolution direction finding method based on propagation operator | |
| CN109655819B (en) | Clutter suppression three-dimensional imaging method based on real-aperture Doppler beam sharpening | |
| CN110850400A (en) | LFMCW radar multi-target detection method based on interferometer direction finding | |
| CN104122548B (en) | Meter wave array radar angle-measuring method swept by machine | |
| CN109884337B (en) | Method for detecting sea surface wind direction by using high-frequency ground wave radar | |
| CN109521418B (en) | Foundation radar angle measurement method based on interference field | |
| CN116148784B (en) | Automatic phase calibration system and method for single-station lightning positioning system | |
| CN106125058B (en) | The improvement mass centre angle-measuring method of metre wave radar is swept based on machine | |
| Zhou et al. | Angle measurement accuracy analysis of sum-difference amplitude-comparison monopulse in onshore or shipborne ISAR |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |