haar_1d.m

《数字信号处理：基于计算机的方法（第三版）》资源下载第三部分。 著者： （美）米特拉著 作译者： 孙洪等译 ISBN号： 7-121-02652-X 出版日期： 2006-06 出版

function [Y] = haar_1D(X)
% Computes the Haar transform of the input vector X
% The length of X must be a power-of-2
% Recursively builds the Haar matrix H
v = log2(length(X)) - 1;
H = [1 1;1 -1];
for m = 1:v,
A = [kron(H,[1 1]);
2^(m/2).*kron(eye(2^m),[1 -1])];
H = A;
end
Y = H*double(X(:));