haar_kron 用克罗内克积求haar 矩阵

function y=haar_kron(N)
%  2^n=N ,n=1,2,3·······
y=1;
a=[1,1];
b=[1,-1];
k=1;
while(k<N)
yy=diag(ones(1,k)); 
y=[kron(y,a);kron(yy,b)];
k=k*2;
end
for i=1:N
y(i,1:end)=y(i,1:end)/sqrt(sum(abs(y(i,1:end)')));
% 
end
end

这里有haar 矩阵的的另一种算法原文地址http://blog.csdn.net/abcd1992719g/article/details/25332211#comments

参考文献：http://zh.wikipedia.org/wiki/%E5%93%88%E7%88%BE%E5%B0%8F%E6%B3%A2%E8%BD%89%E6%8F%9B

请多多指教，谢谢。