Где ошибка в расчете ФАПЧ ?

goshka
Сообщения: 1
Зарегистрирован: 13 фев 2021, 13:57

Где ошибка в расчете ФАПЧ ?

Сообщение goshka »

Расчитываю по методике: http://www.dsplib.ru/content/dpll/dpll.html

Код: Выделить всё

 Fs = 30.72e6; % частота дискретизации (Гц)
 Fc = 500e3; % частота несущей (Гц)
 A  = 511;   % средняя амплитуда входного сигнала 
 M = 32; % Разрядность фазового аккумулятора DDS 
 fres = Fs/(2^M); % Разрешение по частоте DDS (Гц) 
 Ko = fres*2*pi; % Коэфф. усиления DDS (рад/сек * разряд) 
 fprintf('Коэфф. усиления DDS %d\n', Ko);
 Kd = 1.0; % коэффициент Kd фазового детектора
 wp = 2.0*pi*200.0; % резонансная частота контура    (рад/с)
 zeta = 0.7; % демпинг фактор
 T = 1.0/Fs; % интервал дискретизации (сек)
 phi = -10.5; % начальная фазовая расстройка (рад)
 dF = 30e3; % начальная частотная расстройка (Гц)
 fg = Fc-dF; % частота ГУН (Гц)

 %расчет к-тов g1 и g2
 g1 = 2.0*(1.0 - exp(-wp * zeta * T) *cos(wp * sqrt(1.0 - zeta * zeta) * T));
 g2 = exp(-2.0 * wp * zeta * T)-1.0 + g1;

 %расчет к-тов Ki и Kp петлевого фильтра
 ki = g2/(Ko * Kd * A);
 kp = g1/(Ko * Kd * A);
fprintf('Пропорц. коэфф Kp %d\n', kp);
fprintf('Интегр. коэфф Ki %d\n', ki);
Получаю:

Код: Выделить всё

Коэфф. усиления DDS 4.494084e-02
Пропорц. коэфф Kp 2.493761e-06
Интегр. коэфф Ki 7.286227e-11
Реально моделька работает с коэфициентами на 7*10^6 ,больше чем расчетные.
В чем ошибка ?

Суть модельки: подстроиться по частоте на 30 КГц к несущей по частоте и фазе.
Допустим, расстройка по частоте и фазе постоянна.
При АЦП 12 бит (имитируется усилителями на модельке).
При NCO 10 бит (имитируется усилителями на модельке)
Модель simulink_2017- прикрепленный файл.

https://drive.google.com/file/d/1jBjfQ0 ... sp=sharing

https://drive.google.com/file/d/1Y10Izz ... sp=sharing

Ответить