Mean of quantization error
m = errmean(q)
m = errmean(q) returns the mean of a uniformly
distributed random quantization error that arises from quantizing
a signal by quantizer object q.
Note
The results are not exact when the signal precision is close
to the precision of the quantizer.
Find m, the mean of the quantization error
for quantizer q:
q = quantizer; m = errmean(q)
m =
-1.525878906250000e-05Now compare m to m_est,
the sample mean from a Monte Carlo experiment:
r = realmax(q); u = 2*r*rand(1000,1)-r; % Original signal y = quantize(q,u); % Quantized signal e = y - u; % Error m_est = mean(e) % Estimate of the error mean
m_est =
-1.526738835715480e-05