rootuss

filip31Print

Nov 3rd, 2019
208
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. %część 1:
  2.  
  3. N=4000;
  4. fp=1000;
  5. t=0:1/fp:(N-1)/fp;
  6.  
  7. %sygnał - sinusoida o zmiennej częstotliwości
  8. x=chirp(t, 100, 4.0, 400);
  9. subplot(3,2,1);
  10. plot(t,x);
  11. xlabel('czas[s]');
  12. ylabel('sygnal');
  13.  
  14.  
  15.  
  16.  
  17. %moduł widma
  18. Nf=2^12;
  19. N21=Nf/2+1;
  20. v=fft(x,Nf);
  21. wx=abs(v);
  22. f=linspace(0,fp/2,N21);
  23. subplot(3,2,2);
  24. plot(f,wx(1:N21));
  25. xlabel('czestotliwosc[Hz]');
  26. ylabel('modul widma');
  27.  
  28.  
  29.  
  30.  
  31.  
  32.  
  33. %odpowiedź impulsowa
  34. M=201;
  35. h=fir1(M-1,0.5);
  36. th=0:1/fp:(M-1)/fp;
  37. subplot(3,2,3);
  38. plot(th,h);
  39. xlabel('czas[s]');
  40. ylabel('odp. impulsowa');
  41.  
  42.  
  43.  
  44.  
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60. %moduł transmitancji
  61. v=fft(h,Nf);
  62. wh=abs(v);
  63. subplot(3,2,4);
  64. plot(f,wh(1:N21));
  65. xlabel('czestotliwosc[Hz]');
  66. ylabel('modul transmitancji');
  67.  
  68.  
  69.  
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  
  77. %sygnał po filtracji
  78. y=filter(h,1,x);
  79. subplot(3,2,5);
  80. plot(t,y);
  81. xlabel('czas[s]');
  82. ylabel('sygnal po filtracji');
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95. %widmo amplitudy po filtracji
  96. v=fft(y,Nf);
  97. wy=abs(v);
  98. subplot(3,2,6);
  99. plot(f,wy(1:N21));
  100. xlabel('czestotliwosc[Hz]');
  101. ylabel('widmo aplitudy po filtracji');
RAW Paste Data