An iterative convex optimization (ICO) algorithm is proposed to solve pattern synthesis problem under the framework of dual-polarized conformal arrays in this paper. The subproblems of shaping main lobe, optimizing side lobe, and suppressing cross-polarization component are summarized as a joint optimization problem. To solve this problem, the nonconvex constraint about main lobe is rewritten as a convex constraint, which will bring error. And an auxiliary phase function is introduced to correct this error alternatively. Due to the deviation between auxiliary phase and real phase of pattern function, a method minimizing the peak of the synthesis error over observation angles is effectively applied to further improve the performance of the method. Numerical examples show good pattern synthesis ability and convergence performance of the ICO method.
1. Introduction
Benefiting from their incomparable flexibility, conformal arrays have been widely applied to communication systems and aircrafts design, etc [1]. Compared with traditional planar arrays, the inconsistency of the elements’ pattern and polarization feature brings enormous challenges to the pattern synthesis problem for conformal arrays.
The main purpose of conformal array pattern synthesis problem is to obtain desired pattern through optimizing the array weight vector and polarization information. Intelligence optimization algorithms such as genetic algorithm (GA) [2] and simulated annealing [3] are commonly used to solve this problem. However, these articles ignore optimization for cross-polarization component. In [4], the particle swarm optimization (PSO) method is applied to synthesize the elliptical and linear polarization patterns for dual-polarized array with shaped main lobe and low side lobe and cross-polarization level. In [5], an improved differential evolution algorithm is proposed to achieve lower side lobe levels and suppress the cross-polarization. Nevertheless, these stochastic optimization methods are computationally intensive because of their randomness in global optimization.
Other convenient and flexible methods are also utilized to synthesize polarization patterns. In [6, 7], the iterative least-squares method is used to achieve desired main lobe with polarization optimization. In [8], based on the adaptive theory, a numerical pattern synthesis method is proposed, which abandons the traditional method of applying virtual interference but controls the co-polarization and cross-polarization components by two weighting functions. In [9], an enhanced alternate projection algorithm is proposed to synthesize practical conformal array with mutual coupling corrected. In [10], a fast 3D synthesis method is proposed to synthesize copolarized and cross-polarized patterns with the consideration of dynamic range ratio of excitations. Although the above algorithms have a fast synthesis speed, the solution is easy to fall into the local optimal value.
The convex optimization algorithm can simultaneously optimize multiple problems expressed as a series of linear and convex quadratic constraints and can obtain the global optimal solution. Thus, it is also applied to the pattern synthesis problem. In [11], convex optimization starts to be really used for pattern synthesis with interior-point methods. In [12], the authors linearize the nonconvex constraints and divide the power pattern optimization problem into some subproblems, which can be solved by iterative second-order cone programming (SOCP). In [13, 14], it is studied as a convex programming problem to jointly synthesize the polarization and spatial power pattern. In [15], the constraints about constant modulus and main lobe are relaxed as convex constraints to solve. In [16], the authors use semidefinite relaxation to solve the nonconvex constraint to shape beam and control side lobe and cross-polarization. Compressive sensing (CS) theory and convex relaxation are used together to solve sparse conformal array synthesis problem, as described in [17]. It should be noted that the dual-polarized conformal arrays have advantages for the polarization information processing due to the two individual feed points of each dual-polarized antenna. In [18], a dual-polarized patterns optimization problem is formulated as a convex form, and related parameters determination method is presented. However, it lacks effective verification of the ability about main lobe shaping. It is worth studying pattern synthesis problem for dual-polarized conformal arrays.
This paper presents an iterative convex optimization (ICO) algorithm to synthesize pattern for dual-polarized conformal arrays. Main lobe shaping, side lobe optimization, and cross-polarization component suppression are formulated as a joint optimization problem. In the process of finding solution, the nonconvex constraint about main lobe is rewritten as a convex form, which leads to error. In order to reduce the error, an alternating optimization process is designed by introducing auxiliary phase. Since the auxiliary phase can only be approximately equal to the real phase, this approximation would bring error. To further improve the synthesis performance, we find the peak error point over main lobe and increase its weight in the iteration. Finally, numerical examples show the good performance of the proposed method in pattern synthesis.
2. Problem Formulation
Consider a dual-polarized conformal array composed of N elements located on an arbitrary carrier surface, as shown in Figure 1. Ox,y,z and O′x′,y′,z′ represent the global Cartesian coordinate system and the local Cartesian coordinate system, respectively. Without loss of generality, it is assumed that the θ-polarized and ϕ-polarized components are generated by two individual feed points of each dual-polarized antenna.
Coordinate systems.
The two orthogonal components of the array radiation pattern can be written as ϕ-polarization pattern Eϕθ,ϕ and θ-polarization pattern Eθθ,ϕ:(1)Eθθ,ϕ=∑i=1NEiθθ,ϕexpjkviTeθ,ϕ,=aTθ,ϕEθθ,ϕ,Eϕθ,ϕ=∑i=1NEiϕθ,ϕexpjkviTeθ,ϕ,=aTθ,ϕEϕθ,ϕ,where Eiθθ,ϕ and Eiϕθ,ϕ are the contribution of the ith element to the θ-polarized and ϕ-polarized components, respectively. k=2π/λ is the wavenumber, vi=xi,yi,ziT is the position vector of the ith element, eθ,ϕ=sinθcosϕ,sinθsinϕ,cosθT is the unit direction vector, •T denotes the transposition, and aθ,ϕ=ejkv1Teθ,ϕ,…,ejkvNTeθ,ϕT is the ideal steering vector of the array.
Note that the θ-polarized and ϕ-polarized components of each antenna are inconsistent because of the curved structure of conformal array. Therefore, they cannot be superimposed directly. The electromagnetic field excited by each element is decomposed and projected to the unified global coordinate system. And Eiθθ,ϕ and Eiϕθ,ϕ can be obtained though some mathematical transformation, which are described in [19] in detail. The global to local Euler rotation matrix of the ith element is defined as Rgli, and Rlgi denotes the local to global Euler rotation matrix of the ith element. Then, Eiθθ,ϕ and Eiϕθ,ϕ can be written as(2)Eiθθ,ϕEiϕθ,ϕ=Tcsθ,ϕRlgiΤsciθi′,ϕi′wiθ′∗giθ′θi′,ϕi′wiϕ′∗giϕ′θi′,ϕi′,where Tcsθ,ϕ is the transform matrix from Cartesian coordinate to spherical coordinate in global coordinate system and Tsciθi′,ϕi′ is the spherical coordinate to Cartesian coordinate transform matrix in local coordinate system. θi′,ϕi′ is the corresponding angle of θ,ϕ in the ith local coordinate and can been obtained though Euler rotation. wiθ′ and wiϕ′ are excitation weights of two feed points, which control θ-polarized and ϕ-polarized patterns of the ith element, respectively. •∗ denotes conjugation. giθ′θi′,ϕi′ and giϕ′θi′,ϕi′ are θ-polarized and ϕ-polarized patterns of the ith element. Then, for all of the array elements, we have(3)E1θθ,ϕE1ϕθ,ϕ⋮ENθθ,ϕENϕθ,ϕ=T1θ,ϕ⋱TNθ,ϕGθ,ϕw∗=Tθ,ϕGθ,ϕw∗,where Tiθ,ϕ=Tcsθ,ϕRlgiΤsciθi′,ϕi′ and w=w1θ′,w1ϕ′,…,wNθ′,wNϕ′T. Gθ,ϕ=diagg1θ′θ1′,ϕ1′,g1ϕ′θ1′,ϕ1′,…,gNθ′θN′,ϕN′,gNϕ′θN′,ϕN′T. It should be noted that the size of Tiθ,ϕ is 2×2, which means that Eiθ or Eiϕ depends on both giθ′θi′,ϕi′ and giϕ′θi′,ϕi′. Through the above transformation, array pattern can be obtained by(4)Eθθ,ϕEϕθ,ϕT=aTθ,ϕEθθ,ϕaTθ,ϕEϕθ,ϕT=E1θθ,ϕE1ϕθ,ϕ⋮ENθθ,ϕENϕθ,ϕTaTθ,ϕ⊗I2T,=wHGθ,ϕTTθ,ϕaθ,ϕ⊗I2=wHa˜θ,ϕ,where “⊗” denotes the Kronecker product, I2 is an identity matrix with the size of 2×2, and(5)a˜θ,ϕ=Gθ,ϕTTθ,ϕaθ,ϕ⊗I2.
It should be noted that a˜θ,ϕ is a matrix of size 2N×2, which contains two steering vectors of both orthogonal polarization components. Although the above derivation is based on the linear polarization model, it can be easily extended to any generic elliptic polarization condition. For example, in the condition of circular polarization, the corresponding steering vector can be expressed as(6)a˜cpθ,ϕ=Gθ,ϕTTθ,ϕaθ,ϕ⊗P−T,where(7)P=22−22j2222j.
As long as we know the parameters of any polarization state, the steering vector of this polarization can be obtained.
The dual-polarized conformal array pattern synthesis problem aims at finding a weight vector w with the size of 2N×1 to achieve desired pattern shape. Specifically, the ϕ-polarized component is assumed to be the copolarization component, and the requirements can be described as follows:
The main lobe of beam Eϕθ,ϕ is close to the desired shape dθ,ϕ in the main lobe region Ωm.
The side lobe level of beam Eϕθ,ϕ is under a upper bound ρθ,ϕ in the side lobe region Ωs.
The cross-polarization Eθθ,ϕ level cannot exceed the bound qθ,ϕ.
Therefore, this problem can be formulated as the following optimization problem:(8a)minwε,(8b)s.t.Eϕθ,ϕ−dθ,ϕ≤ε,θ,ϕ∈Ωm,(8c)Eϕθ,ϕ≤pθ,ϕ,θ,ϕ∈Ωs,(8d)Eθθ,ϕ≤qθ,ϕ,θ,ϕ∈Ωm∪Ωs,where ε denotes the maximum error between Eϕθ,ϕ and dθ,ϕ over Ωm. And through minimizing this error, we can obtain an array pattern approaching the desired one.
3. Pattern Synthesis Method3.1. The Convex Optimization with Additional Auxiliary Phase
It is known that constraint (8b) is NP-hard because of the quadratic constraint. However, it should be noted that(9)Eϕθ,ϕ−dθ,ϕ≤Eϕθ,ϕ−ejφθ,ϕdθ,ϕ,and if only φθ,ϕ=∠Eϕθ,ϕ, the “=” in formula (9) holds.
Therefore, the constraint (8b) can be rewritten as(10)Eϕθ,ϕ−ejφθ,ϕdθ,ϕ≤ε,θ,ϕ∈Ωm,where dθ,ϕ is often deemed to be real-valued and φθ,ϕ∈−π,π is defined as an auxiliary phase function.
In actual processing, there are many sampling angles over main lobe and side lobe. Thus each constraint in problem (8) contains a series of constraints. To ease the convex optimization programming, norm programming can be introduced. The constraint (8b) can be expressed as(11)Eϕm−ejΨ⊙d2≤ε,where Eϕm=Eϕθm1,ϕn1,…,Eϕθm2,ϕn2 denotes the discrete sampling of copolarization pattern in the sampling angles of main lobe region, Ψ=φθm1,ϕn1,…,φθm2,ϕn2, d=dθm1,ϕn1,…,dθm2,ϕn2, “⊙” denotes the Hadamard product, and •2 denotes L2 norm.
Similarly, other constraints in problem (8) can also be written as sampling vector form. The constraint (8c) can be expressed as(12)Eϕs∞≤ρ,where Eϕs=Eϕθ1,ϕ1,…,Eϕθm1−1,ϕn1−1,Eϕθm2+1,ϕn2+1,…,Eϕθk1,ϕk2 denotes the discrete sampling of copolarization pattern in the sampling angles of side lobe region.
The constraint (8d) can be expressed as(13)Eθ∞≤q,where Eθ=Eθθ1,ϕ1,…,Eθθk1,ϕk2 denotes the discrete sampling of cross-polarization pattern and •∞ denotes L∞ norm. q and ρ are both small real values.
Then problem (8) can be rewritten as(14a)minwε,(14b)s.t.Eϕm−ejΨ⊙d2≤ε,(14c)Eϕs∞≤p,(14d)Eθ∞≤q.
Observe that when φθ,ϕ is fixed, problem (14) is a jointly convex problem with respect to w,ε. And if w is given, the phase φθ,ϕ can be obtained. However, before we obtain the optimal weight vector, we can only know an approximate phase. Therefore, an iterative process is introduced to make the approximate phase gradually approach the real one. Specifically, when w is given, the auxiliary phase φθ,ϕ can be updated as follows:(15)φi+1θ,ϕ=α∠Eϕi+1θ,ϕ+1−αφiθ,ϕ,φi+1θ,ϕ of next iteration is not only simply assigned by ∠Eϕi+1θ,ϕ, but also φiθ,ϕ of last iteration. The parameter α determines the ratio of φiθ,ϕ and ∠Eϕi+1θ,ϕ in the process of obtaining φi+1θ,ϕ, and it is generally ranging from 0 to 1. Therefore, an alternating optimization process can be applied.
3.2. The ICO Algorithm via Minimax Error
Although the algorithm described above can reduce error through additional auxiliary phase, it should be noted that the algorithm still has inevitable error because φθ,ϕ is not exactly equal to ∠Eϕθ,ϕ. Thus it is necessary to find another way to further reduce the error and get smaller ripple of main lobe based on the above method.
For the constraint (14b), the “worst point” can be defined as(16)θ,ϕworsti=argmaxθ,ϕ∈ΩmEϕθ,ϕ−dθ,ϕ.
Compared with [20], the modulus of Eϕθ,ϕ is used because dθ,ϕ is real-valued. In the obtained pattern Eϕ, the point which exceeds the bound d most in the main lobe region would be defined as the “worst point.”
Then the point θ,ϕworst would be added to the observation points of next iteration. The constraint (14b) can be expressed as(17)Eϕmi+1−ejΨi+1⊙di+12≤ε,Eϕmi+1=Eϕmi,Eϕθ,ϕworsti,where Ψi+1 and di+1 are similarly defined. After each iteration, the weight at the angle θ,ϕworst would increase by “1” automatically. Therefore, the error between the obtained pattern and desired pattern can be reduced by accumulating the weights of some observation angles.
Through several iterations, the “worst point” may be at the same angle many times. Note that when the “worst point” are too concentrated, the weights of these points could be too large, which may cause the errors of other points to become larger and then may lead to worse result. In many cases, these worst points are most likely to be concentrated around the transitional region between the main lobe and the side lobe. Therefore, measures need to be carried out to avoid the above condition. Suppose that in the previous i iterations, Li denotes the sum of errors in main lobe. Compare Li+1 with Li, and determine whether w will be updated. It can be specifically described as the following calculation procedure:
Search for the “worst point” θ,ϕworst, and we can obtain Ωmi+1, Ωmi+1=Ωmi,θ,ϕworsti.
Calculate(18)Li=Eϕm−d2.
If Li+1 < Li, update w. Otherwise, w remains unchanged.
After each iteration, when the sum of errors in main lobe is reduced, w would be updated. Accordingly, when w remains unchanged, Li+1 should be substituted by Li. Through the above process, the convergence of the ICO algorithm can be guaranteed.
The initial value L0 should be set reasonably because it has an impact on the convergence of the algorithm. In the preparation stage, we solve problem (14) without auxiliary phase function by convex optimization toolbox. ∠Eϕ0θ,ϕ can be obtained. Then the sum of errors over main lobe is set as L0. And ∠Eϕ0θ,ϕ can be used as φ0θ,ϕ. The flow chart of the proposed algorithm is shown in Figure 2, and the problem (14) can be solved by convex optimization solver in the toolbox, such as the CVX [21].
Algorithm flow chart.
4. Synthesis Examples
The numerical examples of pattern synthesis are presented here to verify the effectiveness of the proposed method. Patch antenna is one of the most frequently used antennas in conformal array. In the following simulation, the lowest order circular patch antenna model [22] is adopted, and its pattern can be written as follows:(19)gθθ′,ϕ′=J2πdλsinθ′−J0πdλsinθ′⋅cosϕ′−jsinϕ′,θ′∈0,π2,gϕθ′,ϕ′=J2πdλsinθ′λ+J0πdλsinθ′⋅cosθ′sinϕ′−jcosϕ′,θ′∈0,π2,gθθ′,ϕ′=gϕθ′,ϕ′=0,Otherwise,where J2 and J0 are the first kind of second-order and zero-order Bessel functions, respectively. A convex optimization method via linear programming (LP) [23] is applied for comparison, where the upper bound and lower bound are set. And the upper bound is defined as the desired pattern in the proposed method. Two typical conformal arrays with half-wavelength-spaced elements are used to verify the effectiveness of the algorithm in the following.
4.1. Semicircular Arc Array
In this subsection, pattern synthesis is carried out based on a semicircular arc array as shown in Figure 3. The goal of the first example is to achieve a flat-top beam of ϕ-polarization components. With the fixed elevation angle θ=90∘, the main lobe region of desired pattern is ϕ∈−30∘,30∘, and the side lobe region is ϕ∈−90∘,−30∘∪30∘,90∘. The cross-polarization level is optimized in the whole observation region. And the value of α is 0.5 in this experiment. Figure 4 shows the synthesis result of copolarization component and cross-polarization component. The initial pattern is obtained by solving convex problem without auxiliary phase, as described in setting the initial values. In the case of achieving similar side lobe level −40 dB and cross-polarization level −50 dB, the proposed ICO algorithm and LP method perform better than the initial pattern on main lobe shaping. By the ICO method, the ripple of gain level in main lobe is less than that of the LP method. The result shows that the optimization performance of the ICO algorithm is improved. And the corresponding elements weights are shown in Figures 5 and 6.
Schematic view of the conformal array with 30 elements.
Flat-top-shaped beam synthesis of a 30-element array. (a) Copolarization component Eϕ at θ=90∘. (b) Cross-polarization component Eθ at θ=90∘.
Normalized amplitudes of the elements weights.
Phases of the elements weights.
After 100 iterations, the count of the “worst points” in Ωmi+1 is plotted in Figure 7. Note that the sum of errors in main lobe stops reducing through 28 iterations. Therefore, there are only 28 “worst points” in Ωmi+1. It can be seen that the “worst points” are centered around ±20∘. The error caused by the transition region is relatively large. Thus the result is reasonable.
Count of the “worst point” (iterations = 100).
The sum of errors in main lobe Li is defined as an index related to convergence. Since Li is not involved in the LP algorithm, it is uniformly calculated by (18). According to the first experiment, the main lobe region is ϕ∈−30∘,30∘, and so is the calculation range of (18). And it should be noted that due to different desired patterns, the sum of LP algorithm (d = (dup + dlow)/2) error appears smaller than that of the proposed ICO method (d=dup), where dup denotes the upper bound and dlow denotes the lower bound. The convergence curve can be seen in Figure 8. After 28 iterations, the ICO algorithm converges, but LP method does not. Therefore, the proposed ICO method performs better on convergence.
Li versus the number of iterations.
Further experiments are carried out to analyse the performance of the proposed method with different α, as shown in Figure 9. The results show that α has an effect on the ripple of main lobe. If we want to get the optimal solution, it is necessary to select that by trial and error according to actual requirements.
Copolarization components with different α.
The second example is to shape a cosecant-square beam over main lobe region defined as 0∘,40∘ with the array shown in Figure 3. Simultaneously, the side lobe level and cross-polarization level should be optimized. Figure 10 shows the optimized patterns with α=0.2. When using the proposed ICO method, the main lobe is similar to the initial pattern. But side lobe level and cross-polarization level are about 2 dB lower than those of initial pattern. Although the LP method can achieve the similar side lobe level and cross-polarization level with the proposed algorithm, its ripple of main lobe is larger. To sum up, the proposed algorithm has a better overall performance on cosecant-square beam shaping, side lobe level, and cross-polarization level suppression.
Cosecant-square-shaped beam synthesis. (a) Copolarization component Eϕ at θ=90∘. (b) Cross-polarization component Eθ at θ=90∘.
4.2. Semicylindrical Array
In this subsection, the synthesis of flat-top beam centered at 20∘ with a semicylindrical array of 2×40 elements is considered, as shown in Figure 11. With the fixed elevation angle θ=90∘, the main lobe and side lobe regions of desired pattern are ϕ∈−10∘,50∘ and ϕ∈−90∘,−10∘∪50∘,90∘, respectively. The cross-polarization component is suppressed in the whole observation region. And the value of α is 0.7. Figure 12 shows that the ICO and LP methods can achieve similar side lobe level −40 dB and cross-polarization level −50 dB with the initial pattern. The proposed ICO and LP methods perform better than the initial pattern on main lobe shaping. In addition, the ICO method can achieve pattern with less ripple of main lobe than LP method. In the case of using different arrays and setting main lobe centered at different directions, the above result shows the good adaptability of the proposed ICO method.
Schematic view of the conformal array with 2×40 elements.
Flat-top-shaped beam synthesis of a 2×40-element array. (a) Copolarization component Eϕ at θ=90∘. (b) Cross-polarization component Eθ at θ=90∘.
Considering all the experiment results, the comparative algorithm can also have a good performance, but the proposed method performs better on ripple of main lobe and cross-polarization level. Thus the algorithm improvement achieves a good effect.
5. Conclusion
This paper proposes an ICO algorithm to synthesize pattern for dual-polarized conformal arrays. A series of constraints are formulated as a joint optimization problem. In the process of solving this nonconvex problem, the convex form for the constraint of main lobe leads to error. Then an auxiliary phase function is introduced to reduce this error through an alternating optimization process. But the auxiliary phase cannot be completely equal to the real phase, and the error occurs in the main lobe. To get better performance, we find the “worst point” with the peak error in main lobe region, and then the performance of main lobe shaping is improved by increasing the weight at “worst point” iteratively. Finally, numerical experiments are carried out, and results show the well convergence and excellent capability of beam shaping of the proposed ICO method.
Data Availability
The image data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
MortonT. E.PasalaK. M.Pattern synthesis of conformal arrays for airborne vehiclesProceedings of the IEEE Aerospace ConferenceMarch 2004Big Sky, MT, USA10301039YangJ. O.YuanQ. R.YangF.ZhouH. J.NieZ. P.ZhaoZ. Q.Synthesis of conformal phased array with improved NSGA-II algorithm200957124006400910.1109/tap.2009.20267142-s2.0-72649086091FerreiraJ. A.AresF.Pattern synthesis of conformal arrays by the simulated annealing technique199733141187118910.1049/el:199708382-s2.0-0031551270LiM.WangX.LiY.DaiH.Polarization beam pattern synthesis based on particle swarm optimization20129211648165310.1587/elex.9.16482-s2.0-84870454920LiR.XuL.YinY.Low-cross-polarisation synthesis of conformal antenna arrays using a balanced dynamic differential evolution algorithm201711131855186210.1049/iet-map.2017.04612-s2.0-85031931923VaskelainenL. I.Iterative least-squares synthesis methods for conformal array antennas with optimized polarization and frequency properties19974571179118510.1109/8.5969122-s2.0-0031192533VaskelainenL. I.Constrained least-squares optimization in conformal array antenna synthesis200755385986710.1109/tap.2007.8918602-s2.0-34047200421ChenS.-L.QinP.-Y.GuoY. J.LiuY.YouP.Generalized 2-D numerical pattern synthesis algorithm for low cross polarization and low sidelobe synthesis2017162578258110.1109/lawp.2017.27348412-s2.0-85028878262ZhaoF.ChaiS.QiH.Hybrid alternate projection algorithm and its application for practical conformal array pattern synthesis201223562563210.1109/jsee.2012.00078ComissoM.VescovoR.Fast co-polar and cross-polar 3D pattern synthesis with dynamic range ratio reduction for conformal antenna arrays201361261462610.1109/tap.2012.22248342-s2.0-84873336697LebretH.BoydS.Antenna array pattern synthesis via convex optimization199745352653210.1109/78.5584652-s2.0-0031101492TsuiK. M.ChanS. C.Pattern synthesis of narrowband conformal arrays using iterative second-order cone programming20105861959197010.1109/tap.2010.20468652-s2.0-77953084142LiM.ChangY.LiY.DongJ.WangX.Optimal polarised pattern synthesis of wideband arrays via convex optimisation20137151228123710.1049/iet-map.2012.06552-s2.0-84890201787FuchsB.FuchsJ. J.Optimal polarization synthesis of arbitrary arrays with focused power pattern201159124512451910.1109/tap.2011.21654922-s2.0-80155123266LiuY.BaiJ.XuK. D.Linearly polarized shaped power pattern synthesis with sidelobe and cross-polarization control by using semidefinite relaxation20186663207321210.1109/tap.2018.28167822-s2.0-85044090827MassaA.RoccaP.OliveriG.Compressive sensing in electromagnetics - a review201557122423810.1109/map.2015.23970922-s2.0-84925870018CaoP.ThompsonJ. S.HaasH.Constant modulus shaped beam synthesis via convex relaxation20171661762010.1109/lawp.2016.25941412-s2.0-85018466040HuW.WangX.LiY.XiaoS.Synthesis of conformal arrays with matched dual-polarized patterns2016151341134410.1109/lawp.2015.25084382-s2.0-84964881936LiuC.DingZ.LiuX.Pattern synthesis for conformal arrays with dual polarized antenna elementsProceedings of the International Congress on Image and Signal ProcessingOctober 2014Dalian, China968973WangF.YangR.FrankC.A new algorithm for array pattern synthesis using the recursive least squares method2003108235238CVX Research, Inc., “CVX: MATLAB software for disciplined convex programming, Version 2.1,” Dec. 2018 [Online]. Available: http://cvxr.com/cvxJamesJ. R.HallP. S.WoodC.1981New York, NY, USAPeter PeregrinusFuchsB.Shaped beam synthesis of arbitrary arrays via linear programming2010948148410.1109/lawp.2010.20512102-s2.0-77953274545