CN115032600A - Circular array secondary radar sectional type weight coefficient processing method based on matrix array - Google Patents
Circular array secondary radar sectional type weight coefficient processing method based on matrix array Download PDFInfo
- Publication number
- CN115032600A CN115032600A CN202210953623.5A CN202210953623A CN115032600A CN 115032600 A CN115032600 A CN 115032600A CN 202210953623 A CN202210953623 A CN 202210953623A CN 115032600 A CN115032600 A CN 115032600A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- matrix
- calibration
- phase
- value
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
The invention discloses a matrix array-based circular array secondary radar sectional weight coefficient processing method, which comprises the following steps of: s1, receiving calibration data of the circular array secondary radar calibration channel and judging the effectiveness of the calibration data; s2, forming adaptive and linear cyclic calculation of weight coefficients for the calibration data through iterative verification to obtain a temporary matrix B 1 (ii) a S3, obtaining an amplitude-phase error table, calculating the circular mapping of matrix arrays and a temporary matrix B 1 Performing one-to-one correspondence on the matrix array B obtained by circular mapping to generate a matrix C; and S4, performing orthogonal IQ signal calculation on the matrix C to generate a final value weight coefficient matrix D for calling a wave control DBF. The matrix array introduced by the invention can quickly and effectively map the N array element wave beam forming standard weighting coefficient value which is required to be changed correspondingly, thereby saving time, saving labor, being accurate and quickly and effectively ensuring the detection performance of the circular array secondary radar.
Description
Technical Field
The invention belongs to the technical field of secondary radar calibration, and particularly relates to a matrix array-based circular array secondary radar sectional weight coefficient processing method.
Background
The circular array secondary radar is still a phased array system essentially, the amplitude-phase characteristics of each channel in the system can be changed in the processes of installation, use, daily maintenance, overhaul and the like of the circular array secondary radar, at the moment, T/R calibration needs to be carried out on the system or DBF (direct digital filter) of a receiving and transmitting channel needs to be compensated according to coefficients, and a formula is introduced
Where Y is the output signal, X is the vector of the input signals, W is the vector of the weighting coefficients, and H is the conjugate matrix. For digital beamforming, the receiver is usually composed of two data channels orthogonal to each other, one being the real part of the signal and the other being the imaginary part of the signal, the signal and the weighting coefficients now being represented in complex form. And calculating a coupling amplitude-phase error table (including transmitting (1030 Mhz) amplitude and phase value of N array elements of an antenna and receiving (1090 Mhz) amplitude and phase value of N array elements) transmitted and received by each channel according to transmitting amplitude and phase value of all N channels and receiving amplitude and phase value of N channels acquired by AD sampling from a TR calibration channel, and calculating a DBF weighting coefficient table. The invention provides a matrix array-based sectional weight coefficient processing method for a circular array secondary radar, according to the characteristics of a circular array, each array element participates in a group array of a plurality of wave positions, a created matrix A represents a coupling amplitude-phase error table transmitted and received by each channel, N array element wave beams are introduced to form a standard weighting coefficient table, a matrix B is created to represent the same, and the amplitude-phase errors of each channel can be quickly compensated through mapping and superposition calculation of the coupling amplitude-phase error tables of the channels and the standard weighting coefficient table, so that the receiving and transmitting wave beams are stably formed and point to the right direction.
Aiming at whether the weight coefficient is variable, digital beam forming can be divided into common digital beam forming and Adaptive Digital Beam Forming (ADBF), and the weight coefficient in the common digital beam forming process is fixed and unchanged as the traditional phased array antenna beam forming, while the weight coefficient matrix B introduced by the invention can be adjusted along with the current working environment of an antenna and a secondary radar system to form a beam with optimal performance.
The method mainly includes that amplitude and phase errors in a phased array system mainly come from antenna phase center error measurement (negligible), a TR channel and an AD/DA channel, and mainly discussed here is that under the condition of system wireless calibration, channel amplitude/phase errors are obtained, and a coupling amplitude and phase error table is established to be used for generating DBF weighting coefficients. When the traditional phased array secondary radar is calibrated, the receiving and transmitting amplitude phases of a TR channel and an antenna array element can be compensated and calibrated only in a darkroom environment, and the calibration is limited by environmental conditions and is not flexible enough; and in the processes of installation, use, daily maintenance, overhaul and the like of the circular array phased array secondary radar, when the secondary deviation of the amplitude and the phase occurs to part of radio frequency channels due to external or internal factors such as service life or temperature, the secondary deviation can not be calibrated and compensated again in an external field, so that the performance of the phased array secondary radar is reduced.
Disclosure of Invention
Aiming at the defects in the prior art, the matrix array-based sectional weight coefficient processing method for the circular array secondary radar solves the problems that the conventional phased array secondary radar is limited by environmental conditions and is not flexible when being calibrated; in the same process of installation, use, daily maintenance, overhaul and the like of the circular array phased array secondary radar, when the radio frequency channel deviates again in amplitude and phase due to external or internal factors such as service life or temperature, the performance of the phased array secondary radar can not be calibrated and compensated again, and the performance of the phased array secondary radar is reduced.
In order to achieve the purpose of the invention, the invention adopts the technical scheme that: a circular array secondary radar sectional type weight coefficient processing method based on a matrix array comprises the following steps:
s1, receiving calibration data of the circular array secondary radar calibration channel and judging the effectiveness of the calibration data;
s2, forming adaptive and linear circulation calculation of weight coefficient for the calibration data through iterative verification to obtain a temporary matrix B 1 ；
S3, obtaining the amplitude-phase error table, calculating the circular mapping of the matrix array, and calculating the temporary matrix B 1 Performing one-to-one correspondence on the matrix array B obtained by cyclic mapping to generate a matrix C;
and S4, performing orthogonal IQ signal calculation on the matrix C to generate a final value weight coefficient matrix D for calling a wave control DBF.
Further: the specific steps of step S1 are:
s11, receiving amplitude and phase sampling values of the circular array secondary radar calibration channel in a receiving and transmitting deviation table, and respectively storing a calibration receiving amplitude value, a calibration receiving phase value, a calibration transmitting amplitude value and a calibration transmitting phase value according to the channels;
s12, judging whether the current calibration is the first calibration of the current process, if so, sending a calibration request, returning to the step S11, and otherwise, entering the step S13;
and S13, respectively carrying out data validity judgment on the calibration receiving amplitude value and the calibration transmitting amplitude value, and the calibration receiving phase value and the calibration transmitting phase value, if the validity conditions are met, entering the step S2, and if the validity conditions are not met, returning to the step S11.
Further, the method comprises the following steps: the validity judgment is as follows: the amplitude value needs to respectively meet the conditions that the receiving amplitude value threshold and the transmitting amplitude value threshold are met, and the sampling difference of the calibration receiving amplitude value and the calibration transmitting amplitude value of the same channel between two times of continuous calibration is within 1 dB; the phase value needs to satisfy the condition that the sampling difference between the calibration receiving phase value and the calibration transmitting phase value of the same channel between two successive calibrations is within 10 degrees.
Further, the method comprises the following steps: the specific steps of step S2 are:
s21, judging whether the calibration receiving amplitude value and the calibration transmitting amplitude value are normal, if all the channels are normal, entering a step S22, otherwise, checking the corresponding channel fault, and ending the process;
s22, starting normal air-to-air work, through analyzing the response data condition of each wave position secondary radar, if the number of each wave position response data is equal to the number of inquiry triggers, entering step S23, if the number of each wave position response data is far less than the number of inquiry triggers, entering step S26;
s23, when the air-to-air work is normal, the sum-difference amplitude of each wave position is normal, namely the sum-amplitude difference amplitude is high, the standard weight coefficient does not need to be modified, the default standard weight coefficient is adopted to carry out linear circulation calculation to obtain a temporary matrix B1, the step S3 is carried out, and if the sum-amplitude difference amplitude of a certain wave position beam _ x is low, the step S24 is carried out;
s24, adjusting the wave position beam _ x central position to be closed corresponding to the channel, and performing linear circulation calculation to obtain a temporary matrix B 1 Entering step S3, starting normal air-to-air operation again, if the secondary radar response data and the difference amplitude of the wave position beam _ x are normal, restoring the channel corresponding to the central position of the wave position beam _ x to be started, adjusting the value of the receiving standard and/or difference phase weight coefficient at the corresponding position through an interface, adjusting by taking 90 degrees or 45 degrees as steps, entering step S27, and if the secondary radar response data and the difference amplitude are still abnormal, entering step S25;
s25, adjusting wave position beam _ x times of closing of the corresponding channel of the central position, and performing linear circulation calculation to obtain a temporary matrix B again 1 Step S3, normal air work is started again, if the secondary radar response data and the difference amplitude are normal, the wave position beam _ x corresponding channel is restored to be started, the corresponding position receiving standard and/or difference phase weight coefficient value is adjusted through the interface, adjustment is carried out by taking 90 degrees as stepping, step S27 is entered, if the secondary radar response data and the difference amplitude are still abnormal, the corresponding channel of the beam _ x wave position is checked through an instrument, and the process is ended;
s26, transmitting standard phase weight coefficient value through interface adjusting wave position beam _ x central position, adjusting by taking 90 degrees as step, and then performing linear circulation calculation to obtain a temporary matrix B 1 The step S3 is entered, normal air-to-air work is started again, and if the wave position beam _ x response data number is equal to the inquiry trigger number, the step S23 is returned; if the number of the wave position response data is still far less than the number of the inquiry triggers, the corresponding channel of the wave position beam _ x is checked through the instrument, and the process is ended;
s27, linear circulation calculation is carried out to obtain a temporary matrix B again 1 And starting normal air work again, observing whether the inquiry detection performance of the wave position beam _ x is normal, if so, entering the step S3, and if not, checking a corresponding channel of the wave position beam _ x through an instrument, and ending the process.
Further: the wave position participates in DBF by 14 channels, and when in transmission: standard amplitude weighting factor of [13,11,8,7,4,2,0,0,2,4,7,8,11,13], standard phase weighting factor of [180,2,168,310,63,141,180,0,321,243,130,348,182,0], receive: the standard and amplitude weighting coefficients are [13,11,8,7,4,2,0,0,2,4,7,8,11,13], the standard deviation amplitude weighting coefficients are [13,11,8,7,4,2,0,0,2,4,7,8,11,13], the standard and phase weighting coefficients are [180,2,168,310,63,141,180,0,321,243,130,348,182,0], and the standard deviation phase weighting coefficients are [0,182,348,130,243,321,0,0,321,243,130,348,182,0 ].
Further, the method comprises the following steps: the specific steps of step S3 are:
s31, calculating a receiving and/difference amplitude error table, namely finding the rxMeasure _ Amp [ minimum ] with the lowest amplitude in N radio frequency channels as a reference amplitude value, and calculating the error of the amplitude of each channel relative to the reference amplitude value, and recording the error as rxMeasure _ AmpError [ i ];
s32, calculating a receiving and/or difference Phase error table, taking rxMeasure _ Phase [0] as a reference Phase value, calculating the error of each channel Phase relative to the reference Phase value, and marking as rxMeasure _ PhaseError [ i ];
s33, calculating a transmission and/or difference amplitude error table and a transmission and/or difference Phase error table, taking the txMeasure _ Amp [0] and the txMeasure _ Phase [0] as reference amplitude values and Phase values, and calculating the errors of the amplitudes of all channels relative to the reference amplitude values and the phases relative to the reference Phase values, which are marked as txMeasure _ Amperer [ i ] and txMeasure _ PhaseError [ i ];
s34, circularly mapping the variable standard weight coefficients to obtain 6 original temporary matrixes B of 14 multiplied by 1 1 Circularly mapping to N wave bits to obtain a matrix array B of 30 multiplied by 14 multiplied by 6;
s35, carrying out wave position cyclic mapping one-to-one corresponding calculation on a matrix A and a variable scale matrix array B by taking the rxMeasure _ Amperer [ i ], the rxMeasure _ PhaseError [ i ], the rxMeasure _ Amperer [ i ] and the txMeasure _ PhaseError [ i ] as basic data, and generating a DBF weight coefficient matrix C.
Further: the specific steps of step S4 are:
s41, obtaining a linear amplitude according to the amplitude/phase value corresponding to the DBF weight coefficient matrix C;
and S42, obtaining the orthogonal I/Q value of the signal through linear amplitude calculation, and generating a final value weight coefficient matrix D for calling a wave control DBF.
The invention has the beneficial effects that: the invention provides a matrix array-based sectional weighting coefficient processing method for a circular array secondary radar, which is characterized in that an NxL matrix A is created by adopting a matrix array to represent a coupling amplitude-phase error table transmitted and received by each channel, an N array element wave beam forming standard weighting coefficient table is introduced at the same time, an NxMxL three-dimensional matrix B is created to represent the same, the amplitude-phase error of each channel can be quickly compensated through cyclic mapping of the matrix A and the matrix B, and a circular array wave position weighting coefficient is generated. The matrix B of the N-element wave beam forming standard weighting coefficient table introduced by the invention can be quickly modified in a reading and writing mode of the interface software chart, so that the real-time compensation and calibration of the channel corresponding to the N-element antenna array are realized, the introduced matrix array can quickly and effectively map the N-element wave beam forming standard weighting coefficient value which corresponds to the N-element wave beam to be changed, the time and labor are saved, the accuracy is realized, and the detection performance of the circular array secondary radar can be quickly and effectively ensured.
With the application of the invention on the circular array secondary radar equipment, when the detection performance of each wave position is reduced, the detection performance of the circular array secondary radar is ensured by the sectional type weight coefficient read-write processing method of the matrix array, and powerful guarantee is provided for ensuring the airspace flight safety.
Drawings
FIG. 1 is a flow chart of the sectional type weight coefficient processing of a circular array secondary radar based on a matrix array according to the present invention;
FIG. 2 is a schematic diagram of a matrix array-based sectional weight coefficient processing calculation for a circular array secondary radar according to the present invention;
FIG. 3 is a flow chart illustrating the design of step S1 relating to the calibration data reception and correctness determination;
FIG. 4 is a block diagram of step S2 of the present invention of the temporary matrix B for variable standard beamforming 1 Processing the associated design flow diagram;
FIG. 5 is a flowchart of the related design of the amplitude-phase DBF weight coefficient matrix C in step S3 according to the present invention;
FIG. 6 is a flow chart of the design of step S4 DBF IQ signal calculation correlation.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
The invention provides a matrix array-based sectional weight coefficient processing method for a circular array secondary radar, which adopts an NxL matrix A established by a matrix array to represent a coupling amplitude-phase error table transmitted and received by each channel, wherein M is the number of array elements participating in array grouping, and the value of L is 6, and represents 6 deviation tables of a receiving and amplitude deviation table, a receiving difference amplitude deviation table, a receiving and phase deviation table, a receiving difference phase deviation table, a transmitting amplitude deviation table and a transmitting phase deviation table. And simultaneously, introducing an N array element beam forming standard weighting coefficient table, similarly creating an NxMxL three-dimensional matrix B for representation, performing superposition calculation on a coupling amplitude-phase error table of the channel and the standard weighting coefficient table, namely performing cyclic mapping on the matrix A and the matrix B to quickly compensate the amplitude-phase error of each channel to obtain a three-dimensional matrix C of a final value of the weighting coefficient, and performing orthogonal calculation on the matrix C to generate an I/Q coefficient matrix D for DBF (direct guided wave) so as to form stable receiving and transmitting beams and point correctly. When the matrix A, namely a coupling amplitude-phase error table of a channel, changes, the matrix B of the N-element beam forming standard weighting coefficient table introduced by the invention can be rapidly modified in a reading and writing mode of an interface software chart, so that the channel corresponding to the N-element antenna array can be compensated and calibrated in real time, and the introduced matrix array can rapidly and effectively map the N-element beam forming standard weighting coefficient value which corresponds to the N-element antenna array and needs to be modified, thereby saving time, labor and accuracy. The invention is explained by taking 30-channel circular array secondary radar as an example, wherein 30 wave bits are provided, each wave bit is participated in DBF by 14 channels, and the calculation process of the transmitting/receiving amplitude-phase weighting coefficient is shown in fig. 1 and fig. 2.
The invention is further described below with reference to the accompanying drawings.
Step 1: and receiving calibration data, judging correctness and calculating error compensation. Starting a calibration mode through a calibration command to obtain amplitude and phase sampling values of receiving (1090 MHz) and transmitting (1030 MHz) frequencies of a circular array secondary radar calibration channel, and judging thresholds, the difference of the receiving/transmitting calibration data of each channel and the difference of the receiving/transmitting calibration data of the channels;
as shown in fig. 3, the specific process is as follows:
step 11: receiving amplitude and phase sampling values of a circular array secondary radar calibration channel at receiving frequency (1090 MHz) and transmitting frequency (1030 MHz), and respectively storing a calibration receiving amplitude value, a calibration receiving phase value, a calibration transmitting amplitude value and a calibration transmitting phase according to channels, wherein the phase takes a channel 1 as a reference channel, so that the receiving phase value and the transmitting phase value of the channel 1 are both 0;
step 12: performing first calibration judgment whether the current process is the current process, if the current process is the first calibration, sending a calibration request, returning to the step 11, and if the current process is not the first calibration judgment of the current process, performing the step 13;
step 13: respectively judging data validity of the calibration receiving amplitude value and the calibration transmitting amplitude value, the calibration receiving phase value and the calibration transmitting phase, wherein the amplitude value needs to respectively meet a receiving amplitude value threshold (the sampling threshold value set by the invention is 90) and a transmitting amplitude value threshold (the sampling threshold value set by the invention is 90), the difference between the receiving amplitude value and the transmitting amplitude value of the same channel between two continuous calibrations is within 1dB (namely the sampling value is within 3), the phase value needs to meet the condition that the sampling difference between the receiving phase value and the transmitting phase value of the same channel between two continuous calibrations is within 10 degrees, if the condition is met, the step 2 is carried out, if the condition is not met, the calibration process is restarted once, and whether the calibration process is carried out for 3 times continuously is judged.
Step 2: and (5) performing adaptive adjustment and linear calculation on the standard beam forming coefficient table. When the amplitude and the phase of the external or internal factors deviate again, the recalibration and compensation cannot be performed, at the moment, the standard beam forming coefficient needs to be adjusted in a sectional mode, and the temporary matrix B of NxMxL is obtained by reading, writing and storing through desktop software 1 ；
As shown in fig. 4, the specific process is as follows:
step 21: judging whether the calibration receiving or calibration transmitting amplitude value in the step 1 is normal, if all channels are normal, entering the step 22, if so, checking the corresponding channel fault, ending the process, and not causing errors in general amplitude sampling, so that special attention needs to be paid to phase sampling;
step 22: starting normal air-to-air work, analyzing the condition of each wave position secondary radar response data, if the number of each wave position response data is equivalent to the number of inquiry triggers, entering step 23, and if the number of each wave position response data is far less than the number of inquiry triggers, entering step 26;
step 23: when the air-to-air system works normally, the sum-difference amplitude of each wave position is normal, that is, the sum-amplitude is higher than the difference amplitude, the standard weight coefficient is not required to be modified, and the default standard weight coefficient or the weight coefficient matrix adjusted in step 26 is adopted to enter step 31. In the invention, a single wave bit is transmitted by 14 channels participating in DBF: standard amplitude weighting factors of [13,11,8,7,4,2,0,0,2,4,7,8,11,13], standard phase weighting factors of [180,2,168,310,63,141,180,0,321,243,130,348,182,0], receive: the standard and amplitude weighting coefficients are [13,11,8,7,4,2,0,0,2,4,7,8,11,13], the standard deviation amplitude weighting coefficients are [13,11,8,7,4,2,0,0,2,4,7,8,11,13], the standard and phase weighting coefficients are [180,2,168,310,63,141,180,0,321,243,130,348,182,0], and the standard deviation phase weighting coefficients are [0,182,348,130,243,321,0,0,321,243,130,348,182,0 ]. If the sum amplitude of a certain wave position (beam _ x) is lower than the difference amplitude, entering the next step;
step 24: adjusting the center position (7 th or 8 th) of a certain wave position (beam _ x) to correspond to the channel switchAfter closing (the attenuation value of the standard amplitude weight coefficient is changed from 0dB to 13dB through interface adjustment), a temporary matrix B is obtained through linear circulation calculation 1 Entering step 3, starting normal air-to-air operation again, if the secondary radar response data and the difference amplitude of the wave position (beam _ x) are normal, restoring the attenuation value of the standard amplitude weight coefficient received by the channel corresponding to the central position (7 th or 8 th) of the wave position (beam _ x) from 13dB to 0dB, adjusting the standard sum or difference phase weight coefficient value received by the corresponding position (7 th or 8 th) through an interface, adjusting by taking 90 degrees or 45 degrees as stepping, entering step 27, and if the secondary radar response data and the difference amplitude are still abnormal, entering step 25;
step 25: adjusting the sub-center position (6 th or 9 th) of a certain wave position (beam _ x) to close the corresponding channel (the method is the same as the step 24), and performing linear cycle calculation to obtain a temporary matrix B 1 Entering step 3, starting normal air-to-air work again, if the secondary radar response data and the difference amplitude ratio are normal, restoring the receiving standard amplitude weight coefficient attenuation value of the channel corresponding to the secondary central position (beam _ x) (the 6 th or the 9 th) from 13dB to 0dB, adjusting the receiving standard and/or difference phase weight coefficient value of the corresponding position (the 6 th or the 9 th) through an interface, generally taking 90 degrees as stepping for adjustment, entering step 27, if the secondary radar response data and the difference amplitude ratio are still abnormal, checking the corresponding channel of the beam _ x wave position through an instrument, and ending the process;
step 26: the standard phase weight coefficient value is transmitted at the position (7 th or 8 th) corresponding to the beam _ x wave position through interface adjustment, and a temporary matrix B is obtained by linear circulation calculation after adjustment is generally carried out by taking 90 degrees as steps 1 Entering the step 3 and the step 4, starting normal opposite-to-empty work again, entering the step 23 if the number of beam _ x wave bit response data is equal to the number of inquiry triggers, and if the number of each wave bit response data is still far less than the number of inquiry triggers, checking a corresponding channel of the beam _ x wave bit through an instrument to finish the process;
step 27: linear cyclic calculation is carried out to obtain a temporary matrix B 1 Starting normal air work again, observing whether the inquiry response detection performance of wave level (beam _ x) is normal or not, and if so, enteringAnd step 31, if the wave position is abnormal, checking a corresponding channel of the beam _ x wave position through an instrument, and ending the flow.
Closing the corresponding channel of the central position of the problematic wave position:
in order to eliminate the serious beam deformation caused by the error compensation of the phase of the channel corresponding to the central position of the wave position in question, the beam forming condition when the amplitude of the channel corresponding to the central position is closed is firstly checked, and the channel corresponding to the central position of the 9 wave position is closed (the attenuation value of the standard amplitude weight coefficient is changed from 0dB to 13dB or lower through interface adjustment) to generate a variable standard weight coefficient for regenerating a final weight coefficient matrix D;
for the empty observation target, the quality of the response signal is obviously improved, if the channel phase corresponding to the center position of the 9 wave bits is abnormal, the amplitude weight coefficient value of the channel of the 9 wave bits is recovered to be changed from 13dB to 0dB according to the step 24.
And (3) compensating the phase processing of the channel corresponding to the central position of the problematic wave position:
in order to eliminate the serious beam deformation caused by the error of phase compensation of the channel corresponding to the central position of the wave position with the problem, after the step of closing the channel corresponding to the central position of the wave position with the problem, the receiving standard and/or the difference phase weight coefficient value of the channel are determined to be modified, firstly, the condition of phase inversion is considered, the receiving and phase weight coefficient value of the channel in the 9 wave position is changed from 0 degree to 180 degrees, and the receiving difference phase weight coefficient value is changed from 180 degrees to 0 degrees.
And for the empty observation target, the quality of the response signal is obviously improved, the target detection performance is recovered to be normal, a variable standard weight coefficient is generated according to the state and is used for regenerating a final weight coefficient matrix D, and if the quality is not improved, the adjustment is continuously carried out according to the step of 90 degrees or 45 degrees.
And step 3: and performing circular mapping calculation on the matrix array, normalizing the received receiving/transmitting calibration data of each channel by taking a reference channel as a reference respectively, acquiring an amplitude-phase error table, and obtaining an NxL matrix A, wherein N is 30, and L represents the number 6 of transmitting/receiving weight coefficient tables, namely obtaining a 30 x 6 matrix array, wherein the receiving/transmitting amplitude reference channel needs to select the channel with the minimum value, and the phase selects the first channel as the reference channel. Acquisition webAfter the error table is compared, the temporary matrix B in the step 2 is passed 1 Performing cyclic mapping calculation on a matrix array B corresponding to the standard beam forming coefficient obtained by mapping to form a one-by-one correspondence relationship, and obtaining a three-dimensional matrix array C, namely a matrix with the size of N multiplied by M multiplied by L =30 multiplied by 14 multiplied by 6;
as shown in fig. 5, the specific process is as follows:
step 31: during calibration, a receiving and/difference amplitude error table is calculated, namely the rxMeasure _ Amp [ minute ] with the lowest amplitude in N radio frequency channels is found as a reference amplitude value, and the error of the amplitude of each channel relative to the reference amplitude value is calculated and recorded as rxMeasure _ AmpError [ i ];
step 32: calculating a receiving and/or difference Phase error table, taking rxMeasure _ Phase [0] as a reference Phase value, and calculating the error of each channel Phase relative to the reference Phase value, which is recorded as rxMeasure _ PhaseError [ i ];
step 33: calculating a transmission amplitude phase error table in the same steps 31 and 32, and respectively recording as txMeasure _ AmpError [ i ] and txMeasure _ PhaseError [ i ];
step 34: circularly mapping the variable standard weight coefficients to obtain 6 original temporary matrixes B of 14 multiplied by 1 1 Circularly mapping to N wave bits to obtain a matrix array B of 30 multiplied by 14 multiplied by 6;
step 35: and performing wave position cyclic mapping one-to-one corresponding calculation on a matrix A taking rxMeasure _ Amperror [ i ], rxMeasure _ PhaseError [ i ], txMeasure _ Amperror [ i ] and txMeasure _ PhaseError [ i ] as basic data and a variable standard weight coefficient matrix array B to generate a DBF weight coefficient matrix C.
And 4, step 4: and calculating the orthogonal IQ signal, and waiting for DBF formation to use. The three-dimensional matrix array C (matrix of 30 × 14 × 6) obtained in step 3 may obtain I/Q values, i.e., real part and imaginary part of the output vector signal, respectively, through an orthogonal formula, and may obtain the transmitted and received DBF weighting coefficients, respectively, for DBF formation.
As shown in fig. 6, the specific process is as follows:
step 41: obtaining a linear amplitude according to the amplitude/phase value corresponding to the DBF weight coefficient matrix C;
step 42: and further obtaining the orthogonal I/Q value of the signal for wave control DBF calling.
Claims (7)
1. A circular array secondary radar sectional type weight coefficient processing method based on a matrix array is characterized by comprising the following steps:
s1, receiving calibration data of the circular array secondary radar calibration channel and judging the effectiveness of the calibration data;
s2, forming adaptive and linear cyclic calculation of weight coefficients for the calibration data through iterative verification to obtain a temporary matrix B 1 ；
S3, obtaining the amplitude-phase error table, calculating the circular mapping of the matrix array, and calculating the temporary matrix B 1 Performing one-to-one correspondence on the matrix array B obtained by circular mapping to generate a matrix C;
and S4, performing orthogonal IQ signal calculation on the matrix C to generate a final value weight coefficient matrix D for calling a wave control DBF.
2. The matrix array-based segmented weight coefficient processing method for the circular array secondary radar according to claim 1, wherein the specific steps of S1 are as follows:
s11, receiving amplitude and phase sampling values of the circular array secondary radar calibration channel in a receiving and transmitting deviation table, and respectively storing a calibration receiving amplitude value, a calibration receiving phase value, a calibration transmitting amplitude value and a calibration transmitting phase value according to the channels;
s12, judging whether the current calibration is the first calibration of the process, if so, sending a calibration request, returning to the step S11, and otherwise, entering the step S13;
and S13, respectively carrying out data validity judgment on the calibration receiving amplitude value and the calibration transmitting amplitude value, and the calibration receiving phase value and the calibration transmitting phase value, if the validity conditions are met, entering the step S2, and if the validity conditions are not met, returning to the step S11.
3. The matrix array-based segmented weight coefficient processing method for circular array secondary radar according to claim 2, wherein the validity judgment is as follows: the amplitude value needs to respectively meet the conditions that the receiving amplitude value threshold and the transmitting amplitude value threshold are met, and the sampling difference of the calibration receiving amplitude value and the calibration transmitting amplitude value of the same channel between two times of continuous calibration is within 1 dB; the phase value should satisfy the condition that the sampling difference between the calibration receiving phase value and the calibration transmitting phase value of the same channel between two successive calibrations is within 10 °.
4. The matrix array-based segmented weight coefficient processing method for the circular array secondary radar according to claim 1, wherein the specific steps of S2 are as follows:
s21, judging whether the calibration receiving amplitude value and the calibration transmitting amplitude value are normal, if all the channels are normal, entering a step S22, otherwise, checking the corresponding channel fault, and ending the process;
s22, starting normal opposite-air work, analyzing the response data condition of each wave position secondary radar, if the number of the response data of each wave position is equal to the number of inquiry triggers, entering step S23, and if the number of the response data of each wave position is far less than the number of inquiry triggers, entering step S26;
s23, when the air-to-air work is normal, the sum-difference amplitude of each wave position is normal, namely the sum-amplitude difference amplitude is high, the standard weight coefficient does not need to be modified, the default standard weight coefficient is adopted to carry out linear circulation calculation to obtain a temporary matrix B1, the step S3 is carried out, and if the sum-amplitude difference amplitude of a certain wave position beam _ x is low, the step S24 is carried out;
s24, adjusting the wave position beam _ x center position to close the corresponding channel, and performing linear circulation calculation to obtain a temporary matrix B 1 Entering step S3, starting normal air-to-air operation again, if the secondary radar response data and the difference amplitude of the wave position beam _ x are normal, restoring the channel corresponding to the central position of the wave position beam _ x to be started, adjusting the value of the receiving standard and/or difference phase weight coefficient at the corresponding position through an interface, adjusting by taking 90 degrees or 45 degrees as steps, entering step S27, and if the secondary radar response data and the difference amplitude are still abnormal, entering step S25;
s25, adjusting beam position beam _ x times of the center position corresponding to the channel to be closed, and performing linear circulation calculation to obtain the temporary matrix B again 1 Step S3 is entered, normal air work is started again, if the secondary radar responds to the dataIf the sum-difference amplitude is normal, recovering the beam _ x wave position to open a corresponding channel at the central position, adjusting the receiving standard and/or difference phase weight coefficient value at the corresponding position through an interface, performing adjustment by taking 90 degrees as steps, entering a step S27, if the secondary radar response data and the difference amplitude are still abnormal, checking the corresponding channel at the beam _ x wave position through an instrument, and ending the flow;
s26, transmitting standard phase weight coefficient values through the center position of the interface adjusting wave position beam _ x, adjusting by taking 90 degrees as steps, and then performing linear circulation calculation to obtain a temporary matrix B 1 The step S3 is entered, normal air-to-air work is started again, and if the wave position beam _ x response data number is equal to the inquiry trigger number, the step S23 is returned; if the number of the wave position response data is still far less than the number of the inquiry triggers, the corresponding channel of the wave position beam _ x is checked through the instrument, and the process is ended;
s27, linear circulation calculation is carried out to obtain a temporary matrix B again 1 And starting normal air work again, observing whether the inquiry detection performance of the wave position beam _ x is normal, if so, entering the step S3, and if not, checking a corresponding channel of the wave position beam _ x through an instrument, and ending the process.
5. The matrix array-based segmented weight coefficient processing method for the circular array secondary radar, according to claim 4, wherein the wave position participates in DBF by 14 channels, and when in transmission: standard amplitude weighting factors of [13,11,8,7,4,2,0,0,2,4,7,8,11,13], standard phase weighting factors of [180,2,168,310,63,141,180,0,321,243,130,348,182,0], receive: the weighting coefficients for the standard sum amplitude are [13,11,8,7,4,2,0,0,2,4,7,8,11,13], the weighting coefficients for the standard difference amplitude are [13,11,8,7,4,2,0,0,2,4,7,8,11,13], the weighting coefficients for the standard sum phase are [180,2,168,310,63,141,180,0,321,243,130,348,182,0], and the weighting coefficients for the standard difference phase are [0,182,348,130,243,321,0,0,321,243,130,348,182,0 ].
6. The matrix array-based segmented weight coefficient processing method for the circular array secondary radar according to claim 1, wherein the specific steps of S3 are as follows:
s31, calculating a receiving and/difference amplitude error table, namely finding the rxMeasure _ Amp [ minimum ] with the lowest amplitude in N radio frequency channels as a reference amplitude value, and calculating the error of the amplitude of each channel relative to the reference amplitude value, and recording the error as rxMeasure _ AmpError [ i ];
s32, calculating a receiving and/or difference Phase error table, taking the rxMeasure _ Phase [0] as a reference Phase value, and calculating the error of each channel Phase relative to the reference Phase value, which is recorded as rxMeasure _ PhaseError [ i ];
s33, calculating a transmission and/or difference amplitude error table and a transmission and/or difference Phase error table, taking the txMeasure _ Amp [0] and the txMeasure _ Phase [0] as reference amplitude values and Phase values, and calculating the errors of the amplitudes of all channels relative to the reference amplitude values and the phases relative to the reference Phase values, which are marked as txMeasure _ Amperer [ i ] and txMeasure _ PhaseError [ i ];
s34, circularly mapping the variable standard weight coefficients to obtain 6 original temporary matrixes B of 14 multiplied by 1 1 Circularly mapping to N wave bits to obtain a matrix array B of 30 multiplied by 14 multiplied by 6;
s35, carrying out wave position cyclic mapping one-to-one corresponding calculation on a matrix A and a variable scale matrix array B by taking the rxMeasure _ Amperer [ i ], the rxMeasure _ PhaseError [ i ], the rxMeasure _ Amperer [ i ] and the txMeasure _ PhaseError [ i ] as basic data, and generating a DBF weight coefficient matrix C.
7. The matrix array-based segmented weight coefficient processing method for the circular array secondary radar according to claim 1, wherein the specific steps of S4 are as follows:
s41, obtaining a linear amplitude according to the amplitude/phase value corresponding to the DBF weight coefficient matrix C;
and S42, obtaining the orthogonal I/Q value of the signal through linear amplitude calculation, and generating a final value weight coefficient matrix D for calling a wave control DBF.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210953623.5A CN115032600B (en) | 2022-08-10 | 2022-08-10 | Circular array secondary radar sectional weight coefficient processing method based on matrix array |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210953623.5A CN115032600B (en) | 2022-08-10 | 2022-08-10 | Circular array secondary radar sectional weight coefficient processing method based on matrix array |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115032600A true CN115032600A (en) | 2022-09-09 |
| CN115032600B CN115032600B (en) | 2022-11-08 |
Family
ID=83130124
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210953623.5A Active CN115032600B (en) | 2022-08-10 | 2022-08-10 | Circular array secondary radar sectional weight coefficient processing method based on matrix array |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115032600B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115586501A (en) * | 2022-11-25 | 2023-01-10 | 四川九洲电器集团有限责任公司 | FPGA-based multichannel baseband data amplitude-phase compensation implementation method |
Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080042896A1 (en) * | 2006-08-21 | 2008-02-21 | Yair Alon | Using radar targets of opportunity to build a monopulse calibration table |
| CN101251597A (en) * | 2008-04-08 | 2008-08-27 | 西安电子科技大学 | Method for self-correction of array error of multi-input multi-output radar system |
| CN102313881A (en) * | 2011-08-01 | 2012-01-11 | 四川九洲电器集团有限责任公司 | Simple realization method for static calibration on radio frequency channel of secondary radar of broadband amplitude monopulse system |
| CN103383450A (en) * | 2013-06-25 | 2013-11-06 | 西安电子科技大学 | Conformal array radar amplitude-phase error correction fast achieving method |
| CN103399305A (en) * | 2013-06-28 | 2013-11-20 | 四川九洲空管科技有限责任公司 | Amplitude and phase calibration method for radio frequency channels of digital array secondary radar and antenna arrays |
| EP2667217A1 (en) * | 2012-05-22 | 2013-11-27 | Rohde & Schwarz GmbH & Co. KG | A method and a system with an automatic calibration of a multistatic antenna array |
| CN103605116A (en) * | 2013-12-04 | 2014-02-26 | 西安电子科技大学 | Online imaging radar channel parameter compensation method based on sparse analysis |
| CN103941242A (en) * | 2014-05-21 | 2014-07-23 | 中国科学院电子学研究所 | Frequency modulated continuous wave linear array amplitude-phase error correction method based on single dominant scattering center |
| CN104020455A (en) * | 2014-05-30 | 2014-09-03 | 中国电子科技集团公司第二十二研究所 | Ground and sky wave array calibration method based on direct waves |
| CN104597433A (en) * | 2015-01-30 | 2015-05-06 | 中国电子科技集团公司第三十八研究所 | Phased-array antenna multi-beam automatic calibration device and method |
| CN105044693A (en) * | 2015-06-25 | 2015-11-11 | 西安电子科技大学 | Microwave associated imaging radar amplitude-phase error correction method based on auxiliary array elements |
| CN105425221A (en) * | 2015-11-09 | 2016-03-23 | 四川九洲电器集团有限责任公司 | Radio frequency channel real-time calibration method and secondary radar |
| CN108254725A (en) * | 2017-12-20 | 2018-07-06 | 雷象科技（北京）有限公司 | Phased-array radar transmission channel recycles calibration method |
| CN108549058A (en) * | 2018-02-28 | 2018-09-18 | 四川九洲电器集团有限责任公司 | Secondary radar broadband active phased array system and dynamic calibration method thereof |
| CN108614247A (en) * | 2018-04-03 | 2018-10-02 | 四川九洲电器集团有限责任公司 | A kind of secondary radar channel calibration method |
| CN108828595A (en) * | 2018-05-28 | 2018-11-16 | 成都雷通科技有限公司 | A kind of control method of S-band phased-array radar |
| CN109901150A (en) * | 2019-03-04 | 2019-06-18 | 四川九洲空管科技有限责任公司 | A kind of multifunction array radar device and its detection method |
| CN110456316A (en) * | 2019-07-05 | 2019-11-15 | 四川九州空管科技有限责任公司 | A kind of multifunction array radar multichannel transmitting-receiving correction system and method |
| CN112083389A (en) * | 2020-06-09 | 2020-12-15 | 成都微联云智科技有限责任公司 | Real-time calibration method for circular array antenna |
| CN112083382A (en) * | 2020-09-15 | 2020-12-15 | 四川九洲空管科技有限责任公司 | High-precision small phased array secondary radar azimuth compensation calculation method |
| CN112269167A (en) * | 2020-09-29 | 2021-01-26 | 四川九洲空管科技有限责任公司 | Intelligent self-adaptive calibration method for secondary radar phased array antenna |
-
2022
- 2022-08-10 CN CN202210953623.5A patent/CN115032600B/en active Active
Patent Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080042896A1 (en) * | 2006-08-21 | 2008-02-21 | Yair Alon | Using radar targets of opportunity to build a monopulse calibration table |
| CN101251597A (en) * | 2008-04-08 | 2008-08-27 | 西安电子科技大学 | Method for self-correction of array error of multi-input multi-output radar system |
| CN102313881A (en) * | 2011-08-01 | 2012-01-11 | 四川九洲电器集团有限责任公司 | Simple realization method for static calibration on radio frequency channel of secondary radar of broadband amplitude monopulse system |
| EP2667217A1 (en) * | 2012-05-22 | 2013-11-27 | Rohde & Schwarz GmbH & Co. KG | A method and a system with an automatic calibration of a multistatic antenna array |
| CN103383450A (en) * | 2013-06-25 | 2013-11-06 | 西安电子科技大学 | Conformal array radar amplitude-phase error correction fast achieving method |
| CN103399305A (en) * | 2013-06-28 | 2013-11-20 | 四川九洲空管科技有限责任公司 | Amplitude and phase calibration method for radio frequency channels of digital array secondary radar and antenna arrays |
| CN103605116A (en) * | 2013-12-04 | 2014-02-26 | 西安电子科技大学 | Online imaging radar channel parameter compensation method based on sparse analysis |
| CN103941242A (en) * | 2014-05-21 | 2014-07-23 | 中国科学院电子学研究所 | Frequency modulated continuous wave linear array amplitude-phase error correction method based on single dominant scattering center |
| CN104020455A (en) * | 2014-05-30 | 2014-09-03 | 中国电子科技集团公司第二十二研究所 | Ground and sky wave array calibration method based on direct waves |
| CN104597433A (en) * | 2015-01-30 | 2015-05-06 | 中国电子科技集团公司第三十八研究所 | Phased-array antenna multi-beam automatic calibration device and method |
| CN105044693A (en) * | 2015-06-25 | 2015-11-11 | 西安电子科技大学 | Microwave associated imaging radar amplitude-phase error correction method based on auxiliary array elements |
| CN105425221A (en) * | 2015-11-09 | 2016-03-23 | 四川九洲电器集团有限责任公司 | Radio frequency channel real-time calibration method and secondary radar |
| CN108254725A (en) * | 2017-12-20 | 2018-07-06 | 雷象科技（北京）有限公司 | Phased-array radar transmission channel recycles calibration method |
| CN108549058A (en) * | 2018-02-28 | 2018-09-18 | 四川九洲电器集团有限责任公司 | Secondary radar broadband active phased array system and dynamic calibration method thereof |
| CN108614247A (en) * | 2018-04-03 | 2018-10-02 | 四川九洲电器集团有限责任公司 | A kind of secondary radar channel calibration method |
| CN108828595A (en) * | 2018-05-28 | 2018-11-16 | 成都雷通科技有限公司 | A kind of control method of S-band phased-array radar |
| CN109901150A (en) * | 2019-03-04 | 2019-06-18 | 四川九洲空管科技有限责任公司 | A kind of multifunction array radar device and its detection method |
| CN110456316A (en) * | 2019-07-05 | 2019-11-15 | 四川九州空管科技有限责任公司 | A kind of multifunction array radar multichannel transmitting-receiving correction system and method |
| CN112083389A (en) * | 2020-06-09 | 2020-12-15 | 成都微联云智科技有限责任公司 | Real-time calibration method for circular array antenna |
| CN112083382A (en) * | 2020-09-15 | 2020-12-15 | 四川九洲空管科技有限责任公司 | High-precision small phased array secondary radar azimuth compensation calculation method |
| CN112269167A (en) * | 2020-09-29 | 2021-01-26 | 四川九洲空管科技有限责任公司 | Intelligent self-adaptive calibration method for secondary radar phased array antenna |
Non-Patent Citations (3)
| Title |
|---|
| 吕江涛等: "航管二次雷达测量精度分析及改善方法", 《现代雷达》 * |
| 张彭: "雷神二次航管雷达的方位校准程序的几点经验", 《空中交通管理》 * |
| 王世民等: "一种二次雷达相控阵系统的自动校准方法", 《国外电子测量技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115586501A (en) * | 2022-11-25 | 2023-01-10 | 四川九洲电器集团有限责任公司 | FPGA-based multichannel baseband data amplitude-phase compensation implementation method |
| CN115586501B (en) * | 2022-11-25 | 2023-03-10 | 四川九洲电器集团有限责任公司 | FPGA-based multichannel baseband data amplitude-phase compensation implementation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115032600B (en) | 2022-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2476163B1 (en) | Antenna failure compensation | |
| WO2022095510A1 (en) | Calibration method and device for array antenna, and storage medium | |
| CN111638494A (en) | Multi-channel amplitude and phase correction method for digital beam synthesis system | |
| CN105137409B (en) | The sane space-time adaptive processing method of echo signal mutually constrained based on width | |
| CN110361705B (en) | Phased array antenna near field iterative calibration method | |
| CN115032600B (en) | Circular array secondary radar sectional weight coefficient processing method based on matrix array | |
| CN110824414A (en) | Device and method for estimating angle of arrival | |
| CN110208757B (en) | Steady self-adaptive beam forming method and device for inhibiting main lobe interference | |
| CN110045338A (en) | Amplitude phase error estimation and bearing calibration based on Hilbert transform | |
| CN113835063B (en) | Unmanned aerial vehicle array amplitude and phase error and signal DOA joint estimation method | |
| CN117310707B (en) | On-orbit azimuth direction diagram extraction method for azimuth multichannel satellite-borne SAR | |
| CN113156380B (en) | Main lobe noise suppression interference self-adaptive suppression method based on transceiving mismatch processing | |
| CN115494462A (en) | Phased array radar online automatic amplitude and phase calibration method and storage medium | |
| CN111817765B (en) | Generalized sidelobe cancellation broadband beam forming method based on frequency constraint | |
| CN109407047B (en) | Amplitude-phase error calibration and direction-of-arrival estimation method based on rank loss root finding | |
| CN112858994B (en) | Amplitude comparison direction finding method based on uniform circular array | |
| CN116679306A (en) | Extended aperture sonar imaging method based on coefficient correction | |
| CN1577967A (en) | Method for reducing the radiation load by a mobile radio terminal and a mobile radio terminal | |
| CN115372925A (en) | Array robust adaptive beam forming method based on deep learning | |
| CN112147590B (en) | Channel equalization method based on response estimation frequency domain fitting | |
| CN114563764A (en) | Self-adaptive beam generation method and system | |
| KR100762218B1 (en) | Apparatus for calibrating transmitters and receivers in array antenna system | |
| CN111431575A (en) | Incoming wave direction sparse reconstruction method based on conventional beam forming | |
| CN115412414A (en) | Digital predistortion feedback system | |
| CN114609579B (en) | Defocus direction finding error correction method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |