5.8  The Gap Approximation

The next sections discuss the so-called gap approximation applied to PAM and QAM.

5.8.1  The PAM 6-dB rule: each extra bit requires 6 dB

PIC
Figure 5.16: a) Graph of Eq. (5.17) for σ2 = 2 and d = 1. b) Corresponding Pe.

An interesting relation obtained from Eq. (5.8) is

b = log2M = log2 (12Ec d2 + 1) = 1 2log2 (12SNRσ2 d2 + 1).
(5.17)

Figure 5.16 a) presents a graph of Eq. (5.17) when varying SNR for σ2 = 2 and d = 1. It indicates that, when b is sufficiently high, each additional bit (from 4 to 5 bits, for example) requires an increase of 6 dB in SNR if one desires to maintain the Pe. Another way of observing this fact is to assume that 12SNRσ2 d2 » 1 and obtain

b ≈ 1 2log2 (12SNRσ2 d2 ) = 1 2log2 (SNR) + cte. = 1 2log2 (100.1SNRdB) + cte.,

which can be expressed as Δb ≈ 0.1661ΔSNRdB or ΔSNRdB ≈ 6.0206Δb.

Figure 5.16 b) indicates how Pe varies with b, which was allowed to assume non-integer values. Because σ2 and the separation d among symbols is kept constant, the parcel Q (d∕(2σ)) is constant and, in this case, approximately 0.36. Hence, Pe tends to 2 × 0.36 as b increases.

It is important to understand that, in spite of Eq. (5.17) looking like the “capacity” equation, this expression does not say anything about probability of error nor if this number b of bits is below or above the capacity C¯ = 1 2 log2 (1 + SNR) (Eq. (4.53), repeated here for convenience) of a discrete-time AWGN channel.

We are interested in understanding the gap approximation for uncoded PAM and QAM:

b¯ = 1 2log2 (1 + SNR Γ ).
(5.18)

This approximation is widely used and the reason is the following task: say one has an estimated SNR (or, equivalently, Ec and σ2) and wishes to find b to operate at a given Pe. In other words, some function b = f(Pe,SNR) is needed.

In this situation, one does not know b, M or dmin that allows to achieve the target Pe. Hence, Eq. (5.17) cannot be used. The capacity

c¯ = 1 2log2 (1 + SNR)

only tells the theoretical (assuming the possibility of infinite delay and complexity) maximum value of b¯, but does not indicate the probability of error when using an uncoded PAM and QAM (or any other code). However, one can use expressions for Pe that take into account SNR and b. In the following paragraph, PAM is assumed, but a similar reasoning applies to QAM.

For PAM, one has

Pe = 2 (1 − 1 M )Q ( 3 M2 − 1SNR).

However, because M (and, equivalently, b = log2M) cannot be expressed in closed form, it is not possible to find an “inverse” b = f(Pe,SNR). Two alternatives are to plot the expression for Pe or even find M numerically. More commonly, one adopts the “gap approximation” that (in this case) consists in ignoring the term (2∕M)Q ( 3 M2−1SNR) and adopt:

PeP^e = 2Q ( 3 M2 − 1SNR).
(5.19)

Figure 5.17 illustrates the approximation for SNR = 10 dB, where it can be seen how the discarded term PeP^e = (2∕M)Q ( 3 M2−1SNR) varies with b.

PIC
Figure 5.17: Comparison between the symbol error probability Pe and the estimate P^e using the gap approximation for PAM with SNR = 10 dB.

Figure 5.18 is similar to Figure 5.17 but adopts SNR = 30 dB. Note that the gap approximation improves when the SNR increases.

PIC
Figure 5.18: Comparison between the symbol error probability Pe and the estimate P^e using the gap approximation for PAM with SNR = 30 dB.

Assuming the gap approximation of Eq. (5.19), to find M that corresponds to the desired Pe, one can now simply isolate M:

Q−1 (Pe 2 ) 3 M2 − 1SNR

M = 2b1 + 3 [Q−1 (Pe 2 )] 2SNR

and write

b = 1 2log2 (1 + SNR Γ ),

where

ΓPAM = 1 3 [Q−1 (Pe 2 )]2
(5.20)

is the gap for PAM in linear scale.

5.8.2  The QAM 3-dB rule: each extra bit requires 3 dB

An expression can be derived for the number of bits per dimension of a SQ QAM:

b¯ = 1 2log2 (12Ec¯ d2 + 1).

In this case, D = 2, and b = 2b¯ such that ΔSNRdB ≈ 3.0103Δb for QAM. This expression is similar to Eq. (5.17).

For square QAM, one has

Pe = 4 (1 − 1 M )Q ( 3 M − 1SNR) − 4 (1 − 1 M )2Q ( 3 M − 1SNR)2.
(5.21)

In this case, the gap approximation is based on the assumption that

Pe ≈ 4Q ( 3 M − 1SNR),
(5.22)

which allows to write

M = 1 + 3SNR [Q−1(Pe∕4)]2

and, alternatively

b = log2 (1 + SNR Γ ),

where

ΓQAM = 1 3 [Q−1 (Pe 4 )]2
(5.23)

is the gap for square QAM.

In summary, the actual supportable bit rate at a given error rate for uncoded QAM and PAM is given by the channel capacity for a modified SNR, which is the “capacity” SNR divided by the so-called SNR gap. The notion of a “gap” is more appropriate when the SNR is given in dB, since the SNR gap is then a fixed additive term for a given error rate:

10log10 (SNR Γ ) = SNRdBΓdB.

Note from Eq. (5.20) and Eq. (5.23) that the gap is independent of the number of bits b, (and so independent of M for M-QAM). It purely depends on the error rate.

Listing 5.8 was used to obtain Table 5.2.

Listing 5.8: MatlabOctaveCodeSnippets/snip_digi_comm_pam_qam_gaps.m
1Pe = [2e-5 1e-5 2e-6 1e-6 2e-7 1e-7 2e-8 1e-8 2e-9 1e-9]; 
2argQ_qam = qfuncinv(Pe/4); %QAM 
3gap_linear_qam = (argQ_qam.^2)/3; 
4gap_db_qam = 10*log10(gap_linear_qam); 
5argQ_pam = qfuncinv(Pe/2); %PAM 
6gap_linear_pam = (argQ_pam.^2)/3; 
7gap_db_pam = 10*log10(gap_linear_pam);

Table 5.2: Gaps in linear and dB scales for PAM and square QAM for a symbol error rate Pe.
QAM
PAM
Pe Γ ΓdB Γ ΓdB
2 ×10−5 6.5038 8.1317 6.0631 7.8269
10−5 6.9458 8.4172 6.5038 8.1317
2 ×10−6 7.976 9.0179 7.5317 8.7689
10−6 8.4213 9.2538 7.976 9.0179
2 ×10−7 9.458 9.758 9.011 9.5477
10−7 9.9056 9.9588 9.458 9.758
2 ×10−8 10.9471 10.393 10.4982 10.2111
10−8 11.3965 10.5677 10.9471 10.393
2 ×10−9 12.4416 10.9488 11.9912 10.7886
10−9 12.8924 11.1033 12.4416 10.9488

For example, assume a QAM that operates with SNR = 20 dB (100 in linear scale) and must have a symbol error rate Pe = 2 × 10−6. According to Table 5.2, in this case the gap is Γ = 7.976, which leads to

b = log2 (1 + SNR Γ ) = log2 (1 + 100 7.976 ) ≈ 3.6865.

Using Eq. (5.22), one obtains Pe ≈ 2 × 10−6 as desired. Recalling that the quantities per dimension for QAM are P¯e = Pe∕2 and b¯ = b∕2, the same project can be done evaluating a PAM for each dimension. In this case, P¯e = 10−6 and, for PAM, Γ = 7.976, such that:

b¯ = 0.5log2 (1 + SNR Γ ) = 0.5log2 (1 + 100 7.976 ) ≈ 1.8795

and, using Eq. (5.19), P¯e = 10−6. This corresponds to b = 2b¯ = 3.6865 bits and Pe = 2 × 10−6, as expected. Therefore, Eq. (5.18) is valid for both PAM and squared QAM.

Figure 5.19 was obtained with SNR = 10 dB for QAM. The abscissa informs the Pe that was specified and the “correct” graph corresponds to b using Eq. (5.21) while the other graph uses Eq. (5.23) to calculate Γ and then b^ using Eq. (5.18). It can be noticed that the number b^ of bits to achieve the specified Pe is underestimated by the gap approximation (b^ < b) .

PIC
Figure 5.19: QAM gap approximation for SNR = 10 dB: comparing the error with respect to the number of bits.

Figure 5.19 assumes that given SNR and Pe, the gap approximation is evaluated by comparing the number of bits.

PIC
Figure 5.20: QAM gap approximation for SNR = 5 dB: comparing the error with respect to the number of bits.

A similar procedure to Figure 5.19 was adopted to generate Figure 5.20 but, in this case, the SNR was decreased to SNR = 5 dB. The impact on the accuracy of the gap approximation can be noticed when the SNR decreases, especially for high Pe.