nulovat=99
load  Bumps;
insig=Bumps;

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

%spravime DCT
spec = dct(insig);
%vynulujeme zelany pocet koeficientov
zeroed=spec;
zeroed(zerobegin:inlen)=zeros(1,inlen-zerobegin+1);
%inverzna transformacia
aprox=idct(zeroed);

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