# An Application of Artificial Intelligence for the Joint Estimation of Amplitude and Two-Dimensional Direction of Arrival of Far Field Sources Using 2-L-Shape Array

Authors: Fawad Zaman; Ijaz Mansoor Qureshi; Junaid Ali Khan; Zafar Ullah Khan
Journal: International Journal of Antennas and Propagation
 (2013)

---

## Abstract

An easy and efficient approach, based on artificial intelligence technique, is proposed to jointly estimate the amplitude, elevation, and azimuth angles of far field sources impinging on 2-L-shape array. In these proposed artificial intelligence techniques, the metaheuristics based on genetic algorithm and simulated annealing are used as global optimizers assisted with rapid local version of pattern search for optimization of the adaptive parameters. The performance metric is employed on a fitness evaluation function depending on mean square error which is optimum and requires single snapshot to converge. The proposed approaches are easy to understand, and simple to implement;  the genetic algorithm specifically hybridized with pattern search generates fairly good results. The comparison of the given schemes is carried out with 1-L-shape array, as well as, with parallel-shape array and is found to be in good agreement in terms of accuracy, convergence rate, computational complexity, and mean square error. The effectiveness and efficiency of the given schemes are examined through Monte Carlo simulations and their inclusive statistical analysis.


---

## Body

## 1. Introduction

Two-dimensional direction of arrival (2D DOA) estimation for multiple signals received at sensors array is one of the indispensable steps for adaptive beamformer to maneuver the main beam in any preferred direction, while allocating nulls in the direction of unwanted signals [1, 2]. From this perspective, 2D DOA estimation has direct applications in radar, sonar, and mobile communications [3]. The most common problem in 2D DOA estimation includes pair matching between elevation and azimuth angles, estimation failure and computational complexity. In the literature, many algorithms have already been discussed in the last two decades to estimate 2-D DOA [4–7], but they have the aforementioned problems. In [8], Wu et al. proposed an algorithm based on propagator method (PM) to overcome the computational load of [7], but it fails to avoid the pair matching problem and estimation failure problem in the range of 1.2217 (radians) to 1.5708 (radians). Moreover, the parallel-shape array is used in [8] which not only needs more sensors but also requires a large number of snapshots to achieve the goal (at least 200 snapshots per sensor are required). Besides that, [8] also failed to estimate the sources amplitudes which also an important parameter to be estimated.In order to overcome these problems, we thank evolutionary computing techniques (ECTs) and L-shape arrays (1-L and 2-L shape). It is well acknowledged in the literature that ECTs such as genetic algorithm (GA), ant colony optimization (ACO), differential evolution (DE), and particle swarm optimization (PSO) have achieved great success in solving different optimization problems especially in the field of array signal processing [9, 10]. These techniques have got attention among the researchers due to their ease in understanding, simplicity in implementation, and less probability to get stuck in the local minima. One of the other most significant aspects of ECT is that their accuracy and reliability increase more when hybridized with any other capable algorithm such as pattern search (PS), interior point algorithm (IPA), and active set algorithm (ASA) [11–15]. In [13], GA is hybridized with IPA for 3D source localization of near field sources impinging on uniform linear array (ULA), and the hybrid GA-IPA approach has shown better performance as compared to GA and IPA alone. In [14], PSO is mixed with PS to estimate the DOA of far field sources, and again the hybrid PSO-PS technique produced better results as compared to PSO [9] and PS alone. GA is hybridized with PS to jointly estimate 3D parameters (amplitude, elevation angle, and azimuth angle) of far field sources arriving at 1-L shape array and again the results of hybrid GA-PS technique were remarkable as compared to GA and PS alone [15]. In [13–15], every time GA and PSO are treated as global optimizers, while PS and IPA are used as local search optimizers and mean square error (MSE) is used as a fitness function.In this paper, an easy and efficient approach based on hybridized GA with PS (GA-PS) and hybrid simulated annealing (SA) technique with PS (SA-PS) are developed for joint estimation of amplitudes and 2-D DOA of far field sources impinging on 2-L-shape array. In these approaches GA and SA are exploited as global optimizers, while PS is incorporated for rapid local search to improve the results further. Mean square error (MSE) is used as an objective function based on an error between actual and estimated signals which is basically derived from maximum likelihood principle [13]. This fitness function needs single snapshot to achieve the results and avoids any uncertainty among the angles that are supplement to one another. The validity and effectiveness of the proposed approaches are tested on a large number of simulations which include different cases based on a different number of sources, a different number of sensors, computational complexity, and the proximity of the sources. Initially, the results of GA, PS, SA, SA-PS, and GA-PS are compared with each other for 2-L-shape array, and then, the entire results are compared with the results of 1-L-shape array [15]. Later on, the results are also compared with parallel shape array which uses PM [8].Throughout the paper, the matrices are denoted by capital bold letters, while vectors are expressed by small bold letters, whereas “T” denotes the transpose of vectors or matrices.The left behind paper is organized as follows:  the problem formulation is discussed in Section2, whereas its proposed methodology is provided in Section 3. Section 4 demonstrates the comprehensive simulations along with the discussion of the results, while conclusion and future directions for the research are revealed in Section 5.


## 2. Problem Formulation

In this section, we addressed the unknown parameters of the received signals. Consider a 2-L-shape array which consists of three uniform linear arrays (ULAs) placed alongx-axis, y-axis, and z-axis. Each ULA contains M-1 elements with the same interelement spacing “d”, while the reference sensor is common for all subarrays as shown in Figure 1 [16]. Consider K narrow band sources impinging on the 2-L-shape array where si,θi, and ϕi are the ith source amplitude, elevation angle, and azimuth angle, respectively.Figure 1
2-L shape array consists of three subarrays.Output Vector along
z
-Axis Subarray. The subarray along z-axis is used to estimate the elevation angle. The signal received at lth element from K sources can be written as
(1)wzl=∑i=1Kbzl(θi)si+αzl,
where l=0,1,2,…M-1.In matrix-vector form, the output of entire subarray can be written as(2)wz=Bz(θ)s+αz,
whereB is called steering matrix, which contains steering vectors; that is,
(3)Bz=[bz0,bz1,…bzM-1]T,
where
(4)bzl=[bzl(θ1),bzl(θ2),…bzl(θk)],bzl(θi)=exp(-jψl,i),ψl,i=2πldcosθiλfori=1,2,…K.Similarly, in (2),
(5)wz=[wz0,wz1,…wzM-1]T,s=[s1,s2,…sK]T,αz=[αz0,αz1,…αzM-1]T,
where αz is called additive white Gaussian noise (AWGN) added at the output of each sensor placed in z-axis subarray.Output Vector along
x
-Axis Subarray. The subarray along x-axis is used to estimate the x-component of azimuth angle and elevation angle. The output at lth sensor placed in the x-axis subarray for K sources can be written as
(6)wxl=∑i=1Kbxl(θi,ϕi)si+αxl.
In vector-matrix form, it can be written as
(7)wx=Bx(θ,ϕ)s+αx,
where
(8)Bx=[bx0,bx1,…bxM-1]T,bxl=[bxl(θ1,ϕ1),bxl(θ2,ϕ2),…bxl(θk,ϕk)],bxl(θi,ϕi)=exp(-jψl,i),ψl,i=2πldsinθicosϕiλ.
Similarly,
(9)wx=[wx0,wx1,…wxM-1]T,αx=[αx0,αx1,…αxM-1]T,
where αx is AWGN added at the output of each sensor in x-axis subarray.Output Vector along
y
-Axis Subarray. In the same manner, the y-axis subarray is used to estimate the y-component of elevation and azimuth angle. The output atlth sensor in this subarray is given as follows:
(10)wyl=∑i=1Kbyl(θi,ϕi)si+αyl,wy=By(θ,ϕ)s+αy,
where
(11)By=[by0,by1,…byM-1]T,
where
(12)byl=[byl(θ1,ϕ1),byl(θ2,ϕ2),…byl(θk,ϕk)],byl(θi,ϕi)=exp(-jψl,i),ψl,i=2πldsinθisinϕiλ.
Similarly,
(13)wy=[wy0,wy1,…wyM-1]T,αy=[αy0,αy1,…αyM-1]T,
where αy is AWGN added at the output of each sensor in y-axis subarray. Here, the problem in hand is very clear which is to estimate the unknown parameters, that is, amplitudes and 2D-DOA (elevation, azimuth) angles for K sources impinging from far field on 2-L-shape array.


## 3. Proposed Techniques

In this section, a brief introduction, parameter setting, and flow diagram are given for GA, PS, SA, GA-PS, and SA-PS.Simulated Annealing. Simulated annealing (SA) technique is used for both discrete and continuous optimization problems. SA works on those problems, which have controlled heating and cooling properties as evident from its name “annealing”, which means “To heat and then cool.” Due to easiness in implementation and the capability to avoid getting stuck in local minima, SA is a well-liked technique over the last few decades and has been successfully used in wide range of engineering problems [17]. In this paper, an MATLAB “built-in optimization tool box” is used for which the parameter setting is shown in Table 1. The best individual result of SA is further passed through PS for further tuning.Table 1
Parameter settings for GA, SA, and PS.
GA
PS
SA
Parameters
Settings
Parameters
Setting
Parameters
Setting
Population size
240
Poll method
GPS positive basis 2 N
Annealing function
Fast
No. of generations
1000
Polling order
Consecutive
Reannealing interval
100
Migration direction
Both ways
Maximum iteration
800
Temperature update function
Exponential temperature update
Crossover fraction
0.2
Function evaluation
16000
Initial temperature
100
Crossover
Heuristic
Mesh size
01
Data type
Custom
Function tolerance
10–12
Expansion factor
2.0
Function tolerance
10–12
Initial range
[0-1]
Contraction factor
0.5
Max iteration
2000
Scaling function
Rank
Penalty factor
100
Max function evaluations
3000 * number of variables
Selection
Stochastic uniform
Bind tolerance
10–03
Temperature update function
Exponential temperature update
Elite count
2
Mesh tolerance
10–06
Mutation function
Adaptive feasible
X tolerance
10–06
Hybrid function call interval
EndPattern Search. Pattern search (PS) was introduced by Hookes and Jeeves in 1961 which does not require the gradient or derivative of the problem and can be used for both local and global optimization problems. Basically, PS works on mesh which is defined according to some specific rules. If no improvement in cost function is achieved at the mesh points of current iteration, then the mesh is polished and the process is repeated. It has already got application almost in every field such as signal processing and soft computing [18].Genetic Algorithm. GA is basically different from previously discussed algorithm (SA and PS) and is applicable to a wide range of optimization problems. GA is a more prominent and proficient algorithm than any other evolutionary computing technique due to its ease in conception, ease in implementation, and more importantly less probability to get stuck in the presence of local optima. GA is being successfully applied to a wide range of problems having application in commerce to the leading scientific research [19, 20]. The structure of GA exploited for the joint estimation is given in Figure 2.Figure 2
Generic structure of the hybrid intelligent approach.The general settings consisting of population size, number of generations of the algorithm, function evaluations, and stoppage criteria are defined in Table1, along with some specific parameter setting/values of the three algorithms based on GA, PS, and SA.The logical steps for GA and GA-PS in the form of pseudo code are given in the following steps.Step 1 (initialization).
Generate randomlyQ number of chromosomes (particles) where the length of each chromosome is 3×K. In each chromosome, the first K genes represent amplitudes, the next K genes contain the elevation angles, while the last K genes represent the azimuth angles. Mathematically, the qth chromosome can be written as(14)dq=[sq,1,sq,2,…sq,K,θq,K+1,…θq,2K,ϕq,2K+1,…ϕq,3K].
Here, sqj∈R:La≤sqj≤Ua where La and Ua are the lower and upper bounds of signal amplitudes ∀q=1,2,…Q and j=1,2,…K. In the same way, the elevation and azimuth angles are
(15)θqj∈R:0≤θqj≤π,∀q=1,2,…Q,j=K+1,K+2,…2K,ϕqj∈R:0≤ϕqj≤2π,∀q=1,2,…Q,j=2K+1,2K+2,…3K.Step 2 (objective function).
MSE is used as an objective evaluation function, which is derived from MLP [13]. For qth chromosome, it can be written as
(16)E(q)=Ex(q)+Ey(q)+Ez(q),
where
(17)Ex(q)=1(M-1)∑x=1M-1|wx-w^xq|2,Ey(q)=1(M-1)∑y=1M-1|wy-w^yq|2,Ez(q)=1(M-1)∑z=1M-1|wz-w^zq|2.wz,wx, and wy are defined in (1), (3), and (7), respectively, while w^xq,w^yq, and w^zq at lth element of each subarray forqth chromosome are given as
(18)w^xlq=∑i=1Kexp[-jπlsin(dq,K+i)cos(dq,2K+i)]dq,i,w^ylq=∑i=1Kexp[-jπlsin(dq,K+i)sin(dq,2K+i)]dq,i,w^zlq=∑i=1Kexp[-jπlcos(dq,K+i)]dq,i,
where dq is defined in (14).Step 3 (termination criteria).
The termination criteria depending on the following conditions are achieved.(a)
If the objective function value is achieved, which is predefined,  that is,ε≤10-12.(b)
The total number of iterations is completed.(c)
The Tol Fun exceeded.(d)
Tol Con occurred.Step 4 (reproduction).
New population is generated by using the operators of crossover, elitism, and mutation selection as provided in Table1.Step 5 (hybridization).
In this step, the finest results got through GA are further given to PS for more improvement. The setting used for PS is also provided in Table1.Step 6 (storage).
Store global best of the current iteration which will be used for comparison and better statistical analysis, and repeat Steps2 to 5 for enough numbers of independent runs.


## 4. Results and Discussions

This section is broadly divided into two parts. In first part of simulations, GA, PS, SA, GA-PS, and SA-PS are examined for 2-L-shape array in terms of estimation accuracy, convergence rate, and at the same time, the results are compared with 1-L-shape array [15]. In the second part, the performance of the best scheme among the above-mentioned five schemes for 2-L-shape array is not only compared with 1-L-shape array but also with parallel-shape array [8]. Throughout the simulations, the distance “d” between two consecutive sensors in each subarray is kept λ/2. All the values of the DOA (elevation and azimuth angle) of sources are taken in radian (rad), and all the results for 2-L-shape array are carried out for threshold value of MSE 10−2. The MATLAB built-in “Optimization Tool Box” is used for GA, SA, and PS for which the settings are provided in Table 1. As in [15], no noise is added to the system for estimation accuracy and convergence rate. All the results are comprehensively examined for 100 independent runs.Case 1.
In this case, the estimation accuracy of GA, PS, SA, GA-PS, and SA-PS is discussed for two sources impinging on 2-L-type array. For better comparison and analysis, all the values of amplitudes and DOA are the same as in [15] for 1-L-shape array, so the amplitudes and DOA values are {s1,s2,θ1,θ2,ϕ1,ϕ2}={1,2,0.5236(rad),0.8727(rad),1.2217(rad),1.9199(rad)} where s1,θ1, and ϕ1 represent the amplitudes, elevation angle, and azimuth angles of first source and so on. The 2-L-shape array is composed of 4 sensors in which 1 sensor is placed along each subarray, while the reference sensor is common for them. As provided in Table 2, all the five mentioned techniques are quite capable of producing fairly good estimate of the actual values. However, among them, the hybrid approach GA-PS produced better results as compared to the other techniques used. The second and third best results are provided by GA and PS, respectively. The results of SA are also improved when hybridized with PS.
The results of 1-L-shape array for 2 sources are provided in Table3 which requires 7 sensors [15]. The 1-L-shape array in [15] also uses the same techniques. One can clearly deduce from the comparison of Tables 2 and 3 that 2-L-shape array produced more accurate results as compared to 1-L-shape array with less number of sensors.Table 2
Estimation accuracy of 2-L-shape array for 2 sources.
Scheme
s
1
s
2
θ
1 (rad)
θ
2 (rad)
ϕ
1 (rad)
ϕ
2 (rad)
Desired values
1.0000
2.0000
0.5236
0.8727
1.2217
1.9199
GA
1.0003
2.0003
0.5240
0.8731
1.2221
1.9203
PS
1.0020
2.0021
0.5259
0.8748
1.2239
1.9220
SA
1.0977
2.0178
0.5385
0.8889
1.2361
1.9284
SA-PS
1.0047
2.0047
0.5278
0.8770
1.2260
1.9242
GA-PS
1.0000
2.0000
0.5235
0.8726
1.2216
1.9198Table 3
Estimation accuracy of 1-L-type array for 2 sources [15].
Scheme
s
1
s
2
θ
1 (rad)
θ
2 (rad)
ϕ
1 (rad)
ϕ
2 (rad)
Desired values
1.0000
2.0000
0.5236
0.8727
1.2217
1.9199
GA
1.0008
2.0008
0.5243
0.8735
1.2225
1.9207
PS
1.0032
2.0033
0.5268
0.8759
1.2249
1.9231
SA
1.0196
2.0195
0.5432
0.8923
1.2413
1.9395
SA-PS
1.0063
2.0063
0.5299
0.8790
1.2281
1.9263
GA-PS
1.0003
2.0002
0.5233
0.8723
1.2220
1.9195Case 2.
In this case, the estimation accuracy of all the above-mentioned five techniques is discussed for three sources impinging on 2-L-shape array. This time the 2-L-shape array consists of 7 sensors; that is, 2 sensors are placed along each subarray, while the reference sensor is common for them. For better comparison with 1-L-type array, the same values of amplitudes and DOA are used as in [15], which are provided in Table 4. In this case, few local minima are observed due to which the performance of all five techniques especially SA, SA-PS, and PS is despoiled. However, again the hybrid GA-PS shows excellence in accuracy even in the presence of local minima. The second best result is shown by GA alone.
Now compare the results with 1-L-type array as provided in Table5 which requires thirteen sensors [15]. Now, one can observe that 2-L-shape array not only requires less number of sensors but also produces more accurate results as compared to 1-L-shape array.Table 4
Estimation accuracy of 2-L-type array for 3 sources.
Scheme
s
1
s
2
s
3
θ
1 (rad)
θ
2 (rad)
θ
3 (rad)
ϕ
1
(rad)
ϕ
2 (rad)
ϕ
3 (rad)
Desired values
1.0000
2.0000
3.0000
0.1745
0.8727
1.3090
0.5236
1.9199
2.4435
GA
1.0043
2.0043
3.0042
0.1788
0.8769
1.3132
0.5278
1.9243
2.4478
PS
1.0189
2.0190
3.0189
0.1934
0.9005
0.8917
0.5427
1.9389
2.4624
SA
1.0509
2.0508
3.0509
0.2254
0.9237
1.3598
0.5744
1.9708
2.4944
SA-PS
1.0342
2.0342
3.0343
0.2087
0.9069
1.3432
0.5581
1.9541
2.4778
GA-PS
1.0003
2.0004
3.0003
0.1748
0.8730
1.3094
0.5240
1.9202
2.4439Table 5
Performance of 1-L-type array for 3 sources [15].
Scheme
s
1
s
2
s
3
θ
1 (rad)
θ
2 (rad)
θ
3 (rad)
ϕ
1 (rad)
ϕ
2 (rad)
ϕ
3 (rad)
Desired values
1.0000
2.0000
3.0000
0.1745
0.8727
1.3090
0.5236
1.9199
2.4435
GA
1.0073
2.0073
3.0073
0.1818
0.8800
1.3163
0.5309
1.9272
2.4508
PS
1.0278
2.0277
3.0277
0.2023
0.9005
1.3368
0.5514
1.9477
2.4713
SA
1.0610
2.0611
3.0610
0.2355
0.9337
1.3700
0.5846
1.9809
2.5045
SA-PS
1.0432
2.0432
3.0432
0.2177
0.9159
1.3522
0.5668
1.9631
2.4867
GA-PS
1.0011
2.0011
3.0011
0.1756
0.8738
1.3101
0.5247
1.9210
2.4446Case 3.
In this case, the estimation accuracy is examined for 4-sources. The 2-L-shape array consists of 10 sensors; that is, 3 elements are placed along each subarray, while the reference element is common for them. As provided in Table6, again the hybrid GA-PS leads the edge over the remaining four techniques in terms of estimation accuracy. The results for 1-L-shape array for the same schemes and the same number of sources are provided in Table 7 which needs fifteen sensors [15]. One can again make out the advantages of 2-L-shape array over 1-L-type array in terms of accuracy and the number of sensors required.Table 6
Performance of 2-L-type array for 4 sources.
Scheme
s
1
s
2
s
3
s
4
θ
1 (rad)
θ
2 (rad)
θ
3 (rad)
θ
4 (rad)
ϕ
1 (rad)
ϕ
2 (rad)
ϕ
3 (rad)
ϕ
4 (rad)
Desired
1.0000
2.0000
3.0000
4.0000
0.2618
0.6109
1.0472
1.4835
1.6581
2.1817
2.7925
3.4034
GA
1.0102
2.0103
3.0103
4.0102
0.2721
0.6212
1.0576
1.4937
1.6683
2.1920
2.8028
3.4137
PS
1.0313
2.0312
3.0313
4.0314
0.2932
0.6423
1.0787
1.5148
1.6895
2.2132
2.8239
3.4348
GA-PS
1.0040
2.0042
3.0043
4.0040
0.2659
0.6151
1.0514
1.4877
1.6624
2.1858
2.7966
3.4077
SA
1.1011
2.1010
3.1012
4.1010
0.3628
0.7122
1.1483
1.5846
1.7594
2.2829
2.8937
3.5044
SA-PS
1.0787
2.0786
3.0785
4.0787
0.3405
0.6895
1.1260
1.5624
1.7369
2.2606
2.8713
3.4821Table 7
Performance of 1-L-type array for 4 sources [15].
Scheme
s
1
s
2
s
3
s
4
θ
1 (rad)
θ
2 (rad)
θ
3 (rad)
θ
4 (rad)
ϕ
1 (rad)
ϕ
2 (rad)
ϕ
3 (rad)
ϕ
4 (rad)
Desired
1.0000
2.0000
3.0000
4.0000
0.2618
0.6109
1.0472
1.4835
1.6581
2.1817
2.7925
3.4034
GA
1.0163
2.0163
3.0162
4.0163
0.2781
0.6272
1.0635
1.4998
1.6744
2.1980
2.8088
3.4197
PS
1.0425
2.0425
3.0426
4.0426
0.3043
0.6534
1.0897
1.5260
1.7006
2.2242
2.8350
3.4468
GA-PS
1.0083
2.0083
3.0083
4.0083
0.2701
0.6192
1.0555
1.4918
1.6664
2.1900
2.8008
3.4117
SA
1.1263
2.1263
3.1263
4.1164
0.3881
0.7372
1.1735
1.6098
1.7844
2.3080
2.9188
3.5306
SA-PS
1.0932
2.0932
3.0932
4.0932
0.3550
0.7041
1.1404
1.5767
1.7513
2.2749
2.8857
3.4966Case 4.
In this case, the convergence rate is evaluated for 2-L-shape array against a different number of sources in the presence of 10 dB noise. As shown in Figure3, the convergence rates of all schemes are degraded with the increase of sources. However, one can notice that the convergence rate of GA is remarkable in case of hybridization with PS. The second best convergence is produced by GA itself alone.
Now, by comparing with 1-L-shape array whose convergence rate is shown in Figure4, one can conclude that all the five schemes have better convergence rate for 2-L-shape array as compared to 1-L-shape array.Figure 3
Convergence rate versus number of sources for 2-L-shape array at 10 dB noise.Figure 4
Convergence rate versus number of sources for 1-L-shape array at 10 dB noise.Case 5.
In this case, we checked the validation of all techniques for more practical scenario; that is, we considered DOA on the reference axis of L-shape arrays. For this, the 2-L-shape array is composed of 12 sensors where the amplitudes and DOA of the four sources are(19){s1=2,θ1=1.5702(rad),ϕ1=0(rad)},{s2=5,θ2=1.2217(rad),ϕ2=6.2832(rad)},{s3=1,θ3=0(rad),ϕ3=4.7124(rad)},{s4=4,θ4=1.0472(rad),ϕ4=1.5702(rad)}.
As provided in Table 8, an obvious degradation can be observed in the estimation accuracy of each scheme for all the elevation angles taken near the reference axis, while for azimuth angle near the reference axis, the accuracy is slightly degraded. However, among all of them the GA-PS technique is less degraded as compared to the others.
Now, if we check the same DOA using 1-L-shape array, then one can see in Table9 that the estimation accuracy of all schemes is degraded more as compared to 2-L-shape array.
Up till now, the performance of GA, PS, SA, SA-PS, and GA-PS is discussed for both 1-L- and 2-L-type arrays, and it has been shown through various cases that GA-PS produced fairly good results for both arrays. So, from now onwards, our focus will be limited to the GA-PS technique only. In the upcoming second part of the simulations, we will compare the GA-PS technique used for 2-L-shape array not only with GA-PS technique used for 1-L-shape array but also with PM which utilizes parallel-shape array [8].Table 8
Performance of 2-L-type array for sources located on reference axis.
Scheme
s
1
s
2
s
3
s
4
θ
1 (rad)
θ
2 (rad)
θ
3 (rad)
θ
4 (rad)
ϕ
1 (rad)
ϕ
2 (rad)
ϕ
3 (rad)
ϕ
4 (rad)
Desired
2.0000
5.0000
1.0000
4.0000
1.5702
1.2217
0.0000
1.0472
0.0000
6.2832
4.7124
1.5702
GA
1.9897
5.0105
0.9897
4.0103
1.6633
1.2321
0.0838
1.0573
0.0192
6.3041
4.6914
1.5900
PS
2.0315
4.9688
0.9687
3.9687
1.7139
1.2485
0.1623
1.0785
0.0314
6.2500
4.7473
1.6090
GA-PS
2.0039
4.9958
1.0044
4.0041
1.6197
1.2259
0.0524
1.0514
0.0105
6.2919
4.7211
1.5830
SA
1.8989
5.1011
1.1010
3.8989
1.3439
1.3231
0.2234
1.1483
0.0489
6.3303
4.7613
1.5237
SA-PS
1.9213
4.9214
1.0783
3.9213
1.7453
1.3001
0.1955
1.1261
0.0401
6.3216
4.7490
1.5341Table 9
Performance of 1-L-type array for sources located on reference axis [15].
Scheme
s
1
s
2
s
3
s
4
θ
1 (rad)
θ
2 (rad)
θ
3 (rad)
θ
4 (rad)
ϕ
1 (rad)
ϕ
2 (rad)
ϕ
3 (rad)
ϕ
4 (rad)
Desired
2.0000
5.0000
1.0000
4.0000
1.5702
1.2217
0.0000
1.0472
0.0000
6.2832
4.7124
1.5702
GA
2.0163
4.9837
1.0162
3.9837
1.7069
1.2381
0.1396
1.0635
0.0357
6.3190
4.7483
1.6059
PS
2.0425
4.9575
1.0426
4.0426
1.7698
1.2643
0.1991
1.0897
0.0598
6.3430
4.7721
1.6299
GA-PS
1.9917
5.0083
1.0083
4.0083
1.4661
1.2300
0.1065
1.0555
0.0375
6.2798
4.7301
1.5910
SA
2.1263
5.1263
0.8737
3.8836
1.8500
1.3480
0.2741
1.1735
0.0792
6.3467
4.7741
1.5067
SA-PS
1.9068
5.0932
0.9068
4.0932
1.3334
1.3150
0.2287
1.1405
0.0668
6.3328
4.7563
1.265Case 6.
In this case, the estimation accuracy in terms of MSE is discussed for GA-PS technique against the increasing number of sources. For this, the 2-L-shape array is composed of 13 sensors, while the 1-L-shape array is consisting of 19 sensors. As shown in Figure5, the MSE is degraded for both L-shape arrays against the increasing number of sources. However, the MSE of GA-PS for 2-L-shape array is better than 1-L-shape array.Figure 5
MSE versus number of sources at 10 dB noise.Case 7.
In this case, Tables10, 11, and 12 list the variances, means, and standard deviations for parallel-shape array with PM method and L-shape arrays with hybrid GA-PS approach. For this, the elevation angle varies between 1.2217 (rad) and 1.5708 (rad) for a fixed azimuth angle of 0.6109 (rad) in the presence of 10 dB noise.
The performance of PM method with parallel-shape array is worse for this range of elevation angle, while the 2-L-shape array with GA-PS technique produced accurate estimation. The second best result is produced by GA-PS for 1-L-shape array. This range of elevation angles is very important for mobile communication, and hence, 2-L-shape array with GA-PS is a good choice to be used. Moreover, the range of elevation and azimuth angles for 2-L-shape array is(0,π) and (0,2π), respectively, while the parallel-shape array and 1-L-shape array have the range of elevation, and azimuth angles are (0,π/2) and (0,2π), respectively, so both arrays cannot be operated for elevation angles beyond π/2. For more details, the readers are encouraged to see [21].Table 10
Means, variances, and standard deviations at 10 dB noise for different elevation angles and fixed azimuth angle by using PM with parallel-shape array [8].
θ in radians forϕ = 0.6109 (rad)
Mean ofθ (rad)
Variance ofθ (rad)
Standard deviation ofθ (rad)
1.2392
1.2227
0.0123
0.0146
1.2915
1.2728
0.0171
0.0167
1.3439
1.3207
0.0228
0.0200
1.3963
1.3696
0.0403
0.0265
1.4486
1.4155
0.0822
0.0379
1.5010
1.4552
0.1250
0.0467
1.5533
1.4784
0.1631
0.0534Table 11
Means, variances, and standard deviations at 10 dB noise for different elevation angles and fixed azimuth angle by using GA-PS with 1-L-shape array [15].
θ in radians forϕ = 0.6109 (rad)
Mean ofθ (rad)
Variance ofθ (rad)
Standard deviation ofθ (rad)
1.2392
1.2394
2.0963
e
-
006
1.9128
e
-
004
1.2915
1.2913
2.2918
e
-
006
2.0000
e
-
004
1.3439
1.3436
2.4819
e
-
006
6.5816
e
-
004
1.3963
1.3965
3.8608
e
-
006
2.5959
e
-
004
1.4486
1.4493
2.1907
e
-
006
1.9554
e
-
004
1.5010
1.5012
3.6759
e
-
006
2.5329
e
-
004
1.5533
1.5535
2.3424
e
-
006
2.0219
e
-
004Table 12
Means, variances, and standard deviations at 10 dB noise for different elevation angles and fixed azimuth angle by using GA-PS with 2-L-shape array.
θ in radians for ϕ = 0.6109 (rad)
Mean ofθ
Variance ofθ
Standard deviation ofθ
1.2392
1.2392
2.6389
e
-
007
6.7866
e
-
005
1.2915
1.2915
1.9565
e
-
007
5.8436
e
-
005
1.3439
1.3439
4.0858
e
-
007
8.4446
e
-
005
1.3963
1.3964
2.5674
e
-
007
6.6940
e
-
005
1.4486
1.4487
2.3073
e
-
007
6.3459
e
-
005
1.5010
1.5011
5.9533
e
-
007
1.0193
e
-
004
1.5533
1.5534
2.0944
e
-
007
6.0460
e
-
005Case 8.
In this case, we discussed the computational complexity of GA-PS used for L-shape arrays and PM with parallel-shape array. The PM required totalO(3×M×T×K) computations where M is the total number of sensors, K is the total number of sources, and T is the number of snapshots which is 200 [8]. On the other hand, the major computations involved in GA-PS with 2-L-shape array include the total number of multiplication in fitness function O(Q2(3+16×K)) plus the multiplications involved in crossover which is approximately 16×Q2 and the multiplication required for PS which is 16×K. So, the total number of major multiplications involved in GA-PS for 2-L-shape array is O(Q2((3+32×K))+16×K). Similarly, the major computations required for GA-PS with 1-L-shape array are O(Q2((3+20×K))+20×K)), where Q is the number of particles which is 12 in this work. The GA-PS technique with 2-L-shape array is less computationally expensive as compared to PM scheme with parallel-shape array [8]. However, it is little more computationally expensive from GA-PS with 1-L-shape array [15].Case 9.
In this simulation, we compared the RMSE of GA-PS used for 2-L-shape array with its counterpart used for 1-L-shape array [15] as well as with PM for parallel-shape array [8]. In this, a single source is considered which has elevation and azimuth angles 1.0472(rad) and 1.9199(rad), respectively. The values of signal to noise ratio (SNR) are ranging from 5 dB to 25 dB. For GA-PS with L-shape arrays, we have used 10 log × (MSE) where MSE is defined in (16). It is quite obvious from Figure 6, that GA-PS with 2-L-shape array maintained minimum values of RMSE for all values of SNR. The second best RMSE is maintained by another GA-PS technique with 1-L-shape array.Figure 6
Root mean square error versus SNR.Case 10.
In Table13, some general properties of parallel-shape array [8], 1-L-shape array [15] and 2-L-shape array are listed. The main drawbacks of parallel-shape array with PM method include estimation failure in the range of 1.2217 (rad) to 1.5708 (rad), its computational complexity, pair matching problem, and need of more sensors. Some of the drawbacks are covered in 1-L-shape array by using GA-PS technique. However, the main disadvantage of 1-L-shape array [15] is the range limitation of elevation angles beyond π/2, and it also requires more sensors. The 2-L-shape array with GA-PS technique is more effective and requires not only a minimum number of sensors but also covers the range of elevation angle beyond π/2 to π which is of great practical importance in mobile communication [21].Table 13
Comparison among 2-L-shape aray, 1-L-shape array [15] and parallel-shape array [8].
Property
Parallel-shape array [8]
1-L-shape array [15]
2-L-shape array
Scheme used
PM
GA-PS
GA-PS
Elevation and azimuth angles range
(0,π/2), (0, 2π)
(0,π/2), (0, 2π)
(0,π), (0, 2π)
Number of estimated sources
2
2
2
Number of sensors required
33
7
4
Number of Snapshots required
M * T
1
1
Pair matching
Required
Not required
Not required
Failure estimation
1.2217 (rad) to 1.5708 (rad)
No failure
No failure


## 5. Conclusion and Future Work Direction

In this work, five techniques, that is, GA, PS, SA, SA-PS, and GA-PS, are used to jointly estimate amplitude, elevation angle and azimuth angle of far field sources impinging on 2-L-shape array. It has been shown through various simulations that the hybrid GA-PS technique is the best technique among them. It was also found that the GA-PS technique used for 2-L shape array performed well in terms of estimation accuracy, convergence rate, RMSE, number of sensors required as compared to 1-L-shape array and PM for parallel-shape array. Hence, GA-PS approach with 2-L-shape array is more practical and less expensive in terms of hardware budget for wireless communications.In future, one can look into same hybrid approach for null steering and sidelobe reduction as well.



---

*Source: 593247-2013-05-29.xml*