Фильтрация фазового шума Опции

[masverter]

Доброго времени суток!
Столкнулся со следующей проблемой. Имеется сигнал ФМ8, отношение сигнал/шум порядка 17 дБ. На 8 гармонике в районе несущей частоты присутствуют довольно сильные помехи:
.

В итоге на приемной стороне имеем нестабильное удержание несущей:
.

Можно ли фильтровать подобные помехи и как делать это наиболее эффективным образом? Хотелось бы узнать опыт сообщества по данному вопросу.
Заранее спасибо.

P.S.: сам пробовал использовать схему из работы  cupo1989.pdf ( 1.14 мегабайт ) Кол-во скачиваний: 408
В данной работе адаптивная фильтрация фазового шума выполняется после ФАПЧ, что не позволяет устранить частые "срывы" в самой ФАПЧ.

[serebr]

Судя по спектру, в принятом сигнале нет особых помех, превышающих полезный сигнал в какой-либо полосе частот. Позже посмотрю поподробнее.
Частота предоставленных в файле I/Q отсчетов для простоты принята равной 1 МГц.



Если я правильно выцедил информацию из предоставленного принятого сигнала, то:
1. Частота I/Q пар равна 4 отсчетам на каждый передаваемый символ
2. Способ модуляции - фазовая модуляция, 8 значений фазы

Могу предположить, что
- есть какой-то фильтр на стороне передатчика для уменьшения межсимвольной интерференции.
- используется цифровой полосовой фильтр на стороне приёмника для подавления внеполосных шумов + аналогичный фильтр на стороне передатчика

Дискретная фазовая модуляция является частным случаем QAM модуляции. Поэтому никаких особых премудростей по поводу фильтрации именно фазового шума я бы не ожидал где-либо встретить. Напрашивается использование адаптивного эквалайзера на строне приемника для подавления межсимвольной интерференции и выравнивания результирующей частотной характеристики (как амплитудной, так и фазовой).

[masverter]

Судя по спектру, в принятом сигнале нет особых помех, превышающих полезный сигнал в какой-либо полосе частот.

При возведении в 8-ю степень проявляются

Код

NFFT = 32768;
NAVG = 10;
nsamples =  NFFT*NAVG;

fin = fopen('rcv_iq16');
N = nsamples * 2;
d = fread(fin, N, 'int16', 0, 'l');
re = d(1:2:N);
im = d(2:2:N);
x = complex(re,im);
fclose(fin);

y = zeros(NFFT,1);
for i=1:NAVG
    y = y + abs(fft(x((i-1)*NFFT+1:i*NFFT).^8));
end
yout = zeros(NFFT,1);
yout(1:NFFT/2) = y(NFFT/2+1:NFFT);
yout(NFFT/2+1:NFFT) = y(1:NFFT/2);
plot(yout);