Оптимальный фильтр: 3. Понятие оптимальной фильтрации — Физический факультет СПбГУ
3. Понятие оптимальной фильтрации — Физический факультет СПбГУ
Оптимальное выделение сигнала из шума можно проводить различными методами, в зависимости от того, какая ставится задача — обнаружение сигнала, сохранение формы сигнала и т.д. В каждом методе оптимальной фильтрации вводится понятие критерия оптимальности, согласно которому строится оптимальный алгоритм обработки сигнала [ 11].
Конкретный алгоритм оптимальной фильтрации будет существенно зависеть от того непериодический или периодический сигнал должен быть выделен (обнаружен) на фоне шума. По отношению к периодическому сигналу далее различаются ситуации: известна или нет его частота повторения. Ниже эти варианты задач будут рассмотрены последовательно.
3.1 Оптимальная фильтрация непериодического (одиночного) сигнала
Оценим возможную эффективность обнаружения непериодического сигнала при его аддитивной смеси с белым шумом.
При формулировке задачи нахождения коэффициента передачи «оптимального» фильтра используются существенные требования относительно сигнала: во-первыхсчитается известной форма сигнала
(3.1)
Т.е. сигнал имеет конечную длительность.
Определение оптимальности фильтра формулируется следующим образом:
Оптимальным фильтром в задаче обнаружения одиночного импульса конечной длительности является фильтр, обеспечивающий максимальное отношение пиковой мощности сигнала к мощности шума в момент окончания импульса.
Комплексный коэффициент передачи такого оптимального фильтра прямо определяется спектром заданного, подлежащего обнаружению сигнала (т.е. его формой и длительностью ) [2 ]
Так, если сигнал имеет спектр
(3.
2)
и длительность его , то функция
есть функция комплексно сопряженная функции спектральной плотности сигнала.
Можно показать [2], что комплексный коэффициент передачи оптимального фильтра, в приведенном выше смысле определяется так:
или
(3.3)
Не воспроизводя выкладки доказательства формулы (3.3) приведенных в ряде источников (например [ 2] ), остановимся на физическом смысле результата.
Замечая, что фазовая характеристика коэффициента передачи в (3.3) есть
видим, что — компенсирует фазовые сдвиги составляющих сдвиги спектра сигнала (3.2), что формирует «пик» импульса на выходе, а линейная функция — обеспечивает задержку этого «пика» на время длительности сигнала, т. е. этот пик приходится на момент окончания сигнала.. Можно сказать, что обеспечивается накопление полезного сигнала на интервале всего времени существования импульса.
Формула (3.3) устанавливает также, что модуль коэффициента передачи должен совпадать с модулем спектральной плотности функции заданного сигнала, т. е. оптимальный фильтр ослабляет спектральные составляющие шума тем сильнее, чем меньше модуль , В результате полная мощность шума на выходе фильтра оказывается меньшей, чем при равномерной АЧХ.
Наконец отметим, что произвольная константа размерна. При безразмерном имеет размерность обратной спектральной плотности сигнала.
3.2. Оценка отношения сигнал/шум при оптимальном фильтре
Будем исходить из приведенного выше выражения (3.3). Заметим, что для сигнала (импульса) сложной формы синтез оптимального фильтра является не простой задачей. По этому искомую оценку отношения сигнал/шум проведем на примере прямоугольного импульса.
(3.
4)
Рис.1
Будем считать, что полезный сигнал представляет собой одиночный прямоугольный импульс (3.4) длительностью и с напряжением , изображенный на рис 1.
Его спектральная плотность описывается функцией
и сопряженная
следовательно, для данного импульса в соответствии с (3.3) имеем:
(3.5).
Такой коэффициент передачи может быть обеспечен схемой рис 2.
Рис.2
Содержащей идеальное интегрирующее звено (), линию задержки () и схему вычитания .
Найдем сигнал на выходе оптимального фильтра для рассматриваемого примера (=1)
(3.
6).
Таким образом, импульс на выходе имеет треугольную форму с основанием 2t и максимальным значением U0 при t=t . Оценим теперь мощность шума на выходе этого оптимального фильтра. Используем формулы (3.3) и (3.5). Положим a =1 Будем считать шум ’ белым’.
Рис.3
Представим модуль оптимального коэффициента передачи в виде
(3.7)
Мощность шума на выходе фильтра в соответствии с (2.2 ) и (2.3) определиться так
(3.8).
Этот табличный интеграл [3] имеет значение: . Таким образом, искомая величина мощности шума равна:
откуда искомое отношение мощности сигнала к мощности шума для данного оптимального фильтра будет:
(3.10).
а для отношения амплитуд сигнал/шум
(3.
11).
В заключение еще раз отметим, что оптимальный фильтр, построенный по указанному выше критерию, жестко связан с полезным входным сигналом: изменение полезного входного сигнала ведет к необходимости изменения коэффициента передачи фильтра.
3.3. Определение оптимальной полосы фильтра нижних частот в задаче выделения (обнаружения) одиночного сигнала на фоне белого шума
Учитывая сложность задачи синтеза оптимального фильтра , в результате которого находится его функция можно подойти к задаче по другому.
Сигналу выбирается тип АЧХ фильтра, сообразуясь с формой сигнала-импульса (точнее модулем его спектра). Например, для рассмотренного выше прямоугольного импульса выбирается фильтр НЧ, для импульса с высокочастотным заполнением- резонансный фильтр и т.д. Далее задача оптимизации ставится относительно выбора параметра фильтра — полосы его пропускания. Следуя этому подходу далее рассматривается возможность выделения полезного сигнала из белого шума не с помощью описанного выше оптимального фильтра, а с помощью линейного RC фильтра нижних частот.
Пусть полный входной сигнал U(t) выражается в виде суммы полезного входного сигнала и белого шума — случайного процесса, у которого спектральная плотность не зависит от частоты
(3.12).
В качестве фильтра нижних частот будем рассматривать интегрирующую цепочку (рис 4) — низкочастотный фильтр первого порядка с постоянной времени и коэффициентом передачи
(3.13).
Рис. 4
При исследовании прохождения шума через линейную систему будем использовать формулу (2.3) ,квадрат модуля коэффициента передачи
(3.14).
где — полоса пропускания рассматриваемого фильтра нижних частот по уровню 0.
707. Требуется найти полосу заданного фильтра нижних частот, обеспечивающую максимальное отношение сигнал/шум на выходе фильтра.
Можно рассматривать прохождение через фильтр нижних частот полезного сигнала и шума раздельно, так как интегрирующая цепочка — линейная схема.
3.3.1. Прохождение полезного сигнала через однозвенный RC фильтр нижних частот
Сигнал на выходе линейной системы может быть найден с помощью спектрального метода.
(3.15).
где коэффициент передачи интегрирующей цепочки определяется формулой (3.14) , а спектральная плотность полезного входного сигнала (3.4) была найдена как интеграл Фурье
(3.16).
Подставив в (3.15) формулы (3.16) и (3.13) и вычислив интеграл, получаем следующее выражение для сигнала на выходе фильтра
при
(3.
17)
Рис 5.
Таким образом, выходной сигнал достигает своего максимального значения в момент окончания входного импульса t=t
(3.18).
Это выражение зависит от соотношения полосы частот фильтра (3.13) и полосы частот, занимаемой полезным сигналом , которая связана с длительностью прямоугольного импульса так 1/t . С учетом этого выражение (3.18) можно преобразовать следующим образом
(3.19).
Если полоса частот, занимаемая спектральной плотностью полезного входного сигнала, меньше полосы частот, определяемой коэффициентом передачи интегрирующей цепочки , то максимальное значение полезного сигнала на выходе интегрирующей цепочки (3.18) равно и не зависит от полосы фильтра при .
Если же полоса частот, занимаемая спектральной плотностью полезного входного сигнала, больше полосы частот, занимаемой коэффициентом передачи интегрирующей цепочки, , то разложив экспоненту в выражении (3.
18) в ряд, получаем следующее максимальное значение полезного сигнала на выходе интегрирующей цепочки
В этом случае амплитуда полезного сигнала на выходе фильтра линейно зависит от полосы фильтра .
Мощность полезного входного сигнала, входящая в энергетическое отношение сигнал/шум, будет пропорциональна, таким образом, квадрату от полосы фильтра
Следовательно, если полоса фильтра перекрывает полосу полезного входного сигнала, то дальнейшее увеличение полосы фильтра не приводит к увеличению полезного выходного сигнала. Если же полоса фильтра уже полосы сигнала, то увеличение полосы фильтра приводит к увеличению мощности полезного выходного сигнала, пропорционально квадрату полосы фильтра.
3.3.2. Прохождение случайного сигнала (белого шума) через фильтр нижних частот
Для мощности шума на выходе фильтра с помощью формул (2.2) и (2.3) может быть получено следующее выражение
(3.
20),
в котором положим — спектральная плотность мощности белого шума, а квадрат модуля коэффициента передачи определен формулой (3.14). Вычислив интеграл , получаем
(3.21),
где — ширина полосы фильтра по уровню 1/.
Отсюда следует, что мощность шума на выходе фильтра линейно зависит от полосы коэффициента передачи интегрирующей цепочки.
Используя полученные выражения для максимального значения выходного полезного сигнала (3.18) и мощности шума на выходе фильтра (3.21), можно получить выражение для энергетического отношения сигнал/шум на выходе фильтра нижних частот (RC-цепочки):
(3.22),
где
и (3.23)
Искомую величину оптимальной полосы для выбранного НЧ фильтра (3.13) и сигнала (3.4), обеспечивающей максимальное отношение сигнал/шум в момент t=t , найдем из условия максимума функции (3.
23), т.е .
(3.24).
Функция имеет пологий максимум, ее график приведен на рис 6.
Рис. 6
Таким образом
И следовательно отношение пиковой мощности сигнала и мощности шума при оптимальной полосе НЧ фильтра равно
(3.25).
Отношение же амплитуды сигнала к «амплитуде» шума будет
(3.26).
Напомним, что использование фильтра с оптимальным коэффициентом передачи(3.5) приводило к отношению сигнал/шум по мощности равному (3.10)
(3.10).
Сравнивая (3.10) с (3.25), видим, что использование RC фильтра НЧ (3.13) с правильно выбранной полосой вместо фильтра с оптимальным коэффициентом передачи приводит к ухудшению соотношения сигнал/шум по мощности на 19 %.
(3.27),
и лишь на 10% по отношению амплитуд сигнал/шум
(3.28).
Т.о. для конкретного сигнала — прямоугольного импульса использование простого RC фильтра НЧ можно считать оправданным (целесообразным).
Качественно такой результат понятен. Если полоса фильтра уже полосы сигнала, то целесообразно увеличивать полосу фильтра, так как при этом мощность полезного сигнала на выходе растет пропорционально квадрату полосы, а мощность шума растет пропорционально первой степени полосы. Если полоса фильтра шире полосы сигнала, то целесообразно уменьшать полосу фильтра, так как при этом мощность полезного сигнала на выходе не меняется, а мощность шума уменьшается пропорционально первой степени полосы.
Далее найдем соотношение сигнал/шум для многозвенного RC — фильтра низкой частоты.
3.
4 Определение оптимальной полосы многозвенного фильтра нижних частот
Рассмотрим теперь задачу определения оптимальной полосы многозвенного фильтра с целью обеспечения максимального отношения сигнал/шум в момент окончания импульса. Импульс будем, как и раньше, считать прямоугольным.(3.4). Конкретно рассмотрим фильтр, собранный идентичных RC — звеньев, разделенных буферными каскадами.(рис 7).
Рис 7.
Коэффициент передачи такого фильтра описывается функцией
(3.29)
Если зафиксировать полосу пропускания этого фильтра на заданном уровне неравномерности , то эти два параметра, как это следует из (3.29), оказываются связанными уравнением
(3.30).
Отсюда очевидно, что для обеспечения постоянства общей заданной полосы фильтра при изменении числа звеньев n, постоянную необходимо изменить следующим образом
(3.
31).
При увеличении числа звеньев n, будет увеличиваться крутизна спада АЧХ в области частот, выше заданной полосы . (рис 8).
Рис8.
Переходная характеристика h(t) для рассматриваемого фильтра (3.29) — реакция на включение ступеньки напряжения на входе определяется так:
(3.32).
где , как отмечено выше, если при увеличении числа каскадов n ставится требование =const ( на заданном уровне неравномерности ), то параметр каждого каскада должен изменятся в соответствии с формулой (3.31). Пример зависимости от n для n=1 и n=5 при одинаковой приведен на рис.8, а зависимость переходной характеристики h(t) также при n=1 и n=5 приведен на рис.9.
Оценим теперь уровень шума на выходе фильтра. Считаем шум на выходе белым, имеющим спектральную плотность мощности S0
(3.
33).
Значение этого табличного неопределенного интеграла (3.36) известно [ 3].
(3.34).
При вычислении определенного интеграла (3.34) следует учесть, что функция равна нулю на верхнем () и на нижнем (-) пределах. Поэтому
(3.35).
Учитывая также необходимое изменение RC каждого каскада фильтра при увеличении n (при требовании =const) получаем интересующий количественный результат. В качестве примера приводим численные данные расчета мощности шума и напряжения шума для фильтров разных порядков (n).
и (3.36).
|
n |
1 |
2 |
3 |
4 |
5 |
6 |
|
0. |
0.78 |
0.85 |
0.9 |
0.95 |
1.01 |
|
|
0.5 |
0.61 |
0.72 |
0.81 |
0.91 |
1.025 |
707
Полученные выше данные о прохождении импульсного сигнала и белого шума через n-звенныйфильтр низкой частоты позволит определить оптимальную полосу фильтра ( при заданной длительности импульса) и соотношение сигнал/шум на выходе рассматриваемого фильтра в момент окончания импульс t=t при оптимальном выборе его полосы .
Как и в случае однозвенного фильтра строим функцию отношения величины полученного сигнала к «амплитуде» шума
(3.37).
Здесь , -коэффициент изменения RC= каждого звена, при изменении порядка фильтра n (3.31)
Отношение , как функция ,имеет пологий max, зависящий от порядка фильтра n.
Так, например, для трехзвенного фильтра max достигается при , и значении . Для пятизвенного фильтра получаем и . Из этих значений определяется оптимальная величина параметра фильтра .
Т.о. искомое отношение амплитуд сигнал/шум с учетом коэффициентов Ki в соответствии с выражением (3.37) дает:
Для трехзвенного фильтра (3.38),
Для пятизвенного фильтра (3.39).
Сравнивая эти результаты с полученными ранее, видим, что повышение порядка фильтра дает худшее отношение сигнал/шум, чем для фильтра первого порядка (3.
13):
и (3.40).
Поэтому, если «оптимальный» фильтр определяемый требованием (см(3.3)) заменяется фильтром RC с оптимально подбираемым параметром, то в рассматриваемом случае прямоугольного импульса лучшим оказывается простейший RC фильтр первого порядка .
Этому предпочтению можно дать следующие объяснения.
Во-первых, АЧХ RC фильтр первого порядка оказывается ближе к модулю спектра прямоугольного импульса, чем АЧХ фильтров более высоких порядков. Напомним, что для «оптимального» фильтра в соответствии с (3.3) оказывается, что его АЧХ совпадает с модулем спектра сигнала.
Во-вторых, как показано выше, значение обобщенного параметра , обеспечивающего наибольшее отношение мощности сигнала к мощности шума в конце импульса, увеличивается с повышением порядка фильтра. Так например, при n=1 значение , при n=3 имеем, что , а при n=5 обобщённый параметр .
Поэтому при заданной длительности импульса t полоса фильтра , требуется большей, для фильтра более высокого порядка, что также приводит к повышению шума.
Физически последнюю зависимость от n можно объяснить ростом группового запаздывания, т. е. требуемое максимальное отношение
к концу импульса достигается для больших n при больших значениях . Что при фиксированной длительности импульса означает большее значение . А это , естественно, увеличивает шум на выходе фильтра.
Приложение 3. Оптимальный линейный фильтр при белом шуме [3]
При обнаружении сигнала на фоне шума наибольшее распространение получил критерий максимума отношения сигнал / помеха на выходе фильтра. Для этого требуется синтезировать фильтр, обеспечивающий в заданный момент времени получение наибольшего возможного отношения пикового значения сигнала к среднеквадратичному значению шума.Сохранение формы сигнала не требуется, так как для обнаружения в шумах форма сигнала значения не имеет.
Будем считать, что сигнал детерминированный , известны его временная функция f(t) и, следовательно, спектр F(). Шум будем считать стационарным процессом с равномерным энергетическим спетром, то есть шум будем считать «белым» со спектральной плотностью мощности .
Задача создания такого оптимального фильтра решается в два этапа :
А) Определение его комплексного коэффициента передачи
Б) Определение структуры (схемы) фильтра и его параметров.
Коэффициент может быть найден для любого физически реализуемого сигнала. Второй этап более сложен и может не иметь решения. Необходимо оценить условия физической осуществимости искомого оптимального фильтра.
Перейдём к конкретному решению первого этапа — определению оптимального
коэффициента передачи фильтра в указанном выше смысле для заданного сигнала и при «белом» шуме.
Представив сигнал в виде интеграла Фурье
(п.3.1),
( где фазовая характеристика спектра сигнала)., а комплексный коэффициент передачи фильтра в виде
Запишем сигнал на выходе фильтра:
(п.3.2)
Пусть пик сигнала на выходе фильтра получается в момент (пока не известный).
(п.3.3)
Теперь оценим мощность шума на выходе этого фильтра.
(п.3.4).
Эффективное, среднеквадратичное напряжение шума равно .
Следовательно, отношение пикового значения сигнала к напряжению шума
(п.3.5).
Далее задача заключается в нахождении коэффициента передачи фильтра,
обеспечивающего максимум правой части выражения (3.
5). Воспользуемся неравенством Шварца. Для числителя выражения (3.5) можно записать:
(п.3.6)
Используя (3.6) , неравенство (3.5) можно записать так :
(п.3.7)
Отсюда видно, что отношение сигнал/шум достигает максимума, когда неравенство (3.7) обращается в равенство. Поскольку правая часть не зависит от , для этого необходимо выполнение двух условий (см. (3.3)):
1) чтобы . Откуда
(п.3.8).
2) чтобы , где А постоянный коэффициент. (п.3.9).
Отсюда следует , что коэффициент передачи оптимального фильтра должен иметь вид :
,
заметим, что есть спектральная функция,
сопряженная спектральной функции сигнала. Поэтому оптимальный коэффициент передачи фильтра будет определяться следующим образом:
.
(п.3.10).
Заметим, что если безразмерная функция, то коэффициент А должен иметь размерность обратную размерности спектра. Из формулы (3.10) следует, что . Оптимальный фильтр должен пропускать спектральные составляющие шума неравномерно, с тем большим ослаблением, чем меньше модуль спектра сигнала. В результате полная мощность шума будет меньшей, чем при равномерной амплитудной характеристике. Слагаемое фазовой характеристики в (3.8) компенсирует фазовые сдвиги в спектре сигнала. Слагаемое обеспечивает появление пика выходного сигнала в момент , то есть задержку выходного сигнала относительно его начала на входе. Только при условии может быть использована вся энергия сигнала ,имеющего длительность Т , для создания наибольшего возможного пика в момент.
Ясно, что дальнейшее увеличение запаздывания не влияет на величину пика.
Компания «Оптимум Фильтрация»
Компания «Оптимум Фильтрация»
| ||||||||||||||||||||||||||||||||||||||||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||
Сайт от Deep Fried Creative | ||||||||||||||||||||||||||||||||||||||||||||||
png»> 
Компания «Оптимум Фильтрация»
Компания «Оптимум Фильтрация»
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
./images/shadowlt.png»> 
Корпуса доступны либо с лентообразным хомутом, либо с поворотным болтом, они просты в использовании и обеспечивают превосходное герметичное уплотнение. Все внутренние детали из нержавеющей стали и верхняя пластина с принудительным сжатием исключают байпас жидкости.
