circbitflip.m

linear time－frequency toolbox

function bitorder=circbitflip(order)
%CIRCBITFLIP  Compute sequence of numbers by flipping one bit.
%   Usage: bitorder=circbitflip(order)
%
%   CIRCBITFLIP(order) computes a sequence of numbers from 0 until
%   order-1 such that each number differs from the previous one by
%   a single flipped bits.

bitorder=[0;1];

if order==2
  return;
end;

for thisorder=2:log2(order)
  bitorder=[bitorder;...
	    2^(thisorder-1)+flipud(bitorder)];
end;