An Iterative Technique for the Synthesis of Active Antenna Oscillator Arrays

A design procedure for the synthesis of a coupled active antenna oscillator array is presented. Such an array is synthesized by deriving two sets of parameters: the radiators’ positions and the oscillators’ outputs. The outputs are used to excite the radiators. Minimization of the mean square error between the desired pattern and the resulting one is made. Synthesis starts from an initial array, which is perturbed iteratively by varying simultaneously the element excitations and positions. In the iteration, the first variation of the cost function is set equal to zero. The final array results from the last iteration, where the stopping criteria are met. The procedure designs simultaneously both the antenna and the attached coupled oscillator array providing viable solutions. The second by properly configuring the tuning parameters through the use of closed-form formulas. The resulting arrays are shown to comply with the desired pattern and the nonlinear dynamics thus proving the validity of our method.


Introduction
It is well known that the angular distribution of the radiated power from an antenna array depends on the geometry and the excitation of the elements [1].In a coupled active antenna oscillator linear array consisting of identical and parallel elements, the pattern is expressed in the following form: where the excitation V k is the complex voltage of the oscillator output, (V k = | V k | amplitude and φ k = ∠ V k the respective phase), f k (θ) is the pattern, and z k is the abscissa of the kth element, respectively.λ is the free-space wavelength at the operating frequency.
The two key sets of parameters that are important in the design are the boxed terms of (1).The antenna array synthesis mainly focuses on the excitation [2][3][4], while efforts aiming at extracting the positions of the array elements are scarce [5].In our present study, we cope either with one of the above or with both sets simultaneously.A perturbation method that provides iteratively the excitation and the placement of the elements in a quest of matching the resulting far-field with the desired one is applied.The whole synthesis is formulated as an optimization problem over a vector containing the excitation and the element positions.The cost function is the mean square error between the desired and the attained pattern.
Apart from attaining the necessary parameters of the antenna array, it is of equal importance the problem of feeding the array by a proper network.In the literature, one can find several approaches to the feeding networks [1], containing corporate feeds and feeding matrices [6].Apart from passive networks, there are others that utilize active circuitry.Among them, the active feeding networks consisting of coupled oscillator arrays [7] are of specific interest.It has been found that such devices allow for the manipulation of the phase distribution without additional phase-shifting circuitry, suggesting a potential for low-cost beam-scanning systems.In the past, most of the studies on coupled oscillator arrays analyze the problem of beam steering [8][9][10] and of null control of the pattern, [11][12][13].In the current work, we extend the functionality of the coupled oscillator arrays to more general pattern synthesis problems.It is noticed that there are cases where a uniform excitation of a linear phased array and a nonuniform placement is enough to solve the synthesis problem.In this case, constraining the oscillators' output voltage to a common value reduces the nonlinear dynamics of the system.Alternatively, one can keep the phases equal and seek a design for different amplitudes.The required excitation distribution for all the simplifying cases, and for the general one (see the design examples), can be readily implemented with proper tune of control parameters, such as the varactor capacitance and the PIN conductance of the oscillators.The relative closed-form formulas are presented for the case of the mutually synchronized coupled oscillator array.

International Journal of Antennas and Propagation
The paper is organized as follows.In Section 2, the coupled oscillator array (COA) theory is reviewed, and the key aspects related to the present study are underlined.The system that contains both, the antenna and the oscillators, is presented in Section 3, where closed-form formulas are derived to design the mutually synchronized oscillator array.The formulation and the solution of the radiation pattern shaping optimization problem for the above active antenna array are given in Section 4. In Section 5, design examples that verify the effectiveness of the method follow.Finally, some concluding remarks are given in Section 6.

Single Oscillator Description.
A model that has been widely used for microwave oscillator circuits is the canonical oscillator model in a parallel current-sum form which leads to the Van der Pol equation for certain nonlinearities [14,15].The equivalent circuit of the oscillator can be extracted from measurements [16] or from simulation [17], and is given in Figure 1.
The oscillator is assembled by three distinct parts: the active device, the tuning tank, and the resonator tank.The active device provides the nonlinear voltage dependent negative conductance.The tuning tank is responsible for the adjustment of the tuning parameters in order to feed the attached antenna element with the proper excitation.Last, the resonator tank includes the parallel resonance circuit of the oscillator and the antenna elements.
In detail, the active device is represented as a negative conductance as follows: where |V | is the amplitude of the output voltage, and a, b are constants characterizing the nonlinear device.The voltage-dependent susceptance of the nonlinear device is denoted by B D (|V |).Usually, the susceptance is independent on the device characteristics and can be neglected.In this case, which is the one that we follow, the current of the nonlinear device obeys ( The reactive energy of the circuit is stored into the LC parallel circuit giving a free-running circular frequency of for the kth oscillator.
In our case, the capacitance C t would be a varactor that is controlled by a voltage.This voltage will be the tuning parameter which for the kth oscillator is V t k .The varactor capacitance follows the well-known equation [18]: where V o = 1.2 V (for a GaAs device), C o,t = 1 pF is the capacitance for V o tuning voltage, and Γ is the varactor exponent (usually taken equal to unity, a choice also followed here).Alternatively, other types of voltage controlled capacitors [19,20] or even voltage controlled inductors [21] can be used instead.The conductance G t represents additional losses in the oscillator circuit.This would be the second (together with the C t,k (V t k )) tuning parameter.In our case, the conductance G t would be a PIN diode variable resistor controlled by its forward current I t k .The PIN conductance follows the wellknown equation [22,23]: Here, we use v o = 20 mV.Alternatively, other types of variable resistors can be used instead [24].
To find the total conductance, one should add to G t the input conductance G A of the radiating element of the array.
Let us now compute the total admittance Y of the circuit: International Journal of Antennas and Propagation 3 We assume that the resonant circuit of the oscillator is of high Q.Thus, the output voltage follows, nearly, a sinusoidal waveform with the following phasor representation: 2.2.Coupled Oscillator Array Description.To complete the analysis, we suppose (see Figure 2) that we have N oscillators which are joined together using a linear coupling network [25].
The governing equation reads where Y i is given by (7), and V i is defined in (8) as the complex output voltage of the ith oscillator.Also, is the mutual admittance between the output nodes of the ith and jth oscillators.Using the above, it can be shown [16][17][18][19][20][21][22][23][24][25][26][27] for the voltage amplitude V k of the kth oscillator that where with Often the following substitution is made: where μ k is a parameter that characterizes the oscillator, while α o is the free-running amplitude of the output voltage.This weak quadratic nonlinearity is the simplest nonlinearity that approximates the behavior of a wide range of active devices.Apart from the time evolution of the voltage amplitudes given in (11) for the phase, we have Equation ( 11) in conjunction with ( 14) fully describes the coupled oscillator array.The above should provide a framework for the design of an antenna oscillator array.

Design of the Mutually Synchronized Array
In Figure 3 [28], an antenna array with the corresponding oscillators is presented.In this mutually synchronized array, each oscillator is bilaterally coupled to the neighboring array elements.The configuration of Figure 3was first proposed by Stephan [8,15] and named as "interinjection-locked" oscillators.Mutual coherence is achieved via the injectionlocking process.The steady-state phase relationships of each oscillator depend on its neighbors' amplitude and phase.
A nonlinear system of simultaneous equations must be solved to determine the self-consistent amplitude and phase relationships [26].In such systems, ( 11) and ( 14) take the form where κ lk are different from zero only when |l − k| ≤ 1 since only the nearest neighbor coupling is assumed.There are various sets of parameters in (15).Those are the oscillator detuning, (ω − ω o,k ), the un-normalized coupling strengths, (κ lk = |G lk + jB lk |), and the coupling phases, (Φ lk = ∠(G lk + jB lk )).There are also the capacitance and the conductance of each one of the oscillators.Since our aim is to manipulate the amplitude and the phase distribution of the active antenna array, any set of parameters could be used for control.Here, we eliminate the coupling phases and we use the conductance of the PIN diodes and the capacitance of the varactors to control the steady state of the amplitudes and phases.Before we present the followed procedure, some issues relative to the coupling network should be addressed.To couple the oscillators, we use resistively loaded lines [29], whose characteristic impedance Z 0 and phase Φ of the electrical length are For such lines, it is where n i = 2 except for the edge elements where it is equal to one.Using this coupling network, closed-form expressions for the tuning parameters can be derived.So, for Y i j given in International Journal of Antennas and Propagation (17), (15), for the steady state, where the time derivatives vanish, reads Using the solution for the tuning parameters,

Oscillator k
Oscillator N Figure 2: The coupled oscillator array.
Before closing the section, the design procedure and check of the stability of the solution for the coupled oscillator array are listed.So, to design the mutually synchronized oscillator array in order to produce a given output voltage distribution (amplitudes and phases), the following steps are implemented.
(i) The given output amplitudes and phases, V k , φ k , are inserted into (19) (15) provides the evolution of the amplitudes and phases with time.
The key step of computing the proper V k , φ k is addressed next.

The Synthesis Method
The system of Figure 3 with the restriction of nonvarying interelement distance, (Δz = z k − z k−1 = const), has been used to control linear phased arrays in the limited cases of beam steering and pattern null in certain directions [8][9][10][11][12][13].
In the current study, we attempt a more general (not limited in certain cases) pattern synthesis.The system is designed by Radiator 1 Radiator N Oscillator k Oscillator N z-axis x-axis y-axis z k providing the proper values for the positions of the radiators and for the voltage outputs of the oscillators.It is obvious that the proposed design offers specific advantages due to the degrees of freedom that make use.
To start the design of the array, we derive the cost function C of the following form: where We define From ( 23), (22) becomes where F d (θ) is the desired pattern, F(θ) is the resulting pattern, and F d − F 2 is the mean square error.By taking into account the output voltage of the coupled oscillators, the resulting pattern can be expressed as where We suppose a perturbation of the excitation of the oscillator outputs and of the position of the elements of the array.Due to the perturbation of the excitation, the new output of the kth oscillator becomes V k +ΔV k ; φ k +Δφ k .Also the new position of the kth element is z k + Δz k .The pattern of the array results in where F 0 (θ) is the initial pattern of the array.It is shown in the Appendix that after assembling the perturbed values of the output voltage, (ΔV k , Δφ k ), of the oscillators and the position, (Δz k ), of the elements in an unknown vector [Δ], the solution of the optimization problem leads to where [A] and [B] are proper matrices that include information relative to the specific problem.Expression (27) gives the appropriate perturbed values of the output voltage of the oscillators and the position of the elements.It is supposed that the desired pattern is not much different than the resulting one.In this case, the perturbation will produce the necessary changes in the excitation and the position of the elements.If F d − F 2 is large, then the perturbation technique is applied iteratively until the mean square error becomes less than a given threshold.For further details on the functionality and the problems relative to the convergence of the proposed iterative method, the reader is referred to [5].

Design Examples
In the design examples that follow, the array can be derived by varying the excitation (amplitude and phase) simultaneously with the element positions.The desired far-field  pattern is a Chebyshev one with HPBW = 15 degrees and SLL = −25 dB.We start from an array of 9 elements with uniform initial excitation and 0.45 λ interelement distance.We apply the synthesis method for three different cases.In case A, the excitations and the positions are to be derived.In case B, the phase of the excitations is set equal to zero, and in case C, the amplitude of the excitation is equal to unit.
For each one of the cases A, B, C, three different oscillators are employed to produce the respective far-field.The Van der Pol parameters of those oscillators are taken from published data (see Table 1).
The parameters of the coupling network are the same for all the examples as follows: 5.1.General Design Example.In this first example, both the element placement and the excitation are allowed to vary.After 20 iterations, the designed antenna array produces a pattern that is given in Figure 4.
The respective element placement and the excitations are given in Table 2.
The tuning parameters, V t k and I t k , of the oscillators are computed after inserting the last two columns of Table 2 into ( 19)- (3).For the present case, the relative results are presented in Table 3. (3) the initial, at t = 0, amplitudes and phases: In Figure 5, it is demonstrated that the designed coupled oscillator array can produce the necessary solution (see Table 2 and Figure 4) and furthermore, this solution is stable.In Figure 5(a), the voltage amplitudes of the nine oscillators versus time are given.Also, in Figure 5(b), the time evolution of the successive phase differences (output voltage phase of the kth oscillator minus the output oscillator phase of the (k − 1)th oscillator) are given.It is clear that the amplitudes and phases converge to the correct values (continuous horizontal lines).

Placement and Amplitude Design Example-Constant
Phases.From Table 2, it is shown that the phase of the elements is close to zero.Also, the final array is nonuniform in the excitation and the geometry.Namely, our solution is more complicated than the classical one for Chebyshev patterns [1,2].Since we try to have a technically simple feeding network, we will solve the same problem with no phase difference between the elements and with only three different excitations.In Table 4, the final placement and the excitation of the elements, after 20 iterations, are given.
Also in Figure 6, the pattern of the array is presented.In the pattern, a difference of less than 0.5 dB between the desired and the resulting one is observed.
Comparing Table 2 with Table 4, one can see that the second solution gives simpler results.
For the present case, the tuning parameters are given in Table 5.  Successive phase differences (degrees) In Figure 7, the time evolution curves for the amplitudes show the stability of the solution.
Note that comparing Figure 7 with Figure 5(a), different converging times are observed.In general, this fact is to be expected since the converging times, apart from the initial conditions V k (t = 0), φ k (t = 0), depend on the   tuning parameter values that are different for the two design examples.

Placement and Phase Design
Example-Constant Amplitudes.In this last case, we provide a design where the amplitudes are equal.The attained solution is given in Table 6.
The produced far-field pattern is given in Figure 8.
It is noted that in this case, the SLL criterion of −25 dB is not fulfilled, (SLL = −21.5 dB).The above shows that the phases of the elements are not enough to solve the problem.For the present case, the tuning parameters are given in Table 7.
For the resulting pattern, the respective evolution curves are given in Figure 9.It is obvious that in the array synthesis, one must tradeoff the simplicity of networks, the size of the array, and the number and type of the different excitations.From the presented cases, it is assessed that the general one is rather complex, while the constant amplitude case fails to comply  with the desired pattern.It seems that the constant phase case presents a better compromise between complexity and performance of the device under study.
Before closing the section, to evaluate the validity of the presented method and the viability of the listed results, some comments are in order.Both phase and amplitude control should be implemented with a relatively easy way.Varactor and PIN diodes are simple electronic elements that should be used.That is because varactor diodes are extensively used as tuning elements in conventional VCOs, and also PIN diodes are employed as variable resistors for vector modulator circuits.The issue that raises concern is the imprecise setting of the values of the tuning parameters and the effect that manufacturing tolerances will have on the overall performance of the system.Fabrication and tuning inconsistencies among individual oscillators in an array lead to differences in the free-running frequencies and amplitudes.This randomness causes phase and amplitude errors that are expected to result in pattern deformation manifested in beam-pointing and SLL errors [1,30].
To provide understanding of the relationships among random tuning error distribution and pattern deformation, a statistical study, focused on the sensitivity analysis, is provided through Monte Carlo simulation.The element placement, the tuning voltage, and the tuning current are treated as random vectors following normal probability distribution with means given by the previous analysis and varying values of variance.In each statistical experiment, after choosing the placement and tuning parameters, C t,k (V t k ) and G t,k (I t k ) are used to extract the time evolution of V k (t) and φ k (t).At the steady state, the far field of the array is computed.The results show that when only I t k is varied and if two or three dBs degradation on SLL performance are acceptable, then the variance should be no more than 10% of the parameters nominal value.This is shown in Figure 10 when all the patterns resulted from a 100-run relative Monte Carlo experiment are depicted.Nearly, for one and two orders of magnitude, less variation is needed for attaining the same results when position and phase are varied, respectively.This means that position and amplitudes setting problems are expected to be less harsh compared to the problems that have already encountered in setting the phase [31], in conventional VCOs.If the previously-stated performance degradation is unacceptable, properly chosen calibration procedure should be considered [32,33].

Conclusion
A design procedure for the synthesis of active antenna arrays is presented.Those arrays are sets of radiator oscillator couples with each element connecting through a coupling network to the others.Such a system presents the designer with two sets of parameters: the radiators' positions and the oscillators' outputs, the second used to excite the first.
The procedure designs simultaneously both the antenna and the attached coupled oscillator array providing viable solutions.In this way, the current work demonstrates how to control the dynamics of coupled oscillator arrays in order to achieve pattern shaping.Whereas previous work has been almost exclusively devoted to beam steering and null formation, this technique extends the versatility of coupled oscillator arrays through manipulation of the aperture amplitude and phase distribution.Furthermore, simultaneously an amplitude space taper distribution, (Δz / = 0), is attained to produce the desired far-field pattern.
The implementability of our method is based first on the actually monitoring of the time evolution of the oscillators' outputs (27).Afterward, the existence of the proper varactors and PINs by using (3) is addressed.It is obvious that the limited range of their values should set some constraints on the resultant pattern.So, in some cases, a proper iterative combination between the capacitance and conductance and the desired pattern is important.

International Journal of Antennas and Propagation
For all design cases, and independently on the complexity of the specific problem, the required excitation distribution can be readily implemented with proper tune of control parameters, such as the varactor capacitance and the PIN conductance of the oscillators.The relative closed-form formulas are presented for the case of the mutually synchronized coupled oscillator array.

Appendix
Following (25) and (26), the matrix form the initial pattern and its first variation can be written as where where [•] T stands for the transpose operation.Following (A.1) and (A.2), the new pattern is written as we have (A.5) In (A.5), it is

Figure 1 :
Figure 1: The equivalent circuit of the oscillator in the current-sum form.

Figure 4 :
Figure 4: The attained pattern for the general design case.

Figure 5 :
Figure 5: The coupled oscillator array response evolution with time for the general design case.(a) The evolution of the amplitudes.(b) The evolution of the phases.

Figure 6 :
Figure 6: The attained pattern for the constant phase design case.

Figure 7 :
Figure 7: The coupled oscillator array response evolution with time for the constant phase design case.The evolution of the amplitudes.

Figure 8 :
Figure 8: The attained pattern for the constant amplitude design case.

Table 6 :φ 7 − φ 6 =Figure 9 :
Figure 9: The coupled oscillator array response evolution with time for the constant amplitude design case.The evolution of the phases.

Figure 10 :
Figure 10: The attained amplitude sensitivity pattern for the general design case.

Table 1 :
Van derPol parameters of the employed oscillators.

Table 2 :
Details on the general design case.

Table 3 :
The tuning parameters for the general design case.

Table 4 :
Details on the constant phase design case.

Table 5 :
The tuning parameters for the constant phase design case.

Table 7 :
The tuning parameters for the constant amplitude design case.