Minimum Variance Distortionless Response Beamformer with Enhanced Nulling Level Control via Dynamic Mutated Artificial Immune System

In smart antenna applications, the adaptive beamforming technique is used to cancel interfering signals (placing nulls) and produce or steer a strong beam toward the target signal according to the calculated weight vectors. Minimum variance distortionless response (MVDR) beamforming is capable of determining the weight vectors for beam steering; however, its nulling level on the interference sources remains unsatisfactory. Beamforming can be considered as an optimization problem, such that optimal weight vector should be obtained through computation. Hence, in this paper, a new dynamic mutated artificial immune system (DM-AIS) is proposed to enhance MVDR beamforming for controlling the null steering of interference and increase the signal to interference noise ratio (SINR) for wanted signals.


Introduction
The evolution of adaptive beamforming was initiated in military and aerospace applications through the employment of electronically steered antennas based on phasedarray technology. Typical applications include long-range surveillance radar, active jammer rejection, and multibeam antennas for space communications [1]. The same antenna array techniques were then assumed suitable for mobile radio communication to solve multipath fading and the cochannel interference problem. A partially adaptive antenna array technology, known as the intermediate frequency side lobe canceller (SLC), was invented by Howells in the late 1950s [2]. This technology incorporates the capability of automatic interference nulling. However, SLC was not fully adaptive because the main beam has a fixed pattern, and the auxiliary array contains only a few controlled elements. This simple adaptive antenna facilitated the development of fully adaptive array by Howells' coworker, Applebaum, in 1965.
The algorithm, commonly known as the Howells-Applebaum algorithm, was developed to maximize the signal-to-noise ratio (SNR) at the output of the beamformer. At the same time, in 1960, another independent research group led by Windrow invented another adaptive array approach based on linear covariance minimum variance (LCMV) [3]. This algorithm, later known as the Windrow-Hoff LMS algorithm, was developed based on the minimum mean square error (MMSE) criterion for the automatic adjustment of array weights. This algorithm is well known for its simplicity but only achieves satisfactory performance under specific operational conditions. The major drawback of this algorithm is its low convergence rate, which refers to the speed by which the mean of the estimated weights approaches the optimal value, thereby making it unsuitable for certain applications.
In 1969, Capon introduced another adaptive beamformer known as minimum variance distortionless response (MVDR), which is a technique capable of resolving signals 2 The Scientific World Journal separated by a fraction of antenna beam width. This optimum beamformer, also known as the LCMV beamformer, requires only the knowledge of the desired signal direction of arrival (DOA) to maximize the SNR. Another important contribution was by Reed, Mallett, and Brennen in 1974. They introduced a fast convergence algorithm known as the sample matrix inversion (SMI) technique, which overcame the problem of slow convergence faced by the LMS algorithm. One of the most important factors in smart antenna processes is beamforming, which refers to the allocation of signals in particular positions and phase angles for each antenna for the corresponding angle of the system [4].
Beamforming technology is a key technique in nulling antenna. By performing some processes with the received array signals, such as weighting and summation, beamforming can help the antenna realize many advanced functions, such as beam shaping, beam scanning, and beam nulling [5]. Reception beamforming is independently achieved at each receiver; however, the transmitter in transmit beamforming has to consider all receivers to optimize the beamformer output [6,7].
One of the beamforming algorithm used in smart antenna is MVDR beamforming, which can calculate the weight vector to determine the desired signal from the interference. Moreover, MDVR maximizes the sensitivity in one direction only [8]. The MVDR beamformer, also known as Capon beamformer, minimizes the output power of the beamformer under a single linear constraint on the response of the array toward the desired signal. The idea of combining multiple inputs in a statistically optimum manner under the constraint of no signal distortion can be attributed to Darlington. Several researchers developed beamformers, in which additional linear constraints are imposed (e.g., Er and Cantoni) [9,10].
Many approaches proposed the use of a mathematical model to improve the robustness of the MVDR beamformer, as presented in [11,12]. Research on the artificial immune system (AIS) and its application has become increasingly important in the field of intelligent information systems [13][14][15]. A new optimization technique was presented for the design of linear antenna arrays. The proposed technique was based on a novel variant of particle swarm optimization (PSO) called Boolean PSO with adaptive velocity mutation. The antenna arrays were optimized based on the requirements for maximizing the power gain at a desired direction and minimizing the side lobe level of the radiation pattern [16]. A complex-valued genetic algorithm for the optimization of beamforming in linear array antennas was proposed. The algorithm was proven to enhance searching efficiency significantly while effectively avoiding premature convergence. Numerical results were presented to illustrate the advantages of the proposed technique over conventional pattern synthesis methods [17].
In this paper, the main goal is to design a beamforming method based on MVDR in corporation with new dynamic mutated artificial immune system (DM-AIS) algorithm in order to enhance the null level at interference sources. By using this method, finding a new mathematical model, changing the filter hardware for signal processing, or changing the design of antenna based on the increased number of elements is no longer necessary. Another reason for proposing the new method is the difficulty in obtaining the optimum value using any normal algorithm without intensification. In this investigation, DM-AIS have been applied in beamforming with uniform linear antenna arrays of 0.5 spacing between adjacent elements and radiating at a frequency of 2.3 GHz. The rest of this paper is organized as follows. Section 2 introduces the basics of adaptive beamforming. Sections 3 and 4 summarized a basis to describe the conventional MVDR beamforming and AIS, respectively. Section 5 shows the incorporation of MVDR with DM-AIS. Simulation results of one user with two interferences and comparison of conventional MVDR with mp-QP MVDR and DM-AIS are reported in Section 6. And finally Section 7 concludes this investigation.

Background of Adaptive Beamforming
Adaptive beamforming is a technique for receiving a signal of interest (SOI) from specific directions while suppressing the interfering signals adaptively in other directions using an array of sensors. This technique can automatically optimize the array pattern by adjusting the elemental control weights until a prescribed objective function is satisfied. This technique provides a means for separating a desired signal from interfering signals.
Beamforming has numerous applications in radar, sonar, seismology, microphone array speech processing, and, more recently, wireless communications. In particular, the use of antenna arrays, in combination with signal processing algorithms at the base station, offers the possibility of exploiting the spatial dimension to separate multiple cochannel users. This approach provided increased channel capacity and wider area coverage. Array beamforming methods in such systems use the spatial dimension to combat interference, noise, and multipath fading of the desired signal [11]. The outputs of the individual sensors were linearly combined after being scaled with the corresponding weights. This process optimizes the antenna array to achieve maximum gain in the direction of the desired signal and nulls in the direction of interferers. For a beamformer, the output at any time , ( ) is given by a linear combination of the data at antennas, with ( ) being the input vector and ( ) being the weight vector, as shown in Figure 1: Weight vector ( ) can be defined as follows: For any algorithm that evades the matrix inverse operation and uses the immediate gradient vector ∇ ( ) for weight The Scientific World Journal vector upgrading, the weight vector at time + 1 can be written as follows: where is the step size parameter, which controls the speed of convergence and lies between 0 and 1. The minimum values of facilitate the sluggish concurrence and highquality estimation of the cost function. Comparatively, the huge values of may direct to a rapid union. However, the constancy over the least value may disappear: An exact calculation of instantaneous gradient vector ∇ ( ) is not possible because prior information on covariance matrix and cross-correlation vector is needed. Thus, an instantaneous estimate of gradient vector is given by By substituting values from (5) into (3), the weight vector is derived as follows: The desired signal can be defined by the following three equations: Numerous algorithms were introduced for the design of an adaptive beamformer [8]. One of the most popular approaches for adaptive beamforming was proposed by Capon [4]. His algorithm leads to an adaptive beamformer with an MVDR. Some constraints, such as the antenna gain being constant in the desired direction, are used to ensure that the desired signals are not filtered out along with the unwanted signals. The MVDR beamformer not only minimizes the array output power but also maintains a distortionless main lobe response toward the desired signal. However, the MVDR beamformer may have an unacceptably low nulling level, which may significantly degrade the performance in the case of unexpected interfering signals. In particular, the performance of MVDR degrades in rapidly moving jammer environments. This degradation occurs because of jammer motion, which may bring jammers out of the sharp notches of the adapted pattern. To achieve high interference suppression and SOI enhancement, an adaptive array must introduce deep and widened nulls in the DOAs of strong interferences while keeping the desired signal distortionless. Thus, the issue of nulling level control is especially important for both deterministic and adaptive arrays [18].

Conventional MVDR Beamforming
When a beamformer has a constant response in the direction of a useful signal, the LCMV algorithm becomes an MVDR algorithm [19]. The MVDR algorithm is capable of suppressing the interference, but with high value in SNR and low noise. At the same time, the MVDR algorithm depends on the steering vectors, which in turn depend on the incident angle of the received signal from the element of the array antenna. The direction of useful signal must be known and the output power subject to a unity gain constraint in the direction of desired signal must be minimized. The array output is given by The output power is as follows: where the covariance matrix should be ( , 1) for the received signal and is the hermitian transpose.
The optimum weights are selected to minimize the array output power MVDR while maintaining unity gain in the look direction ( ), which is the steering vector of the desired signal. The MVDR adaptive algorithm can be written as follows: The steering vector ( ) is given by The Scientific World Journal where is the space between the elements of the antenna, is the desired angle, and is the number of elements, as shown in Figure 2.
The optimization weight vectors can then be acquired using the following formula [20]: These weights are the solution of the optimization problem mentioned at (10) and with the use of four elements of array antenna, four weights as below will be obtained: ] .
Subsequently, the beamformer weights are selected based on minimum mean value of output power according to the number of users inside the coverage area while maintaining unity response in the desired direction. Nevertheless, the restraint ensures that the signal passes through the beamformer undistorted. Consequently, the output signal power is similar to the look direction source power. The total noise, including interferences and uncorrelated noise, is then reduced by the minimization process. Notably, the minimization of the total output noise, while constantly maintaining the output signal, is the same as maximizing the output SINR. However, for the optimal beamformer to perform as described above and to maximize the SINR by cancelling interferences, the number of interferences must be less than or equal to − 2 because an array with elements has − 1 degrees of freedom and has been utilized by the constraint in the look direction. Given that the MVDR beamformer maximizes sensitivity in one direction only, this beamformer is unsuitable for multipath environments, where the desired signal spreads in all directions [21]. The multipath occurs in non-line-of-sight environments such as populated urban areas, where numeorus scatterers are close to the users and the base station. Thus, the MVDR beamformer may have an unacceptably low nulling level, which may significantly degrade performance in the case of unexpected interfering signals. As a result, the beamforming optimization problem is formulated as a multiparametric quadratic programming (mp-QP) [18].

Artificial Immune System
The proposed usage of artificial intelligent system as the enhancement method for the adaptive beamforming technique is based on a framework that is built around the concept of reactive artificial immune system (AIS). AIS, which is inspired by theoretical immunology and observed immune functions, is a branch of the metaheuristic algorithm with promising results in the field of optimization. While AIS resembles some other metaheuristics algorithms such as genetic algorithm (GA), except the recombination operator, the former has code simplicity and is low in computational cost. This AIS system is indeed based upon the normal human immune system in the way that it is reactive towards foreign elements. The immune system is highly robust, adaptive, inherently parallel, and self-organized. It has powerful learning and memory capabilities and presents an evolutionary type of response to infectious foreign elements [22,23].
The main agents of the adaptive immune system are lymphocytes that are are called the B cells which produce antibodies to attack the enemy. Some B cells become "memory cells" which keep molecular records of past invader and minimize the body's response time to an infection. The clonal selection principle is the whole process of antigen recognition, cell proliferation, and differentiation into a memory cell.
The clonal-selection theory proposes that as an antigen enters the immune system certain B cells are selected based on their reaction to this antigen to undergo rapid cloning and expansion. This reaction is often termed the affinity of that B cell (or antibody) for the given antigen. Those B cells are selected, based on their affinity to an antigen, to produce a number of clones to attack or neutralize the invading antigen. Cells that have a higher degree of affinity are allowed to produce more clones. The clonal production develops immune cells that are more adept at recognizing and reacting to the antigen through mutation. B cell offsprings undergo mutation based on an inverse proportionality to their affinity values. Through this process, the affinity of subsequent generations of B cells will have greater reaction to the antigen, and more diversity will also have been added to the system through the wider exploration afforded by the high mutation rates of the cells with lower affinity measures.
The Scientific World Journal 5 The process of a standard clonal selection algorithm can be summarized as follows [24,25].
(1) Generate a random initial population of antibody , given by (2) Compute the fitness of each (3) Generate clones by cloning all cells in the population. The amount of clone is given by where is the total clones generated for each , is the multiplying factor, and is the number of .
The next generation of is obtained using = + selection by choosing the best individuals out of the population.
(7) Steps 3 to 6 are repeated until a predefined stopping condition is reached.

MVDR Beamforming Incorporation with New Dynamic Mutated Artificial Immune System
In this paper, AIS with dynamic mutation (DM) was utilized to enhance the null level of the MVDR beamforming technique. In analogy, the adaptive antenna represents the body of an organism, whereas the interference and noise signal sources represent external harmful attacks toward the organism. Hence, the adaptive antenna system will organize its antibody to protect the body of organism from external antigen attacks. The adaptive antenna system will try to optimize through its AIS iteration process to develop deep null at the DOA of the interference sources to achieve the maximum SINR. In this AIS algorithm, the weight vector will be generated as the system antibody. The algorithm will initiate by generating a population of antibodies, which is represented by weight vectors . The number of generated weight vectors depends on the population size size . For the first iteration, the first set of weight vectors 1 is obtained from the computation of the conventional MVDR weight vector. The weight vectors in every antibody contain an number of weight vectors, depending on the number of sensors or antenna elements used, and can be expressed as follows: In matrix format, the weight vectors in the population of any iteration can be represented by where is weight vectors of total population with sensors in each antenna; mvdr is weight vectors from MVDR beamformer; is number of sensor. In this study, antenna sensor of (4, 1) is used; size is population size; is number of sensor.
Each set of antibodies has amplitude and phase ( ∀ ) to steer the radiation beam toward its target user and place the deep null toward the interference sources to achieve the optimum SINR. The best weight vector is determined according to the fitness value obtained from fitness function as shown below: where User = power of target user, Inter = power of interference, and = number of interference sources.
The new candidate population of antibodies (weight vectors) based on the mutation rate depends on the fitness value of in (22), which is given by where rate = mutation rate and best = fitness of best population weight. 6 The Scientific World Journal The mutation rate of clones is inversely proportional to their antigenic affinity [24,26]. A higher affinity denotes a smaller mutation rate.
The mutations and clones needed to create the new weight are given by The antibodies generate clones proportional to their affinities. The number of cloned antibodies can be computed by where is a multiplying factor with a value of 1 and is the weight vector. Each term of this sum corresponds to the clone size of each selected antibody. The clones are then subjected to hyper mutation and receptor editing. To enable the algorithm to find the best solution, a deep null must be achieved to obtain the maximum SINR. Therefore, the DM was introduced to identify the best solution. The selection of the 10 best solutions depends on the maximum SINR value. Thus, from these 10 best solutions, DM-AIS can select three random values from the 10 best weights to achieve the SINR size. When the scaling factor is (0, 1, 2), the best value of the deep null is derived: where 1 , 2 , 3 = random weight select from the best three solutions and sf = scaling factor from (0, 1, 2). The next step is the clonal operation to obtain the best solution.
Affinity is applied to achieve the maximum power for target and deep null for interference. If this value is the best, it will store and yield the final weight value to stop the calculation. The DM-AIS steps in an adaptive antenna are as follows.
(1) A random initial population of is generated, which is given by (2) The fitness of each is computed: (3) Clones are generated by cloning all the cells in the population. The amount of clone is obtained using (23). (4) The clone population is mutated to produce a mature clone population with number of weight. The rate of mutation is given by Thus, the new is composed of ( ) = { 1 , . . . , }.
The next generation of is obtained using (20). The best best is selected to compose the new new population as follows: (6) The 10 best weight vectors are selected depending on the selection assumption ( size ) to identify the best number of the weight.
(7) The vectors are mutated by selecting three random weights from the 10 best weights using (26).
(8) Steps 3 to 6 are repeated until a predefined stopping condition is reached.
Remarks. Similar to weight vector optimization techniques based on the support vector regression (SVR-AS) [27] and mp-QP MVDR [18], the proposed DM-AIS MVDR technique is also effective in finding the optimal weight vector for an antenna array under specific conditions. The differences among the methodologies are as follows.
(1) The proposed DM-AIS MVDR determines the optimal weight vector of an antenna array through its cloning and mutation process to fine-tune the weight vector toward a radiation pattern that can generate deep null toward the interference source, while maintaining high gain at the target users' directions for an instantaneous scenario.
(2) The mp-AP MVDR produces the optimal weight vector of an antenna array for the optimal problem of the receiving beam based on a multiparametric approach [18].
(3) The SVR-AS requires the radiation pattern as the train data to find the optimal weight vector of an antenna array [27].
The benefits of the proposed approach are as follows.
(2) This approach can guarantee that the nulling levels in the specified areas are better than conventional MVDR.
(3) This approach does not require pretrain data and does not require complex mathematical computation to identify the optimal weight vector. Instead, merely the DOAs of the target user and interference sources are required for DM-AIS MVDR to determine its optimal weight vector for best radiation beamforming.

Simulation and Results
6.1. One User with Two Interferences. In this section, simulations were conducted to validate the proposed approach.  The uniform linear array (ULA) consists of four elements ( = 4) equispaced by half-wavelength. The desired signal and two interference signals are plane waves impinging on the ULA from the directions 50 ∘ , 30 ∘ , and 40 ∘ , respectively. In this simulation, the SNR from 5 dB to 30 dB was used for the desired signal and the two interference signals. The beam pattern nulling level should be below −70 dB. The complex vector of beamformer weights calculated by the aforementioned (27)-(31) is presented in Table 2, whereas the beam patterns generated are plotted in Figure 3. All the beam patterns have nulls at the DOAs of the interference signals and maintain a distortionless response for the SOI. However, DM-AIS place deep nulls (with nulling level equal to −80 dB) at the DOAs of two interference signal sources. The MVDR after 10 iterations reduces nulling levels compared with the MVDR after 20 iterations (Table 1). Figure 3 gives the ratio of SINR 67.9479 dB in 20 iterations while at 10 iterations the SINR is 40.5387 dB of the aforementioned beamformer for the previous scenario. It can be seen that the DM-AIS at 20 iterations shows better ratio which is 122.49% of improvement. If we compare the improvement in SINR between normal AIs and DM-AIS for the same scenario in 10 iterations we obtain the percentage 131.48%.
But an insignificant increase in the SINR of 0.2 is obtained for an increase in the number of iterations. The iterations were taken from 30 till 50, as shown in Figure 4. This implies that there is no need to increase the time required for simulations. The 20 iterations will give the best results by the short time. These results should be saved. The simulated results demonstrate that the DM-AIS with dynamic mutation rate give a good result effectively at a faster speed. The simulated results are dynamically adapting its value when the SINR changes. In comparison, the dynamic mutation rate has shown better results with finer resolutions.
We can then discriminate the difference in results after using the new method with DM-AIS. To justify the results, the new method must be proven to be more effective than any previously suggested method. The new method should   Table 3 and Figure 5 provide details for two users at angles 0 ∘ and 10 ∘ with interference at angles 50 ∘ , −50 ∘ , 30 ∘ , and −30 ∘ .

8
The Scientific World Journal

Comparison of Conventional MVDR with MP-QP MVDR and DM-AIS.
To prove the importance of this project, the results of this work were those of previous work to enhance the MVDR in smart antennas. From the studied literature, no researcher has discussed or addressed beamforming by using four elements in the smart antenna or by applying DM-AIS in the smart antenna. Therefore, the results of this project are compared with robust MVDR beamforming for nulling level control via multiparametric quadratic programming results using 10 elements and with the direction of user at 0 ∘ with interference at 40 ∘ and −40 ∘ [8]. Figure 6 and Table 4 show the results. Figure 6 shows that all beam patterns have nulls in the interference signals and maintain a distortionless response for the SOI. However, DM-AIS MVDR places deep nulls (with nulling level equal to −90 dB) at the two interference signal sources, whereas mp-QP MVDR obtained a null level equal  to −80 dB. Therefore, the DM-AIS with MVDR response presents lower nulling levels compared with the new mp-QP MVDR. The maximum SINRs calculation for each technique is compared with one another. The maximum value obtained from DM-AIS is shown in Table 5. This result shows that the new beamforming using artificial intelligence to determine the desired signal of user is effective and that mathematical equations or filter hardware signal processing are unnecessary. For the proposed approach, the weights value vectors make it difficult to obtain the optimum value by using any other algorithm without intensification.

Conclusion
A new DM-AIS was presented and applied in adaptive beamforming with four elements of linear antenna arrays. The proposed DM-AIS was able to enhance the MVDR technique through further optimization of weight vector, which aimed to control the nulling level of interference and the directionality of the desired signal. Very low levels of interference with good accuracy were achieved even in the case of multiple users or multiple interferences. The results of the proposed approach were compared with those of the mp-QP MVDR and conventional MVDR, and the effectiveness of the proposed approach in minimizing the power of interference and increasing SINR was observed. Finally, the DM-AIS can be useful to antenna engineers The Scientific World Journal 9 for the pattern synthesis of antenna arrays because of its good accuracy and the lack of a requirement for complicated mathematical functions.