
wave0Aver(1:3000,1)=0;
wave6aAver(1:3000,1)=0;
data=[wave0Aver wave1Aver wave2Aver wave3Aver wave4Aver wave5Aver wave6aAver];
x=data;
nCh=size(x,2);
step=2;
csdch=[step+1:nCh-step];
csd=-(x(:,csdch+step) - 2*x(:,csdch) + x(:,csdch-step));
interp_csd=interp2(csd,2,'linear');
x=[0 130]; y=[0 6];
imagesc(x,y,interp_csd');
hold on;
t=(130/3000:130/3000:130);
plot(t,-wave2Aver*50000+3,'k');
