Демодуляция АМ сигнала Опции

[Allregia]

Есть сигнал, несущая около 1кгц, промодулированная (АМ) низкочастотным (0.1-10гц) сигналом.
Сигнал приходит (с АЦП) в виде uint16 /10kHz SR.
Нужно выделить модулирующий сигнал и сделать децимацию на 50гц.

Сейчас я делаю просто: output=IIR_LPF_20Гц (abs(input));
Затем беру каждый 20-й отсчет.

Все работает, но - так как индекс модуляции весьма мал (3-5%), то и сигнал на выходе я получаю весьма маленькой амплитуды (разрядности). Вместе с тем, когда я беру модуль входного сигнала, он содержит большую постоянную составляющую. Если перед тем как делать LPF вычесть константу из этого модуля, то наверное можно повысить разрядность выходного сигнала.
Вопрос первый - так ли это (типа "а не дурк ли я?" ), и второй - как лучше выделить эту константу (постоянную составляющую)? Амплитуда несущей может быть разной, это непредсказуемо.
Выделять по простому, поставить еще один LPF, с частотой среза ниже нижней частоты полезного сигнала, и вычитать его выход их модуля входного сигнала, перед основым LPF и децимацией?

P.S. Все происходит не на ДСП а на обычном МК (АРМ7).

[fontp]

Цифровой детекторный приемник? Круто, мечта начинающего радиолюбителя

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

[Allregia]

Цифровой детекторный приемник?

Ага

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

Несущая не фиксирована и не стабильна.

[sup-sup]

Несущая не фиксирована и не стабильна.

Все-таки, какая-то полоса есть? Если это так, то наверняка, полезен будет полосовой фильтр для отделения сигнала от помех. Таким образом, уйдет постоянная составляющая. А разрядность 'сама повысится' на выходе фильтра.
Можно сразу решить задачу с помощью квадратурного приема огибающей. Для вычисления каждой точки огибающей нужно выполнить две свертки (сразу с нужным прореживанием, т.е. один раз в 20 мс). Делать свертку сигнала (перемноженного на окно эквивалентного фнч, выбранное в зависимости от болтанки несущей + полоса сигнала) с синусом и косинусом центральной частоты. Для получения точки огибающей нужно векторно сложить результаты сверток (сумма квадратов под корнем). Чем узкополоснее будет фильтр (длиннее ядро), тем выше разрядность выходного сигнала. Останется только отбросить младшие разряды.
PS. Скорее всего, нужен будет еще один полосовой фильтр уже для огибающей

Сообщение отредактировал sup-sup - Sep 1 2011, 16:29

[Allregia]

Все-таки, какая-то полоса есть? Если это так, то наверняка, полезен будет полосовой фильтр для отделения сигнала от помех.

Полоса есть конечно, примерно 800-2500 гц. Полосовой фильтр есть - после АЦП стоит HPF 4 порядка на 600гц, без него там вообще все в НЧ гумах и наводках потонуло бы.

Таким образом, уйдет постоянная составляющая. А разрядность 'сама повысится' на выходе фильтра.

В исходном сигнале постоянной составляющей нет - там сигнал через несколко конденсаторов еще до АЦП проходит. Она появляется когда я модуль беру (abs)

Можно сразу решить задачу с помощью квадратурного приема огибающей. Для вычисления каждой точки огибающей нужно выполнить две свертки (сразу с нужным прореживанием, т.е. один раз в 20 мс). Делать свертку сигнала (перемноженного на окно эквивалентного фнч, выбранное в зависимости от болтанки несущей + полоса сигнала) с синусом и косинусом центральной частоты. Для получения точки огибающей нужно векторно сложить результаты сверток (сумма квадратов под корнем). Чем узкополоснее будет фильтр (длиннее ядро), тем выше разрядность выходного сигнала. Останется только отбросить младшие разряды.
PS. Скорее всего, нужен будет еще один полосовой фильтр уже для огибающей

Врядли успею в реалтайме, но попробую.

[bahurin]

Сначала вы говорите

Есть сигнал, несущая около 1кгц, промодулированная (АМ) низкочастотным (0.1-10гц) сигналом.
Сигнал приходит (с АЦП) в виде uint16 /10kHz SR.
Нужно выделить модулирующий сигнал и сделать децимацию на 50гц.

А потом

Полоса есть конечно, примерно 800-2500 гц. Полосовой фильтр есть - после АЦП стоит HPF 4 порядка на 600гц, без него там вообще все в НЧ гумах и наводках потонуло бы.

Я себе очень слабо представляю как модулируя АМ 10 -ти герцовым сигналом можно добиться полосы 2500 Гц да еще и на несущей 1 кГц.

Врядли успею в реалтайме, но попробую.

Да бросьте 10 кГц SampleRate можно даже в столбик на листочке успеть в реал тайме обработать. 7 arm схавает и даже не заметит.

Сообщение отредактировал bahurin - Sep 2 2011, 04:20

[GetSmart]

Сначала вы говорите
...
А потом
...
Я себе очень слабо представляю как модулируя АМ 10 -ти герцовым сигналом можно добиться полосы 2500 Гц да еще и на несущей 1 кГц.

А в промежутке было

Несущая не фиксирована и не стабильна" ... примерно 800-2500 гц.

Разве не очевидно?

10 кГц SampleRate можно даже в столбик на листочке успеть в реал тайме обработать

Это смотря какой алгоритм, сколько у него будет мат.операций на 1 сэмпл.

Сообщение отредактировал GetSmart - Sep 2 2011, 04:28

[bahurin]

А в промежутке было


Разве не очевидно?



Это смотря какой алгоритм, сколько у него будет мат.операций на 1 сэмпл.

Не передергиваете сказано было что около 1 кГц. 2.5 кГц это уже нифига не около согласитесь. Это первое. Во вторых полоса 2500 Гц и нестабильность несущей в пределах 2500 Гц это совершенно 2 разные вещи. Если несущая частота меняется от 800 до 2500 Гц, то надо знать на каком временном интервале это может происходить. Потому что изменение несущей во времени приводит к паразитной частотной модуляции на выходе, которая в свою очередь сопровождается паразитной АМ, которая в свою очередь изничтожит полезную АМ с глубиной 3-5%. Постановка задачи совершенно не понятна. Поэтому я и спросил.

[Allregia]

Не передергиваете сказано было что около 1 кГц. 2.5 кГц это уже нифига не около согласитесь.

В среднем - 1кгц.

Это первое. Во вторых полоса 2500 Гц и нестабильность несущей в пределах 2500 Гц это совершенно 2 разные вещи. Если несущая частота меняется от 800 до 2500 Гц, то надо знать на каком временном интервале это может происходить.

На любом, несущая - псевдошумоподобный сигнал. На 3гц она обрезана сверху антиалиасным фильтром на входе АЦП.

АЦП - 16 бит.

Потому что изменение несущей во времени приводит к паразитной частотной модуляции на выходе, которая в свою очередь сопровождается паразитной АМ, которая в свою очередь изничтожит полезную АМ с глубиной 3-5%. Постановка задачи совершенно не понятна. Поэтому я и спросил.

Нет там никакой ПАМ, это не радиоприемник с контуром на входе. Во всяком случае - я ее не вижу.

Да бросьте 10 кГц SampleRate можно даже в столбик на листочке успеть в реал тайме обработать. 7 arm схавает и даже не заметит.

Их таких 8 каналов, имеющийся алголритм успевает, но загрузка 72мгц проца уже под 80%. т.к.. один из фильтров (LPF) приходится делать во float.

[fontp]

Их таких 8 каналов, имеющийся алголритм успевает, но загрузка 72мгц проца уже под 80%. т.к.. один из фильтров (LPF) приходится делать во float.

В любом случае оптимальное решение - синхронный детектор с ФАПЧ.
Вычислительная сложность его совсем не велика - табличный ГУН, умножитель, пропорционально-интегрирующий фильтр и IIR-ФНЧ на выходе
Другое дело, что ФАПЧ нужно проектировать - так есть литература.

Детекторный приемник - решение далеко не оптимальное. Да и у Вас возможности улучшений почти отсутствуют . Что остается - сделать полосовой фильтр на входе до взятия абсолюта, убирать постоянную составляющую на выходе узкополосным режектором, пожалуй что всё...

[Allregia]

В любом случае оптимальное решение - синхронный детектор с ФАПЧ.

Я не знаю как себя поведет PLL на шумоподобном сигнале.

Детекторный приемник - решение далеко не оптимальное. Да и у Вас возможности улучшений почти отсутствуют . Что остается - сделать полосовой фильтр на входе до взятия абсолюта, убирать постоянную составляющую на выходе узкополосным режектором, пожалуй что всё...

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