q = quantizer('nearest',[4 3]);
[f,x] = errpdf(q);
subplot(211)
plot(x,f)
title('Computed PDF of the quantization error.')
输出绘图显示量化误差的概率密度函数。将此结果与蒙特卡罗试验的样本概率密度函数绘图进行比较:
r = realmax(q);
u = 2*r*rand(10000,1)-r; % Original signal
y = quantize(q,u); % Quantized signal
e = y - u; % Error
subplot(212)
hist(e,20)
gca.xlim = [min(x) max(x)];
title('Estimate of the PDF of the quantization error.')