From: manu Date: Tue, 27 Nov 2007 16:41:27 +0000 (+0000) Subject: testacqui enhancements X-Git-Url: https://git.parisson.com/?a=commitdiff_plain;h=e3a32bcc1b17047cc29c889a46f9c20a650f4b84;p=cnaq.git testacqui enhancements git-svn-id: http://svn.parisson.org/svn/CNAQ/trunk@122 5fc3e0e6-29bc-4d03-b52b-c088cb822bde --- diff --git a/tests/testacqui.m b/tests/testacqui.m index c4fae9a..9f18527 100644 --- a/tests/testacqui.m +++ b/tests/testacqui.m @@ -1,25 +1,33 @@ -lt=5; %s +function testacqui(n) + +if nargin < 1, n = 50; end +if n == 1, n=2; end +%signal generation f_min=20; f_max=20000; -n=50; +lt=5; t=0:1/44100:lt; sig=0.9*chirp(t,20,lt,20000,'logarithmic'); len_sig = length(sig); -% sig=wavread('excitation'); -for i = 1:n - i - pause(0.1) +%playrecord signal +for it = 1:n + pause(0.1) inputbuffer = pa_wavplayrecord(sig',1,[44100],0,[1],[1],1,['asio']); - Co=xcorr(sig',inputbuffer'); -% Co=crosscorr(sig',inputbuffer',length(inputbuffer)-1); - [Y(i),I(i)] = max(Co); - I(i)-len_sig + Co=xcorr(inputbuffer',sig'); + [Y(it),I(it)] = max(Co); + delai(it)=I(it)-len_sig; + end +mfv=mode(delai); +mfvt =num2str(mfv); +sprintf('%s','the most frequent delay is ', mfvt,' samples.') -plot(I); -axis([1 n (min(I)+2) (max(I)-2)]); +%Plot results +plot(delai); +axis([1 n (min(delai)-2) (max(delai)+2)]); figure() plot(Y); axis([1 n min(Y) max(Y)]); +