SciELO - Scientific Electronic Library Online

 
vol.43 número3Carrier frequency offset compensation for ofdma systems using circular banded matricesSet-membership estimation theory for coupled mimo Wiener-like models índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

  • No hay articulos citadosCitado por SciELO

Links relacionados

  • No hay articulos similaresSimilares en SciELO

Compartir


Latin American applied research

versión impresa ISSN 0327-0793versión On-line ISSN 1851-8796

Lat. Am. appl. res. vol.43 no.3 Bahía Blanca jul. 2013

 

An efficient spectral shaping method for ofdm systems using correlated interpolation of symbols

E.M. Lizarraga*†§, A.A. Dowhuszko, V.H. Sauchelli§ and S.H. Gallina

* Laboratorio de Comunicaciones Digitales-CONICET. lizarem@efn.uncor.edu, emlizarraga@conicet.gov.ar
Universidad Nacional de Catamarca, Argentina. sgallina@tecno.unca.edu.ar
Department of Commun. and Networking (Comnet), Aalto University, Finland. alexis.dowhuszko@aalto.fi
§ Universidad Nacional de Córdoba, Argentina. vsauch@com.uncor.edu

Abstract — The generation of suitable orthogonal frequency division multiplexing (OFDM) signals on the grounds of fully digital signal processing is considered. The main objective is to obtain a discrete-time signal with adequate allocation of power emissions in both, in-band portion (i.e., the allocated band for communication) and out-of-band portion (i.e., the band allocated to adjacent channels) of the spectrum. The proposed method prevents the transmitter from using traditional filtering techniques, to keep under control power emissions in the system. In addition, the adaptability feature of our proposal makes its implementation attractive within cognitive radio (CR) and software defined radio (SDR) OFDM-based systems. Our spectral shaping approach is based on an optimum interpolation, obtained from the combination of an inverse fast Fourier transform (IFFT) and a spectral precoding operation, both of them transparent from the perspective of a conventional (legacy) OFDM receiver.

Keywords — Correlated Interpolation; IFFT; Double-Length; OFDM; N-Continuous; Symbol Merging; Spectral Precoding; Spectral Shaping.

I. INTRODUCTION

Orthogonal frequency division multiplexing (OFDM) is an attractive modulation technique that is not only widely used in current commercial systems (e.g., xDSL, DVB, WiFi, and WiMAX), but it also represents the selected candidate to be implemented in the air interface of future high-speed mobile communications standards (i.e., LTE and LTE-Advanced). The available spectrum in an OFDM transmission is occupied by orthogonal subcarriers, which are utilized to convey parallel data streams across non-interfering portions of the channel. The use of OFDM allows to exploit effectively both, the frequency- and time- domain dimensions of a (slowly varying) frequency selective fading channel, as it admits simple single-tap equalization using a cyclic prefix (CP) (Hwang et al., 2009). Another important property is the possibility of implementing the required transmit signal processing operations completely in the digital domain, using the well-known and efficient fast Fourier transform (FFT) algorithm. Nevertheless, the implementation of OFDM in practice comes with many challenges. For example, high levels of out-of-band power emissions may be generated in transmission (due to discontinuities in the time-domain OFDM signal) if no corrective measures are taken. Trying to alleviate this problem, in this paper we focus on the use of a novel interpolation method that we named correlated interpolation (CI), as a simple way to achieve a suitable spectrum shape of the output signal in a practical system implementation.

An OFDM signal is a sequence of OFDM symbols, each one consisting of a collection of modulated orthogonal subcarriers. Since the amplitudes and phases of the subcarriers are often statistically independent, OFDM symbols are considered independent as well. Due to this phenomenon, the concatenation of OFDM symbols introduces discontinuities in the corresponding time-domain signal, or equivalently, high levels of out-of-band power emissions are generated. The simplest solution to this problem consists in using filtering techniques that limit this undesired power leakage out-of-band; nevertheless, the main drawback of such an approach is a reduction in the effectiveness of the CP that is introduced (van de Beek and Berggren, 2009a). Mahmoud and Arslan (2008) proposed an interesting alternative to control this problem, which is basically based on implementing adaptive symbol transitions when generating the OFDM signal. However, the main drawback in this case is the necessity to update the transmitter signal processing on a per-symbol basis, increasing as a consequence the complexity of the system (i.e., its implementation requirements).

This work focus on the generation of OFDM signals that comply with common quality requirements specified by emission masks, such as the ones that are presented in LTE specifications (3GPP 36.211, 2012). In other words, we address the generation of an oversampled digital signal that improves the spectral allocation of power in both, in-band and out-of-band regions. The key idea behind our proposal is based on correlating two OFDM symbols in frequency-domain, to obtain a continuous behavior in time-domain by means of the CI. This is obtained using an inverse discrete Fourier transform (IDFT) representation with double length, when compared to the one that is used in conventional OFDM transmissions. At that time, the concatenation points in the resulting sequence of OFDM symbols will show a continuous behavior in half of the original merging points (that correspond to the whole OFDM transmission). A preliminary version of this merging concept was previously presented in Lizarraga et al. (2011a).

A spectral precoding technique is additionally introduced in this paper, with the ultimate goal of guaranteeing a continuous time-domain behavior in the remaining concatenation points of the OFDM signal (i.e., in those time instants where the previously introduced CI processing is not able to provide a smooth transition). Unlike the scheme given by Mahmoud and Arslan (2008), this procedure consists of a static precoding, and does not necessitate extra spectral resources for its implementation like in alternative proposals (e.g., Chen et al., 2011). Simulation results show an important reduction in the out-of-band power emissions, when compared to alternatives previously reported in the literature (e.g., van de Beek and Berggren, 2009b). Moreover, a lower distortion in frequency domain is observed, when using the error-vector magnitude (EVM) as performance measure to assess this effect.

The rest of this paper is organized as follows. Section II presents the OFDM signal equations that are used throughout this work. The fundamentals of both, CI scheme and reduced distortion precoding proposal are presented in Section III and Section IV, respectively. Section V covers the channel modeling and the system operation at the receiver side, while simulation results are given in Section VI. Finally, Section VII shows the main conclusions.

Notation: Vectors are indicated in bold and lower-case letters, unlike complex or real scalar numbers. Matrices are indicated in bold with capital letters. Superscripts {}T, {}H, {}-1 represent transpose, Hermitian, and inversion, respectively. Notation 0M×N indicates the M×N matrix with all-zero entries, while IM stands for the M×M identity matrix.

II. CONVENTIONAL SIGNAL MODELING

Let us define s(t) as a continuous time OFDM signal composed by the concatenation of successive OFDM symbols si(t), i.e.,

(1)

Thus, the i-th symbol si(t) has a duration T, and derives from the sum of K orthogonal subcarriers belonging to the so-called subcarrier allocation scheme i.e.,

(2)

The amplitude and phase of each subcarrier is determined by the k-th complex symbol di,k (obtained from a prior complex modulator). Complex symbols are grouped into the K-element column vector di. The K-element column vector p(t) contains at each entry the corresponding subcarrier signal for . The interval is defined as [-Tg,Ts), where Tg and Ts represent the duration of the CP and the symbol period, respectively (with T=Tg+Ts). Note that the total symbol rate is determined by the relation K/T.

A sampled version of the i-th symbol without CP can be written as where F={Fq,k} represents the K-point discrete Fourier transform (DFT) by means of a matrix with elements with both q,k=0,..., K-1. Then, this matrix-vector approach can be with elements replaced by an inverse fast Fourier transform (IFFT) calculation (to take advantage of well-known efficient algorithms), i.e.,

(3)

where IFFT(a,b) represents the IFFT calculation with b points applied on the a vector. As stated by balanced subcarrier allocation schemes defined in several standards, we consider the case

(4)

which is included in 3GPP 36.104 (2012). To comply with this subcarrier mapping, can be applied as argument of the IFFT operation, considering and assuming that CSH(a,b) represents a circular up-shift operation of b positions in the column vector a. This allows us to get the output signal

(5)

which is an expression equivalent to (3), but considering (4).

In applications such as software defined radio (SDR) and cognitive radio (CR) (Pawelczak et al., 2011), the bandwidth that is considered to generate the output signal is allowed to be larger than the bandwidth that is assigned to the in-band channel. The aim in this case is to have the power emissions adequately accommodated within certain (and possibly variable) emission mask. To achieve this goal, an oversampled discrete-time output signal is required. Oversampling is related to interpolation, and it does not demand extra bandwidth resources for communication purposes. Let us define η as the oversampling factor. Then, following a similar approach to the one presented in (5), we are able to obtain an oversampled discrete-time output signal from

(6)

where can be seen as a conventional zero-padded version of , defined by

(7)

At this stage, attains the form of a vector with (K+1)η entries, containing the oversampled output signal. The classical CP insertion can be performed by the extension , where represents the CP fraction (i.e., the guard interval fraction). An OFDM symbol with CP is then obtained by means of concatenation, and is given by The corresponding non-oversampled expression results when setting η = 1. Here, we emphasize that the concatenation between the original signal and its cyclic extension generates a continuous signal at the merging point. This follows from the cyclo-stationary property of the IDFT, which shows continuity between the first and last point of the Fourier anti-transformed signal (see Fig. 1 in Lizarraga et al., 2011a).


Figure 1: Block diagram of the proposed OFDM transmitter. The correlation block supports the double-length symbols. Spectral precoding is applied to obtain additional improvement in the performance of the spectral shaping method.

III. CORRELATED INTERPOLATION

The IFFT calculation of a conveniently zero-padded vector produces a continuous interpolated signal, which still shows a smooth transition between the first and the last point. Based on this observation, the avoidance of discontinuities using a convenient IFFT processing when generating a sequence of OFDM symbols arose for first time. As previously stated, our goal is to exploit this property to reduce out-of-band power emissions in an OFDM-based system.

In the first instance we consider an OFDM signal that has no cyclic prefix, by setting the value υ=0. In this situation, it is valid to take into account the concatenation point between successive symbols with indexes i and i+1, for even values in i. Then, our proposal is based on obtaining a discrete-time signal that is equivalent to the one presented in (5) for and but that is jointly obtained from a single IFFT calculation that uses a double number of points. To achieve this representation, we start by defining

(8)

where and F represents the corresponding (K+1)-point DFT matrix. Then, we define the interlacing matrix

(9)

to produce the vector that contains (K+1) entries. With these definitions it is verified that

(10)

If we now focus on the corresponding interpolated representation, the expression originates from

(11)

which is a conventional zero-padded version of Consequently, the analyzed concatenation points follow a continuous signal. Notice that some latency has been introduced in this approach. Since the OFDM symbol has a double length, when assuming that the bandwidth remains unchanged (or equivalently, that the sampling rate in the receiver is not modified), we obtain that the transmission interval may span two channel states. This is because the length of a conventional OFDM symbol is determined according to the stationarity of the channel. Meanwhile, if a frequency-selective fading channel is considered, the insertion of CP is necessary (i.e., v> 0 should be used). Given the equivalence in (10), it is possible to insert prefixes for and in the conventional manner. However, even though the insertion of a prefix in does not affect the characteristics of our proposed scheme, the insertion of a prefix before does alter the continuous behavior that was achieved in To address this problem, we define a cyclic suffix (CS) that is attached to . An additional difficulty appears when considering the double-length interpolation. Since both, CP and CS must be extensions extracted from the single-length symbols, it is possible to show that an extension extracted from implies the concatenation of signals that have no continuity. To overcome this difficulty, we consider the interpolation difference, first introduced by Lizarraga et al. (2011a). It is possible to observe that at the time instants closer to the middle part of the single-length periods, both interpolations (i.e., the double-length and single-length interpolations) are similar. With these results, it is possible to express the output signal as

(12)

In this construction, the cyclic extensions are already included to deal with the channel dispersion in time-domain. Later on, we will discuss the effect that these modifications have in the operation of a conventional (i.e., legacy) OFDM receiver. A block diagram of the proposed transmitter is presented in Fig. 1. In the following section, a complementary proposal based on spectral precoding is presented. Spectral precoding can be used in conjunction with CI, to further improve the performance of the OFDM-based system in the frequency domain.

IV. SPECTRAL PRECODING

According to the results presented in Lizarraga et al. (2011a, b), it is verified that the previous development fosters a performance improvement in terms of out-of-band power emissions. This benefit appears due to the avoidance of discontinuities between consecutive OFDM symbols, with indexes i and i+1 for i even. After this, the interest is transferred to the rest of the concatenation points of the OFDM symbol sequence, that have indexes i and i+1 for i odd (i.e., the concatenation between successive double-length OFDM symbols). Taking into account this situation, we are now ready to develop a precoding scheme that avoids the discontinuities that (might) still exist in the output signal. However, it is important to emphasize that the application of the CI can improve the performance of many other proposals for OFDM system optimization. Therefore, this approach provides a general method that can be applied for other purposes, beyond the out-of-band power emission involved here.

Our development follows a similar reasoning to the one presented by van de Beek and Berggren (2009b), and considers to have a continuous signal in the concatenation of two OFDM symbols with indexes i and i+1 for i odd. When the signal and its first N derivatives are considered, the idea is to verify

(13)

for n =0,...,N. Then, a precoder can be designed using a non-injective linear transformation which, however, allows an acceptable decoding complexity. To express the stated condition in matrix form, we introduce the row vector with entries (-K/2,...,K/2), and a matrix with dimension (N+1)× (K+1), i.e.,

(14)

where superscripts 0,1,...,N in the elements indicate element-wise powers. After that, we define two (K+1)× (K+1) diagonal matrices as

(15)

with Then, the vector brings together two single-length frequency-domain OFDM symbols circularly rotated, and is defined for even values of i. Consequently, considering

(16)

the 2(N+1)×2(K+1) matrix defined by

(17)

can express the constraint in (13) as This restriction holds for a pair of OFDM symbols slightly distorted in frequency domain. This distorted vector is represented as where the vector wi with 2(K+1) elements identifies the distortion term. So, the precoding matrix can be now written as

(18)

This matrix is used to project the non-distorted vector vi onto the null space of yielding distorted vector

A. High-Order Continuity Restrictions

The expression given for in (12) is achieved replacing two fractions of the signal by the adequate fractions of and . Since the interpolations in and are not equal to the corresponding interpolations in (though they are similar as K grows), the concatenation of the fractions taken from and with the fraction of does not guarantee a strictly continuous behavior in the signal. However, due to the aforementioned similarities, the discontinuities that were found do not have large magnitude.

Meanwhile, it has been observed that by specifying continuity orders greater than zero in the precoder design (i.e., N > 0), these discontinuities must be effectively kept under control. This phenomenon can be taken into account at the precoder, defining a pair of new restrictions for time instants An interesting analysis is derived from (10); based on this equality, the new specification can be omitted for n=0. Additionally, simulation results justify the observation that due to interpolation similarities, these additional restrictions can be set only for n=1 regardless of the value in N. That is, focusing solely on the first derivative constraint, we are able to control completely the out-of-band power emissions that are generated in these instants (regardless of higher order constraints that may be defined for the complete precoder). These simplifications allow to reduce the distortion in frequency domain, as well as reduce the complexity of calculating the precoding matrix G. Yet, this operation can be performed even offline.

In this context, we develop some additional matrices to control double-length and simple-length interpolations, i.e.,

(19)
(20)

and two diagonal matrices

(21)

given the angles and Then, two new restrictions matrices can be established according to

(22)

Finally, the precoding matrix can be computed using .

B. Distortion Evaluation

The evaluation of our proposal considers that the decoding is carried out in conventional (legacy) OFDM receiver. Based on simulation results, we will show that no output filter will be required to keep the power emissions within pre-defined power masks (note that this is the main advantage behind our proposal). Nevertheless, the price to pay for the proposed precoding operation is a slight distortion in the frequency domain. To evaluate this phenomenon in quantitative way, the observed levels of distortion are assessed using the performance measure defined by 3GPP 36.104 (2012), i.e.

(23)

It is important to notice that for the precoder presented in (18), the latter expression reduces to

(24)

This value represents distortion which is significantly lower than the one achieved by van de Beek and Berggren (2009b), since the equivalent number of constraints in our precoder is half the number that was previously treated. This enhancement enables the use of higher-order complex modulation schemes, with the corresponding improvement in the achievable data rates of the system.

Accordingly, since in (22) two new constraints are added for the N >0 case the EVM expression becomes

(25)

This result indicates that distortion is increased; nevertheless, it holds a value significantly lower than the one reported in the work by van de Beek and Berggren (2009b). These relationships are shown in Fig. 2 for different values of N.


Figure 2: Error-vector magnitude for various continuity constraints N. Continuous lines with filled markers indicate the obtained distortion for our proposed approach, while dashed lines with void markers show the performance for single-length precoding.

V. RECEIVER OPERATION

The obtained signal has a fundamental difference when compared to conventional OFDM transmission: the conformation of the signal adjusts to scheme that inserts CP before and appends a CS after (where i takes even integer values). At this stage, it is important to highlight three observations. The first one relates to the processing introduced to achieve the signal , which is performed at the transmitter. The second one arises from considering that this processing is transparent to the receiver, assuming conventional sampling rate and proper timing. The third observation is that the receiver must perform both, CP extraction and CS extraction, treating the received OFDM sequence containing double-length symbols in the same way as the sequence would be treated if it contains single-length symbols. Alternatively, we can focus on the operation of the receiver and analyze that for a minimum sampling rate, vector

(26)

is obtained from . Then, we divide this vector into two parts of equal length, i.e., Note that and reduce to and when eliminating the CP and CS, respectively. However, if only CP extraction is applied on both sequences, turns into , while turns into

where

(27)

An important result emerges from previous expression: the extraction of CP in symbols with annexed suffix produces a distortion in the frequency domain, which can be interpreted as a rotation in the phase of each subcarrier. Consequently, if this distortion is taken as part of the effect of the channel, we conclude that this is compensated by the equalizer of a conventional receiver. This important result justify the use of our technique in system with conventional receivers.

VI. SIMULATION RESULTS

The proposed system is evaluated by the aid of numerical simulations. The first analysis uses a configuration with K=600 subcarriers, allocating to each subcarrier an equal fraction of the total average transmit power (fixed to 46 [dBm]) when using 16-QAM modulation. The timing parameters are set to Ts = 1/15 [ms] and Tg=9/128Ts. All these settings are in accordance to 3GPP UTRA/LTE specifications in 3GPP 36.104 (2012). The OFDM transmission is implemented in both, traditional way (i.e., using single-length symbols) and using our new proposal (i.e., implementing double-length symbols). Power spectra are estimated by means of Welch's averaged periodigram method, with a 8192-sample Hanning window and a 1024-sample overlap for an observation time of 1 [s].

The power spectral density for each order of continuity is shown in Fig. 3. The results show a reduction in out-of-band power, keeping a constant gap in the area of adjacent channels when compared to the conventional system. It is worth noticing that this feature is manifested for both, the precoding-free scheme and the scheme with precoding for N =0, 1, 2, 3, 4. It is observed that, at 5 [MHz] from the central frequency of the assigned band, the requirements of the power mask in presence of single-length OFDM transmission are verified using a continuity order N=2. On the other hand, when implementing double-length OFDM transmission, the requirements of the mask are already satisfied with unitary continuity order (i.e., for N=1). This analysis shows that the proposed technique reduces the out-of-band power emission, at the expense of introducing a low distortion (see Fig. 2).


Figure 3: Power spectral density. Filled markers indicate the CI performance, while void markers show the performance for previous schemes presented by van de Beek and Berggren (2009b).

The reduction in out-of-band power emissions that is obtained by using CI presents a constant gap. This result allows a significant reduction in the out-of-band emissions for frequencies near the limits of the allocated band (but outside this band). This behavior is difficult to achieve by conventional filtering. In this context, we analyze the reduction gap that is observed for both, different numbers of subcarriers and various orders of continuity at the precoder. These results are shown in the Fig. 4, where gaps values in the order of 3 [dB] are observed in most cases. We note that a smaller gap was observed when the precoder targets high orders of continuity in presence of a reduced number of subcarriers (e.g., when N = 4 and K=300).


Figure 4: Reduction gap observed for different subcarriers numbers according to various continuity orders in the precoder.

The performance has also been evaluated in terms of symbol error rate, in presence of additive white Gaussian noise (AWGN) channels and Rayleigh fading channels. The first case is considered by assuming a complex-valued zero-mean Gaussian noise vector with covariance matrix and a channel frequency response In the second case, a frequency-selective fading channel is represented using a complex-valued diagonal matrix Hi, corresponding to a hi vector with υ(K+1)+1 elements (i.e., an stochastic vector with zero-mean unit variance complex Gaussian entries). Figure 5 shows the symbol error rate that is achieved for K=600 by both, a conventional single-length transmission and our proposed double-length transmission with/without spectralprecoding (SP). In all cases, a conventional OFDM receiver is used to decode the symbols. The degradation introduced by the precoding is mainly reflected in the error floors. Yet, the minor distortion that our proposed scheme introduces reduces this disturbance significantly in both, AWGN and Rayleigh channel cases.


Figure 5: Symbol error rate for a conventional single-length transmission (dashed lines), and the CI approach proposed in this work (continuous lines). Both, AWGN channel and quasi-stationary Rayleigh fading channels are considered for various distortion settings depending on the precoder order.

To conclude, we would like to mention that the peak-to-averagepowerratio (PAPR) of the system is not affected when using our proposed method. This is because the double-length IFFT calculation is applied on a previously correlated vector (in frequency domain). In this manner, the obtained PAPR is equivalent to the one that is observed when analyzing two independent single-length symbols. So, the number of subcarriers that could be added to produce a peak is only K, rather than 2K. Numerical simulation results were obtained to justify this interpretation; nevertheless, they are not included here since they do not provide additional insights to the main topics already discussed in this work.

VII. CONCLUSIONS

The generation of suitable orthogonal frequency division multiplexing (OFDM) signals was studied in this paper, using a practical spectrum shaping method that enables the mitigation of out-of-band power emissions. To achieve this goal, the use of a correlatedinterpolation (CI) of OFDM symbols was first studied. After that, the use of spectral precoding was also considered. Numerical simulations were carried out to evaluate the performance of our proposed spectrum shaping method. Important improvements have been observed in various performance measures, when compared to the performance results that are observed when using similar methods that were previously reported in the literature. Since our method generates the OFDM signal in a completely digital way, its implementation in practice can be done robustly, without new (and usually costly) hardware modifications at the transmitter side. Due to no modifications are required at the receiver side, the same legacy OFDM processing can be used.

REFERENCES
1. 3GPP 36.211, Evolved Universal Terrestrial Radio Access (EUTRA); Physical Channels and Modulation (Release 8), 3GPP TSG RAN TS 36.211 Std. v11.0.0, Sep. (2012).         [ Links ]
2. 3GPP 36.104, Base Station (BS) Radio Transmission and Reception (Release 8), 3GPP TSG RAN TS 36.104 Std. v11.2.0, Sep. (2012).         [ Links ]
3. Chen, H.M., W.C. Chen and C.D. Chung, "Spectrally Precoded OFDM and OFDMA with Cyclic Prefix and Unconstrained Guard Ratios," IEEE Trans. Wireless Commun., 10, 1416-1427 (2011).         [ Links ]
4. Hwang, T., C. Yang, G. Wu, S. Li and G.Y. Li, "OFDM and its Wireless Applications: A Survey," IEEE Trans. Veh. Technol., 58, 1673-1694 (2009).         [ Links ]
5. Lizarraga, E., A. Dowhuszko and V. Sauchelli, "Out- of-Band Power Reduction in OFDM Systems using Double-length Symbols," IEEE Latin-American Conf. on Commun., Belem, Para, Brazil, 1-6 (2011a).         [ Links ]
6. Lizarraga, E., V. Sauchelli and S. Gallina, "Efficient Implementation of Subcarrier Allocation Schemes in N-Continuous OFDM Transmissions," Reunión de Procesamiento de la Informacón y Control, Oro Verde, Entre Ríos, Argentina, 159-164 (2011b).         [ Links ]
7. Mahmoud, H. and H. Arslan, "Sidelobe Suppression in OFDM-based Spectrum Sharing Systems using Adaptive Symbol Transition," IEEE Commun. Lett., 12, 133-135 (2008).         [ Links ]
8. Pawelczak, P., K. Nolan, L. Doyle, S.W. Oh and D. Cabric, "Cognitive Radio: Ten Years of Experimentation and Development," IEEE Commun. Mag., 49, 90-100 (2011).         [ Links ]
9. van de Beek J. and F. Berggren, "N-continuous OFDM," IEEE Commun. Lett., 13, 1-3 (2009a).         [ Links ]
10. van de Beek, J. and F. Berggren, "EVM-constrained OFDM precoding for reduction of out-of-band emission," IEEEVeh. Tech. Conf. Fall, 1-5 (2009b).         [ Links ]

Received: May 5, 2012
Accepted: October 22, 2012
Recommended by Subject Editor: Gastón Schlotthauer, María Eugenia Torres and José Luis Figueroa. (special issue RPIC 2011)

Creative Commons License Todo el contenido de esta revista, excepto dónde está identificado, está bajo una Licencia Creative Commons