b?=?randperm（8）-1;?x?=?sort（b）;
tic
ffn?=?fft（x）;??
t?=?toc
disp（ffn）;
m=log2（2^nextpow2（length（x）））;
N=2^m;
if?length（x）??x=[xzeros（1N-length（x））];
end
nxd?=?bin2dec（fliplr（dec2bin（[1:N]-1m）））+1;
y?=?x（nxd）;
tic
for?M=1:m