Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase Optimization

*Mykhaylo Andriychuk*

## **Abstract**

In the process of solving the problems of generalized phase optimization the necessity to apply an eigenvalue approach often appears. The practical statement of the optimization problems consists of using the amplitude characteristics of functions that are sought. The usual way of optimization is deriving the Euler equation of the functional, which is used as criterion of optimization. As a rule, such equation is an integral one. It is worth pointing out that the integral equations of the generalized phase optimization are nonlinear ones. The characteristic property of such equations is non-uniqueness of solutions and their branching or bifurcation. The determination of branching solutions leads to the investigation of the corresponding homogeneous equations and the respective eigenvalue problem. This problem is nonlinear because of specificity of the statement of the optimization problem. The study of the above problem allows us to determine a set of points, in which the respective eigenvalues are equal to unity that determines the branching points of solutions. The data of calculations testify to the ability of the approach proposed to determine the solutions of nonlinear equations numerically with not large computations.

**Keywords:** nonlinear optimization, variational approach, radiation characteristic, nonlinear eigenvalue problem, bifurcation of solutions, computational modeling

## **1. Introduction**

The nonlinear eigenvalue approach is used in this chapter for the study of the properties to solutions of the generalized phase problem related to the synthesis of radiation systems through the incomplete data. Such incompleteness is considered here in the example of an indeterminate phase characteristic of function, which characterizes the radiation of the plane antenna arrays.

The problems with an indeterminate phase of the wave field arise in various applications and are widely described in the literature. The most well-known of these is the so-called phase problem (see, for example, [1–4]). It consists in restoring the phase distribution (argument) of the Fourier transform of a finite function by its amplitude (module) given (measured) along the entire real axis. This problem belongs to the classical problems of recovery (identification) and requires the conditions of existence of a unique solution.

In this chapter, another class of inverse problems is considered, and it can be termed optimization (design) problems. In sense of the Fourier transform, this can be, for example, the problem of finding such a finite complex function, the modulus of its Fourier transform satisfies a certain requirement (e.g., is close to a given positive function). As a rule, such requirements are formulated in the variational form, as the minimization of certain functionals. Obviously, such a formulation does not require a uniform solution. On the contrary, the existence of many solutions is often desirable because the above allows many degrees of freedom to determine an appropriate solution. The characteristic applications of such phase optimization problems include the theory of power transmission lines, field converters, antennas and resonators. The first works dealing with nonlinear inverse problems of such type appeared in the second part of the last century (see, for instance [5–10]).

In mathematical terms, problems of this type are reduced to the nonlinear integral equations of Hammerstein type [11–13]. They contain a linear kernel and a nonlinear multiplier that depends on a complex unknown function as an integrand. As a rule, the argument (phase) of this function appears there separately from the module. Similar equations are found in the literature in the context of the mentioned phase problem [14, 15]. They have different solutions, and the study of their structure and process of branching or bifurcation is an interesting mathematical problem [16].

Due to their nonlinearity, the problems under consideration require the development and application of special analytical and numerical methods for their solving. Along with the iterative methods that simulate the physical processes of field formation, the various modifications of Newton's method could be the most promising in this direction [17]. One such modification, which uses solving the nonlinear eigenvalue problems and searching for the zero curves of respective determinants, is proposed in this Chapter. It allows simultaneously with the finding of the branch of solutions to detect the presence of branching points on it and to determine them approximately, provided by this the initial approximations for more accurate calculation.

The nonlinear eigenvalue problems arise in pure and applied mathematics, as well as in the different areas of science that investigate the nonlinear phenomena [18, 19]. A variety of analytical-numerical methods have been elaborated till now for solving the nonlinear problems in acoustics, electrodynamics, fluid dynamics and other areas of applied science [20, 21]. The methods, developed until that time, were focused mainly on solving one-dimensional problems. The difficulties of analytical and computational nature appear if to apply them to a multidimensional problem. The method of implicit function is one of effective tools that been applied for solving the two- and three-dimensional nonlinear eigenvalue problem in the last two decades [22–24]. The extension of this method, which leads to solving the Cauchy problem (21) and (22), we apply in Section 3 to solve the nonlinear two-dimensional eigenvalue problem.

## **2. The operators of direct electrodynamics problem**

In the physical relation, the radiation system represents the plane array with the rectangular or hexagonal placement of radiators. Firstly, we consider the array with the rectangular ordering of separate elements (**Figure 1a**).

Consider a plane rectangular array consisting of *N*<sup>2</sup> � *M*<sup>2</sup> ¼ ð Þ� 2*N* þ 1 ð Þ 2*M* þ 1 identical elements (radiators), which are located in the *xOy* plane of the Cartesian coordinate system equidistantly for each of the coordinates. Since the radiators are identical, it is possible to formulate the synthesis problem not for the whole three*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

#### **Figure 1.**

*The placement of elements for the considered arrays. a) rectangular. b) hexagonal.*

dimensional vector directivity pattern (DP), but only for some complex scalar function *f x*ð Þ 1, *x*<sup>2</sup> that is termed as the array multiplier. This function for a rectangular equidistant array has the form [25]:

$$f(\mathbf{x}\_1, \mathbf{x}\_2) = A\mathbf{I} \equiv \sum\_{n=-N}^{N} \sum\_{m=-M}^{M} I\_{nm} e^{i(c\_1 m \mathbf{x}\_1 + c\_2 m \mathbf{x}\_2)} \tag{1}$$

where **I** ¼ f g *Inm*, �*N* ≤ *n*≤ *N*, �*M* ≤ *m* ≤ *M* is a set of excitations (currents) in the array's elements, *x*<sup>1</sup> ¼ sin *θ* cos *φ=* sin *α*1, *x*<sup>2</sup> ¼ sin *θ* sin *φ=* sin *α*<sup>2</sup> are the generalized angular coordinates, *c*<sup>1</sup> ¼ *kd*<sup>1</sup> sin *α*1, *c*<sup>2</sup> ¼ *kd*<sup>2</sup> sin *α*2, *k* ¼ 2*π=λ* is wave number, *d*<sup>1</sup> and *d*<sup>2</sup> are the distances between radiators along the *Ox* axis and *Oy* axis respectively, *α*<sup>1</sup> and *α*<sup>2</sup> are the angular coordinates, within which the desired power DP *P x*ð Þ 1, *x*<sup>2</sup> is not equal to zero (*P x*ð Þ� 1, *x*<sup>2</sup> 0 outside these angles). The function *f x*ð Þ 1, *x*<sup>2</sup> possesses 2*π=c*1�periodicity with respect to *x*<sup>1</sup> and 2*π=c*2� periodicity with respect to *x*2. Let us denote the region of change of coordinates *x*<sup>1</sup> and *x*<sup>2</sup> on one period as Ω ¼ f g ð Þ *x*1, *x*<sup>2</sup> :j*x*1j≤*π=c*1, j*x*2j≤*π=c*<sup>2</sup> . Below, the function *f x*ð Þ 1, *x*<sup>2</sup> is termed as the DP of array.

A similar formula can be derived for the array with the hexagonal placement of separate elements (**Figure 1b**)

$$f(\mathbf{x}\_1, \mathbf{x}\_2) = A\mathbf{I} \equiv \sum\_{m=-M\_2}^{M\_2} \sum\_{n=-N\_1(m)}^{N\_1(m)} I\_{nm} e^{i(c\_1 n x\_1 + c\_2 m x\_2)} \tag{2}$$

where *M* ¼ 2*M*<sup>2</sup> þ 1 is quantity of the linear subarrays, then *N* ¼ 2*N*1ð Þþ *m* 1 is the number of elements in the *m*�th subarray.

Eqs. (1) and (2) for DP *f x*ð Þ 1, *x*<sup>2</sup> represent the result of a linear operator *A*, which acts on a complex-valued space *HI* <sup>¼</sup> *CN*2�*M*<sup>2</sup> (rectangular case) or *HI* <sup>¼</sup> *<sup>C</sup><sup>N</sup>*0�*<sup>M</sup>* (hexagonal case) to the space of complex functions of two variables defined in the domain Ω. The value *N*<sup>0</sup> determines the number of elements in the central linear subarray in the hexagonal case.

Assume that the desired power DP *P s*ð Þ 1, *s*<sup>2</sup> is not equal to zero in some regions *G* ⊆ Ω, and it is equal to zero outside. The optimization problem is formulated as the minimization problem of the functional

$$\sigma\_a(\mathbf{I}) = \left\| P - \left| A\mathbf{I} \right|^2 \right\|\_f^2 + a \left\| \mathbf{I} \right\|\_I^2 \tag{3}$$

where k k� *<sup>f</sup>* and k k� *<sup>I</sup>* determine the norms in the space of DPs and space of currents respectively, which are defined by the inner products

$$\left\|\boldsymbol{f}\right\|\_{\boldsymbol{f}}^{2} = \left(\boldsymbol{f}\_{1}, \boldsymbol{f}\_{2}\right)\_{\boldsymbol{f}} = \left\|\boldsymbol{f}\_{1}(\boldsymbol{\omega}\_{1}, \boldsymbol{\omega}\_{2})\overline{\boldsymbol{f}}\_{2}(\boldsymbol{\omega}\_{1}, \boldsymbol{\omega}\_{2})d\boldsymbol{\omega}\_{1}d\boldsymbol{\omega}\_{2} \tag{4}$$

$$\left\|\mathbf{I}\right\|\_{I}^{2} = (\mathbf{I}\_{1}, \mathbf{I}\_{2})\_{I} = \frac{4\pi^{2}}{\mathfrak{c}\_{1}\mathfrak{c}\_{2}} \sum\_{n=-N}^{N} \sum\_{m=-M}^{M} I\_{1nm} \overline{I}\_{2nm} \tag{5}$$

here the values *f* <sup>2</sup>ð Þ *x*1, *x*<sup>2</sup> and *I*2*nm* are conjugated to *f* <sup>2</sup>ð Þ *x*1, *x*<sup>2</sup> and *I*2*nm*.

The nonlinear integral equation for the complex vector **I**of currents in space *H***I**, which is derived using the necessary condition of the minimum of functional (3), has the form [26].

$$a\mathbf{I} + \mathbf{2A}^\* \left( |\mathbf{A}\mathbf{I}| \mathbf{A} \mathbf{I} \right) - \mathbf{2A}^\* \left( \mathbf{P} \mathbf{A} \mathbf{I} \right) = \mathbf{0} \tag{6}$$

Here *<sup>A</sup>*<sup>∗</sup> is the operator adjoint to *<sup>A</sup>*, its form is defined by equality ð Þ *<sup>A</sup>***I**, *<sup>f</sup> <sup>f</sup>* <sup>¼</sup> **<sup>I</sup>**, *<sup>A</sup>*<sup>∗</sup> ð Þ*<sup>f</sup> <sup>I</sup>*. Using the inner products (4), (5) and Eq. (1) we obtain.

$$\begin{split} (A^\*f)\_{nm} &= \frac{c\_1 c\_2}{4\pi^2} \left[ \iint (\mathbf{x}\_1, \mathbf{x}\_2) e^{-i(c\_1 n \mathbf{x}\_1 + c\_2 n \mathbf{x}\_2)} d\mathbf{x}\_1 d\mathbf{x}\_2, n = -N, \ -N+1 \ldots, \\ (N-1, N, m = -M, \ -M+1, \ldots, M-1, M. \end{split} \tag{7}$$

If to act by operator *A* on both the parts of (6), we get a nonlinear integral equation of Hammerstein type for the function *f*

$$\text{agf} + 2\text{AA}^\* \left( |f| \mathcal{f} \right) - 2\text{AA}^\* \left( \mathcal{P} \mathcal{f} \right) = \mathbf{0} \tag{8}$$

The kernel of the *AA*<sup>∗</sup> operator for the rectangular array is defined as

$$K(c\_1, c\_2, \mathbf{x}\_1, \mathbf{x}\_1', \mathbf{x}\_2, \mathbf{x}\_2') = K\_1(c\_1, \mathbf{x}\_1, \mathbf{x}\_1') K\_2(c\_2, \mathbf{x}\_2, \mathbf{x}\_2'),\tag{9}$$

where

$$K\_1(\mathbf{x}\_1, \mathbf{x}\_1', c\_1) = \frac{c\_1}{\pi} \frac{\sin \left( N\_2 c\_1 (\mathbf{x}\_1 - \mathbf{x}\_1') / 2 \right)}{\sin \left( c\_1 (\mathbf{x}\_1 - \mathbf{x}\_1') / 2 \right)} \tag{10}$$

$$K\_2(\mathbf{x}\_2, \mathbf{x}\_2', c\_2) = \frac{c\_2}{\pi} \frac{\sin\left(\mathbf{M}\_2 c\_2 (\mathbf{x}\_2 - \mathbf{x}\_2')\right) / 2}{\sin\left(c\_2 (\mathbf{x}\_2 - \mathbf{x}\_2')\right) / 2} \tag{11}$$

The kernel of the *AA*<sup>∗</sup> operator for the hexagonal array is more complicated because we can not to present it in the form of two multipliers

$$K\left(\mathbf{c}\_{1},\mathbf{c}\_{2},\mathbf{x}\_{1},\mathbf{x}\_{1}^{\prime},\mathbf{x}\_{2},\mathbf{x}\_{12}^{\prime}\right) = \frac{\sin\left[c\_{1}(\mathbf{N}\_{1}(\mathbf{0}) - \mathbf{1}/2)\left(\mathbf{x}\_{1} - \mathbf{x}\_{1}^{\prime}\right)\right]}{\sin\left(\mathbf{1}/2c\_{1}(\mathbf{x}\_{1} - \mathbf{x}\_{1}^{\prime})\right)} + $$

$$-2\sum\_{m=1}^{M\_{1}}\cos m c\_{2}\left(\mathbf{x}\_{2} - \mathbf{x}\_{2}^{\prime}\right)\left\{\frac{\sin\left[c\_{1}(\mathbf{N}\_{1}(m) - \mathbf{1}/2)\left(\mathbf{x}\_{1} - \mathbf{x}\_{1}^{\prime}\right)\right]}{\sin\left(\mathbf{1}/2c\_{1}(\mathbf{x}\_{1} - \mathbf{x}\_{1}^{\prime})\right)},\ N\_{1}(m)\text{ is odd},\tag{12}$$

$$2\sum\_{n=1}^{N\_{1}(m)}\cos\left[c\_{1}(n - 1/2)\left(\mathbf{x}\_{1} - \mathbf{x}\_{1}^{\prime}\right)\right],\ N\_{1}(m)\text{ is even}.$$

*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

The kernels (9) and (12) of the integral Eq. (8) are real and degenerate. Since Eqs. (6) and (8) are nonlinear ones, both may have a non-unique solution. The number of solutions and their properties is studied according to the method proposed in [16, 27]. In the practical applications, the solution of Eqs. (6) and (8) is performed by the method of successive approximations. The convergence of the method depends on the parameter *α*, desired DP *P x*ð Þ 1, *x*<sup>2</sup> , as well as the parameters *c*1and *c*2contained in the kernels (9) and (12).

## **3. Search for the bifurcation curves**

We should use the linear integral equation to define the bifurcation curves according to [16]. Based on this equation, we pass to the respective eigenvalue problems, solutions of which allow us to find the characteristic values of parameters *c*<sup>1</sup> and *c*<sup>2</sup> in the kernel of equation, at which the bifurcation appears.

#### **3.1 Description of procedure**

The linear equation

$$q\mathcal{f} = \mathcal{Z}AA^\*\left(\mathcal{P}\mathcal{f}\right) \tag{13}$$

is used to study the properties of Eq. (8).

In contrast to a similar equation for the amplitude DP synthesis problem [25], Eq. (8) does not have a trivial nonzero initial solution *f* <sup>0</sup> for all parameters *c*<sup>1</sup> and *c*2; the trivial solution *f* <sup>0</sup> is zero for it, so in contrast to the problem of synthesis by amplitude DP, we are not talking about the branching of solutions, but about their bifurcation.

The problem of finding bifurcation curves is reduced to the corresponding eigenvalue problem. The equation for eigenfunctions and corresponding eigenvalues, which refers to (13), is

$$\mathbf{g}(\mathbf{x}\_1, \mathbf{x}\_2) = 2\lambda a^{-1} \iint\_{\Omega} \mathbf{g}(\mathbf{x}\_1', \mathbf{x}\_2') K\_1(c\_1, \mathbf{x}\_1, \mathbf{x}\_1') K\_2(c\_2, \mathbf{x}\_2, \mathbf{x}\_2') d\mathbf{x}\_1' d\mathbf{x}\_2' \tag{14}$$

As stated by the branching theory of solutions of the nonlinear equations [16], the bifurcation points can be those values of *c*1and *c*<sup>2</sup> at which Eq. (14) has nonzero solutions.

Using the properties of the degeneracy of the kernel *AA*<sup>∗</sup> , we reduce Eq. (14) to the equivalent system of the linear algebraic equations (SLAE). The coefficients of matrix of this equation depend on the parameters *c*<sup>1</sup> and *c*<sup>2</sup> analytically. To this end, the equations for eigenfunctions corresponding to (13) are written as

$$\mathbf{g}(\mathbf{x}\_1, \mathbf{x}\_2) = \sum\_{n=-N}^{N} \sum\_{m=-M}^{M} \mathbf{x}\_{nm} e^{i(c\_1 n x\_1 + c\_2 m x\_2)} \tag{15}$$

where

$$\propto \chi\_{nm} = \frac{c\_1 c\_2}{4\pi^2} \iint\limits\_{\Omega} (\mathbf{x}\_1', \mathbf{x}\_2') \mathbf{g}\left(\mathbf{x}\_1', \mathbf{x}\_2'\right) e^{-i\left(c\_1 \mathbf{x}\_1' + c\_2 \mathbf{x} \mathbf{x}\_2'\right)} d\mathbf{x}\_1' d\mathbf{x}\_2' \tag{16}$$

Multiplying both the parts of (15) on *P x*<sup>0</sup> 1, *x*<sup>0</sup> 2 � �*e* �*i c*1*kx*<sup>0</sup> <sup>1</sup>þ*c*2*lx*<sup>0</sup> ð Þ<sup>2</sup> at *<sup>k</sup>* ¼ �*N*, � *<sup>N</sup>* <sup>þ</sup> 1, … , *N* � 1, *N*, *l* ¼ �*M*, � *M* þ 1, … , *M* � 1, *M* and integrating over the domain Ω, we obtain a system of linear algebraic equations to determine the quantities *xnm*

$$\begin{aligned} \boldsymbol{x}\_{kl} &= \sum\_{n=-N}^{N} \sum\_{m=-M}^{M} a\_{nm}^{(kl)}(c\_1, c\_2) \boldsymbol{x}\_{nm}, \ \boldsymbol{k} = -\boldsymbol{N}, \ -\boldsymbol{N} + \mathbf{1}, \ \dots, \boldsymbol{N} - \mathbf{1}, \boldsymbol{N}, \\ \boldsymbol{l} &= -\boldsymbol{M}, \ -\boldsymbol{M} + \mathbf{1}, \ \dots, \boldsymbol{M} - \mathbf{1}, \boldsymbol{M}, \end{aligned} \tag{17}$$

where

$$a\_{nm}^{(kl)} = \frac{c\_1 c\_2}{4\pi^2} \left[ \int\_{\Omega} (\varkappa\_1, \varkappa\_2) e^{-i[(c\_1(n-k)\mathbf{x}\_1 + c\_2(m-l)\mathbf{x}\_2)]} d\mathbf{x}\_1 d\mathbf{x}\_2 \right] \tag{18}$$

and matrix of the coefficients *a* ð Þ *kl nm* is self-adjoint and Hermitian.

Thus, we obtained a two-parameter nonlinear spectral problem corresponding to a homogeneous SLAE (17). This problem can be given as

$$(E\_M - A\_M(c\_1, c\_2))\mathbf{x} = \mathbf{0} \tag{19}$$

where *AM* is the matrix of coefficients *a* ð Þ *kl nm* , *EM* is a unit matrix of dimension *N*<sup>2</sup> � *M*2. For the system (19), the equality

$$\Psi(c\_1, c\_2) = \det[E\_M - A\_M(c\_1, c\_2)] = \mathbf{0} \tag{20}$$

must be met to have a non-zero solution.

One can easy to make sure that the function Ψð Þ *c*1,*c*<sup>2</sup> is real. Moreover, since *AM*ð Þ *c*1,*c*<sup>2</sup> is the Hermitian matrix, then *EM* � *AM*ð Þ *c*1,*c*<sup>2</sup> is Hermitian too. The determinant of the Hermitian matrix is a real number [28]. Thus, Ψð Þ *c*1,*c*<sup>2</sup> is a real function of real arguments *c*1and *c*2.

Consequently, the problem to find the eigenvalues of Eq. (14) or to determine the solution of the equivalent SLAE (19) is reduced to finding zeros of function Ψð Þ *c*1,*c*<sup>2</sup> .

If to consider the equation Ψð Þ¼ *c*1,*c*<sup>2</sup> 0 as a problem of determining an implicit function *c*<sup>2</sup> ¼ *c*2ð Þ *c*<sup>1</sup> in the vicinity of some point *c*1, we get Cauchy problem [29].

$$\frac{dc\_2}{dc\_1} = \frac{\Psi\_{c\_1}(c\_1, c\_2)}{\Psi\_{c\_2}(c\_1, c\_2)}\tag{21}$$

$$c\_2 \left( c\_1^{(0)} \right) = c\_2^{(0)} \tag{22}$$

To retrieve the initial conditions (22) we pass to an auxiliary one-dimensional nonlinear spectral problem if to substitute *c*<sup>2</sup> by *c*<sup>2</sup> ¼ *γc*<sup>1</sup> in Eq. (20) with some real parameter *γ*. As a result, we get the one-dimensional eigenvalue problem

$$(E\_M - A\_M(c\_1, \gamma c\_1))\ddot{\mathbf{x}} \equiv \left(E\_M - \tilde{A}\_M(c\_1)\right)\ddot{\mathbf{x}} = \mathbf{0} \tag{23}$$

Eq. (20), which corresponds to Eq. (23), is

$$\Psi(\mathcal{c}\_1, \mathcal{yc}\_1) = \det\left[E\_M - \tilde{A}\_M(\mathcal{c}\_1)\right] = \mathbf{0} \tag{24}$$

*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

Let *c* ð Þ 0 <sup>1</sup> be the solution of the Eq. (24), then *c* ð Þ 0 <sup>1</sup> ,*c* ð Þ 0 2 � � <sup>¼</sup> *<sup>c</sup>* ð Þ 0 <sup>1</sup> , *γc* ð Þ 0 1 � � is the point that corresponds to eigenvalue *λ*<sup>0</sup> ≈1 of Eq. (15). By solving Eqs. (21) and (22) in a small vicinity of point *c* ð Þ 0 <sup>1</sup> ,*c* ð Þ 0 2 � �, we find the spectral curve of the matrix-function *AM*ð Þ *c*1,*c*<sup>2</sup> , which is the curve *c*2ð Þ *c*<sup>1</sup> defining a set of the bifurcation points.

The eigenfunctions of Eq. (14) are defined as the eigenvectors of matrix *AM*ð Þ *c*1,*c*<sup>2</sup> using the resulting solution of the Cauchy problem with the sought solutions Ψð Þ *c*1,*c*<sup>2</sup> . In this procedure, a four-dimensional matrix *AM*ð Þ *c*1,*c*<sup>2</sup> is reduced to a twodimensional one by the relevant renouncement of its elements.

#### **3.2 Defining the area of nonzero solutions**

Due to the peculiarity of the problem statement according to desired power DP *P x*ð Þ 1, *x*<sup>2</sup> , Eq. (8) has zero solution at arbitrary values of the parameters *c*1,*c*2, *α*. From an engineering point of view, this is a significant drawback, but for some desired DPs *P x*ð Þ 1, *x*<sup>2</sup> it is possible to fix an area of parameters *c*1,*c*2, *α* at which a nonzero solution exists. At the small *c*<sup>1</sup> and *c*2, the kernel (9) is given approximately in the form

$$K(c\_1, c\_2, \mathbf{x}\_1, \mathbf{x}\_1', \mathbf{x}\_2, \mathbf{x}\_2') \approx \frac{M\_2 N\_{2^2 1} c\_2}{\pi^2} \tag{25}$$

Assuming that *f x*ð Þ 1, *x*<sup>2</sup> is constant, the integral Eq. (8) can be rewritten as (usually for small *c*<sup>1</sup> and *c*<sup>2</sup> *f x*ð Þ 1, *x*<sup>2</sup> ≈ const).

$$\frac{\pi^2 a}{2M\_2 N\_2 c\_1 c\_2} = \int\_{-1}^{1} \int\_{-1}^{1} P(\mathbf{x}\_1, \mathbf{x}\_2) d\mathbf{x}\_1 d\mathbf{x}\_2 - 4 \left| f(\mathbf{x}\_1, \mathbf{x}\_2) \right|^2 \tag{26}$$

The area of integration Ω in Eq. (8) is reduced in the last formula to the area ½ �� � �1, 1 ½ � 1, 1 because of definition of both the arguments *x*1, *x*<sup>2</sup> and parameters *c*1,*c*2.

Taking into account that j j *f x*ð Þ 1, *x*<sup>2</sup> <sup>2</sup> is positive, we get the following relationship between the function *P x*ð Þ 1, *x*<sup>2</sup> and the parameters *c*1, *c*2, and *α*:

$$\int\_{-1}^{1} \int\_{-1}^{1} P(\mathbf{x}\_1, \mathbf{x}\_2) d\mathbf{x}\_1 d\mathbf{x}\_2 - \frac{\pi^2 a}{2M\_2 N\_2 c\_1 c\_2} > 0 \tag{27}$$

Finally, considering the case *P x*ð Þ� 1, *x*<sup>2</sup> 1, we obtain:

$$c\_1 c\_2 > \frac{\pi^2 a}{8M\_2 N\_2} \tag{28}$$

In fact; inequality (28) determines the area of parameters *c*1,*c*2, *α*, where nonzero solutions exist. In **Figure 2**, the dependence curves *c*<sup>2</sup> ¼ *c*2ð Þ *c*<sup>1</sup> for three different values *M*<sup>2</sup> and *N*<sup>2</sup> are shown. The results are given for array with the number of elements *M*<sup>2</sup> ¼ *N*<sup>2</sup> ¼ 3 (curve 1), *M*<sup>2</sup> ¼ *N*<sup>2</sup> ¼ 5 (curve 2) and *M*<sup>2</sup> ¼ *N*<sup>2</sup> ¼ 11 (curve 3). The area of values *c*<sup>1</sup> and *c*2, where the existence of zero solutions is possible, according to the estimate (28) is located below and to the left of the presented curves.

**Figure 2.** *The curves c*<sup>2</sup> ¼ *c*2ð Þ *c*<sup>1</sup> *at the different M*<sup>2</sup> *and N*2*.*

As can be seen, the area of zero values decreases significantly with increasing *N*2and *M*2. The obtained results testify that the zero solutions of Eq. (8) for a given constant power DP can exist either at a small value *c*1*c*<sup>2</sup> corresponding to low frequencies (at a given size of array), or at the values of *c*<sup>1</sup> that significantly exceeding *c*<sup>2</sup> and vice versa. The last case corresponds to arrays with a large difference in distance of elements along the coordinate axes. Such arrays are usually rarely used in practice.

### **3.3 Determination of bifurcation lines**

### *3.3.1 The case of rectangular array*

The finding of bifurcation lines of the nonlinear Eq. (8) was performed for the array containing *N*<sup>2</sup> ∗ *M*<sup>2</sup> ¼ 11 ∗ 11 ¼ 121 radiators for the desired power DP *P x*ð Þ¼ 1, *x*<sup>2</sup> 1 at Λ*<sup>c</sup>* ¼ f g ð Þ *c*1,*c*<sup>2</sup> , 0< *c*1,*c*<sup>2</sup> ≤2 for the different values of the parameter *α* in (3).

The search for bifurcation lines can be performed directly by investigating the properties of the determinant (20) as a function of the parameters *c*<sup>1</sup> and *c*2. In addition, the function (20) depends on the parameter *α*; so the set of its eigenvalues also depends on this parameter, i.e. the set of spectral curves that separate the areas of zero and nonzero solutions.

The behavior of the corresponding curves when changing the parameter *α* is shown in **Figure 3**. The behavior of the determinant (24) depending on the parameters *c*1and *c*<sup>2</sup> at *α* ¼ 0*:*5 is given in **Figure 3a**; and in **Figure 3b**–**d**, the intersection of this function with a plane Ψð Þ¼ *c*1, *γc*<sup>1</sup> 0 is illustrated at the different *α*. This results in a set of curves that correspond to a set of spectral lines separating the area of zero and nonzero solutions. At a fixed size of array, the area where zero solutions can exist expands if the parameter *α* increases, this area is located below the left of the first curve.

The curves marked by number 1 correspond to the solutions with constant (zero or even) phase DP; curves with number 2 correspond to the solutions with phase DP that *Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

**Figure 3.** *The spectral curves of Eq. (19) at the different α.*

is even with respect to the *Ox*<sup>1</sup> axis, and odd with respect to the *Ox*<sup>2</sup> axis, and curves numbered by 3 correspond to the solutions with a phase DP odd with respect to two coordinate axes. The proposed procedure is quite approximate, it does not allow to separate the curves that correspond to different types of solutions and thus identify the areas where there is a nonzero solution for the synthesized power DP with the specified phase property.

The method of implicit function proposed in [23] and developed for plane array in [30] is devoid of this drawback.

At the first step of this method, a series of one-dimensional eigenvalue problems is solved, by this the different values of parameter *γ* are prescribed by the relation *c*<sup>2</sup> ¼ *γc*1and a one-dimensional problem is solved with respect to *c*1. In **Figures 4** and **5**, the first four eigenvaluesof the problem at *γ* ¼ 1*:*0 and *γ* ¼ 0*:*2 are shown. The values *c* ð Þ*i* <sup>1</sup> ,*c* ð Þ*i* <sup>2</sup> ¼ *c* ð Þ*i* 1 , *<sup>i</sup>* <sup>¼</sup> 1, 2, 3, 4, at which *<sup>λ</sup><sup>i</sup>* <sup>¼</sup> 1, are the bifurcation points in the plane *<sup>c</sup>*1*Oc*2. By this, the set of points ð Þ *<sup>c</sup>*1,*c*<sup>2</sup> at which the eigenvalue *<sup>λ</sup>*ð Þ*<sup>i</sup>* <sup>¼</sup> 1 is determined approximately from the graphical data.

The next step is to refine the values *c* ð Þ*i* <sup>1</sup> ,*c* ð Þ*i* 2 by solving the transcendental Eq. (20), and the point *c* ð Þ*i* <sup>1</sup> ,*c* ð Þ*i* 2 , which is considered as the initial approximation.

**Figure 4.** *The first eigenvalues at the ray c*<sup>1</sup> ¼ *c*2*, α* ¼ 0*:*5*.*

**Figure 5.** *The first eigenvalues at the ray c*<sup>2</sup> ¼ 0*:*2*c*1*, α* ¼ 0*:*5*.*

In the final step, the bifurcation curve in the plane ð Þ *c*1,*c*<sup>2</sup> is determined by solving Eqs. (21) and (22), after specification of the values *c* ð Þ*i* <sup>1</sup> ,*c* ð Þ*i* 2 . In **Figure 6**, the bifurcation curves *c* ð Þ1 <sup>1</sup> ,*c* ð Þ1 2 – *<sup>c</sup>* ð Þ 4 <sup>1</sup> ,*c* ð Þ 4 2 that correspond to the first four eigenvalues are shown. The curve with number 1 corresponds to the solution with the zero (even) phase of the created DP. This curve corresponds to that is marked by 1 in **Figure 3b**.

There are no nonzero solutions with such a phase property for the values *c*<sup>1</sup> and *c*<sup>2</sup> above and to the right of this curve. Curves 2 and 2<sup>0</sup> correspond to solutions in which the phase DP is symmetric about one axis and asymmetric about the other axis (obviously, for a plane array, there are two such curves and they are antisymmetrical). Curve 2 corresponds that is marked by 2 in **Figure 3b**. The curve with number 3

*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

**Figure 6.** *The bifurcation curves corresponding to set c*ð Þ*<sup>i</sup>* <sup>1</sup> ,*c* ð Þ*i* 2 , *<sup>i</sup>* <sup>¼</sup> 1, … , 4*, N*<sup>2</sup> <sup>¼</sup> *<sup>M</sup>*<sup>2</sup> <sup>¼</sup> <sup>11</sup>*.*

corresponds to a solution with a phase DP antisymmetrical (odd) with respect to both the axes. The location of the areas of zero and non-zero solutions is the same as in **Figure 3b**. It should be noted that the problem of refining the roots of Eq. (20) is the most time-consuming in computational relation because refining the roots of this equation requires a series of computational experiments with different values *c* ð Þ*i* 1,0,*c* ð Þ*i* 2,0 of initial parameters close to approximate values.

### *3.3.2 The case of hexagonal array*

Firstly, we consider the procedure of determination of bifurcation curves by finding zero lines of determinant (24). The results, similar to those are presented in **Figure 3a** and **b** for the rectangular array, are shown in **Figures 7** and **8**. One can see that the behavior of function Ψð Þ *c*1,*c*<sup>2</sup> is more complex than in the case of rectangular array. The obtained graphs testify that the solutions with other different behavior of phase argð Þ *f x*ð Þ 1, *x*<sup>2</sup> of the DP appear additionally. One such solution is marked by number 4. Other solutions appear when parameters *c*<sup>1</sup> and *c*<sup>2</sup> increase at the fixed *α*.

Search of the bifurcation curves is carried out similarly to the case of rectangular array. The numerical results are presented for the array with *Ntot* ¼ 61 elements for the desired power DP *N*0ð Þ¼ *x*1, *x*<sup>2</sup> 1 at Λ*<sup>c</sup>* ¼ f g ð Þ *c*1,*c*<sup>2</sup> , 0 <*c*1,*c*<sup>2</sup> ≤2*:*0 for the different values of *α* in the functional (3). At the first step, the one-dimensional eigenvalue problems were solved at the different values of parameter *γ*. In **Figure 9**, the first four eigenvalues are shown at *γ* ¼ 1*:*0, and in **Figure 10**, they are shown at *γ* ¼ 0*:*2. Similar to the case of rectangular array, the points, in which *λ<sup>i</sup>* ¼ 1 are moved to right and the distance between them increases at *γ* ¼ 0*:*2. The values *c* ð Þ*i* <sup>1</sup> ,*c* ð Þ*i* <sup>2</sup> ¼ *γc* ð Þ*i* 1 , where

**Figure 7.** *The surface of determinant (24) values for the hexagonal array, α* ¼ 0*:*5*.*

**Figure 8.** *Zero lines of determinant (24) for the hexagonal array, α* ¼ 0*:*5*.*

*i* ¼ 1, 2, 3, 4, are the bifurcation points in the plane ð Þ *c*1,*c*<sup>2</sup> . The points *c* ð Þ*i* <sup>1</sup> ,*c* ð Þ*i* 2 , for which the eigenvalues *<sup>λ</sup>*ð Þ*<sup>i</sup>* <sup>¼</sup> 1 are determined approximately in this step.

The specification of values *c* ð Þ*i* <sup>1</sup> ,*c* ð Þ*i* 2 by solving Eq. (20) is carried out in the next step, and the points *c* ð Þ*i* <sup>1</sup> ,*c* ð Þ*i* 2 of the graph data from the **Figures 9** and **<sup>10</sup>** are used as initial approximations. The usual numerical half-division method is used for this goal.

The bifurcation points *c* ð Þ*i* <sup>1</sup> ,*c* ð Þ*i* 2 , *<sup>i</sup>* <sup>¼</sup> 1, 2, 3, 4 for the first four eigenvalues in the rays, *c*<sup>2</sup> ¼ *γc*<sup>1</sup> are shown in **Figure 11**. The respective curves of bifurcations, which are obtained by solving Eqs. (21) and (22), are shown in **Figure 12**. As in the case of a rectangular array, to obtain the necessary data, we should carry out precise computations.

*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

**Figure 9.**

*The first eigenvalues for the ray c*<sup>1</sup> ¼ *c*<sup>2</sup> *at α* ¼ 0*:*5*.*

**Figure 10.** *The first eigenvalues for the ray c*<sup>2</sup> ¼ 0*:*2*c*1*, α* ¼ 0*:*5*.*

## **4. The engineering applications**

The results presented in this Section demonstrate how the knowledge about the point of bifurcation obtained as the solutions of the nonlinear eigenvalue problems allows us to understand better the process of bifurcation and how to get the solutions, which are the most optimal in sense of the used criterion of optimization.

### **4.1 The method of successive approximations**

The properties of solutions to Eq. (8) obtained by using the method of successive approximation are related directly with the properties of phase characteristic of the

**Figure 11.** *The bifurcation points at the rays c*<sup>2</sup> ¼ *γc*1*.*

eigenfunctions, which are determined at solving the eigenvalue problem. Prescribing the initial approximation *f* <sup>0</sup> for the iterative process for solving Eq. (8) with the specified property of the phase arg *f* <sup>0</sup>, we could receive the solution of Eq. (8) with the same phase property in the wide range of characteristic parameters *c*<sup>1</sup> and *c*2. This is important for the engineering design of arrays having the fixed phase characteristics of radiation in the defined range of frequencies.

The method of successive approximations

$$\left(f\_{n+1} - \beta f\_n + (\mathbf{1} - \beta)\mathcal{B}\left(f\_n\right) = \mathbf{0}, n = \mathbf{0}, \mathbf{1}, \mathbf{2}, \dots\right) \tag{29}$$

is used for solving Eq. (8) with a set of specific physical parameters of array. In the last formula,

*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

$$B(f) = \frac{2}{a} \left[ AA^\* \left[ (P \cdot f) - \left( |f|^2 \cdot f \right) \right] \right] \tag{30}$$

Parameter *β* ∈½ � 0, 1 in (29) is used to accelerate the convergence of iterative process. To substantiate the condition of convergence of the iterative process (29), we apply Theorem 2.6.2 [26] (p. 133), which states that the operator *AA*<sup>∗</sup> be contraction one. This requirement is met when the inequality

$$a > 2||AA^\* \left[ f\left( N\_0 - \left| f \right|^2 \right) \right]|| \tag{31}$$

met. The results of numerical calculations show that condition (31) is overestimated and for some values of the problem parameters the iterative process (29) converges for values *α*that do not satisfy the estimate (31).

#### **4.2 The case of rectangular array**

In **Figure 13**, the dependence of the convergence of the iterative process (29) on the value of parameter *α* for a desired power DP *P x*ð Þ¼ 1, *x*<sup>2</sup> 1 at the fixed values *β* ¼ 0*:*1, *c*<sup>1</sup> ¼ *c*<sup>2</sup> ¼ 2*:*0, the number of radiators *M*<sup>2</sup> � *N*<sup>2</sup> ¼ 11 � 11 ¼ 121 is shown. The required accuracy *<sup>ε</sup>* <sup>¼</sup> <sup>10</sup>�<sup>3</sup> .

The results of solving the optimizing problem for this desired power DP at *α* ¼ 0*:*5 are shown in **Figure 14**. The approximation quality to a desired DP *P* significantly depends on the parameters *c*<sup>1</sup> and *c*<sup>2</sup> at both fixed *N*<sup>2</sup> and *M*2. The mean-square deviation (MSD) (the first term in (3)) is equal to 0.0847 for *c*<sup>1</sup> ¼ *c*<sup>2</sup> ¼ 1*:*0, and it is equal to 0.0075 for *c*<sup>1</sup> ¼ *c*<sup>2</sup> ¼ 3*:*14.

The synthesized DP ∣*f*∣for larger *c*1,*c*<sup>2</sup> has not only a more optimal mean-square approximation, but it is also closer to the shape of the desired DP *P*. The optimal amplitudes ∣*Inm*∣ of currents in the array's elements are close to constant at such parameters *c*<sup>1</sup> and *c*2.

When solving the optimizing problem for desired power DP of a more complex form, the quality of the approximation significantly depends on both the parameter *α*

#### **Figure 13.**

*The character of convergence of iterative process (29) at the different α, M*<sup>2</sup> ¼ *N*<sup>2</sup> ¼ 11*.*

**Figure 14.** *The created power RP* ∣*f*∣ *(a) and optimal distribution of currents* ∣*Inm*∣ *(b) at c*<sup>1</sup> ¼ *c*<sup>2</sup> ¼ 1*:*0*.*

and the type of initial approximation for the phase of a given DP. The results are shown for the desired power DP

$$P(\varkappa\_1, \varkappa\_2) = |\sin\left(\pi\varkappa\_1\right) \cdot \sin\left(\pi\varkappa\_2\right)|, \ -1 \le \varkappa\_1 \le 1, \ -1 \le \varkappa\_2 \le 1\tag{32}$$

at *c*<sup>1</sup> ¼ *c*<sup>2</sup> ¼ 3*:*14 and *α* ¼ 0*:*2 in **Figure 15**. Despite the fact that the shape of desired DP *P x*ð Þ 1, *x*<sup>2</sup> is more complex that in the previous example, decrease of *α* (from 0.5 to 0.2) and simultaneous increase of *c*<sup>1</sup> and *c*<sup>2</sup> (from 1.0 to 3.14) allows us to get the amplitude ∣*f x*ð Þ 1, *x*<sup>2</sup> ∣ of created DP, which is very close to the *P x*ð Þ 1, *x*<sup>2</sup> . The optimal distribution of currents' amplitudes ∣*Inm*∣ (**Figure 15b**) approaches the shape of the created DP.

We have used an additional optimization parameter *β* in Eq. (29), which, as shown by the results of numerical calculations, accelerates the convergence of iterative process significantly. In **Figure 16**, the results of the study of the influence of this parameter on the rate of convergence at a fixed value of the parameter *α* ¼ 0*:*5 are shown. The results are given for the desired power DP *P x*ð Þ¼ 1, *x*<sup>2</sup> 1 at *c*<sup>1</sup> ¼ *c*<sup>2</sup> ¼ 2*:*0. In order to achieve the accuracy 10�<sup>3</sup> of calculations, one needs 157 iterations at *<sup>β</sup>* <sup>¼</sup> 0*:*01. If parameter *β* increases to a certain value, the number of iterations decreases significantly, so at *β* ¼ 0*:*05, *β* ¼ 0*:*10, *β* ¼ 0*:*15, *β* ¼ 0*:*20, and *β* ¼ 0*:*25 one requires

**Figure 15.** *The created power RP f*j j<sup>2</sup> *at c*<sup>1</sup> <sup>¼</sup> *<sup>c</sup>*<sup>2</sup> <sup>¼</sup> <sup>3</sup>*:*14*.*

*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

**Figure 16.** *The convergence of iterative process (29) for the different β.*

62 iterations, 37 iterations, 28 iterations, 20 iterations, and 16 iterations, respectively. At the subsequent increase, the number of required iterations begins to increase and already at *β* ¼ 0*:*30 the iterative process begins to diverge. Numerical calculations show that the limit valueof *β*, at which the iterative process (29) begins to diverge, significantly depends on the value of the parameter *α*. So, if this parameter decreases, the threshold value of *β* increases. The dependence of the convergence on the array's parameters (*c*1,*c*2, *M*2, *N*2, *d*1, *d*2) is not so significant.

The values of the functional (3) for the created DP with different phases are shown in **Figure 17**. The solid curve corresponds to the phase DP even with respect to two axes, the dotted curve corresponds to phase DP even with respect to one axis and odd with respect to the other, the dashed curve corresponds to the phase DP odd with respect to both the axes.

**Figure 17.** *The values of functional (3) versus the phase of created DP.*

One can see that the values of functional at the fixed *c* (frequency) significantly depend on the phase of the created DP. The value *σ* ¼ 0*:*2 is achieved for the "eveneven" solution at *c* ¼ 0*:*79, for the "even-odd" solution at *c* ¼ 0*:*71 and for the "oddodd" solution at *c* ¼ 0*:*625. That is, within the used criterion, the latter type of solution is 21% better than the first one. From this fact, it follows that at a fixed distance between the radiators for the desired DP *P x*ð Þ¼ 1, *x*<sup>2</sup> 1, the number of array's elements can be reduced by 21% with the same value of MSD. A similar situation is observed for the characteristics of DP at *σ* ¼ 0*:*1, i.e. "odd-odd" solution is better on 19.4% than "even-even".

#### **4.3 The case of hexagonal array**

The results of solution of the optimization problem for two given power DPs *P*1ð Þ� *x*1, *x*<sup>2</sup> 1 and

$$P\_2(\mathbf{x}\_1, \mathbf{x}\_2) = \begin{cases} 2\sqrt{\mathbf{x}\_1^2 + \mathbf{x}\_2^2}\sqrt{1 - \mathbf{x}\_1^2 - \mathbf{x}\_2^2}, \mathbf{x}\_1^2 + \mathbf{x}\_2^2 \le \mathbf{1}, \\ 0, \mathbf{x}\_1^2 + \mathbf{x}\_2^2 > \mathbf{1}, \end{cases} \tag{33}$$

in the form of body of rotation are shown in **Figures 18** and **19** at *α* ¼ 0*:*5.

As previously, the optimization problem consists of solving Eq. (8) by the method of successive approximation (29). The MSD (the value of the first term in (3)) for the first desired DP is equal to 0*:*3774, and it is equal to 0*:*2218 for the second desired DP.

Similar to the case of rectangular array, the approximation quality to the desired DP *P* depends on both the parameters *c*1, *c*2, and *α*. The characteristic of MSD of DPs for *α* at the different *c*<sup>1</sup> on the ray *c*<sup>2</sup> ¼ 1*:*118*c*<sup>1</sup> is shown in **Figures 20** and **21**. The chosen relation between *c*<sup>1</sup> and *c*<sup>2</sup> provides the regularity of the array's geometry, and as the numerical computations have shown, gives the ability to get the close characteristics of radiation in the planes *x*1and *x*2.

The largest MSD for the *P*<sup>1</sup> is achieved at *α* ¼ 1*:*0 for *c*<sup>1</sup> ¼ 0*:*5, and it is equal to 1*:*96443, it diminishes almost linearly if parameter *α* decreases. The largest MSD for DP *P*<sup>2</sup> is equal to 1*:*43685. One should note that the value of MSD diminishes if *α*

**Figure 18.** *The amplitude of created DP f*j j<sup>2</sup> *for P*<sup>1</sup> *at c*<sup>1</sup> <sup>¼</sup> <sup>2</sup>*:*0,*c*<sup>2</sup> <sup>¼</sup> <sup>2</sup>*:*236*.*

*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

#### **Figure 19.**

*The amplitude f*j j<sup>2</sup> *of created DP for P*<sup>2</sup> *at c*<sup>1</sup> <sup>¼</sup> <sup>2</sup>*:*0,*c*<sup>2</sup> <sup>¼</sup> <sup>2</sup>*:*236*.*

### **Figure 20.**

*The MSD versus the parameter α for the DP P*1*.*

decreases, but the norm k k**I** *<sup>H</sup>***<sup>I</sup>** of current growth that is unacceptable from the engineering point of view.

The approximation quality to the desired DP depends also on the type of initial data, which are prescribed for the iterative process (29). The dependence of the values of functional (3) on the parity of phase of the initial approximation *f* <sup>0</sup> for the DP *P*<sup>1</sup> is shown in **Figure 22**. The results are shown for four types of initial approximation: odd with respect to both the axes (dotted curve), even with respect to the *Ox*<sup>1</sup> axis and odd with respect to the *Ox*<sup>2</sup> axis (dashed curve), odd with respect to the *Ox*<sup>1</sup> axis and even with respect to the *Ox*<sup>2</sup> axis (dash-dot curve), and even with respect to both the axis (solid curve). The initial approximation *f* <sup>0</sup>, corresponding to the even phase with respect to both the axis, is optimal for this DP, moreover the values of *σα* for the small values of parameters *c*<sup>1</sup> and *c*<sup>2</sup> differ significantly, but starting from *c*<sup>1</sup> ¼ 0*:*8 this difference does not exceed 10%.

**Figure 21.** *The MSD versus parameter α for DP P*2*.*

**Figure 22.**

*The values of σα versus the initial approximation of initial approximation for the iterative process (29).*

The dependence of convergence of the iterative process (29) on the parameter *α* at the fixed *β* ¼ 0*:*1 is shown in **Figure 23**. As in the case of a rectangular array, the iterative process converges most slowly at *α* ¼ 1*:*5, one needs 67 iterations to achieve the accuracy that is equal to 10�3. The number of iterations decreases if *α* diminishes. For example, one needs 30 iterations to achieve the same accuracy at *α* ¼ 0*:*2.

The dependence of convergence on the parameter *β* is studied too. The necessary number of iterations that needs to achieve the accuracy 10�<sup>3</sup> at *β* ¼ 0*:*01, *β* ¼ 0*:*05, *β* ¼ 0*:*10, *β* ¼ 0*:*15, *β* ¼ 0*:*20, and *β* ¼ 0*:*30 (curves 1–6 respectively) is shown in **Figure 24**. It is substantiated that the iterative process converges most slowly at *β* ¼ 0*:*01, it is necessary 152 iterations to achieve the prescribed accuracy. The most optimal among the considered *β* is *β* ¼ 0*:*20, one needs 20 iterations only to achieve

*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

**Figure 23.**

*The convergence of iterative process (29) versus α.*

**Figure 24.** *The convergence of iterative process (29) versus number of iteration, α* ¼ 0*:*5*.*

the above accuracy. The iterative process becomes slow at the subsequent growth of *β*. For example, one needs 37 iterations to achieve this accuracy; the iterative process (29) becomes convergent at *β* >0*:*40. This testifies that in the process of computations one should to limit by non-large values of *β* (*β* ≤0*:*20) that guarantees the convergence and considerably grows its speed on the contrast with small *β* (*β* ≤0*:*01).

More information about the problem under investigation one can find in [30–34].

## **5. Conclusions**

The problem of finding the solutions to the nonlinear integral equations and their properties is reduced to nonlinear two-dimensional eigenvalue problems that lead to

the subsequent application of an implicit function method for solving the Cauchy problem for the linear differential equation. The area of non-zero solutions to the above equations is determined by involving the solving transcendental equation, which is got by equating to zero of determinant related to the eigenvalue problem. The results of solving the nonlinear eigenvalue problems are applied subsequently for specification of the bifurcation points and obtaining the bifurcation curves. The approach does not depend on the form of operator determining the radiation properties of physical system (plane rectangular and hexagonal arrays). The obtained results are the constructive basis on which a series of practical engineering problems of optimization was solved numerically.

## **Conflict of interest**

The author declares no conflict of interest.

## **Author details**

Mykhaylo Andriychuk1,2

1 Pidstryhach Institute for Applied Problems of Mechanics and Mathematics. NASU, Lviv, Ukraine

2 Lviv Polytechnic National University, Lviv, Ukraine

\*Address all correspondence to: andr@iapmm.lviv.ua

© 2022 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

*Advanced Methods for Solving Nonlinear Eigenvalue Problems of Generalized Phase… DOI: http://dx.doi.org/10.5772/intechopen.103948*

## **References**

[1] Burge RE, Fiddy MA, Greenaway AH, Ross G. The phase problem. Proceedings of the Royal Society of London. 1976; **350**(1661):A350191-A350212. DOI: 10.1098/rspa.1976.0103

[2] Zenkova CY, Gorsky MP, Ryabiy PA, Angelskaya AO. Additional approaches to solving the phase problem in optics. Applied Optics. 2016;**55**:B78-B84. DOI: 10.1364/AO.55.000B78

[3] Taylor G. The phase problem. Acta Crystallographica Section D. 2003; **59**(11):1881-1890. DOI: 10.1107/ S0907444903017815

[4] Hauptman HA, Langs DA. The phase problem in neutron crystallography. Acta Crystallographica Section A. 2003; **59**(3):250-254. DOI: 10.1107/ S010876730300521X

[5] Dumber AS. On the theory antenna beam scanning. Journal of Applied Physics. 1958;**13**(5):31-39

[6] Tartakovskyi LB, Tikhonova VK. Synthesis of a linear radiator with a given distribution of amplitudes. Radiotechnica & Electronica. 1959;**4**(12): 2016-2019 (In Russian)

[7] Chony YI. To the synthesis of an antenna system for a given amplitude radiation pattern. Izvestiya Vysshikh Uchebnykh Zavedenii. Radioelektronika. 1968;**11**(2):1325-1327 (In Russian)

[8] Sklyanin EK. The method of the inverse scattering problem and the quantum nonlinear Schrödinger equation. Doklady Akademii Nauk SSSR. 1979;**244**(6):1337-1341 (In Russian)

[9] Ramm AG. Multidimensional Inverse Scattering Problems. New York: Longman Scientific & Wiley; 1992. p. 385 [10] Ikehata M. Reconstruction of an obstacle from the scattering amplitude at a fixed frequency. Inverse Problems. 1998;**14**:949-954

[11] Colton D, Kirsch A. A simple method for solving inverse scattering problems in the resonance region. Inverse Problems. 1996;**12**:383-393

[12] Precup R. Methods in Nonlinear Integral Equations. Alphen aan den Rijn, Netherlands: Kluwer; 2002

[13] Masujima М. Applied Mathematical Methods in Theoretical Physics. Weinheim, Germany: Wiley-VCH; 2005

[14] Bauschke HH, Combettes PL, Luke DR. Phase retrieval, error reduction algorithm, and Fienup variants: A view from convex optimization. Journal of the Optical Society of America. 2002;**19**(7): 1334-1345

[15] Fienup JR. Phase retrieval algorithms: A comparison. Applied Optics. 1982;**21**:2758-2769

[16] Weinberg MM, Trenogin VA. Branching Theory of Solutions to Nonlinear Equations. Moscow: Nauka; 1969. p. 528 (In Russian)

[17] Deuflhard P. Newton Methods for Nonlinear Problems.Affine Invariance and Adaptive Algorithms.Springer Series Computational Mathematics. Berlin, Heidelberg: Springer-Verlag; 2004. p. 35: xii+424

[18] Güttel S, Tisseur F. The nonlinear eigenvalue problem. ActaNumerica. 2017;**26**:1-94. DOI: 10.1017/ S0962492917000034

[19] Voss H. Nonlinear eigenvalue proble ms. In: Hogben L, editor. Handbook of

Linear Algebra. 2nd ed. Boca Raton, FL: Chapman and Hall/CRC; 2014. p. 1904

[20] Ruhe A. Algorithms for the nonlinea r eigenvalue problem. SIAM Journal on Numerical Analysis. 1973;**10**(4):674-689. DOI: 10.1137/0710059

[21] Mehrmann V, Vos H. Nonlinear eigenvalue problems: A challenge for modern eigenvalue methods. GAMM-Mitteilungen. 2004;**27**(2):121-152. DOI: 10.1002/gamm.201490007

[22] Jittorntrum K. An implicit function theorem. Journal of Optimization Theory and Applications. 1978;**25**(4):575-577. DOI: 10.1007/BF00933522

[23] Savenko PA, Protsakh LP. Implicit function method in solving a twodimensional nonlinear spectral problem. Russian Mathematics. 2007;**51**(11): 41-44. DOI: 10.3103/ s1066369x07110060

[24] Savenko P, Tkach M. Numerical approximation of real finite nonnegative function by the modulus of discrete Fourier transformation. Applied Mathematics. 2010;**1**:41-51

[25] Andriychuk MI, Voitovich NN, Savenko PA, Tkachuk VP. Antenna Synthesis According to the Amplitude Radiation Pattern. Numerical Methods and Algorithms. Kyiv: NaukovaDumka; 1993. p. 256 (In Russian)

[26] Savenko PO. Nonlinear Problems of the Radiation System Synthesis. Theory and Methods of Solution. Lviv: IAPMM of NASU; 2002. p. 320 (In Ukrainian)

[27] Kravchenko VF, Protsakh LP, Savenko PA, Tkach MD. Mathematical peculiarities of plane equidistant array synthesis by given amplitude radiation pattern. Antennas. 2010; **3**(154):34-48

[28] Gantmacher FR. The Theory of Matrices.Volume One. New York: The Chelsea Publishing Company; 1959. p. 276

[29] Savenko P. Computational methods in the theory of synthesis of radio and acoustic radiating systems. Applied Mathematics. 2013;**4**:523-549. DOI: 10.4236/am.2013.43078

[30] Andriychuk MI, Savenko PO, Tkach MD. Synthesis of plane equidistant array according to power radiation pattern. In: Proceedings of XVIIth International Seminar/Workshop on Direct and Inverse Problems of Electromagnetic and Acoustic Wave Theory (DIPED-2012); 24-27 September 2012. Tbilisi, New York: IEEE; 2012. pp. 68-74

[31] Andriychuk MI, Voitovich NN. Antenna synthesis according to power radiation pattern with condition of norm equality. In: Proceeding of 2013 XVIIIth International Seminar/Workshop on Direct and Inverse Problems of Electromagnetic and Acoustic Wave Theory (DIPED); 23–26, September 2013. Lviv, New York: IEEE; 2013. pp. 137-140

[32] Andriychuk MI, Kravchenko VF, Savenko PO, Tkach MD. The plane radiation system synthesis according to the given power radiation pattern. FizicheskiyeProblemyPriborostroeniya. 2013;**2**(3):40-55 (In Russian)

[33] Andriychuk M, Savenko P, Tkach M. Non-linear synthesis problems for plane radiating systems according to the prescribed power directivity pattern. Open Journal of Antennas and Propagation. 2013;**1**:23-34. DOI: 10.4236/ ojapr.2013.12006

[34] Andriychuk MI. Antenna Synthesis through the Characteristics of Desired Amplitude. Newcastle, UK: Cambridge Scholars Publishing; 2019. p. xvi+150

## **Chapter 12**
