CN111965602B - Phased array radar amplitude-phase consistency detection method and system - Google Patents
Phased array radar amplitude-phase consistency detection method and system Download PDFInfo
- Publication number
- CN111965602B CN111965602B CN201911128100.1A CN201911128100A CN111965602B CN 111965602 B CN111965602 B CN 111965602B CN 201911128100 A CN201911128100 A CN 201911128100A CN 111965602 B CN111965602 B CN 111965602B
- Authority
- CN
- China
- Prior art keywords
- phase
- probe
- probes
- vibrators
- amplitude
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 23
- 239000000523 sample Substances 0.000 claims abstract description 213
- 238000000034 method Methods 0.000 claims abstract description 72
- 230000005855 radiation Effects 0.000 claims abstract description 44
- 238000010586 diagram Methods 0.000 claims description 21
- 238000012360 testing method Methods 0.000 claims description 21
- 238000005259 measurement Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 description 15
- 230000008878 coupling Effects 0.000 description 6
- 238000010168 coupling process Methods 0.000 description 6
- 238000005859 coupling reaction Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000012423 maintenance Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000011179 visual inspection Methods 0.000 description 1
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/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
- G01S7/022—Road traffic radar detectors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S2013/0236—Special technical features
- G01S2013/0245—Radar with phased array antenna
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 embodiment of the application discloses a phased array radar amplitude phase consistency detection method and a phased array radar amplitude phase consistency detection system, wherein the method comprises the steps of setting N probes in a preset range, wherein the preset range is determined according to the condition that M vibrators and N probes are mutually positioned in the radiation range of each other; acquiring initial data of M vibrators and N probes, wherein the initial data comprise coordinate positions and pointing angles of the M vibrators, radiation patterns of the M vibrators and the N probes and coordinate positions of a phase center; further, N probes are utilized to respectively measure M vibrators one by one to obtain an amplitude set A i,j And phase set P i,j The method comprises the steps of carrying out a first treatment on the surface of the Further, according to the amplitude set A i,j The phase set P i,j And determining the amplitude and the phase of the M vibrators, the coordinate positions and the pointing angles of the N probes by utilizing the initial data and utilizing a least square method. Therefore, the phase uncertainty caused by the cable motion can be avoided, and the purpose of assisting the phased array radar to maintain the amplitude phase consistency is directly achieved.
Description
Technical Field
The embodiment of the application relates to the technical field of intelligent transportation, in particular to a phased array radar amplitude-phase consistency detection method and system.
Background
All radars need to ensure a certain amplitude consistency in order for the receiver to operate at an optimal state. The large-scale phased array radar has huge vibrator quantity, and problems of individual vibrators, such as vibrator failure or overlarge amplitude phase deviation, are unavoidable in the use process. And therefore require periodic maintenance checks. However, the maintenance mode in the prior art has the problems of serious time and labor consumption and huge cost.
Therefore, when maintaining the amplitude-phase consistency of the phased array radar, a more convenient, time-saving and labor-saving manner is needed.
Disclosure of Invention
Therefore, the embodiment of the application provides a phased array radar amplitude consistency detection method and a phased array radar amplitude consistency detection system, which can assist in maintaining phased array radar amplitude consistency.
In order to achieve the above object, the embodiment of the present application provides the following technical solutions:
according to a first aspect of an embodiment of the present application, there is provided a phased array radar amplitude phase consistency detection method, the method including:
setting N probes in a preset range, wherein the preset range is determined according to the condition that M vibrators and N probes are mutually positioned in the radiation range of each other, and M and N are integers larger than 1;
acquiring initial data of M vibrators and N probes, wherein the initial data comprise coordinate positions and pointing angles of the M vibrators, radiation patterns of the M vibrators and the N probes and coordinate positions of a phase center;
n probes are used for measuring M vibrators one by one to obtain an amplitude set A i,j And phase set P i,j The method comprises the steps of carrying out a first treatment on the surface of the Wherein i is an integer from 1 to N, j is an integer from 1 to M;
according to the amplitude set A i,j The phase set P i,j And determining the amplitude and the phase of the M vibrators, the coordinate positions and the pointing angles of the N probes by utilizing the initial data and utilizing a least square method.
Optionally, the determining the amplitude and the phase of the M vibrators and the coordinate positions and the pointing angles of the N probes by using the least square method is calculated according to the following formula:
A i,j ＝Ar j Ap i abs[Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j )]/|r i,j |
P i,j ＝mod[Pr j +Pp i +angle(Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j ))-k|r i,j |,2π]
wherein Ar is j Representing the amplitude of vibrator j, ap i Representing the amplitude, θ, of the probe i j Represents the pointing angle, eta of vibrator j i Representing the pointing angle of the probe i; zeta type toy i,j Representing the corresponding angle coordinates of the connecting line between the two; dr (θ) j ,ξ i,j ) Representing the value of the radiation pattern of vibrator j at probe i, dp (eta) i ,ξ i,j ) Representative probei, taking a value of the directional diagram at a vibrator j, wherein the product of the directional diagram and the vibrator represents a three-dimensional vector inner product under the condition of an electromagnetic phased array and represents a common product under the condition of an acoustic phased array; pr (Pr) j Representing the phase, pp, of vibrator j i Representing the phase of probe i; r is |r i,j The i represents the distance between the phase center of probe i and the phase center of vibrator j; k represents wave number; abs () represents a modulo function, angle () represents a phase function, and mod () represents a remainder function.
Optionally, the measuring M vibrators with N probes one by one includes:
respectively measuring M vibrators by using N probes at different positions one by one; or alternatively
A single probe is used to move to different positions to measure M vibrators separately.
Optionally, the M transducers and the N probes are located within the radiation range of each other, provided that the position of each probe is within the set beam width of all or part of the transducers, and the probe pointing direction satisfies that all or part of the transducers are within the set beam width of the probe.
Optionally, the position of each probe further includes: the shortest distance from each probe to the vibrator array is more than 2d 2 λ, where d represents the vibrator radiation aperture and λ represents the test wavelength.
According to a second aspect of an embodiment of the present application, there is provided a phased array radar amplitude phase consistency detection system, the system comprising:
the probe setting module is used for setting N probes in a preset range, wherein the preset range is determined according to the condition that M vibrators and N probes are mutually positioned in the radiation range of each other, and M and N are integers larger than 1;
the initial data acquisition module is used for acquiring initial data of the M vibrators and the N probes, wherein the initial data comprise coordinate positions and pointing angles of the M vibrators, radiation patterns of the M vibrators and the N probes and phase center coordinate positions;
the measuring module is used for measuring M vibrators by utilizing N probes one by one to obtain a webDegree set A i,j And phase set P i,j The method comprises the steps of carrying out a first treatment on the surface of the Wherein i is an integer from 1 to N, j is an integer from 1 to M;
an amplitude-phase data calculation module for calculating the amplitude set A i,j The phase set P i,j And determining the amplitude and the phase of the M vibrators, the coordinate positions and the pointing angles of the N probes by utilizing the initial data and utilizing a least square method.
Optionally, the amplitude and phase data calculation module specifically calculates according to the following formula:
A i,j ＝Ar j Ap i abs[Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j )]/|r i,j |
P i,j ＝mod[Pr j +Pp i +angle(Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j ))-k|r i,j |,2π]
wherein Ar is j Representing the amplitude of vibrator j, ap i Representing the amplitude, θ, of the probe i j Represents the pointing angle, eta of vibrator j i Representing the pointing angle of the probe i; zeta type toy i,j Representing the corresponding angle coordinates of the connecting line between the two; dr (θ) j ,ξ i,j ) Representing the value of the radiation pattern of vibrator j at probe i, dp (eta) i ,ξ i,j ) Representing the value of the directional diagram of the probe i at the vibrator j, wherein the product of the directional diagram and the vibrator represents a three-dimensional vector inner product under the condition of an electromagnetic phased array and represents a common product under the condition of the acoustic phased array; pr (Pr) j Representing the phase, pp, of vibrator j i Representing the phase of probe i; r is |r i,j The i represents the distance between the phase center of probe i and the phase center of vibrator j; k represents wave number; abs () represents a modulo function, angle () represents a phase function, and mod () represents a remainder function.
Optionally, the measurement module is specifically configured to: respectively measuring M vibrators by using N probes at different positions one by one; or moved to different positions with a single probe to measure M vibrators separately.
Optionally, the M transducers and the N probes are located within the radiation range of each other, provided that the position of each probe is within the set beam width of all or part of the transducers, and the probe pointing direction satisfies that all or part of the transducers are within the set beam width of the probe.
Optionally, the position of each probe further includes: the shortest distance from each probe to the vibrator array is more than 2d 2 λ, where d represents the vibrator radiation aperture and λ represents the test wavelength.
In summary, according to the phased array radar amplitude phase consistency detection method provided by the embodiment of the application, the N probes are arranged in the preset range, and the preset range is determined according to the condition that the M vibrators and the N probes are mutually in the radiation range of each other; acquiring initial data of M vibrators and N probes, wherein the initial data comprise coordinate positions and pointing angles of the M vibrators, radiation patterns of the M vibrators and the N probes and coordinate positions of a phase center; further, N probes are utilized to respectively measure M vibrators one by one to obtain an amplitude set A i,j And phase set P i,j The method comprises the steps of carrying out a first treatment on the surface of the According to the amplitude set A i,j The phase set P i,j And determining the amplitude and the phase of the M vibrators, the coordinate positions and the pointing angles of the N probes by utilizing the initial data and utilizing a least square method. Therefore, the phase uncertainty caused by the cable motion can be avoided, and the purpose of assisting the phased array radar to maintain the amplitude phase consistency is directly achieved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below. It will be apparent to those of ordinary skill in the art that the drawings in the following description are exemplary only and that other implementations can be obtained from the extensions of the drawings provided without inventive effort.
The structures, proportions, sizes, etc. shown in the present specification are shown only for the purposes of illustration and description, and are not intended to limit the scope of the application, which is defined by the claims, so that any structural modifications, changes in proportions, or adjustments of sizes, which do not affect the efficacy or the achievement of the present application, should fall within the scope of the application.
Fig. 1 is a schematic flow chart of a phased array radar amplitude phase consistency detection method provided by an embodiment of the application;
fig. 2 is a schematic diagram of an embodiment of a phased array radar amplitude phase consistency detection method according to the present application;
FIG. 3 is a schematic diagram of a second embodiment of a phased array radar amplitude phase consistency detection method according to the present application;
fig. 4 is a block diagram of a phased array radar amplitude phase consistency detection system according to an embodiment of the present application.
Detailed Description
Other advantages and advantages of the present application will become apparent to those skilled in the art from the following detailed description, which, by way of illustration, is to be read in connection with certain specific embodiments, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
The antenna array plane of a phased array radar is also composed of a number of radiating elements and receiving elements (called array elements), the number of elements being dependent on the functioning of the radar and can be from hundreds to tens of thousands. The units are regularly arranged on a plane to form an array antenna. The electromagnetic wave coherence principle is utilized, and the phase of the current fed to each radiation unit is controlled by a computer, so that the direction of the beam can be changed to scan, and the scanning is called electric scanning. The radiation unit sends the received echo signals to the host computer to complete the searching, tracking and measuring of the radar on the target.
Each antenna element has a device such as a phase shifter in addition to an antenna element. Different oscillators can be fed with currents of different phases through phase shifters, so that beams of different directivities are radiated in space. The greater the number of elements of the antenna, the more orientations the beam is in space possible. The array control radar equipment is complex and expensive, the beam scanning range is limited, and the maximum scanning angle is 90-120 degrees. When the omnibearing monitoring is needed, 3-4 antenna array planes are needed to be configured. Compared with a mechanical scanning radar, the phased array radar has the advantages of more flexible scanning, more reliable performance and stronger anti-interference capability, and can be rapidly adapted to the change of battlefield conditions.
All radars need to ensure a certain amplitude consistency in order for the receiver to operate at an optimal state. But the amplitude-phase consistency has the greatest influence on the monopulse radar, and the monopulse technology is mainly used in the precise tracking measurement radar. The monopulse radar adopts amplitude or phase sum difference type, the receiver is provided with a plurality of difference channels and auxiliary channels besides sum channels, and in order to improve measurement accuracy and realize automatic tracking of a servo system to a target, the signal amplitude and phase of the difference channels are required to be consistent with those of the sum channels.
In general, a test or pilot signal can be respectively input to the sum and difference channels in a period before each target processing period of the radar so as to measure the errors of the amplitudes and phases of the sum and difference channels, an automatic control system is used for forming control signals of the amplitudes and phases, and the amplitudes and phases of the signals of the sum and difference channels are controlled and adjusted so that the errors of the amplitudes and phases of the signals of the sum and difference channels are reduced, thereby realizing amplitude consistency.
In the prior art, aiming at amplitude-phase consistency calibration, the existing schemes comprise a planar near-field scanning method, an inverse matrix method, a middle-field method relying on accurate positioning, an external-field spherical multi-probe scanning method and a mutual coupling self-calibration method. The planar near field scanning method needs to mechanically scan each vibrator point by point, so that time consumption and cost are high; the inverse matrix method needs to strictly and accurately position the vibrator, and the application of the external field is inconvenient; the middle field method relying on accurate positioning also faces the difficulties of inconvenient application of the external field and low precision; the outfield spherical surface multi-probe scanning method does not need accurate positioning, and the scanning time is acceptable, but the number of probes is too large, the equipment is complex and expensive, and even the import is relied on. The self-calibration method based on mutual coupling of vibrators does not need extra equipment, but the applicable array form is limited, and the precision is not high due to the fact that self coupling between vibrators is small.
The embodiment of the application provides a phased array radar amplitude-phase consistency detection method, which is a probe number controllable calibration system without a probe accurate position, can be used for large-scale phased array radar outfield maintenance and repair, improves detection efficiency and reduces detection cost.
Fig. 1 shows a schematic diagram of a phased array radar amplitude phase consistency detection method provided by an embodiment of the present application, where the method includes the following steps:
step 101: setting N probes in a preset range, wherein the preset range is determined according to the condition that M vibrators and N probes are mutually positioned in the radiation range of each other, and M and N are integers larger than one.
Step 102: initial data of M vibrators and N probes are obtained, wherein the initial data comprise coordinate positions and pointing angles of the M vibrators, radiation patterns of the M vibrators and the N probes and coordinate positions of a phase center.
Step 103: n probes are used for measuring M vibrators one by one to obtain an amplitude set A i,j And phase set P i,j The method comprises the steps of carrying out a first treatment on the surface of the Where i is an integer from 1 to N and j is an integer from 1 to M.
Step 104: according to the amplitude set A i,j The phase set P i,j And determining the amplitude and the phase of the M vibrators, the coordinate positions and the pointing angles of the N probes by utilizing the initial data and utilizing a least square method.
In step 101, the condition that the M transducers and the N probes are in the radiation range of each other is that the position of each probe is within the set beam width of all or part of the transducers, and the probe pointing direction satisfies that all or part of the transducers are within the set beam width of the probe.
In one possible embodiment, the distribution of the probes should be such that all the main radiation lobes of the radar elements are covered by at least one main lobe of the probe and that the probe and the individual elements are in the far-field radiation zone of each other. Arrangement of each probeThe location conditions further include: the shortest distance from each probe to the vibrator array is more than 2d 2 λ, where d represents the vibrator radiation aperture and λ represents the test wavelength.
For example, the total number of the probe channels is at least 2, and specifically, the probe channels can be 1 probe in two modes, or two independent probes, or a plurality of probes. The position of each probe is within the main beamwidth of all or part of the transducers (e.g. 3 db) and the probe pointing should be such that all or part of the transducers are within the main beamwidth of the probe (e.g. 3 db). A pair of a probe and a transducer, which are in the main radiation range of each other, are called the probe effectively covering the transducer. The number of probes is selected according to the principle that at least each vibrator in the array can be effectively covered by at least one probe channel; likewise, a large number of jointly covered vibrators improves accuracy. Each probe channel can cover at least 6 vibrators together with at least one other probe channel.
In step 102, the precise position coordinates and pointing angle (expressed in euler angles) of each transducer in the array are obtained, and the pattern of the transducer in the main radiation direction and the position of the transducer phase center are obtained. The above information can be obtained from the phased array producer and the producer of the probe. In addition, the patterns of the vibrator and the probe can also be obtained by an autonomous measurement mode.
In step 103, the measurement methods include, but are not limited to, the following two methods: respectively measuring M vibrators by using N probes at different positions one by one; or the single probe is used for moving to different positions, and M vibrators are respectively measured. I.e. the probe can be a plurality of independent probes or one probe can mechanically scan and move a plurality of positions. That is, the measurement mode provided by the embodiment of the application may or may not involve mechanical scanning.
The radar array to be tested can be various in form, can be a plane array or a ring array of a multi-probe spherical near field. And measuring the jth vibrator by using the ith probe, wherein all other probes and vibrators are in a closed state during the test until all probes and vibrators are traversed. Obtaining the content amplitudeThe signals of information and phase information, i.e. amplitude set A i,j The phase set P i,j . Assuming that there are N probes and M vibrators in total, at most m×n sets of test data can be obtained. However, in practice, the effective test data volume is less than mxn because not every set of probe and transducer combinations satisfies the principle of effective coverage of each other described above.
The pattern of the radiation electromagnetic field of an antenna distributed over a fixed distance with angular position is called a pattern. The phase is referred to as a phase pattern. And respectively moving the reference origins of the radar vibrator and the probe to respective phase center points to obtain a zero-phase directional diagram, wherein the zero-phase directional diagram is amplitude with no phase value of the directional diagram in each direction.
In step 104, the amplitudes and phases of the M transducers, and the coordinate positions and pointing angles of the N probes are determined by using the least square method, and are calculated according to the following formulas (1) and (2):
A i,j ＝Ar j Ap i abs[Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j )]/|r i,j |………………(1)
P i,j ＝mod[Pr j +Pp i +angle(Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j ))-k|r i,j |,2π]……(2)
wherein Ar is j Representing the amplitude of vibrator j, ap i Representing the amplitude, θ, of the probe i j Represents the pointing angle, eta of vibrator j i Representing the pointing angle of the probe i; zeta type toy i,j Representing the corresponding angle coordinates of the connecting line between the two; dr (θ) j ,ξ i,j ) Representing the value of the radiation pattern of vibrator j at probe i, dp (eta) i ,ξ i,j ) Representing the value of the directional diagram of the probe i at the vibrator j, wherein the product of the directional diagram and the vibrator represents a three-dimensional vector inner product under the condition of an electromagnetic phased array and represents a common product under the condition of the acoustic phased array; pr (Pr) j Representing the phase, pp, of vibrator j i Representing the phase of probe i; r is |r i,j The i represents the distance between the phase center of probe i and the phase center of vibrator j; k represents wave number; abs ()Representing a modulo function, angle () representing a phase function, mod () representing a remainder function.
On the premise that the probe and the oscillator directional diagrams as well as the position coordinates and the pointing angles of the oscillators are known, the amplitude, the phase and the position and the pointing angles of the oscillators can be solved as unknowns through a least square method from the formula (1) and the formula (2). The solution flow using the least square method has various options, and the embodiment of the application exemplarily provides a solution flow, which includes the following steps:
step 1, selecting two probe channels i with the maximum number of commonly covered vibrators 1 And i 2 The average value of the included angles between the two probes and the common covered vibrator connecting line is larger than 20 degrees.
And 2, bringing the test data of the two probes and all the commonly covered vibrators into a formula (1), and completely eliminating the unknown quantity of the amplitudes of all the vibrators in the formula (1) in a mode that the two probes detect the signals of the same vibrator, so as to finally obtain a series of equations (3) about the amplitude ratio, the positions and the Euler angles of the two probes.
The amplitude ratio of the two probes and the position coordinates and Euler angles of the two probes are determined by a least square method. The least square method generally starts from an initial value of an unknown amount and gradually iteratively converges to an optimal solution. In this step, the initial value of the least squares method may be set as follows: the amplitude ratio of the two probes is set to 1, and initial values of the positions and angles of the probes can be obtained by a low-precision, low-cost measuring tool or visual inspection.
Step 3: and (3) bringing the test data of the two probes and all the commonly covered vibrators into a formula (2), completely eliminating the unknown quantity of the amplitudes of all the vibrators in the formula (2) in a mode of subtracting the signals of the same vibrator detected by the two probes, and finally obtaining a series of equations (4) about the phase difference and the position of the two probes.
The term related to the phase of the directional diagram in the formula (4) can be directly calculated by using the data of the probe position and the euler angle in the step 1. And (3) solving the formula (4) by using a least square method, and obtaining an accurate solution of the phase difference and the position of the two probes. The phase difference of the initial values may be set to 0, and the position and euler angle may be set to the result in step 1.
Step 4: using the determined precise positions and euler angles of the two probes, the amplitudes and phases of all the transducers covered by at least one of the two probes are determined by equations (1) and (2). Probe i in solving process 1 The amplitude and phase of (a) can be arbitrarily specified as 1 and 0.
Step 5: and (3) determining the probe with the largest number of vibrators which can cover all vibrators in the steps 1-4 (the union of vibrators covered by two probes) in the rest probes, and adding the probe into a solving range.
Step 6: and (3) solving the position and Euler angle of the probe and the amplitude and phase of the probe by using the oscillator phase and amplitude which are determined in the steps 1-4 and the formulas (1) and (2) through a least square method. The amplitude phase of the initial values was set to 1 and 0, respectively, and the position and the euler angle were measured with low accuracy (plus or minus 50%) or visually.
And 7, according to the position and Euler angle of the probe in the step 6, the amplitude and the phase of all vibrators covered by the probe are obtained.
Step 8: and taking the union of all solved probes and the vibrator covered by the probes as to-be-quantified, and solving equations (1) and (2) by using a least square method to obtain a more accurate solution. The initial value of the least square method can be the result of the probe and the vibrator determined in all steps before the step 7.
Step 9: if the solving time is required, the step 9 can be omitted; or after adding a plurality of probes to the solving range, executing the step 9 once.
And step 10, jumping back to the step 5, adding one probe at a time as a solving object, and continuously solving the rest probes and the rest vibrators until all the probes and vibrators are solved.
Step 11: for the vibrators at the edge of the array, the radiation pattern may deviate from the pattern at the center of the array, so that the actual measurement and correction can be performed on the pattern. Then substituting into the formulas (1) and (2), and solving by using a least square method.
In summary, compared with the near field scanning method, the technical scheme provided by the embodiment of the application omits a near field scanning system with huge cost, and reduces the test time. Compared with an external field method and an inverse matrix method, the accurate positioning step of the probe is omitted, and the testing flow is shortened. The scheme provided by the embodiment of the application can avoid phase uncertainty caused by cable motion, has low requirements on test conditions, and is particularly suitable for testing large-scale outfield arrays.
It should be noted that the method provided by the embodiment of the application is also suitable for the rapid calibration of the amplitude phase consistency of the large-scale antenna array and the sonar array of the communication base station.
In order to further verify and explain the effect of the phased array radar amplitude phase consistency detection method provided by the embodiment of the application, the following is exemplified:
a 20 x 20 acoustic array is constructed, each transducer being an isotropic radiation source and the distance between adjacent transducers being half a wavelength. The oscillator amplitude and phase are determined by random numbers. The linear amplitude is uniformly distributed between 0.5 and 2, and the phase is uniformly distributed between 0 and 360 degrees. Assuming that both probes are isotropic (no euler angle is required as a parameter), the linear amplitude and phase are 1 and 0 degrees respectively. Two probes were placed in planes 5 wavelengths from the array, with internal coordinates (1, 5) and (9, 5), respectively, in wavelength.
Based on the above setting conditions, a test data set A is generated ij And P ij . In order to simulate the actual situation, the test error of plus or minus 0.1db is added on the amplitude, and the test error of plus or minus 1 degree is added on the phase.
Further, according to test data A ij And P ij Andand solving the amplitude and the phase of all the vibrators by a least square method according to the position of each vibrator in the acoustic array.
Further, the error between the result obtained by the comparison calculation and the real value data of the vibrator is small.
Finally, as shown in fig. 3 and 4, the probe positions (1.0015,8.9989,5.0009) and (5.0010,4.9976,4.9972) that are reversely deduced from the noisy simulation test data are coincident with the true values of the two probes. The phase and amplitude errors that are ultimately solved are also within the expected range.
In summary, according to the phased array radar amplitude phase consistency detection method provided by the embodiment of the application, the N probes are arranged in the preset range, and the preset range is determined according to the condition that the M vibrators and the N probes are mutually in the radiation range of each other; acquiring initial data of M vibrators and N probes, wherein the initial data comprise coordinate positions and pointing angles of the M vibrators, radiation patterns of the M vibrators and the N probes and coordinate positions of a phase center; further, N probes are utilized to respectively measure M vibrators one by one, and an amplitude set and a phase set are obtained; further, according to the amplitude set, the phase set and the initial data, the amplitude and the phase of the M vibrators, and the coordinate positions and the pointing angles of the N probes are determined by using a least square method. Phase uncertainty caused by cable motion is avoided, and the purpose of assisting the phased array radar to keep amplitude-phase consistency is directly achieved.
Based on the same technical concept, the embodiment of the application also provides a phased array radar amplitude phase consistency detection system, as shown in fig. 4, wherein the system comprises:
the probe setting module 401 is configured to set N probes within a preset range, where the preset range is determined according to a condition that M vibrators and N probes are located in a radiation range of each other, and M and N are integers greater than one.
The initial data acquisition module 402 is configured to acquire initial data of the M vibrators and the N probes, where the initial data includes coordinate positions and pointing angles of the M vibrators, and radiation patterns and phase center coordinate positions of the M vibrators and the N probes.
A measurement module 403 for measuring M vibrators with N probes one by one to obtain an amplitude set A i,j And phase set P i,j The method comprises the steps of carrying out a first treatment on the surface of the Where i is an integer from 1 to N and j is an integer from 1 to M.
An amplitude-phase data calculation module 404 for calculating an amplitude set A according to the amplitude set A i,j The phase set P i,j And determining the amplitude and the phase of the M vibrators, the coordinate positions and the pointing angles of the N probes by utilizing the initial data and utilizing a least square method.
In one possible implementation, the amplitude-phase data calculation module 404 specifically performs the calculation according to the foregoing formulas (1) and (2). The present embodiment is not described herein.
In one possible implementation, the measurement module 403 is specifically configured to:
respectively measuring M vibrators by using N probes at different positions one by one; or alternatively
And the single probe is moved to different positions to respectively measure M vibrators.
In one possible embodiment, the M transducers and the N probes are within the radiation range of each other, provided that the position of each probe is within the set beam width of all or part of the transducers, and the probe orientation satisfies that all or part of the transducers are within the set beam width of the probe.
In one possible embodiment, the position of each probe further includes: the shortest distance from each probe to the vibrator array is more than 2d 2 λ, where d represents the vibrator radiation aperture and λ represents the test wavelength.
In the present specification, each embodiment of the method is described in a progressive manner, and identical and similar parts of each embodiment are referred to each other, and each embodiment mainly describes differences from other embodiments. For relevance, see the description of the method embodiments.
It should be noted that although the operations of the method of the present application are depicted in the drawings in a particular order, this does not require or imply that the operations be performed in that particular order or that all illustrated operations be performed to achieve desirable results. Additionally or alternatively, certain steps may be omitted, multiple steps combined into one step to perform, and/or one step decomposed into multiple steps to perform.
Although the application provides method operational steps as an example or a flowchart, more or fewer operational steps may be included based on conventional or non-inventive means. The order of steps recited in the embodiments is merely one way of performing the order of steps and does not represent a unique order of execution. When implemented by an apparatus or client product in practice, the methods illustrated in the embodiments or figures may be performed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment, or even in a distributed data processing environment). The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, it is not excluded that additional identical or equivalent elements may be present in a process, method, article, or apparatus that comprises a described element.
The units, devices or modules etc. set forth in the above embodiments may be implemented in particular by a computer chip or entity or by a product having a certain function. For convenience of description, the above devices are described as being functionally divided into various modules, respectively. Of course, when implementing the present application, the functions of each module may be implemented in the same or multiple pieces of software and/or hardware, or a module implementing the same function may be implemented by multiple sub-modules or a combination of sub-units. The above-described apparatus embodiments are merely illustrative, for example, the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, for example, multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
Those skilled in the art will also appreciate that, in addition to implementing the controller in a pure computer readable program code, it is well possible to implement the same functionality by logically programming the method steps such that the controller is in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers, etc. Such a controller can be regarded as a hardware component, and means for implementing various functions included therein can also be regarded as a structure within the hardware component. Or even means for achieving the various functions may be regarded as either software modules implementing the methods or structures within hardware components.
The application may be described in the general context of computer-executable instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, classes, etc. that perform particular tasks or implement particular abstract data types. The application may also be practiced in distributed computing environments where tasks are performed by remote processing devices that are linked through a communications network. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
From the above description of embodiments, it will be apparent to those skilled in the art that the present application may be implemented in software plus a necessary general hardware platform. Based on such understanding, the technical solution of the present application may be embodied essentially or in a part contributing to the prior art in the form of a software product, which may be stored in a storage medium, such as a ROM/RAM, a magnetic disk, an optical disk, etc., including several instructions for causing a computer device (which may be a personal computer, a mobile terminal, a server, or a network device, etc.) to execute the method according to the embodiments or some parts of the embodiments of the present application.
Various embodiments in this specification are described in a progressive manner, and identical or similar parts are all provided for each embodiment, each embodiment focusing on differences from other embodiments. The application is operational with numerous general purpose or special purpose computer system environments or configurations. For example: personal computers, server computers, hand-held or portable devices, tablet devices, multiprocessor systems, microprocessor-based systems, set top boxes, programmable electronic devices, network PCs, minicomputers, mainframe computers, distributed computing environments that include any of the above systems or devices, and the like.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the application, and is not meant to limit the scope of the application, but to limit the application to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the application are intended to be included within the scope of the application.
Claims (8)
1. A phased array radar amplitude phase consistency detection method, the method comprising:
setting N probes in a preset range, wherein the preset range is determined according to the condition that M vibrators and N probes are mutually positioned in the radiation range of each other, and M and N are integers larger than 1;
acquiring initial data of M vibrators and N probes, wherein the initial data comprise coordinate positions and pointing angles of the M vibrators, radiation patterns of the M vibrators and the N probes and coordinate positions of a phase center;
n probes are used for measuring M vibrators one by one to obtain an amplitude set A i,j And phase set P i,j The method comprises the steps of carrying out a first treatment on the surface of the Wherein i is an integer from 1 to N, j is an integer from 1 to M;
according to the amplitude set A i,j The phase set P i,j And the initial data, M are determined by using a least square methodAmplitude and phase of the vibrator, coordinate positions and pointing angles of the N probes;
the method for determining the amplitude and the phase of M vibrators and the coordinate positions and the pointing angles of N probes by using the least square method is calculated according to the following formula:
A i,j ＝Ar j Ap i abs[Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j )]/|r i,j |
P i,j ＝mod[Pr j +P p i +angle(Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j ))-k|r i,j |,2π]
wherein Ar is j Representing the amplitude of vibrator j, ap i Representing the amplitude, θ, of the probe i j Represents the pointing angle, eta of vibrator j i Representing the pointing angle of the probe i; zeta type toy i,j Representing the corresponding angle coordinates of the connecting line between the two; dr (θ) j ,ξ i,j ) Representing the value of the radiation pattern of vibrator j at probe i, dp (eta) i ,ξ i,j ) Representing the value of the directional diagram of the probe i at the vibrator j, wherein the product of the directional diagram and the vibrator represents a three-dimensional vector inner product under the condition of an electromagnetic phased array and represents a common product under the condition of the acoustic phased array; pr (Pr) j Representing the phase, pp, of vibrator j i Representing the phase of probe i; r is |r i,j The i represents the distance between the phase center of probe i and the phase center of vibrator j; k represents wave number; abs () represents a modulo function, angle () represents a phase function, and mod () represents a remainder function;
on the premise that the probe and the oscillator directional diagrams as well as the position coordinates and the pointing angles of the oscillators are known, the amplitude, the phase and the position and the pointing angles of the oscillators can be solved as unknowns through a least square method.
2. The method of claim 1, wherein the measuring M transducers individually with N probes comprises:
respectively measuring M vibrators by using N probes at different positions one by one; or alternatively
A single probe is used to move to different positions to measure M vibrators separately.
3. The method of claim 1, wherein the M transducers and the N probes are within the radiation range of each other provided that the position of each probe is within a set beam width of all or part of the transducers and the probe orientation satisfies that all or part of the transducers are within the set beam width of the probe.
4. The method of claim 3, wherein the location of each probe further comprises: the shortest distance from each probe to the vibrator array is more than 2d 2 λ, where d represents the vibrator radiation aperture and λ represents the test wavelength.
5. A phased array radar amplitude phase consistency detection system, the system comprising:
the probe setting module is used for setting N probes in a preset range, the preset range is determined according to the condition that M vibrators and N probes are mutually positioned in the radiation range of each other, and M and N are integers larger than 1;
the initial data acquisition module is used for acquiring initial data of the M vibrators and the N probes, wherein the initial data comprise coordinate positions and pointing angles of the M vibrators, radiation patterns of the M vibrators and the N probes and phase center coordinate positions;
the measuring module is used for measuring M vibrators by utilizing N probes one by one to obtain an amplitude set A i,j And phase set P i,j The method comprises the steps of carrying out a first treatment on the surface of the Wherein i is an integer from 1 to N, j is an integer from 1 to M;
an amplitude-phase data calculation module for calculating the amplitude set A i,j The phase set P i,j The initial data are used for determining the amplitude and the phase of M vibrators and the coordinate positions and the pointing angles of N probes by using a least square method;
the method for determining the amplitude and the phase of M vibrators and the coordinate positions and the pointing angles of N probes by using the least square method is calculated according to the following formula:
A i,j ＝Ar j Ap i abs[Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j )]/|r i,j |
P i,j ＝mod[Pr j +P p i +angle(Dr(θ j ,ξ i,j )Dp(η i ,ξ i,j ))-k|r i,j |,2π]
wherein Ar is j Representing the amplitude of vibrator j, ap i Representing the amplitude, θ, of the probe i j Represents the pointing angle, eta of vibrator j i Representing the pointing angle of the probe i; zeta type toy i,j Representing the corresponding angle coordinates of the connecting line between the two; dr (θ) j ,ξ i,j ) Representing the value of the radiation pattern of vibrator j at probe i, dp (eta) i ,ξ i,j ) Representing the value of the directional diagram of the probe i at the vibrator j, wherein the product of the directional diagram and the vibrator represents a three-dimensional vector inner product under the condition of an electromagnetic phased array and represents a common product under the condition of the acoustic phased array; pr (Pr) j Representing the phase, pp, of vibrator j i Representing the phase of probe i; r is |r i,j The i represents the distance between the phase center of probe i and the phase center of vibrator j; k represents wave number; abs () represents a modulo function, angle () represents a phase function, and mod () represents a remainder function;
on the premise that the probe and the oscillator directional diagrams as well as the position coordinates and the pointing angles of the oscillators are known, the amplitude, the phase and the position and the pointing angles of the oscillators can be solved as unknowns through a least square method.
6. The system of claim 5, wherein the measurement module is specifically configured to:
respectively measuring M vibrators by using N probes at different positions one by one; or alternatively
A single probe is used to move to different positions to measure M vibrators separately.
7. The system of claim 5, wherein the M transducers and the N probes are within the radiation range of each other provided that the position of each probe is within the set beamwidth of all or a portion of the transducers and the probe orientation satisfies that all or a portion of the transducers are within the set beamwidth of the probe.
8. The system of claim 7, wherein the location of each probe further comprises: the shortest distance from each probe to the vibrator array is more than 2d 2 λ, where d represents the vibrator radiation aperture and λ represents the test wavelength.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911128100.1A CN111965602B (en) | 2019-11-18 | 2019-11-18 | Phased array radar amplitude-phase consistency detection method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911128100.1A CN111965602B (en) | 2019-11-18 | 2019-11-18 | Phased array radar amplitude-phase consistency detection method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111965602A CN111965602A (en) | 2020-11-20 |
| CN111965602B true CN111965602B (en) | 2023-11-10 |
Family
ID=73358344
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911128100.1A Active CN111965602B (en) | 2019-11-18 | 2019-11-18 | Phased array radar amplitude-phase consistency detection method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111965602B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112881992B (en) * | 2021-02-25 | 2022-02-11 | 四川九洲电器集团有限责任公司 | Channel amplitude and phase consistency self-checking method and device and computer equipment |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107064887A (en) * | 2016-12-26 | 2017-08-18 | 中国电子科技集团公司第二十研究所 | Planar Phased Array Antenna is quickly measured and automatic calibrating method |
| CN107765104A (en) * | 2017-09-04 | 2018-03-06 | 华为技术有限公司 | The method and school that a kind of phased array school is surveyed survey device |
| CN109309533A (en) * | 2018-09-04 | 2019-02-05 | 华为技术有限公司 | A kind of calibration method and equipment |
| CN109541330A (en) * | 2018-11-12 | 2019-03-29 | 北京航空航天大学 | A kind of array antenna channel calibration system of plane-wave simulation device |
| CN109581279A (en) * | 2018-11-06 | 2019-04-05 | 中国电子科技集团公司第三十八研究所 | A kind of bearing calibration of ultra wide band multiple-beam system direction finding precision and device |
-
2019
- 2019-11-18 CN CN201911128100.1A patent/CN111965602B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107064887A (en) * | 2016-12-26 | 2017-08-18 | 中国电子科技集团公司第二十研究所 | Planar Phased Array Antenna is quickly measured and automatic calibrating method |
| CN107765104A (en) * | 2017-09-04 | 2018-03-06 | 华为技术有限公司 | The method and school that a kind of phased array school is surveyed survey device |
| CN109309533A (en) * | 2018-09-04 | 2019-02-05 | 华为技术有限公司 | A kind of calibration method and equipment |
| CN109581279A (en) * | 2018-11-06 | 2019-04-05 | 中国电子科技集团公司第三十八研究所 | A kind of bearing calibration of ultra wide band multiple-beam system direction finding precision and device |
| CN109541330A (en) * | 2018-11-12 | 2019-03-29 | 北京航空航天大学 | A kind of array antenna channel calibration system of plane-wave simulation device |
Non-Patent Citations (2)
| Title |
|---|
| GPU-accelerated spherical near-to-far field transformation software for antenna near-field measurements;Yi S. Ding;《 2019 International Workshop on Electromagnetics: Applications and Student Innovation Competition (iWEM)》;1-3 * |
| 相控阵天线平面近场误差分析与校准方法研究;焦禹;《中国优秀硕士学位论文全文数据库 信息科技辑》;第2017卷(第2期);I136-128 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111965602A (en) | 2020-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11309977B2 (en) | Calibration method and communications device | |
| CN107329003B (en) | Optimize the method for SAR Antenna Pattern Measuring Test | |
| Joy et al. | Applications of probe-compensated near-field measurements | |
| CN108445303B (en) | Near-field electromagnetic scattering characteristic simulation method | |
| CN107783086B (en) | Method for diagnosing distorted position of antenna array aperture amplitude phase field | |
| CN107037418B (en) | A kind of bounded space bearing calibration of the super directive property multipole vector array of low frequency | |
| CN104391183A (en) | Near-field-measurement-based rapid calculation method of antenna far-field region characteristic | |
| CN113381187B (en) | Spherical phased array antenna coordinate far and near field comparison and correction method | |
| KR101979770B1 (en) | Apparatus and method for simulating beam-focusing of near-field in aesa radar | |
| CN104375129B (en) | A kind of distributive array coherent synthesis radar emission coherent parameter calibrating method | |
| CN105606906A (en) | Millimeter wave phased array test and calibration method | |
| CN112147593B (en) | Four-dimensional parameter estimation method for high-speed dense explosive fragment target | |
| CN109813969A (en) | A kind of diagnostic method of array antenna, equipment, system and computer readable storage medium | |
| CN117214841A (en) | Phased array radar antenna plane calibration method | |
| CN112711040A (en) | Satellite navigation antenna performance evaluation system and method | |
| CN104200110A (en) | Decoupling-based two-dimensional beam scanning angle measurement method | |
| CN116520035A (en) | Two-dimensional pattern testing method adopting plane near-field rapid inversion method | |
| CN111965602B (en) | Phased array radar amplitude-phase consistency detection method and system | |
| CN113419208B (en) | Interference source real-time direction finding method and device, storage medium and electronic equipment | |
| CN104852775A (en) | Rapid evaluation method for large phased-array antenna wave beam performance | |
| CN117031418A (en) | SAR satellite azimuth agility observation mode scanning and pointing detection method and system | |
| CN109412710B (en) | Antenna transmission performance evaluation method and device | |
| CN109815509B (en) | Diagnostic method, device, system and computer readable storage medium for array antenna | |
| CN109921865A (en) | A kind of the calibration bar approximate simulation calibration system and method for full airspace phased array antenna | |
| CN115542268A (en) | Large-aperture phased array antenna block testing method and system |
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 |