compute_histogram.m

This document is intended to serve as an introduction to Wavelet processing through a set of Matlab

function [h,x] = compute_histogram(M, options)

% compute_histogram - compute the (symmetric) histogram of a vector after threshold.
%
% [h,x] = compute_histogram(M,options);
%
%   If M contains only entry in 1...p, where p=max(M), then h is of size p
%       h(k) = #{i \ M(i)=k} / length(M)
%   If M contains only entry in -p...p, where p=max(abs(M)), then h is of size 2*p+1
%       h(k) = #{i \ M(i)=k-p-1} / length(M)
%   (you can also force this by setting options.symmetrize=1).
%
%   For real valued signals, you can add thresholding by setting options.T
%
%   Copyright (c) 2006 Gabriel Peyr