Advertisement
Guest User

Untitled

a guest
Nov 23rd, 2017
81
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.08 KB | None | 0 0
  1. close all
  2. N = 20;
  3. x = zeros(1,N);
  4.  
  5. x(1,1)=1;
  6. a=0.5;
  7. w=1;
  8.  
  9.  
  10. ha = zeros(1,N);
  11. hb = zeros(1,N);
  12. hc = zeros(1,N);
  13. hd1 = zeros(1,N);
  14. hd2 = zeros(1,N);
  15. he1 = zeros(1,N);
  16. he2 = zeros(1,N);
  17. h = zeros(1,N);
  18.  
  19. for n=1:N %obrobka wektorow systemu
  20. %ha(1,n)=(1/2)^n; %ciag dla systemu stabilnego (1/2 moze byc zastapiona inna wartoscia
  21. %ktora bedzie spelniac warunek |z|<1)Ten ciag bierze sie z Z-Transformaty
  22.  
  23. %hc(1,n)=j^n;%ciag dla systemu "generator" (j moze byc zastapione inna wartoscia
  24. %ktora spelnia warunek |z|=1, najlepiej nie 1 bo to szczegolny
  25. %przypadek
  26.  
  27. ha(1,n) = j^n; %kolo jedn
  28. hb(1,n) = (j*0.5)^n + (-j*0.5)^n; % bieguny sprzezone
  29. hc(1,n) = 25^n + 4^n; %bieguny os rzeczywista
  30.  
  31. hd1(1,n) = 5^n+(-0.5)^n;%poza kolem jedn, blizej srodka ukl wsp
  32. hd2(1,n) = 10^n+(-0.5)^n;%poza kolem jedn, dalej srodka ukl wsp
  33.  
  34. he1(1,n) = 0.1^n+0.2^n;%w kole jedn, blizej sr
  35. he2(1,n) = 0.4^n+0.55^n;%w kole jedn, dalej sr
  36.  
  37. h1(1,n) = exp(a*n) * cos(w*n);
  38. end
  39.  
  40. figure(1);%to dajemy po to zeby nam wykresy w oddzielnych oknach rysowalo
  41. stem(0:1:N-1, x(1,1:1:N), 'r*');
  42. title('nasz sygnal');
  43.  
  44. figure(2);
  45. ya = conv(ha, x);%tu robimy splot
  46. stem(0:1:N-1, ya(1,1:1:N),'r*');
  47. title('na kole jedn');
  48.  
  49. figure(3);
  50. yb = conv(hb, x);
  51. stem(0:1:N-1, yb(1,1:1:N),'r*');
  52. title('bieguny sprzezone');
  53.  
  54. figure(4);
  55. yc = conv(hc, x);
  56. stem(0:1:N-1, yc(1,1:1:N),'r*');
  57. title('bieguny os rzeczywista');
  58.  
  59. figure(5);
  60. yd1 = conv(hd1, x);
  61. yd2 = conv(hd2, x);
  62. semilogy(0:1:N-1, yd1(1,1:1:N),'r.');
  63. hold on;
  64. semilogy(0:1:N-1, yd2(1,1:1:N),'b.');
  65. title('poza kolem jedn, czer-blizej, nieb-dalej srodka ukl wsp');
  66.  
  67. figure(6);
  68. ye1 = conv(he1, x);
  69. ye2 = conv(he2, x);
  70. stem(0:1:N-1, ye1(1,1:1:N),'r.');
  71. hold on;
  72. stem(0:1:N-1, ye2(1,1:1:N),'b.');
  73. title('w kole jedn, czer-blizej, nieb-dalej srodka ukl wsp');
  74.  
  75. hv = fft(h1);
  76.  
  77. % ampl = abs(hvv);
  78. % faza = angle(hvv);
  79.  
  80. figure(7);
  81. %y = conv(h, x);
  82. hold on;
  83. stem(0:1:N-1, abs(hv),'r.');
  84. title('amplituda');
  85. figure(8);
  86. hold on;
  87. stem(0:1:N-1, faza(1,1:1:N),'r.');
  88. title('faza');
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement