EURASIPJournal on Wireless Communications and Networking

RESEARCH

Open Access

Precoded GFDM transceiver with low complexity time domain processing

Maximilian Matthe1*, Luciano Mendes2, Ivan Gaspar1, Nicola Michailow1, Dan Zhang1 and Gerhard Fettweis1

CrossMark

Abstract

Future wireless communication systems are demanding a more flexible physical layer. GFDM is a block filtered multicarrier modulation scheme proposed to add multiple degrees of freedom and to cover other waveforms in a single framework. In this paper, GFDM modulation and demodulation is presented as a frequency-domain circular convolution, allowing for a reduction of the implementation complexity when MF, ZF and MMSE filters are employed as linear demodulators. The frequency-domain circular convolution shows that the DFT used in the GFDM signal generation can be seen as a precoding operation. This new point-of-view opens the possibility to use other unitary transforms, further increasing the GFDM flexibility and covering a wider set of applications. The following three precoding transforms are considered in this paper to illustrate the benefits of precoded GFDM: (i) Walsh Hadamard Transform; (ii) CAZAC transform and; (iii) Discrete Hartley Transform. The PAPR and symbol error rate of these three unitary transform combined with GFDM are analyzed as well.

Keywords: Precoding, GFDM, Flexibility, Low complexity

1 Introduction

Fifth generation (5G) networks will face new challenges that will require a higher level of flexibility from the physical layer (PHY). 3D and 4k video will push the throughput and spectral efficiency. Tactile Internet [1] will demand latencies at least one order of magnitude smaller than what fourth generation (4G) can achieve. Machine type communication (MTC) [2] will require loose synchronization, low power consumption and a massive number of connections. Wireless Regional Area Network (WRAN)

[3] needs to cover wide areas to provide Internet access in low populated regions. The future mobile networks must have a flexible PHY to address several different scenarios.

Generalized Frequency Division Multiplexing (GFDM)

[4] is a recent waveform that can be engineered to address various use cases. GFDM arranges the data symbols in a time-frequency grid, consisting of M subsymbols and K subcarriers, and applies a circular prototype

Correspondence: maximilian.matthe@ifn.et.tu-dresden.de

MATLAB code illustrating the presented techniques can be downloaded from

https://wwwpub.zih.tu-dresden.de/~vf5gdemo/GFDM/

1 Technische Universität Dresden, Vodafone Chair Mobile Communication

Systems, Georg-Schumann-Str. 11, D-01187 Dresden, Germany

Full list of author information is available at the end of the article

filter for each subcarrier. GFDM can be easily configured to cover other waveforms, such as Orthogonal Frequency Division Multiplexing (OFDM) [5] and Single Carrier Frequency Domain Multiple Access (SC-FDMA) [6] as corner cases. The subcarrier filtering can reduce the out-of-band (OOB) emissions, control peak to average power ratio (PAPR) and allows dynamic spectrum allocation. GFDM, with its block-based structure, can reuse several solutions developed for OFDM, for instance, the concept of a cyclic prefix (CP) to avoid inter-frame interference. Hence, frequency-domain equalization can be efficiently employed to combat the effects of multi-path channels prior to the demodulation process. With these features, GFDM can address the requirements of 5G networks.

The main disadvantage of a flexible waveform is the complexity required for its implementation. Reducing the GFDM complexity is essential to bring its flexibility to 5G networks. The understanding of GFDM as a filtered multicarrier scheme leads to a modem implementation based on time-domain circular convolution for each subcarrier, performed in the frequency domain [4].

Springer Open

© 2016 Matthe et al. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a linkto the Creative Commons license, and indicate if changes were made.

This paper follows the well-known polyphase implementation of filter bank modulation [5-7] and applies it to GFDM.

As one main contribution, the investigation in this paper reveals that the Poisson summation formula [8] can be utilized in the demodulation process, considering a per-subsymbol circular convolution and decimation in the frequency domain. This operation can be performed as an element-wise multiplication with subsequent simple M-fold accumulation in the time domain. This simple re-orientation of the data symbol processing allows for a considerable reduction in complexity, which makes GFDM applicable in a wide range of scenarios.

This paper also contributes to reformulate the modulation and demodulation as matrix-vector operations based on M (K)-point discrete Fourier transforms (DFTs). As an outcome, it becomes evident that DFT can actually be considered as a precoding operation. Moreover, they can be replaced by other transformations as well. This observation facilitates arbitrary precoding of the GFDM data, which adds a new level of flexibility to the system, which will be exemplified by the use of three different transforms: (i) Walsh-Hadamard transform (WHT) [9]; (ii) constant amplitude zero auto-correlation (CAZAC) [10] and; (iii) discrete Hartley transform (DHT) [11].

It is worth mentioning that the proposed low complexity signal processing complements the work in [4] because the requirement for block alignment is loosened. It can be usable for supporting pipeline inner receiver implementations, particularly when building synchronization and channel estimation circuits for embedded training sequences [12, 13]. It can also be beneficial to develop future non-linear and recursive detection algorithms.

The remainder of this paper is organized as follows: Section 2 describes the classical approach for GFDM modulation and demodulation and introduce reformulated expressions for a less complex implementation. Section 3 presents the matrix notation to describe the modulation and demodulation processes. Section 4 introduces precoding to the GFDM communication chain. Section 5 shows the complexity analysis of the precoded schemes. Section 6 brings three examples of precoding using different unitary transforms and Section 7 concludes the paper.

2 Classical GFDM description and low-complexity reformulation

Consider a wireless communication system that transmits data in a block-based structure that consists of K subcarriers and M subsymbols. Let N = KM be the total number of data symbols in the block and dk>m, k = 0,..., K — 1,

m = 0,..., M — 1, denote the complex valued data symbol that is transmitted on the //th subcarrier and mth subsymbol. The classical description of GFDM signal generation [14] is given by

[n] = Y^ (g[n] exp | j2n K^n) ) © dk[n]

dk[n] = Y dk,mS[n - mK],

where n = 0,..., N — 1, © describes circular convolution carried out with period N and g[n] denotes the impulse response of the transmit prototype filter.

The subcarrier filtering does not guarantee orthogonality between subcarriers, unless specific prototype filters are used, such as the Dirichlet pulse [15]. Offset Quadrature Amplitude Modulation (OQAM) mapping can also be employed to achieve real-orthogonality [16] and avoid self-interference. In a more general case, the demodulator needs to deal with self-interference. Matched filter (MF) can be combined with successive interference cancellation (SIC) [17] or, alternatively, the self-interference can be mitigated with zero-forcing (ZF) or minimum mean square error (MMSE) filters, for which a simple design based on the Zak transform is available [18].

Using the N-point DFT FN, (1) can be carried out as

c[n] =

Fn \g[n] exp ( j2nk^ Í Fn{dk[n]}

G[f -kM)N ]

FM{dk,()}\f

where <^>N denotes the remainder modulo N and G[f] = Fn{g[n]}. Due to periodicity of with period M and f = 0,...,MK — 1, FM№,(.)} is concatenated K times. This leads to the low-complexity modulator implementation described in [4]. Now, expressing the convolution in (1) explicitly, the transmission equation can be reformulated as [6]

c[n] = Y g [in - mK)N]Y dk,m exp I j2n k=0

K F-1 {d(.),m}\n

The sequence obtained from the Inverse Discrete Fourier Transform (IDFT) is concatenated M times in time-domain due to n = 0,..., MK — 1. Such multiplication in the time-domain expresses the convolution of the subcarrier filter with the data in the frequency domain for each subsymbol.

Assuming a flat and noiseless channel, the convolution of the received signal y[n] = x[n] with a demodulation filter y [n], is given by

dk,m = (V * [<-n)N] © y[n] exp ^ —j2n \n=mK. (5)

With the help of DFT, (5) can be rewritten as

dk,m —Fn

Fn {y*[<-«>n]}

'-.--'

Fn \y[n] exp ( -j2nK-n

Y [<f+kM>N ]

dk,m —

n -1 / kM

— y * [<n - mK >N]y[n] exp ( -j2n kM_n

n—0 ^

— Fn {y* [<n - mK>N]y[n]} |kM,

Fn{u[n] }|kM — Fn \ [<n - mK>N]

I m—0

subcarrier. This constraint is circumvented by performing equalization before demodulation or the presented scheme can be used to initially access the channel state information by demodulation of training data. Also, the proposed scheme facilitates the implementation of odd number of subsymbols [18] avoiding the use of M-point DFT and IDFT required in [4].

Notice that (7) holds when MF, ZF or MMSE filters are used. ZF or MMSE filters can deliver good symbol error rate (SER) performance, while MF needs SIC to remove the SER floor caused by the self-interference. Hence, (7) shall be used to implement ZF and MMSE demodulators, while MF can be considered when the transmit filters are orthogonal or OQAM is employed.

where Fn {y [n]} = r[f ] and Fn {y[n] } = Y[f ]. For example, this leads to the low-complexity demodulator implementation described in [4] for the MF y [n] = g [n], which was later extended to more general filter types in [18]. Now, directly expressing the convolution and sampling in (5) leads to the time-domain multiplication

3 GFDM matrix model

Arranging the data symbols structure leads to the data matrix

in a two-dimensional

D — (d0 d1 ••• d

( d0 \ dï

\dK -J

where the N-point DFT is only evaluated at every Mth sample. According to the Poisson summation formula [8], this can be reduced to a MM = K point DFT

where Acm = {do,m d1m ... dKdenotes the data symbols transmitted in the mth subsymbol.

According to the GFDM description in Section 2, the rows and columns correspond to the time and frequency resources in Fig. 1, respectively. Hence, the /cth row of D represents the data symbols transmitted in the //th subcarrier, d/ = (d/,o d/,1 ... d/,M-1). The time-domain circular convolution of the modulation process can be expressed in the frequency domain as [4]

x — 1£ P^GR^WmdkT

Equation 7 describes a sampled short-time Fourier transform (STFT). This is clear since the GFDM demodulator is actually performing a Gabor transform [18], which can be understood as a sampled STFT. For this reason, (7) does not allow to employ different prototype filters per

For each subcarrier, the data vector is taken to the frequency domain by the M-point DFT matrix WM. The corresponding time domain up sampling operation is realized

in the frequency domain by duplicating the transformed data symbols vectors K times, using the repetition matrix r(k,m) = 1ki ® IM, where IM is an M size identity matrix, 1y is a i x j matrix of ones and <g> is the Kronecker product. Each subcarrier is then filtered with G = diag(WNg), where diag (•) returns a matrix that contains the argument vector on its diagonal and zeros otherwise and g is the vector containing the transmit filter impulse response. An up-conversion of the //th subcarrier to its respective subcarrier frequency is performed by the shift matrix

P(/) = y (p(/)) ® Im, (11)

where y(^) returns the circulant matrix based on the input vector and p(/) is a column vector where the //th element is 1 and all others are zeros. The K subcarriers are summed and transformed back to the time domain with WH to compose the GFDM signal. On the demodulator side, the recovered data symbols for the / th subcarrier are given by

d/ = WH (R(K,M^)T r (p(/))T WNy, (12)

where y is the equalized vector at the input of the demodulator, r = diag(WNy) with y being the demodulation filter impulse response, e.g., MF, ZF, or MMSE filters.

The modulation and demodulation processes can be simplified just by changing the processing order of the data symbols, as derived in (4). In this new matrix model, the GFDM vector is given by

x = J2 P(ra)diag(g)R(MK)wHd^. (13)

Observe that (13) is similar to a polyphase filter structure, but with the difference that cyclic time-shifts are

used here and frequency domain convolution is performed in time domain as element-wise vector multiplication. With this approach, the first step is to obtain a time domain version of dcm by multiplying it with an inverse DFT (IDFT) matrix Wjí. M times upsampling in frequency domain is performed in time by duplicating the transformed data symbols with a repetition matrix R(MK). Each subsymbol is then pulse-shaped with g. P(m) shifts the mth subsymbol to its respective position in time. The GFDM signal is obtained by summing all the pulse-shaped subsymbols, with no need of the N-IDFT domain conversion in (10).

The demodulation operations are also simplified by using the circular-convolution in frequency domain, leading to

dm = Wk (r(M,K))T diag(y) (P(m))T y. (14)

While the new approach describing GFDM modulation and demodulation as convolution in the frequency domain reduces the implementation complexity; in the next section, this principle will be expanded to a more general structure, where domain conversions will be understood as a precoding process.

4 DFT-based precoding for GFDM

Equations (13) and (14) bring a new interpretation of the GFDM chain that further increases the flexibility of this waveform. Let the transmit vector from (13) be redefined as

x = J2 P^diagC^R^V (15)

*m = Td m (16)

is the information coefficient vector to be transmitted and T is a transformation matrix. For classical GFDM, where a time-frequency grid is used to transmit the information, T = WH. However, (16) can be seen as precoding the information vector by a generic matrix T. Different matrices can be used to achieve different requirements. Hence, by transmitting precoded data, GFDM can be even more flexible to address requirements of future wireless networks.

On the receiver side, the recovered coefficients are given by

8m = (r(MK))Tdiag(y) (P(m))Ty. (17)

The estimated data symbols are obtained by reverting the precoding operation on (17), leading to

dm = S8 m. (18)

where S removes precoding. For conventional GFDM,

5 = Wk. A block diagram that describes the proposed GFDM modem is presented in Fig. 2. The demodulation part can be implemented to run continuously, which can be used for tracking specific properties of the estimated information coefficient vector 8m containing training data.

The precoding concept can be enlarged to embrace any transform taking into account the two dimensions of the GFDM data block D, which is defined as a time-frequency grid spanning over several subcarriers and subsymbols, as depicted in Fig. 1. Hence, in the default configuration, the columns of the matrix contain data defined in the frequency domain, while the rows have data in the time domain. The precoded data block A can be generalized to

a = Tc • D • Tr , (19)

K xM K xK K xM MxM

S/P r

m blocks

-ri Cire. Fllt

Cire. Filt

w —>

'Vs, —> rt -^Eb P/S

a) Modulator

Fig. 2 Block diagram of low-complexity modulator and demodulator with precoding

b) Demodulator

where D is now a data block defined in an arbitrary domain, Tc is a K by K precoding matrix applied to the columns and Tr is a M by M precoding matrix applied to the rows of D. For instance, choosing Tc and Tr to be Fourier matrices yields four domains for the GFDM data block,as given in Table 1, where the frequency-time (FT) domain corresponds to standard GFDM.

The choice of the data domain has impact on the robustness of the system regarding time- and frequency- selective fading. The concept is not limited to DFT and, in general, any meaningful transform can be applied to the rows and columns, as exemplified in Section 5.

Taking the concept one step further, (19) can be formulated as

vec(a) = vec(TcDTr). (20)

With the help of the Kronecker product, the right part of (20) can be expressed as

vec(TcDTr) = (tT ® Tc) vec(D). (21)

The resulting matrix

T = (tT ® Tc) (22)

can be seen as an even more general way of precoding that is not restricted to subcarriers or subsymbols but allows arbitrary coupling between any elements of d = vec(D). The precoding is then applied leading to 8 = Td.

5 Out of band, PAPR and complexity analysis

The special precoding examples presented in Section 4 can have specific applications. The FT mode allows for

Table 1 Precoding in time (T) and frequency (F) domains

Domain Operation Number of multiplications

FT A = W^DIM MK log2 (K)

TT A = I*- DIm 0

FF A = WHDWM MK log2 (K) + KM log2 (M)

TF A = I* DWM KM log2 (M)

use of null subsymbols and subcarriers to support non-continuous bands, achieving a low OOB radiation as shown in Fig. 3. The time-time (TT) case leads GFDM to become a single-carrier system occupying the whole system bandwidth. This precoding is especially beneficial for power-limited systems that cannot afford complex signal processing. Another beneficial aspect of the TT precoding is the low PAPR, as can be observed in Fig. 3. In contrast, frequency-frequency precoding (FF) results in a DFT-precoded GFDM system, with a spectrum similar to OFDM. In FF mode, data is defined directly in the frequency domain, according to (3). In this case, one data symbol occupies no more than two frequency bins, assuming a bandlimited transmit filter, such as raised cosine. This property simplifies equalization, as only ICI remains, but on the other hand reduces frequency diversity in channels with strong frequency selectivity. Finally, time-frequency precoding (TF) also corresponds to the single carrier, but with symbols spread in time, creating potential for exploiting time diversity.

The complexity of the precoding operation of the GFDM transceiver based on DFT algorithms are evaluated in Table 1. The number of complex valued multiplications, a costly operation in implementations, is the chosen criterion for the analysis. The comparison is carried out under the assumption that N = KM complex data symbols are transmitted. As the baseline for the comparison in Fig. 3, MN multiplications from (17) are added to the precoding operations and the OFDM algorithm requires naturally the complexity of a DFT operation, with Nlog2(N) multiplications. For purposes of comparison, Table 2 also shows an example assuming a scenario equivalent to Long-Term evolution (LTE) with 20 MHz bandwidth [19]. On the receiver side, the GFDM signal is demodulated with a ZF filter [17], which is a worst-case scenario, since the benefits from a sparse frequency response cannot be harvested.

The GFDM modulator and demodulator proposed in [4] require the complexity of the M-point DFT algorithm with additional K times complex multiplications over repeated chunks of M complex samples, resulting

OFDM e GFDM frequency-time GFDM time-time GFDM frequency-frequency GFDM time-frequency---GFDM time convolution

/ in MHz PAPRo [dB] M

Fig.3 OOB, PAPR and complexity analysis for precoded GFDM. FT/FF: M = 16, Mon = 15, K = 128, Kon = 75,TT/TF: K = 1, M = 1200, OFDM: N = 2048, Non = 1200

from a DFT operation. Even when the frequency response of the demodulation filter is assumed to be sparse, e.g., when roll-off is small and a ZF filter span L can be smaller than the total number of subcarriers, the modem scheme presented in this paper is less complex than the solution proposed in [4]. Another obstacle for the implementation in [4] is that it is based M-point DFT operations and M is optimally chosen to be odd as shown in [18]. Hence, DFT of odd length would need to be implemented which cannot be achieved by radix-2 FFT algorithms.

Recently, an alternative formulation for low-complexity implementation of GFDM modulators and demodulators was shown in [20]. The presented algorithm relies on a block-diagonalization of the transmitter and receiver matrix by using a block-DFT matrix, which can be easily implemented by means of the FFT. The achieved complexity statistics in [20] are comparable to the ones presented in the present paper. However, [20] relies on a specific structure of the transmitter and receiver filters, such that the diagonalization can be employed. Fortunately, at the receiver, this structure is obeyed for MF, ZF and MMSE receiver filters. However, the solution is not general in the sense that any filter at the receiver can be employed with the algorithm in [20]. For example, different filters might be employed when the channel noise is correlated,

polluted by interference or not all subcarriers are allocated. In contrast, the algorithm presented in this paper does not rely on any structure of the receiver filter and hence equally works with any time-domain function y [n] in (5). This is achieved by smartly rearranging the terms in (5) without making any assumptions on y [n].

6 Expanding GFDM features with unitary transform precoding

Several different unitary transforms can be associated with GFDM to achieve specific goals. Three examples based on WHT, CAZAC, and DHT will be explored in the following.

6.1 WHT-GFDM

WHT can be combined with GFDM [9] to increase the system robustness against frequency-selective channels. In this case, the WHT is applied to the data symbols on each subsymbol dcm, leading to

K/2 OK/2

\ Ok/2 -Ok/2 /

Table 2 Number of multiplications of DFT models. Modulation and demodulation complexity are the same with this criterion

Type Matrix model Complexity Resource # op.

OFDM x = N WN d, d = WNy Cofdm = Cdft,n = N log2 (N) K= 2048, M = 1 22,528

GFDM Time conv. K-1 x = NWh L P(k)GR(KM)WMd[ k=0 dk = M Wh (r(K^) t r (P(k))t WNy, vk Cgfdm,zf = cdft,n + L(N + cdft,m) K= 128, M = 15 K 274,202

GFDM Freq. conv. M-1 x = J2 P(m)diagCg)R(MX)WHdm m=0 dm = WK (R(MK))Tdiag(y) (P(m))Ty,Vm Cgfdm = ■■M(N + Cdft,k ) K= 128, M = 15 42,240

with Hi = (1) is the Walsh-Hadamard matrix. Notice that, in this case,

T = 7K wHfflK

and (15) can be used to generate the WHT-GFDM symbol.

On the receiver side, (17) can be used to obtain the recovered coefficients and the reverse precoding matrix

S =-K Ok Wk (26)

is used to obtain the received data vectors, i.e.,

dm = S« m = ^k 2kWK 8 m. (27)

6.2 CAZAC-GFDM

The CAZAC transform is a unitary transform where the coefficients of its matrix are given by [21]

[Ck]y = j= exP jn(i + K)2),

with i = 0,1, ..., K - 1, l = 0,1, ..., K - 1 and CHCk = Ik .

Like the WHT, the CAZAC matrix can be used to pre-code the data symbols in the GFDM system. In this case, the coefficient vector transmitted by the CAZAC-GFDM system is given by [10]

8m = WHCk dm. (29)

Therefore, the precoding matrix for the CAZAC-GFDM is given by

T = WHCk (30)

and the inverse precoding matrix is

The recovered data symbols are given by

S = CHWK.

dm = S8 m = CfWk 8 m.

[0k ]ij = vK ГЧ к) + si4 к)) ■ (33)

Notice that &K &к = Ik .

DHT can be combined with GFDM, leading to the following precoding coefficients

8m = WH©K dm. (34)

The precoding and inverse precoding matrices are given by

t = wHe*

s = e* w*,

respectively. Therefore, the recovered data symbols can be written as

dm = S8 m = e* Wk 8m.

6.4 Precoding performance analysis

The main advantage of the precoding schemes presented in this section is that each subcarrier carries a linear combination of all data symbols from a given subsymbol. This procedure spreads the information over all subcarriers, allowing the receiver to exploit frequency diversity. Since the columns and rows of the precoding matrices are orthogonal to each other, the data symbols can be recovered on the receiver side without self-interference introduced by the precoding operation. Therefore, the precoding GFDM schemes experience a performance gain over conventional GFDM in frequency-selective channels, as can be seen in Fig. 4. The parameters used for simulation are presented in Table 3.

The theoretical symbol error rate for GFDM is given by [17]

' _ 1 \ к-i P =2 I ^ILerfc (jY) -

1 /к - 1

1=0 2 к -1

where k2 = 2^, ¡x is the number ofbits per symbol and 3Rt Es

Yl = ---. (38)

n 2(2x - 1) ^iNo ( )

Es is the average energy per symbol of the quadrature amplitude modulation (QAM) constellation, RT is the throughput efficiency reduction due to CP, N0 is the noise

6.3 DHT-GFDM

The discrete Hartley transform (DHT) is a linear unitary transform with matrix coefficients given by [11]

100 10 10-2 10-3

10-5 10-6 10-7 10

- - - Theo. GFDM

Theo. Precoding GFDM GFDM о WHT-GFDM * CAZAC-GFDM □ DHT-GFDM

10 12 14 16

Es/N„

22 24 26

Fig. 4 Symbol error rate of GFDM with and without precoding over the frequency selective channel presented in Table 3

Table 3 Simulations parameters

Parameter Value

Modulation order 16-QAM

Number of subcarriers 512

Number of subsymbols 5

CP length 16 samples

Transmit filter RC with a = 0.25

Channel impulse response hdB[i] = -2//3 for i = 0,1,.. .,15

spectral density and & is the noise enhancement factor per subcarrier, which is given by

& = MK £

G[ -k] H [k]

with H[k] and G[k] being the channel frequency response and the prototype filter frequency response, respectively.

Equation (37) can also be used to estimate the symbol error rate for precoding GFDM using unitary transform [9] but, in this case,

3\Heq\2 Rt Es

2(2^ - 1) Wo'

\MK ^0 \H[k] \2,

and & must be redefined as

I [A(Im ® T)]-

Another benefit that is obtained by precoding the data symbols before the GFDM modulation is the PAPR reduction, as can be seen in the complementary cumulative density function (CCDF) curves shown in Fig. 5. Since the subcarriers carry a weighted average of all data symbols from a given subsymbol, the probability of achieving high amplitude peaks is reduced, leading to a small PAPR when compared with conventional GFDM. Here, the employed precoding matrix affects the achieved PAPR. The best PAPR reduction is achieved by the CAZAC transform, followed by DHT and WHT, respectively. The explanation for this behavior lies on the aperiodic autocorrelation function [10] of the coefficients 8m. CAZAC transforms lead to coefficients with smaller aperiodic autocorrelation sidelobes. This property leads to a lower PAPR for the CAZAC-GFDM.

7 Conclusions

This paper introduces a GFDM modulator and demodulator based on a low-complexity multiplication in the time domain. The new matrix model allows to consider the coefficients transmitted by the subcarriers as the precoding of the data symbols. This approach broadens the GFDM flexibility, allowing to use different transforms to achieve different goals. Precoding can be used to improve the system SER performance over frequency-selective channels and to reduce the PAPR. The three unitary transforms considered in this paper presented the same SER performance. The CAZAC-GFDM achieved the lowest PAPR among the considered transforms. The computational complexity of the new GFDM modulator and demodulator have been compared with the solution proposed in [4], showing that the modem solution proposed in this paper can reduce the number of complex valued multiplications, even when sparse filters in the frequency domain are employed.

Competing interests

The authors declare that they have no competing interests. Acknowledgements

This work has been performed in the framework of the FP7 ICT-619555 project RESCUE, which is partly funded by the European Union. The work presented in this paper was sponsored by the Federal Ministry of Education and Research within the progarmme "Twenty20 - Partnership for Innovation" under contract 03ZZ0505B - "fast wireless" and within in the framework of the SATURN project with contract no. 100235995 which is funded by the Europaischer Fonds Für regionale Entwicklung (EFRE). This work was partially supported by Finep/Funttel Grant No. 01.14.0231.00, under the Radiocommunication Reference Center (Centro de Referenda em Radicomunicagöes - CRR) project of the National Institute ofTelecommunications - Inatel, Brazil. The computations were performed at the Center for Information Services and High Performance Computing (ZIH) at TU Dresden.

Author details

technische Universität Dresden, Vodafone Chair Mobile Communication Systems, Georg-Schumann-Str. 11, D-01187 Dresden, Germany. 2Instituto

Nacional deTelecomunicagöes, Av. Joäo de Camargo, 510, 37540-000 Santa

Rita do Sapucai, MG, Brazil.

Received: 19 February 2016 Accepted: 8 May 2016 Published online: 25 May 2016

References

1. GP Fettweis, The tactile internet: applications and challenges. IEEE Veh. Technol. Mag. 9(1), 64-70 (2014)

2. P Marsch, B Raaf, A Szufarska, P Mogensen, H Guan, M Farber, S Redana, K Pedersen, T Kolding, Future mobile communication networks: challenges in the design and operation. IEEE Veh. Technol. Mag. 7(1), 16-23 (2012)

3. H Kim, J Kim, S Yang, M Hong, Y Shin, An effective mimo-ofdm system for ieee 802.22 wran channels. IEEE Trans. Circ. Syst. II Express Briefs. 55(8), 821-825 (2008)

4. I Gaspar, et al., in IEEE 77th Vehicular Technology Conference (VTC Spring). Low Complexity GFDM Receiver Based On Sparse Frequency Domain Processing, (2013). http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber= 6692619&tag=1

5. B Farhang-Boroujeny, OFDM versus filter bank multicarrier. IEEE Signal Proc. Mag. 28(3), 92-112 (2011)

6. P Banelli, S Buzzi, G Colavolpe, A Modenini, F Rusek, A Ugolini, Modulation formats and waveforms for 5G networks: who will be the heir of OFDM?: an overview of alternative modulation schemes for improved spectral efficiency. IEEE Signal Proc. Mag. 31 (6), 80-93 (2014)

7. M Renfors, J Yli-Kaakinen, FJ Harris, Analysis and design of efficient and flexible fast-convolution based multirate filter banks. IEEE Trans. Signal Proc. 62(15), 3768-3783 (2014)

8. J Benedetto, G Zimmermann, Sampling multipliers and the poisson summation formula. J. Fourier Anal. Appl. 3(5), 505-523 (1997)

9. N Michailow, L Mendes, M Matthe, I Gaspar, A Festag, G Fettweis, Robust WHT-GFDM for the next generation of wireless networks. IEEE Commun. Lett. 19(1), 106-109 (2015)

10. R Qiu, Z Wu, S Zhu, in 4th International Conference on Wireless Communications, Networking and Mobile Computing. A novel paper reduction method in OFDM systems by using Cazac matrix transform, (2008), pp. 1-4. http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber= 4678106

11. S Bouguezel, MO Ahmad, MNS Swamy, Binary discrete cosine and Hartley transforms. IEEE Trans. Circ. Syst. I Regular Papers. 60(4), 989-1002 (2013)

12. I Gaspar, A Festag, G Fettweis, in 82nd Vehicular Technology Conference (VTCFall). Synchronization using a pseudo-circular preamble for generalized frequency division multiplexing in vehicular communication, (2015), pp. 1-5. http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber= 7391164

13. I Gaspar, A Festag, G Fettweis, in IEEE Global Communications Conference (GLOBECOM2015). An embedded midamble synchronization approach for generalized frequency division multiplexing, (2015), pp. 1-5. http:// ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=7417357

14. N Michailow, I Gaspar, S Krone, M Lentmaier, G Fettweis, in IEEE International Symposium on Wireless Communication Systems (ISWCS). Generalized Frequency Division Multiplexing: analysis of an alternative multi-carrier technique for next generation cellular systems, (2012), pp. 171-175. http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber= 6328352

15. M Matthe, et al., in IEEE International Conference on Communications Workshops (ICC). Influence of pulse shaping on bit error rate performance and out of band radiation of generalized frequency division multiplexing, (2014), pp. 43-48. http://ieeexplore.ieee.org/xpl/articleDetails.jsp? arnumber=6881170&filter%3DAND%28p_IS_Number%3A6881162%29

16. H Lin, P Siohan, Multi-carrier modulation analysis and WCP-COQAM proposal. EURASIP J. Adv. Signal Process. 2014(79), 1371-1430 (2014)

17. N Michailow, et al., Generalized frequency division multiplexing for 5th generation cellular networks. IEEE Trans. Commun. 62(9), 3045-3061 (2014)

18. M Matthe, LL Mendes, G Fettweis, Generalized frequency division multiplexing in a Gabor transform setting. IEEE Commun. Lett. 18(8), 1379-1382 (2014)

19. I Gaspar, et al., in 11th International Symposium on Wireless Communications Systems (ISWCS). LTE-compatible 5G PHY based on Generalized Frequency Division Multiplexing, (2014), pp. 209-213. http://ieeexplore. ieee.org/xpl/login.jsp?tp=&arnumber=6933348&url=http%3A%2F %2Fieeexplore.ieee.org%2Fxpls%2Fabs_all.jsp%3Farnumber%3D6933348

20. A Farhang, N Marchetti, LE Doyle, Low Complexity Modem Design for GFDM. IEEE Trans. Signal Process. PP(99), 1-1 (2015). doi:10.1109/TSP.2015.2502546

21. J Meng, G Kang, in International Conference on Computer Application and System Modeling (ICCASM). A novel ofdm synchronization algorithm based on cazac sequence, vol. 14, (2010), pp. 14-63414637. http://ieeexplore. ieee.org/xpl/login.jsp?tp=&arnumber=5622219&url=http%3A%2F %2Fieeexplore.ieee.org%2Fxpls%2Fabs_all.jsp%3Farnumber%3D5622219

Submit your manuscript to a SpringerOpen journal and benefit from:

► Convenient online submission

► Rigorous peer review

► Immediate publication on acceptance

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article

Submit your next manuscript at ► springeropen.com