+
−%
+
−% Test delay estimation with two cases: injected signal or only noise
+
−% 
+
−% M Nofrarias 13-10-09
+
−%
+
−% VERSION:    $Id: test_ao_delayEstimate.m,v 1.1 2009/10/21 08:39:28 miquel Exp $
+
−%
+
−%
+
−%%
+
−
+
−close all
+
−clear
+
−A =1;
+
−sigma = 0.01;
+
−
+
−% test 1
+
−fs = 2048;
+
−nsecs = 0.75;
+
−sig_f = 256;
+
−sig_d = 2;
+
−
+
−% test 2
+
−% fs = 10;
+
−% nsecs = 10;
+
−% sig_f = 1;
+
−% sig_d = 3;
+
−
+
−signal = 1;
+
−
+
−%%
+
−
+
−
+
−if signal == 1
+
−  % create white noise, given a known seed
+
−  pl = plist('waveform','noise','nsecs',nsecs,'fs',fs,'sigma',sigma);
+
−  a1 = ao(pl);
+
−  a1.setName('noise')
+
−  
+
−  % create sine wave
+
−  pl = plist('waveform','sine wave','A',A,'f',sig_f,'phi',0,'nsecs',nsecs,'fs',fs);
+
−  a2 = ao(pl);
+
−  
+
−  % add both
+
−  a = a1 +a2 ;
+
−  
+
−  % create white noise, given a known seed
+
−  pl = plist('waveform','noise','nsecs',nsecs,'fs',fs,'sigma',sigma);
+
−  b1 = ao(pl);
+
−  b1.setName('noise')
+
−  
+
−  b2 = delay(a2,plist('n',sig_d));
+
−  
+
−  b =  b1 +b2 ;
+
−else
+
−  % create white noise, given a known seed
+
−  pl = plist('waveform','noise','nsecs',nsecs,'fs',fs,'sigma',sigma);
+
−  a1 = ao(pl);
+
−  a1.setName('noise')
+
−  
+
−  
+
−  % add both
+
−  a = a1 ;
+
−  
+
−  % create white noise, given a known seed
+
−  pl = plist('waveform','noise','nsecs',nsecs,'fs',fs,'sigma',sigma);
+
−  b1 = ao(pl);
+
−  b1.setName('noise')
+
−  
+
−  b = delay(a,plist('n',sig_d));
+
−  
+
−end
+
−
+
−iplot(a,b)
+
−
+
−
+
−%%
+
−
+
−pl = plist('weight','scc','plot','true');
+
−d(1) = delayEstimate(a,b,pl);
+
−
+
−
+
−
+
−