Abstract

Orthogonal frequency division multiplexing (OFDM) is a powerful modulation choice for wideband wireless communication systems. However, its high peak-to-average power ratio greatly limits the high power amplifier (HPA) power efficiency. Here, we present the design of an adaptive predistorter to compensate the distortion caused by the HPA. Specifically, we deal with the implementation issue of the proposed predistorter in Lee and de Figueiredo's work (2006). The performance improvement by predistorter is verified by both floating-point simulation and fixed-point simulation, where the latter includes the distortion effects from the hardware. The bit widths for OFDM signals, ADC, and DAC are evaluated, and the bit width of 10 is shown to be sufficient for the hardware design.

1. Introduction

Othogonal frequency division multiplexing (OFDM) has attracted a lot of attention from the modern wireless communication community, because of its several desirable features for high-speed data transmission. In OFDM, a broadband signal is broken down into multiple narrowband subcarriers and implemented efficiently by using the IFFT algorithm [1]. OFDM advantages include: lower intersymbol interference, efficient use of frequency/spectrum through the use of different modulation/coding across subcarriers, and superior narrowband interference suppression capabilities. However, we need to consider the practical hardware limitations of low-cost RF and mixed signal devices when designing OFDM systems for broadband wireless data transmission. One of them is the high power amplifier (HPA) linearity and dynamic range, since OFDM signals have higher peak to average power ratios (PAPRs) than other high-performance modulations, and thus extra care is required. One of the most promising approaches to the mitigation of the PAPR problem is to use a predistorter applied to the OFDM signal prior to its entry into HPA. Its purpose is to compensate the nonlinearity of the HPA and improve the system performance. Many researchers have been investigating OFDM predistorter schemes [2–5]. However, all of these techniques are based on a general approximation form for the nonlinear system, rather than exploiting specific forms gleaned from physical device considerations. Due to this reason, we proposed a closed-form predistorter represented by a few parameters [6].

In this paper, we will apply Rapp's SSPA model [7] for HPA devices and show the implementation plan of the predistorter which was introduced in [6]. Furthermore, this paper also provides the design of a tracking algorithm for the case in which the practical HPA is unknown and varying. Finally, the simulation results are presented to investigate the performance improvement from predistorter, and study the distortion effects caused by saturation, overflow, and quantization with different number of bit widths, since the bit width of OFDM baseband (OFDM BB) and DAC/ADC is limited by cost and design constraints in real systems.

2. System Description

As shown in the block diagram of Figure 1, the proposed OFDM predistorter is placed after the OFDM baseband (BB) block to compensate for the degradation function in HPA.

In Figure 1, π‘Ÿ, π‘ž, and 𝑒 are amplitudes of output of the BB, predistorter, and HPA block, respectively. From the normalized Rapp's SSPA model [7], we have 𝑒[π‘Ÿ]=π‘Ÿ(1+(π‘Ÿ/𝐴0)2𝑝)1/2𝑝,(1) where 𝑒 is the HPA output amplitude, 𝐴0>0 is the maximum output amplitude from HPA, and 𝑝>0 is the parameter which controls the smoothness of the transition from the linear to the saturating region. Please note that the phase distortion for HPA is very small and hence can be neglected [7], and thus, the predistorter mainly focuses on amplitude compensation. In order to compensate the nonlinearity of the HPA, using the predistorter, the HPA output 𝑒 is targeted to be linear to π‘Ÿ. Thus, π‘ž(π‘Ÿ)(1+(π‘ž(π‘Ÿ)/𝐴0)2𝑝)1/2𝑝=π‘Ÿ.(2) Thus, we can derive the following equation [6]: π‘žπ‘Ÿ(π‘Ÿ)=(1βˆ’(π‘Ÿ/𝐴0)2𝑝)1/2𝑝,π‘Ÿ<𝐴0.(3)

When π‘Ÿβ‰₯𝐴0, (3) has no solution, and π‘ž has to be clipped which will be explained in the following section. Figure 2 shows a compensation example of the predistorter [6]. The upper and lower lines pertain to the PD and SSPA model, respectively, and the solid line represents the compensated effect. It shows similar effect with soft envelop limiter.

3. Design Architecture

The closed form expression of predistorter output π‘ž is shown in the previous section. However, the information of HPA parameters 𝐴0 and 𝑝 are unknown and time varying in practice. Figure 3 shows the detailed design of the architecture of our predistorter for a real time-varying environment.

The major block in the OFDM predistorter is an LMS block to update 𝐴0 and 𝑝 estimations and π‘ž estimation block which will be explained next in detail. As for other supporting blocks, R/P (rectangular coordinates transfer to polar coordinates) and P/R (polar coordinates transfer to rectangular coordinates) are employed with a 12-stage cordic algorithm. A look-up table is used to store the precalculated values for the calculation of LMS update and π‘ž estimation.

3.1. LMS Update

The total period is divided into two stages: the training stage covers the start time and periodic pilot time, and the remaining time is the operation stage (see Figure 4).

During the training stage, π‘ž estimation block is off, that is, π‘ž is set to be equivalent to the input signal magnitude of π‘Ÿ, and correspondent HPA output 𝑒 is known. The goal during this stage is to track the solution of time-varying 𝐴0 and 𝑝 by LMS algorithm. We define mean square error as follows:

𝐽𝐴0ξ€ΈβŽ›βŽœβŽœβŽœβŽπ‘ž,𝑝=𝐸1+π‘ž/𝐴0ξ€Έ2𝑝1/2π‘βŽžβŽŸβŽŸβŽŸβŽ βˆ’π‘’2.(4)

From the appendix, we have the expressions of πœ•π½(𝐴0,𝑝)/πœ•π΄0 and πœ•π½(𝐴0,𝑝)/πœ•π‘. For LMS algorithm, expected value is replaced by instantaneous value. Therefore, 𝐴0=𝐴(𝑛+1)0(𝑛)βˆ’πœ‡1β‹…βŽ‘βŽ’βŽ’βŽ’βŽ’βŽ’βŽ£βŽ›βŽœβŽœβŽœβŽœβŽœβŽπ‘ž(𝑛)𝐴1+π‘ž(𝑛)/0(𝑛)2̂𝑝(𝑛)ξ‚Ά1/2̂𝑝(𝑛)βŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽ β‹…ξƒ©βˆ’π‘’(𝑛)π‘ž(𝑛)𝐴0ξƒͺ(𝑛)2̂𝑝(𝑛)+1β‹…βŽ›βŽœβŽœβŽœβŽœβŽœβŽ1𝐴1+(π‘ž(𝑛)/0(𝑛)2̂𝑝(𝑛)ξ‚Ά1/2̂𝑝(𝑛)+1⎞⎟⎟⎟⎟⎟⎠⎀βŽ₯βŽ₯βŽ₯βŽ₯βŽ₯⎦,̂𝑝(𝑛+1)=̂𝑝(𝑛)βˆ’πœ‡2β‹…βŽ§βŽͺβŽͺ⎨βŽͺβŽͺβŽ©βŽ›βŽœβŽœβŽœβŽœβŽœβŽπ‘ž(𝑛)𝐴1+π‘ž(𝑛)/0(𝑛)2̂𝑝(𝑛)ξ‚Ά1/2̂𝑝(𝑛)βŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽ β‹…βŽ›βŽœβŽœβŽœβŽœβŽœβŽβˆ’π‘’(𝑛)π‘ž(𝑛)𝐴1+π‘ž(𝑛)/0(𝑛)2̂𝑝(𝑛)ξ‚Ά1/2̂𝑝(𝑛)βŽžβŽŸβŽŸβŽŸβŽŸβŽŸβŽ β‹…βŽ‘βŽ’βŽ’βŽ’βŽ’βŽ£βŽ›βŽœβŽœβŽœβŽœβŽβˆ’1𝐴̂𝑝(𝑛)β‹…1+π‘ž(𝑛)/0(𝑛)2̂𝑝(𝑛)ξ‚ΆβŽžβŽŸβŽŸβŽŸβŽŸβŽ ξƒ©β‹…lnπ‘ž(𝑛)𝐴0ξƒͺ⋅(𝑛)π‘ž(𝑛)𝐴0ξƒͺ(𝑛)2̂𝑝(𝑛)+𝐴ln1+π‘ž(𝑛)/0(𝑛)2̂𝑝(𝑛)ξ‚Ά2̂𝑝(𝑛)2⎀βŽ₯βŽ₯βŽ₯βŽ₯⎦⎫βŽͺβŽͺ⎬βŽͺβŽͺ⎭,(5) where πœ‡1 and πœ‡2 are step factors which will be defined in the following section, and the initial settings are 𝐴0(0)=1 and 𝑝(0)=1.

3.2. π‘ž Estimation

Figure 4 shows that after 𝑁1 samples in training stage, LMS is turned on once per 𝑁2 samples to update the estimates of 𝐴0 and 𝑝, and turned off during the remaining time (operation stage). Within the operation stage, predistorter is on, and π‘ž estimation is calculated based on the LMS estimations of 𝐴0 and 𝑝. That is, ⎧βŽͺβŽͺ⎨βŽͺβŽͺβŽ©π‘ž(𝑛)=π‘Ÿ(𝑛)𝐴1βˆ’π‘Ÿ(𝑛)/0(𝑛)2̂𝑝(𝑛)ξ‚Ά1/2̂𝑝(𝑛)𝐴,π‘Ÿ(𝑛)<0π‘ž(𝑛),max𝐴,π‘Ÿ(𝑛)β‰₯0(𝑛).(6) Please note that when ξπ΄π‘Ÿ(𝑛)β‰₯0(𝑛), π‘ž(𝑛) is clipped to π‘žmax. In this paper, we set π‘žmax=8, which will be shown to be suitable in the numerical result session.

3.3. Look-up Table

Five sets of function results are required to be stored in the look-up table to calculate the following functions: 𝐹1(π‘₯,𝑦)=[π‘₯2𝑦/(1βˆ’π‘₯2𝑦)]1/2𝑦, 𝐹2(π‘₯)=ln(π‘₯), 𝐹3(π‘₯,𝑦)=ln[1+π‘₯2𝑦], 𝐹4(π‘₯,𝑦)=[π‘₯2𝑦/1+π‘₯2𝑦]1/2𝑦, and 𝐹5(π‘₯,𝑦)=π‘₯2𝑦/1+π‘₯2𝑦.

Nonlinear quantization is applied to save table space. As for 𝐹1(π‘₯,𝑦), 𝑦 is quantized to π‘š1/32 when 33β‰€π‘š1≀96; π‘₯ is quantized to π‘š2/32 when 1β‰€π‘š2≀24, π‘š2/64 when 49β‰€π‘š2≀54, π‘š2/128 when 109β‰€π‘š2≀123, or π‘š2/256 when 247β‰€π‘š2≀255. Thus, table I requires the size of [64βˆ—(24+6+15+9)]Γ—16=3456Γ—16. As for 𝐹2(π‘₯), π‘₯ is quantized to π‘š/512 when 1β‰€π‘šβ‰€60, π‘š/256 when 61β‰€π‘šβ‰€120, π‘š/64 when 31β‰€π‘šβ‰€512 (π‘₯ will be clipped as eight if it exceeds eight). Thus, table II requires the size of (60+60+482)Γ—16=602Γ—16. As for 𝐹3(π‘₯,𝑦), 𝑦 is quantized to π‘š1/32 when 33β‰€π‘š1≀96; π‘₯ is quantized to π‘š2/32 when 1β‰€π‘š2≀31, π‘š2/16 when 16β‰€π‘š2≀31, π‘š2/8 when 16β‰€π‘š2≀64. Thus, table III requires the size of [64βˆ—(31+16+49)]Γ—16=6144Γ—16. As for 𝐹4(π‘₯,𝑦), 𝑦 is quantized to π‘š1/32 when 33β‰€π‘š1≀96; π‘₯ is quantized to π‘š2/64 when 1β‰€π‘š2≀63, π‘š2/32 when 32β‰€π‘š2≀63, π‘š2/16 when 32β‰€π‘š2≀47, π‘š2/8 when 24β‰€π‘š2≀64. Thus, table IV requires the size of [64βˆ—(63+32+16+41)]Γ—16=9728Γ—16. As for 𝐹5(π‘₯,𝑦), 𝑦 is quantized to π‘š1/32 when 33β‰€π‘š1≀96; π‘₯ is quantized to π‘š2/32 when 1β‰€π‘š2≀63, π‘š2/8 when 16β‰€π‘š2≀64. Thus, table V requires the size of [64βˆ—(63+49)]Γ—16=6944Γ—16. The total size for all of these five look-up tables is 26874Γ—16 whose area is less than 0.5mm2 for CMOS18.

3.4. Complexity Evaluation

Based on the precalculated parameters, (5) can be expressed as 𝐴0𝐴(𝑛+1)=0(𝑛)βˆ’πœ‡1⋅𝐴0(𝑛)⋅𝐹4ξƒ©π‘ž(𝑛)𝐴0ξƒͺξƒ­(𝑛),̂𝑝(𝑛)βˆ’π‘’(𝑛)⋅𝐹4ξƒ©π‘ž(𝑛)𝐴0ξƒͺ(𝑛),̂𝑝(𝑛)⋅𝐹4ξƒ©π‘ž(𝑛)𝐴0,(𝑛),̂𝑝(𝑛)ξƒͺξƒ°(7)̂𝑝(𝑛+1)=̂𝑝(𝑛)βˆ’πœ‡2⋅𝐴0(𝑛)⋅𝐹4ξƒ©π‘ž(𝑛)𝐴0ξƒͺ⋅𝐴(𝑛),̂𝑝(𝑛)βˆ’π‘’(𝑛)0(𝑛)⋅𝐹4ξƒ©π‘ž(𝑛)𝐴0ξƒͺβ‹…βŽ‘βŽ’βŽ’βŽ£βˆ’πΉ(𝑛),̂𝑝(𝑛)4ξ‚€ξπ΄π‘ž(𝑛)/0(𝑛),̂𝑝(𝑛)̂𝑝(𝑛)⋅𝐹2ξƒ©π‘ž(𝑛)𝐴0(ξƒͺ+𝐹𝑛)3ξ‚€ξπ΄π‘ž(𝑛)/0(𝑛),̂𝑝(𝑛)2̂𝑝2⎀βŽ₯βŽ₯⎦⎫βŽͺ⎬βŽͺ⎭.(𝑛)(8) After HPA parameters 𝐴0 and 𝑝 estimation, the π‘ž is estimated by ξπ΄Μ‚π‘ž(𝑛)=0(𝑛)⋅𝐹1ξƒ©π‘ž(𝑛)𝐴0ξƒͺ(𝑛),̂𝑝(𝑛).(9) Therefore, the complexity in total includes 5 addition/subtracitons and 15 multiplications which are relative low.

4. Numerical Results

From (1), HPA results in a highly nonlinear situation with high input amplitude, and small distortions vice versa. Therefore, a relative level of power back off is required to reduce HPA distortion. Here, we define input back-off (IBO) as IBO=10log10𝐴20𝑃inξƒͺ,(10) where 𝑃in is input average power of OFDM signal. Next, we will perform the algorithmic level and hardware level (fixed-point) simulations, while the latter include all the distortion effects in hardware such as round-off error and coefficient quantization.

We set the simulation parameters as follows.

(i)16QAM modulated OFDM signal with 64-point IFFT are studied. (ii)The average input back-off power is 6.375 dB, if not being mentioned. (iii)The start training sequence is employed with a length of 𝑁1=160 training samples, and every training sample per 𝑁2=16 OFDM symbols is applied to the following sequence. (iv)𝐴0 and 𝑝 are both assumed to be Gaussian random numbers with mean of 1 and variance of 0.0025. (v)Step factor πœ‡1=1.5 in training stage and 0.5 in the operation stage, while πœ‡2 is set to be as much as six times of πœ‡1. (vi)The bit width of OFDM output 𝐼π‘₯,𝑄π‘₯, the bit width of DAC input 𝐼𝑦,𝑄𝑦, and ADC output are evaluated, since the former is limited by the area cost, and the latter is limited by DAC/ADC design. (vii)The channel is assumed as AWGN with variance of πœ‚0/2.

From Figure 5, it shows that the bit width of 10 is recommended for OFDM BB output, DAC input, and ADC output, since there is not much improvement to increase bit width beyond 10. The proposed implementation plan of predistorter is shown to improve system performance even including degradation effect from hardware.

5. Conclusion

In this paper, we have provided an implementation plan of the proposed predistorter in [6] to compensate the nonlinear distortion of SSPA. We used an LMS algorithm for time-varying environment, which we have shown to be capable of tracking SSPA parameters. Finally, a fixed-point simulation including hardware degradation factor was performed to verify the superior performance of the proposed implementation scheme.

Appendix

A. Derivative of 𝐽(𝐴0,𝑝) w.r.t 𝐴0

ξ€·π΄πœ•π½0ξ€Έ,π‘πœ•π΄0=ξ‚ƒπœ•πΈπ‘ž/(1+(π‘ž/𝐴0)2𝑝)1/2π‘ξ‚„βˆ’π‘’2πœ•π΄0π‘ž=2𝐸(1+(π‘ž/𝐴0)2𝑝)1/2𝑝ξƒͺβ‹…ξ‚΅βˆ’π‘žβˆ’π‘’ξ‚Άβ‹…ξ€Ί2𝑝1+(π‘ž/𝐴0)2π‘ξ€»βˆ’1/2π‘βˆ’1β‹…π‘ž2𝑝𝐴⋅(βˆ’2𝑝)β‹…0ξ€Έβˆ’2π‘βˆ’1ξƒ­π‘ž=2𝐸(1+(π‘ž/𝐴0)2𝑝)1/2𝑝ξƒͺβ‹…ξ‚΅π‘žβˆ’π‘’π΄0ξ‚Ά2𝑝+1⋅1(1+(π‘ž/𝐴0)2𝑝)1/2𝑝+1.ξƒͺξƒ­(A.1)

B. Derivative of 𝐽(𝐴0,𝑝) w.r.t 𝑝

ξ€·π΄πœ•π½0ξ€Έ,𝑝=ξ‚»πœ•π‘πœ•πΈ[π‘ž/(1+(π‘ž/𝐴0)2𝑝)1/2π‘βˆ’π‘’]2ξ‚Όπ‘žπœ•π‘=2𝐸(1+(π‘ž/𝐴0)2𝑝)1/2π‘ξƒ­β‹…βˆ’π‘’πœ•[π‘ž/(1+(π‘ž/𝐴0)2𝑝)1/2𝑝]ξƒ°π‘žπœ•π‘=2𝐸(1+(π‘ž/𝐴0)2𝑝)1/2π‘ξƒ­ξƒ¬ξ‚΅βˆ’1βˆ’π‘’β‹…π‘žβ‹…expξ‚Άξƒ©ξ‚΅π‘ž2𝑝⋅ln1+𝐴0ξ‚Ά2𝑝⋅ξƒͺξƒ­πœ•[(βˆ’1/2𝑝)β‹…ln(1+(π‘ž/𝐴0)2𝑝)]ξƒ°βŽ§βŽͺ⎨βŽͺβŽ©ξƒ¬π‘žπœ•π‘=2𝐸(1+(π‘ž/𝐴0)2𝑝)1/2π‘ξƒ­ξƒ¬ξ‚΅π‘žβˆ’π‘’β‹…π‘žβ‹…1+𝐴0ξ‚Ά2π‘ξƒ­βˆ’1/2π‘β‹…βˆ’1ξ‚Έξ‚΅ξ‚Άβ‹…2π‘πœ•[ln(1+(π‘ž/𝐴0)2𝑝)]ξƒ©ξ‚΅π‘žπœ•π‘+ln1+𝐴0ξ‚Ά2𝑝ξƒͺβ‹…πœ•(βˆ’1/2𝑝)⎧βŽͺ⎨βŽͺβŽ©ξƒ¬π‘žπœ•π‘ξƒ­ξƒ°=2𝐸(1+(π‘ž/𝐴0)2𝑝)1/2π‘ξƒ­ξƒ¬ξ‚΅π‘žβˆ’π‘’β‹…π‘žβ‹…1+𝐴0ξ‚Ά2π‘ξƒ­βˆ’1/2π‘β‹…βˆ’1ξ‚Έξ‚΅ξ‚Άβ‹…12𝑝(1+(π‘ž/𝐴0)2𝑝)β‹…ξ€·πœ•[π‘ž/𝐴0ξ€Έ2𝑝]ξƒ©ξ‚΅π‘žπœ•π‘+ln1+𝐴0ξ‚Ά2𝑝ξƒͺβ‹…12𝑝2⎧βŽͺ⎨βŽͺβŽ©ξƒ¬π‘žξƒ­ξƒ°=2𝐸(1+(π‘ž/𝐴0)2𝑝)1/2π‘ξƒ­ξƒ¬ξ‚΅π‘žβˆ’π‘’β‹…π‘žβ‹…1+𝐴0ξ‚Ά2π‘ξƒ­βˆ’1/2π‘β‹…βˆ’1ξ‚Έξ‚΅ξ‚Άβ‹…12𝑝(1+(π‘ž/𝐴0)2𝑝)β‹…πœ•ξ€Ίξ€·ξ€·exp2𝑝⋅lnπ‘ž/𝐴0ξ€Έξ€Έξ€»ξƒ©ξ‚΅π‘žπœ•π‘+ln1+𝐴0ξ‚Ά2𝑝ξƒͺβ‹…12𝑝2⎧βŽͺ⎨βŽͺβŽ©ξƒ¬π‘žξƒ­ξƒ°=2𝐸(1+(π‘ž/𝐴0)2𝑝)1/2π‘ξƒ­ξƒ¬ξ‚΅π‘žβˆ’π‘’β‹…π‘žβ‹…1+𝐴0ξ‚Ά2π‘ξƒ­βˆ’1/2π‘β‹…βˆ’1ξ‚Έξ‚΅ξ‚Άβ‹…π‘ž2𝑝(1+(π‘ž/𝐴0)2𝑝)ξ‚΅ξ‚΅π‘žβ‹…exp2𝑝⋅ln𝐴0β‹…πœ•ξ€·ξ‚Άξ‚Άξ€Ίξ€·2𝑝⋅lnπ‘ž/𝐴0ξ€Έξ€Έξ€»ξ‚€ξ€·πœ•π‘+ln1+π‘ž/𝐴0ξ€Έ2𝑝⋅12𝑝2π‘žξ‚Ήξ‚Ό=2𝐸(1+(π‘ž/𝐴0)2𝑝)1/2𝑝ξƒͺβ‹…ξƒ©π‘žβˆ’π‘’(1+(π‘ž/𝐴0)2𝑝)1/2𝑝ξƒͺβ‹…βˆ’1𝑝⋅(1+(π‘ž/𝐴0)2𝑝)ξƒͺξ‚΅π‘žβ‹…ln𝐴0ξ‚Άβ‹…ξ‚΅π‘žπ΄0ξ‚Ά2𝑝+ln(1+(π‘ž/𝐴0)2𝑝)2𝑝2.ξƒ­ξƒ°(B.1)