Дискриминатор частоты с временным сдвигом квадратурных компонент

Описание дискриминатора

Дискриминатор использует отсчеты коррелятора с текущего и предыдущего такта работы.

$u_{d \omega } = I_k(\widetilde{\tau}_k,\widetilde{\omega}_{d\,k})Q_{k-1}(\widetilde{\tau}_{k-1},\widetilde{\omega}_{d\,k-1}) - Q_k(\widetilde{\tau}_k,\widetilde{\omega}_{d\,k})I_{k-1}(\widetilde{\tau}_{k-1},\widetilde{\omega}_{d\,k-1})$,

где
$I_k(\widetilde{\tau}_k,\widetilde{\omega}_{d\,k}) = \sum_{l=1}^{L}y(t_{k,l})$ $I_k(\widetilde{\tau}_k,\widetilde{\omega}_{d\,k}) = \sum_{l=1}^{L}y(t_{k,l})h_{дк}(t_{k,l}-\widetilde{\tau}_k)cos(\omega_0t_{k,l}+\widetilde{\omega}_{d\,k}(l-1)T_d))$
$Q_k(\widetilde{\tau}_k,\widetilde{\omega}_{d\,k}) = \sum_{l=1}^{L}y(t_{k,l})h_{дк}(t_{k,l}-\widetilde{\tau}_k)sin(\omega_0t_{k,l}+\widetilde{\omega}_{d\,k}(l-1)T_d))$

В процессе выполнения магистерской работы возникла необходимость использовать данный частотный дискриминатор. Для этого были проверены его характеристики, указанные в книге ГЛОНАСС (стр. 262 - 264).