Revision as of 07:15, 10 September 2008 by Asabesan (Talk)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

This is the original code:


F0 = 13;
T0 = 1/F0;
Ts = 0.07;
t  = 0:Ts:13*T0;
x  = real(exp(j*(2*pi*F0*t-pi/2)));
plot(t,x)
<\pre>

This code is wrong because the sampling frequency, Ts, is to large to get an accurate recreation of the
signal. By Nyquist's theorem, the descrete sampling frequency must be twice the continuous frequency in order to avoid unwanted artifacts. In other words Ts=.5*T0.

<pre>
F0 = 13;
T0 = 1/F0;
Ts = .5*T0;
t  = 0:Ts:13*T0;
x  = real(exp(j*(2*pi*F0*t-pi/2)));
plot(t,x)
<\pre>
Explain what the bug is, and modify the above code to fix this bug. Post your answer on a Rhea page.

Alumni Liaison

Ph.D. on Applied Mathematics in Aug 2007. Involved on applications of image super-resolution to electron microscopy

Francisco Blanco-Silva