act_wavelet='bior2.2'
decomp_levels=9;
nulovat=99
load Bumps;
insig=Bumps;

inlen=length(insig)
zerobegin=ceil((1-nulovat/100)*inlen);

%spravime DWT
[spec,l] = wavedec(insig,decomp_levels,act_wavelet);
%vynulujeme zelany pocet koeficientov
zeroed=spec;
zeroed(zerobegin:length(spec))=zeros(1,length(spec)-zerobegin+1);
%inverzna transformacia
aprox=waverec(zeroed,l,act_wavelet);

subplot(4,1,1); plot(insig);  Ylabel('Vstupny signal');

subplot(4,1,2); plot(spec);   Ylabel('Spektrum');

mytitle= ['Nulované spektrum(',int2str(nulovat),'%)'];
subplot(4,1,3); plot(zeroed); Ylabel(mytitle)

subplot(4,1,4); plot(aprox);  Ylabel('Rekonštruovaný signál');

%vypocitame MSE
MSE=sum((insig-aprox)^2)/inlen