Advertisement
Guest User

Untitled

a guest
Aug 19th, 2017
78
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.71 KB | None | 0 0
  1. %izničim stereo -> delam z y(:,1)
  2. %[C,L] = wavedec(y(:,1),5,'db4');
  3.  
  4. [a1,d1] = dwt(y(:,1),'db4');
  5. [a2,d2] = dwt(a1,'db4');
  6. [a3,d3] = dwt(a2,'db4');
  7. [a4,d4] = dwt(a3,'db4');
  8. [a5,d5] = dwt(a4,'db4');
  9. clear a1 a2 a3 a4
  10.  
  11. d1 = resample(d1,100,50);
  12. d2 = resample(d2,100,25);
  13. d3 = resample(d3,100,12);
  14. d4 = resample(d4,100,6);
  15. d5 = resample(d5,100,3);
  16. a5 = resample(a5,100,3);
  17.  
  18. %porežem signal
  19. d1 = d1(1:size(y,1));
  20. d2 = d2(1:size(y,1));
  21. d3 = d3(1:size(y,1));
  22. d4 = d4(1:size(y,1));
  23. d5 = d5(1:size(y,1));
  24. a5 = a5(1:size(y,1));
  25.  
  26. d1 = ovojnica(d1);
  27. d2 = ovojnica(d2);
  28. d3 = ovojnica(d3);
  29. d4 = ovojnica(d4);
  30. d5 = ovojnica(d5);
  31. a5 = ovojnica(a5);
  32.  
  33. %seštevek signalov
  34. A = d1 + d2 + d3 + d4 + d5 + a5;
  35. clear d1 d2 d3 d4 d5 a5
  36. A = A';
  37. %figure; plot(A);
  38.  
  39. %avtokorelacija
  40. acorr = xcorr([A A A],length(A),'coeff');
  41. clear A
  42. %figure; plot(acorr);
  43.  
  44. %odstranim vrhe pod 1/40 max-imuma
  45. acorrMax = max(acorr) / 40;
  46. % A=[]; k = 1;
  47. % for n = 1:length(acorr)
  48. % if acorr(n) >acorrMax
  49. % A(k) = acorr(n);
  50. % k = k+1;
  51. % end
  52. % end
  53.  
  54. %iskanje vrhov
  55. [pks,locs] = findpeaks(acorr,'minpeakdistance',5);
  56. rez = diff(locs);
  57.  
  58. razmerje = length(y) / length(acorr);
  59. razmerje = razmerje * Fs;
  60.  
  61. rez = rez * razmerje;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement