Цифровая ФАПЧ, Расчет контура Опции

[PavPro]

Доброго времени суток.
Знаю что данная тема поднималась на разных ветках форума.
Тем не менее нигде не смог найти ответ на ряд вопросов касающихся расчета цифровой ФАПЧ.
Ниже приводится алгоритм расчета контура ФАПЧ по Незами:

fclk = 10e6; % Частота дискретизации (Гц)
T = 1/fclk; % Время выборки (секунды)
M = 32; % Разрядность фазового аккумулятора NCO
fres = fclk/(2^M); % Разрешение по частоте NCO (Гц)
Ko = fres*2*pi; % Коэфф. усиления NCO (рад/сек * разряд)
Kd = 0.5; % Коэфф. усиления детектора (хрен знает в чем измеряется??)
Bn = 10e3; % Натуральная полоса пропускания ФАПЧ (Гц)
wn = Bn*2*pi; % Круговая полоса пропускания (рад/с)
zeta = 0.707; % Коэфф. демпфирования (затухания)
BL = 2*pi*Bn*(zeta+(1/(4*zeta))); % Шумовая полоса ФАПЧ
Kp = (T*(wn^2)/(2*Ko*Kd))*(1+(1/tan(T*wn/(4*zeta)))); % Пропорциональный коэфф.
Ki = T*(wn^2)/(Ko*Kd); % Интегральный коэфф.

Вопрос первый: правильно ли выполнен расчет в целом. Насколько можно доверять данной методике.
Дело в том что если считать по методике с http://www.dsplib.ru/content/dpll/dpll.html получаются совсем другие результаты
при одних и тех же изходных данных. Точнее числа получаются такими же, но они отличаются порядком запятой.
Тут как я понимаю, дело в масштабировании??

Вопрос второй. Как правильно выражается коэффициент фазового детектора?
В данном случае я принял равным Kd=0.5, изсходя из следующих соображений - при перемножении cos*sin в режиме слежения ФАПЧ,
амплитуда сигнала на выходе детектора становится ровно в половину меньше. Так ли это.
Заранее большое спасибо всем, кто поможет разобраться в данных вопросах!