Pulse Shaping Filter

From the previous articles, we have seen how a Tx generates a signal from the numbers and the Rx recovers those numbers from that signal. In time domain, everything looks nice and perfect.

Let us investigate the system characteristics in frequency domain. As is clear from PAM block diagram, the main component that defines the spectral contents of the signal is the pulse shape at the Tx. We start with our attention towards the rectangular pulse shape used so far.

Frequency domain view

It is straightforward to find the DFT of many signals from DFT definition and the article on DFT examples, or a software routine. The discrete frequency plot as a function of k, however, can only show a limited amount of frequency content due to distinct frequency domain samples. In the spirit of simplicity we have adopted in this text, we will interpolate many of the future graphs for a finer resolution. All we have to do is choose a large value of when plotting.

As discussed in the article The Discrete Fourier Transform (DFT), for a reasonably band-limited signal, we can approximate the continuous frequency by discrete frequency . Then, does a larger imply that we can cover a larger frequency range in reality? No, recall from frequency relation that and

(1)

A large just makes the resolution finer and finer. The range stays within to and so does between and , until the sample rate itself changes. When we illustrate the frequency domain beyond this range, then we want to look into aliasing replicas of the spectrum as well.

Spectrum of a Rectangular Pulse

In time domain, a rectangular pulse has a clear benefit that Tx symbols do not naturally invade into each other’s territory. In case there are no reflections or dispersions in the channel, there are no adjacent symbols interfering into each other at the Rx as well — a phenomenon called Inter-Symbol Interference (ISI). ISI simplifies the detector in Rx design because everything related to a particular th symbol remains within the symbol interval . Or in other words, what happens in a symbol duration stays within that symbol duration.

As for the frequency domain, we derived an analytical expression for the DFT of a rectangular sequence in the article on DFT examples as

In the graph of this sinc function, nulls should occur wherever the numerator is zero, except at where the peak exists. Now when its argument is an integer multiple of . Thus, equating the argument of the numerator to returns the first null position.

Thus, the positions of nulls are with respect to discrete frequency . In Hz, with , the null positions are , , because

Finally, the range is

which for samples/symbol produces

The spectrum of a rectangular pulse shape shows these deductions in Figure below.

A great disadvantage of a rectangular pulse shape is now visible: the sidelobe level is just dB below the mainlobe, a very high value. Furthermore, the subsequent sidelobes are not decaying fast enough. If there is a communications link which occupies the drawn spectrum and the neighboring spectrum is allocated to someone else, then there will be a lot of interference between the two parties due to their sidelobe energies interfering with each other.

Just like real estate, radio spectrum is a very precious resource and must be utilized judiciously. This is why there are spectrum regulatory authorities in every country who impose strict restrictions to comply with a spectral mask. Even for wired channels, there is always a natural bandwidth of the medium (copper wire, coaxial cable, optical fiber) that imposes upper limits on its utilization.

Not only that a rectangular pulse shape is a poor choice due to its large spectral occupancy, another side effect of having a wider bandwidth pulse shape is that the matched filter at the Rx also has a similar bandwidth. Viewed in frequency domain, the larger the bandwidth of the Rx filter, the more noise and interference it allows to enter the system.

Randomness of bit stream

It is important to remember that such a spectrum approximately appears at the output of the transmitter only if there is a sufficient randomness in the bit stream, i.e., bits 0 and 1 occur with equal probability and so there are enough transitions between them. Any pattern in the bit stream can alter the output spectrum significantly.

As an example, imagine a constant bit stream of 1’s that is then modulated to a symbol level . After pulse scaling, the output is nothing but times an all-ones sequence that has a Fourier Transform of an impulse at bin (DC) and nothing else (see this Figure). In real applications, there is enough randomness in the modulated sequence and hence the spectrum will closely resemble the spectrum of a recangular pulse. We will see such an example a little later.

Reducing the bandwidth

Having known the disadvantages of a rectangular pulse shape, a superior pulse shaping filter needs to be looked for that should help in placing multiple channels adjacent to each other while minimizing inter-channel interference between them as well as noise bandwidth at the Rx. However, such advantages should not come at a price of introducing ISI in the system.

Nyquist No-ISI Criterion in Time Domain

To get that superior pulse shape, we trace the following sequence of steps through the help of Figure below that plots the auto-correlation of both a rectangular and a conceptual pulse shape. Auto-correlation of a pulse shape was defined earlier as

and plots are shown in continuous-time to make the related steps below more understandable.

• Remember that we discussed that signals that are narrow in time domain are wide in frequency domain, and vice versa. Accordingly, a consequence of narrow time span of the pulse shape is that it gets expanded in frequency. That is why a rectangular pulse has a large bandwidth. To decrease the bandwidth, the length of the pulse shape needs to be extended in time domain — much more than a single symbol duration .
• The effect of extending the length of the pulse shape is that every symbol now interferes with a number of symbols occurring both before and after that symbol. The dreaded ISI appears! Assuming a usual even shaped pulse, every symbol contributes towards ISI in symbols towards its left and symbols towards its right on time axis.
• To strike peace with this ISI, we exploit a key property of digital communication waveforms explained before: the critical samples at the output of matched filter occur at the end of integer multiples of symbol duration . For symbol detection, only these -spaced samples are needed and the rest can be discarded. This is what we saw at the output of downsampler in PAM block diagram.

If all but one sample during each symbol interval are discarded, they can be assigned any value without any effect on the detector output. These are the “don’t care” samples we give up for ISI to play with in the expectation that we get to shape our desired spectrum in return, and with the result that our desired -spaced samples remain intact. Keep in mind that we are talking about the output of the matched filter — and hence auto-correlation of the pulse shape — here. The reason is to elaborate the nature of signal that is eventually downsampled to yield symbol estimates. The actual pulse shape will be extracted from this auto-correlation, as we will shortly see.

• In the light of above discussion, three pulse shape auto-correlations are illustrated in Nyquist no-ISI Figure.
1. A rectangular pulse shape has a duration of seconds, or samples. Hence, its auto-correlation extending over seconds is shown as dotted red line in correlation of rectangular pulse as well. Notice that it has a maximum value at the desired current sampling instant and zero for adjacent two symbols. The time span, however, is too short giving rise to an unreasonably large spectrum.
2. A pulse shape auto-correlation that is expanded over many symbols is shown in blue curve, where it is scaled by a symbol . Note that at all the sampling instants (integer multiples of ), the effect of this particular symbol is zero. We intend to assign values to its “don’t care” samples in shaded regions to achieve a compact spectrum.
3. An adjacent symbol scales the same pulse shape at time . Observe the same zero-ISI effect on the current symbol at time and on all other symbols.

Not shown in the above figure are all other symbols shaped by pulses at , , and so on. But it is clear that sum total of ISI from all adjacent symbols is zero at sampling instant of every single symbol. For a concrete formulation of this zero-ISI criterion, we turn towards its mathematical foundation.

Just as in the case of rectangular pulse in the article on Pulse Amplitude Modulation (PAM), the input signal is

Although this equation was derived for a pulse with duration , it is also true for longer pulses because

This is because convolution between any signal and a unit impulse results in the same signal. Therefore, even for longer pulses, scaling each individual symbol with its own pulse shape is equivalent to convolution of symbol stream with a pulse shaping filter.

When this signal is input to matched filter , the output is written as

To generate symbol decisions, -spaced samples of the matched filter output at are

(2)

In the above equation, the first term is the currently desired symbol and the second term is inter-symbol interference (ISI). It is obvious that ISI can be zero if the pulse auto-correlation satisfies the condition

(3)

This is Nyquist no-ISI criterion in time domain — a mathematical expression of the same concept we explained above: to obtain zero ISI, the pulse auto-correlation should pass through zero for all integer multiples of before and after the current symbol.

Nyquist filter

A filter with impulse response coefficients satisfying Nyquist no-ISI criterion is called a Nyquist filter. In most texts, this is usually called a Nyquist pulse but we use the term Nyquist filter to avoid confusion, as this is actually not the pulse shape but the auto-correlation of the underlying pulse shape. The coefficients of the pulse shape itself will be derived later.

When Nyquist no-ISI criterion in time domain is satisfied, plugging Eq (3) into Eq (2) gives

i.e., the downsampled matched filter output maps back to the Tx symbol in the absence of noise. This is the crux of digital communications theory.

Nyquist no-ISI Criterion in Frequency Domain

Throughout the text, we have been avoiding complicated mathematical derivations. We will do the same here and instead of proving Nyquist no-ISI theorem which connects time and frequency domain properties of pulse auto-correlation, we again take the intuitive route.

As long as we operate at a sample rate , the spectral replicas are spaced at apart. But since we downsample by to get , our sample rate changes from at the output of matched filter to after the downsampler. To see what happens in frequency domain, observe that

• In discrete domain, in Eq (3) is nothing but a single impulse at time . This can be seen as the blue dots in first part of Figure below. Hence, its DFT is an all-ones rectangular sequence in frequency domain, which was derived in the article on DFT examples as

• Moreover, we know from the post on sample rate conversion that a consequence of downsampling by is appearance of new spectral aliases apart from each other.

Combining the above two facts produces second part of Figure below illustrating this relationship between symbol-spaced pulse auto-correlation and its Fourier transform.

We can also see that the signal bandwidth in this particular case is (never confuse a rectangular pulse shape in time domain with a pulse auto-correlation that is a rectangle in frequency domain)

This is the maximum bandwidth after which aliasing will occur from spectral replicas. This also sets the fundamental limit between the bandwidth and symbol rate as

(4)

Hence, with zero ISI, a maximum symbol rate equal to symbols/second can be supported within a bandwidth .

Remember that the point of this whole exercise was to limit the enormous bandwidth occupied by a time domain rectangular pulse. In other words, the quest is to design a pulse shape such that it complies with the spectral mask within a channel bandwidth . Nyquist no-ISI criterion in frequency domain helps here.

Nyquist showed that a necessary and sufficient condition for zero ISI condition of Eq (3) is that the Fourier transform of the pulse auto-correlation satisfies the condition

(5)

Its proof is not much complicated but we skip it anyway and rely on intuition. A hint about Nyquist no-ISI criterion in frequency domain comes from time and frequency views of Nyquist no-ISI: the spectrum of the pulse auto-correlation should be a constant. All it says is that for the DFT of pulse auto-correlation,

• Draw the primary spectrum whose range is .
• Draw its shifted replicas at frequencies . Repeat the same for all integer multiples of from to .
• Add all these replicas together.
• If this sum results in a constant, only then the time domain pulse auto-correlation will satisfy the no-ISI condition Eq (3).

Figure above shows two spectra, one of which satisfies the Nyquist no-ISI criterion while the other does not. The spectrum on the left is constant because the bandwidth is as wide as allowed by sampling theorem (which is for sample rate ). This not only avoids aliasing but owing to a flat spectrum within produces a flat spectrum in range , , and so on. Consequently, the cumulative spectrum as in Eq (5) is a constant function of frequency from to .

Interestingly, this maximum bandwidth allowed by sampling theorem is also the minimum bandwidth allowed by Nyquist no-ISI criterion in frequency domain. Any smaller than this, and it becomes impossible to obtain an overall flat spectrum. That is shown on the right of Figure above where there are spectral “holes” that cannot be filled regardless of the spectral shape.

Coefficients of an Ideal Pulse Auto-correlation

Above, we said that a filter whose impulse response coefficients satisfy Nyquist no-ISI criterion is a Nyquist filter, which is actually the auto-correlation of the underlying pulse shape. From Figure above and the discussion in the last section, we find our ideal Nyquist filter: a pulse auto-correlation whose spectrum is a rectangle from to . This filter is shown in Figure on time and frequency views of Nyquist no-ISI and left of Figure on minimum Nyquist bandwidth.

To derive the time domain representation of such a filter, remember that a rectangle in frequency domain is a sinc in time domain due to time frequency duality. This rectangle covers the entire continuous frequency range from to (or discrete frequency range to ) yielding in DFT of a rectangular signal (this is not samples/symbol, just the sequence length in that equation). There is a scaling factor of , furthermore, from the definition of iDFT. Combining these facts,

For our purpose with a sample rate of so far,

It satisfies Eq (3) because it is equal to at and zero for and shown as blue dots in Figure on time and frequency views of Nyquist no-ISI.

Using the identity for small , the term in the denominator is quite small as compared to numerator due to division by . Then, the above equation can be written as

This is the time domain waveform sampled at rate . As we are interested in time domain coefficients of this ideal pulse auto-correlation, we sample the waveform at rate instead of just , an increase by a factor of . Correspondingly, a decrease in sample time yields

(6)

These are the coefficients of our ideal pulse auto-correlation, which extends from to in time domain as shown in Figure below.

There are two main problems with this ideal pulse auto-correlation, both arising due to a very slow rate of decay for the tail.

[Truncation Error] Any practical system cannot use a filter with an infinite number of coefficients like the one shown in Figure of coefficients for an ideal pulse auto-correlation. Hence, the filter impulse response must be truncated at a few symbols to the left and a few symbols to the right of the current symbol. To ensure that this truncated filter closely approximates the ideal impulse response, only very small values can be ignored. That generates a large filter length .

[Timing Errors] We know that optimal timing instant coincides with symbol boundaries. Timing synchronization block in a Rx is responsible for extracting this symbol-aligned clock. However, even a small error in timing synchronization output causes a significant increase in ISI due to large value of samples at the tail interacting with neighbouring modulated pulses. The quicker the tail decay is, the lesser the errors arising from timing jitter when sampling adjacent pulses.

To address these issues, we take the reverse route now. We started with the aim of reducing the bandwidth of a rectangular pulse shape. We found an ideal pulse auto-correlation with a rectangular spectrum that satisfies Nyquist no-ISI criterion. To overcome the limitations in time domain associated with it, we seek to trade-off some bandwidth with desirable time domain behaviour.

Raised Cosine (RC) Filter

Remember that we observed regions of acceptable ISI in auto-correlation Figure of a conceptual pulse and deduced that these “don’t care” samples can be given up for ISI to play with in the expectation that we get to shape our desired spectrum in return, while our -spaced samples remain intact. If there is a relaxation and a restriction in time domain, then intuitively there should be a corresponding relaxation and restriction in frequency domain as well. We have found the restriction to be a flat spectrum across the whole frequency range. The relaxation is that how this spectrum becomes flat is a matter of no concern. The actual spectral shape is a “don’t care” case, as long as the sum of all shifted replicas is flat.

It is clear that we cannot reduce the bandwidth below the fundamental limit , see Figure on minimum Nyquist bandwidth. The only way to alter the spectral shape is by expanding the bandwidth beyond in such a way that the sum total of the replicas again becomes flat. For this purpose, an odd symmetry around the point is required. The advantage of such odd symmetry is that the spectral magnitudes before and after are rotated versions of each other. In other words, the spectral copies at fold around into the original bandwidth. That supplies the additional amplitude required to bring the spectrum in a flat shape, as illustrated in Figure below.

Higher frequency components of a signal arise from abrupt changes in time domain, which was the reason of large spectral occupancy of a rectangular pulse shape. Similarly, long tails in a time domain signal arise from abrupt transition in the flat spectrum as in Figure on time and frequency views of Nyquist no-ISI. Also recall from the post on FIR filters that the length of an FIR filter depends on its transition bandwidth.

Having located the abrupt transition bandwidth around as the root cause, we can extend the bandwidth of the pulse auto-correlation in any shape as long as it has odd symmetry around the points . The purpose is to smooth out its spectrum.

The smoothest spectral shape one can imagine is a sinusoid. If such a spectral taper is convolved with the ideal rectangular spectrum, the discontinuity in the spectrum can be removed. Since a half-cosine is an even symmetric shape, it is shown to be convolved with the rectangular spectrum in Figure below. The width of the half-cosine is where which forms the transition bandwidth of the resultant spectrum and its even symmetric shape preserves the odd symmetry around . As a consequence, this is also a Nyquist filter as it satisfies Nyquist no-ISI criterion and is called a Raised Cosine (RC) filter (it has nothing to do with an RC circuit consisting of a resistance and a capacitor). The name raised cosine comes from the fact that the transition bandwidth of the spectrum consists of a cosine raised by a constant to make it non-negative, see Eq (7) and Eq (8 later).

It can be observed that the effect of spectral convolution is an increase in bandwidth from to . Since specifies the bandwidth beyond the minimum bandwidth (Nyquist frequency or folding frequency), it is called excess bandwidth or roll-off factor. For example, when , the total bandwidth is more than the minimum and when , the total bandwidth is exactly more than, or twice, the minimum. Typical values of range from to .

In time domain, this is equivalent to product of a sinc signal with an even signal of time (since the spectrum is real and even, the time domain signal is also real and even). The sinc signal guarantees the zero crossings as they cannot be moved by multiplication, and the even signal dampens the long tails in time. The higher the bandwidth, the narrower this signal is and hence faster the decay.

First, consider , the widest bandwidth case. Here, the convolution between a rectangle of spectral width gets convolved with a cosine of width containing one full period. The convolution of the resultant spectrum has length and it is in cosine shape from to . This is plotted in Figure below, where smooth frequency transition can be seen with an odd symmetry around .

Let us explore its mathematical expression. Due to the way it is usually written in literature, the formula looks more complicated that it actually is. Recall that a sinusoid in time domain is written as

A frequency domain sinusoid will just interchange these roles, with inverse period and independent variable defined in terms of frequency. A half-cosine from to implies that the period in frequency of this cosine is . Also, it is raised by a constant to make it non-negative and scaled by to get a unity gain. Thus, the expression for this half-cosine spectrum is given as

The overall spectral shape from to can be given as

(7)

where the negative half of the spectrum is the same as the positive half, and hence the term .

We obtain its time domain expression through the following steps:

• In frequency domain, the constant term and cosine in the above expression do not have a frequency support from to , but only from to .
• In frequency domain, that is equivalent to multiplication with a rectangular window of the same width.
• Multiplication in frequency domain is convolution in time domain.
• In time domain, the constant term translates to an impulse at time location , and results in two impulses with half that amplitude at time locations (inverse period) .
• In time domain, the rectangular window of width in frequency is a sinc signal with zero crossings at integer multiples of .
• Convolution of a signal with an impulse is the signal itself. When convolved with sinc arising from the window, this generates three sinc signals: one at time while two at time .

The above sequence of steps result in time domain signal of an extended bandwidth pulse auto-correlation shown in Figure below. The reason of choosing a frequency support is now clear. With in frequency, the sinc resulting from cosine adds out of phase at with the sinc from the constant term. The signs of their tails are opposite to each other and that is what brings down the levels of long tails of an ideal pulse auto-correlation. Furthermore, sincs from both constant term and cosine have zero crossings at integer multiples of , as marked through the ellipse in the figure. Nevertheless, a time shift of for the sinc from cosine causes the cumulative signal to pass through zero at integer multiples of . This is how this pulse auto-correlation fulfills the Nyquist no-ISI criterion in time domain. Finally, sampling the pulse auto-correlation at samples/symbol produces the discrete-time coefficients that can be used for pulse shaping in digital domain.

Compare this discrete-time signal with the ideal Nyquist filter in Figure of coefficients for an ideal pulse auto-correlation. Now we can truncate the filter to a small length without much loss in accuracy. In addition, a sampling time misalignment at the Rx will have a minor influence on the current symbol and not much. Remember that the cost of this relief is the wider bandwidth . Observing this trend, we can think of striking a trade-off between the filter length and bandwidth expansion. The excess bandwidth comes into play here.

For this purpose, the excess bandwidth can be reduced from to a lower value. By reducing the bandwidth, the following sequence occurs.

• When , the period in frequency of the half-cosine in Raised Cosine spectrum with excess bandwidth equal to 1 decreases, while still maintaining odd symmetry around .
• As shown in Figure below, a decrease in frequency period results in three distinct regions of the spectrum: a half-cosine in the positive part of spectrum

a half-cosine in the negative portion of the spectrum

and a constant term

Note that when , the bandwidth was which can be written as . A decrease in results in new bandwidth as drawn in Figure above. Also observe that case coincides with a rectangular spectrum of ideal Nyquist filter as in Figure on time and frequency views of Nyquist no-ISI and Figure on minimum Nyquist bandwidth. This was the minimum bandwidth allowed by Nyquist no-ISI criterion.

• A decrease in bandwidth results in multiple windowed cosines and a constant. To avoid complication, its detailed analysis is not necessary. We just note that as gets smaller and smaller, the sincs in time domain resulting from windowed cosine in frequency move away from each other, thus causing less tail suppression as compared to case. In conclusion, the excess bandwidth gives the system designer a trade-off between reduced bandwidth and time domain tail suppression.

Taking clue from Eq (7), the general expression for a Raised Cosine can be written as
\small

(8)

\normalsize

The above equations looks intimidating but it is not. There are only three minor differences from Eq (7).

1. The middle term shows a reduced period of half-cosine and a spectral shift of cosine center from to . This is the edge of the passband now and is called the passband frequency.
2. The first term is a constant occupying the bandwidth left behind by the half-cosine.
3. The final term shows no bandwidth occupied from onwards. This is the start of the stopband and is called the stopband frequency.

Now it is clear why an oversampling factor of samples/symbol – a sampling rate of – is used in digital communication systems. Applying sampling theorem to the stopband frequency , the sampling theorem sets the minimum sampling rate at .

To obtain its time domain expression, recall Figure on spectral convolution where we found that the resultant waveform in time domain is equivalent to product of a sinc signal with an even signal of time.

(9)

This is drawn in Figure below for different values of . Note the simultaneous zero crossings of all waveforms at integer multiples of . Also, plugging above produces the coefficients of a sinc signal, the same ideal Nyquist filter in Eq (6).

The above equation becomes indeterminate for and . It can be shown that (we skip the derivation and use a mathematical technique called L’H\^{o}pital’s rule)

Resulting Pulse Shape

Until now in this section, everything we have discussed so far was about pulse auto-correlation that satisfies Nqyuist no-ISI criteria and is also called a Nyquist filter. However, how to derive the actual pulse shape is still not known. This is what we intend to find out next.

We start with a simple question: where should a Raised Cosine filter be placed in the system? There are only three possible choices.

[Receiver] The most straightforward way to incorporate it is to have it at the Rx. However, then a significant disadvantage is that there would be no mechanism to control the spectral sidelobes at the Tx. The spectral shaping required to minimize out-of-band energy cannot be avoided.

[Transmitter] If the spectrum is fully shaped at the Tx side, then any additional filtering at the Rx will not be “matched” to the incoming signal causing ISI. In the imaginary case of no filtering at all at the Rx, the noise and energy from adjacent channels will enter in the Rx in addition to unavoidable in-band noise. This significantly reduces the SNR as well as causing other issues such as increasing the interference sensitivity and dynamic range (adjacent channel energy can be quite higher than the desired band). So the Rx filter must be as compact as possible around the Tx spectrum so that maximal amounts of noise and adjacent channel interference can be eliminated.

[Both] The solution then is to split the Raised Cosine spectrum into two parts, one at the Tx to control the spectrum and the other at the Rx to reject noise and adjacent channels but still have zero ISI cumulative response. Remember that the pulse auto-correlation can be given by convolution formula as

The pulse shape can be placed at the Tx and its matched filter at the Rx. Above, the response of two filters in cascade is the convolution of their impulse responses, which implies multiplication of their frequency responses.

From above equation, the frequency response of the actual pulse shape can be derived as

(10)

A significant advantage of this arrangement is that the matched filter simultaneously maximizes the Rx SNR.

In the case of Raised Cosine pulse auto-correlation, the underlying pulse shape is called Square-Root Raised Cosine (SR-RC) pulse or filter, where the square-root is in frequency domain. Referring to Eq (8), taking the square-root does not affect and values. Using the identity and taking the square-root on both sides adjusts the middle term as
\small

(11)

\normalsize

For various values of excess bandwidth , this is shown in Figure below. For , there is no difference between a Raised Cosine and a Square-Root Raised Cosine filter due to a rectangular spectrum. Also notice from Eq (11) that the transition band of a Square-Root Raised Cosine pulse is a quarter cycle of a cosine as compared to a half-cosine of a Raised Cosine filter. It has implications which we shortly see.

Without mathematical derivation, the time domain waveform is given by

(12)

Figure below plots these time domain waveforms of Square-Root Raised Cosine for different values of . Observe again that plugging produces a sinc signal, the same ideal Nyquist filter in Eq (6). Through the red ellipse in the figure, the zero crossings are seen at integer multiples of only for . For and , Square-Root Raised Cosine does not satisfy Nyquist no-ISI criterion. This is not surprising because it has to fulfill that criterion only after matched filtering with another Square-Root Raised Cosine at the Rx to form the cumulative Raised Cosine shape.

In a software routine, the above equation can be used to find the coefficients of the Square-Root Raised Cosine pulse shape used for Tx filtering. Just like in Raised Cosine case, the denominator in this equation becomes zero for and . Pulse coefficients at these values can be given as (again, L’H\^{o}pital’s rule and skipping exact derivation}

(13)

Being a band-limited signal, the Square-Root Raised Cosine has is infinitely long in time that cannot be implemented in real systems. Therefore, it must be truncated to symbols to the left and symbols to the right. This is equivalent to in samples, thus resulting in total filter length . This time span of symbols or samples is called group delay (Group delay is a more general concept but this definition is good enough here).

Figure below compares the spectra of a Square-Root Raised Cosine pulse with that of a rectangular pulse (rectangular in time, not frequency). The Square-Root Raised Cosine pulse is generated using , samples/symbol and for a total filter length of samples. The duration of rectangular pulse is obviously samples. A huge improvement in sidelobe suppression is fairly visible.

Example

As an example, Wideband Code-Division Multiple-Access (WCDMA) – the main technology behind 3rd-generation (3G) cellular systems – implements a Square-Root Raised Cosine pulse shape with excess bandwidth , which translates the signaling rate of MHz to a bandwidth of MHz (the factor 2 arises for the RF bandwidth – we will discuss that in a later post). Accounting for the guard-bands to minimize interference between neighboring channels, the signal bandwidth in WCDMA systems is 5 MHz.

In the meantime, if you found this article useful, you might want to subscribe to my email list below to receive new articles.

PAM Revisited

Let us now replace the Square-Root Raised Cosine filter above into our basic PAM system of PAM block diagram in place of a rectangular pulse. For a -PAM modulation system with symbols , the process unfolds as follows. Keep comparing the signals thus generated with those from rectangular pulse shape in PAM block diagram.

• A source generates Tx bits to be sent to a destination. These bits are input to a look-up table that maps bits to symbols according to a chosen modulation scheme. Next, the Tx symbols are upsampled by samples/symbol. An example created of bits is illustrated in Figure below.

• The pulse shaping block with samples/symbol, excess bandwidth and group delay symbols filters the waveform to output a smooth waveform which then is passed to the analog portion of the Tx. A DAC produces the continuous output illustrated as dashed red line in Figure below. The bit information behind the waveform is also shown.

• In the best case scenario, i.e., a channel that only adds AWGN to the Tx signal, the Rx signal is drawn as in Figure below.

• The Rx signal is sampled by an ADC to produce and then input to a matched filter which is the same pulse shaping filter as at the Tx. The matched filter output at samples/symbol is shown in Figure below. As long as the the combination of Tx and Rx filters (i.e., overall Raised Cosine filter) obeys Nyquist criterion, there is no ISI and one out of every samples at can be preserved to form a symbol estimate while the remaining samples are thrown away. If there was zero noise, these downsampled values directly map to the transmitted symbols without any error, from which the sequence of Rx bits can be found with the help of an LUT.

To decide which one sample to keep out of every samples is the job of timing synchronization subsystem.

• When the noise is present which is always the case, the minimum distance rule is employed to produce symbol estimates according to its shortest distance from a constellation point, as illustrated in Figure below. Depending on numerous factors in system design, a proportion of bits will eventually end up in error, i.e., they will be different than Tx bits. The ratio of number of Rx bits in error to the total number of bits is called Bit Error Rate (BER).

Remember that in FIR filters, we said that every FIR filter causes the output signal to be delayed by a number of samples given by

(14)

and as a result, samples of the output at the start — when the filter response is moving into the input sequence — can be discarded. Similarly, the last samples — when the filter is moving out of the input sequence — should be discarded as well.

In summary, when the source information bits are filtered by a square-root Nyquist filter, the sharp edges visible for a rectangular pulse are smoothed out by a considerable margin, as shown in Figure below for 20 2-PAM symbols and excess bandwidth . This smoothness in time domain actually limits the bandwidth. Also observe that in the absence of noise, the values at optimum sampling locations are not all the same and exhibit ISI in a square-root Nyquist case. This is because the zero crossings of a square-root Nyquist shape are not necessarily at integer multiples of symbol times. However, after complete Nyquist filtering and no noise, all values coincide with the same optimum symbol amplitudes shown as red asterisks in Figure below.

Drawbacks of Square-Root Raised Cosine Pulse

As discussed above, Square-Root Raised Cosine pulse is much better than a rectangular pulse in shaping the spectrum but it has two major drawbacks.

[Insufficient Sidelobe Attenuation] There is a limit to the sidelobe suppression that a Square-Root Raised Cosine pulse can achieve. The sidelobe levels are reasonably higher than realistic spectral mask requirements imposed by regulatory authorities of attenuating out-of-band energy up to 60 to 80 dB. Looking back, recall that the transition band of a Raised Cosine pulse is half cycle of a cosine. Therefore, the transition band of a Square-Root Raised Cosine is a quarter cycle of a cosine, see Eq (11) and spectrum of Square-Root Raised Cosine. Its abrupt termination at the stopband results in a discontinuity causing a relatively poor sidelobe response.

[Increase in ISI] Looking at spectral comparison of Square-Root Raised Cosine with a rectangular pulse, an important question arises at this stage: The spectrum of Square-Root Raised Cosine should be precisely beyond Hz but why is it not? This is because as a consequence of truncation, the pulse is no more absolutely band-limited within and assumes infinite support in frequency in the form of sidelobes. Remember that truncation means multiplication by a rectangular window. This multiplication between Square-Root Raised Cosine pulse and rectangular window in time domain is convolution between Square-Root Raised Cosine spectrum and a sinc signal in frequency domain (the sidelobes and in-band ripple are inherited from that oscillating sinc signal and are a function of excess bandwidth and the pulse extension in symbols).

As a result of this truncation in time domain and subsequent convolution in frequency domain, the half amplitude values are moved away from the odd symmetry points of half symbol rate, or as illustrated in Figure below (compare to the case in Figure on odd symmetry). As a result, the spectral replicas now do not add up to exactly a constant value. In other words, Nyquist no-ISI criterion is only approximately satisfied by the truncated pulse shape, thus giving rise to increased ISI. The maximum magnitude of for is called peak ISI, shown in Figure below to be equal to for excess bandwidth and group delay .

The reader is encouraged to plot the spectra for different values of and to appreciate their role in spectral shaping.

Square-Root Raised Cosine and Raised Cosine are good starting points for pulse shape design and capture the involved details fairly well. Moreover, they have closed-form mathematical expressions that are good for analytical purposes. Nevertheless, the quarter cycle cosine of the Square-Root Raised Cosine results in high in-band ripples and insufficient out-of-band attenuation levels. There are other pulse shape design procedures that produce a Nyquist filter with high sidelobe attenuation and preserve Nyquist no-ISI criteria as well.

A Frequency Domain Window Based Pulse

A spectral shape is a discrete-time signal and hence just a sequence of numbers. This sequence of numbers in time domain can be generated through iDFT of a carefully designed discrete-time frequency response. Recall that a Raised Cosine filter was generated in frequency domain through convolution of an ideal rectangular spectrum with a half-cosine taper, see Figure on spectral convolution. Here, we replace the half-cosine with an alternative taper that is an improved spectral window. The span of this spectral window is determined by the excess bandwidth and the iDFT length.

The criteria for this spectral window design are

1. Narrow width of the mainlobe to effect a smoother transition band in resulting pulse spectrum. The width of the transition band remains unchanged
2. Small sidelobe levels that get inherited in the pulse spectrum

One such candidate is a Kaiser window in frequency domain (this technique is devised by Fred Harris in his paper “An Alternate Design Technique for Square-Root Nyquist Shaping Filters”) which has a minimum mainlobe spectral width for specified sidelobe levels. For a Kaiser window of a particular length, the sidelobe height is controlled by a parameter . In Matlab, a Kaiser window of length with parameter can be generated through the command kaiser(N,beta). Figure below compares this taper with a cosine taper with the same transition bandwidth and highlights the difference of smoothness between the two.

In frequency domain, the convolution of this Kaiser window taper with an ideal rectangular spectrum generates the desired pulse shape. This frequency domain result is shown in Figure below for excess bandwidth where it is compared with a Square-Root Raised Cosine pulse with the same excess bandwidth. This Figure depicts both the improved Nyquist filter and a Raised Cosine prior to square-root operation. Notice that a narrower mainlobe width of the improved taper has produced a smoother transition band in the improved Nyquist filter as compared to a Raised Cosine. This smoothness gets inherited by square-root Nyquist filter and it does not exhibit a similar kind of discontinuity as a Square-Root Raised Cosine, both of which are shown here. This discontinuity gives rise to high levels of time response even after many symbol intervals. When such an impulse response is truncated, it generates high sidelobes as well as in-band ripple in the spectrum of a Square-Root Raised Cosine. This is evident where an order of magnitude improvement in sidelobe suppression is visible. For example, an attenuation of dB is obtained through this procedure as compared to dB through Square-Root Raised Cosine.

Not shown here is the fact that low levels of in-band ripple accompanies the low levels of sidelobes because they are always equal in a window based design (arising from convolution of the same spectral taper sliding through both stopband and passband). For reasons that are beyond the scope of this text, the low in-band ripple generates less ISI after matched filtering of the square-root Nyquist filter. Therefore, the additional advantage of sidelobe suppression comes with a benefit, not at a cost.

Finally, there are other pulse shape design procedures as well that employ iterative techniques to convert an initial lowpass filter to a Nyquist filter with high sidelobe attenuation while preserving Nyquist no-ISI criteria as well. For the purpose of this text, we continue using the Square-Root Raised Cosine pulse and Raised Cosine pulse auto-correlation to keep things simple and assume that a reader who goes on to implement the system will use the better alternatives discussed above.

If you found this article useful, you might want to subscribe to my email list below to receive new posts.