Efficient Preamble Design Technique for Millimeter-Wave Cellular Systems with Beamforming

The processing time for beam training in millimeter-wave (mmWave) cellular systems can be significantly reduced by a code division multiplexing (CDM)-based technique, where multiple beams are transmitted simultaneously with their corresponding Tx beam IDs (BIDs) in the preamble. However, mmWave cellular systems with CDM-based preambles require a large number of cell IDs (CIDs) and BIDs, and a high computational complexity for CID and BID (CBID) searches. In this paper, a new preamble design technique that can increase the number of CBIDs significantly is proposed, using a preamble sequence constructed by a combination of two Zadoff-Chu (ZC) sequences. An efficient technique for the CBID detection is also described for the proposed preamble. It is shown by simulations using a simple model of an mmWave cellular system that the proposed technique can obtain a significant reduction in the complexity of the CBID detection without a noticeable performance degradation, compared to the previous technique.


Introduction
Mobile traffic has been increasing markedly owing to the growth of smart devices such as smart phones, tablets, and laptops. By 2020, the mobile traffic volume is expected to increase by more than 1000-fold because of the increase in multimedia and IoT services. In order to support such a considerable increase in the mobile traffic, the 30-300 GHz millimeter-wave (mmWave) frequency band supporting a wide bandwidth is considered as a possible candidate for 5G mobile communication systems with antenna arrays for directional beamforming at both the base station (BS) and the mobile station (MS) [1,2]. The array gain obtained with the directional beamforming extends the coverage of the mmWave systems that suffer from a high path loss owing to a high carrier frequency [3]. Because of the small wavelength of the mmWaves, antenna arrays can be easily installed in the MSs.
Owing to significant costs in the RF chain and the power consumption, analog beamforming is currently preferred over digital beamforming for mmWave communication systems [4]. In mmWave communication systems with analog beamforming, switched beamforming techniques with a set of predefined angles are usually used for transmit-receive (Tx-Rx) beamforming. In these systems, a maximum array gain is obtained when the Tx and Rx beams are perfectly aligned. A Tx-Rx beamforming technique using the 60 GHz unlicensed spectrum has already been standardized for the wireless LAN (IEEE 802.11ad and IEEE 802.15c) to provide a multigigabit-per-second data rate [5][6][7]. In the IEEE 802.11ad standard, a two-level scheme composed of a sector level sweep (SLS) phase and a beam refinement protocol (BRP) phase is adopted for beam training [5]. During the SLS phase, a sector level beam sweep is performed, where the best transmit sector beam is selected with a quasi-omni antenna pattern at the receiver. During the BRP phase, a refined beam selection is carried out with the receive beam at a high resolution to compensate for the imperfections in the quasi-omni antenna pattern. However, in an mmWave cellular system, it is not possible to use the quasi-omni antenna pattern on the receive side because it requires a large transmission range unlike the wireless LANs [8]. Directional beamforming antennas are also necessary at both the BS and the MS to maintain a link budget in the beam-training period. The beam training technique in the IEEE 802.11ad standard cannot be used for cellular communication systems because it is mainly targeted for indoor communication, without a handover capability. The beam training protocol in cellular systems requires a considerable training time and network resources because it needs to select a serving BS with the best beam pair after checking the link qualities of all the possible beam pairs for all the neighbor BSs.
The best beam pair in the mmWave cellular system can be found by using a time-division beam switching (TDBS) during the beam-training period. In the TDBS scheme, single Tx beams are individually transmitted from the BS until all the Tx beams are transmitted. A Tx beam ID (BID) is distinguished by the location of the time slot where the corresponding Tx beam is transmitted. The Rx beam sweep is performed at the MS for each Tx beam for measuring the SNR of every Tx-Rx beam pair. The measurement of the SNRs for all the possible Tx-Rx beam pairs must be performed for all the neighbor BSs to select a serving BS with the best beam pair. The processing time required for beam training in the TDBS increases proportionally to the product of the number of Tx beams, the number of Rx beams, and the number of neighbor BSs. This considerable processing time will create a significant overhead for a moving MS because beam training should be performed periodically for a possible handover or for beam tracking.
In order to reduce the processing time for beam training in mmWave cellular systems, a code division multiplexing (CDM)-based technique was proposed in [9], where multiple beams with their corresponding Tx BIDs are transmitted simultaneously during a beam-training or a beam-tracking period. In [9], the cell ID (CID) as well as the corresponding BID is transmitted in the preamble so that the BID can be detected in a multicell environment. The BID is mapped to a single Zadoff-Chu (ZC) sequence in association with its CID because the ZC sequence is suitable for preamble designs in OFDM systems owing to its good correlation property and low peak-to-average power ratio (PAPR) [10][11][12][13]. However, the computational complexity required for the CID and BID (CBID) detection in [9] increases significantly because the correlation operation required for the CBID detection increases proportionally to the product of the number of CIDs and BIDs. In addition, the number of CIDs in mmWave cellular systems is expected to increase significantly because the cells should be deployed in a denser manner than in the existing systems. Therefore, a new preamble design technique is needed for mmWave cellular systems that can provide a larger number of CIDs and permit a CBID detection with a lower computational complexity.
In this paper, a CDM-based preamble design technique for the mmWave cellular communication systems is proposed using a preamble sequence constructed by a combination of two ZC sequences. It is shown that the proposed preamble design technique can increase the number of CBIDs significantly, compared to the previous technique generated using a single ZC sequence. It is also shown that the proposed technique can reduce the computational complexity in the CBID detection without a noticeable performance degradation, using a specific mapping rule and a parameter selection method. Because of the special structure used in preamble generation, the proposed preamble is shown to have a low PAPR characteristic similar to the one obtained using a single ZC sequence.
The rest of this paper is organized as follows: in Section 2, a preamble design technique using two ZC sequences is proposed for OFDM-based mmWave cellular systems with beamforming. In Section 3, an efficient technique for the CBID detection is described for the proposed preamble. Parameter selection methods are also discussed to avoid the performance degradation caused by a preamble sequence constructed using a combination of two ZC sequences. The performance of the proposed technique is evaluated by computer simulation in Section 4. The conclusions are drawn in Section 5.

Proposed Preamble Design Technique
The two ZC sequences used to generate a preamble sequence in the proposed technique are given by [14]: Z r,d s pmq " Z r s pm`d s q " e´j πr s pm`d s qpm`d s`1 q{L , GCD pr s , Lq " 1 (1) where m P t0, 1,¨¨¨, L´1u , s P t0, 1u. L, s, and GCD denote a prime number representing the length of the ZC sequence, a parameter indicating one of the two different ZC sequences, and the greatest common divisor, respectively. In addition, r s P t1,¨¨¨, L´1u and d s P t0, 1,¨¨¨, L´1u denote the root index and the cyclic shift of the s-th ZC sequence, respectively. The number of preamble sequences that can be generated by a combination of the root indices and the cyclic shifts of the two ZC sequences, is given by pL´1q 2 L 2 . Therefore, the number of preamble sequences that can be generated by combining two ZC sequences is considerably larger than the ones generated by a single ZC sequence. For example, the number of available preamble sequences is 3660 for a single ZC sequence with a length of 61 and 864,900 for two ZC sequences with lengths of 31 each. Therefore, the number of preamble sequences that can be generated by combining two ZC sequences with length of 31 is 236 times greater than the ones generated by a single ZC sequence with a length of 61. The proposed preamble is generated by allocating two ZC sequences with lengths, L, in the frequency domain, as follows: where r X r,d s ppk`N 0 q %N q " r X r,d s pkq: Here, N denotes the number of subcarriers (FFT size) in an OFDM symbol. Figure 1 shows the concept of the preamble generation in the proposed technique. As shown in Figure 1a,b, r X r,d 0 pkq and r X r,d 1 pkq can be viewed as the two repetitions of Z r,d 0 pmq and Z r,d 1 pmq, respectively. Here, null subcarriers generated by the difference between the FFT size and the sequence length are ignored. As given in Equation (2), X r,d pkq is the sum of r X r,d 0 pkq and r X r,d 1 pkq multiplied by e´j 2πk{N . The proposed preamble in the time domain can be obtained by performing an IFFT of X r,d pkq, as follows:
The preamble in Equation (3) can be rewritten as: where q P t0, 1u.   x n [14]. The upper bound for the sum of the correlation values of the two sequences, , r d C  , that will be used for the sequence detection in Section 3, can be derived using the triangular inequality as follows: becomes 1 / 2 as given in Property 1 of Equation (5). Here, % denotes the modulo operation. The parameter, q, is used to distinguish the even samples (q " 0) from the odd samples (q " 1). As shown in Figure 1c, N 0 -point FFTs of the even and odd samples of x r,d pnq are given by r X r,d 0´p kq %N 0¯a nd r X r,d 1´p kq %N 0¯, respectively. In other words, even and odd samples of x r,d pnq consist of Z r,d 0 pmq and Z r,d 1 pmq, respectively. Because of this property, the proposed preamble can have a low PAPR characteristic similar to the one obtained using a single ZC sequence.
The correlation property for each sequence (composed of either odd or even samples) can be easily obtained from x r,d pnq [14]. The upper bound for the sum of the correlation values of the two sequences, r C r,d , that will be used for the sequence detection in Section 3, can be derived using the triangular inequality as follows: ere, {r 1 s , d 1 s } and {r 2 s , d 2 s } are used to distinguish the parameters of {r s , d s } for the two different preamble sequences, respectively. δp¨q denotes the Kronecker delta function. From Equation (5), it can be seen that if r s " 0, r C r,d becomes one when d s " 0 and becomes zero when d s ą 0. In this case, we can obtain an ideal correlation property. Further, if r s ą 0, r C r,d becomes 1{ ? L for all values of d s . In this case, a relatively good correlation property can be obtained. However, when r s " 0 and C r,d becomes 1{2 as given in Property 1 of Equation (5). When ? L`1˘) as given in Property 3 of Equation (5). These properties (Property 1 and Property 3) can lead to a performance degradation in the sequence detection. In the proposed technique, we select a set of root indices and cyclic shifts such that the case satisfying the condition in Property 1 and Property 3 does not occur or is minimized. The method for selecting the parameters will be described in the next section.

Efficient CBID Detection Technique for mmWave Cellular Systems
In this section, an efficient detection technique that can reduce the computational complexity of the CBID search using the preamble sequence described in the previous section is proposed for mmWave cellular systems. Figure 2 shows an example of an mmWave cellular system with two BSs and one MS. It is assumed that the switched beamforming technique with a set of predefined angles is used at both the BS and the MS. It is also assumed that the BS has multiple antenna arrays, whereas the MS has only one antenna array. That is, multiple beams are transmitted simultaneously from the BS, while a single switched beam is used at the MS. The sidelobes of each beam are omitted in this figure for simplicity. Details on the parameters for the preamble generation, beamforming, and the CBID information will be described in Section 4.
as given in Property 3 of Equation (5). These properties (Property 1 and Property 3) can lead to a performance degradation in the sequence detection. In the proposed technique, we select a set of root indices and cyclic shifts such that the case satisfying the condition in Property 1 and Property 3 does not occur or is minimized. The method for selecting the parameters will be described in the next section.

Efficient CBID Detection Technique for mmWave Cellular Systems
In this section, an efficient detection technique that can reduce the computational complexity of the CBID search using the preamble sequence described in the previous section is proposed for mmWave cellular systems. Figure 2 shows an example of an mmWave cellular system with two BSs and one MS. It is assumed that the switched beamforming technique with a set of predefined angles is used at both the BS and the MS. It is also assumed that the BS has multiple antenna arrays, whereas the MS has only one antenna array. That is, multiple beams are transmitted simultaneously from the BS, while a single switched beam is used at the MS. The sidelobes of each beam are omitted in this figure for simplicity. Details on the parameters for the preamble generation, beamforming, and the CBID information will be described in Section 4.
respectively. In Figure 3, 'SP' represents a synchronization preamble from which the symbol timing offset (STO) and carrier frequency offset (CFO) are estimated. The beam ID associated with the cell ID is acquired from the 'BP' (BID preamble). In this paper, we assume that the synchronization procedure has been completed with the SP using the same approach as in an LTE system with synchronization signals or in a Mobile WiMAX with preambles [12,15]. In this paper, we will focus only on the design of the BP. As shown in Figure 3,    Figure 3 shows only the preamble segment, ignoring the data transmission segment in the frame structure. In this figure, N B and N M , denote the number of Tx beams (BIDs) and Rx beams, respectively. The variables for the set of BIDs and Rx beam indices are defined as b P t0, 1,¨¨¨, N B´1 u and i P t0, 1,¨¨¨, N M´1 u, respectively. In Figure 3, 'SP' represents a synchronization preamble from which the symbol timing offset (STO) and carrier frequency offset (CFO) are estimated. The beam ID associated with the cell ID is acquired from the 'BP' (BID preamble). In this paper, we assume that the synchronization procedure has been completed with the SP using the same approach as in an LTE system with synchronization signals or in a Mobile WiMAX with preambles [12,15]. In this paper, we will focus only on the design of the BP. As shown in Figure 3, N B beams are simultaneously transmitted from the BS and the Rx beams are swept over in this period. N B beams are transmitted repeatedly, N M times, until one round of the Rx beam sweep is completed. The Rx beam switching takes place for every 2 nd symbol such that the MS can receive a pair of preambles (SP and BP) for each beam. Note that the processing time for the beam training when the CDM-based scheme is used is N B times smaller than in the case of the TDM-based scheme.  If the CDM-based technique is used for multiple beam transmission in cellular communication systems, the BID as well as the corresponding CID needs to be transmitted in the preamble because the BID needs to be detected in a multicell environment. Therefore, the information (both the CID and the BID) assigned to the beam should be transmitted in the BP. A hierarchical preamble design concept is used to carry the BID information as well as the CID information in the BP. The information on the CID ( c ) and BID ( b ) are mapped to the root indices ( 0 r , 1 r ) and the cyclic shifts ( 0 d , 1 d ) given in Equation (1) as follows: Here, g and v denote the group ID (GID) and sequence ID (SID), respectively. The CID ( c ) is expressed by a combination of g and v . C N , g N , and v N denote the number of CIDs, GIDs, and SIDs, respectively. P L denotes the scaling factor for a phase change.

Spatial domain
Time domain If the CDM-based technique is used for multiple beam transmission in cellular communication systems, the BID as well as the corresponding CID needs to be transmitted in the preamble because the BID needs to be detected in a multicell environment. Therefore, the information (both the CID and the BID) assigned to the beam should be transmitted in the BP. A hierarchical preamble design concept is used to carry the BID information as well as the CID information in the BP. The information on the CID (c) and BID (b) are mapped to the root indices (r 0 , r 1 ) and the cyclic shifts (d 0 , d 1 ) given in Equation (1) as follows: where Z g s pmq " e´j πr g s pmqpm`1q{L , S c,b s " e´j πd c,b s pd c,b s`1 q{L Here, g and v denote the group ID (GID) and sequence ID (SID), respectively. The CID (c) is expressed by a combination of g and v. N C , N g , and N v denote the number of CIDs, GIDs, and SIDs, respectively. L P denotes the scaling factor for a phase change. r g s , d c,b s , and O v s denote the root index for the GID (g), the cyclic shift corresponding to the CID (c) and BID (b), and the phase rotation offset for the SID (v), respectively, for the s-th ZC sequence. In Equation (6)  . The mapping rule for the GID is given by: where r g 1 " g`r g{ pr g`1q {2 Here, N r denotes the number of root indices used in the preamble design. An example of Equation (7) when N g " 12 and (N r " 12, N r " 7, N r " 5) is given in Table 1. In this table, the root index pairs corresponding to the GIDs are listed. For example, 12 root indices are used when N r " 12, whereas five root indices are used when N r " 5. Table 1. Example of the GID mapping: (a) N r " 12; (b) N r " 7 ; (c) N r " 5. As given in Equation (6), v is mapped to a pair of phase rotation offsets, , N C can be expressed by w pN r´p w´1q {2q N v . The maximum value of N C is obtained when N r " L´1 and is given by w pL´pw`1q {2q N v . Therefore, the proposed technique can increase the number of available CIDs compared to the previous technique with a single ZC sequence. For example, when the length of a ZC sequence is 1021, the number of available CIDs is 1020 in the previous technique and 6084 in the proposed technique, when L, w, and N v are set to 509, 3, and 4, respectively. The number of CIDs in the proposed technique is six times larger than the one in the previous technique. The values of r d v,b s are selected such that the conditions r are satisfied for two different preambles with pv 1 , b 1 q and pv 2 , b 2 q. Here, 0 ď v 1 , v 2 ď N v´1 and 0 ď b 1 , b 2 ď N B´1 . This condition is required to remove the poor correlation properties (Property 1 and 3) in Equation (5). Note that the correlation property of the proposed preamble can be changed depending on the value of N r . As N r decreases, the number of side peaks satisfying the condition in Property 3 increases. The maximum number of side peaks is given by 2 pw´1q. In the example for the GID mapping rule in Table 1(b/c), the maximum number of side peaks satisfying the condition in Property 3 is two/four because the value of w is two/three. For example, the root index 2/3 appears twice/four times for different GIDs, producing high correlation values (side peaks). Thus, as N r decreases, the number of side peaks producing high correlation values at incorrect positions can increase. However, the computation complexity for the CBID detection is reduced as N r decreases because the number of available root indices decreases. In the CBID detection, it is assumed that N 0 -point FFT is used. However, because the length of the ZC sequence, L, is usually different from the FFT size, N 0 , P v,b s pmq in Equation (6) can be rewritten as: where ρ denotes the ratio of N 0 to L and is a rational number. The non-integer factor, ρ, usually degrades the performance of the CBID detection. In order to compensate the effect of the fractional part in ρ on the CBID detection, we multiply Z c,b s pmq by a polyphase sequence, D v,b s pmq, as follows: s . In the following, Z r,d s , r X r,d s , X r,d , and x r,d , used in Equations (2)-(4), are replaced by r Z c,b s , r X c,b s , X c,b , and x c,b , respectively, to formulate the equations for the CBID detection in mmWave cellular systems.
In addition, there may exist a residual STO even after the synchronization procedure has been completed with the SP. Because the residual STO causes a phase rotation in the frequency domain, the BP needs to be designed to be robust to the residual STO [16]. Next, we analyze the effect of the residual STO on the CBID detection. When the residual STO (χ) exists, the decimated version (even and odd samples) of x c,b ppn´χq %N q with q P t0, 1u is given by: where r m " N 0 {2´tL{2u`m, χ " 2µ l`l U l " tu l | 0 ď u l ď tχ Max {2u´l pχ Max`1 q %2 u l " pχq %2 , χ P t0, 1,¨¨¨, χ Max u Here, l denotes a parameter indicating that χ is an even number (l " 0) or an odd number (l " 1). It can be seen from Equation (1) that s " pq`lq %2 changes depending upon the value of the STO for the same q. When χ is an even number, the root index of the ZC sequence of a decimated version with q " 0 in Equation (10) corresponds to the case, s " 0. However, when χ is an odd number, the root index of a ZC sequence of a decimated version with q " 0 corresponds to the case, s " 1. Hence, an ambiguity in the GID detection may occur because the order of the indices in a root pair corresponding to the GID can be reversed depending upon the value of the STO in the CBID detection. For example, N 0 -point FFTs of the even and odd samples of x c,b ppn´χq %N q are expressed as r X c,b 0´p kq %N 0¯e´j 2πu 0 k{N 0 and r X c,b 1´p kq %N 0¯e´j 2πu 0 k{N 0 , when χ is an even number, but as r X c,b 1´p kq %N 0¯e´j 2πpu 1`1 qk{N 0 and r X c,b 0´p kq %N 0¯e´j 2πu 1 k{N 0 , when χ is an odd number. In order to avoid ambiguity, the value of r d v,b 0 needs to be selected such that it is not equal to r should be selected such that´r Next, we consider an OFDM-based mmWave cellular system with a Tx-Rx beamforming. When the proposed preamble with the CID, c, and the BID, b, is transmitted from the b-th Tx beam of the BS, the received signal at the i-th Rx beam in the MS is given in the frequency domain as follows: Here, it is assumed that the initial STO and CFO synchronizations have been completed and only the residual STO exists. The parameters η, µ, H, G, and W denote the Tx array gain, Rx array gain, channel frequency response, total gain including the Tx and Rx array gains, and the additive white Gaussian noise (AWGN), respectively. The detection of the CID, BID, and the residual STO using the received signal can be performed by correlating the received signal in Equation (11) with the preamble sequence in the frequency domain as follows: Here,ĝ,v,ĉ, andî denote the detected GID, SID, CID, and Rx beam index, respectively.χ,l, andûl, denote the estimated parameters for the residual STO. Y q i´p kq %N 0¯r epresents the frequency domain version of y i p2n`qq. Here, y i p2n`qq denotes the decimated version of the received signal in the time domain with q P t0, 1u. Because r g s P t1, 2,¨¨¨, N r u is always satisfied for the root indices in ! r g 0 , r g 1 ) selected by Equation (7), the IFFT operation required for the CBID detection needs to be performed N r times for each r Y 0 i pmq and r Y 1 i pmq. Therefore, the computational complexity for the CBID detection decreases as N r decreases.

Simulation
In this section, the performance of the proposed preamble design technique is evaluated by computer simulation with a simple setup of an OFDM-based mmWave cellular system, as shown in Figure 2. The center frequency, bandwidth, FFT size, CP size, and subcarrier spacing are set to 28 GHz, 250 MHz, 1024, 128, and 270 kHz, respectively [17]. A Rician channel model consisting of one LOS path and one NLOS path is used for an mmWave channel [3] and is programmed with a 3-dimensional spatial channel model (3D-SCM) [18]. The number of rays in the NLOS path is set to 20, and the K-factor is set to 10 dB in our simulation.
In the simulation, only BS1 is considered for a one-cell environment and both BS1 and BS 2 are considered for a two-cell environment. The cell radius is set to 500 m. In a two-cell environment, it is assumed that the MS with a switched beam is located 400 m away from BS1. The residual STO (χ) is generated uniformly in the range 0-4. A uniform linear array (ULA) with eight antenna elements is used at the BS. A uniform circular array (UCA) with eight antenna elements is used at the MS to cover a 360˝azimuth angle around the MS. The parameters for the preamble generation, beamforming, and the CBID information used in the simulation are listed in Table 2. Table 2. Parameters for the simulation.

Parameters for Preamble Generation
(1020,255,4,12,12)  Different notations are used for the length of ZC sequence, L, to distinguish between the parameters used for each technique. The notation, L, is used for the previous technique [9] and r L is used for the proposed technique. Here, the previous technique refers to the case where the CBID information is mapped to the preamble sequence generated by a single ZC sequence. In the proposed technique, a preamble sequence is generated by a combination of two ZC sequences. Here, L and 2 r L are set to 1021 and 1018, respectively. The number of CIDs, N C , is set to 1020 (L´1) for both the previous and the proposed techniques. The values of L P and O v 0 , O v 1 ( are selected such that r d v,b s satisfies the conditions described in Section 3 and are listed in Table 2. In the simulation, two types of the proposed preamble sequences (Type-1 and Type-2) with different values of N r are considered. The value of N r is set to 255 for Type-1 and 86 for Type-2, respectively. The value of N g is set to 255 for both Type-1 and Type-2. Figure 4 shows the correlation properties of the proposed preamble sequences (Type-1 and Type-2) when CID is 278 and BID is six. It can be seen from Figure 4a that the preamble sequence of Type-1 has a good correlation property because the number of side-peaks is zero. Note that the number of side-peaks (satisfying the condition in Property 3) is given by 2 pw´1q, where w " P N g {N r T . In the case of the Type-2 preamble, there are side-peaks with a value of 0.5 at four different CIDs because the value of 2 pw´1q is four, as can be seen in Figure 4b. The root index pairs corresponding to the CIDs (618, 622, 954, 962) are ({69, 70}, {70, 71}, {68, 70}, {70, 72}), whereas the root index pair corresponding to the CID (278) is {70,70}. A high correlation occurs at these four CIDs with BID = 6 because the same root index, 70, is used.
proposed preamble sequences (Type-1 and Type-2) when CID is 278 and BID is six. It can be seen from Figure 4a that the preamble sequence of Type-1 has a good correlation property because the number of side-peaks is zero. Note that the number of side-peaks (satisfying the condition in Property 3) is given by   In the case of the Type-2 preamble, there are sidepeaks with a value of 0.5 at four different CIDs because the value of   2 1 w  is four, as can be seen in Figure 4b. The root index pairs corresponding to the CIDs (618, 622, 954, 962) are ({69, 70}, {70, 71},  {68, 70}, {70, 72}), whereas the root index pair corresponding to the CID (278) is {70,70}. A high correlation occurs at these four CIDs with BID = 6 because the same root index, 70, is used.  Figure 5 shows the success probability of the CBID detection in one-cell and two-cell environments. CBID detection using the preambles simultaneously transmitted from the BS(s), is performed at the MS using Equations (12)- (14). The success of the CBID detection refers to the case where the CID and BID detected at the MS match the ones used in the BS for the beam transmission.
Simulations are performed for three different cases (1 × 1, 8 × 1, 8 × 8), where the first and second terms denote the number of antenna elements in the BS and MS, respectively. When the number of antenna elements is equal to one, an omnidirectional antenna is used. From this figure, it can be seen that an array gain of 9 dB is obtained when the number of antenna elements increases from one to eight. The success probability of the CBID detection decreases in a two-cell environment compared to a one-cell environment because of the interference transmitted from the adjacent BS. It can be seen from this figure that the proposed techniques (both Type-1 and Type-2) are robust to the STO and achieve almost the same performance as the previous technique. The reason that the proposed technique, Type-2, does not suffer from a performance degradation is that the parameters that can   Figure 5 shows the success probability of the CBID detection in one-cell and two-cell environments. CBID detection using the preambles simultaneously transmitted from the BS(s), is performed at the MS using Equations (12)- (14). The success of the CBID detection refers to the case where the CID and BID detected at the MS match the ones used in the BS for the beam transmission.
Simulations are performed for three different cases (1ˆ1, 8ˆ1, 8ˆ8), where the first and second terms denote the number of antenna elements in the BS and MS, respectively. When the number of antenna elements is equal to one, an omnidirectional antenna is used. From this figure, it can be seen that an array gain of 9 dB is obtained when the number of antenna elements increases from one to eight. The success probability of the CBID detection decreases in a two-cell environment compared to a one-cell environment because of the interference transmitted from the adjacent BS. It can be seen from this figure that the proposed techniques (both Type-1 and Type-2) are robust to the STO and achieve almost the same performance as the previous technique. The reason that the proposed technique, Type-2, does not suffer from a performance degradation is that the parameters that can avoid the condition in Property 3 and ambiguities in the (GID, SID, BID) detection are selected at the cell planning stage. In the simulation, the CIDs of BS1 and BS2 are selected as 278 and 378. The root index pairs corresponding to the CIDs are {70, 70} and {9, 10}, respectively, where no common root index is used. The performance degradation caused by the intercell interference (high correlation) in the Type-2 preamble can be avoided by selecting appropriate root index pairs at the cell planning stage. Note that the computational complexity required for the CBID detection in the Type-2 preamble is reduced to 33.3%, compared to the Type-1 preamble because N r in the Type-2 preamble is 1/3 that of the Type-1 preamble.
index pairs corresponding to the CIDs are {70, 70} and {9, 10}, respectively, where no common root index is used. The performance degradation caused by the intercell interference (high correlation) in the Type-2 preamble can be avoided by selecting appropriate root index pairs at the cell planning stage. Note that the computational complexity required for the CBID detection in the Type-2 preamble is reduced to 33.3%, compared to the Type-1 preamble because r N in the Type-2 preamble is 1/3 that of the Type-1 preamble.  Figure 6 shows the bit error rate (BER) performance of an 8 × 8 mmWave system in Figure 2. The BER is measured during the data transmission period after completing the cell and beam search period. A quadrature phase shift keying (QPSK) signal is used for the input and the channel impulse response is assumed to be known. In this figure, analytic curves for AWGN and Rician channel (kfactor = 10 dB) with an 8 × 8 beamforming are included for comparison. From this figure, it can be seen that the BER performance of the beamforming system (one-cell) is better than the analytic (Rician) curve because the channel can be approximated as an AWGN when beamforming is   Figure 6 shows the bit error rate (BER) performance of an 8ˆ8 mmWave system in Figure 2. The BER is measured during the data transmission period after completing the cell and beam search period. A quadrature phase shift keying (QPSK) signal is used for the input and the channel impulse response is assumed to be known. In this figure, analytic curves for AWGN and Rician channel (k-factor = 10 dB) with an 8ˆ8 beamforming are included for comparison. From this figure, it can be seen that the BER performance of the beamforming system (one-cell) is better than the analytic (Rician) curve because the channel can be approximated as an AWGN when beamforming is performed in the direction of an LoS path. The performance is slightly worse than the AWGN case because a small power transmitted through the non-LOS path of the Rician channel will not be received in the MS. The previous and proposed (Type-1 and Type-2) techniques exhibit the same performance because the success probability of the CBID detection is one in the range of SNR, >´30 dB, as shown Figure 5. In a one-cell environment, the BER performance of the proposed (or previous) technique with an 8ˆ8 beamforming is similar to the AWGN case. In a two-cell environment, the BER performance is degraded by approximately 4 dB compared to the one-cell case because of the interference from the adjacent BS. performed in the direction of an LoS path. The performance is slightly worse than the AWGN case because a small power transmitted through the non-LOS path of the Rician channel will not be received in the MS. The previous and proposed (Type-1 and Type-2) techniques exhibit the same performance because the success probability of the CBID detection is one in the range of SNR, >−30 dB, as shown Figure 5. In a one-cell environment, the BER performance of the proposed (or previous) technique with an 8 × 8 beamforming is similar to the AWGN case. In a two-cell environment, the BER performance is degraded by approximately 4 dB compared to the one-cell case because of the interference from the adjacent BS. Figure 6. BER performance in one-cell and two-cell environments. Table 3 shows the number of complex multiplications required for CBID detection at the MS for each Rx beam direction. From this table, it can be seen that the numbers of complex multiplications required for the Type-1 and Type-2 preambles are reduced to 1.9% and 0.6%, respectively, compared to the previous technique, when C N is 1020. A significant reduction in the computational complexity can be achieved using the Type-2 preamble. Figure 7 shows the number of complex multiplications required for the CBID detection when C N varies. When C N is 1020, the number of multiplications required for the Type-1 and Type-2 preambles are 1,175,040 and 396,288, respectively. When C N is 24, the number of multiplications required for the previous technique is 1,470,240 that is larger than that for the proposed technique (Type-1 and Type-2), when C N is equal to 1020.
Therefore, the proposed technique (Type-1 and Type-2) can provide an increase in the number of CIDs by a factor of 43, compared to the previous technique, when the same number of complex multiplications is used for the CBID detection.   Table 3 shows the number of complex multiplications required for CBID detection at the MS for each Rx beam direction. From this table, it can be seen that the numbers of complex multiplications required for the Type-1 and Type-2 preambles are reduced to 1.9% and 0.6%, respectively, compared to the previous technique, when N C is 1020. A significant reduction in the computational complexity can be achieved using the Type-2 preamble. Figure 7 shows the number of complex multiplications required for the CBID detection when N C varies. When N C is 1020, the number of multiplications required for the Type-1 and Type-2 preambles are 1,175,040 and 396,288, respectively. When N C is 24, the number of multiplications required for the previous technique is 1,470,240 that is larger than that for the proposed technique (Type-1 and Type-2), when N C is equal to 1020. Therefore, the proposed technique (Type-1 and Type-2) can provide an increase in the number of CIDs by a factor of 43, compared to the previous technique, when the same number of complex multiplications is used for the CBID detection.

Conclusions
In this paper, an efficient preamble design technique is proposed for OFDM-based mmWave cellular systems with beamforming, by combining two ZC sequences in the frequency domain. The parameter selection methods for the proposed preamble are also described for avoiding the performance degradation caused by the side-peaks (high correlation), STO, and ambiguities in the (GID, SID, BID). It was shown by simulation that the proposed techniques (Type-1 and Type-2) can obtain a significant complexity reduction in the CBID detection without a noticeable performance degradation compared to the previous technique. The number of complex multiplications required for the Type-2 preamble is reduced to 0.6% compared to the previous technique. The performance degradation caused by the side-peaks in the Type-2 preamble can be avoided by selecting appropriate root index pairs in the cell planning stage. With a similar computational complexity for the CBID detection, the proposed technique can provide an increase in the number of CIDs by a factor of 43, compared to the previous technique.

Conclusions
In this paper, an efficient preamble design technique is proposed for OFDM-based mmWave cellular systems with beamforming, by combining two ZC sequences in the frequency domain. The parameter selection methods for the proposed preamble are also described for avoiding the performance degradation caused by the side-peaks (high correlation), STO, and ambiguities in the (GID, SID, BID). It was shown by simulation that the proposed techniques (Type-1 and Type-2) can obtain a significant complexity reduction in the CBID detection without a noticeable performance degradation compared to the previous technique. The number of complex multiplications required for the Type-2 preamble is reduced to 0.6% compared to the previous technique. The performance degradation caused by the side-peaks in the Type-2 preamble can be avoided by selecting appropriate root index pairs in the cell planning stage. With a similar computational complexity for the CBID detection, the proposed technique can provide an increase in the number of CIDs by a factor of 43, compared to the previous technique.