dr=dir('*.wav');
ss=wavread(dr(1).name);
size(ss)
ss=ss(:,1);
plot(ss(1:100))
plot(ss(1:1000))
plot(ss(1:10000))
plot(ss((1:10000)+1000))
plot(ss((1:10000)+10000))
plot(ss((1:10000)+100000))
plot(ss((1:100000)+100000))
plot(ss((1:100000)+1000000))
plot(ss((1:100000)))
plot(ss((1:1000000)))
plot(ss((1:1000000)+2e4))
plot(ss((1:1000000)+2e5))
plot(ss((1:10000)+2e5))
plot(ss((1:10000)+3e5))
plot(ss((1:100000)+3e5))
plot(ss((1:100000)+3.1e5))
plot(ss((1:100000/5)+3.1e5))
sss=ss(1:10:end);
plot(sss((1:10000/5)+3.1e4))
sss=ss(1:100:end);
plot(sss((1:1000/5)+3.1e3))
plot(conv2(sss,sss((1:1000/5)+3.1e3),'same'))
plot(conv2(sss,fliplr(sss((1:1000/5)+3.1e3)),'same'))
plot(conv2(sss,flipud(sss((1:1000/5)+3.1e3)),'same'))
plot([conv2(sss,flipud(sss((1:1000/5)+3.1e3)),'same') sss])
plot([conv2(sss(1:10000),flipud(sss((1:1000/5)+3.1e3)),'same') sss(1:10000)])
plot([conv2(sss(1:10000),fliplr(sss((1:1000/5)+3.1e3)),'same') sss(1:10000)])
plot([conv2(sss((1:1000/5+3.1e3)),fliplr(sss((1:1000/5)+3.1e3)),'same') conv2(sss((1:1000/5+3.1e3))])
plot([conv2(sss((1:1000/5+3.1e3)),fliplr(sss((1:1000/5)+3.1e3)),'same') sss((1:1000/5+3.1e3))])
plot(sss((1:1000/5)+3.1e3))
plot([conv2(sss((1:1000/5+3.1e3)),fliplr(sss((1:1000/5)+3.1e3)),'same') sss((1:1000/5+3.1e3))])
sss((1:1000/5+3.1e3))
plot([conv2(sss((1:1000/5)+3.1e3),fliplr(sss((1:1000/5)+3.1e3)),'same') sss((1:1000/5+3.1e3))])
plot(sss((1:1000/5)+3.1e3))
plot(conv2(sss((1:1000/5)+3.1e3),sss((1:1000/5)+3.1e3),'same'))
plot(conv2(sss((1:1000/5)+3.1e3),fliplr(sss((1:1000/5)+3.1e3)),'same'))
plot(conv2(sss((1:1000/5)+3.1e3),flipud(sss((1:1000/5)+3.1e3)),'same'))
plot(conv2(sss((1:10000/5)+3.1e3),flipud(sss((1:1000/5)+3.1e3)),'same'))
plot(sss((1:1000/5)+3.1e3))
plot(sss((1:200)+3.1e3))
plot(sss((1:200)+3.1e3+200))
plot(sss((1:200)+3.1e3+400))
plot(sss((1:200)+3.1e3))
plot(ss((1:20000)+3.1e5))
plot(ss((1:20000)+3.1e5+20000))
plot(ss((1:20000)+3.1e5))
plot(conv2(ss((1:20000)+3.1e5),ss((1:20000)+3.1e5),'same'))
plot(conv2(ss((1:20000)+3.1e5),flipud(ss((1:20000)+3.1e5)),'same'))
plot(conv2(ss((1:200000)+3.1e5),flipud(ss((1:20000)+3.1e5)),'same'))
plot(conv2(ss((1:20000)+3.1e5),flipud(ss((1:20000)+3.1e5)),'same'))
plot(conv2(ss((1:40000)+3.1e5-20000),flipud(ss((1:20000)+3.1e5)),'same'))
plot(conv2(ss((1:60000)+3.1e5-20000),flipud(ss((1:20000)+3.1e5)),'same'))
plot(conv2(ss((1:60000)+3.1e5-20000),conv2(flipud(ss((1:20000)+3.1e5)),ones(1,3),'same'),'same'))
plot(conv2(ss((1:60000)+3.1e5-20000),conv2(flipud(ss((1:20000)+3.1e5)),ones(1,1),'same'),'same'))
plot(conv2(ss((1:60000)+3.1e5-20000),conv2(flipud(ss((1:20000)+3.1e5)),ones(3,1),'same'),'same'))
plot(conv2(ss((1:60000)+3.1e5-20000),conv2(flipud(ss((1:20000)+3.1e5)),ones(5,1),'same'),'same'))
plot(conv2(ss((1:60000)+3.1e5-20000),conv2(flipud(ss((1:20000)+3.1e5)),ones(7,1),'same'),'same'))
plot(conv2(ss((1:60000)+3.1e5-20000),conv2(flipud(ss((1:20000)+3.1e5)),ones(9,1),'same'),'same'))
plot(conv2(ss((1:60000)+3.1e5-20000),conv2(flipud(ss((1:20000)+3.1e5)),ones(1,1),'same'),'same'))
plot(conv2(ss((1:20000)+3.1e5),flipud(ss((1:20000)+3.1e5)),'same'))
plot(conv2(ss((1:60000)+3.1e5-20000),flipud(ss((1:20000)+3.1e5)),'same'))
(conv2(ss((1:60000)+3.1e5-20000),flipud(ss((1:20000)+3.1e5)),'same'));
ssp=0*ss;
tmp=(conv2(ss((1:60000)+3.1e5-20000),flipud(ss((1:20000)+3.1e5)),'same'));
ssp((1:20000)+0)=tmp(20001:40000);
NN=length(ss)/20000
NN=fix(length(ss)/20000)
for n=1:NN
tmp=(conv2(ss((1:60000)+3.1e5+20000(n-2)),flipud(ss((1:20000)+3.1e5)),'same'));
for n=1:NN
tmp=(conv2(ss((1:60000)+3.1e5+20000*(n-2)),flipud(ss((1:20000)+3.1e5)),'same'));
ssp((1:20000)+(n-1)*20000)=tmp(20001:40000);
plot(tmp);
drawnow;
n
end
plot(ssp)
size(ssp)(
clf
size(ssp)
size(ss)
plot(ss((1:20000)+3.1e5))
plot(sp((1:20000)))
plot(ssp((1:20000)))
plot(ssp((1:40000)))
plot(ssp((1:100000)))
plot(ssp((1:1000000)))
plot(ssp((1:10000000)))
plot(abs(ssp((1:10000000))))
plot(abs(ssp((1:1000000))))
plot(ssp((1:40000)))
plot(ssp((1:100000)))
cplot(ssp((1:100000)).*exp(pi*i*(1:100000)))
cplot(ssp((1:100000)).*exp(pi*i*(1:100000)'))
cplot(ssp((1:100000)).*exp(0*pi*i*(1:100000)'))
cplot(ssp((1:100000)).*exp(pi*i*(1:100000)'))
plot(ssp((1:100000)).*exp(pi*i*(1:100000)'))
plot(real(ssp((1:100000)).*exp(pi*i*(1:100000)')))
plot(real(ssp((1:100000)).*exp(0*pi*i*(1:100000)')))
plot(real(ssp((1:100000)).*exp(pi*i*(1:100000)')))
plot(real(ssp((1:100000)).*exp(pi*i*(1:100000)'+.1)))
plot(real(ssp((1:100000)).*exp(pi*i*(1:100000)'+.2)))
plot(real(ssp((1:100000)).*exp(pi*i*(1:100000)'+.01)))
plot(real(ssp((1:100000)).*exp(pi*i*(1:100000)')))
plot(exp(pi*i*(1:100000)'))
cplot(exp(pi*i*(1:100000)'))
cplot(exp(pi*i*(1:100000)))
plot(real(ssp((1:100000)).*exp(pi*i*(1:100000)'+.01)))
plot(real(ssp((1:100000)).*exp(pi/2*i*(1:100000)')))
cplot(exp(pi*i*(1:100000)/2))
cplot(exp(pi*i*(1:100000)))
plot(real(ssp((1:100000)).*exp(pi/2*i*(1:100000)')))
plot(real(ssp((1:100000)).*exp(pi*i*(1:100000)')))
plot(real(ssp((1:100000)).*exp(pi*i*(1:100000)'))); axis([1e4-100 1e4+100 -1 1])
plot(real(ssp((1:100000)).*exp(0*pi*i*(1:100000)'))); axis([1e4-100 1e4+100 -1 1])
plot(real(ssp((1:100000))); axis([1e4-100 1e4+100 -1 1])
plot(real(ssp((1:100000)))); axis([1e4-100 1e4+100 -1 1])
plot(real(ssp((1:100000))));
plot(real(ssp((1:100000))),'.');
plot(real(ssp((1:100000))),'.'); axis([1e4-100 1e4+100 -1 1])
plot(real(ssp((1:100000))),'.-'); axis([1e4-100 1e4+100 -1 1])
plot(real(ss((1:100000))),'.-'); axis([1e4-100 1e4+100 -1 1])
plot(real(ss((1:100000)+3.1e5)),'.-'); axis([1e4-100 1e4+100 -1 1])
plot(ss((1:100000)+3.1e5),'.-'); axis([1e4-100 1e4+100 -1 1])
plot(ss((1:100000)+3.1e5-20000),'.-'); axis([1e4-100 1e4+100 -1 1])
plot(ss((1:100000)+3.1e5-20000),'.-'); axis([1e4-100 1e4+100 -.1 .1])
plot(ss((1:100000)+3.1e5-20000),'.-'); axis([1e4-100 1e4+100 -.01 .01])
plot(ss((1:20000)+3.1e5))
plot(ss((1:100000)+3.1e5-20000),'.-'); axis([1e4-100 1e4+100 -.01 .01])
plot(conv2(ss((1:60000)+3.1e5+20000*(n-2)),flipud(ss((1:20000)+3.1e5)),'same'));
plot(conv2(ss((1:60000)+3.1e5+20000*(1-2)),flipud(ss((1:20000)+3.1e5)),'same'));
plot(conv2(ss((1:60000)+3.1e5+20000*(1-2)),flipud(ss((1:20000)+3.1e5)),'same')); axis([3e4-100 3e4+100 -.01 .01])
plot(conv2(ss((1:60000)+3.1e5+20000*(1-2)),flipud(ss((1:20000)+3.1e5)),'same')); axis([3e4-100 3e4+100 -1 1])
plot(conv2(ss((1:60000)+3.1e5+20000*(1-2)),flipud(ss((1:20000)+3.1e5+1)),'same')); axis([3e4-100 3e4+100 -1 1])
plot(abs(fftshift(fft(tmp))))
plot(-30000:30000,abs(fftshift(fft(tmp))))
plot(-30000:30000-1,abs(fftshift(fft(tmp))))
plot((-30000:30000-1)/60000,abs(fftshift(fft(tmp))))
plot((-30000:30000-1)/60000*44100,abs(fftshift(fft(tmp))))
plot(tmp)
plot(ss((1:100000)+3.1e5-20000),'.-');
plot(ss((1:100000)+3.1e5-20000));
plot(tmp)
plot(tmp,'.-')
44100/17
1/(44100/17)
for n=1:NN
tmp=(conv2(ss((1:60000)+3.1e5+20000*(n-2)),flipud(ss((1:20000)+3.1e5)),'same'));
ssp((1:20000)+(n-1)*20000)=tmp(20001:40000);
n
end
plot(tmp)
plot(abs(fftshift(fft(fftshift(tmp)))))
ftmp=abs(fftshift(fft(fftshift(tmp))));
plot(ftmp)
plot(ftmp(1:end/2-100))
plot(ftmp(1:end/2-1000))
max(ftmp(1:end/2-1000))
[mx pp]=max(ftmp(1:end/2-1000));
pp
pp-60000/2
plot(tmp)
plot([tmp exp(2*pi*i*7811/60000*(1:60000))])
plot([tmp exp(2*pi*i*7811/60000*(1:60000)')])
plot(real([tmp exp(2*pi*i*7811/60000*(1:60000)'))])
plot(real([tmp exp(2*pi*i*7811/60000*(1:60000)')]))
plot(real([tmp exp(2*pi*i*7811/60000*(1:60000)')/10]))
plot(real([tmp.*exp(2*pi*i*7811/60000*(1:60000)')/10]))
plot(abs([tmp.*exp(2*pi*i*7811/60000*(1:60000)')/10]))
plot(abs([tmp.*exp(2*pi*i*7811/60000*(1:60000)')/10])); axis([30540-100 30540+100 -1 1])
plot(abs([tmp.*exp(2*pi*i*7811/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(abs([tmp.*exp(2*pi*i*(-7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(abs([tmp.*exp(pi*i*(-7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(abs([tmp.*exp(0*pi*i*(-7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(abs([tmp.*exp(2*pi*i*(-7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
clf
plot(abs([tmp.*exp(2*pi*i*(-7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(abs([tmp exp(2*pi*i*(-7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(abs([tmp; exp(2*pi*i*(-7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(abs([tmp; exp(2*pi*i*(-7811)/60000*(1:60000)')]));
plot(abs([tmp; exp(0*pi*i*(-7811)/60000*(1:60000)')]));
plot(abs([tmp exp(0*pi*i*(-7811)/60000*(1:60000)')]));
plot(([tmp exp(0*pi*i*(-7811)/60000*(1:60000)')]));
plot(real([tmp exp(2*pi*i*(-7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(real([tmp exp(2*pi*i*(7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(real([tmp.*exp(2*pi*i*(7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(imag([tmp.*exp(2*pi*i*(7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(real([tmp exp(2*pi*i*(7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(imag([tmp.*exp(2*pi*i*(7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(real([tmp exp(2*pi*i*(7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(imag([tmp.*exp(2*pi*i*(7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(real([tmp exp(2*pi*i*(7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(real([tmp.*exp(2*pi*i*(7811)/60000*(1:60000)')])); axis([30540-100 30540+100 -1 1])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(9,1),'same')])); axis([30540-100 30540+100 -1 1])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(9,1),'same')])/9); axis([30540-100 30540+100 -1 1])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(11,1),'same')])/9); axis([30540-100 30540+100 -1 1])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(13,1),'same')])/9); axis([30540-100 30540+100 -1 1])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(17,1),'same')])/9); axis([30540-100 30540+100 -1 1])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(19,1),'same')])/19); axis([30540-100 30540+100 -1 1])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(21,1),'same')])/21); axis([30540-100 30540+100 -1 1])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(19,1),'same')])/19); axis([30540-100 30540+100 -1 1])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(19,1),'same')])/19); axis([30540-100 30540+100 -inf inf])
plot(imag([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(19,1),'same')])/19); axis([30540-100 30540+100 -inf inf])
plot(real([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(19,1),'same')])/19); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(19,1),'same')])/19); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(19,1),'same')])/19);
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(19,1),'same')])/19); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(23,1),'same')])/23); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),ones(29,1),'same')])/29); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),([1:9 8:-1:1],'same')])/9); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),[1:9 8:-1:1],'same')])/9); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),[1:9 8:-1:1]','same')])/9); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),([1:9 8:-1:1].^2)','same')])/9); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),[1:19 18:-1:1]','same')])/9); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),[1:19 18:-1:1]','same')])/9);
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),[1:19 18:-1:1]','same')])/9); axis([30540-100 30540+100 -inf inf])
tmp=(conv2(ss((1:60000)+3.1e5+20000*(n-2)),flipud(ss((1:20000)+3.1e5)),'same'));
n
n=100
tmp=(conv2(ss((1:60000)+3.1e5+20000*(n-2)),flipud(ss((1:20000)+3.1e5)),'same'));
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),[1:19 18:-1:1]','same')])/9); axis([30540-100 30540+100 -inf inf])
plot(abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),[1:19 18:-1:1]','same')])/9);
for n=1:NN
tmp=(conv2(ss((1:60000)+3.1e5+20000*(n-2)),flipud(ss((1:20000)+3.1e5)),'same'));
tmp=abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),[1:19 18:-1:1]','same')]));
for n=1:NN
tmp=(conv2(ss((1:60000)+3.1e5+20000*(n-2)),flipud(ss((1:20000)+3.1e5)),'same'));
tmp=abs([conv2(tmp.*exp(2*pi*i*(7811)/60000*(1:60000)'),[1:19 18:-1:1]','same')]);
ssp((1:20000)+(n-1)*20000)=tmp(20001:40000);
plot(tmp);
drawnow;
n
end
NN
save
clear all
trackprep000
clear all
trackprep000
edit trackprep000.m
trackprep000
clear all
trackprep000
plot(mean(aimg,2))
plot(mean(aimg,2)/max(mean(aimg,2))
plot(mean(aimg,2)/max(mean(aimg,2)))
trackprep000
simage(ledimg); caxis([0 6]); hold on; plot(ly,lx,'w.'); hold off
simage(ledimg); caxis([0 6]);
simage(aimg)
simage(mn(box(1):box(2),box(3):box(4)))
plot(sum(aimg(:,y2+bd:end))/max(sum(aimg(:,y2+bd:end))))
simage(aimg(:,y2+bd:end)))
simage(aimg(:,y2+bd:end))
simage(smx)
trackprep000
simage(ledimg); caxis([0 6]);
simage(ledimg)
simage(smx)
simage(aimg(:,y2+bd:end))
colormap(jet)
simage(ledimg)
simage(smx)
simage(ledimg)
simage(ledimg>10)
simage(ledimg>100)
simage(ledimg>1)
simage(ledimg)
colorbar
simage((1./chiimg(smx/max(smx(:)),ipf(abs(rr),D,.8),[])))
simage((1./chiimg(smx/max(smx(:)),ipf(abs(rr),D,.8),[]))); colorbar
simage((1./chiimg(smx/max(smx(:)),ipf(abs(rr),D,.8),[]))<10); colorbar
simage((1./chiimg(smx/max(smx(:)),ipf(abs(rr),D,.8),[]))>10); colorbar
simage((1./chiimg(smx/max(smx(:)),ipf(abs(rr),D,.8),[]))>6); colorbar
simage((1./chiimg(smx/max(smx(:)),ipf(abs(rr),30,.8),[]))>6); colorbar
simage((1./chiimg(smx/max(smx(:)),ipf(abs(rr),40,.8),[]))>6); colorbar
simage((1./chiimg(smx/max(smx(:)),ipf(abs(rr),50,.8),[]))>6); colorbar
simage((1./chiimg(smx/max(smx(:)),ipf(abs(rr),50,.8),[]))); colorbar
ledimg=0;
fprintf('Finding LED centers');
for D=30:2:50
ledimg=max(ledimg,(1./chiimg(smx/max(smx(:)),ipf(abs(rr),D,.8),[])));
fprintf('.');
end
simage(ledimg)
colorbar
simage(ledimg)
colorbar
[Np lx ly]=findpeaks(ledimg,1,6,20); Np
simage(ledimg>6)
simage(ledimg); hold on; plot(ly,lx,'w.'); hold off
simage(ledimg)
simage(ledimg); hold on; plot(ly,lx,'w.'); hold off
[Np lx ly]=findpeaks(ledimg,1,6,30); Np
simage(ledimg); hold on; plot(ly,lx,'w.'); hold off
[Np lx ly]=findpeaks(ledimg,1,6,40); Np
simage(ledimg); hold on; plot(ly,lx,'w.'); hold off
[Np lx ly]=findpeaks(ledimg,1,7,40); Np
simage(ledimg); hold on; plot(ly,lx,'w.'); hold off
trackprep000
%-- 4/7/10  6:15 PM --%
load('matlab.mat','ssp');
plot(ssp(1:100000))
plot(ssp(1:100000)>80)
plot(ssp(ssp(1:100000)>80))
plot(ssp(ssp>80))
i1=find(diff(ssp>80)==1);
i2=find(diff(ssp>80)==-1);
size(i1)
size(i2)
plot((i1+i2)/2)
plot(diff((i1+i2)/2))
plot(diff((i1+i2)/2)/44100)
plot(n,diff((i1+i2)/2)/44100)
n=1:1770;
plot(n,diff((i1+i2)/2)/44100,n,1)
plot(n,diff((i1+i2)/2)/44100,n,1+n*0)
plot(n,diff((i1+i2)/2)/44100)
size(i1)
plot(n',diff((i1+i2)/2)/44100,n,1+n*0)
plot(n,diff((i1+i2)/2)/44100,n,1+n*0)
n=1:1770-1;
plot(n,diff((i1+i2)/2)/44100,n,1+n*0)
tt=(i1+i2)/2/44100;
plot(tt(n),diff((i1+i2)/2)/44100,n,1+n*0)
plot(tt(n),diff((i1+i2)/2)/44100,tt(n),1+n*0)
hist(diff(tt),100)
hist(diff(tt),1000)
hist(diff(tt),10000)
plot(tt(n),diff((i1+i2)/2)/44100,tt(n),1+n*0)
hist(diff(tt),10000)
hist(diff(tt),100000)
hist(diff(tt),1000000)
2^32
2^64
(2^64)/1e9
(2^64)/1e9/3600
(2^64)/1e9/3600/24
(2^64)/1e9/3600/24/365