Bias Correction to Antenna Frequency Response for Wideband Polarimetric Phased Array Radar

The frequency response of wideband polarimetric phased array radar antenna biases polarization measurement. A bias correction method named the iterative frequency division (IFD) algorithm is proposed in this paper to eliminate the effect of frequency response of the antenna system. The method is built on a new measurement bias model in which the effects of frequency response on the backscattering covariance matrix are analyzed. To demonstrate the effectiveness of the proposed method, a wideband microstrip patch array element is designed. Simulation results show that the biases of differential reflectivity ( Z D R ) and linear depolarization ratio ( L D R ) can be decreased to less than 0.1 dB and − 40 dB even for wide beam direction. The proposed method maintains a favorable polarization measurement accuracy for wideband polarimetric phased array radar.


Introduction
Recently, polarimetric phased array radar (PPAR) has attracted broad interest, owing to its potential in weather observation, air traffic control, and air surveillance [1][2][3][4][5][6][7]. However, some challenges, especially the requirement of high accuracy for polarimetric measurement, still hinder the development of the PPAR technology [7][8][9]. To be specific, the electric fields radiated from the horizontal (H) and vertical (V) ports are not orthogonal to each other caused by beam steering off the principal plane [3]. This leads to a cross-polar component where the polarization bias is generated. To address the problem, some calibration/correction methods have been studied. The bias correction technique for polarization measurement in PPAR was firstly put forward in [3], simulated in [9][10][11][12], and further updated in [13], and demonstrated in [5,14] using the large-scale testbeds. Thus far, most bias corrections focus on the narrow-band pulse. If the transmitted/received signal is wideband from a wideband polarimetric phased array radar (WPPAR), the pulse will be distorted due to the bandwidth effect of the antenna system [15]. The bias correction methods aforementioned cannot cover the wide bandwidth and inevitably cause the polarization measurement bias. Therefore, the wideband effect needs to be considered for accurate polarimetric measurement applications.
The wideband signal promises a higher resolution and detailed polarization profiles [15][16][17]. However, the frequency-dependent power amplifier (PA) and antenna system distort the output waveform in the WPPAR [18,19]. To circumvent the problem from PA, some digital pre-distortion schemes have been explored in [20][21][22][23]. Another challenge is from the transmitting/receiving antenna. The relative amplitude and phase of the far-field radiated by the wideband antenna are measured as a function of frequency and beam direction. The high fluctuation of amplitude and phase distort the wideband pulse [15,24]. Consequently, the distortion of received signals deteriorates the polarization measurement because the magnitude and phase of the signal lead to the bias of the scattering matrix. This is the problem we would like to address in this paper.
The polarimetric characteristic of the practical antenna varies with frequency, and the polarization non-orthogonality of H and V ports is also frequency-dependent. The polarimetric characteristic and polarization non-orthogonality are mutually coupled in the WPPAR. Hence, the frequency response of the wideband antenna is more complex for deteriorating the waveform fidelity for the transmitted/received pulse. In this paper, a polarization bias correction method named iterative frequency division (IFD) is proposed to address the above-mentioned coupled factors for the WPPAR. We analyze the process of transmission and reception of the pulse and derive the source of the polarization bias. The bias is reduced by a novel subband method involving the division of frequency band and the non-uniform selection of frequency point. The performance of the proposed method is verified by a designed wideband antenna element. Numerical simulations versus bandwidths and beam directions are investigated. The polarimetric variables, such as Z DR and LDR, are evaluated for the accurate polarization measurement.
The rest of the paper is organized as follows. In Section 2, the frequency response of the transmission process is deduced. In Section 3, the IFD algorithm for the wideband signal is provided to correct the scattering matrix. In Section 4, the polarimetric characteristics of the wideband microstrip antenna are analyzed, and the bias is quantified by the polarimetric variables. Simulation results and discussions are presented. Conclusions are presented in Section 5.

Theoretical Formulation of Bias
A typical signal transmission procedure between the antenna system and target, as illustrated in Figure 1, is considered. Horizontal (H) and Vertical (V) ports radiate and receive the individual electric field components. For the polarization measurement to a point target, the received voltage equation of the receiver can be expressed as [3,25] where matrix V is the received open-circuit voltage. matrix R and matrix T, describing the spacial polarization characteristics of the antenna, denote the transmitting pattern and receiving pattern, respectively. Generally, R is equal to T according to reciprocity principle. S (b) is the intrinsic backscatter matrix of the target. R T is the transpose of R. Rewriting Formula (1) in the form of a matrix gives (2) Formula (2) can only be applied to the narrow-band signal or the situation where the bandwidth effect of the antennas is neglected. For the wideband case, the equation becomes where f R is the radio frequency. The detailed derivation of Formula (2) is shown in Appendix A. For simplicity, the amplitude and phase imbalance of H and V channels, and the imperfect polarization isolation between H and V ports are neglected or idealized.
The voltage V( f R ) induced on the open-circuit terminal includes not only the scattering characteristic of the target but also the wideband properties of the transmitting and receiving antennas. It is frequency-dependent. We describe the characteristics of an antenna as a system transfer function mathematically, assuming that S (b) ( f R ) is constant. Given that the amplitudes of the transmitted signals are unit amplitudes and the relative phases are 0°, V( f R ) represents the amplitude and phase modulations to the signal's spectrum with the input frequency. The waveform information is carried by the emission pulse and is not related to V( f R ). Thus, it can be seen that the waveform of the received signals is modulated in the frequency domain. For each channel, V( f R ) is a complex number and denotes the relation between the source and received signal. The relation is expressed as follows: where X r ( f R ) is the output signal of the receiving antenna and X i ( f R ) is the spectrum of the input signal x(t). The overbar on symbol denotes that the signal is on a carrier, i.e., it has not yet been demodulated. It is worth mentioning that the transmitted/received signal should be amplified in practical radar systems by the high power amplifier (HPA) in transmitting antenna and the low noise amplifier (LNA) in receiving antenna. The nonlinear behavior of the amplifier introduces distortion in the output pulse. Dunn et al. used the digital pre-distortion technique to minimize the distortion of amplified signals that appear to be linearly amplified [20][21][22][23]. With the linear assumption on the amplifier in the WPPAR system, this paper focuses on the frequency response of the antenna systems. When the antenna system operates in the alternate transmission and simultaneous reception (ATSR) mode, all four components of the backscattering covariance matrix can be estimated from two or more pulses. Hence, the received signals can be expressed as For simplicity, Equation (5) is rewritten as It can be seen that the waveforms of received signals in four channels are modulated by transfer function V( f R ) in the frequency domain. The transfer function is determined by the impedance matching, gain, polarization matching et al. of the transmitting/receiving antenna.
The signal representations in the time and frequency domain are symmetrical. To facilitate the following derivation, the measurement bias is conducted in the frequency domain. In addition, the internal thermal noise and the external environmental noise are ignored herein. In a modern radar system, the carrier frequency of the received signal should be removed firstly [26], and the echo signal after downconversion is referred to as a baseband signal. The spectra of the baseband signals in matrix form can be written as where f is the frequency of the baseband signal. The four signals are sent to matched filters in the frequency domain, and the outputs of four channels from the matched filters are where X * ( f ) is the impulse response of matched filters. The output of the matched filter in the time domain is where IFT() denotes the Inverse Fourier Transform. It is clear in Equation (9) that both the amplitude and phase of four output signals are modulated by the complex matrix through the transmitting, scattering and receiving process. The modulation finally results in the distortion of echo signals, thereby deteriorating the measurement accuracy of the backscattering matrix. The maximum values of the filters' output are obtained and the measured scattering matrix can be expressed as where S If the projection matrix method in [3] is utilized to correct the measured scattering matrix, only one frequency point is considered for correction. For example, the projection matrix corresponding to the center frequency f 0 is selected as the correction matrix, and the expression is written as where C is the correction matrix, C = T −1 .
In fact, the correction matrix C involves the polarization characteristics including the cross-polar, polarization non-orthogonality. For the narrow-band case, C is related to the beam direction. However, for the wideband case, C is related to the frequency as well and can be denoted as C( f R ). Therefore, Equation (11) should be rewritten for the wideband antenna as Theoretically, each frequency point corresponds to one correction matrix, and each correction matrix corresponds to one corrected scattering matrix. The bias of polarization measurement varies with the input frequency, so we should select the frequency point that could minimize the correction bias. However, this expected frequency point that could satisfy the accurate polarization measurement does not necessarily exist because the selection of the correction matrix is related to the smoothness of the system transfer function V( f R ). Therefore, a bias correction method for the wideband case is proposed in this paper for the accuracy requirement of polarization measurement for WPPAR.

Correction Algorithm
The scheme of bias correction for the frequency response is displayed in Figure 2. The correction is implemented in the digital domain (the coloring regions in Figure 2). The correction matrix C( f R ) is frequency-dependent and is diverse from each frequency point. To reduce the number of frequency points, the correction can be done by using some specific points in the corresponding subbands, namely the frequency band division method. The measured scattering matrix S (p) is corrected by the correction matrices corresponding to each given subband. The sum of the estimated values is taken as the correction result. This method will improve correction effectiveness and maintain the polarization measurement accuracy.  The output of the matched filter in the frequency domain is denoted as X m [k] in the digital signal processing, where k = 1, 2, · · · , K, K is the total number of sampling points in the frequency band. We divide the bandwidth (specifically, 3 dB bandwidth herein) of X m [k] to multiple subbands. To be specific, the emphasis is on how to divide the frequency into subband and how to select one frequency point in each subband to obtain the corresponding correction matrix. The frequency point in each subband is selected to obtain the minimal estimation bias. Therefore, the problem in this paper can be described as a frequency point selection problem. The number of the selected frequency points is defined as N, the same as the number of subbands. The frequency division strategies and their computational complexity are shown in Figure 3. Several alternative ways are available.  • Strategy 1: Selecting the non-uniform frequency point in the individual non-uniform subband is equivalent to selecting the frequency point randomly. The number of permutation and combination for selecting n frequency points in K sampling points is C n K . To compare the effectiveness of the various permutation of the frequency point selection, the correction to the scattering matrix is implemented for each permutation. Thus, the correction times for n frequency points is nC n K . To obtain N optimal frequency points iteratively, the total correction times is ∑ N n=1 nC n K . The computation complexity has an exponential growth with the increasing number of frequency points. • Strategy 2: If the non-uniform point is selected in the individual uniform subband, the correction times K n , that is, the sampling points in each subband is K/N. The correction times in the total frequency band remains unchanged with N and is K. The total correction times is NK before the iteration operation is completed. Obviously, the computation complexity of this strategy is considerably lower than that of Strategy 1. If the frequency response does not fluctuate dramatically, the way dividing the transfer function V[k] uniformly with an equal bandwidth interval is reasonable. • Strategy 3: Uniform point selection in the uniform subband describes the case that the frequency point position is invariant for each subband, such as the center frequency. The strategy has fewer correction times. However, the frequency point is not necessarily optimal in each subband, and the number of needed frequency points will be greater than the two strategies mentioned above.
represents the nth subband of the matched filter output. According to Equations (8)-(10), the measured scattering matrix S (p) n of the nth subband can be derived. When the correction is done within the subband, the estimated scattering matrix is given by where C n [k] denotes the kth correction matrix for the nth subband, corresponding to the k frequency point. Now, the problem has become a choice matter of the frequency point, thereby obtaining the global optimal correction matrix. The iteration calculation is a candidate as a way to get the optimal results mathematically. If the frequency is divided iteratively, the number of subbands increases gradually and the number of frequency points has a continuous increase. The operation is repeated until the bias is corrected acceptably. The method is referred to as the iterative frequency division (IFD) algorithm. Using the IFD method, the bias ofŜ n [k] could be as low as possible. The formula in (13) is expanded for analyzing the polarization bias as follows: In addition, the four elements of the corrected scattering matrix are written as The sum of N corrected matrices is referred to as the scattering matrix of the wideband polarization measurement, where S n = S nhh S nhv S nvh S nvv .
Therefore, the estimated polarization scattering matrix (PSM) of the observed target can be computed.
To evaluate the performance of the correction matrix for polarization bias correction, we use the polarimetric variables including the differential reflectivity (Z DR ) and linear depolarization ratio (LDR) [25]. The scattering phenomenon is modeled as a reflection from a single point scatterer. The echo signal is received after the two-way propagations, that is, before and after scattering. The scattering matrix can be estimated through the echo by using correction methods, and the corrected Z DR scaled in dB is given by We describe the bias of differential reflectivity with The expected requirement for ∆Z DR is less than 0.1 dB for the phased array radar system [14]. The scattering matrix of the target in this research is a unit matrix. Thus, S hh = S vv and Z DR = 0. Hence, the bias of differential reflectivity is equal to the differential reflectivity after correcting, namely, Consider the assumption that the intrinsic scattering matrix of the target is a unit matrix, namely, S hv = S vh = 0, the measured LDR is strictly the bias of LDR. Hence, the ∆LDR of H channel after correcting scaled in dB is given by We assume that the constraint condition of the polarization variable measurement is ∆LDR h ≤ −40 dB. Two constraint conditions are utilized to determine how many subbands it needs to be divided into, and the two polarization variables imply the relative relations of four components in the scattering matrix. To increase the data comparability, we normalize the ∆Z DR and ∆LDR h , and the sum of normalized data ∆Z nor DR and ∆LDR nor h is written as The framework of the IFD algorithm for bias correction is displayed in Algorithm 1. The input of the algorithm is the matched filter output X m [k], system transfer function V[k]. The output of the algorithm is the minimal number of the subbands N. The iterative operations start with n = 1. Firstly, we divide X m [k] to N sections uniformly, such as X m1 [k], X m2 [k], · · · , X mN [k]. Secondly, a time-frequency conversion is made for the nth subband X mn [k] by using the IFT. The scattering matrix S (p) n can be measured through detecting the peak value. Thirdly, S (p) n is corrected by multiplying the correction matrices corresponding to K n frequency points in the current subband.

Algorithm 1: IFD for bias correction
We select the frequency point which could minimize the δ as the expected one, and take the N corrected results as the total corrected results. Lastly, N increases progressively and repeats the operations above until ∆Z DR ≤ 0.1 dB and ∆LDR h ≤ −40 dB. Therefore, the corrected S of the observed target is computed and the number N of the subbands is determined.

Simulation Results and Analyses
In this section, we provide the polarimetric characteristics with the frequency and off-axis deviation angle based on a designed wideband polarimetric antenna. The numerical results and comparisons present the correction to the polarization measurement bias. Furthermore, the reasons for the number of frequency points selection are discussed.

The Polarization Characteristics of a Wideband Polarimetric Antenna
A wideband microstrip patch antenna is designed by High Frequency Structure Simulator (HFSS) software to demonstrate the frequency response. The design of this wideband antenna is relatively mature and applied abroad, and its radiation characteristics can be used to test the influence of frequency response on polarization measurement. This is an X-band microstrip patch antenna which adopts two layers of substrates. The dielectric constant of the substrates r = 4.3, and the medium between the upper substrate and the lower substrate is air. The feeding configuration is a coaxial probe as displayed in Figure 4. Considering the antenna's voltage standing wave ratio (VSWR), as shown in Figure 5a, the bandwidth is up to 1.77 GHz when VSWR ≤ 2.  The electric fields of H and V ports create cross-polar coupling (i.e., polarization non-orthogonality) when the beam direction is away from the broadside (x-axis) of the antenna array. The off-axis deviation angle, the beam deviating from the boresight direction, is indicated with (φ, θ). φ is the azimuth angle. θ is the angle subtended from the z-axis. Without loss of generality, we choose (φ, θ) in the three-dimensional space and assume that a complementary angle relationship is restrained between φ and θ, which is We use HFSS to compute the radiation fields E φ and E θ from the H and V ports which are activated alternately. Figure 6 shows the polarization state of the radiated fields from the H port. The points distribution on the Poincare sphere describes the variation of the polarization state. The investigation has manifested that polarization state changes with the bandwidth, that is, the amplitude of E φ and E θ is not invariable or flat and the phase of E φ and E θ is not linear. Meanwhile, the state changes across with the deviation angle.
The results reveal that polarization state changes with not only the frequency but also the deviation angle.
The polarization state of fields from V port is similar to that of the H port. In addition to the variation of the polarization state in each frequency point, the non-orthogonality state of fields from H and V ports also changes over the bandwidth and beam direction separately. The intercoupling between the polarization state and the non-orthogonality makes it more complicated in the bias terms, and both of them contribute to the frequency response. In fact, the transmitting and receiving process all have the frequency response, so it will affect the measurement twice. The influence of the complex polarization characteristics on the polarization measurement accuracy is investigated in the next subsection.

Bias Correction with Bandwidth and Deviation Angle
The simulations and tests are conducted on the array antenna element and the tests are carried out in the ATSR mode. The emission waveform is the linear frequency modulation (LFM) signal and the central frequency is 9.5 GHz. The bandwidth B increases gradually from 100 MHz to 1700 MHz with a 100 MHz interval, and the pulse duration τ from 1 µs to 17 µs with 1 µs interval. The chirp rate B/τ remains invariant, so the waveforms are consistent for pulses with different bandwidths. The trends of the ∆Z DR and ∆LDR h are investigated based on different bandwidths and beam directions.
To manifest that the superiority of the non-uniform frequency point selection in the subband (Strategy 2)is better than the uniform selection (Strategy 3), a comparison is conducted. Situations when n = 1 and n = 2 are presented in Figure 7. The upper mesh surface in Figure 7a describes the bias after polarization correction using Strategy 3 across with the bandwidth B and deviation angle φ, and the lower mesh surface using Strategy 2. It can be seen that the polarization variable ∆Z DR is lower when the frequency point selection is optimized. Similar results have been observed when N = 2, as illustrated in Figure 7b. It is worth noting that Strategy 3 when N = 2 does not result in lower correction bias than that when N = 1, so the effect of Strategy 3 is doubtful. However, Strategy 2 can reduce the ∆Z DR considerably. In addition to ∆Z DR , a similar situation exists for ∆LDR h and differential phase (Φ DP ). Therefore, it is essential to search the optimal frequency point, thereby obtaining the corresponding correction matrix. Thus, the polarization measurement bias could be reduced. In the following analysis and discussion, only Strategy 2 is considered. Comparison of the corrected results for different N is displayed in Figure 8. Figure 8a shows that the measured ∆Z DR is 0.56 dB greater than that when the system frequency response acts on the polarization scattering matrix without correction. When N = 1, namely a single matrix corresponding to a single frequency point is used for correction, the biases decrease significantly. When N is larger than 2, the change is slight. However, the biases increase gradually with the gain of bandwidth, which proves the effect of bandwidth on the polarimetric measurement accuracy. Similarly, the ∆LDR h in Figure 8b increases with bandwidth. As manifested in Figure 8, the corrected ∆Z DR is less than 0.1 dB over the whole bandwidth if N ≥ 2, and ∆LDR h is less than −40 dB when the bandwidth is up to 1700 MHz. If the frequency band is divided into more than two parts, the goal of accurate measurement can be achieved. In addition, the influence of the deviation angle on the bias correction is investigated. The polarization characteristics of the wideband array element, especially the non-orthogonal state of the dual-polarized antenna, are very complicated. The co-polarization and the cross-polarization vary with the deviation angle. Taking the bandwidth B = 1700 MHz, for example, the corresponding results are shown in Figure 9. Without any correction, the ∆Z DR is very sensitive to deviation angle (fluctuating from 0 dB to 0.7 dB), and the ∆LDR h keeps large values (around −10°). With IFD correction, the ∆Z DR and ∆LDR h decrease significantly. Additionally, the correction biases decrease with the increase of N. The correction using the matrix corresponding to a single frequency point cannot meet the accuracy requirement. The ∆Z DR is greater than 0.2 dB at several angles and ∆LDR h reaches −24.97 dB. The corrected ∆Z DR is less than 0.1 dB over the whole bandwidth if N ≥ 2, and the ∆LDR h is below −40 dB when N = 3. Thus, N = 3 should be considered for this case. ∆Z DR and ∆LDR h are designated as the constraints of bias correction. The simulations present the ∆Z DR and ∆LDR h , which denote the relationship of the amplitudes between co-polarization and cross-polarization components. To test the validity of the proposed algorithm, the phase of the corrected PSM is investigated. This paper focuses on the differential phase between the two elements of diagonal terms, and it can be corrected similarly, as displayed in Figure 10. The Φ DP after bias correcting is less than 3.15 • under different bandwidths and deviation angles. The bias will be so small that it can be ignored when the corrections are conducted within multiple frequency subbands.

Discussion
Numerical results and comparisons show that the greater the number of divided subbands is, the smaller the corrected bias will be. However, the computational complexity grows with the increase of the number of divided subbands. The simulations are carried out using MATLAB 2017b (MathWorks, Natick, United States) on a laptop with Intel Core CPU i7-8550U (Intel, Santa Clara, United States) at 1.99 GHz and 8 GB RAM. The computation time is 5.62 minutes, 8.02 minutes, 9.56 minutes and 10.86 minutes for N = 1, 2, 3 and 4, respectively. According to the given constrained polarization variables, three non-uniform frequency points are enough to make a bias correction. In addition, the time consumed for N = 3 is 13.6% lower than that for N = 4. Therefore, we divide the frequency band to three subbands and select three correction matrices corresponding to three frequency points.
In this case, the number of frequency points seems to be small. The reasons may be as follows. On one hand, some idealizations in bias modeling are made, such as the nonlinearity of power amplifier, the amplitude and phase imbalance of H and V channels. On the other hand, the parameters, such as electric field, gain and impedance matching are the function of frequency. These factors contribute to the frequency response of the antenna system which affects the bias correction of the polarization measurement, thereby determining the number of divided subbands. Therefore, N depends on the antenna system parameters. The frequency response of the designed wideband antenna is relatively smooth, as shown in Figure 5b,c. Furthermore, the antenna effective length of an ultrawideband antenna is presented in the frequency-domain in [27], whose amplitude is flat and the phase is linear approximately. Hence, the bias correction operation by several frequency points is enough in the ideal situation.

Conclusions
This paper presents a key consideration for the measured PSM bias that arises from the frequency response of the antenna system with wide bandwidth. We propose an IFD method to correct the polarization measurement bias. The measured PSM can be corrected by iteratively multiplying the correction matrices in each subband. The simulation results show that the iterative non-uniform selection of frequency point in each subband is more effective than uniform selection on obtaining the corresponding correction matrix. When the frequency band is divided into three or more parts, it can meet the measurement requirement when the deviation angle varies from 0 • to 45 • and bandwidth from 100 MHz to 1700 MHz. That is, ∆Z DR is less than 0.1 dB; ∆LDR h is less than −40 dB. The ramification of this paper is that the frequency band division is an effective way to improve the polarization measurement accuracy, especially in the context of wideband antenna system in WPPAR.
An isolated phased array element is analyzed in this paper, without involving the mutual coupling. Future research would investigate the whole array with multiple elements, and the mutual coupling between the elements is necessary to be considered for correcting the bias of polarization measurement.
functions of frequency and can also be written as E t1 ( f c + f ), E t2 ( f c + f ); f c is the carrier frequency and f is the baseband signal frequency. Let f R = f c + f for simplicity herein. Project E t1 ( f R ) and E t2 ( f R ) onto the polarization base (a φ ,a θ ), and the projection matrix T( f R ) can be expressed as The projection matrix T( f R ) is not necessarily a unit matrix, especially when the beam is away from the broadside of the antenna plane. The cross-polarization is even on the order of co-polarization. Hence, the incident filed at the target when H port is activated is Simultaneously, the incident filed when the V port is activated is The backscattered electric field in the back scatter alignment (BSA) convention [25] can be written as where S (b) is the intrinsic backscattering matrix of the target. r is the distance between antenna and target. k is the wave number, and k = 2πf R /c, c is the electromagnetic wave velocity. S (b) is assumed to be a unit matrix, and the receiving antenna is assumed to be the same as the transmitting array element. When the H port is energized, the backscattered electric field E s1 ( f R ) at the receiving array element can be written as As expressed in [25], the radiation electric field can be formulated in terms of the transmitter power P t , antenna gain G, the unit polarization vector e, and the intrinsic impedance Z 0 of free space as The equivalent length of the antenna is defined as Hence, h can be derived as If the H port with a vector effective length h 1 (assumed that h 1 = h) is used to receive pulse radiated from H port, the open-circuit voltage V 11 ( f R ) of the receiving antenna can be expressed as [25] V The cofactor c/(2Z 0 √ P t ) is used to normalize the open-circuit voltage. V( f R ) can be written as Similarly, the open-circuit voltage V 21 ( f R ) in the V port can be written as Thus, the open-circuit voltage V 1 ( f R ) of the received element can be written as where T T ( f R ) is the tranpose of T( f R ).
Simultaneously, the open-circuit voltages V 12 ( f R ) and V 22 ( f R ) when V port is activated are The open-circuit voltage V 2 ( f R ) is The total voltage V( f R ) of the dual-polarization element can be reduced as below: