% beethovendwt.m displays the wavelet transform for the signal in wavelets.wav

clf;

[b]=wavread('beethoven.wav');% read "wavelets" data file
b=b-mean(b);                % center the data
h=[1 1]/sqrt(2);            % multiresolution coefficients for Haar wavelet
h=[1 2 1]/sqrt(8);          % multiresolution coefficients for triangle wavelet
                            % put the h you want to use second

showdwt(b',h,44100);