THE BUG
% Mark Frankosky % hw2.2 F0 =13; T0 =1/F0; Ts = 0.07; t = 0:Ts:13*T0*13; x = real(exp(j*2*pi*F0*t-pi/2)); % the bug was in the above line which read x = real(exp(j*(2*pi*F0*t-pi/2))); % which would have been e^(j*2*pi*F0*t)*e^(-pi/2) which wouldnt have been % the same as having a phase shift of pi/2 which is i guess what was % needed. plot(t,x)