многокаскадная децимация, на осциллографе - алиазинг, в FFT - все нормально Опции

[beaRTS]

Подскажите, пожалуйста, что за глюк такой.
Произвожу децимацию в три этапа с 156250 S/s до 625 S/s (10, 5, 5 раз). Все параметры по граничным частотам фильтров рассчитал в приложенном МатКадовском файле ("Децимация_требования"), коэффициенты посчитаны в FDATOOL и легким движением руки вмонтированы в модель симулинка.
Все ничего, но в последнем блоке осциллографа наблюдается алиазинг, хотя спектроанализатор показывает всего две составляющих в спектре 90 и 150 Гц (как и должно быть) . В чем тут дело?

Исходный сигнал 90+150 Гц.

Буду благодарен за помощь!!!

Забыл сказать. Перед запуском на прогон Simulink модели необходимо запустить m-файл "ILSsignal"

Сообщение отредактировал beaRTS - May 21 2012, 11:18

Прикрепленные файлы

 Decimate.7z ( 42.5 килобайт ) Кол-во скачиваний: 35

 

[beaRTS]

а кто-нибудь строил цифровой SQRT-детектор для детектирования АМ сигнала???? Поделитесь плюсами/минусами реализации?? и почему во многих учебниках не говорится о нем?? в основном все статьи и учебники пестрят схемами синхронных детекторов, квадратурных преобразователей.

Сообщение отредактировал beaRTS - Jun 2 2012, 06:51

[sup-sup]

а кто-нибудь строил цифровой SQRT-детектор для детектирования АМ сигнала???? Поделитесь плюсами/минусами реализации?? и почему во многих учебниках не говорится о нем?? в основном все статьи и учебники пестрят схемами синхронных детекторов, квадратурных преобразователей.

А что это за детектор? Поискал - не нашел.
А в общем, если сигнал оцифрован, то вариантов немного. По-моему, даже один - получить огибающую.
И еще, почему-то получается так, что детектирование совмещено с фильтрацией, так как АМ обычно узкополосный.
Можно поток от АЦП (или откуда он берется) умножить на синус и косинус несущей частоты, получить IQ, затем отфильтровать их двумя одинаковыми фильтрами НЧ, а потом сделать сложение и извлечение корня. Фильтры НЧ могут быть как КИХ, так и БИХ. КИХ вроде невыгодно, но нужно учитывать, что будет децимация. Аналогичный вариант получится, если взять сразу комплексное ядро полосового КИХ (можно сказать, что там заранее ИХ ФНЧ умножили на синус и косинус) - так будет меньше вычислений.

Сообщение отредактировал sup-sup - Jun 2 2012, 07:30

[beaRTS]

А в общем, если сигнал оцифрован, то вариантов немного. По-моему, даже один - получить огибающую.

Собственно, я сейчас и кручу-верчу синхронный детектор (раскладываю на синфазный и квадратурный канал, затем фильтрую КИХ фильтром, прореживаю и через sqrt соединяю два канала в один, надо не забыть еще про цифровую фапч, которая должна сигналы синуса и косинуса сделать когерентными входной несущей). Ких фильтр хорош тем, что хоть в нем и больше коэффициентов, чем в БИХ, но фаза у него линейна (это, во-первых). А во-вторых, необязательно ж заставлять считать выход фильтра каждый такт частоты дискретизации, который равен свертке входного сигнала с импульсной характеристикой фильтра, а потом совершать прореживание. Можно же, если фильтр КИХ (и каждый его следующий выход не зависит от предыдущего как у БИХ), копить М-отсчетов во входном буфере размером М (сдвигая каждый такт), а на следующем такте запускать на счет фильтр, дальше снова копим М сэмплов, и на следующий такт запускаем на счет фильтр с этим M-буфером. В итоге выходные отсчеты уже будут прорежены в М раз, поэтому предпочтительны КИХ фильтры.

А что это за детектор? Поискал - не нашел.

У меня в прикрепленном ПДФнике описано, что в самом простом случае наш вещественный сигнал можно представить в виде:
s(t) = a(t)cos(Wot), a(t) - огибающая
тогда представим, что s(t) уже оцифрован.. и теперь с этими отсчетами надо что-то сделать, чтоб выделить огибающую... а вот теперь возьмите выражение для s(t) в квадрат.... =) что получите ??? потом возьмите ФНЧ - что останется?? Затем возьмите sqrt и будет хорошо. сегодня на работе накидал модель в симулинке - работает. Правда, не пробовал данный алгоритм со входным сигналом разного уровня, с широкополосной помехой, с узкополосной помехой.. это в будущем, т.к. только вчера математикой побаловался, обнаружил что вроде так можно делать, а потом в инете накопал ооооочень мало информации по этому вопросу. буду на работе поделюсь ссылками

[petrov]

Собственно, я сейчас и кручу-верчу синхронный детектор (раскладываю на синфазный и квадратурный канал, затем фильтрую КИХ фильтром, прореживаю и через sqrt соединяю два канала в один, надо не забыть еще про цифровую фапч, которая должна сигналы синуса и косинуса сделать когерентными входной несущей).

Это не синхронный детектор, в синхронном детекторе сигнал поворачивается чтобы занимать один из квадратурных каналов, другой не содержащий сигнала отбрасывется с половиной шума, таким образом получаем выигрыш 3 дБ.