02.10.2019 Views

UploadFile_6417

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Sampling and Reconstruction of Analog Signals 85<br />

which is a real-valued function since x a(t) isareal and even signal. To evaluate<br />

X a(jΩ) numerically, we have to first approximate x a(t) byafinite-duration<br />

grid sequence x G(m). Using the approximation e −5 ≈ 0, we note that x a(t)<br />

can be approximated by a finite-duration signal over −0.005 ≤ t ≤ 0.005 (or<br />

equivalently, over [−5, 5] msec). Similarly from (3.36), X a(jΩ) ≈ 0 for Ω ≥<br />

2π (2000). Hence choosing<br />

∆t =5× 10 −5 ≪ 1 =25× 10−5<br />

2(2000)<br />

we can obtain x G(m) and then implement (3.35) in MATLAB.<br />

% Analog Signal<br />

>> Dt = 0.00005; t = -0.005:Dt:0.005; xa = exp(-1000*abs(t));<br />

% Continuous-time Fourier Transform<br />

>>Wmax = 2*pi*2000; K = 500; k = 0:1:K; W = k*Wmax/K;<br />

>>Xa = xa * exp(-j*t’*W) * Dt; Xa = real(Xa);<br />

>>W = [-fliplr(W), W(2:501)]; % Omega from -Wmax to Wmax<br />

>>Xa = [fliplr(Xa), Xa(2:501)]; % Xa over -Wmax to Wmax interval<br />

>>subplot(2,1,1);plot(t*1000,xa);<br />

>>xlabel(’t in msec.’); ylabel(’xa(t)’)<br />

>>title(’Analog Signal’)<br />

>>subplot(2,1,2);plot(W/(2*pi*1000),Xa*1000);<br />

>>xlabel(’Frequency in KHz’); ylabel(’Xa(jW)*1000’)<br />

>>title(’Continuous-time Fourier Transform’)<br />

Figure 3.11 shows the plots of x a(t) and X a(jΩ). Note that to reduce the number<br />

of computations, we computed X a(jΩ) over [0, 4000π] rad/sec (or equivalently,<br />

over [0, 2] KHz) and then duplicated it over [−4000π, 0] for plotting purposes.<br />

The displayed plot of X a(jΩ) agrees with (3.36).<br />

□<br />

□ EXAMPLE 3.19 To study the effect of sampling on the frequency-domain quantities, we will<br />

sample x a(t) inExample 3.18 at 2 different sampling frequencies.<br />

a. Sample x a(t) atF s = 5000 sam/sec to obtain x 1(n). Determine and plot<br />

X 1(e jω ).<br />

b. Sample x a(t) atF s = 1000 sam/sec to obtain x 2(n). Determine and plot<br />

X 2(e jω ).<br />

Solution<br />

a. Since the bandwidth of x a(t) is2KHz, the Nyquist rate is 4000 sam/sec,<br />

which is less than the given F s. Therefore aliasing will be (almost) nonexistent.<br />

Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).<br />

Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!