1. Technical Field
Embodiments of the present disclosure relate to signal processing, and more particularly to systems and devices and a method of crest factor reduction utilizing amplitude gain compression.
2. Description of Related Art
Modern communication systems use wideband transmission waveforms that contain multiple carriers and/or sub carriers. These transmission waveforms are not constant in amplitude, but can have very fast amplitude variations wherein the peak amplitude is substantially larger than the root mean square (rms) amplitude. The multi carrier signal is transmitted using a multi carrier power amplifier (MCPA) to transmit the signal at a very high efficiency, while maintaining acceptable signal quality and a high adjacent channel power ratio to meet the spectral emission mask requirements.
To maximize efficiency of a power amplifier (PA), it is desirable to transmit the signal near the saturation power levels PSAT. However, if the peak amplitude of the signal reaches PSAT then the signal will suffer distortion which may accumulate to cause high inter-modulation distortion. As the levels of inter-modulation distortion may be very high, the digital pre-distortion may fail to fully correct the amplitude distortion, and the inter-modulation (IMD) may be higher than that allowed by the required spectral emission mask.
To avoid this problem, the signal may be transmitted at lower levels to avoid the peak signal reaching PSAT. But this simple back-off approach results in lower PA efficiency, which is a key performance parameter.
The disclosure is illustrated by way of example and not by way of limitation in the accompanying drawings in which like references indicate similar elements. Various embodiments illustrate different features of the disclosure. It should be noted that references to “an” or “one” embodiment in this disclosure are not necessarily to the same embodiment, and such references can mean “at least one.” Descriptions of components in the embodiments are given for the purpose of illustrating rather than limiting.
An electronic device can comprise an MCPA transmission system. The proposed CFR method compresses the peak amplitudes of signals in the CFR processor 12. Embodiments of systems, devices, and CFR methods using amplitude compression are given. The disclosed systems and communication devices may be implemented as standalone devices, or integrated into various network gateway devices or network terminal devices. The various network gateway devices comprise base stations, bridges, routers, switches, or hot spots or access points for wireless networking. The network terminal devices comprise set-top boxes, cell phones, tablet personal computers, laptop computers, multimedia player, digital cameras, personal digital assistants, navigation devices, or mobile internet devices.
A multi-carrier combiner 11 combines multiple carriers on a common spectrum, wherein each carrier has an gain of γi at a frequency fi to produce a composite signal expressed as follows:
wherein γi and fi are the carrier gain and frequency of the individual carrier, respectively with a subscript of variable i, and Nc is a total number of the multiple carriers.
The process of CFR reduces the PAR or crest factor of the multi-carrier signal in the digital domain, and the resulting out-of-band IMD is then filtered by a digital filter.
A digital predistortion processor 13 performs digital predistortion on signals outputted by the CFR processor 12 to compensate for distortion which may be caused by a power amplifier 15 based on feedback provided by an analog feedback unit 17. The digital predistortion processor 13 transmits post predistortion signals to the power amplifier 15 through an analog transmission unit 14.
With reference to
Block 100 Up-Sampler:
The up-sampler 100 is operable to up-sample the multi-carrier signal x(t) having a sampling rate of Rs by a factor of M, to produce the output signal y(t) having the sampling rate M·Rs. This up-sampling has two main purposes:
(1) to prevent aliasing of the noise caused by the subsequent amplitude compression from being folded back into the bandwidth of the fundamental signals, which would increase the error vector magnitude (EVM). The inter-modulation distortion products that are not aliased will be filtered without increasing the EVM; and
(2) to prevent the signal from overshooting, such as the Gibb's phenomenon, at the signal transitions, which may increase the PAR. The higher the up-sampling factor M, the lower the resulting EVM, and the lower the PAR which can be achieved for a given required minimum EVM.
The signal x(t) will typically have a sampling rate Rs that is slightly higher than the Nyquist sampling rate. The up-sampling by M will generate the signal y(t) having a sampling rate of M·Rs. It is desirable to have the sampling rate M·Rs to be about 2 to 3 times the signal bandwidth of the signal x(t) in order to reduce the undesired IMD from folding back to the fundamental signal bandwidth.
This implies that the higher the bandwidth the lower the signal degradation due to aliasing. However, increasing the bandwidth increases the complexity, especially in the CFR filter 600.
Block 200 Compute Input Amplitude
The absolute value processor 200 is designed to compute the amplitude of the input signal:
a(t)=|x(t)|=√{square root over (Ix(t)2+Qx(t)2)}{square root over (Ix(t)2+Qx(t)2)} (1)
the variables Ix and QX are respectively the I and Q components of the signal x(t).
Block 300 amplitude compression processing:
The special amplitude compressor 300 is operable to suppress the amplitude in a fashion that minimizes the in-band distortion, and also minimizes the sample regrowth in the subsequent filter process.
The amplitude compressor 300 executes an amplitude compression method that computes the compression gain g1(t) based on the amplitude a(t),
The method of computing g1(t) based on a(t) and generating the desired amplitude is shown in Table 1, wherein the amplitude signal a(t) is abbreviated as “a”, and “c” represents a desired amplitude of the signal a(t):
In Zone Z1, because the amplitude and the phase are unchanged, no distortion will occur, there is no EVM degradation for this section of the transfer function, and the PAR is unchanged. Referring to
In Zone Z2, a very small amplitude range is limited to T1, and some low-level distortion will occur. If T2 and T1 are kept close to each other, very small changes of the EVM and PAR result. Referring to
In Zone Z3, the large signal is compressed in a gradually heavier fashion, that is, the larger the signal amplitude, the heavier the compression which is applied. Referring to
The compression gain g1(a) is a function of the amplitude a(t) and can be generalized by the formulation of a non-linear function that will minimize the PAR for a given EVM, and this function can be searched for the optimum non-linear gain function (or profile) that will be unity gain in Zone 1, about 1/a gain in Zone 2, and increasingly hard compression in Zone 3. The function g1(a) can be modeled as a function that yields the desired EVM, and have the lowest PAR. The function g1(a) can be represented as a polynomial or a polynomial ratio:
The variables p(a) and q(a) are two polynomial functions of amplitude “a” with order m−1 and n−1 respectively, and m, n, pi and qi are the polynomial coefficients selected to yield the desired EVM with the lowest PAR.
Alternatively, the function g1(t)=ƒ(a(t)) may be implemented as a look-up table (LUT). This gain function or profile can be stored in an LUT for use in the actual allocation of the gain compression.
The LUT may comprise L pairs of amplitude and gain points according to the function g1(t)=ƒ(a(t)), and with different scale methods as summarized in Table 2.
As a possible illustration without loss of generality, for Option 1, the method to process the amplitude gain compression is as follows:
Step 1: initiate with L1 amplitude thresholds Ti, wherein the “i” is an integer variable ranging from 1 to L1, that is i=1, 2, . . . , or L1 in decibels (dB) above the RMS of the signal level, wherein
T1<T2<T3<T4< . . . <TL1 (3)
These amplitude thresholds do not have to be equally spaced, and T1 and TL1 are respectively the smallest and largest possible values of a(t) in dB.
Step 2: select the gain set gi in dB corresponding to Ti as shown in
Step 3: perform interpolation of the L1 pairs {(Ti,gi)} using linear interpolation, log-based or loglog-based interpolation, or any other non-linear interpolation method to create a larger set of amplitude gain pairs {(aai, gg1)}, i=1:L2 (L2 is the number of points in this interpolation) wherein aai is one of a plurality of equally spaced signal amplitudes with small step sizes in
Step 4: convert aai and gg1, into a linear scale:
Step 5: create L3 samples in the LUT. The interpolation is repeated to provide {(aaai,gggi)}, where i=1:L3, in which aaai covers the range of the signal amplitude and are uniformly spaced with very small step sizes. Typically L3 minimizes the quantization noise by about 1000, 0=aaa1>aaa2>aaa3 . . . > . . . >, aaaL2=Max (amplitude), (aaai+1−aaai)=Δaaa=Max (amplitude)/(L3−1), and gggi is a gain value being the interpolated gain value in the interpolation process.
Step 6: store {(aaai,gggi)} in an LUT, referred to as a compression gain LUT, where aaai is used to address the LUT and gggi are gain values outputted by the LUT.
Step 7: for a value a(t), the processor finds the closest location ai that is smaller than a(t) in the compression gain LUT.
Step 8: based on the LUT, the amplitude compression processor computes the value g1(t) via linear interpolation, such as:
wherein Δa=(ai+1−ai) (7)
Step 9: apply g1(t) for input amplitude a(t).
Table 2 shows 3 different gain compression options. For Option 1, the conversions (4) and (5) in Step 4 are skipped. For Option 2, conversion (4) in Step 4 is skipped. For Option 3 conversion (5) in Step 4 is skipped.
To obtain the best performance of the amplitude compression for CFR, the value needs to be determined empirically. One method is by searching by adjusting Ti and gi until the lowest PAR at complementary cumulative distribution function (CCDF)=10−6 is achieved with a given EVM.
Block 400 Amplitude Multiplier:
The multiplier 400 multiplies gain g1(t) with the gain correction g2(t) to generate gain g(t) according to equation (8), to perform signal amplitude compression:
g(t)=g1(t)·g2(t) (8)
wherein g1(t) is the gain value extracted from the gain compression LUT, and g2(t) is the correction gain value (>1), to compensate for the signal rms power loss due to peak amplitude compression.
Block 500 Signal Multiplier:
The signal multiplier obtains signal z(t) by multiplying signal y(t) with the gain g(t) to perform signal amplitude compression
zI(t)=g(t)yI(t) (9)
zQ(t)=g(t)yQ(t) (10)
A signal in the formulas (9) and (10), such as one of zI(t), zQ(t), yQ(t), and yI(t), have subscripts I and Q which respectively represent I and Q components of the signal, such as z(t) or y(t). This signal compression reduces the I and Q samples with the same gain g(t), thus keeping the phase constant, and only suppresses the high amplitude signals as defined by the gain LUT profile. This eliminates the phase distortion in the peak limiting process, and consequently the EVM is not reduced by this factor.
Block 600 Crest Factor Reduction Filter:
The CFR filter 600 suppresses the undesired IMD outside of the fundamental carriers. This process can be accomplished by various approaches.
First approach: composite multi carrier filter.
The first approach comprises down-sampling the signal z(t) to a lower sampling rate, performing multi-carrier finite impulse response (FIR) filtering on the down-sampled signal, and then up-sampling the filtered signal to the desired sampling rate to generate a signal u(t). The block diagram of this filter is shown in
One disadvantage of this filter is that if the signal is decimated by a large ratio then there are large signal overshoots in the multi-carrier filter output because the ratio of the sampling rate over the occupied bandwidth is low. This would result in a high PAR. However, if the down sampling ratio K is small, then the complexity of the multi-carrier filter is large in terms of the required number of coefficients. A second approach outlined below can help with this compromise.
Second approach: multi carrier sub-band CFR filter.
The signal z(t) is the composite of many carriers, each at different center frequency, bandwidth, and different transition bandwidth requirements, as tabulated in Table 3.
The spectrum of the multi-carrier signal z(t) is as shown in
Referring to
that is slightly larger than the Nyquist bandwidth. The down-sampler may comprise one of the down samplers 6201, 6202, . . . and 620m.
The carrier is filtered by a low complexity FIR filter with Ni taps. Ni is computed approximately as:
wherein
Rs is the input signal sampling rate of the composite signal x(t)
M is the interpolation factor in Block 100
Ki is the decimation factor in Block 620, where i is an index
Δƒi is the required transition bandwidth of the ith carrier
R is the required filter rejection sidelobe in dB.
The FIR filter may comprise one of carrier filters 6301, 6302, . . . and 630m.
After carrier filtering, the signals are then subject to an amplitude, phase, and time (APT) correction by an APT processor, such as one of blocks 6401, 6402, . . . and 640m. Each of the blocks 6401, 6402, . . . and 640m comprises a complex multiplier providing amplitude adjustment, a phase correction device, and a delay buffer providing time delay adjustment. The delay adjustment can be performed with a delay buffer at an up-sampler, such as 6501, 6502, . . . or 650m to output higher resolution signals. The outputs of the blocks 6401, 6402, . . . and 640m are respectively up-sampled to the rate M·Rs by up samplers 6501, 6502, . . . and 650m, and are then respectively frequency shifted by frequency shifters 6601, 6602, . . . and 660m back to the desired frequency fi. The required NCOs at the front and back of the CFR filter are synchronized to avoid unwanted time and phase offsets. After converting to the correct frequencies, a sub-band combiner 680B recombines the carriers outputs by frequency shifters 6601, 6602, . . . and 660m to create the desired signals that have the inter-modulation energy outside of the required signal bands reduced.
Block 700 Compute Output Amplitude.
Another absolute value processor 700 computes the amplitude b(t) of the output signal u(t)
b(t)=|u(t)|=√{square root over (Iu(t)2+Qu(t)2)}{square root over (Iu(t)2+Qu(t)2)} (12)
The Iu(t) and Qu(t) are respectively the I and Q components of the signal u(t).
Block 800 Gain Correction.
A gain correction unit 800 produces the compensation gain correction to maintain a constant power gain through the CFR processor 12. The peak amplitude compression by the multipliers 500 produces a power reduction of the signal y(t). The gain correction unit 800 produces gain correction factor g2(t) to adjust the average power of the output signal u(t) to be substantially the same as that of the input signal y(t). The process of this function is shown in
The process starts by an integration and dump (ID) unit 810 computing the square of the amplitude ai2, and an integration and dump (ID) unit 811 computing bi2. The ID unit 810 obtains a sum of the squared amplitudes Py=Σ1Nai2 over N samples, and the ID unit 811 obtains a sum of the squared amplitudes Pu=Σ1Nbi2 over N samples. The ID units 810 and 811 then feed the Py and Pu to a divider 820 to compute the short term power ratio
The short term power ratios are then passed through a filter 830, such as a filter of ID, infinite impulse response (IIR), or finite impulse response (FIR), to compute the average power ratio Gm. Finally, a square root unit 840 obtains the square root of the power ratio Gm to determine the gain correction factor g2(t).
The gain correction factor g2(t) is multiplied with the compression gain g1(t) to produce the instantaneous amplitude gain g(t).
As shown in
Note that for the same peak PAR6, the amplitude compression technique has a higher SNR compared to the amplitude clipping technique.
Referring to
Referring to
The foregoing disclosure of various embodiments has been presented for purposes of illustration and description. It is not intended to be exhaustive or to limit the disclosure to the precise forms disclosed. Many variations and modifications of the embodiments described herein will be apparent to one of ordinary skill in the art in light of the above disclosure. The scope of the disclosure is to be defined only by the claims appended hereto and their equivalents.
Number | Name | Date | Kind |
---|---|---|---|
7239257 | Alexander et al. | Jul 2007 | B1 |
7308041 | Schrader et al. | Dec 2007 | B1 |
20030108214 | Brennan et al. | Jun 2003 | A1 |
20040218689 | Akhtman | Nov 2004 | A1 |
20050181754 | Wu et al. | Aug 2005 | A1 |
20080095284 | Hori et al. | Apr 2008 | A1 |
20090180527 | Asami | Jul 2009 | A1 |
20120281786 | Lindemann et al. | Nov 2012 | A1 |
20140192913 | Bogdan | Jul 2014 | A1 |