function X=shift_beat(delta,X)		
	beatpos=X(:,4);
	if any(diff(beatpos)<0)
		disp('WARNING: metrical position not monotonic');
	end
	X(:,4)=beatpos+delta;
end

