Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- %izničim stereo -> delam z y(:,1)
- %[C,L] = wavedec(y(:,1),5,'db4');
- [a1,d1] = dwt(y(:,1),'db4');
- [a2,d2] = dwt(a1,'db4');
- [a3,d3] = dwt(a2,'db4');
- [a4,d4] = dwt(a3,'db4');
- [a5,d5] = dwt(a4,'db4');
- clear a1 a2 a3 a4
- d1 = resample(d1,100,50);
- d2 = resample(d2,100,25);
- d3 = resample(d3,100,12);
- d4 = resample(d4,100,6);
- d5 = resample(d5,100,3);
- a5 = resample(a5,100,3);
- %porežem signal
- d1 = d1(1:size(y,1));
- d2 = d2(1:size(y,1));
- d3 = d3(1:size(y,1));
- d4 = d4(1:size(y,1));
- d5 = d5(1:size(y,1));
- a5 = a5(1:size(y,1));
- d1 = ovojnica(d1);
- d2 = ovojnica(d2);
- d3 = ovojnica(d3);
- d4 = ovojnica(d4);
- d5 = ovojnica(d5);
- a5 = ovojnica(a5);
- %seštevek signalov
- A = d1 + d2 + d3 + d4 + d5 + a5;
- clear d1 d2 d3 d4 d5 a5
- A = A';
- %figure; plot(A);
- %avtokorelacija
- acorr = xcorr([A A A],length(A),'coeff');
- clear A
- %figure; plot(acorr);
- %odstranim vrhe pod 1/40 max-imuma
- acorrMax = max(acorr) / 40;
- % A=[]; k = 1;
- % for n = 1:length(acorr)
- % if acorr(n) >acorrMax
- % A(k) = acorr(n);
- % k = k+1;
- % end
- % end
- %iskanje vrhov
- [pks,locs] = findpeaks(acorr,'minpeakdistance',5);
- rez = diff(locs);
- razmerje = length(y) / length(acorr);
- razmerje = razmerje * Fs;
- rez = rez * razmerje;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement