Advertisement
squelch

EE110B Lab2 in progress

Jan 23rd, 2017
124
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.00 KB | None | 0 0
  1. % LAB 2
  2.  
  3. % TASK 1
  4.  
  5. % Create a random sequence.
  6. x = randi(10,1,100);
  7.  
  8. % Part (a)
  9. n = 1:101;
  10. h1 = 0.9.^(n-1);
  11. h2 = [0,ones(1,99),0];
  12. h = h1.*h2;
  13. y1a = conv(h,x);
  14. stem(y1a)
  15.  
  16. % Part (b)
  17. y1b = zeros(1,101);
  18. y1b(1)=0;
  19. for i = 1:100
  20. y1b(i+1)=y1b(i)+x(i);
  21. end
  22. stem(y1b)
  23.  
  24. % Part (c)
  25. % No code needed. Discuss results.
  26.  
  27. % TASK 2
  28.  
  29. % Create a new random sequence.
  30. x = randi(10,1,101);
  31.  
  32. % Part (a)
  33. n = 1:101;
  34. h1 = (0.9.^(n-1)).*cos((pi/5)*(n-1));
  35. h2 = [0,ones(1,99),0];
  36. h = h1.*h2;
  37. y2a = conv(h,x);
  38. stem(y2a)
  39.  
  40. % Part (b)
  41. y2b = zeros(1,102);
  42. y2b(1)=0;
  43. y2b(2)=0;
  44. for i = 1:100
  45. y2b(i+2)=(1.8).*cos(pi/5).*y2b(i+1) - 0.81.*y(i) + x(i+1) - (0.9).*cos(pi/5).*x(i);
  46. end
  47. stem(y2b)
  48.  
  49. % TASK 3
  50.  
  51. % Sequence (a)
  52. h1a = 0.9.^(n-1);
  53. h2a = [0,ones(1,99),0];
  54. ha = h1a.*h2a;
  55.  
  56. % Sequence (b)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement