% Wavelet Analyse des lokalen Hoelder-Exponenten von 
% 1-|x| auf dem Intervall [-1,1]

x = linspace(-1,1-1/256,512);
f = 1-abs(x);

c = cwt(f,1:64,'sym2','abslvl',[100 400]);
% Abschneiden der Raender bei der Darstellung der cwt

u=diag(1./log(abs([1:64]/512)))*log(abs(c))-.5;
% Berechnung in Bemerkung 2.24

u=u(2:2:end,2:2:end);
% Matlab produziert fuer ungerade shifts und ungerade
% Skalen relativ grosse Rundungsfehler, diese 
% werden aussortiert

figure(1)
mesh(u)

figure(3)
for i=1:8
subplot(4,2,i)
plot(u(:,125+i))
end

figure(3)
for i=1:8
subplot(4,2,i)
plot(u(:,50+i))
end