if (~exist('filter_spike')) 
    load filter_spike;
end
    
f = 0.995; 
in = 1:1200; 
filter100 = sin(f*(in)*3.14159/20);
filter300 = sin(3*f*(in)*3.14159/20);
filter1000 = sin(10*f*(in)*3.14159/20);

for i = 1:3
    b(:,i) = A(:,i) - mean(A(:,i));
    b100(:,i) = conv(b(:,i), filter100); 
    b300(:,i) = conv(b(:,i), filter300); 
    b1000(:,i) = conv(b(:,i), filter1000); 
    cspike(:,i) = conv(b(:,i), filter_spike);
end;

figure; 
plot(A)

figure; 
subplot(4,1,1); plot(b100, '.');
subplot(4,1,2); plot(b300, '.');
subplot(4,1,3); plot(b1000, '.');
subplot(4,1,4); plot(cspike);