希尔伯特曲线及其生成

希尔伯特曲线及其生成

今天在做MATLAB社区（MATLAB cody）习题时，学习了数学中的Hilbert曲线，其能够遍历一个正方形中的所有点，是一个连续但不可导的曲线。在计算机中，它能够遍历一个(2X4^n) X (2X4^n)的二维网格，在图像处理中可以据此把一个图像映射为一个向量，后来人们将遍历的方法扩展到任意尺寸的矩形。我查阅了文献《Hilbert曲线扫描矩阵的生成算法及其MATLAB程序代码》（王笋2006），并将其中的代码做了简单实现。实现后提交到MATLAB cody，得到了在cody的2个badges，分别是speed demon和leader，分别表示第一个解决某个问题的，以及这个问题的最小size解法。

代码：

function H = HilbertScan(k)
H0 = [1 2; 4 3];
flag = 0;
for ii = 1:k-1
    if flag ==0
        H = [ H0,               (4^(ii+1)+1)*ones(2^ii) - fliplr(H0);
            (4^ii)*ones(2^ii)+H0' ,(3*4^ii+1)*ones(2^ii)-fliplr(H0')];
    else
        H = [ H0,               (4^ii)*ones(2^ii)+H0';
            (4^(ii+1)+1)*ones(2^ii)-flipud(H0) ,(3*4^ii+1)*ones(2^ii)-(fliplr(H0))'];
    end
    flag = 1 - flag;
end

[20170426]