Pulse Shaping Filter, Правильно ли понимаю? Опции

[Politeh]

Здравствуйте!

Хотелось бы узнать работу PSF(Pulse Shaping Filter).

В модеме 8PSK в моём случае используется PSF преподнятый косинус. Не могу точно понять как он работает. В книге К.Феер "Беспроводная цифровая связь" написано, что подавая импульсы на PSF преподнятого косинуса, на выходе получаем отклик - импульсную характеристику. Это понятно, но конкретно в модеме тоже нужно модулировать несущую такими же откликами(импульсной характеристикой), или просто обычные импульсы(NRZ) пропускаем через этот PSF, а потом этим резльтатом модулируем несущие в квадратурных каналах?

Спасибо.

Эскизы прикрепленных изображений

 

[Politeh]

В матлабе поэкспериментировал с функциями. Вот что получилось. Верхний рисунок - импульсная характеристика. Нижний - реакция фильта при подаче на него трёх последовательных единиц (1 1 1) с символьной скоростью. Теперь я так понимаю такими импульсами и модулируются квадратурные несущие: косинус и синус. Естественно что что они должны подаваться в нужном соотношении(после последовательно параллельного преобразователя который отражает группы из трёх бит в координаты векторов). После фильтра ставится ЦАП и перемножается в аналоговом виде. Только не понимаю, если у меня подаётся два одинаковых импульса(+1) на каждый канал, то после модуляции должна быть фаза в 45 градусов, но тогда и амплитудная модуляция будет присутствовать, так как не могу же сделать я амплитуды импульсов 0.707... .

Эскизы прикрепленных изображений

 

[petrov]

Естественно что что они должны подаваться в нужном соотношении(после последовательно параллельного преобразователя который отражает группы из трёх бит в координаты векторов). После фильтра ставится ЦАП и перемножается в аналоговом виде. Только не понимаю, если у меня подаётся два одинаковых импульса(+1) на каждый канал, то после модуляции должна быть фаза в 45 градусов, но тогда и амплитудная модуляция будет присутствовать, так как не могу же сделать я амплитуды импульсов 0.707... .

А что это за координата такая (1,1), нету в 8PSK такой точки, на выходе таблицы у вас должно быть как раз (0.707, 0.707) для 45 градусов.

[Politeh]

А что это за координата такая (1,1), нету в 8PSK такой точки, на выходе таблицы у вас должно быть как раз (0.707, 0.707) для 45 градусов.

Понял. Спасибо.

В общем да, но в реальных системах, работающих в эфире или с частотным уплотнением, есть есть требование к уровню помех в соседнем канале. Например для КВ диапазона это -40дБ, насколько помню.
Чтобы это требование обеспечить, используют цифровые фильтры для модулирующих сигналов. Подавать просто единицы и нули на вход модулятора нельзя. Уровень боковых лепестков этого фильтра должен быть соответствующим.

Неправильно, на ЦАП подаются отсчеты уже прошедшие фильтрацию цифровым фильтром. В реальном времени это делать ни к чему, обычно все просчитывается заранее и записывается в ПЗУ.
За один битовый отсчет нужно выбрать из ПЗУ и подать на ЦАП число отсчетов равное порядку фильтра. Поэтому удобнее выбирать порядок как степень двух.

Почему не можете?
Все это можно учесть и выбирать для такой комбинации соответствующие части прошивок.
Здесь лежат ссылки на программы расчета цифровых фильторов.
Я из них пользовался только FilterExpress, работает.

Это я уже промоделировал в Матлабе с использованием фильтра. Правда сейчас у меня почему то отклики налазят друг на друга, странно - слишком растянутая импульсная характеристика получилась(ширина её при порядке равном 9 выходит около 0.4 мс), хотя параметры выбрал из стандарта: коэффициент сглаживания = 0.6. Частота симвовлов 10500 Гц, период соот-но 0.095 мс.

Вот код:

fs = 10500; - частота символов


Ts = 1/fs;

fd = 2*fs; - частота дискретизации

[num,den] = rcosine(fs,fd,'fir',0.6, 2); - расчет фильтра.
[H T]= impz(num); - импульсная хар-ка
T = T.*(1/fd); - ось времени
subplot(2,1,1)
stem(T, H) - рисуем
grid on;
% [A w] = freqz(num, den); - рисуем АЧХ
% subplot(2,1,2)
% w*fd;
% plot(w*fd/pi, 20*log10(abs(A)))
% ylim([-100 10])
% grid on


[y TT] = rcosflt([0.7 ,0.7],fs/6, fd,['filter', 'fir'],num); пропускаем два импульса по 0.7 через фильтр.
subplot(2,1,2)
stem(TT,y) - рисуем.

[alexkok]

Понял. Спасибо.
Это я уже промоделировал в Матлабе с использованием фильтра. Правда сейчас у меня почему то отклики налазят друг на друга, странно - слишком растянутая импульсная характеристика получилась(ширина её при порядке равном 9 выходит около 0.4 мс), хотя параметры выбрал из стандарта: коэффициент сглаживания = 0.6. Частота симвовлов 10500 Гц, период соот-но 0.095 мс.

По поводу Матлаба ничего сказать не могу, не пользуюсь.
Но частота отсчетов должна быть больше символьной в N раз, где N - порядок фильтра.

[petrov]

Но частота отсчетов должна быть больше символьной в N раз, где N - порядок фильтра.

Никак не связан порядок фильтра с частотой отсчётов.

[alexkok]

Никак не связан порядок фильтра с частотой отсчётов.

Смотря что понимать под порядком фильтра.
Для FIR фильтра порядок это число коэффициентов и частота отсчетов не может быть меньше чем символьная частота умноженная на число коэффициентов. Больше - может.

[petrov]

Смотря что понимать под порядком фильтра.
Для FIR фильтра порядок это число коэффициентов и частота отсчетов не может быть меньше чем символьная частота умноженная на число коэффициентов. Больше - может.

Ерунда какая-то...
Допустим символьная частота 1, 8 отсчётов на символ, т е частота отсчётов 8, коэффициентов 200, вот меньше запросто.

[alexkok]

Ерунда какая-то...
Допустим символьная частота 1, 8 отсчётов на символ, т е частота отсчётов 8, коэффициентов 200, вот меньше запросто.

Ну и промоделируйте какой уровень боковых лепестков будет в соседнем канале для обсуждаемого случая:
символьная частота - 10кГц, т.е. полоса канала 25кГц,
тип фильтра - приподнятый косинус,
плюс Ваши условия.
Что будет за этой полосой?
Я навскидку оцениваю в - 15..-18дБ.