b The plot of the cosine wave with the noise just barely visible.. This mul-tiplication produces the baseband phase noise output VT n and a sine wave of amplitude 0.5 at twice the freque
Trang 1V1
V2
V1
V2
n N
n N
V1
V2
n N
2 )
1 N
k N
126
Trang 2At the upper left is a noise-free discrete sine wave V1(n) at frequency
f = 1.0, amplitude 1.0, in 128 positions of (n) A discrete cosine wave
V2(n), amplitude 1.0 at the same frequency, has some phase noise added,
0.1·[rnd(1) − 0.5] The rnd(1) function creates a random number from
0 to + 1 at each position of (n) The value 0.5 is subtracted, so the
random number is then between−0.5 and +0.5 The index of the phase
modulation is 0.1
(a) The plot of the noise-free sine wave
(b) The plot of the cosine wave with the noise just barely visible (c) We now multiply the sine wave and the cosine wave This
mul-tiplication produces the baseband phase noise output VT (n) and a
sine wave of amplitude 0.5 at twice the frequency of the two input waves We subtract this unwanted wave so that only the phase noise
is visible in part (c) This is equivalent to a lowpass Þlter that rejects the times 2 frequency Note the vertical scale in the graph of part (c) that shows the phase noise greatly ampliÞed
(d) We next use the DFT to get the noise spectrum VT (k) in dB format.
At this point we also perform two 3-point smoothing operations on
VT (k), Þrst to get VT1(k) and then to get VT2(k) This operation smoothes the spectrum of VT (k) so that VT (k) in the graph in part (e) is smoothed to VT2(k) in the graph in part (f) This is
postdetection Þltering that is used in spectrum analyzers and many
other applications to get a smoother appearance and reduce noise peaks; it improves “readability” of the noise shelf value
(e) Also in part (d) we perform lowpass Þltering [−20 log(1 + k2)]
(a Butterworth lowpass Þlter) to get VN (k) This result is also smoothed two times and the comparison of VN (k) and VN2(k) is
seen in the graphs of parts (e) and (f)
(f) Note that in parts (e) and (f) the upper level of the phase noise plot
VT (k) and VT2(k) is >53 dB below the 0-dB reference level of the test signal V2 at frequency k = 1 This is called the relative noise
shelf for the noisy test signal that we used It is usually expressed
as a dBc number (dB below the carrier, in this case, >53 dBc) This
noise shelf is of great signiÞcance in equipment design It deÞnes the ability to reject interference to and from closely adjacent signals
Trang 3128 DISCRETE-SIGNAL ANALYSIS AND DESIGN
and also to analyze unwanted phase disturbances on input signals
The lowpass Þlter to get VN2(k) greatly improves phase noise, but
only at frequencies somewhat removed from the signal frequency Still, it is very important that wideband phase noise interference is greatly reduced
In conclusion, there are many advanced applications of the cross power spectrum that we cannot cover in this book but that can be explored using various search engines and texts
REFERENCES
Carlson, A., 1986, Communication Systems, 3rd ed., McGraw-Hill, New York Dorf, R C., and R H Bishop, 2005, Modern Control Systems, 10th ed., Prentice
Hall, Englewood Cliffs, NJ
Gonzalez, G., 1997, Microwave Transistor AmpliÞer Analysis and Design,
Pren-tice Hall, Upper Saddle River, NJ
Oppenheim, A V., and R W Schafer, 1999, Discrete-Time Signal Processing,
2nd ed., Prentice-Hall, Upper Saddle River, NJ, p 189
Papoulis, A., 1965, Probability, Random Variables, and Stochastic Processes,
McGraw-Hill, New York
Sabin, W E., 1988, Envelope detection and noise Þgure measurement, RF
Design, Nov., p 29.
Sabin, W E., and E O Schoenike, 1998, HF Radio Systems and Circuits,
SciTech, Mendham, NJ
Schwartz, M., 1980, Information Transmission, Modulation and Noise, 3rd ed.,
McGraw-Hill, New York, Chap 5
Shearer, J L., A T Murphy, and H H Richardson, 1971, Introduction to System
Dynamics, Addison-Wesley, Reading, MA.
Trang 4The Hilbert
This Þnal chapter considers a valuable resource, the Hilbert transform (HT), which is used in signal-processing systems to achieve certain prop-erties in the time domain and the frequency domain The DFT, IDFT, FFT, and Hilbert transform work quite well together with discrete signals
if certain problem areas to be discussed are handled correctly
Example of the Hilbert Transform
Figure 8-1 shows a two-sided square wave x(n) time sequence, and
we will walk through the creation of an HT for this wave Design the
two-sided square-wave time sequence using N = 128 The value at n = 0
is zero, which provides a sloping leading edge for better plot results
Val-ues from 1 to N /2 − 1 = + 1.0 Set the N /2 position to zero, which has
been found to be important for successful execution of the HT because
N /2 is a special location that can cause problems because of its small
but non-zero value Values from N /2 + 1 to N − 1 = − 1 At N the wave
returns to zero
(a) Plot the two-sided square wave from n = 0 to N − 1.
Discrete-Signal Analysis and Design, By William E Sabin
Copyright 2008 John Wiley & Sons, Inc.
129
Trang 5130 DISCRETE-SIGNAL ANALYSIS AND DESIGN
(b) Execute the DFT to get X (k), the two-sided, positive-frequency
Þrst-half and negative-frequency second-Þrst-half phasor spectrum This spec-trum is a set of sine waves as shown in Fig 2-2c Each sine wave consists of two imaginary components
(c) Multiply the Þrst-half positive-frequency X (k) values by –j to get
a 90◦ phase lag Multiply the second-half negative-frequency X (k)
values by+ j to get a 90◦phase lead As in step (a), be sure to set the
N /2 value to zero This value can sometimes confuse the computer
unless it is forced to zero This step (c) is the ideal Hilbert transform
(d) Part (d) of the graph shows the two-sided spectrum XH (k) after the
phase shifts of part (c) This spectrum is a set of negative cosine waves as shown in Fig 2-2b Each cosine wave consists of two real negative components, one at +k and one at N − k.
(e) Use the IDFT to get the two-sided xh(n) time response Note that
xh(n) is a sequence of real numbers because x(n) is real.
(f) The original square wave and its HT are both shown in the graph Note the large peaks at the ends and in the center These peaks are
characteristic of the HT of an almost square wave The perfect square
wave would have inÞnite (very undesirable) peaks
(g) Calculate two 3-point smoothing sequences described in Chapter 4 for the sequence in part (f) This smoothing is equivalent to a lowpass Þlter, or at radio frequencies to a narrow band-pass Þlter
(h) Plot the Þnal result The sharp peaks have been reduced by about
2 dB Further smoothing is usually required in narrowband circuit design applications, as we will see later
The three peaks in part (f) are usually a problem in any peak-power-limited system (which is almost always the practical situation) The smoothing in part (h) thus becomes important Despite the peaks, the rms voltage in part (f) is the same for both of the waveforms in that diagram (nothing is lost)
For problems of this type, the calculation effort becomes extensive, and the use of the FFT algorithm, with its greater speed, would ordinarily
be preferred The methods of the Mathcad FFT and IFFT functions are
described in the User Guide and especially in the online Help In this