# Research on Quantization Error Influence of Millimeter-Wave Phased Array Antenna

Authors: Haigen Yang; Linqun Zhu; Zhun Xia; Yanqing Chen; Luohao Dai; Ruotian Xu; YuanHao Chen; Luyang Li; GuiYing Sun; Hongyang Yu; Wenting Xu
Journal: International Journal of Antennas and Propagation
 (2021)

---

## Abstract

The millimeter-wave phased array antenna is a higher integration system that is composed of different subarray modules, and in actual engineering, the existing amplitude, phase errors, and structural errors will change the performance of the array antenna. This paper studies the influence of the random amplitude and phase errors of the antenna array in the actual assembly process and the actual position errors between the subarrays on the electrical performance of the antenna. Based on the planar rectangular antenna array-electromagnetic coupling model, we propose a method of verifying the effect of random errors on the phased array antenna. The simulation result shows that the method could obtain the critical value of the error generated by the antenna subarray during processing and assembly. To reduce the error factor, it is necessary to ensure that the random phase and amplitude error should not exceed10°,0.5dB. The error in the X-direction during assembly should be ≤0.05λ, and the error in the Y-direction should be ≤0.1λ. When symmetrical deformation occurs, the maximum deformation should be less than 0.05λ.


---

## Body

## 1. Introduction

The antenna is an important part of radar, and the advancement of radar technology is inseparable from reliable and stable radar signals, which will directly affect the detection effect of radar [1]. However, the directivity of a single antenna is limited, to use antennas for electrical scanning in space, several antennas can be arranged together regularly to produce a directional pattern, which is called an antenna array [2]. The millimeter wave has a short wavelength and has the comprehensive advantages of microwave and light waves. The antenna array using millimeter wave has the advantages of extremely wide bandwidth, small size, compact system structure, and electromagnetic energy focusing [3], which is especially suitable for radar and other equipment.Driven by the rapid development of mobile communication systems, in order to better apply this multi-antenna array element structure to radar equipment, phased array systems and multiple input multiple output (MIMO) systems have emerged. Among them, the phased array system has been widely used on ground radars as early as the 1850s, and the MIMO concept was first introduced into radar technology in 2004, mainly using beam diversity to analyze and study the angle of arrival of the signal. Although the theories are similar, the application scenarios are quite different. MIMO radar is mostly used in civilian applications, such as the Industrial Internet of Things (IoT), due to its high packet loss rate and decreased multiplexing gain in high-speed conditions, and it is used in conjunction with the 5G technology being deployed to meet the needs of industrial communication systems [4]. In order to improve its performance, the model and algorithm need to be continuously improved. The so-called quaternion noncircular MUSIC (QNC-MUSIC) algorithm was proposed to improve the accuracy of DOA estimation in [5]. Shi et al. constructed a generalized tensor model in [6] and optimized the tensor by maximizing the number of detection targets, and finally derived the Cramér–Rao Bound (CRB) of nested radars, which proves the superiority of the method. A bistatic co-prime EMVS-MIMO radar framework was proposed in [7], and the work in [8] proposed a new closed-form estimation algorithm for EMVS-MIMO radar, constructed a new rotation invariant characteristic, and achieved better estimation than existing algorithms.Compared with MIMO radar, phased array radar has a huge cost, larger volume and weight, a long history of development, and relatively complete technology, and it is mostly used in large-scale shipborne and airborne radars [9]. The antenna element spacing at the transmitting end of a phased array radar is usually on the order of wavelength. In order to prevent grating lobes from appearing, the spacing is usually set to half a wavelength. In MIMO radar, the sensor spacing should not exceed half a wavelength to avoid phase ambiguity [10]. Compared with conventional arrays, active phased array antennas have the advantages of multifunction, high reliability, and high detection and tracking capabilities [11]. Vollbracht determined the optimal phase excitation distribution through the study of the single-feed antenna array, and it was confirmed on the 4×8 antenna subarray [12]. Sharma also proposed a new phased array composed of antijamming antennas [13]. Recently, a highly efficient polarized 8×8 millimeter-wave antenna array in the 60GHz frequency band was proposed in [14]. Compared with the traditional array, its impedance bandwidth, stable gain, and radiation within the bandwidth are improved. Ortiz et al. used a mathematical model based on diffraction theory to evaluate the influence of mutual coupling between polarized antennas [15].Conventional phased array antennas have problems such as greater complexity and high cost [16]; we can use digital technologies splice digital array unit formed of different sizes and form fronts to reduce design difficulty [17]. However, in practical engineering, there are extensive excitation errors caused by amplitude and phase changes [18], as well as deformation errors caused by processing and using. Therefore, in order to achieve the expected antenna sidelobe requirements and obtain a stable array design, it is necessary to perform error analysis on the assembled antenna array [19]. Wang analyzed the influence of random feeding errors such as the failure rate of the array element and the feeding amplitude and phase error on the electrical performance of the phased array antenna, which analyzed the bowl surface deformation and bending deformation of both systematic errors impact the electrical properties of a phased array antenna [20]. Chen and Zhou proposed a mathematical model of the thermal deformation error of the active phased array to analyze the influence of such errors on the antenna pattern [21].The current research on the performance of phased array radar mainly focuses on the influence of various factors on the electrical performance of the antenna when the radar is used, and the errors caused by the production process and assembly accuracy are often ignored. This article is mainly based on the error theory modeling of the millimeter-wave array antenna. The model is composed of two16×8 rectangular subarrays and works in the 24.25 GHz∼27.5 GHz frequency band. Based on the array antenna algorithm of this model, we can introduce random amplitude and phase errors of each channel and structural errors caused by array deformation [22]. Observe its influence on the performance of antenna gain, sidelobe level, beam width, etc. [23–25], realize the rapid analysis of the influence of the two types of errors on the electrical performance of the antenna array, establish the interval model of the relevant input parameters, and give the fluctuation interval of the relevant index, to provide support for the robust design of the millimeter-wave phased array antenna array.


## 2. Research of Active Phased Array Antennas

As a special type of antenna, the phased array antenna is composed of many identical independent antenna elements to form an antenna array. By controlling the radiation energy and phase relationship between each unit and using different array arrangements for each antenna unit in the array, the excitation feedback relationship affects the radiation field of the entire array, and accurate and predictable radiation patterns and beam directions are obtained, to improve the gain, scanning frequency, and anti-interference of the array [26].Compared with the traditional passive array, the active array adopts a distributed feedback structure, that is, each unit in the antenna array has a complete transmit/receive (T/R) unit to achieve high-power amplification and high-sensitivity reception [27]. Compared with other arrays, the active phased array has flexibility in beam direction, can form and independently control multiple beams, and search, identify, and track multiple targets. At the same time, because the active phased array antenna adopts millimeter wave with high power, its anti-interference ability is greatly increased [28].Phased array antennas are usually divided into two types: linear array and area array. The linear array only has one-dimensional scanning capability. If you want to have two-dimensional scanning capability, you need to combine multiple one-dimensional linear arrays to form a planar array antenna.Figure1 shows a rectangle of M×N antenna elements of the grid planar array antenna; the center position of the antenna element is rm,n=x^mdx+y^ndy. Assuming that the antenna element patterns are the same, and the analysis method is similar to a linear array, the array factor of the planar array antenna can be obtained as(1)Fθ,ϕ=∑m,nam,nexpjkmdxu−u0+ndyv−v0,where u=sinθcosϕ, v=sinθsinϕ, u0=sinθ0cosϕ0, and v0=sinθ0sinϕ0. The pencil beam is a typical feature of the planar antenna array pattern. When the excitation amplitude of each antenna element is equal, a uniform planar array is obtained. Normalize the above formula to get the matrix factor of the uniform plane array:(2)Fθ,ϕ=sinMπdxu−u0/λMsinπdxu−u0/λ+sinNπdyv−v0/λNsinπdyv−v0/λ.Figure 1
Planar array antenna diagram [29].


## 3. Error Modeling Theory

In the actual design process, the influence of various errors on the active phased array antenna must be considered. These errors have both random errors and systematic errors, which may be caused by component defects, or they may be caused by the feeder network, or other factors. These errors may cause the antenna gain to decrease, the sidelobe level to increase, the beam width to expand, the radiation efficiency to decrease, and even affect the antenna’s beam pointing, which largely determines the performance of active phased array radar [30]. The errors mentioned above are all random errors, which cannot be eliminated due to their contingency and unpredictability, and the performance of all large antenna arrays is affected by random errors such as element position and amplitude and phase changes [31].


### 3.1. Planar Array Antenna Modeling

In order to facilitate the simulation calculation of the two-dimensional phased array antenna pattern, using the coordinate system is usually chosen to express the direction in the simulation. In order to facilitate multidimensional process analysis and experimental verification, the representation method adopted is as follows: the two-dimensional pattern is represented by the antenna coordinate system, and the three-dimensional pattern is represented by the sinusoidal coordinate system [32].This article revolves around two-dimensional array pattern modeling. In the antenna coordinate system, the coordinate ofR can be characterized by the parameter θ and the parameter ϕ. Among them, θ is the angle from the z-axis to the point R00≤θ≤1800, ϕ, and ϕ is the angle between the x-axis and the projection line of R on the XY-plane. Normalize the amplitude of R, and the coordinate value corresponding to R0 is sinθcosϕ,sinθsinϕ,cosθ. The antenna coordinate system of the array direction diagram is shown in Figure 2.Figure 2
Antenna coordinate system [29].The sine space is a hemispherical mapping from a three-dimensional space to a two-dimensional plane, represented by three variablesU, V, and W, which can more intuitively characterize the antenna coordinate system. The transformation formula from antenna coordinate system to sine space is as follows:(3)U=sinθcosϕ,V=sinθsinϕ,W=cosθ.For the visible space of a two-dimensional planar phased array antenna, the value range ofU and V is −1,1, and the value range of W is 0,1. For an ideal planar array, there is no z-component in the position function of the array element, so there is no need to consider the influence of W on the array.Assuming the two-dimensional planar phased array hasM×N radiation channels [33], the rectangular grid arrangement, the channel spacing is dx and dy, respectively, the phase center position of the radiation unit is rm,n=x^mdx+y^ndy. The array arrangement can be seen in Figure 1.Then, the array factor pattern of the planar array antenna is(4)AF=∑m,nM,NAmnexpjkmdxu+ndyv,k=2πλ.Among them,M and N, respectively, represent the number of azimuth and elevation dimensions of the planar phased array antenna, and u and v, respectively, represent the coordinates of the UV-plane in the sine space representation, where the corresponding point of the incident angle is located, u=sinθcosϕ, v=sinθsinϕ. θ and ϕ are two parameters representing the spatial angle in the antenna coordinate system; Amn is the complex excitation signal of radiation channels m,n:(5)Amn=Imnexp−jkmdxu0+ndyv0,u0=sinθ0cosϕ0,v0=sinθ0sinϕ0.Imn is the weighted amplitude of radiation channels m,n, and θ0,φ0 is the beam direction of the antenna. Then, the pattern of the two-dimensional array antenna can be expressed as(6)AF=∑m,nImnexpjkmdxu−u0+ndyv−v0.According to the principle of pattern product, considering the influence of the element pattern on the array pattern, the expression of the two-dimensional array antenna can be expressed as(7)Eθ,ϕ=EP∗AF=fθ,ϕ∗AF,where EP=fθ,ϕ is the pattern of the antenna elements in the array.


### 3.2. Random Error of Amplitude and Phase

In actual engineering, there will be some errors in the amplitude and phase of the array element. Assume that the amplitude error of array elementm,n is ΔImn and the phase error is Δϕmn, where Δϕmn is a smaller amount after being converted to the radian system. Then, the array factor of the array antenna becomes(8)AF=∑m,nImn1+Δδmnexpjkmdxu−u0+ndyv−v0expjΔϕmn.According to Taylor’s expansion,(9)expjΔϕmn=1+jΔϕmn+1+j∗oΔϕmn.After finishing the expression of the antenna array factor, we can get(10)AF=∑m,nImn1+Δδmnexpjkmdxu−u0+ndyv−v0∗1+jΔϕmn,=∑m,nImnexpjkmdxu−u0+ndyv−v0+∑m,nΔδmnexpjkmdxu−u0+ndyv−v0+jΔϕmn∑m,nImnexpjkmdxu−u0+ndyv−v0+jΔϕmn∑m,nΔδmnexpjkmdxu−u0+ndyv−v0=AF0+ΔAF1+ΔAF2+ΔAF3.Among them, the elements are(11)ΔAF1=∑m,nΔδmnAF0,ΔAF2=jΔϕmnAF0,ΔAF3=jΔϕmn∑m,nΔδmnexpjkmdxu−u0+ndyv−v0.SinceΔAF3 is the product of two minimal errors, it is ignored here. The array factor can be simplified as(12)AF=AF0+ΔAF1+ΔAF2,where AF0 is the array factor of the two-dimensional array antenna under the ideal amplitude and phase distribution; ΔAF1 is the amount of change caused by the amplitude error; and ΔAF2 is the amount of change caused by the phase error.(13)Eθ,ϕ=EP∗AF=fθ,ϕ∗AF.The ideal power lobe function is(14)P0=E0u∗E0∗u.Theoretically, the sidelobes ofE0u and P0u can be designed to be arbitrarily low, but due to the existence of amplitude and phase errors, the reduction of the sidelobe level is limited. When there is amplitude and phase error between each unit, column and column, the lobe function is(15)Eθ,ϕ=fθ,ϕ∑m,nImn1+δmnexpjΔφmn∗expjkmdxu−u0+ndyv−v0,where δmn and φmn are the amplitude and phase distributions of radiating element m,n in the antenna array. The amplitude and phase error can be expressed by Gaussian distribution, the mean value is 0, and the variance is δmn2, φmn2. According to the central limit theorem, it can be proved that the sidelobe level R after considering the error obeys the Rician distribution; namely,(16)PR=Rδ2=exp−R2+Sm2δ2∗I0RSmδ2,I0x=12π∫−ππexcosφdφ,δE2=δae2+δφe2MNη,η=∑∑Imn2MN∑∑Imn2.I0 is the zero-order modified Bessel function; Sm is the sidelobe level value under ideal conditions; δ2=δE2/2, δE2 is the variance of the lobe, which represents the degree of agreement between the actual lobe and the theoretical design band; and η characterizes the aperture efficiency of the array weight. Therefore, the probability that the sidelobe is lower than the given value RT is(17)PR<RT=∫0RTPRdR,limit sidelobe level≈11dB+101gδE22.In practical applications, in order to facilitate the implementation of the project, according to the modular design requirements of the antenna array, the array is usually divided into multiple small subarrays or modules, which are cascaded through the feed network. In the cascading process, random errors are generated in units of subarrays. The following formula is the expression of the antenna pattern when the array is used as a subarray to construct the array. The division of other subarrays is similar:(18)Eθ,ϕ=fθ,ϕ∑n1+δnexpjΔφn∗∑mImn1+δmnexpjΔφmnexpjkmdxu−u0+ndyv−v0,where δn and φn are the amplitude and phase distribution of the nth column of radiating elements in the antenna array.


### 3.3. Structural Error

The structural error of the antenna array includes the processing error of the array, the installation error of the radiating element, the antenna frame and the subarray, and the structural deformation error caused by the deformation of the antenna. In the active phased array, processing, assembly, and other links will cause the deformation of the array and generate random errors. In the actual working environment, the factors such as vibration, impact, and high and low temperature will also cause the deformation of the planar array and finally change the position of the element to produce deformation errors, which will reduce the electromagnetic performance of the antenna [34]. As a result, the active phased array has raised sidelobes, decreased gain, and worsened pointing accuracy. Therefore, the study of the relationship between the coupling active phased array structure and electromagnetic analysis antenna electrical properties varies with structural changes in the error active phased bursts of the same size curve surface [35, 36], to obtain the critical value wavefront deformation and random error combined, can provide quantitative theoretical guidance for structural design and reasonable allocation of tolerances.


## 3.1. Planar Array Antenna Modeling

In order to facilitate the simulation calculation of the two-dimensional phased array antenna pattern, using the coordinate system is usually chosen to express the direction in the simulation. In order to facilitate multidimensional process analysis and experimental verification, the representation method adopted is as follows: the two-dimensional pattern is represented by the antenna coordinate system, and the three-dimensional pattern is represented by the sinusoidal coordinate system [32].This article revolves around two-dimensional array pattern modeling. In the antenna coordinate system, the coordinate ofR can be characterized by the parameter θ and the parameter ϕ. Among them, θ is the angle from the z-axis to the point R00≤θ≤1800, ϕ, and ϕ is the angle between the x-axis and the projection line of R on the XY-plane. Normalize the amplitude of R, and the coordinate value corresponding to R0 is sinθcosϕ,sinθsinϕ,cosθ. The antenna coordinate system of the array direction diagram is shown in Figure 2.Figure 2
Antenna coordinate system [29].The sine space is a hemispherical mapping from a three-dimensional space to a two-dimensional plane, represented by three variablesU, V, and W, which can more intuitively characterize the antenna coordinate system. The transformation formula from antenna coordinate system to sine space is as follows:(3)U=sinθcosϕ,V=sinθsinϕ,W=cosθ.For the visible space of a two-dimensional planar phased array antenna, the value range ofU and V is −1,1, and the value range of W is 0,1. For an ideal planar array, there is no z-component in the position function of the array element, so there is no need to consider the influence of W on the array.Assuming the two-dimensional planar phased array hasM×N radiation channels [33], the rectangular grid arrangement, the channel spacing is dx and dy, respectively, the phase center position of the radiation unit is rm,n=x^mdx+y^ndy. The array arrangement can be seen in Figure 1.Then, the array factor pattern of the planar array antenna is(4)AF=∑m,nM,NAmnexpjkmdxu+ndyv,k=2πλ.Among them,M and N, respectively, represent the number of azimuth and elevation dimensions of the planar phased array antenna, and u and v, respectively, represent the coordinates of the UV-plane in the sine space representation, where the corresponding point of the incident angle is located, u=sinθcosϕ, v=sinθsinϕ. θ and ϕ are two parameters representing the spatial angle in the antenna coordinate system; Amn is the complex excitation signal of radiation channels m,n:(5)Amn=Imnexp−jkmdxu0+ndyv0,u0=sinθ0cosϕ0,v0=sinθ0sinϕ0.Imn is the weighted amplitude of radiation channels m,n, and θ0,φ0 is the beam direction of the antenna. Then, the pattern of the two-dimensional array antenna can be expressed as(6)AF=∑m,nImnexpjkmdxu−u0+ndyv−v0.According to the principle of pattern product, considering the influence of the element pattern on the array pattern, the expression of the two-dimensional array antenna can be expressed as(7)Eθ,ϕ=EP∗AF=fθ,ϕ∗AF,where EP=fθ,ϕ is the pattern of the antenna elements in the array.


## 3.2. Random Error of Amplitude and Phase

In actual engineering, there will be some errors in the amplitude and phase of the array element. Assume that the amplitude error of array elementm,n is ΔImn and the phase error is Δϕmn, where Δϕmn is a smaller amount after being converted to the radian system. Then, the array factor of the array antenna becomes(8)AF=∑m,nImn1+Δδmnexpjkmdxu−u0+ndyv−v0expjΔϕmn.According to Taylor’s expansion,(9)expjΔϕmn=1+jΔϕmn+1+j∗oΔϕmn.After finishing the expression of the antenna array factor, we can get(10)AF=∑m,nImn1+Δδmnexpjkmdxu−u0+ndyv−v0∗1+jΔϕmn,=∑m,nImnexpjkmdxu−u0+ndyv−v0+∑m,nΔδmnexpjkmdxu−u0+ndyv−v0+jΔϕmn∑m,nImnexpjkmdxu−u0+ndyv−v0+jΔϕmn∑m,nΔδmnexpjkmdxu−u0+ndyv−v0=AF0+ΔAF1+ΔAF2+ΔAF3.Among them, the elements are(11)ΔAF1=∑m,nΔδmnAF0,ΔAF2=jΔϕmnAF0,ΔAF3=jΔϕmn∑m,nΔδmnexpjkmdxu−u0+ndyv−v0.SinceΔAF3 is the product of two minimal errors, it is ignored here. The array factor can be simplified as(12)AF=AF0+ΔAF1+ΔAF2,where AF0 is the array factor of the two-dimensional array antenna under the ideal amplitude and phase distribution; ΔAF1 is the amount of change caused by the amplitude error; and ΔAF2 is the amount of change caused by the phase error.(13)Eθ,ϕ=EP∗AF=fθ,ϕ∗AF.The ideal power lobe function is(14)P0=E0u∗E0∗u.Theoretically, the sidelobes ofE0u and P0u can be designed to be arbitrarily low, but due to the existence of amplitude and phase errors, the reduction of the sidelobe level is limited. When there is amplitude and phase error between each unit, column and column, the lobe function is(15)Eθ,ϕ=fθ,ϕ∑m,nImn1+δmnexpjΔφmn∗expjkmdxu−u0+ndyv−v0,where δmn and φmn are the amplitude and phase distributions of radiating element m,n in the antenna array. The amplitude and phase error can be expressed by Gaussian distribution, the mean value is 0, and the variance is δmn2, φmn2. According to the central limit theorem, it can be proved that the sidelobe level R after considering the error obeys the Rician distribution; namely,(16)PR=Rδ2=exp−R2+Sm2δ2∗I0RSmδ2,I0x=12π∫−ππexcosφdφ,δE2=δae2+δφe2MNη,η=∑∑Imn2MN∑∑Imn2.I0 is the zero-order modified Bessel function; Sm is the sidelobe level value under ideal conditions; δ2=δE2/2, δE2 is the variance of the lobe, which represents the degree of agreement between the actual lobe and the theoretical design band; and η characterizes the aperture efficiency of the array weight. Therefore, the probability that the sidelobe is lower than the given value RT is(17)PR<RT=∫0RTPRdR,limit sidelobe level≈11dB+101gδE22.In practical applications, in order to facilitate the implementation of the project, according to the modular design requirements of the antenna array, the array is usually divided into multiple small subarrays or modules, which are cascaded through the feed network. In the cascading process, random errors are generated in units of subarrays. The following formula is the expression of the antenna pattern when the array is used as a subarray to construct the array. The division of other subarrays is similar:(18)Eθ,ϕ=fθ,ϕ∑n1+δnexpjΔφn∗∑mImn1+δmnexpjΔφmnexpjkmdxu−u0+ndyv−v0,where δn and φn are the amplitude and phase distribution of the nth column of radiating elements in the antenna array.


## 3.3. Structural Error

The structural error of the antenna array includes the processing error of the array, the installation error of the radiating element, the antenna frame and the subarray, and the structural deformation error caused by the deformation of the antenna. In the active phased array, processing, assembly, and other links will cause the deformation of the array and generate random errors. In the actual working environment, the factors such as vibration, impact, and high and low temperature will also cause the deformation of the planar array and finally change the position of the element to produce deformation errors, which will reduce the electromagnetic performance of the antenna [34]. As a result, the active phased array has raised sidelobes, decreased gain, and worsened pointing accuracy. Therefore, the study of the relationship between the coupling active phased array structure and electromagnetic analysis antenna electrical properties varies with structural changes in the error active phased bursts of the same size curve surface [35, 36], to obtain the critical value wavefront deformation and random error combined, can provide quantitative theoretical guidance for structural design and reasonable allocation of tolerances.


## 4. Simulation Analysis of Random Error in Amplitude and Phase

### 4.1. Random Error Simulation Analysis

The amplitude and phase errors in the simulation analysis process of amplitude and phase random errors are random values. The random error probability distribution obeys the Gaussian distribution, and the error value range is defined by the variance.In the simulation process, taking a16×16 planar phased array antenna as an example, the element spacing dx=0.43λ, dy=0.52λ, and the antenna element pattern uses Gaussian beams. In the process of random error analysis, the error source includes quantization error. The minimum quantization step size of the numerically controlled attenuator is 0.5dB, and the digital phase shifter adopts a 6-phase shifter.Let us take the random phase and amplitude errors with variances of0°,0dB\0.3dB\0.5dB\1dB\2dB and 10°,0dB\0.3dB\0.5dB\1dB\2dB as an example to simulate the array pattern under the influence of random error. The simulation results of random amplitude and phase antenna without scanning and scanning −45° are shown in Figures 3 and 4. In the same way, the simulation result data under the conditions of 5°,0dB\0.3dB\0.5dB\1dB\2dB and 20°,0dB\0.3dB\0.5dB\1dB\2dB are plotted and analyzed with the above experimental data.Figure 3
Simulation results of random amplitude and phase antenna without scanning pattern: (a) phase is0°; (b) phase is 10°.
(a)(b)Figure 4
Random amplitude and phase antenna scanning−45° pattern simulation results: (a) phase is 0°; (b) phase is 10°.
(a)(b)Table1 shows the statistical results of technical indicators under the abovementioned random error conditions. Random errors have a significant impact on antenna gain, sidelobe level, and beam width and have little impact on beam pointing. As the error value increases, the antenna gain gradually decreases, the sidelobe level gradually increases, and the beam width gradually expands. When the error value is large, the beam direction will also change to a certain extent.Table 1
Technical index statistics under the condition of random amplitude and phase error.
Random phase error (variance (°))Scan angleRandom amplitude error (variance (dB))GainSidelobe level3 dB beam widthBeam pointing00028.22−20.258.000.000.328.24−19.618.000.000.528.17−20.778.000.00128.05−20.228.000.00227.51−17.928.200.00−45026.08−18.9511.10−44.400.326.02−18.5611.20−44.400.525.94−18.1711.30−44.40125.66−17.5611.40−44.40225.66−17.5611.40−44.4050028.20−20.088.000.000.328.18−20.028.000.000.528.17−19.828.000.00128.02−18.868.200.00227.70−19.368.100.00−45026.07−18.5511.20−44.400.326.04−19.3811.30−44.400.526.00−18.8211.30−44.40125.98−17.1711.10−44.40225.62−17.0311.00−44.50100028.14−20.048.00−0.100.328.13−19.618.000.000.528.08−19.658.00−0.10127.96−18.148.100.00227.67−19.187.900.00−45026.00−18.9111.20−44.500.325.99−19.9211.20−44.400.525.98−18.1211.20−44.40125.91−16.2711.10−44.50225.30−17.1611.50−44.40200027.84−17.838.000.100.327.90−20.038.00−0.100.527.84−17.778.100.10127.76−16.428.000.00227.27−19.628.000.10−45025.80−19.0711.20−44.500.325.81−18.8711.20−44.400.525.72−18.8411.30−44.40125.65−17.1911.10−44.30225.27−19.2211.30−44.30According to the simulation results in Table1, under the influence of random errors, the gain, sidelobe level, etc., change, which affects the array antenna gain and overall performance. As the random error variance increases, the average sidelobe level and the first sidelobe level of the array pattern also increase. When the amplitude random error variance exceeds 0.5dB and the phase random error variance exceeds 10°, not only will the energy of the array pattern become more dispersed, the randomness of the sidelobe positions will increase, and the level of the first sidelobe will increase significantly compared to ideal conditions. In addition, the effect of random error on antenna beam scanning is like that when it is not scanning.


### 4.2. Verification of Multiple Random Results

The following is to verify the conclusion that the amplitude and phase random error is not greater than0.5dB,10°. The simulation result of the antenna pattern with five random amplitude and phase errors is shown in Figure 5.Figure 5
Simulation results of antenna pattern with five random amplitude and phase errors. (a) Azimuth plane. (b) Elevation plane.
(a)(b)Table2 counts the main technical indicators of antenna azimuth and elevation antennas with 5 different random errors. After comparison, the gain drop is less than or equal to 0.2dB, the sidelobe level is less than or equal to 1dB on average, the beam width changes up to 0.2°, and the influence of beam pointing is small.Table 2
Five-order random amplitude and phase error antenna technical index statistics.
No.Directional map sectionRandom timesGainSidelobe level3 dB beam widthBeam pointing1Azimuth planeIdeal situation28.22−20.348.000.002128.07−18.767.800.003228.13−19.407.800.004328.11−19.848.000.005428.12−19.737.800.006528.11−18.898.000.007Elevation planeIdeal situation28.22−20.296.400.008128.09−20.566.600.009228.15−19.256.400.0010328.14−19.446.600.0011428.12−19.976.600.0012528.11−18.886.600.00


## 4.1. Random Error Simulation Analysis

The amplitude and phase errors in the simulation analysis process of amplitude and phase random errors are random values. The random error probability distribution obeys the Gaussian distribution, and the error value range is defined by the variance.In the simulation process, taking a16×16 planar phased array antenna as an example, the element spacing dx=0.43λ, dy=0.52λ, and the antenna element pattern uses Gaussian beams. In the process of random error analysis, the error source includes quantization error. The minimum quantization step size of the numerically controlled attenuator is 0.5dB, and the digital phase shifter adopts a 6-phase shifter.Let us take the random phase and amplitude errors with variances of0°,0dB\0.3dB\0.5dB\1dB\2dB and 10°,0dB\0.3dB\0.5dB\1dB\2dB as an example to simulate the array pattern under the influence of random error. The simulation results of random amplitude and phase antenna without scanning and scanning −45° are shown in Figures 3 and 4. In the same way, the simulation result data under the conditions of 5°,0dB\0.3dB\0.5dB\1dB\2dB and 20°,0dB\0.3dB\0.5dB\1dB\2dB are plotted and analyzed with the above experimental data.Figure 3
Simulation results of random amplitude and phase antenna without scanning pattern: (a) phase is0°; (b) phase is 10°.
(a)(b)Figure 4
Random amplitude and phase antenna scanning−45° pattern simulation results: (a) phase is 0°; (b) phase is 10°.
(a)(b)Table1 shows the statistical results of technical indicators under the abovementioned random error conditions. Random errors have a significant impact on antenna gain, sidelobe level, and beam width and have little impact on beam pointing. As the error value increases, the antenna gain gradually decreases, the sidelobe level gradually increases, and the beam width gradually expands. When the error value is large, the beam direction will also change to a certain extent.Table 1
Technical index statistics under the condition of random amplitude and phase error.
Random phase error (variance (°))Scan angleRandom amplitude error (variance (dB))GainSidelobe level3 dB beam widthBeam pointing00028.22−20.258.000.000.328.24−19.618.000.000.528.17−20.778.000.00128.05−20.228.000.00227.51−17.928.200.00−45026.08−18.9511.10−44.400.326.02−18.5611.20−44.400.525.94−18.1711.30−44.40125.66−17.5611.40−44.40225.66−17.5611.40−44.4050028.20−20.088.000.000.328.18−20.028.000.000.528.17−19.828.000.00128.02−18.868.200.00227.70−19.368.100.00−45026.07−18.5511.20−44.400.326.04−19.3811.30−44.400.526.00−18.8211.30−44.40125.98−17.1711.10−44.40225.62−17.0311.00−44.50100028.14−20.048.00−0.100.328.13−19.618.000.000.528.08−19.658.00−0.10127.96−18.148.100.00227.67−19.187.900.00−45026.00−18.9111.20−44.500.325.99−19.9211.20−44.400.525.98−18.1211.20−44.40125.91−16.2711.10−44.50225.30−17.1611.50−44.40200027.84−17.838.000.100.327.90−20.038.00−0.100.527.84−17.778.100.10127.76−16.428.000.00227.27−19.628.000.10−45025.80−19.0711.20−44.500.325.81−18.8711.20−44.400.525.72−18.8411.30−44.40125.65−17.1911.10−44.30225.27−19.2211.30−44.30According to the simulation results in Table1, under the influence of random errors, the gain, sidelobe level, etc., change, which affects the array antenna gain and overall performance. As the random error variance increases, the average sidelobe level and the first sidelobe level of the array pattern also increase. When the amplitude random error variance exceeds 0.5dB and the phase random error variance exceeds 10°, not only will the energy of the array pattern become more dispersed, the randomness of the sidelobe positions will increase, and the level of the first sidelobe will increase significantly compared to ideal conditions. In addition, the effect of random error on antenna beam scanning is like that when it is not scanning.


## 4.2. Verification of Multiple Random Results

The following is to verify the conclusion that the amplitude and phase random error is not greater than0.5dB,10°. The simulation result of the antenna pattern with five random amplitude and phase errors is shown in Figure 5.Figure 5
Simulation results of antenna pattern with five random amplitude and phase errors. (a) Azimuth plane. (b) Elevation plane.
(a)(b)Table2 counts the main technical indicators of antenna azimuth and elevation antennas with 5 different random errors. After comparison, the gain drop is less than or equal to 0.2dB, the sidelobe level is less than or equal to 1dB on average, the beam width changes up to 0.2°, and the influence of beam pointing is small.Table 2
Five-order random amplitude and phase error antenna technical index statistics.
No.Directional map sectionRandom timesGainSidelobe level3 dB beam widthBeam pointing1Azimuth planeIdeal situation28.22−20.348.000.002128.07−18.767.800.003228.13−19.407.800.004328.11−19.848.000.005428.12−19.737.800.006528.11−18.898.000.007Elevation planeIdeal situation28.22−20.296.400.008128.09−20.566.600.009228.15−19.256.400.0010328.14−19.446.600.0011428.12−19.976.600.0012528.11−18.886.600.00


## 5. Simulation Analysis of Structural Error

Figure6 is a schematic diagram of a millimeter-wave antenna array and the center number of the array unit, where the x-axis is the azimuth direction, and the y-axis is the distance direction. The related data of the array is as follows:(a)
Working frequency: 24.25 GHz–27.5 GHz(b)
Model boundary size: 100 mm (azimuth direction)∗ 112 mm (range direction)(c)
Rectangular array unit arrangement: 16 (azimuth direction)∗ 16 (azimuth direction)(d)
Antenna unit center spacing: 5 mm (azimuth direction)∗ 6 mm (azimuth direction)(e)
Substrate material: Rogers 4350B. The thickness is 2.5 mm.Figure 6
Schematic diagram of a millimeter-wave antenna array.


### 5.1. Subarray Assembly Error

In the actual antenna design,16×8=128 units are used as a subarray, and errors in the position of the subarray must be considered during the assembly process.This paper takes16×16 planar phased array antenna, that is, two subarrays, as an example. The element spacing is dx=0.43λ，dy=0.52λ, the antenna element pattern uses Gaussian beam, and the element pattern function is EP=cosθEF/2. Among them, EF determines the gain of the radiating element. Here, 1.2 is selected, and the gain is about 6.4 dB. The position error of the left and right half arrays is simulated by MATLAB programming simulation, and the influence of the subarray error in the XY-direction on the antenna performance is analyzed.


#### 5.1.1. Subarray Error in theX-Direction

When there is an0.2λ error in the X-direction of the left and right half of the array, the array element grid is shown in Figure 7. The experiment simulated the changes of the antenna pattern with an error value of 0:0.05λ:0.2λ.Figure 7
Position error of the left and right half arrays in theX-direction.Figure8 shows the projection of the three-dimensional pattern in the UV space when there is a 0.2λ error in the left and right half of the antenna. Figure 9 shows the simulation results of the directional pattern of the antenna azimuth plane when the left and right half arrays have different errors in the X-direction. The elevation plane has almost no effect. Figure 10 shows the changes in the technical indicators of the antenna azimuth plane when the left and right half arrays have different errors in the X-direction.Figure 8
Change of the position error pattern of the left and right half arrays in theX-direction 0.2λ.Figure 9
The change of the azimuth pattern caused by the position error of the left and right half arrays in theX-direction.Figure 10
Changes in technical indicators caused by the position error of the left and right half arrays in theX-direction. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)


#### 5.1.2. Subarray Error in theY-Direction

When there is a0.2λ error in the Y-direction of the left and right half of the array, the array element grid is shown in Figure 11. The experiment simulated the changes of the antenna pattern with an error value of 0:0.05λ:0.2λ.Figure 11
Position error of the left and right half arrays in theY-direction.Figure12 shows the projection of the three-dimensional pattern in the UV space when there is a 0.2λ error in the left and right half of the antenna. Figure 13 shows the simulation results of the antenna elevation plane pattern when the left and right half arrays have different errors in the Y-direction. The azimuth plane has little influence. Figure 14 shows the changes in the technical indicators of the antenna azimuth plane when the left and right half arrays have different errors in the Y-direction.Figure 12
Change of the position error pattern of the left and right half arrays in theY-direction 0.2λ.Figure 13
Change of the azimuth pattern caused by the position error of the left and right half arrays in theY-direction.Figure 14
Changes in technical indicators caused by the position error of the left and right half arrays in theY-direction. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)


#### 5.1.3. Subarray Error inXY-Direction

Figure15 shows the projection of the three-dimensional pattern in the UV space when there is a 0.1λ error in the left and right half of the antenna.Figure 15
Position error pattern changes of the left and right half arrays in theXY-direction 0.1λ.Table3 shows the statistics of technical indicators such as antenna gain and sidelobe level when the XY-directions change at the same time.Table 3
XY-direction changes at the same time technical indicator statistics.
No.Directional map sectionX-axisY-axisGainSidelobe level3 dB beamwidthBeam pointing1Azimuth plane0028.22−20.338020.0528.22−20.338030.128.22−20.338040.05028.26−19.438050.0528.26−19.438060.128.26−19.438070.1028.30−18.627.8080.0528.30−18.627.8090.128.30−18.627.8010Elevation plane0028.22−20.296.60110.0528.22−20.306.60120.128.22−20.316.60130.05028.26−20.296.60140.0528.26−20.306.60150.128.26−20.316.60160.1028.30−20.296.60170.0528.30−20.306.60180.128.30−20.316.60From the above simulation results, it can be concluded that the gaps produced by the left and right half arrays in theX-direction will affect the distribution of antenna radiation power in space, and the impact on the sidelobe level is more obvious. As the gap increases, the electrical aperture of the planar array becomes larger, the gain becomes larger, the beam width becomes narrower, and the beam direction does not change. The gaps produced by the left and right half arrays in the Y direction will affect the distribution of antenna radiation power in space, causing tilt rotation, which has a significant impact on the far-area sidelobe level. As the gap increases, the electrical aperture of the planar array becomes larger, the gain becomes larger, the beam width becomes narrower, and the beam direction does not change. When the XY-direction changes at the same time, it is a combination of the above two situations.In practical applications, according to the principle that the sidelobe level deterioration is not greater than 1 dB, the error in theX-direction should be ≤0.05λ, and the error in the Y-direction should be ≤0.1λ.


### 5.2. Subarray Deformation Error

Assuming that a vibration load is applied during the operation of the radar, under the constraints of the four corners, the antenna array will undergo symmetrical and asymmetrical deformation. For the planar rectangular active phased array antenna, the deformation of its subarray is usually saddle-shaped deformation. The constraints are shown in Figure16. A typical saddle-shaped deformation is listed below, the effect of antenna deformation on the antenna pattern is simulated and analyzed, and the main technical indicators of the antenna are counted.Figure 16
Matrix constraint condition.


#### 5.2.1. Saddle-Shaped Deformation of Arrays

Figure17 shows the saddle-shaped deformation of the antenna array. The normalized deformation data of the center of each antenna element is shown in Table 4.Figure 17
Saddle deformation.Table 4
Normalized deformation data at the center of the antenna element.
Unit no.Normalized data10.055920.2130.38440.5650.72260.85570.9518191100.951110.855120.722130.56140.384150.21160.055917−0.065180.0901190.257200.423210.575220.701230.791240.837250.837260.791270.701280.575290.423300.257310.090132−0.06533−0.21734−0.0543350.113360.275370.422380.543390.629400.673410.673420.629430.543440.422450.275460.11347−0.054348−0.21749−0.38150−0.20951−0.0381520.124530.269540.388550.472560.515570.515580.472590.388600.269610.12462−0.038163−0.20964−0.38165−0.53866−0.35867−0.18268−0.0176690.128700.245710.328720.371730.371740.328750.245760.12877−0.017678−0.18279−0.35880−0.53881−0.67282−0.48583−0.30584−0.138850.0083860.126870.209880.251890.251900.209910.126920.008393−0.13894−0.30595−0.48596−0.67297−0.76998−0.57899−0.395100−0.225101−0.07811020.04021030.1231040.1661050.1661060.1231070.0402108−0.0781109−0.225110−0.395111−0.578112−0.769113−0.821114−0.627115−0.442116−0.271117−0.123118−0.00481190.07811200.1211210.1211220.0781123−0.0048124−0.123125−0.271126−0.442127−0.627128−0.821129−0.821130−0.627131−0.442132−0.271133−0.123134−0.00481350.07811360.1211370.1211380.0781139−0.0048140−0.123141−0.271142−0.442143−0.627144−0.821145−0.769146−0.578147−0.395148−0.225149−0.07811500.04021510.1231520.1661530.1661540.1231550.0402156−0.0781157−0.225158−0.395159−0.578160−0.769161−0.672162−0.485163−0.305164−0.1381650.00831660.1261670.2091680.2511690.2511700.2091710.1261720.0083173−0.138174−0.305175−0.485176−0.672177−0.538178−0.358179−0.182180−0.01761810.1281820.2451830.3281840.3711850.3711860.3281870.2451880.128189−0.0176190−0.182191−0.358192−0.538193−0.381194−0.209195−0.03811960.1241970.2691980.3881990.4722000.5152010.5152020.4722030.3882040.2692050.124206−0.0381207−0.209208−0.381209−0.217210−0.05432110.1132120.2752130.4222140.5432150.6292160.6732170.6732180.6292190.5432200.4222210.2752220.113223−0.0543224−0.217225−0.0652260.09012270.2572280.4232290.5752300.7012310.7912320.8372330.8372340.7912350.7012360.5752370.4232380.2572390.0901240−0.0652410.05592420.212430.3842440.562450.7222460.8552470.951248124912500.9512510.8552520.7222530.562540.3842550.212560.0559Based on the normalized deformation, the maximum deformation of the antenna array was set, and the azimuth and elevation patterns of the antenna were simulated when the maximum deformation became0:0.025λ:0.4λ. Figure 18 shows the changes in the pattern of the array azimuth plane under saddle-shaped deformation. Figure 19 shows the changes in the main technical indicators of the array azimuth under saddle-shaped deformation. Figure 20 shows the changes in the pattern of the array elevation plane under saddle-shaped deformation. Figure 21 shows the changes in the main technical indicators of the array elevation plane under saddle-shaped deformation. Table 5 shows the changes in the main technical indicators of the antenna under the saddle-shaped deformation of the array.Figure 18
Azimuth plane of saddle-shaped deformation.Figure 19
Changes in technical indicators of azimuth plane caused by saddle-shaped deformation. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)Figure 20
Elevation plane of saddle-shaped deformation.Figure 21
Changes in technical indicators of elevation plane caused by saddle-shaped deformation. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)Table 5
Statistics of saddle-shaped deformation technical indicators.
No.Directional map sectionMax deformationGainSidelobe level3 dB beam widthBeam pointing1Azimuth plane028.22−20.338.000.0020.02528.20−20.098.000.0030.0528.14−19.408.000.0040.07528.05−18.378.000.0050.127.92−17.128.200.0060.12527.75−15.738.200.0070.1527.55−20.908.200.0080.17527.31−20.758.400.0090.227.02−20.548.400.00100.22526.71−20.248.600.00110.2526.35−19.868.800.00120.27525.95−19.379.000.00130.325.51−18.789.200.00140.32525.04−18.099.600.00150.3524.52−17.3010.000.00160.37523.96−16.4310.800.00170.423.37−15.4711.800.0018Elevation plane028.22−20.296.600.00190.02528.20−20.066.600.00200.0528.14−19.396.600.00210.07528.05−18.406.600.00220.127.92−17.216.800.00230.12527.75−15.916.800.00240.1527.55−21.056.800.00250.17527.31−20.996.800.00260.227.02−20.897.000.00270.22526.71−20.717.000.00280.2526.35−20.447.200.00290.27525.95−20.077.400.00300.325.51−19.597.400.00310.32525.04−18.997.600.00320.3524.52−18.298.000.00330.37523.96−17.498.400.00340.423.37−16.618.800.00From Figures18–21 and Table 5, the saddle-shaped deformation has a great influence on the azimuth plane and the elevation plane pattern of the planar array. As the amount of deformation increases, the antenna gain gradually decreases, the sidelobe level gradually increases, and the beam width gradually expands, which has almost no effect on the antenna beam direction. When the deformation amount is greater than or equal to 0.15λ, the azimuth plane and the elevation plane pattern begin to defocus, the pattern is no longer focused, and the first zero point disappears.


#### 5.2.2. Analysis of Deformation Error Impact

In order to verify whether the saddle-shaped deformation will affect the scanning beam of the antenna, the azimuth plane and the elevation plane, respectively, scan the directional patterns of −30° under different maximum deformation variables, as shown in Figures22 and 23. It can be seen from the figure that the influence of saddle-shaped deformation on antenna beam scanning is basically the same as when the antenna is not scanning.Figure 22
Saddle-shaped deformation azimuth scanning−30° direction map.Figure 23
Saddle-shaped deformation elevation scanning−30° direction map.Integrating the technical indicators of the azimuth plane and the elevation plane, the maximum deformation error of the saddle-shaped deformation is≤0.05λ and the deterioration of the technical indicators is acceptable. In engineering practice, the maximum deformation error value should be controlled to be ≤0.05λ.


## 5.1. Subarray Assembly Error

In the actual antenna design,16×8=128 units are used as a subarray, and errors in the position of the subarray must be considered during the assembly process.This paper takes16×16 planar phased array antenna, that is, two subarrays, as an example. The element spacing is dx=0.43λ，dy=0.52λ, the antenna element pattern uses Gaussian beam, and the element pattern function is EP=cosθEF/2. Among them, EF determines the gain of the radiating element. Here, 1.2 is selected, and the gain is about 6.4 dB. The position error of the left and right half arrays is simulated by MATLAB programming simulation, and the influence of the subarray error in the XY-direction on the antenna performance is analyzed.


### 5.1.1. Subarray Error in theX-Direction

When there is an0.2λ error in the X-direction of the left and right half of the array, the array element grid is shown in Figure 7. The experiment simulated the changes of the antenna pattern with an error value of 0:0.05λ:0.2λ.Figure 7
Position error of the left and right half arrays in theX-direction.Figure8 shows the projection of the three-dimensional pattern in the UV space when there is a 0.2λ error in the left and right half of the antenna. Figure 9 shows the simulation results of the directional pattern of the antenna azimuth plane when the left and right half arrays have different errors in the X-direction. The elevation plane has almost no effect. Figure 10 shows the changes in the technical indicators of the antenna azimuth plane when the left and right half arrays have different errors in the X-direction.Figure 8
Change of the position error pattern of the left and right half arrays in theX-direction 0.2λ.Figure 9
The change of the azimuth pattern caused by the position error of the left and right half arrays in theX-direction.Figure 10
Changes in technical indicators caused by the position error of the left and right half arrays in theX-direction. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)


### 5.1.2. Subarray Error in theY-Direction

When there is a0.2λ error in the Y-direction of the left and right half of the array, the array element grid is shown in Figure 11. The experiment simulated the changes of the antenna pattern with an error value of 0:0.05λ:0.2λ.Figure 11
Position error of the left and right half arrays in theY-direction.Figure12 shows the projection of the three-dimensional pattern in the UV space when there is a 0.2λ error in the left and right half of the antenna. Figure 13 shows the simulation results of the antenna elevation plane pattern when the left and right half arrays have different errors in the Y-direction. The azimuth plane has little influence. Figure 14 shows the changes in the technical indicators of the antenna azimuth plane when the left and right half arrays have different errors in the Y-direction.Figure 12
Change of the position error pattern of the left and right half arrays in theY-direction 0.2λ.Figure 13
Change of the azimuth pattern caused by the position error of the left and right half arrays in theY-direction.Figure 14
Changes in technical indicators caused by the position error of the left and right half arrays in theY-direction. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)


### 5.1.3. Subarray Error inXY-Direction

Figure15 shows the projection of the three-dimensional pattern in the UV space when there is a 0.1λ error in the left and right half of the antenna.Figure 15
Position error pattern changes of the left and right half arrays in theXY-direction 0.1λ.Table3 shows the statistics of technical indicators such as antenna gain and sidelobe level when the XY-directions change at the same time.Table 3
XY-direction changes at the same time technical indicator statistics.
No.Directional map sectionX-axisY-axisGainSidelobe level3 dB beamwidthBeam pointing1Azimuth plane0028.22−20.338020.0528.22−20.338030.128.22−20.338040.05028.26−19.438050.0528.26−19.438060.128.26−19.438070.1028.30−18.627.8080.0528.30−18.627.8090.128.30−18.627.8010Elevation plane0028.22−20.296.60110.0528.22−20.306.60120.128.22−20.316.60130.05028.26−20.296.60140.0528.26−20.306.60150.128.26−20.316.60160.1028.30−20.296.60170.0528.30−20.306.60180.128.30−20.316.60From the above simulation results, it can be concluded that the gaps produced by the left and right half arrays in theX-direction will affect the distribution of antenna radiation power in space, and the impact on the sidelobe level is more obvious. As the gap increases, the electrical aperture of the planar array becomes larger, the gain becomes larger, the beam width becomes narrower, and the beam direction does not change. The gaps produced by the left and right half arrays in the Y direction will affect the distribution of antenna radiation power in space, causing tilt rotation, which has a significant impact on the far-area sidelobe level. As the gap increases, the electrical aperture of the planar array becomes larger, the gain becomes larger, the beam width becomes narrower, and the beam direction does not change. When the XY-direction changes at the same time, it is a combination of the above two situations.In practical applications, according to the principle that the sidelobe level deterioration is not greater than 1 dB, the error in theX-direction should be ≤0.05λ, and the error in the Y-direction should be ≤0.1λ.


## 5.1.1. Subarray Error in theX-Direction

When there is an0.2λ error in the X-direction of the left and right half of the array, the array element grid is shown in Figure 7. The experiment simulated the changes of the antenna pattern with an error value of 0:0.05λ:0.2λ.Figure 7
Position error of the left and right half arrays in theX-direction.Figure8 shows the projection of the three-dimensional pattern in the UV space when there is a 0.2λ error in the left and right half of the antenna. Figure 9 shows the simulation results of the directional pattern of the antenna azimuth plane when the left and right half arrays have different errors in the X-direction. The elevation plane has almost no effect. Figure 10 shows the changes in the technical indicators of the antenna azimuth plane when the left and right half arrays have different errors in the X-direction.Figure 8
Change of the position error pattern of the left and right half arrays in theX-direction 0.2λ.Figure 9
The change of the azimuth pattern caused by the position error of the left and right half arrays in theX-direction.Figure 10
Changes in technical indicators caused by the position error of the left and right half arrays in theX-direction. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)


## 5.1.2. Subarray Error in theY-Direction

When there is a0.2λ error in the Y-direction of the left and right half of the array, the array element grid is shown in Figure 11. The experiment simulated the changes of the antenna pattern with an error value of 0:0.05λ:0.2λ.Figure 11
Position error of the left and right half arrays in theY-direction.Figure12 shows the projection of the three-dimensional pattern in the UV space when there is a 0.2λ error in the left and right half of the antenna. Figure 13 shows the simulation results of the antenna elevation plane pattern when the left and right half arrays have different errors in the Y-direction. The azimuth plane has little influence. Figure 14 shows the changes in the technical indicators of the antenna azimuth plane when the left and right half arrays have different errors in the Y-direction.Figure 12
Change of the position error pattern of the left and right half arrays in theY-direction 0.2λ.Figure 13
Change of the azimuth pattern caused by the position error of the left and right half arrays in theY-direction.Figure 14
Changes in technical indicators caused by the position error of the left and right half arrays in theY-direction. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)


## 5.1.3. Subarray Error inXY-Direction

Figure15 shows the projection of the three-dimensional pattern in the UV space when there is a 0.1λ error in the left and right half of the antenna.Figure 15
Position error pattern changes of the left and right half arrays in theXY-direction 0.1λ.Table3 shows the statistics of technical indicators such as antenna gain and sidelobe level when the XY-directions change at the same time.Table 3
XY-direction changes at the same time technical indicator statistics.
No.Directional map sectionX-axisY-axisGainSidelobe level3 dB beamwidthBeam pointing1Azimuth plane0028.22−20.338020.0528.22−20.338030.128.22−20.338040.05028.26−19.438050.0528.26−19.438060.128.26−19.438070.1028.30−18.627.8080.0528.30−18.627.8090.128.30−18.627.8010Elevation plane0028.22−20.296.60110.0528.22−20.306.60120.128.22−20.316.60130.05028.26−20.296.60140.0528.26−20.306.60150.128.26−20.316.60160.1028.30−20.296.60170.0528.30−20.306.60180.128.30−20.316.60From the above simulation results, it can be concluded that the gaps produced by the left and right half arrays in theX-direction will affect the distribution of antenna radiation power in space, and the impact on the sidelobe level is more obvious. As the gap increases, the electrical aperture of the planar array becomes larger, the gain becomes larger, the beam width becomes narrower, and the beam direction does not change. The gaps produced by the left and right half arrays in the Y direction will affect the distribution of antenna radiation power in space, causing tilt rotation, which has a significant impact on the far-area sidelobe level. As the gap increases, the electrical aperture of the planar array becomes larger, the gain becomes larger, the beam width becomes narrower, and the beam direction does not change. When the XY-direction changes at the same time, it is a combination of the above two situations.In practical applications, according to the principle that the sidelobe level deterioration is not greater than 1 dB, the error in theX-direction should be ≤0.05λ, and the error in the Y-direction should be ≤0.1λ.


## 5.2. Subarray Deformation Error

Assuming that a vibration load is applied during the operation of the radar, under the constraints of the four corners, the antenna array will undergo symmetrical and asymmetrical deformation. For the planar rectangular active phased array antenna, the deformation of its subarray is usually saddle-shaped deformation. The constraints are shown in Figure16. A typical saddle-shaped deformation is listed below, the effect of antenna deformation on the antenna pattern is simulated and analyzed, and the main technical indicators of the antenna are counted.Figure 16
Matrix constraint condition.


### 5.2.1. Saddle-Shaped Deformation of Arrays

Figure17 shows the saddle-shaped deformation of the antenna array. The normalized deformation data of the center of each antenna element is shown in Table 4.Figure 17
Saddle deformation.Table 4
Normalized deformation data at the center of the antenna element.
Unit no.Normalized data10.055920.2130.38440.5650.72260.85570.9518191100.951110.855120.722130.56140.384150.21160.055917−0.065180.0901190.257200.423210.575220.701230.791240.837250.837260.791270.701280.575290.423300.257310.090132−0.06533−0.21734−0.0543350.113360.275370.422380.543390.629400.673410.673420.629430.543440.422450.275460.11347−0.054348−0.21749−0.38150−0.20951−0.0381520.124530.269540.388550.472560.515570.515580.472590.388600.269610.12462−0.038163−0.20964−0.38165−0.53866−0.35867−0.18268−0.0176690.128700.245710.328720.371730.371740.328750.245760.12877−0.017678−0.18279−0.35880−0.53881−0.67282−0.48583−0.30584−0.138850.0083860.126870.209880.251890.251900.209910.126920.008393−0.13894−0.30595−0.48596−0.67297−0.76998−0.57899−0.395100−0.225101−0.07811020.04021030.1231040.1661050.1661060.1231070.0402108−0.0781109−0.225110−0.395111−0.578112−0.769113−0.821114−0.627115−0.442116−0.271117−0.123118−0.00481190.07811200.1211210.1211220.0781123−0.0048124−0.123125−0.271126−0.442127−0.627128−0.821129−0.821130−0.627131−0.442132−0.271133−0.123134−0.00481350.07811360.1211370.1211380.0781139−0.0048140−0.123141−0.271142−0.442143−0.627144−0.821145−0.769146−0.578147−0.395148−0.225149−0.07811500.04021510.1231520.1661530.1661540.1231550.0402156−0.0781157−0.225158−0.395159−0.578160−0.769161−0.672162−0.485163−0.305164−0.1381650.00831660.1261670.2091680.2511690.2511700.2091710.1261720.0083173−0.138174−0.305175−0.485176−0.672177−0.538178−0.358179−0.182180−0.01761810.1281820.2451830.3281840.3711850.3711860.3281870.2451880.128189−0.0176190−0.182191−0.358192−0.538193−0.381194−0.209195−0.03811960.1241970.2691980.3881990.4722000.5152010.5152020.4722030.3882040.2692050.124206−0.0381207−0.209208−0.381209−0.217210−0.05432110.1132120.2752130.4222140.5432150.6292160.6732170.6732180.6292190.5432200.4222210.2752220.113223−0.0543224−0.217225−0.0652260.09012270.2572280.4232290.5752300.7012310.7912320.8372330.8372340.7912350.7012360.5752370.4232380.2572390.0901240−0.0652410.05592420.212430.3842440.562450.7222460.8552470.951248124912500.9512510.8552520.7222530.562540.3842550.212560.0559Based on the normalized deformation, the maximum deformation of the antenna array was set, and the azimuth and elevation patterns of the antenna were simulated when the maximum deformation became0:0.025λ:0.4λ. Figure 18 shows the changes in the pattern of the array azimuth plane under saddle-shaped deformation. Figure 19 shows the changes in the main technical indicators of the array azimuth under saddle-shaped deformation. Figure 20 shows the changes in the pattern of the array elevation plane under saddle-shaped deformation. Figure 21 shows the changes in the main technical indicators of the array elevation plane under saddle-shaped deformation. Table 5 shows the changes in the main technical indicators of the antenna under the saddle-shaped deformation of the array.Figure 18
Azimuth plane of saddle-shaped deformation.Figure 19
Changes in technical indicators of azimuth plane caused by saddle-shaped deformation. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)Figure 20
Elevation plane of saddle-shaped deformation.Figure 21
Changes in technical indicators of elevation plane caused by saddle-shaped deformation. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)Table 5
Statistics of saddle-shaped deformation technical indicators.
No.Directional map sectionMax deformationGainSidelobe level3 dB beam widthBeam pointing1Azimuth plane028.22−20.338.000.0020.02528.20−20.098.000.0030.0528.14−19.408.000.0040.07528.05−18.378.000.0050.127.92−17.128.200.0060.12527.75−15.738.200.0070.1527.55−20.908.200.0080.17527.31−20.758.400.0090.227.02−20.548.400.00100.22526.71−20.248.600.00110.2526.35−19.868.800.00120.27525.95−19.379.000.00130.325.51−18.789.200.00140.32525.04−18.099.600.00150.3524.52−17.3010.000.00160.37523.96−16.4310.800.00170.423.37−15.4711.800.0018Elevation plane028.22−20.296.600.00190.02528.20−20.066.600.00200.0528.14−19.396.600.00210.07528.05−18.406.600.00220.127.92−17.216.800.00230.12527.75−15.916.800.00240.1527.55−21.056.800.00250.17527.31−20.996.800.00260.227.02−20.897.000.00270.22526.71−20.717.000.00280.2526.35−20.447.200.00290.27525.95−20.077.400.00300.325.51−19.597.400.00310.32525.04−18.997.600.00320.3524.52−18.298.000.00330.37523.96−17.498.400.00340.423.37−16.618.800.00From Figures18–21 and Table 5, the saddle-shaped deformation has a great influence on the azimuth plane and the elevation plane pattern of the planar array. As the amount of deformation increases, the antenna gain gradually decreases, the sidelobe level gradually increases, and the beam width gradually expands, which has almost no effect on the antenna beam direction. When the deformation amount is greater than or equal to 0.15λ, the azimuth plane and the elevation plane pattern begin to defocus, the pattern is no longer focused, and the first zero point disappears.


### 5.2.2. Analysis of Deformation Error Impact

In order to verify whether the saddle-shaped deformation will affect the scanning beam of the antenna, the azimuth plane and the elevation plane, respectively, scan the directional patterns of −30° under different maximum deformation variables, as shown in Figures22 and 23. It can be seen from the figure that the influence of saddle-shaped deformation on antenna beam scanning is basically the same as when the antenna is not scanning.Figure 22
Saddle-shaped deformation azimuth scanning−30° direction map.Figure 23
Saddle-shaped deformation elevation scanning−30° direction map.Integrating the technical indicators of the azimuth plane and the elevation plane, the maximum deformation error of the saddle-shaped deformation is≤0.05λ and the deterioration of the technical indicators is acceptable. In engineering practice, the maximum deformation error value should be controlled to be ≤0.05λ.


## 5.2.1. Saddle-Shaped Deformation of Arrays

Figure17 shows the saddle-shaped deformation of the antenna array. The normalized deformation data of the center of each antenna element is shown in Table 4.Figure 17
Saddle deformation.Table 4
Normalized deformation data at the center of the antenna element.
Unit no.Normalized data10.055920.2130.38440.5650.72260.85570.9518191100.951110.855120.722130.56140.384150.21160.055917−0.065180.0901190.257200.423210.575220.701230.791240.837250.837260.791270.701280.575290.423300.257310.090132−0.06533−0.21734−0.0543350.113360.275370.422380.543390.629400.673410.673420.629430.543440.422450.275460.11347−0.054348−0.21749−0.38150−0.20951−0.0381520.124530.269540.388550.472560.515570.515580.472590.388600.269610.12462−0.038163−0.20964−0.38165−0.53866−0.35867−0.18268−0.0176690.128700.245710.328720.371730.371740.328750.245760.12877−0.017678−0.18279−0.35880−0.53881−0.67282−0.48583−0.30584−0.138850.0083860.126870.209880.251890.251900.209910.126920.008393−0.13894−0.30595−0.48596−0.67297−0.76998−0.57899−0.395100−0.225101−0.07811020.04021030.1231040.1661050.1661060.1231070.0402108−0.0781109−0.225110−0.395111−0.578112−0.769113−0.821114−0.627115−0.442116−0.271117−0.123118−0.00481190.07811200.1211210.1211220.0781123−0.0048124−0.123125−0.271126−0.442127−0.627128−0.821129−0.821130−0.627131−0.442132−0.271133−0.123134−0.00481350.07811360.1211370.1211380.0781139−0.0048140−0.123141−0.271142−0.442143−0.627144−0.821145−0.769146−0.578147−0.395148−0.225149−0.07811500.04021510.1231520.1661530.1661540.1231550.0402156−0.0781157−0.225158−0.395159−0.578160−0.769161−0.672162−0.485163−0.305164−0.1381650.00831660.1261670.2091680.2511690.2511700.2091710.1261720.0083173−0.138174−0.305175−0.485176−0.672177−0.538178−0.358179−0.182180−0.01761810.1281820.2451830.3281840.3711850.3711860.3281870.2451880.128189−0.0176190−0.182191−0.358192−0.538193−0.381194−0.209195−0.03811960.1241970.2691980.3881990.4722000.5152010.5152020.4722030.3882040.2692050.124206−0.0381207−0.209208−0.381209−0.217210−0.05432110.1132120.2752130.4222140.5432150.6292160.6732170.6732180.6292190.5432200.4222210.2752220.113223−0.0543224−0.217225−0.0652260.09012270.2572280.4232290.5752300.7012310.7912320.8372330.8372340.7912350.7012360.5752370.4232380.2572390.0901240−0.0652410.05592420.212430.3842440.562450.7222460.8552470.951248124912500.9512510.8552520.7222530.562540.3842550.212560.0559Based on the normalized deformation, the maximum deformation of the antenna array was set, and the azimuth and elevation patterns of the antenna were simulated when the maximum deformation became0:0.025λ:0.4λ. Figure 18 shows the changes in the pattern of the array azimuth plane under saddle-shaped deformation. Figure 19 shows the changes in the main technical indicators of the array azimuth under saddle-shaped deformation. Figure 20 shows the changes in the pattern of the array elevation plane under saddle-shaped deformation. Figure 21 shows the changes in the main technical indicators of the array elevation plane under saddle-shaped deformation. Table 5 shows the changes in the main technical indicators of the antenna under the saddle-shaped deformation of the array.Figure 18
Azimuth plane of saddle-shaped deformation.Figure 19
Changes in technical indicators of azimuth plane caused by saddle-shaped deformation. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)Figure 20
Elevation plane of saddle-shaped deformation.Figure 21
Changes in technical indicators of elevation plane caused by saddle-shaped deformation. (a) Gain. (b) Sidelobe level. (c) Beam width. (d) Beam pointing.
(a)(b)(c)(d)Table 5
Statistics of saddle-shaped deformation technical indicators.
No.Directional map sectionMax deformationGainSidelobe level3 dB beam widthBeam pointing1Azimuth plane028.22−20.338.000.0020.02528.20−20.098.000.0030.0528.14−19.408.000.0040.07528.05−18.378.000.0050.127.92−17.128.200.0060.12527.75−15.738.200.0070.1527.55−20.908.200.0080.17527.31−20.758.400.0090.227.02−20.548.400.00100.22526.71−20.248.600.00110.2526.35−19.868.800.00120.27525.95−19.379.000.00130.325.51−18.789.200.00140.32525.04−18.099.600.00150.3524.52−17.3010.000.00160.37523.96−16.4310.800.00170.423.37−15.4711.800.0018Elevation plane028.22−20.296.600.00190.02528.20−20.066.600.00200.0528.14−19.396.600.00210.07528.05−18.406.600.00220.127.92−17.216.800.00230.12527.75−15.916.800.00240.1527.55−21.056.800.00250.17527.31−20.996.800.00260.227.02−20.897.000.00270.22526.71−20.717.000.00280.2526.35−20.447.200.00290.27525.95−20.077.400.00300.325.51−19.597.400.00310.32525.04−18.997.600.00320.3524.52−18.298.000.00330.37523.96−17.498.400.00340.423.37−16.618.800.00From Figures18–21 and Table 5, the saddle-shaped deformation has a great influence on the azimuth plane and the elevation plane pattern of the planar array. As the amount of deformation increases, the antenna gain gradually decreases, the sidelobe level gradually increases, and the beam width gradually expands, which has almost no effect on the antenna beam direction. When the deformation amount is greater than or equal to 0.15λ, the azimuth plane and the elevation plane pattern begin to defocus, the pattern is no longer focused, and the first zero point disappears.


## 5.2.2. Analysis of Deformation Error Impact

In order to verify whether the saddle-shaped deformation will affect the scanning beam of the antenna, the azimuth plane and the elevation plane, respectively, scan the directional patterns of −30° under different maximum deformation variables, as shown in Figures22 and 23. It can be seen from the figure that the influence of saddle-shaped deformation on antenna beam scanning is basically the same as when the antenna is not scanning.Figure 22
Saddle-shaped deformation azimuth scanning−30° direction map.Figure 23
Saddle-shaped deformation elevation scanning−30° direction map.Integrating the technical indicators of the azimuth plane and the elevation plane, the maximum deformation error of the saddle-shaped deformation is≤0.05λ and the deterioration of the technical indicators is acceptable. In engineering practice, the maximum deformation error value should be controlled to be ≤0.05λ.


## 6. Conclusions

Taking a 16 × 16 two-dimensional array antenna as an example, this paper constructed a millimeter-wave antenna array error theoretical model. The influence of random errors caused by production,XY-direction errors caused by assembly, and saddle-shaped deformation of the front on the electrical performance of the antenna under the condition that the size of the front is unchanged is analyzed. Through a large amount of data calculation, we have given the quantitative relationship between the influence of the error and the electrical performance of the antenna, drawn the influence relationship curve, and finally obtained the critical value. Engineers can refer to the analysis methods and conclusions in the article to estimate the antenna performance and put forward reasonable requirements for the processing tolerance of the antenna array within the allowable range of electrical performance. In follow-up research, we will further analyze the impact of other forms of errors on the electrical performance of the antenna and find the best correction method so that engineers can get a more optimized antenna design and assembly plan. The specific indicators obtained from the experiment are as follows:(1)
In engineering practice, we should try to eliminate the possible factors that produce random errors. At the same time, strictly control the production process and flow during the production and development process, so that the variance of the random phase and amplitude error does not exceed10°,0.5dB, to ensure a small impact on the array antenna pattern and meet the performance requirements of the array antenna.(2)
The errors caused by the assembly of the left and right half arrays, the gaps in theX-direction, have a more obvious impact on the sidelobe level, and the gaps in the Y-direction will affect the distribution of antenna radiation power in space, causing tilt rotation. When the XY-direction changes at the same time, it is a combination of the two situations. According to the sidelobe level deterioration being not more than 1 dB, the error in the X-direction should be ≤0.05λ, and the error in the Y-direction should be ≤0.1λ.(3)
Symmetrical deformations such as saddle-shaped deformations mainly affect the gain and energy distribution of the antenna and will not affect the beam direction. In engineering practice, to ensure that the sidelobe level deterioration is better than1dB, the maximum deformation should be better than 0.05λ.



---

*Source: 1874537-2021-09-25.xml*