CN104199004A - Phase correction method for radar transmission channels - Google Patents

Phase correction method for radar transmission channels Download PDF

Info

Publication number
CN104199004A
CN104199004A CN201410105470.4A CN201410105470A CN104199004A CN 104199004 A CN104199004 A CN 104199004A CN 201410105470 A CN201410105470 A CN 201410105470A CN 104199004 A CN104199004 A CN 104199004A
Authority
CN
China
Prior art keywords
radar
transmission channel
rotation
phase
channel
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
Application number
CN201410105470.4A
Other languages
Chinese (zh)
Other versions
CN104199004B (en
Inventor
杨明磊
陈伯孝
夏碧君
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201410105470.4A priority Critical patent/CN104199004B/en
Priority claimed from CN201410105470.4A external-priority patent/CN104199004B/en
Publication of CN104199004A publication Critical patent/CN104199004A/en
Application granted granted Critical
Publication of CN104199004B publication Critical patent/CN104199004B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/40—Means for monitoring or calibrating
    • G01S7/4004—Means for monitoring or calibrating of parts of a radar system
    • G01S7/4008—Means for monitoring or calibrating of parts of a radar system of transmitters

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention belongs to the technical field of consistence correction of radar transmission channels, and discloses a phase correction method for radar transmission channels. The phase correction method for the radar transmission channels comprises the following steps: acquiring I/Q signals echoed by a plurality of transmission channels of a radar; converting the I/Q signals into a plurality of vectors, and rotating the plurality of vectors to a first quadrant or a fourth quadrant if the plurality of vectors are positioned in a second quadrant or a third quadrant; performing vector rotation based on one-way search till a set accuracy is reached, and determining a final phase corresponding to each transmission channel; performing phase correction on each transmission channel according to the final phase corresponding to each transmission channel.

Description

The bearing calibration of a kind of radar emission channel phases
Technical field
The invention belongs to radar transmit-receive passage consistency alignment technique field, particularly a kind of radar emission channel phases bearing calibration.
Background technology
Modern radar adopts digital array antenna more, and this makes the coverage of radar wider, and antijamming capability is stronger.But being activated at of the each array element of digital array antenna receiving-sending passage (being transmission channel and receiving cable) exists inconsistency in phase place.The reason that causes these phase place inconsistencies is a lot, for example, bay damages, the error of each functional module of transmission channel or receiving cable inside itself, amplitude phase error between the each array element of feeder line, the impedance variation of the bay that the mutual coupling between bay causes etc.Phase place inconsistency can make a big impact to the minor level of antenna lobe, antenna gain and beam position etc.Therefore the passage that transmits and receives that must pair array antenna carries out phase correction.
Phase correction is divided into transmitting and proofreaies and correct and receive correction, the point-frequency signal of being launched successively same time length by transmission channel to be corrected is proofreaied and correct in transmitting, and obtain i/q signal (wherein I signal represents real component, and Q signal is imaginary part component) by same correction channel reception signal by quadrature sampling.Receive to proofread and correct and launch the point-frequency signal of corresponding time span by a corrected transmission channel, each receiving cable to be corrected receives signal simultaneously, obtains i/q signal by quadrature sampling.Receive to proofread and correct or transmitting correction no matter be, the treatment scheme obtaining after i/q signal is all the same, the i/q signal that utilizes sampling to obtain calculates the phase place of each passage to be corrected, a selected passage to be corrected is reference channel again, other each passage to be corrected is all obtained phase differential taking this reference channel as benchmark, and this phase differential is the correction coefficient of each passage to be corrected.In the process of channel correcting, there is a considerable link, utilize i/q signal to ask phase place.
Traditional i/q signal that utilizes asks phase method for position to comprise look-up table, polynoimal approximation, table look-up and polynomial expression combined techniques, digit-by-digit method.But these methods all can not meet the demands in speed and precision and hardware is realized difficulty.Therefore rotation of coordinate numerical calculation method (being cordic algorithm) arises at the historic moment.Sense of rotation and the step-length of traditional C ORDIC algorithm all change with iterations.The sense of rotation of coordinate becomes according to vectorial position, if vector is in the second quadrant, first by vector around initial point dextrorotation turnback, if vector in third quadrant, first vector is rotated counterclockwise to 180 degree around initial point.So just all vectors are all concentrated on to first, fourth quadrant, after this, if being in first quartile, vector rotates clockwise, if be in fourth quadrant by being rotated counterclockwise, the anglec of rotation is successively decreased by 2 integral number power, until meet predetermined number of revolutions.Ask phase method for position with respect to traditional utilize I, Q signal, the precision of this cordic algorithm and speed all obviously improve, and implementation complexity also obviously reduces, but this kind of each iteration of cordic algorithm all needs to judge sense of rotation, have increased the complexity realizing; The anglec of rotation is all successively successively decreased at every turn, has reduced search speed; Thereby reduce the efficiency that radar transmit-receive channel phases is proofreaied and correct.In addition, the iterations of this kind of cordic algorithm is fixed, and possible maximum error is also just more fixing, for some low-angle selections, there will be the situation that error rate is higher.
Summary of the invention
The object of the invention is to propose the bearing calibration of a kind of radar emission channel phases.Can the in the situation that of faster speed, obtain higher phase correction precision, reduce implementation complexity simultaneously, make the transceiver channel Concordance of radar can be more accurate sooner.
For realizing above-mentioned technical purpose, the present invention adopts following technical scheme to be achieved.
The bearing calibration of a kind of radar emission channel phases, is characterized in that, comprises the following steps:
S1: utilize the i/q signal of the echo of i/q signal to the N transmission channel of the 1st transmission channel echo of arbitrary receiving cable acquisition radar of radar, the number of the transmission channel that N is radar;
S2: final phase place corresponding to each transmission channel of determining radar; The final phase place corresponding to each transmission channel of described definite radar comprises the following steps:
S21: in plane right-angle coordinate, according to the phase place of the i/q signal of the echo of j transmission channel and amplitude, by corresponding complex vector located ω for the i/q signal of the echo of j transmission channel jrepresent, j gets 1 to N; According to ω jplace quadrant, by ω jbe rotated ω around initial point jafter rotation, become complex vector located ξ j, draw rotation start-phase z jand ξ jcorresponding plural imaginary part y 0;
S22: parameters d, if ξ jbe positioned at first quartile, make d=-1; Otherwise, make d=1; Variable i and k are set, and i gets odd number, and the initial value of i is 1; K gets 1,2, and 3 In the time of k=1, establish complex vector located ξ j, i, k=ξ j, and establish z j, i, k=z j; ξ j, i, kcorresponding plural imaginary part y j, i, kfor y 0;
S23: by ξ j, i, kbe rotated around initial point, its anglec of rotation is θ i, θ i=arctan (2 -i), sense of rotation is by the decision of d value, and in the time of d=-1, sense of rotation is clockwise direction; In the time of d=1, sense of rotation is counterclockwise; ξ j, i, kafter rotation, become complex vector located ξ j, i, k+1; Draw ξ j, i, k+1corresponding rotatable phase z j, i, k+1, z j, i, k+1=z j, i, k-d θ i; Draw ξ j, i, k+1corresponding plural imaginary part y j, i, k+1;
S24: if | y j, i, k+1| be less than setting accuracy ε, obtain the final phase place α that j transmission channel is corresponding j, α j=z j, i, k+1; If | y j, i, k+1| be more than or equal to setting accuracy ε, judge y j, i, k+1with y j, i, kwhether jack per line; If y j, i, k+1with y j, i, knot jack per line, adds 2 by i value, is back to step S23; If y j, i, k+1with y j, i, kjack per line, adds 1 by k value, is back to step S23;
S3: using the 1st transmission channel of radar as with reference to passage, obtain the j transmission channel of radar and the phase difference α of reference channel j: Δ α j=α 1-α j; Then, according to the j transmission channel of radar and the phase difference α of reference channel j, the j transmission channel of radar is carried out to phase correction.
Feature of the present invention and further improvement are:
In step S1, first set 0 degree direction of the wave beam of radar emission, the 1st transmission channel to the N transmission channel of radar successively launching beam is oriented to the signal of 0 degree, the number of the transmission channel that N is radar; Received successively the echoed signal of each transmission channel by arbitrary receiving cable, the echoed signal of each transmission channel is sampled and orthogonal transformation, obtain the i/q signal of each transmission channel echo of radar.
In step S21, if ω jbe positioned at the second quadrant, by it around initial point dextrorotation turnback, ω jafter rotation, become complex vector located ξ j, draw rotation start-phase z j, z j=π; If ω jbe positioned at third quadrant, it be rotated counterclockwise to 180 degree, ω around initial point jafter rotation, become complex vector located ξ j, draw rotation start-phase z j, z j=-π; If ω jbe positioned at first quartile or fourth quadrant, by it around initial point 0 degree that turns clockwise, ω jafter rotation, become complex vector located ξ j, draw rotation start-phase z j, z j=0; Drawing complex vector located ξ jafterwards, calculate ξ jcorresponding plural imaginary part y 0.
In step S23, ξ j, i, k+1the imaginary part y of corresponding plural number j, i, k+1for:
y j,i,k+1=y j,i,k+d·2 -i·x j,i,k。
In step S24, setting accuracy ε gets 0.005 to 0.15.
In step S1, the initial phase of the transmitted waveform of the j transmission channel of radar is φ j;
In step S3, the j transmission channel of radar is carried out to phase correction and comprise the following steps: obtaining Δ α jafterwards, by the phase compensation Δ α of the transmitted waveform of the j transmission channel of radar j, make the phase place of the transmitted waveform of the j transmission channel of radar become φ j+ Δ α j.
Beneficial effect of the present invention is:
1) implementation complexity reduces
Traditional C ORDIC algorithm is in the time searching for, and the direction of search and the anglec of rotation are all changing.The direction of search changes with the position of current vector, in the time that vector is positioned at first quartile, rotate clockwise, when vector during in fourth quadrant by being rotated counterclockwise, step-size in search (anglec of rotation) successively decreases successively.Therefore traditional C ORDIC algorithm should be controlled sense of rotation, control again the anglec of rotation, more complicated while realization, therefore increase the difficulty that passage consistency is proofreaied and correct.But in the present invention, the one dimension individual event search adopting, according to the direction of the residing location positioning rotation of vector, once the direction of search just can not change after determining again, therefore implement comparatively simply, only need to control the anglec of rotation, iteration all will judge sense of rotation at every turn, therefore doing field of radar passage consistency timing, implementing comparatively simple.
2) computational accuracy improves:
The iterations of traditional C ORDIC algorithm is determined, therefore computational accuracy is subject to the impact of iterations, iterations revocable in the present invention, but determine in advance required precision, then continue iteration until meet default precision, therefore in the situation that calculated amount is suitable, the phase accuracy of trying to achieve with the present invention is apparently higher than the phase accuracy of trying to achieve with traditional C ORDIC algorithm, so improved the precision of radar emission channel correcting.
3) computing velocity improves:
What traditional C ORDIC algorithm adopted is bidirectional iteration method, iteration step length successively decreases successively, when make vector be rotated to another quadrant because iteration step length is oversize, in order to make vector finally infinitely approach X-axis, sense of rotation must change in good time, therefore just easily form a kind of being similar to and come and go the searching route of running, search speed is declined.What the present invention adopted is single direction rotation, once sense of rotation just no longer changes after determining, has avoided the situation of round race, and has reduced initial rotation angle degree,, in the time that the anglec of rotation need to reduce, has reduced without efficient search meanwhile, has improved computing velocity.
Brief description of the drawings
Fig. 1 is the schematic flow sheet of a kind of radar emission channel phases of the present invention bearing calibration;
The arc-tangent value that Fig. 2 is the present invention draws in emulation experiment one arc-tangent value, traditional C ORDIC algorithm draws and the contrast schematic diagram of theoretical value;
Fig. 3 is the partial enlarged drawing of Fig. 2;
Fig. 4 is traditional C ORDIC algorithm and the error of calculation contrast schematic diagram of the present invention in the time calculating arc-tangent value in emulation experiment one;
The arc-tangent value that Fig. 5 is the present invention draws in emulation experiment two arc-tangent value, traditional C ORDIC algorithm draws and the contrast schematic diagram of theoretical value;
Fig. 6 is the partial enlarged drawing of Fig. 5;
Fig. 7 is traditional C ORDIC algorithm and the error of calculation contrast schematic diagram of the present invention in the time calculating arc-tangent value in emulation experiment two.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described:
With reference to Fig. 1, it is the schematic flow sheet of a kind of radar emission channel phases of the present invention bearing calibration.This radar emission channel phases bearing calibration comprises the following steps:
S1: utilize the i/q signal of the echo of i/q signal to the N transmission channel of the 1st transmission channel echo of arbitrary receiving cable acquisition radar of radar, the number of the transmission channel that N is radar.Be described as follows:
First the 0 degree direction of setting the wave beam of radar emission, the 1st transmission channel to the N transmission channel of radar successively launching beam is oriented to the signal of 0 degree, the number of the transmission channel that N is radar; Received successively the echoed signal of each transmission channel by arbitrary receiving cable (special correcting circuit passage), the echoed signal of each transmission channel is sampled and orthogonal transformation, obtain the i/q signal of the echo of i/q signal to the N transmission channel of the 1st transmission channel echo of radar.In the embodiment of the present invention, the initial phase of the transmitted waveform of the j transmission channel of radar is φ j, j gets 1 to N.
S2: final phase place corresponding to each transmission channel of determining radar; Final phase place corresponding to each transmission channel of determining radar comprises the following steps:
S21: in plane right-angle coordinate, according to the phase place of the i/q signal of the echo of j transmission channel and amplitude, by corresponding complex vector located ω for the i/q signal of the echo of j transmission channel jrepresent; According to ω jplace quadrant, by ω jbe rotated ω around initial point jafter rotation, become complex vector located ξ j, draw rotation start-phase z jand ξ jcorresponding plural imaginary part y 0.Be described as follows:
In this plane right-angle coordinate, taking X-axis forward as 0 degree direction, taking X-axis negative sense as 180 degree directions, taking Y-axis forward as 90 degree directions, taking Y-axis negative sense as 270 degree directions.If the i/q signal of the echo of j transmission channel comprises real component I and imaginary part component Q, real component I and complex vector located ω jhorizontal ordinate correspondence, imaginary part component Q and complex vector located ω jordinate correspondence.If I<0, Q>0, ω jin the second quadrant; If I<0, Q<0, ω jin third quadrant; If I>0, Q<0, ω jin fourth quadrant; If I>0, Q>0, ω jin fourth quadrant.
If ω jbe positioned at the second quadrant, by it around initial point dextrorotation turnback, ω jafter rotation, become complex vector located ξ j, draw rotation start-phase z j, z j=π.
If ω jbe positioned at third quadrant, it be rotated counterclockwise to 180 degree, ω around initial point jafter rotation, become complex vector located ξ j, draw rotation start-phase z j, z j=-π.
If ω jbe positioned at first quartile or fourth quadrant, by it around initial point 0 degree (namely not rotating) that turns clockwise, ω jafter rotation, become complex vector located ξ j, draw rotation start-phase z j, z j=0; Drawing complex vector located ξ jafterwards, calculate ξ jcorresponding real x 0with imaginary part y 0.
S22: parameters d, if ξ jbe positioned at first quartile, make d=-1; Otherwise, make d=1; Variable i and k are set, and i gets odd number, and the initial value of i is 1; K gets 1,2, and 3 In the time of k=1, establish complex vector located ξ j, i, k=ξ j, and establish z j, i, k=z j; ξ j, i, kcorresponding real x j, i, kfor x 0, ξ j, i, kcorresponding plural imaginary part y j, i, kfor y 0.
S23: by ξ j, i, kbe rotated around initial point, its anglec of rotation is θ i, θ i=arctan (2 -i), sense of rotation is by the decision of d value, and in the time of d=-1, sense of rotation is clockwise direction; In the time of d=1, sense of rotation is counterclockwise; ξ j, i, kafter rotation, become complex vector located ξ j, i, k+1; Draw ξ j, i, k+1corresponding rotatable phase z j, i, k+1, z j, i, k+1=z j, i, k-d θ i; Draw ξ j, i, k+1corresponding real x j, i, k+1with imaginary part y j, i, k+1:
x j,i,k+1=x j,i,k-d·2 -i·y j,i,k，y j,i,k+1=y j,i,k+d·2 -i·x j,i,k。
S24: if | y j, i, k+1| (y j, i, k+1absolute value) be less than setting accuracy ε, obtain the final phase place α that j transmission channel is corresponding j, α j=z j, i, k+1.
If | y j, i, k+1| be more than or equal to setting accuracy ε, judge y j, i, k+1with y j, i, kwhether jack per line (if y j, i, k+1with y j, i, kbe positive negative two numbers, y j, i, k+1with y j, i, kjack per line, otherwise y j, i, k+1with y j, i, knot jack per line), if y j, i, k+1with y j, i, knot jack per line, illustrate that the last time anglec of rotation is too large, complex vector located another quadrant that has been rotated to (is rotated to fourth quadrant from first quartile, or rotate to first quartile from fourth quadrant), for realizing one direction search, avoid coming and going the searching route of race formula, therefore judge that this rotation is invalid, while being again rotated, the anglec of rotation need to be diminished, therefore i value is added 2(by i+2 assignment to i), be back to step S23; Because i value becomes large, anglec of rotation θ idiminish.
If y j, i, k+1with y j, i, kjack per line, illustrates that the last time anglec of rotation is suitable, ξ j, i, k+1and ξ j, i, kin same quadrant, now, can continue to be rotated according to equidirectional, while needing to continue rotation, the anglec of rotation is identical with the anglec of rotation of rotation last time.Now, by k value add 1(by k+1 assignment to k), be back to step S23.In the embodiment of the present invention, setting accuracy ε gets 0.005 to 0.15.
Because j gets 1 to N, therefore, after step S2, just determine the final phase place corresponding to each transmission channel of radar.
S3: using the 1st transmission channel of radar as with reference to passage, obtain the j transmission channel of radar and the phase difference α of reference channel j: Δ α j=α 1-α j; Then, according to the j transmission channel of radar and the phase difference α of reference channel j, the j transmission channel of radar is carried out to phase correction.Particularly, the j transmission channel of radar being carried out to phase correction comprises the following steps: obtaining Δ α jafterwards, by the phase compensation Δ α of the transmitted waveform of the j transmission channel of radar j, make the phase place of the transmitted waveform of the j transmission channel of radar become φ j+ Δ α j.
It should be noted that, the present invention also can be applicable to radar receiving cable phase correction, in the time carrying out radar receiving cable phase correction, be oriented to the signal (needing equally 0 degree direction of the wave beam of setting in advance radar emission) of 0 degree by a transmission channel launching beam, then there is each receiving cable of radar to receive the signal of this transmission channel, each receiving cable, by the signal receiving is sampled and orthogonal transformation, obtains corresponding i/q signal; After obtaining i/q signal, corresponding receiving cable phase correction and radar emission channel phases of the present invention bearing calibration are similar, no longer repeat at this.
Generally speaking, technical thought of the present invention is: the each passage of radar successively launching beam is oriented to the signal of 0 degree, receives successively each transmission channel echo by receiving cable, and the echo receiving is sampled and orthogonal transformation, obtains i/q signal; Judge corresponding vectorial position according to i/q signal again, will rotate to first, fourth quadrant in the vector of second and third quadrant; Then determine the direction of rotation according to the residing quadrant of rotating vector, the angle of each rotation and vectorial reference position are to judge according to the rotation result of last time, even after rotation last time, vector has been rotated to another quadrant, the reference position of this rotating vector is identical with the vectorial reference position of rotation last time, and the anglec of rotation reduces; Otherwise this rotating vector reference position is updated to the position of the rear vector of rotation last time, the anglec of rotation remains unchanged, and so circulation, until meet setting accuracy; Obtain respectively afterwards the final phase place of each transmission channel; Finally select a wherein transmitting passage as with reference to passage, to obtain the phase differential of each transmission channel and reference channel, when radar signal, each transmission channel is compensated respectively to a phase differential, make each transmission channel realize phase equalization.
Effect of the present invention further illustrates by following emulation experiment:
1) simulated conditions:
Get at random the angle of 20-180 °～180 °, be distributed in four quadrants; Emulation experiment is divided into two large classes, and accuracy requirement is higher and accuracy requirement is lower; When accuracy requirement is higher, setting accuracy of the present invention ε=0.01, the iterations of traditional C ORDIC algorithm is 7 times; When accuracy requirement is lower, setting accuracy of the present invention ε=0.1, the iterations of traditional C ORDIC algorithm is 4 times.In following emulation experiment, use the calculating of arc-tangent value as the judgment criteria of the present invention and traditional C ORDIC algorithm.
2) emulation content
Emulation experiment one: in the time that computational accuracy is had relatively high expectations, use this method and traditional C ORDIC algorithm to carry out the correction of radar emission channel phases.In emulation experiment one, get at random 20 angles between-180 °～180 °, be randomly distributed in four quadrants, precision parameter of the present invention ε=0.01, the iterations of traditional C ORDIC method is 7 times.With reference to Fig. 2, be the present invention draws in emulation experiment one arc-tangent value, arc-tangent value that traditional C ORDIC algorithm draws and the contrast schematic diagram of theoretical value.As can be seen from Figure 2, now the result of calculation of two kinds of methods has all approached theoretical value very much.Fig. 3 is the partial enlarged drawing of Fig. 2, and from then on Fig. 3 can find out, although now the result of calculation of two kinds of methods has all approached theoretical value very much, the present invention more approaches theoretical value than traditional C ORDIC algorithm.With reference to Fig. 4, it is the error of calculation contrast schematic diagram in the time calculating arc-tangent value of traditional C ORDIC algorithm and the present invention in emulation experiment one.As can be seen from Figure 4, the error of calculation of the present invention is less than the error of calculation of traditional C ORDIC algorithm, when actual computation, average error of the present invention is 0.2555 °, total calculation times is 132 times (be often called once and calculate to step S23), the average error of traditional C ORDIC algorithm is 0.5166 °, total calculation times is 140 times, therefore in the time that accuracy requirement is higher, the present invention not only computational accuracy higher than traditional C ORDIC algorithm, and calculation times is also less than traditional C ORDIC algorithm, therefore solving high precision anyway when function, with the present invention calculate precision and speed be all better than traditional C ORDIC algorithm.
Emulation experiment two, when computational accuracy requirement is lower, is used this method and traditional C ORDIC algorithm to carry out the correction of radar emission channel phases.In emulation experiment two, get at random 20 angles between-180 °～180 °, be randomly distributed in four quadrants, precision parameter of the present invention ε=0.1, the iterations of traditional C ORDIC method is 4 times.With reference to Fig. 5, be the present invention draws in emulation experiment two arc-tangent value, arc-tangent value that traditional C ORDIC algorithm draws and the contrast schematic diagram of theoretical value; With reference to Fig. 6, it is the partial enlarged drawing of Fig. 5.As can be seen from Figure 5, the result of calculation of two kinds of methods and theoretical value have certain error, then in conjunction with Fig. 6 observe, can find result of calculation of the present invention with respect to traditional C ORDIC algorithm more close to theoretical value.With reference to Fig. 7, it is the error of calculation contrast schematic diagram in the time calculating arc-tangent value of traditional C ORDIC algorithm and the present invention in emulation experiment two.As can be seen from Figure 7, the error of calculation entirety of traditional C ORDIC algorithm is greater than the error of calculation of the present invention; When actual computation, average error of the present invention is 1.5069 °, total calculation times is 66 times (every execution step S23 is called once calculating), the average error of traditional C ORDIC algorithm is 4.1628 °, and total calculation times is 80 times, therefore in the time that accuracy requirement is not high, the present invention not only computational accuracy higher than traditional C ORDIC algorithm, and calculation times is also less than traditional C ORDIC algorithm, therefore in the time of phase place rough estimate, be all better than traditional C ORDIC algorithm by precision and speed that the present invention calculates.
In sum, when the present invention solves echoed signal phase place for the correction of field of radar passage consistency, no matter be requirement Exact Solution or rough estimate, computational accuracy of the present invention and computing velocity are all higher than the result of calculation of using traditional C ORDIC algorithm, and owing to having adopted unidirectional search, be that in search procedure, a needs is controlled the step-length of search and do not need to control sense of rotation, therefore implementation complexity of the present invention is also less than traditional C ORDIC algorithm, so the present invention can effectively improve the effect of the transmission channel phase correction of radar.
Obviously, those skilled in the art can carry out various changes and modification and not depart from the spirit and scope of the present invention the present invention.Like this, if these amendments of the present invention and within modification belongs to the scope of the claims in the present invention and equivalent technologies thereof, the present invention is also intended to comprise these changes and modification interior.

Claims (6)

1. the bearing calibration of radar emission channel phases, is characterized in that, comprises the following steps:
S1: utilize the i/q signal of the echo of i/q signal to the N transmission channel of the 1st transmission channel echo of arbitrary receiving cable acquisition radar of radar, the number of the transmission channel that N is radar;
S2: final phase place corresponding to each transmission channel of determining radar; The final phase place corresponding to each transmission channel of described definite radar comprises the following steps:
S21: in plane right-angle coordinate, according to the phase place of the i/q signal of the echo of j transmission channel and amplitude, by corresponding complex vector located ω for the i/q signal of the echo of j transmission channel jrepresent, j gets 1 to N; According to ω jplace quadrant, by ω jbe rotated ω around initial point jafter rotation, become complex vector located ξ j, draw rotation start-phase z jand ξ jcorresponding plural imaginary part y 0;
S22: parameters d, if ξ jbe positioned at first quartile, make d=-1; Otherwise, make d=1; Variable i and k are set, and i gets odd number, and the initial value of i is 1; K gets 1,2, and 3 In the time of k=1, establish complex vector located ξ j, i, k=ξ j, and establish z j, i, k=z j; ξ j, i, kcorresponding plural imaginary part y j, i, kfor y 0;
S23: by ξ j, i, kbe rotated around initial point, its anglec of rotation is θ i, θ i=arctan (2 -i), sense of rotation is by the decision of d value, and in the time of d=-1, sense of rotation is clockwise direction; In the time of d=1, sense of rotation is counterclockwise; ξ j, i, kafter rotation, become complex vector located ξ j, i, k+1; Draw ξ j, i, k+1corresponding rotatable phase z j, i, k+1, z j, i, k+1=z j, i, k-d θ i; Draw ξ j, i, k+1corresponding plural imaginary part y j, i, k+1;
S24: if | y j, i, k+1| be less than setting accuracy ε, obtain the final phase place α that j transmission channel is corresponding j, α j=z j, i, k+1; If | y j, i, k+1| be more than or equal to setting accuracy ε, judge y j, i, k+1with y j, i, kwhether jack per line; If y j, i, k+1with y j, i, knot jack per line, adds 2 by i value, is back to step S23; If y j, i, k+1with y j, i, kjack per line, adds 1 by k value, is back to step S23;
S3: using the 1st transmission channel of radar as with reference to passage, obtain the j transmission channel of radar and the phase difference α of reference channel j: Δ α j=α 1-α j; Then, according to the j transmission channel of radar and the phase difference α of reference channel j, the j transmission channel of radar is carried out to phase correction.
2. a kind of radar emission channel phases as claimed in claim 1 bearing calibration, it is characterized in that, in step S1, first set 0 degree direction of the wave beam of radar emission, the 1st transmission channel to the N transmission channel of radar successively launching beam is oriented to the signal of 0 degree, the number of the transmission channel that N is radar; Received successively the echoed signal of each transmission channel by arbitrary receiving cable, the echoed signal of each transmission channel is sampled and orthogonal transformation, obtain the i/q signal of each transmission channel echo of radar.
3. a kind of radar emission channel phases as claimed in claim 1 bearing calibration, is characterized in that, in step S21, if ω jbe positioned at the second quadrant, by it around initial point dextrorotation turnback, ω jafter rotation, become complex vector located ξ j, draw rotation start-phase z j, z j=π; If ω jbe positioned at third quadrant, it be rotated counterclockwise to 180 degree, ω around initial point jafter rotation, become complex vector located ξ j, draw rotation start-phase z j, z j=-π; If ω jbe positioned at first quartile or fourth quadrant, by it around initial point 0 degree that turns clockwise, ω jafter rotation, become complex vector located ξ j, draw rotation start-phase z j, z j=0; Drawing complex vector located ξ jafterwards, calculate ξ jcorresponding plural imaginary part y 0.
4. a kind of radar emission channel phases as claimed in claim 1 bearing calibration, is characterized in that, in step S23, and ξ j, i, k+1the imaginary part y of corresponding plural number j, i, k+1for:
y j,i,k+1=y j,i,k+d·2 -i·x j,i,k。
5. a kind of radar emission channel phases as claimed in claim 1 bearing calibration, is characterized in that, in step S24, setting accuracy ε gets 0.005 to 0.15.
6. a kind of radar emission channel phases as claimed in claim 1 bearing calibration, is characterized in that, in step S1, the initial phase of the transmitted waveform of the j transmission channel of radar is φ j;
In step S3, the j transmission channel of radar is carried out to phase correction and comprise the following steps: obtaining Δ α jafterwards, by the phase compensation Δ α of the transmitted waveform of the j transmission channel of radar j, make the phase place of the transmitted waveform of the j transmission channel of radar become φ j+ Δ α j.
CN201410105470.4A 2014-03-20 A kind of radar emission channel phases bearing calibration Expired - Fee Related CN104199004B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410105470.4A CN104199004B (en) 2014-03-20 A kind of radar emission channel phases bearing calibration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410105470.4A CN104199004B (en) 2014-03-20 A kind of radar emission channel phases bearing calibration

Publications (2)

Publication Number Publication Date
CN104199004A true CN104199004A (en) 2014-12-10
CN104199004B CN104199004B (en) 2016-11-30

Family

ID=

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107229046A (en) * 2016-03-25 2017-10-03 松下电器产业株式会社 Radar installations and radar method
CN109032562A (en) * 2018-06-14 2018-12-18 浙江大学 The cordic algorithm device and algorithm of low time delay high frequency single direction rotation
CN109076040A (en) * 2016-06-17 2018-12-21 Ntt 电子株式会社 Phase compensation device, phase compensating method and communication device
CN110289920A (en) * 2019-03-22 2019-09-27 北京理工雷科电子信息技术有限公司 A kind of phase alignment of the multi-channel wide with receive-transmit system
CN114137474A (en) * 2021-11-30 2022-03-04 天津光电通信技术有限公司 Rotation method and device suitable for broadband direction finding receiver

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0959522A1 (en) * 1998-05-19 1999-11-24 Toyota Jidosha Kabushiki Kaisha Method for determining phase correction values in radar apparatus
CN101957444A (en) * 2010-09-30 2011-01-26 中国船舶重工集团公司第七二三研究所 Multichannel radar amplitude and phase automatic correcting method and device
CN103399307A (en) * 2013-08-05 2013-11-20 中国科学院电子学研究所 Method and device for correcting channel errors

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0959522A1 (en) * 1998-05-19 1999-11-24 Toyota Jidosha Kabushiki Kaisha Method for determining phase correction values in radar apparatus
CN101957444A (en) * 2010-09-30 2011-01-26 中国船舶重工集团公司第七二三研究所 Multichannel radar amplitude and phase automatic correcting method and device
CN103399307A (en) * 2013-08-05 2013-11-20 中国科学院电子学研究所 Method and device for correcting channel errors

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ERAN FISHLER ET AL.: "MIMO RADAR:AN IDEA WHOSE TIME HAS COME", 《IEEE NATIONAL RADAR CONFERENCE》 *
N.OZBEN ONHON ET AL.: "A Sparsity-Driven Approach for Joint SAR Imaging and Phase Error Correction", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 *
党晓方等: "双基地MIMO雷达发射阵列幅相误差的旋转阵列校正法", 《***工程与电子技术》 *
杨明磊等: "多载频MIMO雷达幅相误差校正", 《***工程与电子技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107229046A (en) * 2016-03-25 2017-10-03 松下电器产业株式会社 Radar installations and radar method
CN109076040A (en) * 2016-06-17 2018-12-21 Ntt 电子株式会社 Phase compensation device, phase compensating method and communication device
CN109032562A (en) * 2018-06-14 2018-12-18 浙江大学 The cordic algorithm device and algorithm of low time delay high frequency single direction rotation
CN109032562B (en) * 2018-06-14 2020-12-08 浙江大学 CORDIC algorithm device and algorithm with low time delay and high frequency unidirectional rotation
CN110289920A (en) * 2019-03-22 2019-09-27 北京理工雷科电子信息技术有限公司 A kind of phase alignment of the multi-channel wide with receive-transmit system
CN114137474A (en) * 2021-11-30 2022-03-04 天津光电通信技术有限公司 Rotation method and device suitable for broadband direction finding receiver

Similar Documents

Publication Publication Date Title
CN109298388B (en) Direct over-the-horizon target geographic coordinate estimation method based on azimuth information
CN105738865B (en) The direct localization method of multiple target under signal waveform known conditions
CN103135083A (en) Electromagnetic vector sensor array amplitude and phase error self-correcting method based on array rotation
TWI663849B (en) Beamforming calibration system and method
CN105158759B (en) HRWS SAR channel phases deviation correction methods based on clutter cancellation
US9300518B2 (en) Pulse radar device and control method therefor
CN105487063A (en) Direct positioning method based on external radiation source time delay and Doppler frequency
WO2017181763A1 (en) Carrier phase recovery method and apparatus, and storage medium
CN105262527B (en) Antenna axis calibration device and method
JP6178245B2 (en) Synthetic aperture radar apparatus and image processing method thereof
CN103760522A (en) Method and system used for time difference estimation and multi-station clock error correction
CN105738891A (en) Method for tracking weak maneuvering target angle through airborne digital array radar
CN107085202A (en) A kind of Capon robust adaptive beamforming method and system for mutual coupling error
CN103780294B (en) A kind of maximum signal to noise ratio output phased array antenna method of weighting
US20160259036A1 (en) Radar apparatus and signal generating apparatus
CN110350990B (en) Phased array network calibration method, device, equipment and storage medium
CN104199004A (en) Phase correction method for radar transmission channels
CN109412984B (en) Aitken acceleration method-based blind signal-to-noise ratio estimation method in multi-antenna scene
CN109917384B (en) Medium and long distance frequency domain fast imaging method and device for cantilever scanning foundation SAR
CN104199004B (en) A kind of radar emission channel phases bearing calibration
CN108445445A (en) A kind of distribution passive location method and device
CN210243826U (en) Radar multichannel signal preprocessing device and pulse compression unit thereof
CN110943949B (en) Broadband angle tracking method based on SC-FDE transmission system
CN107395547A (en) Compensation method and system under the conditions of a kind of IQ branch roads imbalance
JP6373408B2 (en) Channel equalization tracking apparatus, method and receiver

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20161130