blp.m

Double Density Wavelet Soft

function h = blp(k)
% h = blp(k)
% Binomial lowpass filter
% H(z) = [(1/z+1)/2]^k

h = 1;
for i = 1:k
  h = ([h 0] + [0 h])/2;
end