Lab 5 Primer#

Slides

The range of complexity for a pulse amplitude modulation (PAM) transceiver system can vary drastically. A working version that demonstrates the core principles can be simulated with just a few lines of code. In practice, many additional subsystems need to be added for reliable operation.

Simulation#

In simulation, pulse amplitude modulation is relatively straightforward, and can be achieved with just a handful steps.

Mapping data to pulses#

The first step in any digital communication system is serialization, the process of converting a useful data structure into a stream of bits. For example, the following MATLAB code serializes the ascii string “meatball” into an array of 56 bits.

>> serialized_data = reshape((dec2bin(uint8('meatball'))'),[1,56])

serialized_data =

    '11011011100101110000111101001100010110000111011001101100'

In 2-PAM, we map a bit of ‘0’ to a pulse with an amplitude \(-d\) and a bit of one to a pulse with an amplitude of \(+d\). A common pulse shape is the raised cosine

upsampled_data = upsample(str2num(serialized_data')*2-1,16);
pulse_shaping_filter = rcosdesign(0.8,4,16,'normal');
sequence_of_pulses = conv(upsampled_data,pulse_shaping_filter,'same')
figure; plot(sequence_of_pulses); hold on; stem(upsampled_data.*max(pulse_shaping_filter));

Polyphase filter bank for pulse shaping#

Recall that we upsampled the data by a factor of \(L\) before applying the pulse shaping filter. Using the standard FIR filter structure, this would mean that a large fraction, \(\frac{L-1}{L}\), of the multiply-accumulate operations used in the filter are wasted on zeros! The polyphase filterbank structure solves this problem.

The diagram above depicts the polyphase structure to implement an 8-tap FIR filter with impulse response \(g[n]\) after upsampling by a factor of \(L=4\). The filter coefficients \(g[n]\) are divided into \(L\) smaller filters of length \(\frac{8}{L}=2\). Even though the upsampling operation quadrupled the sampling rate of the output, these smaller filters are allowed to run at the original rate, reducing the computational cost by a factor of roughly \(L=4\). By allowing the commutator the cycle through these four smaller filters, we produce the same output signal as if we had upsampled first, then applied an FIR filter at a higher rate.

Modulation#

The signal above is what we call the baseband signal. Before transmitting, we want shift the spectrum to higher frequencies in a way that can be reversed at the transmitter. A simple method is with sinusoidal modulation and demodulation

f0 = 12000; fs = 48000; w0 = 2*pi*f0/fs; n = 1:length(sequence_of_pulses);
modulated = sequence_of_pulses.' .* cos(w0*n);
figure; plot(sequence_of_pulses); hold on; plot(modulated);

Demodulation#

Demodulation can be achieved by once again multiplying by the carrier. The second multiplication by cosine introduces frequency components at \(\pm f_c\), which can be removed by a low pass filter.

demodulated = modulated .* cos(w0*n);
recovered_pulses = conv(demodulated,pulse_shaping_filter,'same');
figure; plot(sequence_of_pulses); hold on; plot(recovered_pulses);

Mapping pulses to data#

We end up with a signal that can be downsampled and quantized to recover the original data.

>> recovered_data = downsample(recovered_pulses,16) > 0

recovered_data =

    '11011011100101110000111101001100010110000111011001101100'

Receiver#

The receiver is where most of complexities occur. In the simulation above, we assumed that the receiver’s carrier exactly matches the transmitter (both in frequency and phase). We assumed that the signal traveled from transmitter to receiver with no distortion or time delay. We also assumed that the transmitter and receiver have exactly the same sampling rates so that every symbol was separated by exactly 16 samples. Though these assumptions are rarely true in practice, each type of error can be addressed by adding an appropriate subsystem.

Costas loop for carrier recovery#

The Costas loop is one way to track the phase \(\theta\) using an iterative algorithm. When processing the \(k\)th sample at the receiver, the Costas loop produces the next phase estimate \(\theta[k+1]\) based on the following rule:

\[ \theta[k+1] = \theta[k] - \mu \text{LPF} \left\{ r[k] \cos(\omega_0 k + \theta[k]) \right\} \times \text{LPF} \left\{ r[k] \sin(\omega_0 k + \theta[k]) \right\}\]

where \(r[k]\) is the received signal and \(w_0\) is the carrier frequency.

For an in-depth discussion of the Costas loop, see Software Receiver Design by C. Richard Johnson, Jr., William A. Sethares and Andrew Klein, section 10.4.

Symbol Clock Recovery#

The carrier recovery process will correct a frequency offset so that the signal can be properly demodulated.

The demodulated signal is downsampled and quantized to produce a stream of bits at the receiver. This process is sensitive to any timing offset. For example, the red and black stem plots show two possible ways that the demodulated signal might be downsampled. In one case (red), the values are recovered exactly. In the other case (black) many of the samples are close to zero and might easily be overwhelmed by noise.

A symbol clock recovery subsystem should adapt the sampling times to minimize errors.

Symbol Clock Recovery#

One simple clock recovery method consists of two bandpass filters and a squaring block.

The first bandpass filter passes any oscillations that occur at a rate \(f_{sym}/2\), which is the highest frequency that can occur in the baseband signal. Squaring the result will double the frequency. The second bandpass filter will then produce an estimate of the symbol clock.

Error probability#

It is common to model each sample of noise an independent draw from Gaussian distribution with variance \(\sigma^2\).

As the noise variance increases, it is more likely that a symbol is incorrectly decoded. However, as the symbol spacing \(A\) increases, a larger amount of noise is required to cause a bit error.

As a result, the bit error probability will depend on the ratio \(A/\sigma\). Since the signal power is proportional to \(A^2\) and the noise power is \(\sigma^2\), this can be expressed as

\[\frac A \sigma = \sqrt{\frac{A^2}{\sigma ^2}} = \sqrt{\text{SNR}}\]

Power maximization for timing recovery#

For an in-depth discussion of timing recovery via output power maximization, see Software Receiver Design by C. Richard Johnson, Jr., William A. Sethares and Andrew Klein, section 12.4.