|
|
Я проверил (матлабовский код см. в конце). Ошибка получается 3 младших бита. Если взять синус на частоте 1.8 кГц, то амплитуда синуса вместо входной 20*32=640 будет равна 0.015*640=9.6 (0.015 - это подавление на 1.8 кГц). Он будет как раз чуть выше шума. А про 0.18 кГц я вообще молчу.
clear;
c_N = 200;
FR = 17;
yp = zeros(1, c_N);
yp(1) = 1;
x = zeros(1, c_N);
round_factor = 0.5;
a=131071/131072;
b=117340/131072;
for n=2:c_N,
yp(n)=(floor(((a*(x(n)-x(n-1))+b*yp(n-1))*2^FR) + round_factor))/2^FR;
end;
|
|
