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 c4fae9a3f2666ad010ed25c004658ecc3ec0b3ed..9f18527a7d388dbc5c623762cded22a7cd06a24a 100644 (file)
--- 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)]);
+