<< Chapter < Page Chapter >> Page >
J S D ( θ ) = 1 4 LPF { ( r p ( k T s ) - cos ( 4 π f 0 k T s + 2 θ ) ) 2 } .

Substituting r p ( k T s ) from [link] and assuming f c = f 0 , this can be rewritten

= 1 4 LPF { ( cos ( 4 π f 0 k T s + 2 Φ ) - cos ( 4 π f 0 k T s + 2 θ ) ) 2 } .

Expanding the square gives

= 1 4 LPF { cos 2 ( 4 π f 0 k T s + 2 Φ ) - 2 cos ( 4 π f 0 k T s + 2 Φ ) · cos ( 4 π f 0 k T s + 2 θ ) + cos 2 ( 4 π f 0 k T s + 2 θ ) } .

Using the trigonometric formula [link] for the square of a cosine and the formula [link] for the cosine angle sum (i.e., expand cos ( x + y ) with x = 4 π f 0 k T s and y = 2 Φ , and then again with y = 2 θ ) yields

J S D ( θ ) = 1 8 LPF { 2 + cos ( 8 π f 0 k T s + 4 Φ ) - 2 cos ( 2 Φ - 2 θ ) - 2 cos ( 8 π f 0 k T s + 2 Φ + 2 θ ) + cos ( 8 π f 0 k T s + 4 θ ) } .

By the linearity of the LPF,

= 1 4 + 1 8 LPF { cos ( 8 π f 0 k T s + 4 Φ ) } - 1 4 LPF { cos ( 2 Φ - 2 θ ) } - 1 4 LPF { cos ( 8 π f 0 k T s + 2 Φ + 2 θ ) } + 1 8 LPF { cos ( 8 π f 0 k T s + 4 θ ) } .

Assuming that the cutoff frequency of the lowpass filter is less than 4 f 0 , this simplifies to

J S D ( θ ) 1 4 ( 1 - cos ( 2 Φ - 2 θ ) ) ,

which is shown in the top plot of [link] for Φ = - 0 . 8 . The algorithm [link] is initialized with θ [ 0 ] at some point on the surface of the undulating sinusoidal curve. At each iterationof the algorithm, it moves downhill. Eventually, it will reach one of the nearby minima, which occur at θ = - 0 . 8 ± n π for some n . Thus, [link] provides evidence that the algorithm can successfully locate the unknown phase, assuming that the preprocessed signal r p ( t ) has the form of (10.3).

The error surface Equation 17 for the SD phase tracking algorithm is shown in the top plot. Analogous error surfaces for the phase locked loop Equation 22 and the Costas loop Equation 27 are shown in the middle and bottom plots. All have minima (or maxima) at the desired locations (in this case -0.8) plus nπ offsets.
The error surface [link] for the SD phase tracking algorithm is shown in the top plot. Analogous error surfaces forthe phase locked loop [link] and the Costas loop [link] are shown in the middle and bottom plots. All have minima (or maxima) at the desired locations (in thiscase - 0.8) plus n π offsets.

[link] shows the algorithm [link] with the averaging operation replaced by the more general LPF. In fact, this provides a concrete answer to [link] ; the averaging, the LPF, and the integral block allact as lowpass filters. All that was required of the filtering in order to arrive at [link] from [link] was that it remove frequencies above 4 f 0 . This mild requirement is accomplished even bythe integrator alone.

A block diagram of the squared difference phase tracking algorithm Equation 11. The input r_p(kT_s) is a preprocessed version of the received signal as shown in Figure 10-3. The integrator block Σ has a lowpass character, and is equivalent to a sum and delay as shown in Figure 7-5.
A block diagram of the squared difference phase tracking algorithm [link] . The input r p ( k T s ) is a preprocessed version of the received signal as shown in [link] . The integrator block Σ has a lowpass character, and is equivalent to a sum and delay as shown in [link] .

The code in pllsd.m is simplified in the sense that the received signal rp contains just the unmodulated carrier. Implement a more realistic scenario by combining pulrecsig.m to include a binary message sequence, pllpreprocess.m to create rp , and pllsd.m to recover the unknown phase offset of the carrier.

Using the default values in pulrecsig.m and pllpreprocess.m results in a ψ of zero. [link] provided several situations in which ψ 0 . Modify pllsd.m to allow for nonzero ψ , and verify the code for the cases suggested in [link] .

Investigate how the SD algorithm performs when the received signal contains pulse shaped 4-PAM data.Can you choose parameters so that θ Φ ?

Consider the sampled cosine wave

x ( k T s ) = cos ( 2 π f 0 k T s + α )

where the frequency f 0 is known but the phase α is not. Form

v ( k T s ) = x ( k T s ) cos ( 2 π f 0 k T s + β ( k ) )

using the current (i.e. at the k th sample instant) estimate β ( k ) of α and define the candidate cost function

J ( β ) = LPF { v 2 ( k T s ) }

where the cutoff frequency of the ideal lowpass filter is 0 . 8 f 0 .

  1. Does minimizing or maximizing J ( β ) result in β = α ? Justify your answer.
  2. Develop a small stepsize gradient descent algorithm for updating β ( k ) to estimate α . Be certain that all of the signals needed to implement thisalgorithm are measurable quantities. For example, x is directly measurable, but α is not.
  3. Determine the range of the initial guesses for β ( k ) in the algorithm of part (b) that will lead to the desired convergence to α given a suitably small stepsize.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Software receiver design. OpenStax CNX. Aug 13, 2013 Download for free at http://cnx.org/content/col11510/1.3
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Software receiver design' conversation and receive update notifications?

Ask