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)
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.
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)
Explain what the bug is, and modify the above code to fix this bug. Post your answer on a Rhea page.