We are limited to signal sequences that are discrete in both time and frequency domains from 0 to N − 1, which makes things a little easier.. The Power Spectrum We have learned that a ti
Trang 1for it is Eq (6-16), where we use E (expectation) to mean the same as
“averaging”, assuming that many repetitions have been performed:
ρxy = E {[x(n) − E(x(n))][y(n) − E(y(n))]}√
= E {[x(n) − E(x(n))][y(n) − E(y(n))]}
σXσY After many repetitions and averaging of ρxy, the numerator is the
expected value of the cross-covariance of x(n) and y(n) [Eq (6-15)],
and the denominator is the square root of the product of the variances of
x(n) and y(n), or more simply, justσxσy V (x(n)) and V (y(n)) or (σx and
σy) must both be greater than 0.0 This equation can be simpliÞed as
ρxy = E[x(n)y(n)]√ − E[x(n)]E[y(n)]
= E[x(n)y(n)] − E[x(n)]E[y(n)]
σXσY
If x(n) and y(n) are independent then the numerator of Eq (6-17) is zero:
E[x(n)y(n)] = E[x(n)]E[y(n)] (6-18) and ρxy= 0.0 However, there are some cases, not to be explored here,
where x(n) and y(n) are not independent, yetρxy is nevertheless equal to zero So “uncorrelated” and “independent” do not always coincide Look-ing at Eq (6-18), we can guess that this might happen For further insight about the correlation coefÞcient, see [Meyer, 1970, Chap 7]
As an example we will calculate ρxy of Fig 6-5 using Eq (6-17) and the time-averaged values instead of expected values because Eq (6-17)
is assumed to be noise-free:
ρXY = (xy) − xy
= 0.096 0.344 − (0.31 · 0.277)
· 0.286 = 0.099
The same calculation on Fig 6-4 produces a value of 1.00
Trang 2112 DISCRETE-SIGNAL ANALYSIS AND DESIGN
This brief introduction to correlation and variance is no more than a
“get acquainted” starting point for these topics and is not intended as a substitute for more advanced study and experience with probability and statistical methods We are limited to signal sequences that are discrete in
both time and frequency domains from 0 to N − 1, which makes things
a little easier Mathcad calculates very easily all of the equations in this chapter
REFERENCES
Carlson, A B., 1986, Communication Systems, 3rd ed., McGraw-Hill, New York Meyer, P L., 1970, Introductory Probability and Statistical Methods,
Addison-Wesley, Reading, MA
Oppenheim, A V., and R W Schafer, 1999, Discrete-Time Signal Processing,
2nd ed., Prentice Hall, Upper Saddle River, NJ
Schwartz, M., 1980, Information Transmission, Modulation and Noise, 3rd ed.,
McGraw-Hill, New York
Zwillinger, D., Ed., 1996, CRC Standard Mathematical Tables and Formulae,
CRC Press, Boca Raton, FL
Trang 3The Power Spectrum
We have learned that a time-domain discrete sequence x(n) that extends
from 0≤ n ≤ N − 1 can be considered as two-sided, positive-time for the Þrst half and negative-time for the second half Each sample x(n), considered by itself, is just a magnitude (see Chapter 1) It also has
a time-position attribute but none other, such as frequency or phase or
properties such as real or imaginary In other words, x(n) is not a phasor.
It is what we see on an ordinary oscilloscope
On the other hand, the x(n) sequence (the entire scope screen dis-play) can consist of a set of complex-valued voltage or current waveforms
applied to a complex load network of some kind However, time-domain analysis of complex signals combined with complex loads requires math methods that we will not explore in this book [Oppenheim and Schafer, 1999; Carlson, 1986; Schwartz, 1980; Dorf and Bishop, 2005; Shearer
et al., 1971], so we prefer to convert the time sequence x(n) to the fre-quency X (k) domain using the DFT After processing the signal in the
frequency domain we can, if we wish, use the IDFT to get the time
domain x(n) sequence representation of the processed discrete signal.
Discrete-Signal Analysis and Design, By William E Sabin
Copyright 2008 John Wiley & Sons, Inc.
113
Trang 4114 DISCRETE-SIGNAL ANALYSIS AND DESIGN
This is a simple and very useful approach that is widely used, especially
in computer-aided design
In this chapter we are interested in power We are also interested in phasors The problem is that any phasor that has constant amplitude has zero average power, so it makes no sense to talk about average phasor power Therefore, we will combine the positive- and negative-frequency phasors coherently, using the methods described in Fig 2-2 and employed elsewhere, to get a positive-frequency sine wave or cosine wave at
fre-quency (k) and phase θ(k) from 1 ≤ k ≤ N/2 − 1 We then have a true signal that has average power at frequency (k), and we can look at its
power spectrum
There is another approach available The real or imaginary part of the
phasor Me j ωt is a sinusoidal wave that has a peak value M The rms value
of this sinusoidal wave, considered by itself, is M· 0.7071 In our Mathcad examples the method of the previous paragraph, where we combine both sides of the phasor spectrum coherently, is an excellent and very simple approach that takes into account the two complex-conjugate phasors that are the constituents of the true sine or cosine signal
FINDING THE POWER SPECTRUM
We will use voltage values, but current values apply equally well, using the Norton source transformation [Shearer et al., 1971] The discrete Fourier
transform DFT [Eq (1-2)] of an x(n) signal sequence leads to a dis-crete two-sided X (k) steady-state spectrum of complex voltage phasors oscillating at frequency (k) from 1 to N − 1 with amplitude X (k) and
relative phase θ(k) At each (k) and (N − k) we will combine a pair of complex-conjugate phasors to get a positive-side sine or cosine V (k) from
k = 1 to k = N /2 − 1.
In order to keep the analysis consistent with circuit realities, assume that
V (k) is an open-circuit voltage generator whose internal impedance is for
now a constant resistance R g , but a complex Z g (k) can very easily be used.
V (k) is then the steady-state open-circuit voltage at frequency (k) The
voltage VL(k) across the load Y (k) (see Fig 7-1) is found independently
Trang 5V(k) ± jB(k)
G(k)
G B (k)
R G
Y(k)
VL(k) P(k)
PL(k):= (XE(k))2+ (XO(k))2
(1+ R·Y(k))2
·Y(k)
VL(k):= V(k)
1+ R G ·Y(k)
Figure 7-1 Equivalent circuit of frequency-domain power spectrum at
frequency position k.
at each frequency (k) as
V L(k)= 1 V (k)
+ RG Y (k) = 1 V (k)
+ ZG (k)Y (k) (7-1)
For a given V (k) the steady-state value of VL(k) depends only on the present value of Y (k) and Z G (k) or R G, and not on previous or future
val-ues of (k) The complex load admittance Y (k) = [G(k) + GB (k)] ± jB(k)
siemens which we have pre-determined by calculation or measurement at
each frequency (k), is driven by the complex signal voltage V L(k)=
Re[V L(k)] ± jIm[V L(k)], and the complex power to the load is
PL(k) = V L(k)2Y (k) watts and vars (7-2)
Figure 7-1 shows the equivalent circuit for Eq (7-1) PL(k) is the
power spectrum with real part (watts), imaginary part (vars), and phase angle θ(k), that is delivered to the complex load admittance Y (k) =
G(k) + GB (k) ± jB(k) If B(k) is zero, the power Pl(k) is in phase with