CN112147588B - Rapid calculation method for asymmetric radar irradiation area - Google Patents

Rapid calculation method for asymmetric radar irradiation area Download PDF

Info

Publication number
CN112147588B
CN112147588B CN202011092555.5A CN202011092555A CN112147588B CN 112147588 B CN112147588 B CN 112147588B CN 202011092555 A CN202011092555 A CN 202011092555A CN 112147588 B CN112147588 B CN 112147588B
Authority
CN
China
Prior art keywords
zone
irradiation area
equation
antenna
area
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
Application number
CN202011092555.5A
Other languages
Chinese (zh)
Other versions
CN112147588A (en
Inventor
张浙东
张玉石
黎鑫
许心瑜
李清亮
夏晓云
张金鹏
赵鹏
万晋通
尹志盈
朱秀芹
余运超
李善斌
尹雅磊
李慧明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute of Radio Wave Propagation CETC 22 Research Institute
Original Assignee
China Institute of Radio Wave Propagation CETC 22 Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Institute of Radio Wave Propagation CETC 22 Research Institute filed Critical China Institute of Radio Wave Propagation CETC 22 Research Institute
Priority to CN202011092555.5A priority Critical patent/CN112147588B/en
Publication of CN112147588A publication Critical patent/CN112147588A/en
Application granted granted Critical
Publication of CN112147588B publication Critical patent/CN112147588B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01S7/4004—Means for monitoring or calibrating of parts of a radar system
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Abstract

The invention discloses a method for quickly calculating an asymmetric radar irradiation area, which comprises the following steps: step 1, establishing an asymmetric transmitting and receiving antenna irradiation area equation: and 2, establishing a self-adaptive cross irradiation area rapid numerical calculation method. The invention provides a method for quickly calculating asymmetrical radar irradiation areas, which comprises the steps of establishing a rectangular coordinate system by taking a radar receiving antenna as a reference, then respectively establishing radar irradiation area equations with unequal long axes and equal short axes of a near area and a far area of the radar receiving and transmitting antenna according to the position relation erected by the radar receiving and transmitting antenna, incident beams and azimuth beams, finally obtaining calculation step lengths corresponding to different incident angles through a step length correction method, and solving the cross irradiation areas of the receiving and transmitting antenna corresponding to different incident angles. Compared with the original method, the method improves the calculation speed and the calculation precision of the radar irradiation area of the dual-antenna system, and increases the collection amount of effective radar echo data in unit time.

Description

Rapid calculation method for asymmetric radar irradiation area
Technical Field
The invention belongs to the field of radar measurement, and particularly relates to a method for quickly calculating radar irradiation area for measuring radar receiving and transmitting antenna separation in the field.
Background
The radar is generally divided into a pulse system and a continuous wave system in terms of system, wherein the pulse system radar adopts a mode of sharing a transmitting and receiving antenna, and the continuous wave system radar adopts a mode of separating the transmitting and receiving antenna. When the radar collects echo data to the ground and sea surface background, the data is generated by each scattering element in the irradiation area of the radar background. The method for calculating the irradiation area is relatively simple for a pulse system, and when the depression angle is small, the irradiation area is calculated
Figure GDA0003479001310000011
Wherein R and phiaRespectively the distance from the radar to a target and the width of an antenna beam on an azimuth plane, c is the propagation speed of electromagnetic waves in vacuum, tau is the pulse width, and theta is a depression angle; when the depression angle is close to 90 degrees, the irradiation area
Figure GDA0003479001310000012
Wherein phicIs the width of the antenna beam in the elevation plane. For a radar with separate transmitting and receiving antennas, the illumination area belongs to the intersection area of the illumination areas of the transmitting and receiving antennas, as shown in fig. 1, O1And O2Respectively belong to the areas where the central beams of the transmitting antenna and the receiving antenna are located, O1-T1Belonging to the nearly 3dB spot beam irradiation area of the transmitting antenna, O1-T2Belongs to the far 3dB spot beam irradiation area of the transmitting antenna, O2-R1Belonging to a receiving antenna3dB spot beam irradiation area, O2-R2Belonging to the far 3dB spot beam irradiation area of the receiving antenna. As can be seen from fig. 1, the irradiation area of the near 3dB spot beam is smaller than the irradiation area of the far 3dB spot beam, and the difference increases with the increase of the incident angle, and similarly, the irradiation area of the intersection region of the transmitting antenna and the receiving antenna also has a similar phenomenon, and the irradiation area O-M of the near 3dB spot intersection region is obviously smaller than the irradiation area O-N of the far 3dB spot intersection region.
In the process of calculating the irradiation area of the intersection region, the conventional method is to respectively and equivalently equalize the irradiation area of the transmitting antenna and the irradiation area of the receiving antenna into an ellipse, and the semimajor axis and the semiminor axis of the ellipse of the transmitting antenna are respectively aTAnd bTThe semi-major axis and the semi-minor axis of the receiving antenna are respectively aRAnd bR
Figure GDA0003479001310000021
bT＝O1A
Figure GDA0003479001310000022
bR＝O1B
Wherein, aT＝aR，bT＝bR. Then, the irradiation area of the cross region is determined by an integral method, see the related articles of "derivation and calculation of radar irradiation area" from "radio wave and antenna" at 2 nd stage 1983 and "calculation of scatterometer equivalent irradiation area" from "radio wave science bulletin" at 8 st stage 1993 ". The problems of this method are: firstly, the cross irradiation area of the double antennas obtained by calculation is larger than that of the actual situation; secondly, an integral mode is adopted, so that the calculation real-time performance is poor, the sizes of the irradiation areas of the cross areas corresponding to different heights, different incidence angles and different wave bands need to be calculated in advance, and then a table look-up mode is adopted for use.
Disclosure of Invention
The invention aims to solve the technical problem of providing a method for quickly calculating the asymmetrical radar irradiation area with high precision and high speed.
The invention adopts the following technical scheme:
the improvement of the method for quickly calculating the irradiation area of the asymmetrical radar is that the method comprises the following steps:
step 1, establishing an asymmetric transmitting and receiving antenna irradiation area equation:
step 11, transmitting antenna irradiation area equation:
step 111, calculating semi-major axis a of an elliptic equation of near zone irradiation area of the transmitting antennajT
Figure GDA0003479001310000023
Step 112, calculating the semimajor axis a of the ellipse equation of the irradiation area of the far zone of the transmitting antennayT
Figure GDA0003479001310000024
Step 113, calculating the semiminor axis b of the ellipse equation of the near zone and the far zone of the transmitting antennaT：
bT＝H·tan(θT)/cos(β) (4)
Step 114, writing an ellipse equation of the near zone and the far zone of the transmitting antenna, and the central point x of the ellipse equation of the near zone and the far zone of the transmitting antenna0＝0,y0＝T0R0L, where L is the erection horizontal distance between the receiving antenna and the transmitting antenna, and the elliptic equation is as follows:
near zone ellipse equation:
Figure GDA0003479001310000031
wherein x is less than or equal to 0 and 5;
far zone ellipse equation:
Figure GDA0003479001310000032
wherein x is more than or equal to 0 and 6;
step 12, receiving an antenna irradiation area equation:
step 121, calculating a semimajor axis a of an elliptic equation of the near-zone irradiation area of the receiving antennajR：
Figure GDA0003479001310000033
Step 122, calculating the semimajor axis a of the receiving antenna far zone irradiation area elliptic equationyR；
Figure GDA0003479001310000034
Step 123, calculating the semi-minor axis b of the ellipse equation of the near zone and the far zone of the receiving antennaR；
bR＝H·tan(θR)/cos(β) (9)
Step 124, writing out an ellipse equation of the near zone and the far zone of the receiving antenna, and the central point x of the ellipse equation of the near zone and the far zone of the receiving antenna0＝0,y0The ellipse equation is 0 as follows:
near zone ellipse equation:
Figure GDA0003479001310000035
wherein x is less than or equal to 0 (10);
far zone ellipse equation:
Figure GDA0003479001310000041
wherein x is more than or equal to 0 and 11;
because the receiving antenna and the transmitting antenna are the same in performance index and are consistent in erection height, the receiving antenna and the transmitting antenna have the following properties:
Figure GDA0003479001310000042
step 2, establishing a self-adaptive cross illumination area rapid numerical calculation method:
step 21, setting a frame for radar background measurementSetting height H, incident angle beta, distance L between antennas, and half beam width of incident angle
Figure GDA0003479001310000043
Azimuth half beam width theta, wherein the number of incident angles is greater than or equal to 4, and the range of the incident angles is
Figure GDA0003479001310000044
Step 22, calculating the length of the cross area MN according to the parameters set in step 21
Figure GDA0003479001310000045
While
Figure GDA0003479001310000046
Near-zone elliptical crossing region L of transmitting-receiving antennajMNElliptical cross region L with far zone of transmitting-receiving antennayMNThe composition is as follows:
Figure GDA0003479001310000047
LjMNcan be calculated from the equations (5) and (10), and the expression is as follows:
Figure GDA0003479001310000048
LyMNcan be calculated from equations (6) and (11), and the expression is as follows:
Figure GDA0003479001310000049
from formula (15) and formula (16)
Figure GDA00034790013100000410
A series of enlarged LMNValue of
LMN1,LMN2,...,LMNn-1,LMNn，n≥4 (17)
Step 23, determining the maximum step length m according to the maximum crossing area corresponding to the maximum incident anglemsnThe step length corresponding to other incident angles is based on the correction function f (x), and the numerical calculation step length m of a certain incident angle is obtained in a mapping modesnThe correction function f (x) is related to the incident angle and the size of the cross area, and the correction factor is calculated as follows:
f(x)＝1/(1+k*exp(-x/10)) (1)
wherein x represents an incident angle, n is total, n is not less than 4, and k is (MN (x)n-1)-MN(x2))/(xn-1-x2) N is more than or equal to 4, and MN represents the length of the cross area in the X-axis direction;
the slope k in equation (1) is obtained as follows:
Figure GDA0003479001310000051
step 24, according to the equations (1) and (18), the step length of the cross-illumination area numerical calculation corresponding to different incident angles can be obtained, assuming that L isMNMaximum step number of mmsnThen the number of steps m for different incident anglessnThe formula (c) is as follows:
Figure GDA0003479001310000052
step 25, obtaining the step number corresponding to different incident angles according to the formula (19)
Figure GDA0003479001310000053
Dx and dy in (1), the formula is as follows:
Figure GDA0003479001310000054
formula (II)
Figure GDA0003479001310000055
In, XN、XMIs represented by the following formula:
Figure GDA0003479001310000061
Figure GDA0003479001310000062
and Y isR、YTRelating to the irradiation area of the transmitting and receiving antennas, divided into a near zone and a far zone,
Figure GDA0003479001310000063
the near zone is represented by the formula:
Figure GDA0003479001310000064
in the formula (I), the compound is shown in the specification,
Figure GDA0003479001310000065
Figure GDA0003479001310000066
the distal region is represented by the formula:
Figure GDA0003479001310000067
in the formula (I), the compound is shown in the specification,
Figure GDA0003479001310000068
the invention has the beneficial effects that:
the invention provides a method for quickly calculating asymmetrical radar irradiation areas, which comprises the steps of establishing a rectangular coordinate system by taking a radar receiving antenna as a reference, respectively establishing a radar irradiation area equation with unequal semi-major axes and equal semi-minor axes of a near area and a far area of the radar receiving and transmitting antenna according to the position relation of the erection of the radar receiving and transmitting antenna, obtaining the calculation step lengths corresponding to different incidence angles by a step length correction method, and solving the cross irradiation areas of the receiving and transmitting antenna corresponding to different incidence angles. Compared with the original method, the method improves the calculation speed and the calculation precision of the radar irradiation area of the dual-antenna system, and increases the collection amount of effective radar echo data in unit time.
Drawings
FIG. 1 is a top view of a dual antenna crossover region;
FIG. 2 is a schematic diagram of a radar illuminated area with separate transmit and receive antennas;
FIG. 3 is a schematic diagram of a rectangular coordinate system of the cross area of the dual antenna;
FIG. 4 is a schematic flow chart of building transmit receive antenna ellipse equations;
fig. 5 is a graph showing the variation of MN length with incident angle when H is 15 m;
fig. 6 is a graph showing the variation of MN length with incident angle when H is 10 m;
FIG. 7 is a schematic diagram of adaptive step size as a function of incident angle;
FIG. 8 is a graph showing the variation of the illuminated area with the incident angle.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Embodiment 1, this embodiment discloses a method for quickly calculating an asymmetric radar irradiation area, and the improvement is that the method includes the following steps:
step 1, establishing an asymmetric transmitting and receiving antenna irradiation area equation:
the irradiation area of the dual-antenna system radar is the cross area of the irradiation area of the transmitting antenna and the irradiation area of the receiving antenna, the near-zone irradiation area shape of the antenna is different from the far-zone irradiation area shape of the antenna, and the far-zone irradiation area is far larger than the near-zone irradiation area along with the increase of the incident angle.
As shown in fig. 2, the radar transmits a signal from the transmitting antenna T, and after the signal is reflected by each scattering unit in the irradiation area of the cross region, the signal is transmitted back to the receiving antenna R, so that the irradiation area of the cross region must exist on the premise that the receiving antenna R can receive the signal. First, the irradiation areas of the transmitting antenna T and the receiving antenna R are divided into a near-zone irradiation area and a far-zone irradiation area, and the near-zone irradiation area and the far-zone irradiation area are respectively approximated to two ellipses having different semi-major axes and equal semi-minor axes. A near zone irradiation area elliptic equation and a far zone irradiation area elliptic equation are established for the transmitting antenna T and the receiving antenna R in an asymmetric mode, and a schematic diagram of a double-antenna cross area rectangular coordinate system is shown in fig. 3. The semimajor axis and semiminor axis of the ellipse can pass through the erection height H of the antenna, the incidence angle beta of the antenna, the distance L between the antennas and the half-beam width of the incidence direction of the antenna
Figure GDA0003479001310000071
And calculating parameters such as the azimuth angle direction half beam width theta of the antenna. The process of establishing the ellipse equation is shown in figure 4. The method specifically comprises the following steps:
step 11, transmitting antenna irradiation area equation:
step 111, calculating semi-major axis a of an elliptic equation of near zone irradiation area of the transmitting antennajT
Figure GDA0003479001310000072
Step 112, calculating the semimajor axis a of the ellipse equation of the irradiation area of the far zone of the transmitting antennayT
Figure GDA0003479001310000081
Step 113, calculating the semiminor axis b of the ellipse equation of the near zone and the far zone of the transmitting antennaT：
bT＝H·tan(θT)/cos(β) (4)
Step 114, write the near zone of the transmitting antennaEllipse equation with far zone, as can be seen from fig. 3, the central point x of the ellipse equation of the near zone and the far zone of the transmitting antenna0＝0,y0＝T0R0L, where L is the erection horizontal distance between the receiving antenna and the transmitting antenna, and the elliptic equation is as follows:
near zone ellipse equation:
Figure GDA0003479001310000082
wherein x is less than or equal to 0 and 5;
far zone ellipse equation:
Figure GDA0003479001310000083
wherein x is more than or equal to 0 and 6;
step 12, receiving an antenna irradiation area equation:
step 121, calculating a semimajor axis a of an elliptic equation of the near-zone irradiation area of the receiving antennajR：
Figure GDA0003479001310000084
Step 122, calculating the semimajor axis a of the receiving antenna far zone irradiation area elliptic equationyR；
Figure GDA0003479001310000085
Step 123, calculating the semi-minor axis b of the ellipse equation of the near zone and the far zone of the receiving antennaR；
bR＝H·tan(θR)/cos(β) (9)
Step 124, writing out the ellipse equation of the near zone and the far zone of the receiving antenna, as can be seen from fig. 3, the central point x of the ellipse equation of the near zone and the far zone of the receiving antenna0＝0,y0The ellipse equation is 0 as follows:
near zone ellipse equation:
Figure GDA0003479001310000091
wherein x is less than or equal to 0 (10);
far zone ellipse equation:
Figure GDA0003479001310000092
wherein x is more than or equal to 0 and 11;
because the receiving antenna and the transmitting antenna are the same in performance index and are consistent in erection height, the receiving antenna and the transmitting antenna have the following properties:
Figure GDA0003479001310000093
step 2, establishing a self-adaptive cross illumination area rapid numerical calculation method:
as shown in FIG. 3, the shadow area between M-N is the required irradiation area of the dual antenna, and its size depends on the height H of the antenna and the half-beam width of the incident direction of the antenna
Figure GDA0003479001310000094
The difference in the azimuth direction half beam width θ of the antenna changes. H15 m, L1 m, theta 3 deg,
Figure GDA0003479001310000095
For example, as shown in fig. 5, when the incident angle is less than 60 degrees, MN does not change much, and when the incident angle exceeds 60 degrees, MN length increases sharply, and the irradiation area also increases sharply. The irradiation area is an irregular area, and the size of the area is difficult to be calculated by an equation, but can be calculated by numerical method, namely
Figure GDA0003479001310000096
The method for calculating the irradiation area specifically comprises the following steps:
step 21, setting the erection height H, the incidence angle beta, the distance L between the antennas and the half-beam width of the incidence angle of the radar background measurement
Figure GDA0003479001310000097
Azimuth half beam width θWherein the number of the incident angles is more than or equal to 4, and the range of the incident angles is
Figure GDA0003479001310000098
Step 22, calculating the length of the cross area MN according to the parameters set in step 21
Figure GDA0003479001310000099
While
Figure GDA00034790013100000910
Near-zone elliptical crossing region L of transmitting-receiving antennajMNElliptical cross region L with far zone of transmitting-receiving antennayMNThe composition is as follows:
Figure GDA0003479001310000101
LjMNcan be calculated from the equations (5) and (10), and the expression is as follows:
Figure GDA0003479001310000102
LyMNcan be calculated from equations (6) and (11), and the expression is as follows:
Figure GDA0003479001310000103
from formula (15) and formula (16)
Figure GDA0003479001310000104
A series of enlarged LMNValue of
LMN1,LMN2,...,LMNn-1,LMNn，n≥4 (17)
Step 23, since the irradiation area is a dynamic process with different measurement parameters, a difficulty in engineering implementation is how to divide the step length of the numerical calculation, i.e. dx anddy, too fine granularity of division will result in long calculation time and low efficiency, too coarse granularity of division, high calculation efficiency, but the calculation accuracy will be reduced. Therefore, a self-adaptive cross-illumination area fast numerical calculation method is provided, and the principle of the method is that the maximum step length m is determined according to the maximum cross-illumination area corresponding to the maximum incident anglemsnThe step length corresponding to other incident angles is based on the correction function f (x), and the numerical calculation step length m of a certain incident angle is obtained in a mapping modesnThe correction function f (x) is related to the incident angle and the size of the cross area, and the correction factor is calculated as follows:
f(x)＝1/(1+k*exp(-x/10)) (1)
wherein x represents an incident angle, n is total, n is not less than 4, and k is (MN (x)n-1)-MN(x2))/(xn-1-x2) N is not less than 4, x is takenn-1And x2Eliminating the influence of the maximum value and the minimum value, wherein MN represents the length of the cross area in the X-axis direction;
by establishing a step length correction function, the step length of numerical calculation corresponding to different incidence angles is changed, and under the condition that the calculation precision is not changed, the calculation efficiency and the calculation real-time performance are improved.
The slope k in equation (1) is obtained as follows:
Figure GDA0003479001310000111
step 24, according to the equations (1) and (18), the step length of the cross-illumination area numerical calculation corresponding to different incident angles can be obtained, assuming that L isMNMaximum step number of mmsnThen the number of steps m for different incident anglessnThe formula (c) is as follows:
Figure GDA0003479001310000112
step 25, obtaining the step number corresponding to different incident angles according to the formula (19)
Figure GDA0003479001310000113
Dx and dy in (1), the formula is as follows:
Figure GDA0003479001310000114
formula (II)
Figure GDA0003479001310000115
In, XN、XMIs represented by the following formula:
Figure GDA0003479001310000116
Figure GDA0003479001310000117
and Y isR、YTRelating to the irradiation area of the transmitting and receiving antennas, divided into a near zone and a far zone,
Figure GDA0003479001310000118
the near zone is represented by the formula:
Figure GDA0003479001310000119
in the formula (I), the compound is shown in the specification,
Figure GDA00034790013100001110
Figure GDA0003479001310000121
the distal region is represented by the formula:
Figure GDA0003479001310000122
in the formula (I), the compound is shown in the specification,
Figure GDA0003479001310000123
therefore, the specific implementation process of the calculation of the irradiation area of the asymmetric radar is completed.
The height H of the radar is 10 m, the horizontal distance L between the double antennas is 1.2 m, and the half-wave beam width in the incident direction of the transmitting and receiving antenna
Figure GDA0003479001310000124
The half beam width θ in the azimuth direction is 3.5 °, the incident angle is 5 ° to 85 °, the interval is 5 °, and the maximum step size is 2000 as an example. The variation curve of the length of the irradiation area MN of the intersection region with the incident angle is shown in FIG. 6. The slope k is 0.6445, and the curve of the adaptive step size along with the incident angle is shown in fig. 7. The curve of the calculated variation of the irradiation area with the incident angle is shown in fig. 8. Fourthly, calculating efficiency: the step-size mode requires 68000000 calculation cycles; the adaptive step size approach requires 32500000 calculation cycles;
Figure GDA0003479001310000125
the improvement is more than 1 time.

Claims (1)

1. A method for rapidly calculating an asymmetrical radar irradiation area is characterized by comprising the following steps:
step 1, establishing an asymmetric transmitting and receiving antenna irradiation area equation:
step 11, transmitting antenna irradiation area equation:
step 111, calculating semi-major axis a of an elliptic equation of near zone irradiation area of the transmitting antennajT
Figure FDA0003479001300000011
Step 112, calculating the semimajor axis a of the ellipse equation of the irradiation area of the far zone of the transmitting antennayT
Figure FDA0003479001300000012
Step 113, calculating the semiminor axis b of the ellipse equation of the near zone and the far zone of the transmitting antennaT：
bT＝H·tan(θT)/cos(β) (4)
Step 114, writing an ellipse equation of the near zone and the far zone of the transmitting antenna, and the central point x of the ellipse equation of the near zone and the far zone of the transmitting antenna0＝0,y0＝T0R0L, where L is the erection horizontal distance between the receiving antenna and the transmitting antenna, and the elliptic equation is as follows:
near zone ellipse equation:
Figure FDA0003479001300000013
wherein x is less than or equal to 0 and 5;
far zone ellipse equation:
Figure FDA0003479001300000014
wherein x is more than or equal to 0 and 6;
step 12, receiving an antenna irradiation area equation:
step 121, calculating a semimajor axis a of an elliptic equation of the near-zone irradiation area of the receiving antennajR：
Figure FDA0003479001300000015
Step 122, calculating the semimajor axis a of the receiving antenna far zone irradiation area elliptic equationyR；
Figure FDA0003479001300000021
Step 123, calculating the semi-minor axis b of the ellipse equation of the near zone and the far zone of the receiving antennaR；
bR＝H·tan(θR)/cos(β) (9)
Step 124, writing the ellipse equation of the near zone and the far zone of the receiving antenna, the near zone and the far zone of the receiving antennaCentral point x of the far zone elliptic equation0＝0,y0The ellipse equation is 0 as follows:
near zone ellipse equation:
Figure FDA0003479001300000022
wherein x is less than or equal to 0 (10);
far zone ellipse equation:
Figure FDA0003479001300000023
wherein x is more than or equal to 0 and 11;
because the receiving antenna and the transmitting antenna are the same in performance index and are consistent in erection height, the receiving antenna and the transmitting antenna have the following properties:
Figure FDA0003479001300000024
step 2, establishing a self-adaptive cross illumination area rapid numerical calculation method:
step 21, setting the erection height H, the incidence angle beta, the distance L between the antennas and the half-beam width of the incidence angle of the radar background measurement
Figure FDA0003479001300000025
Azimuth half beam width theta, wherein the number of incident angles is greater than or equal to 4, and the range of the incident angles is
Figure FDA0003479001300000026
Step 22, calculating the length of the cross area MN according to the parameters set in step 21
Figure FDA0003479001300000027
While
Figure FDA0003479001300000028
By near-zone elliptical cross-over zone of transmitting and receiving antennasDomain LjMNElliptical cross region L with far zone of transmitting-receiving antennayMNThe composition is as follows:
Figure FDA0003479001300000031
LjMNcan be calculated from the equations (5) and (10), and the expression is as follows:
Figure FDA0003479001300000032
LyMNcan be calculated from equations (6) and (11), and the expression is as follows:
Figure FDA0003479001300000033
from the formulae (15) and (16), it is possible to obtain a compound which is obtained by reacting a compound of formula (I) with H, beta,
Figure FDA0003479001300000034
a series of L's where theta becomes largeMNValue of
LMN1,LMN2,...,LMNn-1,LMNn，n≥4 (17)
Step 23, determining the maximum step length m according to the maximum crossing area corresponding to the maximum incident anglemsnThe step length corresponding to other incident angles is based on the correction function f (x), and the numerical calculation step length m of a certain incident angle is obtained in a mapping modesnThe correction function f (x) is related to the incident angle and the size of the cross area, and the correction factor is calculated as follows:
f(x)＝1/(1+k*exp(-x/10)) (1)
wherein x represents an incident angle, n is total, n is not less than 4, and k is (MN (x)n-1)-MN(x2))/(xn-1-x2) N is more than or equal to 4, and MN represents the length of the cross area in the X-axis direction;
the slope k in equation (1) is obtained as follows:
Figure FDA0003479001300000035
step 24, according to the equations (1) and (18), the step length of the cross-illumination area numerical calculation corresponding to different incident angles can be obtained, assuming that L isMNMaximum step number of mmsnThen the number of steps m for different incident anglessnThe formula (c) is as follows:
Figure FDA0003479001300000041
step 25, obtaining the step number corresponding to different incident angles according to the formula (19)
Figure FDA0003479001300000042
Dx and dy in (1), the formula is as follows:
Figure FDA0003479001300000043
formula (II)
Figure FDA0003479001300000044
In, XN、XMIs represented by the following formula:
Figure FDA0003479001300000045
Figure FDA0003479001300000046
and Y isR、YTRelating to the irradiation area of the transmitting and receiving antennas, divided into a near zone and a far zone,
the near zone is represented by the formula:
Figure FDA0003479001300000047
in the formula (I), the compound is shown in the specification,
Figure FDA0003479001300000048
the distal region is represented by the formula:
Figure FDA0003479001300000049
in the formula (I), the compound is shown in the specification,
Figure FDA00034790013000000410
CN202011092555.5A 2020-10-14 2020-10-14 Rapid calculation method for asymmetric radar irradiation area Active CN112147588B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011092555.5A CN112147588B (en) 2020-10-14 2020-10-14 Rapid calculation method for asymmetric radar irradiation area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011092555.5A CN112147588B (en) 2020-10-14 2020-10-14 Rapid calculation method for asymmetric radar irradiation area

Publications (2)

Publication Number Publication Date
CN112147588A CN112147588A (en) 2020-12-29
CN112147588B true CN112147588B (en) 2022-03-01

Family

ID=73951708

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011092555.5A Active CN112147588B (en) 2020-10-14 2020-10-14 Rapid calculation method for asymmetric radar irradiation area

Country Status (1)

Country Link
CN (1) CN112147588B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102112892A (en) * 2007-05-29 2011-06-29 剑桥顾问 Radar system and method
CN107005092A (en) * 2014-09-05 2017-08-01 索雷斯能源公司 Radio field electric power transmission system, method and its transmitter and receiver
CN109270529A (en) * 2018-12-07 2019-01-25 电子科技大学 Forward sight array SAR high-resolution imaging method and system based on virtual-antenna
CN110736971A (en) * 2019-11-05 2020-01-31 中国电波传播研究所（中国电子科技集团公司第二十二研究所） real-time recognition method for non-cooperative targets in sea clutter measurement area of shore-based radar

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IT1275349B (en) * 1994-11-25 1997-08-05 Alenia Spazio Spa ANTENNA WITH ROTARY ELLIPTICAL BEAM WITH POSSIBILITY OF RECONFIGURATION AND BEAM ZOOM
JP2000241529A (en) * 1999-02-22 2000-09-08 Honda Motor Co Ltd Radar apparatus and abnormality detecting method for radar apparatus
US20060082989A1 (en) * 2004-10-18 2006-04-20 Yeong Long Electric Industry Co., Ltd. Asymmetric LED in-ground lamp
FR2967787B1 (en) * 2010-11-23 2012-12-14 Thales Sa METHOD OF ESTIMATING THE ANGULAR POSITION OF A TARGET BY RADAR AND RADAR DETECTION USING THE METHOD
WO2013013468A1 (en) * 2011-07-26 2013-01-31 深圳光启高等理工研究院 Offset feed radar antenna
CN112947793A (en) * 2015-05-22 2021-06-11 触觉实验室股份有限公司 Transmitting and receiving system and method for bidirectional orthogonal signaling sensor
CN108363046B (en) * 2018-02-09 2020-03-27 周口师范学院 Passive radar direction finding method and system based on double-receiver elliptic scanning
JP6968978B2 (en) * 2018-03-13 2021-11-24 古野電気株式会社 Precipitation particle discrimination device, precipitation particle discrimination system, precipitation particle discrimination method, and precipitation particle discrimination program
CN110596701B (en) * 2019-07-01 2021-09-28 杭州电子科技大学 Non-level-flight double-station SAR frequency domain FENLCS imaging method based on quadratic ellipse model
CN110954877B (en) * 2019-11-05 2022-04-26 中国电波传播研究所（中国电子科技集团公司第二十二研究所） Method for generating terrain of measurement area of scattering characteristics of vehicle-mounted radar

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102112892A (en) * 2007-05-29 2011-06-29 剑桥顾问 Radar system and method
CN107005092A (en) * 2014-09-05 2017-08-01 索雷斯能源公司 Radio field electric power transmission system, method and its transmitter and receiver
CN109270529A (en) * 2018-12-07 2019-01-25 电子科技大学 Forward sight array SAR high-resolution imaging method and system based on virtual-antenna
CN110736971A (en) * 2019-11-05 2020-01-31 中国电波传播研究所（中国电子科技集团公司第二十二研究所） real-time recognition method for non-cooperative targets in sea clutter measurement area of shore-based radar

Also Published As

Publication number Publication date
CN112147588A (en) 2020-12-29

Similar Documents

Publication Publication Date Title
CN102288944B (en) Super-resolution height measuring method based on topographic matching for digital array meter wave radar
CN110764059B (en) Method for transmitting and receiving vertical beam three-coordinate phased array radar
WO2022262303A1 (en) Millimeter wave vehicle-mounted mimo radar antenna array design method and apparatus
CN103885054A (en) Meter-wave radar low-elevation height measuring method based on distributed source reflection model
WO2019071917A1 (en) Satellite tracking method
CN111398914B (en) Radar tracking beam angle measurement method
CN113532437A (en) Ultra-wideband indoor target positioning method based on multipath utilization
CN112147588B (en) Rapid calculation method for asymmetric radar irradiation area
CN108089161B (en) Antenna array synthetic beam spatial domain selection method based on power feedback
CN106291556B (en) A kind of high squint SAR imaging algorithm based on hyperbolic model
CN111649803B (en) Three-dimensional radar level meter based on vertical linear array and design method thereof
CN110176665A (en) A kind of phased array antenna and phased-array radar
CN115113156B (en) Calibration method and system for dual-polarized phased array meteorological radar
CN110907925B (en) Weight positioning method under high-frequency ground wave radar double-station model
CN111381212A (en) Virtual ultrashort baseline positioning method based on subarray division
CN109901131A (en) Multipath based on oblique projection utilizes coherent beam formation method
CN110334480B (en) Design method of secondary surface extended curved surface of double-offset antenna for reducing noise temperature
CN114910944A (en) Phased array antenna tracking precision measuring method
CN103913731B (en) Based on the target elevation method of estimation of prior imformation
CN110532631B (en) 6G communication antenna array element position tolerance determination method based on channel capacity sensitivity
Yang et al. A single-layer reflectarray antenna for sub-thz applications
CN103913726B (en) A kind of radar frame height scaling method of based target echoed signal
CN217544916U (en) Deformed residual cutting shaped millimeter wave antenna
CN113904127B (en) Ultra-wideband high-gain direction-finding antenna based on side lobe suppression antenna feed source
CN109001704B (en) Target echo signal calculation method for sea surface low-angle radar

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