[Waso 1]

Имеется радиосигнал, модулированный прямоугольными импульсами переменной скважности и немного зашумленный. Сигнал подвергается 2048-точечному БПФ и по полученному спектру я должен определить текущую скважность сигнала. Если смотреть на спектр чистого импульсного сигнала, то первый минимум огибающей (который рядом с максимумом) соответствует частоте, обратной величине ширины импульсов. Я пока написал дубовый алгоритм, который перебирает кратные гармоники и приблизительно определяет этот минимум. Но когда в сигнал добавляется шумок, спектр начинает дрожать. Основная частота определяется четко, а вот по скважности моя програмулина начинает давать шнягу...

 

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

[AlexeyW 0]

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

Вообще, в зашумленных задачах хороши интегральные методы. Например, я как-то использовал цифровой осциллограф для измерения АЧХ. Сначала примитивно - прога считывает с осциллографа цифровые осциллограммы и находит максимум. Но намного точнее и устойчивее было, когда производилось интегрирование квадрата функции и т.д. (не забыв о сдвиге нуля).

 

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

Огибающая ожидаемой функции плавает вместе со скважностью, являясь ее функцией. Соответственно, при корелляционном методе нужно искать максимум корелляции F(экспериментальной) с F(теоретической как функции от скважности).

[Waso 1]

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

 

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

[тау 29]

Очень медленно получается. Или есть способ попроще? Я просто не очень знаком с корелляционными методами...

тогда просто проинтегрируйте огибающую перед фурье - шум уменьшится.

[AlexeyW 0]

Кстати, а для чего Вам первоначальное Фурье-преобразование? Почему не хотите использоватьпросто сигнал, и работать с ним без особых сложностей?

Изменено 14 декабря, 2010 пользователем AlexeyW

[Waso 1]

тогда просто проинтегрируйте огибающую перед фурье - шум уменьшится.

Это интересный вариант! Спасибо.

Кстати, а для чего Вам первоначальное Фурье-преобразование? Почему не хотите использоватьпросто сигнал, и работать с ним без особых сложностей?

Ну оно так по ТЗ надо. Да и по другому легче не получится. Потому что сигнал плавает по амплитуде, шумы иногда бывают и посильнее... Так что тут как-то просто так фронты не выцепишь... Вернее выцепишь много лишнего. А так - знаешь примерно свою частоту. Ищешь в спектре там палку, от нее определяешь первый минимум и готово. Научиться чистить спектр и будет счастье!

[AlexeyW 0]

тогда просто проинтегрируйте огибающую перед фурье - шум уменьшится.

А я, честно признаться, не совсем понял. До преобразования? А как именно?

[Гость Савинова Анна]

Не поняла: так есть только спектр и приблизительный вид сигнала?

[Fast 0]

А я, честно признаться, не совсем понял. До преобразования? А как именно?

суммируем на интервале T

 

а можно проинтегрировать по спектральным компонентам после FFT, усреднив несколько раз мгновенный спектр, я б так сделал

 

 

[Waso 1]

Не поняла: так есть только спектр и приблизительный вид сигнала?

Вообще нужно анализировать только спектр. Но в принципе доступ к данным АЦП тоже есть.

 

[тау 29]

А я, честно признаться, не совсем понял. До преобразования? А как именно?

да это я чушь предложил, не поможет.

 

 

Вообще нужно анализировать только спектр.

поскольку гармоники спектра соответствуют 2E/(pi*n)*sin(pi*n * t/T) , то организовав цикл по длительности (и по амплитуде Е если она не известна) предполагаемого импульса t в разумных пределах, ищете сумму среднеквадратических разностей между квадратами значений спектра идущих с интервалом 1/T и квадратом по формуле . Определяете наилучшее совпадение по минимуму суммы для всех итераций. Чтобы промежуточные бины спектра не мешались - в этих точках вычислений не производите. Можно к примеру отсеять по амплитуде >0.1 макс.

T (период) определите усреднением по "высоким" бинам , а длительность t по наилучшему совпадению . и затем скважность.

[AlexeyW 0]

суммируем на интервале T

увеличит такой метод шумовую погрешность или уменьшит (по сравнению с прямым обнаружением фронта и спада для измерения длительности), напрямую зависит от длительности фронта "прямоугольного импульса".

 

то организовав цикл по длительности (и по амплитуде Е если она не известна)

не исключено, что можно обойтись и без цикла. Есть же такое понятие - поиск экстремали функционала. Кажется, тут именно то - поиск значения параметра, от которого зависит корелляционная функция, при котором она имеет экстремум.

[Andrey_1 0]

увеличит такой метод шумовую погрешность или уменьшит (по сравнению с прямым обнаружением фронта и спада для измерения длительности), напрямую зависит от длительности фронта "прямоугольного импульса".

 

 

не исключено, что можно обойтись и без цикла. Есть же такое понятие - поиск экстремали функционала. Кажется, тут именно то - поиск значения параметра, от которого зависит корелляционная функция, при котором она имеет экстремум.

 

Измерьте спектр шума в отсутствии сигнала и тупо вычтете из спектра сигнал+шум

[AlexeyW 0]

Измерьте спектр шума в отсутствии сигнала и тупо вычтете из спектра сигнал+шум

Если б все было так просто, мы бы шум вообще победили :)

[Andrey_1 0]

А я, честно признаться, не совсем понял. До преобразования? А как именно?

 

ФНЧ