t = linspace(0, 1, 44000);
signaali = sin(2*pi*100*t);
kvantisoitu_signaali = round(signaali);
kvantisointivirhe = signaali-kvantisoitu_signaali;
ditheroity_signaali = round(signaali+rand(1,44000)-0.5);
dither_kvantisointivirhe = signaali-ditheroity_signaali;
dither_rms = sum(dither_kvantisointivirhe.^2)
pcm_rms = sum(kvantisointivirhe.^2)
subplot(2,1,1);
plot(t, signaali,'b','Linewidth',3);
hold on
plot(t, kvantisoitu_signaali,'r','Linewidth',3);
plot(t, ditheroity_signaali,'k');
legend({'Signaali','PCM','PCM+dither'});
xlim([0 0.01]);
subplot(2,1,2);
v1 = abs(fft(kvantisointivirhe).^2);
v2 = abs(fft(dither_kvantisointivirhe).^2);
hold on
plot(v1,'r');
plot(v2,'k');
legend({'PCM','PCM+dither'});
title('Kvantisointikohinan taajuusspektri');
xlim([0 5000]);