Выделение сигнала на фоне помех - Форум разработчиков электроники ELECTRONIX.ru

Выделение сигнала на фоне помех

Марк_Я Просмотр профиля	Oct 28 2009, 13:27 Сообщение #1
Частый гость Группа: Свой Сообщений: 178 Регистрация: 19-09-07 Из: М.О. г.Фрязино Пользователь №: 30 656	Господа! Рискну сформулировать вопрос в основной теме, поскольку начинающие могут его проигнорировать.. Имеем периодический импульсный сигнал. Медленно меняющий форму и расположение самого импульса в периоде, но с очень стабильным периодом. Назовем семплом оцифрованный один период. Всего семплов анализируем восемь штук. Помеха выглядит как хаотично появляющийся в любом месте периода импульс с формой подобной сигналу. Анализ провожу, дифференцируя семпл на отсчеты АЦП, и рассматривая одноименные отсчеты АЦП всех 8 семплов сигнала. Далее строю гистограмму по каждой точке внутри семпла и нахожу весовой центр распределения этой гистограммы. Извините за вольную терминологию... Вопрос в следующем. Я никак не могу найти или придумать компактный алгоритм нахождения этого весового центра. То, что в лоб имеет следующий алгоритм: Берем нулевой массив из 256 ячеек для 8-разрядного преобразования АЦП (у меня такое). Инкрементируем элемент массива с адресом равным отсчету АЦП. Получаем гистограмму этих 8 одноименных (в периоде) отсчетов. Затем двигаемся воротами анализа от 0 до 255 и, как только в воротах анализа окажется не менее, скажем 5 единиц (сумма значений элементов массива в воротах анализа), так центр ворот и будет весовым центром гистограммы. Если нигде не будет обнаружено пика, значит помеха слишком велика и детектирования сигнала не произойдет. Поругайте мой "деревянный велосипед с квадратными колесами"и направьте на путь истинный... Спасибо.

Ответов

Марк_Я Просмотр профиля	Nov 1 2009, 09:14 Сообщение #2
Частый гость Группа: Свой Сообщений: 178 Регистрация: 19-09-07 Из: М.О. г.Фрязино Пользователь №: 30 656	Специально для алекскока. Имеем три резонансных контура со сдвигом перекрывающим девиацию возбуждающих их ОТДЕЛЬНЫХ генераторов с ЧМ модуляцией. Все три генератора работают на одних частотах, но частота модуляции слегка у всех расходится. Принимает сумму всех трех сигналов на одну магнитную антенну. Получаем сплошную кашу. Теперь подаем на приемник синхронизирующий модулирующий ОДИН из генераторов сигнал. Получаем синхронизацию только одного из трех откликов. Теперь давайте, умненький Вы наш, коррелируйте... Языком молоть я тоже умею, и при непреодолимом желании (дождетесь...) могу обложить вежливо и непринужденно - мама не балуй... Еще раз для "жирафов" обожающих все вокруг себя коррелировать: весь фрейм забит помехами. контура могут иметь произвольную внутри полосы анализа частоту. калибрующих интервалов нет. добротности контуров произвольны. Вопрос конкретный: чего с чем сворачивать для корреляции? Пока тут языком "корреляторы" мололи впустую, я спробовал один вариант стат обработки. Поделюсь "назло врагам" с друзьями. Мож оценят-помогут-скорректируют. Итак. Каждый массив одноименных выборок (помним, у нас их ровно 8) усредняем (сиречь тупо складываем и сдвигаем результат вправо на три разряда). Далее пробегаем по всем восьми выборкам и сравниваем с полученным средним. Ежели отклонение от среднего выходит за некоторый порог (находим опытным путем), то заменяем такие выборки средним. Кроме того, во время такого "причесывания" считаем количество скорректированных отсчетов. Если таких "уклонистов" более 2-х(или 3-х, зависит от интенсивности помех, определяем эмпирически) - метим всю эту "братию" из 8 отсчетов как пораженную точку, в противном случае находим новое среднее... Все... После окончательного расчета (включая квадратурное детектирование амплитуды) получаем результирующий массив и вектор пораженных точек. Можно эти точки при желании экстраполировать, сохраняя три последних результирующих фрейма... Черт возьми... работает и в скользящем окне... и на 70-рублевом камне... Тем не менее жду отзывов и иных вариантов... Спасибо... Медиана пригодилась... Thanks... ЗЫ. Насчет "оторвать яйца". Ежели такой студент-"коррелятор" ко мне на практику придет - всенепременно причинное место отсеку, аки барану... ЗЗЫ. DRUID3, это круто! Я не великий спец в мат.методах, но случайно попал точно в метод по Вашей ссылке. Только среднеквадратичное отклонение вычислять нет возможности (не поспеваю извлекать корень), потому довольствовался просто средним отклонением... Сообщение отредактировал Марк_Я - Nov 1 2009, 09:30

sup-sup Просмотр профиля	Nov 1 2009, 09:36 Сообщение #3
Знающий Группа: Участник Сообщений: 674 Регистрация: 26-08-05 Пользователь №: 7 997	Цитата(Марк_Я @ Nov 1 2009, 13:14) Каждый массив одноименных выборок (помним, у нас их ровно 8) усредняем (сиречь тупо складываем и сдвигаем результат вправо на три разряда). Далее пробегаем по всем восьми выборкам и сравниваем с полученным средним. Ежели отклонение от среднего выходит за некоторый порог (находим опытным путем), то заменяем такие выборки средним. Кроме того, во время такого "причесывания" считаем количество скорректированных отсчетов. Если таких "уклонистов" более 2-х(или 3-х, зависит от интенсивности помех, определяем эмпирически) - метим всю эту "братию" из 8 отсчетов как пораженную точку, в противном случае находим новое среднее... Все... После окончательного расчета (включая квадратурное детектирование амплитуды) получаем результирующий массив и вектор пораженных точек. Можно эти точки при желании экстраполировать, сохраняя три последних результирующих фрейма... Черт возьми... работает и в скользящем окне... и на 70-рублевом камне... Тем не менее жду отзывов и иных вариантов... Спасибо... Медиана пригодилась... Thanks... ЗЫ. Насчет "оторвать яйца". Ежели такой студент-"коррелятор" ко мне на практику придет - всенепременно причинное место отсеку, аки барану... ... прибежал, взволнован крайне, сообщеньем нас потряс ... ( :-) - уравновешивающий эпиграф). Еще есть анекдот про верблюда в пустыне ("подержи верблюда" - не про автора темы, а просто ассоциации. И еще, в книжке про то как физики шутят есть определение "точки зрения" Ж-) В любом случае, если дело движется - это хорошо. Сообщение отредактировал sup-sup - Nov 1 2009, 09:55

Сообщений в этой теме

Марк_Я Выделение сигнала на фоне помех Oct 28 2009, 13:27

SFx сам я тоже особо не в теме, но в теории вроде для ... Oct 28 2009, 17:01

Марк_Я Это не символы. Определять начало не требуется, по... Oct 28 2009, 19:24

alexkok Цитата(Марк_Я @ Oct 28 2009, 17:27) Берем... Oct 28 2009, 20:20

Марк_Я Не понял. Я ничего не селектирую. Я пытаюсь найти ... Oct 28 2009, 20:57

alexkok Цитата(Марк_Я @ Oct 28 2009, 23:57) Не по... Oct 29 2009, 05:52

rezident Раз период стабильный и фаза импульса медленно изм... Oct 28 2009, 21:40

Марк_Я Я специально выделил в названии слово ПОМЕХА. Зада... Oct 29 2009, 04:18

des00 Цитата(Марк_Я @ Oct 28 2009, 22:18) Я спе... Oct 29 2009, 06:17

sup-sup Допустим, отбрасывать по два максимальных и миним... Oct 29 2009, 07:11

DRUID3 Цитата(Марк_Я @ Oct 29 2009, 06:18) Я спе... Oct 30 2009, 03:33

bahurin Цитата(DRUID3 @ Oct 30 2009, 06:33) Т.е. ... Oct 30 2009, 05:25

Марк_Я Ни максимальные ни минимальные значения не являютс... Oct 29 2009, 09:09

sup-sup Цитата(Марк_Я @ Oct 29 2009, 12:09) Ни ма... Oct 29 2009, 17:40

анатолий Это типичная задача когерентного накопления или ав... Oct 29 2009, 11:52

Марк_Я По слогам: ТАМ ШУМА НЕТ... Совсем... Как класса.... Oct 29 2009, 12:36

mdmitry На основе фазовой автоподстройки: подстраивается п... Oct 29 2009, 16:43

SFx почитал топик, может вам нейронную сеть обучить? Oct 29 2009, 19:32

rezident Что-то я не понял, почему корреляционный способ не... Oct 29 2009, 19:34

AndeyP Цитата(Марк_Я @ Oct 28 2009, 17:27) Берем... Oct 29 2009, 19:37

Марк_Я Цитата(DRUID3 @ Oct 30 2009, 06:33) Есть ... Oct 30 2009, 09:10

анатолий Не надо нервничать. Товарищ предлагает делать корр... Oct 30 2009, 09:23

fontp Цитата(анатолий @ Oct 30 2009, 12:23) Не ... Oct 30 2009, 09:42

Марк_Я Цитата(анатолий @ Oct 30 2009, 12:23) Не ... Oct 30 2009, 10:26

fontp Цитата(Марк_Я @ Oct 30 2009, 13:26) Госпо... Oct 30 2009, 10:32

Марк_Я Цитата(fontp @ Oct 30 2009, 12:42) Да-да-... Oct 30 2009, 10:34

fontp Цитата(Марк_Я @ Oct 30 2009, 13:34) Они в... Oct 30 2009, 10:36

Марк_Я Цитата(fontp @ Oct 30 2009, 13:32) Хорошо... Oct 30 2009, 10:46

fontp Цитата(Марк_Я @ Oct 30 2009, 13:46) Гисто... Oct 30 2009, 10:50

Марк_Я Вооо! Это уже по теме... Смысл моего вопроса в... Oct 30 2009, 11:36

анатолий Цитата(Марк_Я @ Oct 30 2009, 14:36) Вооо... Oct 30 2009, 12:12

fontp Если это всё скользит во времени, скользящее окно ... Oct 30 2009, 11:54

Марк_Я Ну скользящее окно на последние 2х8х64 это просто ... Oct 30 2009, 12:07

fontp Цитата(Марк_Я @ Oct 30 2009, 15:07) Ну ск... Oct 30 2009, 12:13

Марк_Я "аппаратная циклическая" есть - это назы... Oct 30 2009, 12:29

fontp Цитата(Марк_Я @ Oct 30 2009, 15:29) ... Oct 30 2009, 12:32

Марк_Я Цитата(fontp @ Oct 30 2009, 15:13) Не-е з... Oct 30 2009, 12:36

fontp Цитата(Марк_Я @ Oct 30 2009, 15:36) Трудн... Oct 30 2009, 12:39

Марк_Я Так думаю... всю голову уже сломал. Вот сюда обра... Oct 30 2009, 12:45

pokos Чото я не пойму метаний. Как уже отметил анатолий,... Oct 30 2009, 13:10

Марк_Я Специально для отстающих... Сов.секретно. Есть д... Oct 30 2009, 13:24

sup-sup Цитата(Марк_Я @ Oct 30 2009, 17:24) Специ... Oct 30 2009, 16:03

Марк_Я Даже и не знаю что сказать... Если коррелятор мож... Oct 30 2009, 18:17

sup-sup Цитата(Марк_Я @ Oct 30 2009, 22:17) Даже ... Oct 30 2009, 22:12

alexkok Цитата(sup-sup @ Oct 31 2009, 02:12)... Oct 30 2009, 22:26

sup-sup Цитата(alexkok @ Oct 31 2009, 01:26) Я эт... Oct 30 2009, 22:56

Марк_Я Господа! Я уже практически рассказал больше, ч... Oct 31 2009, 08:43

alexkok Цитата(Марк_Я @ Oct 31 2009, 12:43) Госпо... Nov 1 2009, 03:55

fontp Я вот чтот подумал: А зачем Вам гистограмма, Марк?... Oct 31 2009, 11:40

alexkok Цитата(fontp @ Oct 31 2009, 14:40) Ну не ... Nov 1 2009, 00:58

fontp Цитата(alexkok @ Nov 1 2009, 03:58) Приме... Nov 1 2009, 06:49

DRUID3 Цитата(alexkok @ Nov 1 2009, 02:58) Приме... Nov 1 2009, 06:56

Марк_Я Специально для алекскока. Имеем три резонансных ко... Nov 1 2009, 09:14

sup-sup Цитата(Марк_Я @ Nov 1 2009, 13:14) Каждый... Nov 1 2009, 09:36

alexkok Цитата(Марк_Я @ Nov 1 2009, 12:14) Специа... Nov 1 2009, 19:50

анатолий Цитата(Марк_Я @ Nov 1 2009, 12:14) Еще ра... Nov 7 2009, 15:34

Марк_Я Цитата(sup-sup @ Nov 1 2009, 12:36) ... Nov 1 2009, 11:10

sup-sup Цитата(Марк_Я @ Nov 1 2009, 15:10) Снисхо... Nov 1 2009, 12:20

Марк_Я Эка Вы любезный куда загнули? Преподленькая тракт... Nov 1 2009, 13:00

sup-sup Цитата(Марк_Я @ Nov 1 2009, 17:00) Эка Вы... Nov 1 2009, 14:50

Марк_Я Что, обидно? В ноги не упал... не бился в конвуль... Nov 1 2009, 15:23

sup-sup Цитата(Марк_Я @ Nov 1 2009, 18:23) Что, о... Nov 1 2009, 16:59

Марк_Я За два часа ни одного. Больше не смотрел. Это при ... Nov 1 2009, 19:52

alexkok Цитата(Марк_Я @ Nov 1 2009, 23:52) Что ка... Nov 4 2009, 21:25

Марк_Я Увы, нет. Это совсем не RFID... Даже не тепло... Nov 6 2009, 06:47

Марк_Я Еще один популяризатор... Вы хоть цитату прочтите ... Nov 8 2009, 03:48

анатолий Цитата(Марк_Я @ Nov 8 2009, 06:48) Чем ан... Nov 10 2009, 13:04

Марк_Я 1. С конца. Ни на что не претендую и ни на кого не... Nov 10 2009, 14:15

анатолий В этом случае линейная теория учит, что: - надо ра... Nov 10 2009, 16:57

fontp Цитата(анатолий @ Nov 10 2009, 19:57) В э... Nov 10 2009, 17:40

Марк_Я Цитата(анатолий @ Nov 10 2009, 19:57) В э... Nov 10 2009, 19:18

alexkok Цитата(Марк_Я @ Nov 10 2009, 23:18) - Рас... Nov 10 2009, 20:08

Марк_Я Цитата(alexkok @ Nov 10 2009, 23:08) А от... Nov 11 2009, 04:20

« Предыдущая тема · Алгоритмы ЦОС (DSP) · Следующая тема »

1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)

Пользователей: 0