[dimitriy91 0]

Доброго времени суток. Подскажите пожалуйста как рассчитать фильтры НЧ (не петлевой) петли Костаса для BPSK. Где то слышал что частота среза должна быть равна бодовой скорости, но мне кажется это бред. Бодовая скорость теоретически может быть равна двум полосам сигнала (по Найквисту). Не могу понять чему равна наивысшая частота спектра сигналов:

cos( fi(t) ) и sin(fi(t) ), где fi(t) - принимает значение 0 или 180 в бодовые интервалы, предполагаем, что fi(t) случайно и равномерно принимает 0 или 180.

 

http://www.google.ru/imgres?imgurl=https:/...ZBQQQMwghKAUwBQ

Изменено 30 ноября, 2015 пользователем abraziv

[litv 0]

не петлевой петли :laughing: .

 

Были модели by Petrov для simulink на этом форуме.

Еще было:

sampling_freq = частота дискретизации;

capture_band = полоса захвата;

AMP = средняя амплитуда входного сигнала (видимо ошибки);

koef = коэффициент, действительно принятый за 0.707;

 

Вычисляем след. коэффициенты:

T (период) = 1/sampling_freq;

wn = (capture_band*PI)/koef;

 

g1 = 1 - exp(-2.0*koef*wn*T);

g2 = 1 + exp(-2.0*koef*wn*T) - 2.0*exp(-koef*wn*T)*cos(wn*T*sqrt(1-(koef*koef)) );

 

G1 = g1/AMP;

G2 = g2/AMP;

 

Z = 0;

 

формула фильтрации

 

out = inp*G1 + (inp+Z)*G2; - выходное значение

Z+=inp; // накопительный регистр

 

 

 

Ведет себя очень здорово.

Использовал для тактовой и фазовой синхронизации:

При частоте входного сигнала 2400 Гц :

полоса захвата 1 Гц - при стабильной работе;

полоса захвата 10 - 100 Гц - когда надо быстро что-нибудь выкрутить;

 

Кстати, расчитываемые коэффициенты можно менять во время работы, чтобы регулировать скорость реакции системы.

 

[petrov 6]

...

 

Просто согласованный фильтр. У вас на рисунке прямоугольные импульсы, вот в I и Q ветках по согласованному фильтру с прямоугольным импульсом должно стоять.

[gosha 0]

On 11/30/2015 at 2:18 PM, litv said:

Были модели by Petrov для simulink на этом форуме.

sampling_freq = частота дискретизации;

capture_band = полоса захвата;

AMP = средняя амплитуда входного сигнала (видимо ошибки);

koef = коэффициент, действительно принятый за 0.707;

Вычисляем след. коэффициенты:

T (период) = 1/sampling_freq;

wn = (capture_band*PI)/koef;

g1 = 1 - exp(-2.0*koef*wn*T);

g2 = 1 + exp(-2.0*koef*wn*T) - 2.0*exp(-koef*wn*T)*cos(wn*T*sqrt(1-(koef*koef)) );

G1 = g1/AMP;

G2 = g2/AMP;

 

 Fs = 30720000; % частота дискретизации (Гц)
 f0 = 16000.0; % частота сигнала (Гц)
 kd = 1.0; % коэффициент Kd фазового детектора
 ko = 1000.0; % коэффициент Ko ГУН
 wp = 2.0*pi*16000.0; % резонансная частота контура    (рад/с)
 zeta = 0.7; % демпинг фактор
 T = 1.0/Fs; % интервал дискретизации (сек)
 phi = -1.5; % начальная фазовая расстройка (рад)
 df = 4.0; % начальная частотная расстройка (Гц)
 fg = f0-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);
 kp = g1/(ko*kd);
fprintf('Пропорц. коэфф Kp %d\n', kp);
fprintf('Интегр. коэфф Ki %d\n', ki);

Коэфициент Kd фазового детектора как определяется ?

Коэфициент Ko VCO == Input sensitivityThis value scales the input voltage and, consequently, the shift from the Quiescent frequency value. The units of Input sensitivity are Hertz per volt. ?

 

https://uk.mathworks.com/help/comm/ref/discretetimevco.html

Если задача- выделить поднесущую 16КГц, для вычисления резонансной частоты подставляем 16КГц в эту формулу? :  wp = 2.0*pi*16000.0; % резонансная частота контура    (рад/с) ?

 

pll_model.slx