f = 0.995; 
in = 1:1200; 
filter1000 = sin(1000*f*(in)*3.14159/2000);
filter316 = sin(316*f*(in)*3.14159/2000);
filter100 = sin(100*f*(in)*3.14159/2000);

for rowrange = 1:2
	rows = [1:30000]; 
	% rows = [30001:size(A,1)];
	% columns = [4:6];  
	columns = [4:10];
	% columns = [10];
	
	%figure; 
	%plot(A(rows,columns))
	
	for i = columns
        b(:,i) = A(rows,i);
        b1000(:,i) = conv(b(:,i), filter1000); 
        b316(:,i) = conv(b(:,i), filter316); 
        b100(:,i) = conv(b(:,i), filter100); 
	end;
	
	figure; 
	col = length(columns);
	for(i = columns)
        1+(i-columns(1))*col
        subplot(3,col,i-columns(1)+1); plot(b1000(:,i), '.');
        2+(i-columns(1))*col
        subplot(3,col,i-columns(1)+1+col); plot(b316(:,i), '.');
        3+(i-columns(1))*col
        subplot(3,col,i-columns(1)+1+2*col); plot(b100(:,i), '.');
	end;
	
	clear b;
	clear b1000;
	clear b316;
	clear b100; 
end