Optimum low complexity filter bank for generalized orthogonal frequency division multiplexing

EURASIP Journal on Wireless Communications and Networking, Jan 2018

Generalized frequency division multiplexing (GFDM) is one of the multicarrier modulation candidates proposed for the 5th generation of wireless networks. Among GFDM linear receivers, GFDM MMSE receiver achieves the best error performance for multipath fading channels at the cost of high numerical complexity. Hence, the combination of GFDM match filter (MF) receiver and double-side successive interference cancellation (DSIC) method is used instead. However, there is a significant gap between the error performance of GFDM MMSE and DSIC/MF receivers for the case of employing modern channel coding. Recently, we have proposed a new multicarrier scheme based on GFDM architecture called generalized orthogonal frequency division multiplexing (GOFDM). This study derives an optimized cyclic tree-structured perfect reconstruction-quadrature mirror filter (PR-QMF) bank for GOFDM transceiver and then introduces a novel method for implementation of the optimum filter bank in the frequency domain. Employing such a fast and optimum filter bank provides several advantages for GOFDM transceiver. GOFDM transmitter mitigates out-of-band spectrum leak to the level of that of GFDM. In addition, choosing an appropriate configuration of filter bank yields lower peak to average power ratio in transmit signal of GOFDM compared to that of OFDM. On the other hand, while GOFDM MMSE receiver has lower numerical complexity compared with GFDM DSIC/MF receiver, its coded bit error rate curve is close to that of GFDM MMSE receiver. The aforementioned advantages envision GOFDM as a competitive candidate to be employed in the physical layer of new wireless applications.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1186%2Fs13638-017-1017-x.pdf

Optimum low complexity filter bank for generalized orthogonal frequency division multiplexing

Abbaszadeh et al. EURASIP Journal on Wireless Communications and Networking Optimum low complexity filter bank for generalized orthogonal frequency division multiplexing Mohammad Hadi Abbaszadeh 2 Babak H. Khalaj 0 1 Afrooz Haghbin 2 0 School of Computer Science, Institute for Research in Fundamental Sciences (IPM) , Tehran , Iran 1 Department of Electrical Engineering, Sharif University of Technology , Tehran , Iran 2 Department of Electrical Engineering, Science and Research Branch, Islamic Azad University , Tehran , Iran Generalized frequency division multiplexing (GFDM) is one of the multicarrier modulation candidates proposed for the 5th generation of wireless networks. Among GFDM linear receivers, GFDM MMSE receiver achieves the best error performance for multipath fading channels at the cost of high numerical complexity. Hence, the combination of GFDM match filter (MF) receiver and double-side successive interference cancellation (DSIC) method is used instead. However, there is a significant gap between the error performance of GFDM MMSE and DSIC/MF receivers for the case of employing modern channel coding. Recently, we have proposed a new multicarrier scheme based on GFDM architecture called generalized orthogonal frequency division multiplexing (GOFDM). This study derives an optimized cyclic tree-structured perfect reconstruction-quadrature mirror filter (PR-QMF) bank for GOFDM transceiver and then introduces a novel method for implementation of the optimum filter bank in the frequency domain. Employing such a fast and optimum filter bank provides several advantages for GOFDM transceiver. GOFDM transmitter mitigates out-of-band spectrum leak to the level of that of GFDM. In addition, choosing an appropriate configuration of filter bank yields lower peak to average power ratio in transmit signal of GOFDM compared to that of OFDM. On the other hand, while GOFDM MMSE receiver has lower numerical complexity compared with GFDM DSIC/MF receiver, its coded bit error rate curve is close to that of GFDM MMSE receiver. The aforementioned advantages envision GOFDM as a competitive candidate to be employed in the physical layer of new wireless applications. 1 Introduction Spectrum is a limited media for growing demand of radio communications, and hence, spectrum efficiency is one of the important parameters to be taken into account for any modulation proposed as the air interface of the 5th generation (5G) wireless networks. Low computational complexity of the transceiver algorithm is another important parameter especially for use cases that require low battery consumption. OFDM has been widely used in wireless transmission due to its low computational complexity and good performance for multipath channels [ 1, 2 ]. In order to provide the possibility of using a simple one-tap frequency domain equalizer (FDE) at the receiver side of OFDM, cyclic prefix (CP) is generally added to OFDM symbols. Since CP length should, at least, be equal to the channel length, CP overhead ratio is high for large channel spreads. Moreover, OFDM suffers from high out-of-band spectrum leak [3] which forces many subcarriers to be left vacant in order to fit the OFDM spectrum in an emission spectrum mask. Therefore, OFDM achieves low spectral efficiency due to large guard time and guard band and thus could not be recommended as an appropriate choice for new wireless applications. Among the alternative filtered multicarrier schemes [ 4, 5 ], general frequency division multiplexing (GFDM) is known as a versatile multicarrier modulation being able to address many of the requirements of new wireless applications [6]. GFDM was first introduced in [ 7 ] and then has been widely investigated in the 5GNOW project [ 8, 9 ]. Three common linear receivers, i.e., zero-forcing (ZF), matched filter (MF), and minimum mean square error (MMSE) receivers, are evaluated for the receiver side of GFDM. Among GFDM linear receivers, GFDM MMSE receiver yields the lowest error rate at the cost of the highest complexity which makes it unusable for many applications [10]. Therefore, matched filter (MF) receiver in combination with the double-side successive interference cancellation (DSIC) algorithm is used instead to provide a trade-off between computational complexity and bit error rate (BER) performance [ 11 ]. On the basis of GFDM architecture, we have recently introduced a generic cyclic filtered multicarrier scheme called extended-GFDM. Employing perfect reconstructionquadrature mirror filter (PR-QMF) bank in the generic scheme yields generalized orthogonal frequency division multiplexing (GOFDM). GOFDM benefits from architectural advantages of GFDM including flexibility of the structure, high spectral efficiency, and low latency. We have shown GOFDM MMSE outperforms GFDM ZF receiver in terms of error performance on the assumption of perfect synchronization. In addition, we have proved GOFDM MMSE receiver has a simple structure that consists of onetap MMSE-FDE prior to PR-QMF receiver filter bank [ 12 ]. The first contribution of this paper is to propose an optimized tree-structured PR-QMF bank in order to achieve low out-of-band (OOB) spectrum emission. Simulation results show that employing proposed optimized PR-QMF bank in GOFDM structure results in OOB spectrum leak as low as that of GFDM transmitter. In addition, it is shown GFDM and GOFDM transmitters which employ an appropriate configuration of filter bank can improve peak to average power ratio (PAPR) compared with that of OFDM transmit signal. As the second contribution, a novel method for implementing the cyclic tree-structured filter bank in the frequency domain is introduced. Such implementation results in a lower computational complexity for GOFDM MMSE receiver compared with existing implementation techniques proposed for GFDM DSIC/MF receiver under multipath fading channels. Further, it will be seen GOFDM MMSE receiver yields better error performance compared with OFDM for multipath fading channel where modern error control coding is presented in the setup. In the similar conditions, GFDM encounters complexity and error performance challenges for MMSE and DSIC/MF receivers, respectively. The remaining parts of the paper are organized as follows. In Section 2, the extended-GFDM scheme is described. In Section 3, the first optimum low-complexity PR-QMF bank is derived and then the GOFDM transceiver block diagram based on such filter bank is introduced. Simulation results presented in Section 4 compare GOFDM and GFDM in terms of spectrum localization, computational complexity, error performance, and PAPR considering OFDM as a benchmark. The paper is concluded in Section 5. 2 Extended-GFDM scheme Figure 1 illustrates a general form of baseband multicarrier modulation inspired from GFDM architecture; therefore, it is called extended-GFDM scheme. At the transmitter side, xti(n) for i = 1 : K are K blocks of M complex symbols generated by a QAM mapper. Up-sampling xti(n) sequences by the factor of K generates xteiðnÞ sequences of N = MK symbols whose symbols are given as 8 0 n 1 xteiðnÞ ¼ >>< xti@K A; n ¼ 0; K ; 2K ; …; ðM−1ÞK ð1Þ >>: 0 otherwise: Assume gi(n) and hi(n) for i = 1 : K to be impulse responses of transmitter and receiver filters with length less than or equal to N, we have yðnÞ ¼ XK N i¼1xteiðnÞ ⊛ giðnÞ; where ⊛ N denotes N-point circular convolution. Substituting for xteiðnÞ symbols from (1) yields y(n) symbols as yðnÞ ¼ XK i¼1 XjM¼−01xtiðjÞg~iðn−jK Þ for n ¼ 0 : N −1; where g~iðn−jK Þ is achieved through jK cyclic shift of gi(n) after extending to N symbols through zero padding. Next, inserting CP at the beginning of y(n) sequence generates yt(n) sequence. At the receiver side, removing CP symbols yields yr(n) sequence that may be expressed as the cyclic convolution of y(n) and the impulse response of multipath fading channel. Employing FDE for yr(n) yields yEq(n), where FDE needs an N-point DFT to transfer symbols to frequency domain and an N-point IDFT to return them ð2Þ ð3Þ back to time domain. Subsequently, xreiðnÞ for i = 1 : K are N-symbol sequences obtained as N xreiðnÞ ¼ yEqðnÞ ⊛ hiðnÞ; and then K-factor down-sampling operations yield xri(n) symbols as xriðnÞ ¼ xreiðnK Þ for n ¼ 0 : M−1: Since each gk(n) is a band-passed filter located at a separate part of frequency domain, it can be considered as a subcarrier waveform. Assuming the sampling time to be Ts, (3) shows that it will take KTs seconds to transmit one symbol of each subcarrier, and thus, M symbols of each branch are carried within M time slots taking the total duration of MKTs seconds. Therefore, the scheme is K-subcarrier and M time slot modulation with block based output of N = MK symbols. This result is in contrast with OFDM block (OFDM symbol) in which each subcarrier transmits one symbol during one time slot. The generic scheme gains some architectural advantages that are described as follows: Low latency: Since equalization is processed within one block, extended-GFDM has the advantage of low delay. In addition, cyclic filtering used in extended-GFDM scheme prevents to generate tail symbols, and thus, it has lower latency compared to linear filtered multicarrier modulations. Spectral efficiency: Extended-GFDM improves guard time overhead ratio compared with OFDM through sharing one CP among multiple time slots. Flexible structure: In wireless cellular networks, smallM, large-K configurations can be employed for downlink to provide multiple access and large-M, small-K configurations may be used for uplink to mitigate PAPR. Assume p(n) to be the impulse response of a prototype filter. Substituting giðnÞ ¼ pðnÞ exp j2π ði−K1Þn for i = 1 : K in Fig. 1 creates a modulated filter bank; y(n) is thus obtained as ð4Þ ð5Þ ð6Þ ð7Þ In GFDM transmitter, data of i branch is filtered by use of a prototype filter and shifted in the frequency domain through multiplying by exp j2π ði−1Þn factor [ 13 ]; K y(n) is thus given as XiK¼1 xteiðnÞ⊛N pðnÞ It is straightforward to show (6) and (7) are equal, and therefore, employing giðnÞ ¼ pðnÞ exp j2π ði−1Þn K extended-GFDM scheme leads to GFDM structure. in the 3 Filter bank design for GOFDM transceiver GOFDM emerges through employing synthesis and analysis PR-QMF banks at the transmitter and receiver sides of extended-GFDM scheme, respectively. Since GOFDM follows extended-GFDM structure, it takes the common architectural advantages of extended-GFDM described in previous section. In this section, we will address the design of PR-QMF bank to achieve low OOB spectrum leak and low order of computational complexity when employed in the GOFDM structure. 3.1 Background As seen in Fig. 2, PR-QMF bank consists of synthesis filter bank as multiplexing part and analysis filter bank as demultiplexing part. Connecting the output of the synthesis part directly to the input of the analysis part will completely reconstruct the original inputs of the filter bank. Furthermore, in PR-QMF bank, synthesis filters are matched to corresponding analysis filters [ 14 ]. In addition, perfect reconstruction property is preserved even when the convolutions are performed circularly, and thus, PR-QMF banks can also be classified as cyclic orthogonal filter banks [ 15 ]. Two-band PR-QMF banks are the simplest form of PRQMF banks which are easily designed by using a prototype filter. Assume h(n) to be L-tap prototype filter-derived such that L is an even number and its Z transform satisfies the following condition: a b HðzÞH z−1 þ Hð−zÞH −z−1 ¼ Constant: ¼ RðzÞ þ Rð−zÞ Perfect reconstruction property yields [ 16 ] h1ðnÞ ¼ hðnÞ; h2ðnÞ ¼ ð−1Þðnþ1ÞhðL−1−nÞ; g1ðnÞ ¼ hðL−1−nÞ; and g2ðnÞ ¼ ð−1ÞnhðnÞ; ð8Þ ð9Þ ð10Þ ð11Þ ð12Þ where h1(n) and h2(n) are low-pass and high-pass filters at the synthesis part and g1(n) and g2(n) are low-pass and high-pass filters at the analysis part of 2-band PRQMF bank, respectively. Based on (8), considering R(z) = H(z)H(z−1), the L-tap prototype filter of 2-band PRQMF bank can be obtained as the causal part of 2L-tap half-band filter R(z). The half-band filter R(z) is a linear phase FIR filter, and thus, we may use the ParksMcClellan algorithm to derive R(z) so as to provide the desired features for prototype filter H(z) [ 17 ]. Using a uniform tree structure of log2K levels, 2-band filter banks can be used for building K-band filter banks for the case where K is a power of two. By using the same prototype filter for generating 2-band filter banks in all levels of the tree, we may obtain the conventional tree-structured filter bank. However, 2-band filter banks used in different levels of a tree, even at each node of the same tree level, can be different for the most general case. Figure 3 shows 4-band tree-structured filter bank in which the filters with superscripts 1 and 2 belong to the first and second levels of the tree, respectively. 3.2 Optimized tree-structured PR-QMF bank In this part, we introduce the tree-structured PR-QMF filter bank with localized frequency response sub-band filters. Our approach is the design of non-conventional tree structure wherein different levels of tree have different 2-band filter banks but 2-band filter banks of the same level are similar. Let M be the size of input sequences at the synthesis part of K-band tree-structured filter bank and hi(n) be the impulse response of the prototype filter employed to derive 2-band filter banks used in ith level of the tree. The OOB spectrum leak depends on the behaviour of the filters within the passband and stopband. The optimization is thus aimed in designing hi(n) for i = 1 : log2K such that the following optimization problem is satisfied min Ei Subject to Ei ¼ R f ∉passband Hiðf Þ df variables Hiðf Þ RiðzÞ þ Rið−zÞ ¼ Constant RiðzÞ ¼ HiðzÞHiðz−1Þ ð13Þ where Hi(f ) and Hi(z) are frequency response and Z transform of hi(n), respectively. Substituting hi(n) for h(n) in (9–12) yields gi1ðnÞ; gi2ðnÞ; hi1ðnÞ and hi2ðnÞ and, subsequently, the optimized K-band tree-structured filter will be established when the 2-band filter banks used in the ith level of tree are obtained by use of the derived impulse responses gi1ðnÞ; gi2ðnÞ; hi1ðnÞ and hi2ðnÞ for i = 1 : log2K. In our approach, first, the Parks-McClellan algorithm is used to design Ri(n) where the designing parameters are length, cut-off frequency and maximum ripple of the filter and then, hi(n) is derived as the causal part of Ri(n). We determine the designing parameters of Ri(n) so that Ei is minimized for hi(n). Noting that going from each level to the successive higher one, the sequence size is duplicated, and hence, the size of input sequences at the ith level of tree is 2i − 1M which is changed to 2iM after up-sampling by the order of two. In order to preserve the size of filtered sequences through cyclic filtering, the filters length has to be less than or equal to the size of filtered sequences, and thus, we have ð14Þ Li ≤ 2iM; where Li is the length of hi(n). Since more spread in the time domain permits designing of more localized sub-band filters in the frequency domain, the length of hi(n) is chosen to have the maximum possible value, i.e., Li = 2iM. As a consequence, the order of respective half-band filter Ri(n) will be 2 × 2iM − 1. In addition, cut-off frequency is set at the middle of frequency space (i.e., 0.5 for normalized frequency space). Ri(n) is then derived by use of Parks-McClellan algorithm and on the conditions that the order and cut-off frequency are fixed and the maximum ripple is numerically determined such that the achieved hi(n) minimizes Ei. Figure 4 compares frequency response of 32-tap optimum prototype filter (case 2) versus two non-optimum 32-tap prototype filters of cases 1 and 3 with larger and smaller maximum-ripples, respectively. The filters are derived as the causal parts of 64-tap half-band filters with the order and cut-off frequency equal to 61 and 0.5, respectively. For case 2, the maximum-ripple of respective half-band filter is numerically derived as 0.0121 while the maximum-ripples of cases 1 and 3 are chosen as 0.1 and 0.001, respectively. The values of Ei for cases 1, 2, and 3 are 0.0944, 0.064, and 0.0712, respectively. The results show optimum filter with lowest Ei value may have larger maximum-ripple compared to non-optimum prototype filters. 3.3 Low-complexity tree-structured PR-QMF bank In this section, we propose a practical implementation algorithm in order to reduce the computational complexity of the optimum cyclic tree-structured filter bank. In order to reach this goal, the cyclic tree-structured filter bank is implemented in the frequency domain using Discrete Fourier Transform (DFT) and Inverse DFT (IDFT), where M-point DFT and IDFT are defined as: 1 XðkÞ ¼ pffiMffiffiffiffi XM−1 n¼0 xðnÞW kMn for k ¼ 0 : M−1 1 xðnÞ ¼ pffiMffiffiffiffi XM−1 k¼0 XðkÞW −Mkn for n ¼ 0 : M−1; ð15Þ ð16Þ where W M ¼ e−jð2MπÞ . Time domain operations for a 2band cyclic filter bank are cyclic convolution, upsampling, and down-sampling by the order of two. Throughout the paper, we will denote the frequency domain operations equivalent to the time domain factor 2 up-sampling and down-sampling as 2-Copy and 2-Fold, respectively, where the definitions of 2-Copy and 2-Fold are derived in the following lemmas. Lemma 1 If the M-symbol sequence x(n) is up-sampled by the factor of 2 to generate xe(n), then Xe(k)=2-Copy{X(k)}, where X(k) is a sequence of M-point DFT coefficients of x(n), Xe(k) is a sequence of 2M-point DFT coefficients of xe(n) and 1 pffi2ffi XðkÞ f or k ¼ 0 : M−1 XeðkÞ ¼ 2−CopyfXðkÞg ¼ f pffi2ffi Xðk−MÞ f or k ¼ M : 2M−1 1 . Lemma 2 If 2M-symbol sequence x(n) is down-sampled by the factor of 2 to generate xd(n), then Xd(k)= 2-Fold {X(k)}, where X(k) is a sequence of 2M-point DFT coefficients of x(n) and Xd(k) is a sequence of M-point DFT coefficients of xd(n) and 2−FoldnXðkÞg ¼ p1ffi2ffi ½XðkÞ þ Xðk þ MÞ for k ¼ 0; 1; …; M. Proof In the earlier study, we have expressed up- and down-sampling by the factor K in the frequency domain by use of DFT coefficients [ 11 ]. Substituting K with 2 in the derived results yields the results of Lemmas 1 and 2. Further, we know that cyclic convolution of the Msymbol sequence x(n) with the M-coefficient filter g(n) is equivalent to the multiplication of DFT coefficients of x(n) and frequency domain coefficients of g(n) given by XM−1 n¼0 gðnÞW kMn for k ¼ 0 : M−1: ð17Þ Substituting up-sampling, down-sampling, and cyclic convolution operations of K-band tree-structured filter bank with 2-Copy, 2-Fold, and multiplication operations yields a new structure that we call it as frequency domain tree structure (FDTS) filter bank. Consequently, in order to completely implement the filter bank in the frequency domain, at the first stage of synthesis part, a set of K blocks of M-point DFT are required to transfer input blocks to the frequency domain. The K-band FDTS synthesis filter bank is the subsequent stage, and then, N-point IDFT is required to return these symbols from the frequency domain back to the time domain. At the analysis part, the first stage is an N-point DFT. The K-band FDTS analysis filter bank and K blocks of M-point DFT are the subsequent stages. Figure 5 shows how to implement a 4-band treestructured filter bank in frequency domain where xiðnÞ and x^ iðnÞ for i = 1 : 4 are M-symbol sequences and y n ð Þ and y^ ðnÞ are 4M-symbol sequences, respectively. In general case, based on the results of the previous section derived for optimum filter bank, when the inputs at the synthesis part are M-symbol sequences, gi1ðnÞ; gi2ðnÞ; hi1 ðnÞ and hi2ðnÞ are 2iM-coefficient filters. Substituting gi 1 ðnÞ; gi2ðnÞ; hi1ðnÞ and hi2ðnÞ for g(n) in (17) yields the frequency domain coefficients of the filters to be employed at the ith level of FDTS filter bank, i.e., Gi1ðkÞ; Gi2ðkÞ; H i1 ðkÞ and Hi2ðkÞ. Finally, it should be noted that since cyclic filtering creates a clockwise circular shift, outputs of the analysis part of PR-QMF bank are not exactly equal to the corresponding inputs of synthesis part. In order to compensate such shifts and obtain full reconstruction cyclic filter bank, it is enough to change frequency domain coefficients of the synthesis filters for i = 1 : log2K as i G^ 1ðkÞ ¼ Gi1ðkÞW 22iiMMðk−1Þ for k ¼ 0 : 2iM−1 i G^ 2ðkÞ ¼ Gi2ðkÞW 22iiMMðk−1Þ for k ¼ 0 : 2iM−1: ð18Þ ð19Þ 3.4 Complexity analysis In this part, the complexities of frequency and time domain implementations of tree-structured filter bank are compared. The complexity is calculated in terms of the number of needed basic operations of additions and multiplications to generate one output block. The filter bank is the optimal K-band tree-structured filter bank where the inputs at the synthesis part are sequences of M symbols. In general, the coefficients of the filters may be complex numbers, and so, all multiplications and additions are generally complex instructions. We assume that the filter coefficients and the DFT coefficients of the filters are known at the receiver side and they are registered in the buffer in the initialization step. 3.4.1 Frequency domain implementation Figure 6a, b shows a typical 2-band frequency domain synthesis and analysis filter banks used at the ith level of the synthesis and analysis parts of optimum FDTS filter bank, respectively. Let V1(k), V2(k) be 2i − 1M-symbol input sequences of 2-band synthesis filter bank. Based on Lemma 1, the next 2-Copy operations generate 2iM-symbol sequences whose symbols are given as 1VeðkÞ ¼ V e1 k þ 2i−1M ¼ V 1ðkÞ for k ¼ 0 : 2i−1M−1 2VeðkÞ ¼ V e2 k þ 2i−1M ¼ V 2ðkÞ for k ¼ 0 : 2i−1M−1; and, subsequently, the output symbols of synthesis part are obtained as ð20Þ ð21Þ a b For 2-band analysis filter bank, letting W^ ðkÞ be 2iMsymbol input sequence yields V^ 1ðkÞ ¼ W^ ðkÞ V^ 2ðkÞ ¼ W^ ðkÞ H i1ðkÞ for k ¼ 0 : 2iM−1 H i2ðkÞfor k ¼ 0 : 2iM−1; subsequently, based on Lemma 2, output symbols of 2band analysis filter bank are obtained as 1 V^ 1ðkÞ þ V^ 1 k þ 2i−1M V^ 1dðkÞ ¼ 2 1 V^ 2ðkÞ þ V^ 2 k þ 2i−1M V^ 2dðkÞ ¼ 2 for k ¼ 0 : 2i−1M−1 for k ¼ 0 : 2i−1M−1 Based on (22), the 2-band frequency domain synthesis filter banks used at the ith level of FDTS synthesis filter bank require 2 × 2iM multiplications and 2iM additions. The K-band FDTS synthesis filter bank consists of log2K levels, and there are 2Ki 2-band filter banks at the ith level. Therefore, the number of multiplications and additions performed in the synthesis part of FDTS is Pi¼lo1g2K 2Ki 2 2iM¼2KM log2K Pi¼lo1g2K 2Ki 2iM¼KM log2K In order to implement K-band synthesis filter bank in the frequency domain, the total basic operations needed consists of basic operations for the calculation of K blocks of M-point DFT, K-band FDTS, and N-point IDFT. K is always a power of two and N=MK , thus assuming that M to be a power of two yields N to be also a power of two. Therefore based on the radix-2 FFT algorithm, the overall number of basic operations is obtained as ð22Þ ð23Þ ð24Þ ð25Þ ð26Þ mults adds: ð27Þ On the other hand, referring to Eqs. (23–26), 2-band analysis filter banks used at ith level of analysis part need 2 × 2iM multiplications for two multiplication operators and 2 × 2i − 1M additions for two 2-Fold operators. Therefore, the 2band analysis and synthesis filter banks used at the ith level of FDTS filter bank need equal number of basic mathematical calculations, and consequently, the total number of basic mathematical operations of synthesis and analysis parts of frequency domain implementation are similar. 3.4.2 Time domain implementation The inputs of the ith level of synthesis part are blocks of 2i − 1M symbols, and thus, after up-sampling, the size of input blocks is changed to 2iM where half of them are zeroes. Two-band filter banks used at the ith level of analysis part require two 2iM-point cyclic convolutions and one 2iM-point addition at the end. Therefore, number of basic operations performed for the 2-band synthesis filter banks used at ith level of synthesis part is 8 < 22i M2 : 2iM 2iM−1 mults adds: The tree is built of log2K levels where the ith level at the synthesis part consists of 2Ki 2-band synthesis filter bank. Therefore, the total number of basic operations of the synthesis part yields as (P log2K Ki 22i M2 ¼ NMX log2K 2i Pi¼lo1g2K 2Ki 2iMð2iM−1Þ¼NPi¼lio¼1g21K ð2iM−1Þ i¼1 2 In a similar approach, it can be shown that the numbers of basic operations at the analysis part is equal to mults adds: ð28Þ ð29Þ mults adds: ð30Þ 3.5 GOFDM system model based on optimum lowcomplex filter bank Deploying the optimum low-complexity treestructured PR-QMF bank in the generic scheme of Fig. 1 yields the block diagram of Fig. 7 wherein lower case letters denote time domain symbols and higher case letters denote frequency domain symbols. As stated before, the advantage of implementing treestructured PR-QMF filter bank in frequency domain instead of time domain is reducing the computational complexity of the structure. This advantage is more highlighted at the receiver side wherein the existing FDE implies that the received symbols are transferred to the frequency domain even in case of time domain implementation. 4 Simulation results 4.1 Simulation setup In what follows, we compare GOFDM with GFDM in terms of spectrum localization, error performance, numerical complexity, and peak to average power ratio. Optimum tree-structured PR-QMF bank is employed for GOFDM. In the GFDM structure, the prototype filter is root-raised cosine (RRC) filter with roll-off factor of 0.3. Setting up the configuration parameters of filter banks needs some consideration. Provided that M and K being even numbers GFDM transmit matrix becomes singular [ 13 ] but PR-QMF bank does not face such a problem and we are permitted to use M and K both as a power of 2. In the simulations, we evaluate two configurations where the number of subcarriers in GFDM and GOFDM schemes for configuration cases 1 and 2 are 256 and 16, respectively. Therefore, we encounter the GFDM limitation of using odd number for M values. Hence, M is chosen as 9 and 129 for GFDM scheme and the nearest power 2 numbers, i.e., 8 and 128 are chosen for GOFDM scheme in the configuration cases 1 and 2, respectively. In two GOFDM configurations, M × K = 2048 and the number of CP symbols is equal to 64; thus, both configurations have similar spectral efficiency. When OOB emission is analyzed, configuration with large number of subcarriers, i.e., case 1 is taken into account while PAPR is analyzed for configuration with small number of subcarriers, i.e., case 2. Both configurations are evaluated in terms of error performance and numerical complexity. The channel profile and other simulation parameters are presented in Table 2. 4.2 OOB leakage In order to reduce OOB emission, the normal approach is to use empty subcarriers at the corner sides of assigned bandwidth. For block-based modulations, abrupt changes of signal between successive blocks result in high OOB emission. Windowing is thus exploited at the transmitter side by multiplying output blocks with a window function to smooth abrupt changes and combat the OOB emission. For the case of windowed waveforms in addition to CP, cyclic suffix (CS) is also inserted. Through the simulations, G prefix denotes windowed waveform and GW prefix denotes windowing and guard subcarriers are simultaneously applied for a waveform. According to LTE standard for the case of 256 subcarriers, 150 subcarriers are occupied as allocated bandwidth. We consider six empty subcarriers of each side as guard band and, consequently, power spectral densities (PSDs) leak to other empty subcarriers may be considered as OOB emission. For the case of GOFDM, despite what might be expected, frequency responses of sub-band filters of the K-band tree-structured filter OFDM 64 16 16 samples in both sides Not used 256 1 256 1 16-QAM bank are disorderly arranged in the spectrum for the case where K is higher than 2. The 4-band frequency split of the spectrum is illustrated in Fig. 8 for ideal band-pass filters. As explained, 4-band tree-structured filter bank has two stages, stage 1 consists of two 2band filter banks and stage 2 consists of one 2-band filter bank. It is seen in Fig. 8 that filter bank of stage 2 reverses the arrangement of two higher sub-bands and the frequency response of sub-band 3 comes before sub-band 4. Therefore, the order of 4-band tree-structured filter bank in the spectrum may be shown as {1, 2, 4, 3}. Similarly, 8band tree-structured filter bank may be considered as a stage of two 4-band tree-structured filter banks prior to one 2-band filter bank. Combining two 4-band treestructured filter banks using next 2-band filter bank reverses the order of higher 4-band filter bank, and thus, the sub-bands are ordered in the spectrum as {1,2,4,3,7,8,6,5}. Following such approach, one can find the order of subband filters of K-band tree-structured filter bank in the spectrum is matched to a sequence of gray numbers. In our case study, the 256-number gray sequence starting from 1 shows the order of subcarriers in the spectrum. Figure 9 compares PSDs of OFDM, GFDM, and GOFDM. In order to evaluate OOB emission, OOB radiation parameter is defined as the ratio of the energy that is emitted into the OOB range and the amount of energy within the allocated bandwidth, i.e., jBWj :R f ∈OOBPðf Þdf OOB radiation ¼ jOOBj R f ∈BWPðf Þdf ð31Þ where BW and OOB are the set of frequencies that are considered as in-band and out of band, respectively. The OOB radiations of G-OFDM, G-GFDM, and G-GOFDM are calculated as − 25.7, − 34.2, and − 33.9 dB, respectively. In case of applying both windowing and guard subcarriers, the OOB radiations of GW-OFDM, GW-GFDM, and GW-GOFDM are reduced to − 35.6, − 50.5, and −49.8 dB, respectively. As expected, the OOB emissions of cyclic multicarrier modulations are significantly less than that of OFDM. Furthermore, the optimized PR-QMF bank causes GOFDM OOB radiations to be only 0.3 and 0.7 dB more than that of GFDM for the cases of inserting guard subcarriers and employing window and inserting guard subcarriers together, respectively. We now compare the bandwidth efficiency for three modulations. Noting that the overhead symbols are CP and CS guard symbols in the time domain and empty guard subcarriers in the frequency domain, we define the spectral efficiency as the ratio of useful symbols to the total transmitted symbols during each output block, and then, the spectral efficiency of OFDM, GOFDM, and GFD25M6 8are 256þ26546þ16 ¼ 0:76; 256 2956þ694þ32 ¼ 0:96 and 256 8þ64þ32 ¼ 0:95 , respectively. The results show that GOFDM transmitter yields almost the same level of OOB leakage and spectral efficiency as those of GFDM transmitter. 4.3 Error performance In the case of GFDM MMSE receiver, an inverse matrix has to be calculated according to changes in noise power. Therefore, GFDM MMSE receiver suffers from high complexity, and DSIC/MF receiver is suggested to use instead. In the earlier study, we have proved, on the assumption of perfect synchronization, GOFDM MMSE receiver which simply consists of one-tap MMSE-FDE and analysis PR-QMF bank is the best GOFDM linear receiver in terms of error performance [ 12 ]. On the other hand, GFDM DSIC/MF receiver is the combination of ZF-FDE, matched filters and DSIC algorithm. One-tap coefficients of MMSE-FDE and ZF-FDE for k = 1:N are defined as Cð1kÞ and jCðCkÞjð2kþÞσ2w , respectively, where C(k) and σ2w are N-point DFT coefficients of channel impulse response and the power of white Gaussian noise, respectively. In this study, OFDM, GFDM DSIC/MF with four iterations of DISC algorithm, and GFDM MMSE and GOFDM MMSE receivers are compared in terms of error performance for transmitting windowed signals through multipath fading channel when complete synchronization is assumed. At the first step, uncoded BER is analyzed. As seen in Fig. 10a for the case of configuration 1, while GOFDM MMSE slightly outperforms OFDM, the BER curves of GFDM MMSE and DSIC/MF receivers are tightly close to that OFDM. For the configuration case 2 shown in Fig. 10b, GOFDM and GFDM MMSE receivers outperform OFDM in terms of error performance for signal to noise ratio (SNR) larger than 12 and 15 dB, respectively, and the situations are reversed for SNR smaller than these thresholds. However, there is a considerable gap between GFDM DSIC/MF and OFDM error performance curves which means the negative effect of self-interference at the receiver side of GFDM MF receiver is more serious for configuration case 2 compared with configuration case 1. At the second step, the parallel concatenated convolutional code (PC-CC) with code rate R = 1/3 is employed for transmitter side, and a soft demapper in combination with a turbo decoder with 10 iterations is employed at the receiver side. It is known that the performance of employed advanced receiver depends on the mean square error in the received symbols. Therefore, as seen in Fig. 10c, d, while GOFDM MMSE outperforms OFDM, there are significant gaps between coded BER curves of DSIC/MF GFDM and OFDM in both configurations. On the other hand, the coded error performance of GOFDM MMSE receiver compared to GFDM MMSE receiver depends on the configuration of filter bank and SNR value. It can be verified that BER curve of GOFDM MMSE is slightly lower/higher than that of GFDM MMSE for SNR smaller/larger than 9.7 and 7.5 dB for the configuration cases 1 and 2, respectively. 4.4 Computational complexity Among the methods introduced to implement GFDM receivers, authors in [ 18 ] proposed the lowest complexity modem for GFDM zero-forcing and matched filter receivers. However, the proposed method does not support DSIC/MF receiver for the case of multipath fading channel. Therefore, sparse frequency domain processing method introduced in [ 19 ] still leads to lowest complexity for our case study. Based on this method, the number of multiplications in GFDM DSIC/MF receiver with the configuration of K subcarriers and M time slots ignoring the equalizer complexity is reduced to N ½ð2J þ 1Þ log2M þ ðJ þ IÞ þ log2N ; ð32Þ where I is the number of adjacent sub-bands wherein the frequency response of prototype filter has significant values and J is the number of iterations for DSIC algorithm. I depends on the frequency response of prototype filter and usually is small; however, J depends on the order of QAM constellation and for the case of large QAM constellation, large J is required to achieve good error performance [ 19 ]. Ignoring the equalizer part, the complexity of GOFDM receiver with the configuration of K subcarriers and M time slots is equivalent to the complexity of K-band analysis filter bank with an input block of N = MK symbols. When M and K are a power of two, numerical complexity for such configuration is derived in (28). We compare numerical complexity of GFDM DSIC/ MF and GOFDM MMSE receivers in terms of the number of multiplications. In the case of GOFDM, K and M are power of 2 numbers where K is ranged from 16 to 2048 and M is ranged from 4 to 128. For the case of GFDM, K has the same range but M values increased by 1, i.e., M = 5, 9, 17, 33, 65, 129. The values of M in GFDM and GOFDM configurations are not equal; thus, we consider GOFDM configuration parameters as the reference and calculate its complexity relative to GFDM complexity with equal K but M values increased by 1. As seen in Fig. 11a, b, GOFDM MMSE receiver achieves less numerical complexity compared to GFDM DSIC/ MF receiver especially for the case of large M and large J values. For example, in configuration case 1, GOFDM to GFDM complexity ratio for J = 4 and J = 8 are 0.4474 and 0.2724, respectively, and in configuration case 2, this 4.5 PAPR The PAPR of yt(n) sequence is defined as Larger PAPR values demand a higher linear range amplifier which is a challenge for designing low-cost terminals. Complementary cumulative density function (CCDF) of PAPR, i.e., the probability that the PAPR exceeds a certain value, is the typical measure to compare different systems in terms of PAPR. PAPR value of multicarrier schemes depends on the number of subcarriers, and the schemes with a few number of subcarriers perform better than the schemes with large number of subcarriers [ 20 ] in terms of PAPR. Authors in [ 21 ] showed considering GFDM multicarrier scheme for the uplink of a cellular network wherein only one subcarrier assigned to each user yields smaller PAPR compared to OFDM scheme. Therefore, we consider GFDM and GOFDM with the configuration case 2 for an uplink scenario where one subcarrier assigned to each of 16th users. For the case of OFDM, 16 contiguous subcarriers are allocated to each user. As the benchmark, single-carrier FDM system is considered. Then, a pulse-shaping raised cosine filter with roll-off = 0.5 is applied in all schemes. The average values of CCDF of PAPR of transmit signals of the waveforms are shown in Fig. 12. It is seen, OFDM is significantly outperformed by GFDM and GOFDM in terms of PAPR. It can also be verified that GOFDM is superior to single-carrier FDM but it is outperformed by GFDM in terms of PAPR. According to results of previous parts, GFDM transmitter achieves such superiority at the costs of more computational complexity or lower error performance of the receiver side. 5 Conclusions GOFDM encompasses the architectural advantages of GFDM including high spectral efficiency, low latency, and flexible structure. The novelty of current work relies on exploiting a fast and optimized PR-QMF bank to be employed in the GOFDM structure. As a consequence, GOFDM transmitter yields the same level of OOB leakage as that of GFDM transmitter and GOFDM receiver yields lower computational complexity compared with GFDM DSIC/MF receiver. Further, it was shown that GOFDM MMSE receiver is superior to OFDM in terms of error performance when modern channel coding is used in the setup. In the similar conditions, GFDM DSIC/MF receiver could not perform adequately, and thus, it was required to use high complexity MMSE receiver for GFDM to achieve acceptable error performance. Furthermore, it was shown employing the configurations of small numbers of subcarriers, large numbers of time slots in GOFDM architecture preserves good spectral efficiency and, besides, yields PAPR even lower than single-carrier FDM scheme. For such configurations, although GFDM transmitter is superior to GOFDM transmitter in terms of PAPR, GFDM MMSE and DSIC/MF receivers suffer from significant cost of numerical complexity and error performance, respectively. In spite of aforementioned advantages, GOFDM has to be precisely evaluated for the case of asynchronous transmitters which is a challenging issue for new 5G applications. Abbreviations 5G: 5th generation; CCDF: Complementary cumulative density function; CP: Cyclic prefix; CS: Cyclic suffix; DSIC: Double-side successive interference cancellation; FDE: Frequency domain equalizer; FDTS: Frequency domain tree structure; GFDM: Generalized frequency division multiplexing; GOFDM: Generalized orthogonal frequency division multiplexing; MF: Matched filter; MMSE: Minimum mean square error; OOB: Out-of-band; PAPR: Peak to average power ratio; PR-QMF: Perfect reconstructionquadrature mirror filter; PSD: Power spectral density; RRC: Root-raised cosine; ZF: Zero-forcing Acknowledgements This work was supported in part by a grant from IPM and in part by a grant from the Iran National Science Foundation under grant 95824827. Funding There is no funding for the current research. Availability of data and materials The custom codes are available via email to corresponding author. Authors’ contributions MHA and BHK have been responsible for providing the idea of the paper. MHA has also been responsible for implementation of ideas and performing simulations. MHA, BHK, and AH have been involved in writing the manuscript and proofreading it. Authors’ information Mohammad Hadi Abbaszadeh received his B.Sc. and M.Sc. degrees from Shiraz University and Iran University of Science and Technology, respectively. He is currently pursuing the Ph.D. degree in Tehran Science and Research Branch, Islamic Azad University, Tehran, Iran. His research interests include MIMO wireless communications and multicarrier modulations. Babak Hossein Khalaj received his M.Sc. and Ph.D. degrees from Stanford University. At Stanford, he has been among the pioneering team working on adoption of multi-antenna arrays in mobile networks. He has also been a senior member of Advanced Communications Research Institute (ACRI) at Sharif University, Tehran, Iran, and the recipient of Alexander von Humboldt Fellowship in 2007–2008. Afrooz Haghbin received her M.Sc. and PhD degrees from Tehran University and Tarbiat Modares University, respectively. She is currently with the electrical and computer department of Science and Research Branch in Islamic Azad University, Tehran, Iran, as an assistant professor. Her research interests include MIMO wireless communications, channel coding, precoding, multicarrier modulation, and estimation theory. Competing interests The authors declare that they have no competing interests. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 1. SA Fechtel , A Blaickner , Efficient FFT and equalizer implementation for OFDM receivers . IEEE Trans. Consum. Electron . 45 ( 4 ), 1104 - 1107 ( 1999 ) 2. M Mirahmadi , A Al-Dweik , A Shami , BER reduction of OFDM based broadband communication systems over multipath channels with impulsive noise . IEEE Trans. Commun . 61 ( 11 ), 4602 - 4615 ( 2013 ) 3. JVD Beek , F Berggren , Out-of-band power suppression in OFDM . IEEE Commun. Lett . 12 ( 9 ), 609 - 611 ( 2008 ) 4. B Farhang-Boroujeny , OFDM versus filter bank multicarrier . IEEE Signal Process. Mag . 28 ( 3 ), 92 - 112 ( 2011 ) 5. V Vakilian , T Wild , F Schaich , St Brink, JF Frigon, in 2013 IEEE Globecom Workshops (GC Wkshps) . Universal-filtered multi-carrier technique for wireless systems beyond LTE ( 2013 ) 6. M Van Eeckhaute , A Bourdoux , P De Doncker , F Horlin , Performance of emerging multi-carrier waveforms for 5g asynchronous communications . EURASIP J. Wirel. Commun. Netw . 2017 ( 1 ), 29 ( 2017 ) 7. G Fettweis , M Krondorf, S Bittner, in VTC spring 2009 - IEEE 69th Vehicular Technology Conference . GFDM-generalized frequency division multiplexing ( 2009 ) 8. N Michailow , M Matthe , IS Gaspar , AN Caldevilla, LL Mendes , A Festag , G Fettweis , Generalized frequency division multiplexing for 5th generation cellular networks . IEEE Trans. Commun . 62 ( 9 ), 3045 - 3061 ( 2014 ) 9. G Wunder , P Jung , M Kasparick, T Wild , F Schaich , Y Chen , ST Brink , I Gaspar , N Michailow , A Festag , L Mendes , N Cassiau , D Ktenas , M Dryjanski , S Pietrzyk , B Eged , P Vago , F Wiedmann, 5gnow : non-orthogonal, asynchronous waveforms for future mobile applications . IEEE Commun. Mag . 52 ( 2 ), 97 - 105 ( 2014 ) 10. N Michailow , S Krone , M Lentmaier , G Fettweis , in 2012 IEEE Vehicular Technology Conference (VTC Fall) . Bit error rate performance of generalized frequency division multiplexing ( 2012 ) 11. R Datta , N Michailow , M Lentmaier , G Fettweis , in 2012 IEEE Vehicular Technology Conference (VTC Fall) . GFDM interference cancellation for flexible cognitive radio PHY design ( 2012 ) 12. MH Abbaszadeh, BH Khalaj, A Haghbin, Error performance analysis for generalized orthogonal frequency division multiplexing . Computers & Electrical Engineering 61 , 139 - 150 ( 2017 ) 13. M Matthe , LL Mendes , G Fettweis , Generalized frequency division multiplexing in a Gabor transform setting . IEEE Commun. Lett . 18 ( 8 ), 1379 - 1382 ( 2014 ) 14. PP Vaidyanathan , Multirate Systems and Filter Banks (Prentice- Hall , Englewood Cliffs, 1993 ) 15. PP Vaidyanathan , A Kirac, in 1997 IEEE International Conference on Acoustics, Speech, and Signal Process.. Theory of cyclic filter banks ( 1997 ) 16. AN Akansu , R.A.H. , Multiresolution Signal Decomposition: Transforms, Subbands, and Wavelets (Academic pres, Inc, 1992 ) 17. AV Oppenheim, RWSchafer, JR Buck, Discrete-Time Signal Processing Upper. Saddle River , N.J. , Prentice Hall , (18th ed.). ( 1999 ) 18. A Farhang , N Marchetti , LE Doyle, Low-complexity modem design for GFDM . IEEE Trans. Signal Process . 64 ( 6 ), 1507 - 1518 ( 2016 ) 19. I Gaspar , N Michailow , A Navarro , E Ohlmer , S Krone , G Fettweis, in 2013 IEEE 77th Vehicular Technology Conference (VTC Spring) . Low complexity GFDM receiver based on sparse frequency domain processing ( 2013 ) 20. H Seung Hee , L Jae Hong , An overview of peak-to-average power ratio reduction techniques for multicarrier transmission . IEEE Wirel. Commun . 12 ( 2 ), 56 - 65 ( 2005 ) 21. N Michailow , G Fettweis, in 2013 International Symposium on Intelligent Signal Processing and Communication Systems . Low peak-to-average power ratio for next generation cellular systems with generalized frequency division multiplexing ( 2013 )


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1186%2Fs13638-017-1017-x.pdf

Mohammad Hadi Abbaszadeh, Babak H. Khalaj, Afrooz Haghbin. Optimum low complexity filter bank for generalized orthogonal frequency division multiplexing, EURASIP Journal on Wireless Communications and Networking, 2018, 12, DOI: 10.1186/s13638-017-1017-x