Homework Seven: Difference between revisions

From Class Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(3 intermediate revisions by the same user not shown)
Line 22: Line 22:
[[Image:signal_sampled.jpg |thumb|center|upright=2|Figure 4: Signal being sampled at a rate of 1/T.]]
[[Image:signal_sampled.jpg |thumb|center|upright=2|Figure 4: Signal being sampled at a rate of 1/T.]]


In order to get the original signal, we simply need to create a lowpass filter that will essentially encompass the the original (desired) signal:
In order to get the original signal, we simply need to create a lowpass filter that will essentially encompass the the original (desired) signal and filter out any high frequency components:


[[Image:signal_lowpass.jpg |thumb|center|upright=2|Figure 5: The use of a lowpass filter to obtain the original  signal as shown in Figure 1.]]
[[Image:signal_lowpass.jpg |thumb|center|upright=2|Figure 5: The use of a lowpass filter to obtain the original  signal as shown in Figure 1.]]


So what have we accomplished? We have taken a signal (Figure 1), prepared it to be sampled (Figure 2 & 3), sampled it at a sampling rate of <math>\textstyle f_s = \frac{1}{T} </math> (Figure 4), and used a lowpass filter to collect the original (now sampled) signal (Figure 5).
[[Image:original.jpg |thumb|center|upright=2|Figure 6: Original signal obtained using a lowpass filter.]]
 
So what have we accomplished? We have taken a signal (Figure 1), prepared it to be sampled (Figure 2 & 3), sampled it at a sampling rate of <math>\textstyle f_s = \frac{1}{T} </math> (Figure 4), and used a lowpass filter to collect the original (now sampled) signal (Figure 5 & 6).
<br/>
<br/>
 
Our next task is to find an expression of x(t) in terms of x(nT). To accomplish this, recall from class that the transfer function of the lowpass filter is
<br/>
<br/>
 
<math>H(f)=\begin{cases} T,  & \mbox{ } |f| < \frac{1}{2T} \\ 0, & \mbox{ else } \end{cases}\!</math><br><br>
 
<br/>
 
So to find the transfer function in the time-domain we simply need to find the inverse Fourier Transform of H(f):
 
<math> h(t) = \mathcal{F}^{-1}[H(f)] = \int_{-\frac{1}{2T}}^{\frac{1}{2T}}Te^{j2 \pi ft} \,df  = \frac{Te^{j2 \pi ft}}{j2 \pi t} \Bigg|_{-\frac{1}{2T}}^{\frac{1}{2T}} =  \frac{T}{\pi t}\left[ \frac{e^{j \pi t/T}-e^{-j \pi t/T}}{j2}\right] </math>
 
By Euler's Identiy, however, we can see that
<br/>
<br/>
<math> \frac{T}{\pi t}\left[ \frac{e^{j \pi t/T}-e^{-j \pi t/T}}{j2}\right] = \frac{T}{\pi t} sin \left( \frac{ \pi t}{T} \right) </math>
<br/>
<br/>
 
Recall that <math>\textstyle sinc(\theta ) = \frac{sin(\pi \theta )}{\pi \theta }</math>. Thus,
<br/>
<br/>
<math> \frac{sin \left( \frac{ \pi t}{T} \right)}{\frac{\pi t}{T}} = sinc\left( \frac{t}{T} \right) \Longrightarrow h(t) = sinc\left( \frac{t}{T} \right)</math>
<br/>
<br/>
 
We learned in class that <math> \textstyle x'(t) = \sum_{n=-\infty }^{\infty}x(nT)\delta (t-nT)</math> and, thus, we must convolve <math>\textstyle x'(t) \mbox{ with } h(t)</math> to find the function of <math>\textstyle x(t) </math> desired:
<br/>
<br/>
<math> x(t) = x'(t)*h(t) = \sum_{n=-\infty }^{\infty}x(nT)\delta (t-nT)*\left[ sinc\left( \frac{t}{T} \right) \right]
= \sum_{n=-\infty }^{\infty}x(nT) \mbox{ } sinc\left( \frac{t-nT}{T} \right)
</math>
<br/>
<br/>
Therefore, the original signal (in terms of <math> \scriptstyle x(nT)</math>) is
<br/>
<br/>
<math> x(t)= \sum_{n=-\infty }^{\infty}x(nT) \mbox{ } sinc\left( \frac{t-nT}{T} \right) = \sum_{n=-\infty }^{\infty}\frac{x(nT) \mbox{ } sin \left( \frac{\pi (t-nT)}{T} \right)}{\frac{\pi (t-nT)}{T}}</math>

Latest revision as of 18:17, 29 November 2009

Figure out what happens if your sampled signal, x(t), has frequency components only for 12fs<f<fs. Can you recover the original signal from it? If so, find the expression for x(t) in terms of x(nt).


Nick Christman


Imagine we have an original signal that has the following frequency plot, centered about the "zero frequency" origin:

Error creating thumbnail: File missing
Figure 1: Original signal for which we want to analyze.

In order to analyze this signal we must look at both the positive and negative frequency aspects -- therefore, we will split the signal into two parts, positive and negative:

Error creating thumbnail: File missing
Figure 2: Signal separated into positive and negative aspects.

To "send the signal" we must essentially move it to a higher frequency and, remembering that x(t) has frequency components for 12fs<f<fs, we get the following signal, X(f):

Error creating thumbnail: File missing
Figure 3: Preparing signal to be sent and sampled.

Finally, we need to sample the signal at a rate of fs=1T. Theoretically, this leads to the following frequency plot:

Error creating thumbnail: File missing
Figure 4: Signal being sampled at a rate of 1/T.

In order to get the original signal, we simply need to create a lowpass filter that will essentially encompass the the original (desired) signal and filter out any high frequency components:

Error creating thumbnail: File missing
Figure 5: The use of a lowpass filter to obtain the original signal as shown in Figure 1.
Error creating thumbnail: File missing
Figure 6: Original signal obtained using a lowpass filter.

So what have we accomplished? We have taken a signal (Figure 1), prepared it to be sampled (Figure 2 & 3), sampled it at a sampling rate of fs=1T (Figure 4), and used a lowpass filter to collect the original (now sampled) signal (Figure 5 & 6).

Our next task is to find an expression of x(t) in terms of x(nT). To accomplish this, recall from class that the transfer function of the lowpass filter is

H(f)={T, |f|<12T0, else 


So to find the transfer function in the time-domain we simply need to find the inverse Fourier Transform of H(f):

h(t)=1[H(f)]=12T12TTej2πftdf=Tej2πftj2πt|12T12T=Tπt[ejπt/Tejπt/Tj2]

By Euler's Identiy, however, we can see that

Tπt[ejπt/Tejπt/Tj2]=Tπtsin(πtT)

Recall that sinc(θ)=sin(πθ)πθ. Thus,

sin(πtT)πtT=sinc(tT)h(t)=sinc(tT)

We learned in class that x(t)=n=x(nT)δ(tnT) and, thus, we must convolve x(t) with h(t) to find the function of x(t) desired:

x(t)=x(t)*h(t)=n=x(nT)δ(tnT)*[sinc(tT)]=n=x(nT) sinc(tnTT)

Therefore, the original signal (in terms of x(nT)) is

x(t)=n=x(nT) sinc(tnTT)=n=x(nT) sin(π(tnT)T)π(tnT)T