Haar wavelet transformation matrix implementation

Create Haar wavelet transformation matrix

Create Haar wavelet transformation matrix H for the matrix vector
multiplication implimentation of Haar wavelet transformation.
This function uses the following nice formula to create the Haar
transformation matrix:
H_n=1/sqrt(2)[H_(n/2) kron (1 1)
I_(n/2) kron (1 -1)],
where 'kron' denotes the kronecker product.
The iteration starts with H_1=[1]. The normalization constant 1/sqrt(2)
ensure that H_n^T*H_n=I, where I is identity matrix. Haar wavelets are the
rows of H_n.