Код: Выделить всё
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
В чем ошибка ?
Суть модельки: подстроиться по частоте на 30 КГц к несущей по частоте и фазе.
Допустим, расстройка по частоте и фазе постоянна.
При АЦП 12 бит (имитируется усилителями на модельке).
При NCO 10 бит (имитируется усилителями на модельке)
Модель simulink_2017- прикрепленный файл.
https://drive.google.com/file/d/1jBjfQ0 ... sp=sharing
https://drive.google.com/file/d/1Y10Izz ... sp=sharing