Optimal Location of the Access Points for MIMO-UWB Systems

Featured Application: Provide a secure, fast, and adaptable wireless transmission network for manufacturing systems. Abstract: A multiple-input and multiple-output ultra-wideband (MIMO-UWB) system provides a higher data rate. However, the multipath effect of the intersymbol interference (ISI) increases the bit error rate (BER) and outage probability of the MIMO-UWB system. For this paper, the authors applied the real orthogonal design (ROD) to an MIMO-UWB system to improve the efﬁciency of that system. A ray-tracing technique and an inverse fast Fourier transform were used to get the impulse response of the indoor environment. In addition, a rake receiver was used to increase the strength of the received signal to minimize the multipath effect. For this paper, two cases of an indoor wireless MIMO-UWB system were studied: case (A) used different antenna arrays, whereas case (B) placed antenna arrays in different locations to ﬁnd the best position of the transmitter. In case (A), three different shapes of antenna arrays, namely L-shape, circular-shape, and Y-shape, were used for the transmitter and receiver. The BER performance for these arrays in the UWB frequency of 3.1–10.6 GHz was examined. Numerical results showed that the outage probability of the circular array was better than that of the other two arrays. In case (B), the transmitter used was an array with two antenna elements. The optimal location for the transmitter was found by using both asynchronous particle swarm optimization (APSO) and self-adaptive dynamic differential evolution (SADDE). The numerical results indicated that the performance of APSO was better than that of SADDE.

In an IVWSN, the vehicle control systems and energy budgets of the sensor nodes in an in-car channel contain a large number of metal reflectors because of the short distance. Specifically, the distance between the transmitter and the receiver in an in-car channel is less than 10 m.
Smart antenna technology allows the device to collect more energy to improve the signal-to-noise ratio (SNR) at the receiver, while the multiple-input and multiple-output (MIMO) system has the advantage of high capacity gain [7][8][9][10][11]. MIMO systems are well-suited for working in diverse systems. In rich multipath channels, the probability of having more than one independent path increases. The properties of an MIMO channel, especially its capacity, are highly dependent on the spatial structure of the channel. Different diversity techniques have been applied to communication systems to improve their quality of communication [12][13][14][15].
UWB is mainly applied to indoor environments characterized by the multipath effect. The UWB technique has a plurality of diversity techniques and has almost exclusively been applied to MIMO to improve its quality of communication. The MIMO-UWB system has the potential for short-range, high-speed transmission requirements, and is especially suitable for complex indoor environments. It is a wireless solution that performs better in the current system design. MIMO-UWB has become a compelling candidate for extremely high data-rate communications [16][17][18][19].
This paper applies the real orthogonal design (ROD) [20,21], a diversity technique, in an MIMO-UWB system. Using an ROD to increase the transmit diversity scheme, the transmitter does not require any feedback from receivers. The computation complexity of the ROD is equal to maximal ratio receive combining.
Applying an ROD in an MIMO-UWB system addresses the market demands for quality and efficiency without requiring a significant redesign of existing UWB systems. The bit error rate (BER) performance is better when applying a high-order ROD. However, the intersymbol interference (ISI) effect is proportional to the size of the ROD. The multipath effect of the ISI, which increases the BER and outage probability, degrades the performance of an MIMO-UWB system. The rake receiver combines signals arriving on multiple propagation paths, and the tapped delay line (TDL) is used in this matched filter [22][23][24]. Therefore, a rake receiver is deployed to not only suppress the ISI multipath effect, but also utilize the energy present in multiple propagation paths.
Ray-tracing techniques are used to calculate the channel state information for different antenna radiation patterns [25]. The process can be carried out through simulation to overcome these drawbacks as long as the accuracy of the simulation tools can be proven. It conceptually assumes that a pincushion of rays is sent out from the transmitter, and is traced over the environments until either the ray has reached the receiver or the ray has lost enough power and its contribution to the received signal is negligible. The transmitting ray will be transmitted, reflected, or diffracted to the receiver. The receiver collects a narrowband signal around a given frequency range and converts it to the complex equivalent baseband signals.
For this paper, the authors used three different shapes of 8 × 8 antenna arrays (i.e., L-shape, circular-shape, and Y-shape arrays) for the transmitter and receiver. The BER performances for these arrays in MIMO-UWB communication systems are presented. To the best of the authors' knowledge, there is still no published work using asynchronous particle swarm optimization (APSO) and self-adaptive dynamic differential evolution (SADDE) to search for transmission locations to minimize the outage probability of indoor wireless communication channels. In this paper, two different optimization algorithms are proposed and compared to find the optimal transmitter location for the MIMO-UWB communication system [26][27][28]. The outage probability of the MIMO-UWB systems is investigated. The MIMO-UWB system and its bit error rate calculation formula are described in Section 2.1. APSO and the SADDE evolution algorithms are described in Section 2.2. The results of applying different shape of antenna arrays and finding an optimal transmitter location by applying APSO and SADDE are shown and discussed in Section 3.

System Description
UWB combined with MIMO is a viable way to achieve high data rates in wireless communication. The authors focused on the application of RODs, because they are key to achieving low BER and cost reduction in MIMO-UWB systems. A plurality of complex orthogonal designs for narrowband is available in the literature [20,24]. Figure 1 shows the basic idea of a binary pulse amplitude modulation (BPAM) for a UWB system consisting of N T transmitting and N R receiving antennas with N F rake fingers. Each information symbol is being conveyed through monocycle pulses with a period of T d and with different polarities. A frame is associated with a set of slots. Multi-antenna systems provide a diversity order without any feedback signal from the receiver to the transmitter. Diversity techniques are used to mitigate degradation in the error performance of the radio channel without using additional radiated power or spectral bandwidth [29]. The transmit symbols are applied to the ROD [20,21], where each symbol with a slot time Ts seconds is pulse-shaped by the second order of Gaussian waveform in the transmitters.
communication. The authors focused on the application of RODs, because they are key to achieving low BER and cost reduction in MIMO-UWB systems. A plurality of complex orthogonal designs for narrowband is available in the literature [20,24]. Figure 1 shows the basic idea of a binary pulse amplitude modulation (BPAM) for a UWB system consisting of NT transmitting and NR receiving antennas with NF rake fingers. Each information symbol is being conveyed through monocycle pulses with a period of Td and with different polarities. A frame is associated with a set of slots. Multi-antenna systems provide a diversity order without any feedback signal from the receiver to the transmitter. Diversity techniques are used to mitigate degradation in the error performance of the radio channel without using additional radiated power or spectral bandwidth [29]. The transmit symbols are applied to the ROD [20,21], where each symbol with a slot time Ts seconds is pulse-shaped by the second order of Gaussian waveform in the transmitters.
The received signal Y(t) passes through a matching filter whose impulse response is hij(t) and the sampled signal in the first frame is divided into two branches. These signals are fed to the rake receiver and combined with the MRC algorithm, producing signals r1(t) and r2(t) respectively. Finally, the transmitted symbols S1(t) and S2(t) can be independently decoded from the output r(t). The same procedure is repeated for the next frame. The signals Ŝ(t) are summed up to collect all the diversities provided by multiple receive antennas, multi-time slot, and multipath. The figure shows that the inherent diversity of UWB signals in the temporal domain is collected by the rake receiver, while the diversity in the spatial domain is harvested by the use of the space-time code (STC). It is beneficial to employ an ROD code for collecting more diversities. An ROD of size n is an n by n orthogonal square matrix with entries in the indeterminate ±S; the rows of the matrix represent antennas, whereas its columns represent transmission time slots [24,30]. The ROD of size 2 is as shown below: S1 and S2 are two information symbols transmitted simultaneously from two transmitted antennas (NT = 2). Then, the authors attempted to detect the information symbols S1 and S2 based on the consecutive time slot 1 and time slot 2 (Nslot = 2). Finally, they combined all the spatial diversities provided by the receiving antennas (NR = 2) and the signal components that propagated through the channel by different paths. The received signal Y(t) passes through a matching filter whose impulse response is h ij (t) and the sampled signal in the first frame is divided into two branches. These signals are fed to the rake receiver and combined with the MRC algorithm, producing signals r 1 (t) and r 2 (t) respectively. Finally, the transmitted symbols S 1 (t) and S 2 (t) can be independently decoded from the output r(t). The same procedure is repeated for the next frame. The signalsŜ(t) are summed up to collect all the diversities provided by multiple receive antennas, multi-time slot, and multipath. The figure shows that the inherent diversity of UWB signals in the temporal domain is collected by the rake receiver, while the diversity in the spatial domain is harvested by the use of the space-time code (STC). It is beneficial to employ an ROD code for collecting more diversities.
An ROD of size n is an n by n orthogonal square matrix with entries in the indeterminate ±S; the rows of the matrix represent antennas, whereas its columns represent transmission time slots [24,30]. The ROD of size 2 is as shown below: S 1 and S 2 are two information symbols transmitted simultaneously from two transmitted antennas (N T = 2). Then, the authors attempted to detect the information symbols S 1 and S 2 based on the consecutive time slot 1 and time slot 2 (N slot = 2). Finally, they combined all the spatial diversities provided by the receiving antennas (N R = 2) and the signal components that propagated through the channel by different paths. The ROD codes were developed, and the decoding algorithm was given. Based on the Hurwitz-Radon theorem, the full rate can be achieved by using RODs of size 8, that is, with eight transmitted antennas [20,24]. The ROD of size 8 is as shown below: For this paper, an impulse radio with a BPAM was applied. Time-modulated UWB communication is based on the discontinuous emission of a very short waveform as a train of pulses. Each pulse has ultra-wide spectral requirements in the frequency domain. Because the pulse propagates well in the radio channel, it does not require additional carrier modulation. This radio concept is referred to as impulse radio.
The transmitted PAM signals are as follows: where S j is the transmitted signal of transmitter j, T d is the transmitting signal duration, and the value of T d is usually much larger than the pulse duration. E s is the average transmitted energy per symbol. In a BPAM, each symbol takes only two values S = {+1,−1} and is assumed to be independent and identically distributed (i.i.d.). The second derivative of the Gaussian waveform is used here as the transmit waveform p(t) at the nanosecond scale. The second derivative of Gaussian waveform p(t) can be described by the following expression: where σ is the standard deviation of the Gaussian wave The channel state information (CSI) is available at the receiver, which is the channel impulse response between a receiving antenna and a transmitting antenna. A ray-tracing channel model was developed to calculate the channel state information for the UWB-MIMO system.
The ray-tracing technique can deal with radio wave propagation in any environment. A raytracing approach based on frequency domain was applied to the deterministic simulation for the indoor channel. A pincushion of rays was sent out from the transmitter and traced over a given environment until either it had reached the receiver or it had lost enough power, and its contribution to the received signal was negligible. During transmission, the ray was transmitted, reflected, or diffracted before reaching the receiver. The receiver collected a narrowband signal around the given frequency range and converted the signal to a complex equivalent baseband signal. The goal of channel modeling was to determine the impulse response for any given transmitter-receiver location in the system.
where h ij is the channel impulse response between a receiving antenna i and a transmitting antenna j, L is the path index, β ij,L and τ ij,L are the amplitude and the time delay of L-th path between receiving antenna i and transmitting antenna j. N L is the number of paths. The impulse response with input transmitted signal S(t) and received signal Y(t) is described by convolution. The corresponding received signals are as follows: where * denotes the convolution operator.
After receiving the signals through the channel, they are processed with a correlator in receivers with rake finger N F . With rake receivers, the ability to improve the received signal can be increased in multi-path fading channels. The different signal paths propagated by different paths are combined by rake receivers [31][32][33]. These programs produce time diversity. Combining different signal components increases the signal-to-noise ratio (SNR) and also improves the link performance. The rake receiver is designed to compensate the UWB channel's multipath. It is assumed that the receiver noise n(t) at the corresponding time instant, including the effect of correlation, is Gaussian white noise with zero mean and variance N 0 /2. The probability of BER can be represented as [24]: where α ij,f is the corresponding receiving amplitude for Sj at receiver antenna i with rake finger f, and erfc is the complementary error function. Making the generalization for the case with the ROD of size 8 and rake receiver, the BER formula is as follows: Alternatively, the outage probability can also be used for evaluating the system. The bit error rate (BER) is an essential factor when calculating the outage probability in an MIMO-UWB system. The receiver is defined as the "outage point" as its BER > 10 −6 . The definition of the outage probability is as follows [25,34]

Evolution Algorithms
The following two algorithms were used to minimize the value of the objective function (OF), which is defined in Equation (7).

Self-Adaptive Dynamic Differential Evolution
The SADDE is based on the DDE scheme with the ability to automatically adjusting scaling factors without increasing the time complexity [28,35]. The SADDE algorithm starting from the initial population consists of a randomly generated set of individual coordinates that represent each location of the transmitter antenna. The SADDE algorithm contains the following steps: Step 1: According to the problem, initialize the population of a set of D-dimensional vectors in the starting population, parameter space is X i : i = 1, 2, · · · , M p , where D is the number of parameters to be optimized, and M p is the population size.
Step 2: Evaluate the objective function (7) for each individual in the population.
Step 3: The mutation operation is performed by an arithmetic combination of individuals. Each trial vector V i is generated from the parent's generation parameter vector X i according to the following equation: where α g i and β g i are the scaling factors associated with the vector differences (X g best − X g i ) and (X g m − X g n ), respectively. The disturbance vector V of the mutation mechanism consists of the parameter vector X g i , the best particle X g best , and two randomly selected vectors, α g i and β g i , which are adjusted automatically.
Control parameters in SADDE continue to evolve from generation to generation. The new vector is generated by using the evolution values of the control parameters. During the selection process, these new carriers have a better chance of surviving and passing improved control parameters to the next generation. Each generation's control parameters for each individual here are self-adjusting according to the following scheme: where α l and β l are the lower limits of α and β, their values being set at 0.1. α u and β u are the upper limits of α and β, and their values are set at 0.9 [28,36]; rand 1 , rand 2 , rand 3 , and rand 4 are random numbers with values uniformly distributed between 0 and 1.
Step 4: Perform crossover operations to increase the diversity of parameter vectors. The cross vector u i is replaced by the current vector X g i or trial vector V i based on the probability of crossover Crs g i . It can be expressed as follows: where κ(j) is the random number generated uniformly between 0 and 1. Crs g i is the crossover probability, Crs g i ∈ (0, 1). Rand 5 and rand 6 are random numbers with the values uniformly distributed between 0 and 1. Crs g i is the crossover probability, Crs g i ∈ (0, 1). The bigger the value of Crs g i , the more rapidly the termination criteria can be reached. However, a smaller value of Crs g i increases the diversity of the parameter vectors.
Step 5: The selection operation is used to generate offspring. Compare the parent vector X g i with the cross vector u g+1 i and select the vector with the smaller objective function value as the next generation member. The selection operation is shown below: SADDE is a self-adaptive evolutionary algorithm. Specifically, if the offspring has a better objective function value than its parent, then this offspring individual can be used to replace each parent. If there are individuals in the offspring that meet the termination criteria, stop this process and get the best individual; otherwise, return to step 2 to continue the entire cycle.

Asynchronous Particle Swarm Optimization
Asynchronous particle swarm optimization (APSO), proposed by Clerc, is an alternative to the particle swarm optimization (PSO) algorithm [36]. APSO uses a different velocity update equation, which utilizes a parameter called constriction factor. The purpose of the constriction factor is to improve the convergence. As APSO updates the global best particle position as soon as one is found that is less than the value of the objective function. APSO will converge more quickly than PSO. For this paper, APSO contains six steps in minimizing the value of the objective function, as shown below.
The number of particles X i (i = 1, 2, . . . , N p ) with D dimensional vectors is chosen, where D is the number of parameters to be optimized, and N p is the population size: Step 1: Randomly set the initial particle position based on the search range.
Step 2: Evaluate the value of the objective function for each particle and record the particle information, that is, the particle position which has a value less than the current value of the objective function, and the new and smaller value of the objective function.
Step 3: Update the global best particle position x gb and best particle position x pb for the gth generation based on the data recorded in Step 2.
Step 4: Ten percent of x gb will be mutated by the following equation: where c 1 and c 2 are the scaling parameters. Both ϕ 1 and ϕ 2 are random real numbers from zero to one. g (g = 1, 2, 3, ..., g max ) is the current iteration number. The space limit of the environment is between x max and x min . If the objective function of the mutation particle is the best, update the x gb with x mu .
Step 5: Calculate the particle velocity using the following equation: where ij is the ith particle and jth dimension, ∆ = Step 6: Stop the algorithm and output the final data when the termination criteria have been reached.  Figure 2a is the indoor environment layout, and the actual office was simulated for this paper. In this simulation, the dimension of the office was 5 m long, 10 m wide, and 4.5 m high. A partition, a wooden desk, and an iron cabinet were put inside the office. All the floors, walls, and ceilings were made of 20 cm thick concrete. The transmitter was placed at the center of the office at axis (2.5, 5.0, 2.5) m, and 105 receivers were uniformly distributed in the office, as shown in Figure 2b The height of the receivers was 1.5 m above the ground. The number of fingers for the rake receiver was 4 (N F = 4).  Figure 2a is the indoor environment layout, and the actual office was simulated for this paper. In this simulation, the dimension of the office was 5 m long, 10 m wide, and 4.5 m high. A partition, a wooden desk, and an iron cabinet were put inside the office. All the floors, walls, and ceilings were made of 20 cm thick concrete. The transmitter was placed at the center of the office at axis (2.5, 5.0, 2.5) m, and 105 receivers were uniformly distributed in the office, as shown in Figure 2b The height of the receivers was 1.5 m above the ground. The number of fingers for the rake receiver was 4 (NF = 4). In a real environment, the dielectric constant and conductivity change with the frequency. Therefore, the frequency dependence for the dielectric constant and conductivity was also considered in the simulation. Since there were many receivers, the average BER and the outage probability were used for evaluating the system.

Case (A): Antenna Array
Both the transmitter and the receiver used an eight-antenna array arranged in different shapes (i.e., L-shape, circular-shape, and Y-shape), as shown in Figure 3. The transmitter was located at the center of the office at 2 m height. In a real environment, the dielectric constant and conductivity change with the frequency. Therefore, the frequency dependence for the dielectric constant and conductivity was also considered in the simulation. Since there were many receivers, the average BER and the outage probability were used for evaluating the system.

Case (A): Antenna Array
Both the transmitter and the receiver used an eight-antenna array arranged in different shapes (i.e., L-shape, circular-shape, and Y-shape), as shown in Figure 3. The transmitter was located at the center of the office at 2 m height. Appl. Sci. 2018, 8, x FOR PEER REVIEW 9 of 13 (a) (b) (c) The receivers were evenly distributed in the office at 1.5 m height. The performance for each shape of the antenna array in UWB frequency of 3.1-10.6 GHz was evaluated. The average BER versus SNR, the ratio of transmitting signal power to the noise power, for the 105 receivers are shown in Figure 4a, whereas the outage probability versus SNR for these receivers is shown in Figure 4b. In both figures, for most of the SNR values, the circular array has the smallest average BER and outage probability. For example, in Figure 4a, when the SNR is 68 dB, the BER values for the circular array, the L array, and the Y array are 0.18%, 0.61%, and 0.87%, respectively. In Figure 4b, when the SNR is 68 dB, the outage probability values for the circular array, the L array, and the Y array are 0.95%, 1.90%, and 1.90%, respectively. It is clear that the BER and the outage probability of the circular-shape array are better than those of the other two arrays in most of the circumstances.

Case (B): Optimal Location by Evolution Algorithm
Case (B) used SADDE and APSO to find the optimal transmitter location. To save computing time, both the transmitter and the receiver used an array of two antennas. The height of the transmitter was 2 m, and that of the receiver was 1.5 m. The outage probability was calculated as defined in Equation (7).
Initially, the optimal transmitter location was set at axis (2.5, 5.0, 2.0) m, the center of the office. The parameter values are shown in Table 1. As shown in Figure 5, when the SNR was 56 dB, there were 11 outage points. To improve the communication quality and reduce the outage probability [26], SADDE and APSO were used to find the optimal location of the transmitter. In Figure 6a, after using the SADDE to adjust the optimal transmitter location to axis (2.81, 5.08, 2.0) m, the outage points dropped from 11 to 10. After using the APSO to adjust the optimal transmitter location to axis (2.61, 3.95, 2.0) m, the outage points dropped from 11 to 7, as shown in Figure 6b. When the SNR increases The receivers were evenly distributed in the office at 1.5 m height. The performance for each shape of the antenna array in UWB frequency of 3.1-10.6 GHz was evaluated. The average BER versus SNR, the ratio of transmitting signal power to the noise power, for the 105 receivers are shown in Figure 4a, whereas the outage probability versus SNR for these receivers is shown in Figure 4b. In both figures, for most of the SNR values, the circular array has the smallest average BER and outage probability. For example, in Figure 4a, when the SNR is 68 dB, the BER values for the circular array, the L array, and the Y array are 0.18%, 0.61%, and 0.87%, respectively. In Figure 4b, when the SNR is 68 dB, the outage probability values for the circular array, the L array, and the Y array are 0.95%, 1.90%, and 1.90%, respectively. It is clear that the BER and the outage probability of the circular-shape array are better than those of the other two arrays in most of the circumstances. The receivers were evenly distributed in the office at 1.5 m height. The performance for each shape of the antenna array in UWB frequency of 3.1-10.6 GHz was evaluated. The average BER versus SNR, the ratio of transmitting signal power to the noise power, for the 105 receivers are shown in Figure 4a, whereas the outage probability versus SNR for these receivers is shown in Figure 4b. In both figures, for most of the SNR values, the circular array has the smallest average BER and outage probability. For example, in Figure 4a, when the SNR is 68 dB, the BER values for the circular array, the L array, and the Y array are 0.18%, 0.61%, and 0.87%, respectively. In Figure 4b, when the SNR is 68 dB, the outage probability values for the circular array, the L array, and the Y array are 0.95%, 1.90%, and 1.90%, respectively. It is clear that the BER and the outage probability of the circular-shape array are better than those of the other two arrays in most of the circumstances.

Case (B): Optimal Location by Evolution Algorithm
Case (B) used SADDE and APSO to find the optimal transmitter location. To save computing time, both the transmitter and the receiver used an array of two antennas. The height of the transmitter was 2 m, and that of the receiver was 1.5 m. The outage probability was calculated as defined in Equation (7).
Initially, the optimal transmitter location was set at axis (2.5, 5.0, 2.0) m, the center of the office. The parameter values are shown in Table 1. As shown in Figure 5, when the SNR was 56 dB, there were 11 outage points. To improve the communication quality and reduce the outage probability [26], SADDE and APSO were used to find the optimal location of the transmitter. In Figure 6a, after using the SADDE to adjust the optimal transmitter location to axis (2.81, 5.08, 2.0) m, the outage points dropped from 11 to 10. After using the APSO to adjust the optimal transmitter location to axis (2.61, 3.95, 2.0) m, the outage points dropped from 11 to 7, as shown in Figure 6b. When the SNR increases

Case (B): Optimal Location by Evolution Algorithm
Case (B) used SADDE and APSO to find the optimal transmitter location. To save computing time, both the transmitter and the receiver used an array of two antennas. The height of the transmitter was 2 m, and that of the receiver was 1.5 m. The outage probability was calculated as defined in Equation (7).
Initially, the optimal transmitter location was set at axis (2.5, 5.0, 2.0) m, the center of the office. The parameter values are shown in Table 1. As shown in Figure 5, when the SNR was 56 dB, there were 11 outage points. To improve the communication quality and reduce the outage probability [26], SADDE and APSO were used to find the optimal location of the transmitter. In Figure 6a, after using the SADDE to adjust the optimal transmitter location to axis (2.81, 5.08, 2.0) m, the outage points dropped from 11 to 10. After using the APSO to adjust the optimal transmitter location to axis (2.61, 3.95, 2.0) m, the outage points dropped from 11 to 7, as shown in Figure 6b. When the SNR increases from 56 dB to 64 dB, the number of the outage points for SADDE was reduced from 10 to 1, whereas the number for APSO dropped from 7 to 0. This result demonstrated that APSO was better than SADDE in finding an optimal transmitter location, as shown in Table 2. from 56 dB to 64 dB, the number of the outage points for SADDE was reduced from 10 to 1, whereas the number for APSO dropped from 7 to 0. This result demonstrated that APSO was better than SADDE in finding an optimal transmitter location, as shown in Table 2.

Parameter
Value search target (m) = (x, y, z) location of the transmitter (x, y, 2.0) fitness function outage probability of the receivers frequency range 3.1-10.6 GHz SNR 40 dB to 80 dB

Conclusions
For this paper, the authors applied an ROD in an MIMO-UWB system to achieve low BER and to reduce cost. Moreover, a rake receiver was used to increase the strength of the received signals to minimize the multipath effect. Three different shapes of antenna arrays, namely L-shape, circularshape, and Y-shape, were used for both the transmitter and the receiver, and the BER performance and outage probability of each shape of the antenna array were compared. The BER and the outage probability of the circular array were better than those of the other two arrays. Two-antenna element arrays were used for the transmitter and receiver, and asynchronous particle swarm optimization (APSO) and self-adaptive dynamic differential evolution (SADDE) were used to find the optimal location of the transmitter. The results showed that APSO performs better than SADDE.