Алгоритм Герцеля при обнаружении короткого импульса
Re: Гитберт при обнаружении короткого импульса
у меня главная проблема в Fs/N
есть неясность по поводу сигнал.шум . Попробую обьяснить
Допустим есть RF импульс 1мк , его полоса ~1мгц
Если взять окно длинной 1мкс получим при частоте 25мпс N=25
и полосу бина 1мгц
Вроде все замечательно, вкладываемся как раз в спектр(~1/T=1мкс)
Теперь вопрос если при том же входном сигнале я увеличиваю окно и соответсвенно N в 10 раз
Я получаю N=100 и ширину бина уже 100кгц, при этом сигнал будет физически воздействовать только 1 мкс из 10.
По факту я получу десять 100кгц бинов с нашим спектром.
Теперь самое главное что интересует, если сравнить сигнал _шум для нашего 1мгц бина и для суммы из десяти 100кгц бинов то сигнал_шум будет одинаковый в обоих случаях?
Поясняю зачем это мне надо. Кроме основного пульса в этой полосе возможно наличие постоянной несущей в виде помехи. Если я буду использовать 1мкс оптимальное окно(N=25) то я не получу разрешения по частоте меньше 1мгц
с другой стороны если я возьму 10 бинов и в одном из них будет помеха например на 300 кгц выше ценра я вырежу этот бин и оставлю 9 без помехи .
Просто не вижу как можно по другому убрать помеху с одной стороны и сохранить чувствительность в данной полосе с другой
есть неясность по поводу сигнал.шум . Попробую обьяснить
Допустим есть RF импульс 1мк , его полоса ~1мгц
Если взять окно длинной 1мкс получим при частоте 25мпс N=25
и полосу бина 1мгц
Вроде все замечательно, вкладываемся как раз в спектр(~1/T=1мкс)
Теперь вопрос если при том же входном сигнале я увеличиваю окно и соответсвенно N в 10 раз
Я получаю N=100 и ширину бина уже 100кгц, при этом сигнал будет физически воздействовать только 1 мкс из 10.
По факту я получу десять 100кгц бинов с нашим спектром.
Теперь самое главное что интересует, если сравнить сигнал _шум для нашего 1мгц бина и для суммы из десяти 100кгц бинов то сигнал_шум будет одинаковый в обоих случаях?
Поясняю зачем это мне надо. Кроме основного пульса в этой полосе возможно наличие постоянной несущей в виде помехи. Если я буду использовать 1мкс оптимальное окно(N=25) то я не получу разрешения по частоте меньше 1мгц
с другой стороны если я возьму 10 бинов и в одном из них будет помеха например на 300 кгц выше ценра я вырежу этот бин и оставлю 9 без помехи .
Просто не вижу как можно по другому убрать помеху с одной стороны и сохранить чувствительность в данной полосе с другой
- Бахурин Сергей
- Администратор
- Сообщения: 1114
- Зарегистрирован: 05 окт 2010, 19:55
- Контактная информация:
Re: Гитберт при обнаружении короткого импульса
ну да так можно пытаться бороться с помехой. Только надо будет потом просуммировать остальные 9 бинов, потому что при увеличении N бинов станет 10 но каждый будет ниже по уровню
Re: Гитберт при обнаружении короткого импульса
у меня главная проблема в Fs/N
есть неясность по поводу сигнал.шум . Попробую обьяснить
Допустим есть RF импульс 1мк , его полоса ~1мгц
Если взять окно длинной 1мкс получим при частоте 25мпс N=25
и полосу бина 1мгц
Вроде все замечательно, вкладываемся как раз в спектр
Теперь вопрос если при том же входном сигнале я увеличиваю окно и соответсвенно N в 10 раз
Я получаю N=100 и ширину бина уже 100кгц, при этом сигнал будет физически воздействовать только 1 мкс из 10.
По факту я получу десять 100кгц бинов с нашим спектром.
Теперь самое главное что интересует, если сравнить сигнал _шум для нашего 1мгц бина и для суммы из десяти 100кгц бинов то сигнал_шум будет одинаковый в обоих случаях?
Поясняю зачем это мне надо. Кроме основного пульса в этой полосе возможно наличие постоянной несущей в виде помехи. Если я буду использовать 1мкс оптимальное окно(N=25) то я не получу разрешения по частоте меньше 1мгц
с другой стороны если я возьму 10 бинов и в одном из них будет помеха например на 300 кгц выше ценра я вырежу этот бин и оставлю 9 без помехи .
или как можно по другому убрать узколосную помеху с одной стороны и сохранить чувствительность в данной полосе с другой
есть неясность по поводу сигнал.шум . Попробую обьяснить
Допустим есть RF импульс 1мк , его полоса ~1мгц
Если взять окно длинной 1мкс получим при частоте 25мпс N=25
и полосу бина 1мгц
Вроде все замечательно, вкладываемся как раз в спектр
Теперь вопрос если при том же входном сигнале я увеличиваю окно и соответсвенно N в 10 раз
Я получаю N=100 и ширину бина уже 100кгц, при этом сигнал будет физически воздействовать только 1 мкс из 10.
По факту я получу десять 100кгц бинов с нашим спектром.
Теперь самое главное что интересует, если сравнить сигнал _шум для нашего 1мгц бина и для суммы из десяти 100кгц бинов то сигнал_шум будет одинаковый в обоих случаях?
Поясняю зачем это мне надо. Кроме основного пульса в этой полосе возможно наличие постоянной несущей в виде помехи. Если я буду использовать 1мкс оптимальное окно(N=25) то я не получу разрешения по частоте меньше 1мгц
с другой стороны если я возьму 10 бинов и в одном из них будет помеха например на 300 кгц выше ценра я вырежу этот бин и оставлю 9 без помехи .
или как можно по другому убрать узколосную помеху с одной стороны и сохранить чувствительность в данной полосе с другой
- Бахурин Сергей
- Администратор
- Сообщения: 1114
- Зарегистрирован: 05 окт 2010, 19:55
- Контактная информация:
Re: Гитберт при обнаружении короткого импульса
При увеличении длительности выборки отношение сигнал шум для белого аддитивного шума не изменится, потому что шум также размажется по большему числу бинов. А вот амплитуда узкополосного шума возрастет, но станет локализована по частоте и вы сможете вырезать часть узкополосного спектра с этой помехой.
После того как вы удалите УП помеху можно обнулить часть длинной выборки после 25 отсчётов. Тем самым вы существенно повысите отношение сигнал шум до 10 раз, потому что будете использовать временную локализацию полезного сигнала. Но делать так можно только после удаления УП помехи
После того как вы удалите УП помеху можно обнулить часть длинной выборки после 25 отсчётов. Тем самым вы существенно повысите отношение сигнал шум до 10 раз, потому что будете использовать временную локализацию полезного сигнала. Но делать так можно только после удаления УП помехи
Re: Алгоритм Герцеля при обнаружении короткого импульса
извините за повтор
Можно немного подробнее по последней теме. Идея понятна но я не совсем понимаю как это сделать правильно
1. Беру N=250 Герцелем
2. Анализирую 10 бинов в полосе 1мгц
3. Моя помеха например в 7 бине
4. удаляю его.
5. Как дальше получить из этих 9 бинов исходный сиглал без помехи что бы взять 25 первых выборок?
Можно немного подробнее по последней теме. Идея понятна но я не совсем понимаю как это сделать правильно
1. Беру N=250 Герцелем
2. Анализирую 10 бинов в полосе 1мгц
3. Моя помеха например в 7 бине
4. удаляю его.
5. Как дальше получить из этих 9 бинов исходный сиглал без помехи что бы взять 25 первых выборок?
- Бахурин Сергей
- Администратор
- Сообщения: 1114
- Зарегистрирован: 05 окт 2010, 19:55
- Контактная информация:
Re: Алгоритм Герцеля при обнаружении короткого импульса
2 варианта:
1. Поставить временной фильтр который будет убирать эти частоты и потом делать FFT. Так делать хорошо когда частота помехи известна.
2. FFT - убрать уп помеху - IFFT - Ограничение по времени - анализировать энергию во времени и принимать решение об обнаружении.
1. Поставить временной фильтр который будет убирать эти частоты и потом делать FFT. Так делать хорошо когда частота помехи известна.
2. FFT - убрать уп помеху - IFFT - Ограничение по времени - анализировать энергию во времени и принимать решение об обнаружении.
Re: Алгоритм Герцеля при обнаружении короткого импульса
про первый вариант думал сразу. Не получиться увы. Все та же проблемм короткий импульс.
У нас минимум несколько десятков первых отсчетов а наверное не меньше сотни уйдет на заполнения фильтра а это и есть наш пульс, мы его просто потеряем.
ихмо получаеться толькоти иФФТ.
Есть еще мысль по поду оконных функций.
Если применять окно на выборку больше пульса в 10 раз то окно подавит наш пульс в начале .
А что если взять всю выборку . поедлить ее на две части и переставить их что бы пульс был в средине нового общего окна
пульсу это не повредит а для гармонической помехи это будет эквивалентно скачку фазы . не очень хорошо но возможно но зато можно применинять окна
У нас минимум несколько десятков первых отсчетов а наверное не меньше сотни уйдет на заполнения фильтра а это и есть наш пульс, мы его просто потеряем.
ихмо получаеться толькоти иФФТ.
Есть еще мысль по поду оконных функций.
Если применять окно на выборку больше пульса в 10 раз то окно подавит наш пульс в начале .
А что если взять всю выборку . поедлить ее на две части и переставить их что бы пульс был в средине нового общего окна
пульсу это не повредит а для гармонической помехи это будет эквивалентно скачку фазы . не очень хорошо но возможно но зато можно применинять окна
- Бахурин Сергей
- Администратор
- Сообщения: 1114
- Зарегистрирован: 05 окт 2010, 19:55
- Контактная информация:
Re: Алгоритм Герцеля при обнаружении короткого импульса
Есть возможность записать 256 отсчётов из которых первые 25 это импульс? Желательно чтобы присутствовала УП помеха.
Re: Алгоритм Герцеля при обнаружении короткого импульса
да, без проблем могу гдето 1000 минимум сделать
- Бахурин Сергей
- Администратор
- Сообщения: 1114
- Зарегистрирован: 05 окт 2010, 19:55
- Контактная информация:
Re: Алгоритм Герцеля при обнаружении короткого импульса
Сделайте я посмотрю данные.