We analyse some of the candidates for modulations for 5G: FBMC-OQAM, GFDM-OQAM, and WCP-COQAM. Unlike most of the related bibliographies, which are oriented to mobile communications, our research is focused on 5G in cognitive radio based industrial wireless communications. According to the ultrareliability and low-latency requirements of industrial communications, we simulate the aforementioned modulations in low-band transmissions (carrier frequencies below 6 GHz and a bandwidth narrower than 100 MHz) through large indoor spaces and severe multipath channels that emulate industrial halls. Moreover, we give detailed information about WCP-COQAM and how the windowing affects the protection against multipath effect and reduces spectral efficiency compared to GFDM-OQAM. We also compare the aforementioned filtered multicarrier techniques and OFDM in terms of robustness against multipath channels, power spectral density, and spectral efficiency. Based on these results, we aim at providing an approximate idea about the suitability of 5G MCM candidates for industrial wireless communications based on CR.

1. Introduction

Multicarrier modulation (MCM) schemes are widely adopted in wireless communication systems, orthogonal frequency division multiplexing (OFDM) being the most extended one among them. Despite the features and benefits that OFDM offers in wireless communications, it also presents significant drawbacks such as high out-of-band (OOB) radiation. Considering this fact and the upcoming applications and needs in wireless communications, a lot of research on new waveforms for the physical (PHY) communication layer of the fifth generation (5G) has been carried out in order to overcome OFDM’s limitations [14]. As a result, several pulse-shaping techniques exist, which allow subcarrier filtering in MCM systems. These techniques provide good frequency localization and reduce the OOB radiation. Among the main PHY candidates for 5G, the most extended MCM schemes in the bibliography are filter-bank multicarrier offset quadrature amplitude modulation (FBMC-OQAM) [5], generalized frequency division multiplexing (GFDM) [6], and universal filtered multicarrier (UFMC) [7]. FBMC-OQAM and GFDM consist in filtering each subcarrier in order to get good frequency localization in each subcarrier; UFMC, on the contrary, consists in filtering groups of subcarriers in order to reduce only the OOB radiation.

Although most of the research carried out on PHY proposals for 5G is oriented to mobile communications, 5G is actually meant to cover a wide range of application fields [8]. In order to cope with the diversity of services and the flexibility that 5G will demand, in [9, 10] innovative frameworks allowing multiservice subband filtered multicarrier systems are introduced. These frameworks present a structure in which the available spectrum is divided into subbands dedicated to different services. The division and interference avoidance between different services is provided by means of subband filtering.

Among all the application fields that 5G will cover, our research focuses on wireless industrial communications [11]. Specifically, we investigate robust PHY layers designed for industrial wireless communications based on cognitive radio (CR) [12, 13]. Since the existing industrial wireless communication systems based on the IEEE 802.15.4 standard lack robustness against these harsh channel conditions [14, 15], new schemes for the PHY layer in industrial wireless communications must be proposed. This kind of industrial scenarios is characterized by the presence of high impulsive noise and interference from other wireless systems and fading channels due to severe multipath effect and shadowing caused by the numerous, large, and metallic obstacles within the propagation channel. In this work, we simulate transmissions through large indoor spaces and severe multipath channels in order to model some of these characteristics.

Besides the aforementioned harsh environments, requirements for industrial communications suppose an extra challenge for wireless systems, especially in some factory automation (FA) applications. Since 5G is also meant to cover these industrial communications’ ultrareliability and low-latency requirements (up to packet error rates, near 100 s over-the-air transmission time, and 1 ms end-to-end latency) [1618], in this work, we simulate and evaluate some of the MCM 5G candidates in low-band transmissions (carrier frequencies below 6 GHz and a bandwidth narrower than 100 MHz).

In this article, among the aforementioned industrial environment characteristics, we focus on severe multipath effect. We analyse FBMC-OQAM, GFDM-OQAM, and also windowed cyclic prefix-circular OQAM (WCP-COQAM) modulation systems. Although WCP-COQAM [1921] is not as extended as the other two schemes in the bibliography, it shares several similarities with GFDM, which make it another candidate for 5G [22] and an interesting MCM scheme to be analysed. We evaluate their performance in terms of bit error rate (BER), power spectral density (PSD), and spectral efficiency (SE) under highly dispersive channels. We compare these results with traditional OFDM and we use the SE to represent the cost of using cyclic prefix (CP) extensions, subcarrier filtering, and windowing schemes in each MCM system.

Based on those analyses and results, we aim at providing a notion about the suitability of 5G MCM candidates for industrial wireless communications based on CR.

From the analysis and the assessment of the aforementioned MCM systems, we provide the following contributions:(1)This article complements the work in [1921] by bringing additional details about how windowing affects the protection against multipath channels and the SE in WCP-COQAM. We state the conditions to provide full orthogonality in multipath channels and we show how windowing reduces SE in WCP-COQAM compared to GFDM-OQAM. We explain all these details in Sections 2.3, 2.4, and 3 and provide some results and discussion in Section 4.3.(2)While most research about 5G is focused on mobile communications, here we simulate and assess some of the main 5G modulation candidates under different conditions. We simulate low-band transmissions through large indoor spaces and severe multipath channels in order to model some of the wireless industrial environments’ characteristics. Under these conditions and by means of BER, PSD, and SE analysis, we evaluate the suitability of the aforementioned MCM systems for wireless industrial communications based on CR.

The paper is organised as follows: in Section 2, we explain the theoretical background of the MCM systems we are going to analyse; in Section 3, we give more detailed explanations about WCP-COQAM and make some clarifications with respect to the bibliography; in Section 4, we present the results obtained from our simulations, based on which we compare and discuss the performances of the considered MCM systems; in Section 5, we state the conclusions we obtained from our research; finally, in Section 6, we explain which will be the next steps in our investigation.

2. Multicarrier Modulations

The MCM systems we analyse are based on OQAM instead of classical IQ modulations (e.g., N-PSK and N-QAM). The reason why we focus on OQAM-based schemes is that MCM systems cannot simultaneously keep good time-frequency localization (TFL), Nyquist symbol rate, and orthogonality between transmitted symbols if conventional IQ modulations are used, as it is stated by the Balian-Low Theorem (BLT) [23]. According to BLT, intersymbol interference (ISI) and intercarrier interference (ICI) introduced by the prototype filters of these MCM systems make it impossible to achieve perfect reconstruction of complex valued symbols transmitted at Nyquist rate. Consequently, in case matched filter receivers are used, FBMC and GFDM perform worse than OFDM in terms of BER. One solution to overcome the constraints defined by the BLT is to send alternately real and imaginary valued symbols at twice the Nyquist rate instead of complex symbols. This technique is known as OQAM [5]. This strategy allows the filtered data symbols to remain orthogonal and keep both good TFL and Nyquist symbol rate, all at the same time.

Hence, we first introduce the OQAM technique and next we explain the rest of the MCM schemes.

2.1. OQAM

OQAM consists, basically, in splitting one complex symbol into two semisymbols, one real and one imaginary. After this operation, the whole OQAM symbol’s duration remains equal to the duration of the original complex symbol, while the duration of each semisymbol is half the whole symbol, as shown in Figure 1. Thus, OQAM symbol structure’s sampling rate is twice the original complex symbol structure’s. It is worth mentioning that this oversampling does not carry any change in the total bandwidth of the transmitted signal when OQAM is used in MCM systems. As we explain in Section 2.2, when the multicarrier OQAM signal is serialized, its symbol rate is the same as it would be in a conventional IQ modulation.

The OQAM scheme discussed up to now can reduce ISI caused by pulse-shaping processes. Since imaginary and real semisymbols are transmitted alternately, interference coming from adjacent semisymbols can be ignored in the receiver (i.e., interference caused by a real semisymbol to its adjacent imaginary semisymbol can be ignored and vice versa).

Similarly, in frequency domain, pulse-shaping makes adjacent subcarriers overlap each other, thus causing severe ICI. If those subcarriers are organised so that real and imaginary ones are alternated one next to the other, ICI from adjacent subcarriers will be avoided because interference coming from them can be ignored in the receiver. Figure 2 shows the block diagram of an OQAM MCM system. Figure 3, on the other hand, shows the demodulation process to recover the original complex data symbols.

More details about OQAM can be found in [5].


FBMC consists in filtering each subchannel in order to get well-localized subcarriers. Figure 4 shows the basic representation of a FBMC-OQAM transmultiplexer. The main blocks are the OQAM preprocessing or modulator block, the synthesis filter bank, the analysis filter bank, and OQAM postprocessing or demodulator.

The key components of the FBMC-OQAM scheme are the synthesis and analysis filter banks, which are formed by one filter per subchannel. The filter of the th subchannel is defined aswhere is the prototype filter of length and . The exponential factor corresponds to the th subcarrier, where is the total number of subcarriers.

Note that the FBMC-OQAM scheme works as a multirate system. Before filtering the OQAM data symbols, they are upsampled by using a zero padding process to get the oversampled signal :The reason why the upsampling factor is and not is that OQAM preprocessing already introduces an upsampling factor of 2. Thus, the FBMC-OQAM signal can be expressed as

At the receiver side, the FBMC-OQAM signal is demodulated so that, in first place, each subchannel is filtered in order to get only its corresponding signal:where .

Then, is sampled in the time domain, so that the estimation of the original OQAM signal is obtained:

Transmission channel, noise, equalization, and other issues have been omitted in order to focus just on FBMC-OQAM waveform modulation and demodulation. More details about FBMC-OQAM can be found in [5].


The structure of a GFDM-OQAM transmultiplexer is similar to that of the FBMC-OQAM transmultiplexer shown in Figure 4. GFDM-OQAM is also based on OQAM modulation and subchannel filtering. The difference is that in GFDM-OQAM filtering is performed by a circular convolution instead of the linear convolution used in FBMC-OQAM. This way, the modulation system now adopts a block-based signal structure, so that a CP can be added to provide orthogonality against multipath channels without affecting signal’s TFL.

We consider a GFDM-OQAM system with subchannels and complex valued data symbols in each block (i.e., OQAM symbols in each block). Thus, if we consider as one GFDM-OQAM block, it can be expressed aswhere and filters are obtained by the periodic repetition of filters defined in (1). So .

At the receiver side, the GFDM-OQAM signal is demodulated so that, in first place, each subchannel is filtered in order to get only its corresponding signal. Once again, the difference between FBMC-OQAM and GFDM-OQAM receivers consists in the circular convolution:where .

Then, is sampled in the time domain, so that the estimation of the original OQAM signal is obtained, just as it is shown in (5).

The insertion of a CP, transmission channel, noise, equalization, and other issues have been omitted in order to focus just on GFDM-OQAM waveform modulation and demodulation. More details about GFDM-OQAM can be found in [19, 20]. In those works, the idea of GFDM-OQAM is introduced, although the authors refer to this scheme as FBMC-COQAM or just as COQAM.


WCP-COQAM is introduced in [21]. This technique can be considered as a GFDM-OQAM system with a windowing process, which improves the PSD with respect to GFDM-OQAM.

In WCP-COQAM, both the CP and the windowing are related. The CP is divided into two parts: the guard interval (GI) and the window interval (WI) (referred to as RI in [21]). So the total length of the CP will be equal to the sum of the lengths of the GI and WI parts:where GI is the part of the CP aimed at protecting against multipath channels and WI is the part that will be used for windowing.

Let us consider a GFDM-OQAM signal as a queue of several GFDM-OQAM blocks, each with its CP extension of length . Equally, we will consider to be the total number of subcarriers and to be the number of symbols in each block. Then the th block of a WCP-COQAM signal is defined aswhere ; , and the window function is defined byIn our work, we use a Hamming window for the window coefficients. Figure 5 shows the block diagram that represents this whole process.

It must be noted from (9) and (10) that the samples of the CP corresponding to the WI part and also the last samples of the block are multiplied by the window coefficients and they are overlapped with the WI regions of the previous and next blocks, respectively. By overlapping the WI regions, these extra samples are prevented from reducing SE, as shown in Figure 6.

The receiver structure is depicted in Figure 7. It is worth noticing that Figure 7 only shows the recovery process of one GFDM-OQAM block from a received WCP-COQAM signal. Once the GFDM-OQAM signal is recovered, there are several algorithms for its demodulation. In Section 2.3 of this paper, the basic demodulation scheme for GFDM-OQAM is explained and in [21] the authors introduce a computationally efficient demodulation scheme.

In [21], for this scheme, perfect synchronization and an ideal channel are assumed. So, in the first place, the first samples of one WCP-COQAM block are taken. Next, the first (which correspond to the CP extension) samples are removed, while the last samples are left to process them within the next block. Finally, in order to get all the samples back to their original positions within the original GFDM-OQAM block, a cyclic shift is carried out, which can be expressed aswhere ; is the estimation of the original signal, and is the signal after the removal of the windowing effect and the CP.

3. Orthogonality Condition and Spectral Efficiency Analysis of WCP-COQAM

Several analyses on 5G modulation candidates are available in the bibliography [24, 25], some of them even including WCP-COQAM modulation scheme [26, 27]. However, to the best of our knowledge, none of them address the implications of windowing in the WCP-COQAM waveform as we do here.

In this section, we contribute some additional details regarding how windowing affects the orthogonality against multipath channels and the SE in WCP-COQAM. Specifically, we state the condition that the GI region of the CP extension must fulfil in order to provide full orthogonality in the presence of multipath channels. Besides, we also compare WCP-COQAM and GFDM-OQAM (due to the similarities between both MCM systems) in terms of SE and we show that, in this aspect, the former is outperformed by the latter if equal multipath protection is assumed.

3.1. Orthogonality Condition

In WCP-COQAM, multipath channel protection is directly related to the CP extension and also to the windowing process. In order to explain this fact, we show in detail the steps to form a WCP-COQAM block starting from a basic GFDM-OQAM block (as seen in Section 2.3).

Once a GFDM-OQAM block is formed, the CP extension is added to it. Figure 8(a) shows the details of this first stage in the construction of WCP-COQAM signal, where is the total number of subcarriers; is the number of symbols per block; WI is the region of length devoted to windowing and overlapping; GI is the region of length that protects against multipath channel effect; and the CP extension has a total length of . , , and are those last samples of the block which are copied in order to form the CP extension. As is represented in Figure 8(a), before the windowing is applied, this signal can be considered as a GFDM-OQAM block with CP extension.

Figure 8(b) represents the WCP-COQAM block right after windowing is applied. As shown in Figure 6 and in (9) and (10), the window function is applied on the first and the last samples, so that these samples are somehow distorted. Because of that, CP, WI, and GI are no longer equal to CP′, WI′, and GI′, respectively.

Since the last samples are affected by the window, we define a new region for them which we call . Thus, now the original GI and regions are divided into GI1 and GI2 and GI1′ and , respectively, so that , , and . As for the lengths of these regions, and . Note that GI1 and are now the only really redundant parts. Because of that, only GI1 and not the whole GI (as it is explained in [21]) acts as a real guard interval against multipath channels when frequency domain equalization (FDE) is performed. Hence, must be fulfilled in order to maintain full orthogonality through a multipath channel, where is the transmission channel length. Thus, considering that and so , we can define the full orthogonality condition for WCP-COQAM asIt is worth mentioning that the FDE should be performed for samples between GI2 and , inclusive, because that is the section corresponding to the cyclic convolution between GI1 and with the transmission channel.

3.2. Spectral Efficiency Analysis

Figure 9 shows how overlapping between adjacent blocks is performed in WCP-COQAM. Every windowed region is overlapped with the next or the previous block, so that the WI region of each block is overlapped with the region of the previous block. Thus, the samples within WI regions do not suppose an overload and they do not affect the SE. Based on this fact, in [21], the authors claim that the use of the window function has no effect on SE. However, this statement is not precise. Although the samples in WI do not directly affect SE, because they overlap with the samples of of the previous block, windowing implies adding the GI2 region in order to fulfil the orthogonality condition defined in (12). Note that the use of the window function and overlapping of WI and regions distort those samples in region. This and the fact that FDE is performed between GI2 and , inclusive, imply that the GI2 region is necessary in order to recover those samples originally placed in the region. At this point, we conclude that both GI1 and GI2 are essential parts within the GI region of WCP-COQAM.

In order to make a fair analysis of SE in WCP-COQAM, we compare it to GFDM-OQAM, because both modulation schemes only differ in windowing and overlapping processes. In order to maintain a homogeneous notation between GFDM-OQAM and WCP-COQAM, we assume the following considerations: (for the sake of readability, from this point on, parameters corresponding to OFDM will be expressed with the subindex , parameters corresponding to GFDM-OQAM will be expressed with the subindex , and parameters corresponding to WCP-COQAM will be expressed with the subindex ).

The criterion we applied for this comparison is to provide equal protection against multipath channels for both modulations, so that . Thus

From this analysis, it is clear that if full orthogonality is to be maintained, the windowing process brings an extra GI2 region to the CP which has effects on the SE of a WCP-COQAM signal. In particular, SE is reduced by a factor of compared to a windowless system with equal multipath channel protection.

4. Simulation Results and Discussion

In this section, we show and discuss the results obtained from the simulations of the systems previously described. Based on these results, we assess the performances of the simulated systems in terms of BER, PSD, and SE. In our analysis, we prioritize robustness against highly dispersive channels, so that we use the BER of transmissions through Rayleigh fading multipath channels in order to assess the robustness of the MCM systems. In addition to BER, we also take into consideration the PSD of the transmitted signals, for it is an important factor when it comes to CR systems. Regarding SE, we use it to represent the cost of the CP extensions, windowing process, and convolution tails caused by the prototype filters in FBMC-OQAM, so that we can obtain a fair comparison between different MCM systems.

Since we are evaluating MCM systems, we consider it appropriate to compare the analysed FMC modulation techniques with an OFDM reference system in order to assess their performance. Table 1 shows the configuration parameters we used in most of our simulations. We choose parameters suitable for low-band transmissions in ultrareliable and low-latency indoor scenarios. For that, we consider aspects like subcarrier bandwidth and transmission channel’s coherence bandwidth, data block duration, and window length for PSD improvement.

4.1. Robustness against Multipath Channels

In this subsection, we analyse the robustness of each MCM system against multipath channels. In order to make a fair comparison, we simulate every system with equal protection against multipath channels, matching effective GI lengths (, , and ). These parameters are selected as shown in Table 1.

The equalization technique we use in our simulations is one-tap zero-forcing (ZF) FDE for all the MCM systems. While OFDM, GFDM-OQAM, and WCP-COQAM employ CP extensions to prevent ISI caused by the transmission channel, FBMC-OQAM could be considered a special case in this aspect. Since it is not a blockwise modulation and no CP extension is used, one-tap FDE is not a straightforward equalization technique for FBMC-OQAM. However, due to the subcarrier bandwidth and symbol duration near to transmission channel’s coherence bandwidth and time, respectively, we consider the one-tap FDE as a suitable equalization method for FBMC-OQAM in this case. An in-depth analysis about the doubly dispersive channels’ impact on FBMC systems is given in [28].

We simulated transmissions through a Rayleigh fading channel model with exponential power delay profile (PDP) whose root mean square (rms) delay spread and channel length are  ns and  taps. Besides, we simulate uncoded and coded communications, with turbo codes, soft decoding, and a code rate of . In Figure 10, we compare the BER versus of the MCM systems over these channel models, assuming perfect synchronization and full channel state information (CSI). We got these BER curves by simulating transmissions through Monte Carlo iterations for each MCM system and each channel model. We transmit 12288 random data bits and we generate new random channels at each Monte Carlo iteration.

Here we show that, for uncoded transmissions, in terms of error rate, FBMC-OQAM performs slightly worse than the rest of MCM systems under the simulated conditions. While at low values every MCM system provides similar BER, an error floor appears at high values because of the lack of GI in FBMC-OQAM. On the other hand, although OFDM, GFDM-OQAM, and WCP-COQAM provide similar performances in terms of BER providing equal multipath effect protection, OFDM slightly outperforms the other two MCM schemes due to its perfect orthogonality.

Finally, it is worth noting that, for coded transmissions, all the simulated MCM systems present similar BER curves, so that they provide similar robustness against multipath channels.

4.2. Power Spectral Density Analysis

PSD is a relevant characteristic in order to design a CR communication system suitable for scenarios with high spectrum occupation. In this subsection, we compare the spectrums of the signals of the analysed MCM systems in order to assess their potential suitability for CR based systems.

Figure 11 shows the spectrums of OFDM, GFDM-OQAM, WCP-COQAM, and FBMC-OQAM signals one on top of the other. For this test, we only activate the 32 central subcarriers in order to make the power leakage into adjacent subcarriers visible.

As shown in Figure 11, the FB based MCM systems carry significant reduction of OOB radiation compared to OFDM. However, even between these MCM schemes, there is still a remarkable difference. In this sense, FBMC-OQAM outperforms both GFDM-OQAM and WCP-COQAM. On the other hand, thanks to the windowing, WCP-COQAM shows considerably lower OOB radiation than GFDM-OQAM. It is worth mentioning that the longer the WI region is, the narrower WCP-COQAM signal’s spectrum gets. Since extra samples within the GI region are needed for windowing, as described in Section 3, PSD improvement of WCP-COQAM over GFDM-OQAM comes at the cost of lower SE. More detailed analysis on SE is given in Section 4.3.

4.3. Spectral Efficiency Analysis

In this subsection, we analyse SE in order to evaluate the cost of obtaining higher robustness against dispersive channels and better PSD by means of CP extensions and windowing. For this analysis, we consider the parameters defined in Table 1, so that effective multipath channel protection is equal for OFDM, GFDM-OQAM, and WCP-COQAM. Equations (15) show the SE values of the three blockwise MCM systems:Here we show a particular case with some specific parameters. On one hand, GFDM-OQAM presents higher SE than OFDM because its GI precedes symbols, while in OFDM the GI precedes only one symbol. On the other hand, because of the windowing operation, the CP extension of WCP-COQAM, specifically the GI part, has extra samples with respect to GFDM-OQAM and that is the reason why WCP-COQAM also presents lower SE than GFDM-OQAM. However, under these conditions, with a high number of subcarriers, we can conclude that the differences in the SE between OFDM, GFDM-OQAM, and WCP-COQAM are not really significant.

Regarding FBMC-OQAM, once again, we consider it different from the rest of MCM systems because it is not a blockwise modulation and it uses no GI protection for multipath effect. So its SE depends on the length of the transmitted frames with respect to the convolution tails. For long data transmissions, its SE will tend to 1, while for short data transmissions it should tend to 0.5, assuming a minimum amount of data samples equivalent to a GFDM-OQAM or WCP-COQAM block. This calculation is simple if we introduce 4 complex data symbols (8 OQAM symbols) in (2).and we introduce these upsampled data in (3). The resultant signal will contain 4095 samples, of which half will correspond to actual data and the rest to convolution tails from subcarrier filtering.

Considering the low-latency scenarios our research focuses on, SE might suppose a significant drawback for FBMC-OQAM compared to GFDM-OQAM and WCP-COQAM during short data transmissions.

5. Conclusions

In this article, we make a comparison between three candidates for PHY layer modulation in 5G and a reference OFDM system in order to assess their suitability for industrial wireless communications based on CR. For that, we simulate low-band communications and highly dispersive indoor channel scenarios. Under these conditions, we show the performances and features of these MCM systems from different points of view.

Regarding robustness against multipath channels, we conclude that all the analysed MCM systems provide similar performances under the simulated conditions. Although for uncoded transmissions FBMC-OQAM presents an error floor at high values, this error floor is corrected when coding is used.

In terms of PSD and OOB radiation, these filtered MCM schemes provide more restrained spectrum than OFDM. In particular, FBMC-OQAM outperforms the rest of the modulation schemes in this aspect, so it could be a suitable candidate for the CR scenarios considered in our research. On the other hand, although not as good as FBMC-OQAM, WCP-COQAM still has considerably better PSD than GFDM-OQAM, which only provides less than 10 dB of difference in OOB radiation with respect to OFDM. Therefore, we conclude that circular filtering only is not sufficient and windowing must be applied in order to ensure more efficient use of the spectrum. Moreover, we consider that, due to its also restrained spectrum, WCP-COQAM might be another suitable modulation scheme for CR applications.

As for the SE analysis, as discussed in Section 4.3, the difference between OFDM, GFDM-OQAM, and WCP-COQAM is not really significant for the parameters we used in our simulations. On the other hand, FBMC-OQAM might suffer from severe degradation of SE during short data transmissions in the low-latency scenarios we consider in our research.

6. Future Work

In this work, we considered perfect synchronization and full channel state information at the receiver, but it is also important to consider more realistic simulation conditions in the future. So the next step in our investigation will be to analyse and propose time-frequency synchronization and channel estimation methods for GFDM-OQAM and WCP-COQAM.

The other key point in our research will be to simulate channel models that fairly emulate industrial wireless channel conditions. For that, we will simulate severe multipath conditions, time-varying channels, impulsive noise, and interference from other wireless systems.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


This work was partly supported by TIPOTRANS (ETORTEK) and CIIRCOS (PC2013-68) projects of the Basque Government (Spain) and by COWITRACC (TEC2014-59490-C2-1-P/2-P) project of the Spanish Ministry of Economy and Competitiveness.