1. close all
  2. clear all
  3. clc
  4.  
  5. x = inline('upulse(t-0.5,0,0.5,0) - upulse(t,0,0,0.5)');
  6. t = linspace(0,1,1001);
  7. xat = -.55*sin(2*pi*t)-0.13*sin(4*pi*t);
  8. fs = 5; T = 1/fs;
  9. tn = 0.2:T:1; xn = inline('upulse(tn-0.5,0,0.5,0) - upulse(tn,0,0,0.5)');
  10. subplot(121), plot(t, xat,t, x(t), tn, xn(tn), 'o')
  11. grid on; xlabel('t (sec)');ylabel('x(t),xa(t)'); title('fs = 5 Hz')
  12.  
  13. % part e
  14.  
  15. xat10 = -.6155*sin(2*pi*t)-0.2753*sin(4*pi*t)-.1453*sin(6*pi*t)-.065*sin(8*pi*t);
  16. fs10 = 10; T10 = 1/fs10;
  17. tn10 = 0.1:T10:1; xn10 = inline('upulse(tn-0.5,0,0.5,0) - upulse(tn,0,0,0.5)');
  18. subplot(122), plot(t, xat10,t, x(t), tn10, xn10(tn10), 'o')
  19. grid on; xlabel('t (sec)');ylabel('x(t),xa(t)'); title('fs = 10 Hz')
  20.  
  21. a = [sin(pi/5) sin(2*pi/5) sin(3*pi/5) sin(4*pi/5);
  22.     sin(2*pi/5) sin(4*pi/5) sin(6*pi/5) sin(8*pi/5);
  23.     sin(3*pi/5) sin(6*pi/5) sin(9*pi/5) sin(12*pi/5);
  24.     sin(4*pi/5) sin(8*pi/5) sin(12*pi/5) sin(16*pi/5)];