Пенсионер
Нет
УДК 53.087, 681.2
Введение.
Предлагаемый подход предусматривает использование метода нелинейной фильтрации сигнала, когда после дискретизации и оцифровки, его периодические неоднородности существенно сглаживаются, а разовые аномалии усиливаются.
Метод основан на разрешении противоречия между необходимостью сгладить изменения сигнала, обусловленные структурой его неоднородностей, фоновыми шумами, и необходимостью усилить изменения сигнала, обусловленные присутствием в среде аномальной (не свойственной ее равновесному состоянию) неоднородности путем определенной нелинейной фильтрации.
Цель статьи.
Представить результаты моделирования оригинальной фильтрации временной структуры сложных, негармонических сигналов для выделения аномалий, полученных от различных непрерывных сред.
Научная новизна.
Предложен метод фильтрации сигнала, позволяющей выявить аномалии в поведении частично организованной структуры сигнала от протяженной среды.
Методология
Алгоритм такой фильтрации предполагает следующие этапы:
1) Деление пополам текущего временного интервала регистрации импульсного сигнала от среды, задаваемого с определенным шагом дискретизации.
2) Вычисление логарифма отношения сигналов для таких интервалов, при нарастании интервалов регистрации сигнала в процессе его распространения, с шагом равным удвоенному шагу дискретизации принимаемого сигнала.
3) Умножение полученного ряда дискретных значений сигнала на отношение числа отсчетов сигнала от среды к числу отсчетов на взятом интервале.
4) Выделение неоднородности среды, как аномалии в поведении сигнала, путем соотнесения провалов во временном ходе отфильтрованного сигнала с границами такой аномалии, когда чем уже и глубже такой провал, тем контрастней аномалия.
Для формирования модельной структуры неоднородностей среды по выделенной аномалии сигнала, соответствующей конкретной ситуации используется следующая последовательность действий 59-61:
Из дискретных отсчетов сигнала формируют две новые последовательности дискретных отсчетов путем их суммирования, для первой, каждый раз начиная с первого отчета и до заданного (тип А), с нарастанием от начала и до конца выборки; для второй, каждый раз начиная с последнего отсчета и до заданного, с нарастанием от конца до начала выборки (тип В), при шаге приращения равном требуемому разрешению,
Для каждого сформированного таким образом отчета формируют дополнительный отсчет (тип AD и тип BD соответственно), путем
суммирования такого же числа последующих (в первом случае) или предшествующих исходных отсчетов (во втором случае),
Делят меньшее значение из такой пары отсчетов на большее (AD/A, B/BD).
результаты деления логарифмируют и делят на число исходных значений сигнала, соответствующее каждой паре, получая таким образом две новые последовательности отсчетов (IAn и IBn), число
которых равно числу исходных значений сигнала, а их значения пропорциональны средним коэффициентам ослабления среды на соответствующих участках трассы σ(An), σ(Bn),
Берут отношение для пар отсчетов разного типа IA1/IBn, IA2/IBn…, IAn/IB1, когда им соответствующие участки среды имеют общую часть Δ, равную половине шага необходимой дискретизации сигнала.
В сформированной последовательности находят глобальный минимум IAk/IBn–k + 1 = min и число исходных значений сигнала до общей части Δ соответствующих участков трассы
Результаты
Первичная проверка эффективности специальной фильтрации по описанному алгоритму проводилась в численном эксперименте при сравнении ее с линейной и медианной фильтрациями [1,2,3]. Отличие исследуемой фильтрации - использование при построении алгоритма априорной информации о возможном скачке сигнала обратного рассеяния на шаге его дискретизации (на основе, знания характера изменения оптических характеристик реальной атмосферы).
Численное моделирование проводилось с применением интерактивной системы спектрально-статистической обработки зашумленного модельного сигнала. Использовался экспоненциально-спадающий сигнал с перепадом в центре выборки из 50 отсчетов. Результаты численной обработки такого модельного сигнала приведены ниже.
Рис.1. Модельный сигнал с перепадом и наложенной структурой (верхний график) и результаты фильтрации выявляющие такой перепад (нижний график), где по оси ординат отложена величина сигнала, а по оси абсцисс – число отсчетов.
Критерием эффективности фильтрации служило отношение дисперсий отфильтрованного и исходного (без шумов) сигналов как для всей выборки, так и ее частей до и после перепада. Пример специальной фильтрации другого модельного сигнала программным образом приведен ниже.
Рис.2. Модельный сигнал с перепадом и наложенной структурой (верхний график) и результаты фильтрации выявляющие такой перепад (нижний график), где по оси ординат отложена величина сигнала, а по оси абсцисс – число отсчетов.
Следует отметить, что в результате фильтрации граница перепада четко выделяется, сигнал после перепада выше его исходного значения.
В этом случае линейная и медианная фильтрации не эффективны и не позволяют отождествить введенный перепад модельного сигнала.
Общие оценки эффективности фильтрации представлены в таблице 1, где введены следующие обозначения: LINF – линейный фильтр, MEDF - медианный фильтр, GRAN(0), GRAN(1) – специальная фильтрация (0,1 – число повторных поисков максимальных перепадов отфильтрованного сигнала), DISP1, DISP2 – дисперсии до и после введенного перепада сигнала, SIG1, 2, 3 – сигналы при отношениях между максимальным значением, перепадом и шумом: 10:1:1; 10:1:5; 10:1:10 соответственно.
|
SIG 1 |
DISP1 |
DISP2 |
LINF MEDF GRAN(0) GRAN (1) |
1,021 1,008 33,62 3,738 |
0,990 0,992 15,90 2,000 |
1,198 1,055 1,055 2,016 |
|
SIG 2 |
DISP1 |
DISP2 |
LINF MEDF GRAN(0) GRAN (1) |
1,131 1,070 2,438 4,651 |
1,198 1,148 2,420 3,148 |
2,888 2,492 1,420 |
|
SIG 3 |
DISP1 |
DISP2 |
LINF MEDF GRAN(0) GRAN (1) |
1,487 1,400 5,193 4,789 |
1,494 1,383 6,908 4,320 |
109 2,893 3,969 4,355 |
В таблице значения дисперсии нормированы относительно дисперсии исходного сигнала. Для специальной фильтрации DISP1 и DISP2 существенно выше единицы, что указывает на нелинейный характер фильтрации, высокую чувствительность к отдельным перепадам сигнала.
Рис.3. Результаты сравнения предлагаемой обработки приходящего сигнала рассеяния (нелинейной фильтрации - GRAN) с методами линейной и медианной фильтраций (LINF, MEDF).
Сопоставления
Поиск аналогов введенной фильтрации привел к методам цифровой обработки сигналов в оптике и голографии, конкретнее к методам нелинейной фильтрации и робастным системам.
Такую фильтрацию можно отождествить с рекурсивной фильтрацией при элементах робастности (нечувствительности к влиянию небольшого числа ложных срабатываний и малым изменениям входных данных). Она имеет общие черты со следующими методами фильтрации:
Однако предлагаемая фильтрация имеет следующие преимущества:
На основе проведенного анализа можно сделать заключение, что данный тип фильтрации обеспечивает синтез преимуществ разных типов фильтрации, при избавлении от многих ограничений им свойственных.
Выводы
Результаты численного моделирования, с генерацией шумовой компоненты, или дополнительной синусоидальной компоненты, с периодом соизмеримым с протяженностью взятой выборки сигнала, подтверждают перспективность применения предлагаемого подхода по нелинейной фильтрации полученного от среды сигнала, который позволяет успешно выделять перепады сигнала, обусловленные аномальной неоднородностью среды, на фоне определенным образом организованной структуры неоднородностей, или в условиях значительных фоновых шумов, когда известные методы не приносят успеха.
Дополнительными преимуществами такого анализа являются слабая зависимость результатов от длины выборки протяженного сигнала, простота реализации, высокое быстродействие, возможность органично использовать цифровую технику при аппаратурной реализации.
Рецензии:
11.11.2022, 9:46 Ашрапов Улугбек Товфикович
Рецензия: В статье "Метод выявления аномальной неоднородности в структуре шумоподобного сигнала" автором проведены исследования которые могут иметь применение его подхода по нелинейной фильтрации сигнала, который позволяет выделять перепады сигнала, обусловленные аномальной неоднородностью среды, на фоне организованной структуры неоднородностей, или в условиях значительных фоновых шумов. Статью рекомендую к публикации.