Оптимальная линейная фильтрация непрерывных сигналов

Рассмотрим теорию оптимальной линейной фильтрации стационарных процессов – фильтр Колмогорова-Винера или оптимальный линейный фильтр (ОЛФ).

Пусть на входе линейного фильтра с передаточной функцией H (j w) действует сумма полезного сигнала s (t ) и помехи n (t ). Отклик фильтра на это действие – восстановленный полезный сигнал (оценка сигнала s (t )). Будем считать, что s (t ) и n (t ) – стационарные взаимнонекоррелированные процессы с известными спектральными плотностями мощности (СПМ) G s (f ) и G n (f ). Нужно найти такую функцию H (j w), которая обеспечивает минимум среднего квадрата ошибки восстановления сигнала

. (22.1)

Иначе говоря, критерием оптимальности фильтра является минимум среднего квадрата ошибки восстановления сигнала. В такой постановке задача была решена А.Н. Колмогоровым (1939 г.) для дискретных случайных последовательностей и Н. Винером (1941 г.) для непрерывных процессов. Поэтому оптимальный (в указанном смысле) линейный фильтр называется фильтром Колмогорова-Винера.

Действующие на входе фильтра сигнал s (t ) и помеха n (t ) проходят через фильтр независимо и создают на выходе фильтра соответственно фильтрованные сигнал s вых (t ) и помеху n вых (t ). С учетом этого ошибка восстановления запишется

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

. (22.3)

Величина линейных искажений полезного сигнала фильтром зависит от степени отличия АЧХ фильтра от постоянной величины и степени отличия ФЧХ от линейной зависимости. Средний квадрат шума на выходе фильтра зависит только от АЧХ фильтра. Для того чтобы линейные искажения полезного сигнала были минимальными, ФЧХ фильтра должна быть линейной

j(w) = –wt 0 , (22.4)

где t 0 – задержка сигнала в фильтре. Ясно, что с учетом задержки соотношение (22.1) имеет вид . Это уточнение не влияет на критерий оптимальности, поскольку в системах связи и вещания ожидаемая задержка сигнала в фильтре несущественная.

Перейдем к определению АЧХ фильтра. Для этого определим спектральные плотности мощностей левой и правой частей соотношения (22.2)

Выразим СПМ помехи на выходе фильтра через СПМ помехи n (t ) и искомую АЧХ фильтра:

. (22.6)

По определению СГП ергодичного процесса

, (22.7)

где S D s (w) – амплитудный спектр ошибки за счет линейных искажений Ds (t );

Т – длительность реализации Ds (t ).

Поскольку амплитудный спектр сигнала s в ых (t ) определяется как S (w)×H (w), где S (w) – амплитудный спектр сигнала s (t ), то

Переходя от амплитудных спектров к СПМ, получим

После подстановки соотношений (22.6) и (22.9) в (22.5) получим

Средний квадрат ошибки восстановления (средняя мощность) вычисляется

. (22.11)

Поскольку функция G e (w) ³ 0 на всех частотах, то, обеспечив min G e (w) на всех частотах, достигнем минимума величины . Искомую АЧХ H (w) определим из условия экстремума функции G e (w):

После решения уравнения (22.13) получим выражение для АЧХ фильтра

. (22.14)

На рис. 22.1 иллюстрируется АЧХ ОЛФ, определенная соотношением (22.14).

Из рис. 22.1 видно особенности АЧХ ОЛФ:

На частотах, где G n (f ) = 0, значение АЧХ H (f ) = 1 – в этих областях частот фильтр не должен вносить искажений;

На частотах, где G s (f ) = 0, значение АЧХ H (f ) = 0 – в этих областях частот фильтр должен полностью ослаблять составляющие помехи;

На частотах, на которых G s (f ) = G n (f ), АЧХ H (f ) = 0,5;

На других частотах значения АЧХ определяются вычислениями по формуле (22.14).

Подставим выражение (22.14) в соотношение (22.10) для определения СПМ ошибки:

(22.15)

При подстановке соотношения (22.15) в выражение (22.11) можно вычислить средний квадрат ошибки восстановления сигнала .

Из (22.15) видно, что ошибка равна нулю только в том случае, когда G s (f )G n (f ) = 0, т.е. когда спектры сигнала и помехи не перекрываются (хотя бы один с сомножителей равен нулю).

Перепишем соотношение (22.14) в виде

. (22.16)

Из последнего соотношения видно, что коэффициент передачи оптимального фильтра на каждой из частот тем меньше, чем больше отношение G n (f )/G s (f ) на этой частоте.

Следует отметить, что оптимальные линейные фильтры, обеспечивающие минимум ошибки , существенным образом отличаются от согласованных фильтров, рассмотренных ранее. Если основное назначение рассмотренных здесь фильтров состоит в наилучшем воспроизведении формы сигнала, то задача согласованных фильтров состоит в формировании максимального отношения сигнал/шум в момент отсчета.

При использовании ОЛФ в аналоговых системах связи и вещания выявляется такая особенность. Имеет место высокое отношение спектральных плотностей сигнала и шума: G s (f )/G n (f ) >> 1. Выражение для АЧХ ОЛФ (22.14) в случае полосовых сигналов переходит в следующее

(22.17)

где f min и f max – граничные частоты спектра сигнала. В случае НЧ сигналов

(22.18)

где F max – максимальная частота спектра сигнала.

Таким образом, оптимальные линейные фильтры в системах связи и вещания имеют П-образную АЧХ.

Контрольные вопросы

1. Сформулируйте критерий оптимальности ОЛФ.

2. При каком условии ошибку восстановления сигнала можно свести к нулю?

3. Объясните отличие ОЛФ и согласованного фильтра с точки зрения ослабления помех.

23. Сравнение помехоустойчивости оптимальных
демодуляторов сигналов аналоговых видов модуляции

Мы выяснили, что в условиях слабых помех демодулятор должен содержать: фильтр додетекторной обработки, детектор, фильтр последетекторной обработки. Для того чтобы демодулятор был оптимальным, фильтры должны быть оптимальными. В условиях слабой помехи АЧХ фильтров должны быть П-образными:

– фильтр додетекторной обработки – полосовой фильтр, граничные частоты полосы пропускания которого совпадают с граничными частотами спектра модулированного сигнала;

– фильтр последетекторной обработки – ФНЧ, частота среза которого совпадает с максимальной частотой спектра первичного сигнала F max .

Помехоустойчивость определим в условиях действия АБГШ. Анализ помехоустойчивости состоит в определении выигрыша демодулятора в отношении сигнал/шум

Для определения выигрыша нужно определить 4 величины: P b , P e , P s , P n .

Оптимальныйдемодулятор сигнала балансовой модуляции . Математически сигнал БМ записывается

s БМ (t ) = A 0 b (t )cos2pf 0 t . (23.2)

Полосовой фильтр имеет полосу пропускания DF БМ = 2F max . Для детектирования БМ сигнала необходимо использовать синхронный детектор (рис. 23.1). ФНЧ, являющийся обязательным элементом схемы синхронного детектора, используется как фильтр последетекторной обработки, т.е. частота среза фильтра равна F max , а АЧХ в полосе пропускания постоянная и равна 1.

P s = = = 0,5 P b . (23.3)

Средняя мощность помехи на входе демодулятора P n .

На выходе перемножителя за счет сигнала имеем

A 0 b (t )cos2pf 0 t ×2cos2pf 0 t = A 0 b (t ) + A 0 b (t ) cos2p2f 0 t . (23.4)

ФНЧ пропускает первый компонент, а второй ослабляет. С учетом умножения на 1/A 0 на выходе демодулятора получим b (t P b .

Помеху на входе синхронного детектора (как полосовой процесс) представим квадратурными составляющими

n (t ) = N C (t )cos2pf 0 t + N S (t )sin2pf 0 t . (23.5)

Мощность помехи делится поровну между квадратурными составляющими, мощность каждой из них P n /2. Квадратурная составляющая помехи не проходит через синхронный детектор, и на выходе демодулятора получим

e(t ) = N C (t )/A 0 . (23.6)

Поскольку = P n /2, а левая часть равенства равняется 0,5 , то = P n , а

P e = P n / . (23.7)

g БМ = = 2. (23.8)

Оптимальныйдемодулятор сигнала однополосной модуляции . Математически сигнал ОМ записывается

s ОМ (t ) = A 0 b (t )cos2pf 0 t ± A 0 sin2pf 0 t . (23.9)

Оптимальный демодулятор сигнала ОМ выполняется по схеме, приведенной на рис. 23.1. Полосовой фильтр имеет полосу пропускания DF ОМ = F max .

Средняя мощность модулированного сигнала

P s = = P b . (23.10)

Средняя мощность помехи на входе синхронного детектора P n .

Синхронный детектор не пропускает квадратурную составляющую сигнала ОМ, поэтому на основе анализа демодуляции сигнала БМ на выходе демодулятора сигнала ОМ получим сигнал b (t ). Его средняя мощность равняется P b . Прохождение помехи через синхронный детектор проанализирован выше и получено значение мощности помехи на выходе демодулятора (23.7).

Определим выигрыш демодулятора

g ОМ = = 1. (23.11)

Оптимальныйдемодулятор сигнала амплитудной модуляции на базесинхронного детектора . Математически сигнал АМ записывается

s АМ (t ) = A 0 (1 + m АМ ×b (t ))cos2pf 0 t . (23.12)

Схема демодулятора сигнала АМ приведена на рис. 11.5. Полосовой фильтр имеет полосу пропускания DF АМ = 2F max . Исходя из приведенного выше анализа, очевидно, что на выходе ФНЧ за счет сигнала получим A 0 + A 0 m АМ ×b (t ). Фильтр верхних частот ослабляет постоянную составляющую A 0 и пропускает вторую составляющую A 0 m АМ ×b (t ).

Средняя мощность модулированного сигнала

P s = = 0,5 (1 + m АМ P b ). (23.13)

Для дальнейшего анализа удобно учесть (см. модуль 1), что P b = 1/ , где K A – коэффициент амплитуды сигнала b (t ).

Прохождение шума через синхронный детектор проанализирован выше. Мощность шума на выходе демодулятора сигнала АМ

P e = . (23.14)

Определим выигрыш демодулятора

= . (23.15)

Демодулятор сигнала амплитудной модуляции на базедетектора огибающей . Схема такого демодулятора приведена на рис. 23.3. Она отличается от схемы демодулятора на рис. 23.2 типом амплитудного детектора – синхронный детектор заменен детектором огибающей с целью упрощения схемы демодулятора. Выходной сигнал детектора огибающей пропорциональный огибающей входного сигнала .

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

(23.16)

Оптимальныйдемодулятор сигнала фазовой модуляции . Математически сигнал ФМ записывается

s ФМ (t ) = A 0 cos(2pf 0 t + Δj д ∙b (t )), (23.17)

где Δj д – девиация фазы сигнала, которую часто называют индексом фазовой модуляции m ФМ.

Схема оптимального демодулятора сигнала ФМ приведена на рис. 23.4: ФНЧ1 и ФНЧ2 имеют частоты среза F max (m ФМ + 1); ФНЧ3 – фильтр последетекторной обработки с частотой среза F max ; АЧХ фильтров постоянная в полосах пропускания и равна 1; ФД – фазовый детектор.

Средняя мощность модулированного сигнала

. (23.18)

Средняя мощность исходного сигнала P b = 1/ .

Средняя мощность помехи на входе ФД P n . Прохождение помехи через ФД анализируют при отсутствии модулирующего сигнала, т.е. b (t ) º 0, и s (t ) = A 0 cos2pf 0 t . Помеху представляют квадратурными составляющими в виде (23.5). Тогда

В условиях слабой помехи |N c (t )| << A 0 , и помеха на выходе ФД N s (t )/(A 0 m ФМ), а ее мощность равна P n /(A 0 m ФМ) 2 .

Полоса пропускания ФНЧ3 в раз меньше, чем полоса пропускания ФНЧ2 и ширина спектра помехи N s (t ).

Поскольку спектр помехи равномерный, то мощность помехи e(t ) в m ФМ + 1 раз меньше, чем мощность помехи N s (t )/(A 0 m ФМ) и определяется

. (23.20)

Определим выигрыш демодулятора

Оптимальныйдемодулятор сигнала частотной модуляции. Математически сигнал ЧМ записывается

где Δf д – девиация частоты. Для последующего изложения удобно использовать индекс ЧМ m ЧМ = Δf д /F max .

Схема оптимального демодулятора сигнала ЧМ приведена на рис. 23.5. Схема отличается от схемы демодулятора ФМ сигнала (рис. 23.4) наличием дифференциатора; ЧД – частотный детектор.

Средняя мощность модулированного сигнала

. (23.23)

Средняя мощность сигнала на выходе демодулятора P b =1/ .

Средняя мощность помехи на входе ЧД P n .

Прохождение помехи через ЧД аналогично прохождению помехи через ФД, следует рассмотреть прохождение помехи N s (t )/(A 0 2pΔf д) через дифференциатор. Поскольку помеха N s (t ) – квазибелый шум в полосе частот F max (m ЧМ + 1), то спектральная плотность мощности этой помехи на входе дифференциатора

Поскольку передаточная функция дифференциатора j ω, то спектральная плотность мощности помехи на выходе дифференциатора

Определим мощность помехи e(t )

Определим выигрыш демодулятора

При анализе мы выявили, что спектральная плотность мощности помехи e(t ) имеет параболическую зависимость – формула (23.24):

G e (f ) = kf 2 , (23.26)

где k – коэффициент пропорциональности. Эта особенность спектра часто учитывается при разработке систем передачи с ЧМ.

Сравнениеаналоговых систем передачи. Основными параметрами, по которым сравниваются системы передачи, является выигрыш демодулятора в отношении сигнал/шум g и коэффициент расширения полосы частот при модуляции a = ΔF s /F max. Для рассмотренных методов модуляции эти параметры сведены в табл. 23.1.

Проведем сравнение числовых значений параметров при некоторых исходных данных: K А = 5; m ЧМ = m ФМ = 6; m АМ = 1.

Вычисления дают: g АМ = 0,038; g БМ = 2; g ОМ = 1; g ЧМ = 60,5; g ФМ = 20,2.

Сравнение числовых значений выигрыша показывает, что самой низкая помехоустойчивость присуща системе передачи с АМ: выигрыш g АМ << 1, что логически назвать проигрышем. Однако АМ используется в системах радиовещания, где этот недостаток компенсируется простотой демодулятора на основе детектора огибающей (огромное количество более простых радиоприемников и один радиопередатчик с большей мощностью, чем при использовании БМ и ОМ).

Наибольшая помехоустойчивость присуща системе передаче с ЧМ. «Платой» за высокую помехоустойчивость является широкая полоса частот сигнала. Так, при F max = 3,4 кГц ΔF ЧМ = 47,6 кГц, в то время как полоса частот сигнала ОМ ΔF ОМ = 3,4 кГц.

Таблица 23.1 – Основные параметры аналоговых систем передачи

Метод модуляции	g	a	Примечания
АМ			Синхронный детектор
		Детектор огибающей
БМ
ОМ
ЧМ	×a	2(m ЧМ + 1)	r вх > r пр
ФМ	×a	2(m ФМ + 1)	r вх > r пр

Порогпомехоустойчивости демодулятора сигнала ЧМ. Из соотношения для выигрыша демодулятора ЧМ (23.25) вытекает, что, чем больше индекс m ЧМ, тем больше выигрыш (правда, ценой увеличения полосы частот сигнала). Может показаться, что это дает возможность работать демодулятору со слабым сигналом (низким отношением сигнал/шум). Но, когда отношение сигнал/шум на входе демодулятора r вх меньше порогового отношения сигнал/шум r пр, то выигрыш демодулятора резко уменьшается (рис. 23.6). Такое явление резкого уменьшения величины выигрыша называют порогом помехоустойчивости приема сигнала ЧМ.

Пороговое отношение сигнал/шум r пр несколько зависит от значения m ЧМ (рис. 23.6). Считают, что демодулятор по схеме стандартного частотного детектора характеризуется ориентировочным значением r пр = 10. Область значений r вх, когда r вх < r пр, – это нерабочая область.

Были предложены так называемые порогопонижающие схемы демодуляторов сигналов ЧМ, которые получили названия:

Демодулятор со следящим фильтром;

Демодулятор с обратной связью по частоте;

Демодулятор на основе синхронно-фазового детектора.

Схемы этих демодуляторов описаны в специальной литературе. Демодуляторы, которые выполнены по таким схемам, характеризуются пороговым отношением сигнал/шум r пр = 5...7 дБ (в зависимости от исходных данных на систему передачи). Снижение r пр разрешает:

1) работать демодулятору с более низким отношением сигнал/шум;

2) увеличить выигрыш, так как

,

и если допустить уменьшение отношения сигнал/шум r вх, то можно увеличить индекс m ЧМ, а увеличение m ЧМ приводит к увеличению выигрыша.

Литература

Основная

1. Стеклов В.К., Беркман Л.Н. Теорія електричного зв’язку: Підручник для ВНЗ під редакцією В.К. Стеклова. – К.: Техніка, 2006.

2. Теория электрической связи: Учебник для вузов / А.Г. Зюко, Д.Д. Кловский, В.И. Коржик, М.В. Назаров; Под ред. Д.Д. Кловского. – М.: Радио и связь, 1998.

Условие оптимальности фильтра. Фильтр Колмогорова-Винера является оптимальным фильтром формирования из входного сигнала x(t) выходного сигнала z(t) при известной форме полезного сигнала s(t), который содержится во входном сигнале в сумме с шумами. В качестве критерия его оптимизации используется среднее квадратическое отклонение сигнала y(t) на выходе фильтра от заданной формы сигнала z(t). Подставим уравнение свертки (12.2.1) в раскрытой форме весового суммирования в выражение (12.2.2") и получим отклонение  2 выходного сигнала y(k) = h(n)③x(k-n) от заданной формы выходного сигнала z(k) по всем точкам массива данных:

В частном случае воспроизведения формы полезного сигнала в качестве функции z(k) принимается функция s(k). Минимум выражения (12.3.1) определяет значения коэффициентов h(n) оптимального фильтра. Для нахождения их значений продифференцируем выражение (12.3.1) по коэффициентам фильтра и приравняем полученные уравнения нулю. В итоге получаем:

где
- корреляционная функция входного сигнала,
- взаимная корреляционная функция сигналов z(k) и x(k). Отсюда:

h n R(m-n) = B(m), n = m = 0,1,2, ... , M. (12.3.2)

В краткой форме символической записи:

h(n) ③ R(m-n) = B(m). (12.3.3)

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

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

m=0: h o R(0)+ h 1 R(1)+ h 2 R(2)+ h 3 R(3)+ ...+ h M R(M) = B(0), (12.3.3")

m=1: h o R(1)+ h 1 R(0)+ h 2 R(1)+ h 3 R(2)+ ...+ h M R(M-1) = B(1),

m=2: h o R(2)+ h 1 R(1)+ h 2 R(0)+ h 3 R(1)+ ...+ h M R(M-2) = B(2),

..............................................................................................................

m=M: h o R(M)+ h 1 R(M-1)+ h 2 R(M-2)+ .... + h M R(0) = B(M).

Решение данной системы уравнений относительно значений h i дает искомый оператор фильтра.

При расчете каузальных (односторонних) фильтров выходной сигнал z(t) следует задавать со сдвигом вправо по оси координат таким образом, чтобы значимая часть функции взаимной корреляции B(m) полностью располагалась в правой части системы уравнений (12.3.3").

Отметим, что R(m) = R s (m)+R q (m), где R s - функция автокорреляции сигнала, R q - функция автокорреляции шума, а B(m) = B zs (m)+B zq (m), где B zs - функция взаимной корреляции сигналов z(k) и s(k), B zq - функция взаимной корреляции сигнала z(k) и помех q(k). Подставляя данные выражения в (12.3.3), получаем:

h(n) ③ = B zs (m)+B zq (m). (12.3.4)

Частотная характеристика фильтра находится преобразованием Фурье левой и правой части уравнения (12.3.4):

H() = W zs ()+W zq (),

H() = / , (12.3.5)

где W s ()  R s (m) и W q ()  R q (m) - энергетические спектры (плотности мощности) сигнала и помех, W zs ()  B zs (m) - взаимный энергетический спектр входного и выходного сигналов, W zq ()  B zq (m) - взаимный энергетический спектр выходного сигнала и помех.

В геофизической практике обычно имеет место статистическая независимость полезного сигнала, а, следовательно, и сигнала z(k), от шумов, при этом B zq = 0 и фильтр называют оптимальным по сглаживанию шумов при заданной форме выходного сигнала:

H() = W zs () / , (12.3.6)

Фильтр (12.3.6) оптимален в том смысле, что максимизирует отношение мощности сигнала к мощности шума по всему интервалу сигнала, но не в каждой индивидуальной точке.

Выражения (12.3.5-12.3.6), как правило, всегда имеют решение. Однако это не означает возможность формирования фильтром любой заданной формы выходного сигнала. Из чисто практических соображений можно сразу предполагать, что если спектр заданного сигнала z(t) больше значимой части спектра полезного сигнала s(t), то оператор фильтра попытается сформировать требуемые высокие частоты заданного сигнала из незначимых частот спектра полезного сигнала, что может потребовать огромных коэффициентов усиления на этих частотах, под действие которых попадут и частотные составляющие шумов. Результат такой операции непредсказуем. Эти практические соображения можно распространить и на все частотные соотношения входного и выходного сигналов линейных фильтров: значимые гармоники спектров выходных сигналов должны формироваться из значимых гармоник спектров входных сигналов.

Если заданная форма сигнала z(k) совпадает с формой полезного сигнала s(k), то B(m) = B ss = R s и фильтр называют фильтром воспроизведения полезного сигнала :

H() = W s () / , (12.3.7)

Выражения (12.3.6-12.3.7) достаточно наглядно демонстрируют физический смысл формирования передаточной функции фильтра. При воспроизведении сигнала частотная функция взаимной корреляции входного сигнала с выходным W zs (плотность взаимной мощности) повторяет частотную функцию автокорреляции W s (плотность мощности сигнала). Плотность мощности статистических шумов W q распределена по частотному диапазону равномерно, в отличие от плотности мощности сигнала W s , которая, в зависимости от формы сигнала, может занимать любые частотные интервалы спектрального диапазона. Частотная передаточная функция фильтра воспроизведения сигнала формируется отношением W s ()/. На частотах, где сосредоточена основная энергия сигнала, имеет место W s ()>>W q () и H()  1 (как минимум, больше 0.5). Там, где значение W s () становится меньше W q , коэффициент передачи фильтра становится меньше 0.5, и в пределе H()=0 на всех частотах, где полностью отсутствуют частотные составляющие сигнала.

Аналогичный процесс имеет место и при формировании произвольного сигнала z(t) на выходе фильтра, только в этом случае из частот входного сигнала устанавливаются на выделение и усиление частоты взаимной мощности входного и выходного сигнала, необходимые для формирования сигнала z(t), причем коэффициент на этих частотах может быть много больше 1, а подавляться могут не только шумы, но и частоты основного сигнала, если их нет в сигнале z(t).

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

Задание мощности шумов. Следует внимательно относиться к заданию функции шумов Wq(). При полной неопределенности шума необходимо, как минимум, выполнять оценку его дисперсии  2 и распространять на весь частотный диапазон с соответствующей нормировкой на его величину (2Wq() d =  2), т.е. считать его белым шумом. При известной функции полезного сигнала в зарегистрированной реализации значение дисперсии шумов в первом приближении может быть оценено по разности дисперсий реализации и функции полезного сигнала. Можно выполнить и выделение шумов из входного сигнала в отдельный шумовой массив, например, вейвлетным преобразованием. Однако использовать выделенный шум непосредственно для вычисления функции Wq() допустимо только по достаточно представительной реализации при условии стационарности и эргодичности шума. В противном случае функция Wq() будет отображать только распределение шумов в зарегистрированной реализации сигнала, а соответственно фильтр будет оптимален только к этой реализации, что не гарантирует его оптимальности к любой другой реализации. Но для обработки единичной зарегистрированной реализации сигнала такой метод не только вполне допустим, но и может существенно повысить точность формирования выходного сигнала.

Эффективность фильтра. Из выражений (12.3.5-12.3.7) следует, что с позиции минимального искажения полезного сигнала при максимальном подавлении шумов фильтр Колмогорова-Винера эффективен в тем большей степени, чем больше отношение сигнал/шум на входе фильтра. В пределе, при W q ()<>W s () имеем H() 0 и сигнал будет сильно искажен, но никакой другой фильтр лучшего результата обеспечить не сможет.

Пример. Расчет оптимального фильтра воспроизведения сигнала. Выполняется в среде Mathcad .

Форма входного сигнала считается известной и близка к гауссовой. На входной сигнал наложен статистический шум с равномерным распределением мощности по всему частотному диапазону (белый шум), некоррелированный с сигналом, и функцию Wzq принимаем равной нулю. Для наглядного просмотра связи параметров фильтра с параметрами сигнала задаем модели сигналов в двух вариантах:

K:= 1000 k:= 0 .. K A:= 50

s1 k:= A·exp[-0.0005·(k-500) 2 ] s2 k:= A·exp[-0.00003·(k-500) 2 ]  информационные сигналы

Q:= 30 q k:= rnd(Q) – Q/2 x1 k:= s1 k + q k x2 k:= s2 k + q k  входные сигналы

Рис. 12.3.1. Модельные сигналы.

В качестве выходных сигналов задаем те же самые функции s1 и s2. Быстрым преобразованием Фурье вычисляем спектры сигналов и формируем спектры плотности мощности:

S1:= CFFT(s1) S2:= CFFT(s2) Q:= CFFT(q)  спектры сигналов

 спектры мощности

Ds1:= var(s1) Ds2:= var(s2) Dx1:= var(x1) Dx2:= var(x2) Dq:= var(q)  дисперсии

Ds1 = 124.308 Ds2 = 310.264 Dx1 = 202.865 Dx2 = 386.78 Dq = 79.038  информация

mean(Wq) = 0.079 Wq1:= (Dx1 – Ds1)/(K+1) Wq1 = 0.078  информация

Wq2:= (Dx2 – Ds2)/(K+1) Wq2 = 0.076  информация

Wq k:= Wq1  замена на постоянное распределение

Для воспроизведения сигналов вычисления функций Wzs не требуется, т.к. Wzs = Ws. Вычисление Wq также имеет только информативный характер.

Передаточные функции оптимальных фильтров (приведены на рис. 12.3.2):

Рис. 12.3.2. Передаточные функции оптимальных фильтров

в сопоставлении с нормированными модулями спектров сигналов

Как следует из рисунка 12.3.2, для плавных монотонных функций, основная энергия которых сосредоточена в низкочастотной части спектра, передаточные функции оптимальных фильтров, по существу, представляют собой низкочастотные сглаживающие фильтры с автоматической подстройкой граничной частоты пропускания под основные частоты входного сигнала. Операторы фильтров можно получить обратным преобразованием Фурье:

h1:= ICFFT(H1)/(K+1) h2:= ICFFT(H2)/(K+!)  обратное преобразование Фурье

Рис. 12.3.3. Импульсные отклики фильтров.

Оператор фильтра, в принципе, бесконечен. В данном случае, при использовании БПФ максимальное число отсчетов равно К/2 = 500. Усечение размеров оператора может выполняться по типовым методам с применением весовых функций (в расчете операторы нормируются к 1, весовые функции не применяются).

N1:= 160 n1:= 0 .. N1 N2 ;= 500 n2:= 0 .. N2  размеры и нумерация операторов

hm1:= h1 0 + 2·h1 n 1 hm1=0.988 h1:= h1/hm1  нормировка

hm2:= h2 0 + 2·h2 n 2 hm2=1.001 h2:= h2/hm2  нормировка

 свертка

Рис. 12.3.4. Проверка действия оптимальных фильтров.

Коэффициент усиления дисперсии помех  Kd:= (h 0) 2 + 2·h n Kd1=0.021 Kd2= 0.0066

Среднеквадратическое отклонение воспроизведения сигнала:

e1= 1.238 e2 = 0.701

Проверка действия оператора фильтра приведена на рис. 12.3.4.

Особую эффективность оптимальный фильтр имеет при очистке от шумов сигналов, имеющих достаточно сложный спектральный состав. Оптимальный фильтр учитывает конфигурацию спектра сигнала и обеспечивает максимальное подавление шумов, в том числе внутри интервала основного частотного диапазона сигнала. Это наглядно видно на рис. 12.3.5 для сигнала, близкого к прямоугольному, спектр которого имеет кроме основной низкочастотной части затухающие боковые осцилляции. Расчет фильтра выполнялся по методике, приведенной в примере 1.

Рис. 12.3.5. Оптимальная фильтрация сигнала со сложным спектральным составом.

Рис. 12.3.6. Оптимальная фильтрация радиоимпульса.

На рис. 12.3.6 приведен пример фильтрации оптимальным фильтром радиоимпульса. Основной пик спектра радиоимпульса находится в области несущей частоты, а боковые полосы определяются формой модулирующего сигнала, в данном случае – прямоугольного импульса. На графике модулей сигнала и передаточной функции фильтра можно видеть, что оптимальный фильтр превратился в полосовой фильтр, при этом учитывается форма боковых полос спектра сигнала.

Фильтры прогнозирования и запаздывания. Если в правой части уравнения (12.3.3) желаемым сигналом задать входной сигнал со сдвигом на величину kt, то при этом B(m) = R(m+k) и уравнение принимает вид:

h(n) ③ R(m-n) = R(m+k). (12.3.8)

При k > 0 фильтр называется фильтром прогнозирования и вычисляет будущие значения сигнала по его предшествующим значениям. При k < 0 фильтр является фильтром запаздывания. Реализация фильтра заключается в решении соответствующих систем линейных уравнений для каждого заданного значения k. Фильтр может использоваться для интерполяции геофизических полей, в том числе в наперед заданные точки, а также для восстановления утраченных данных.

Курсовая работа 18-07. Разработка программы расчета оптимального фильтра приема данных в коде Манчестер-II и фильтрации цифровых данных.

Курсовая работа 19-07. Разработка программы расчета фильтра синхронизации приема данных в коде Манчестер-II и фильтрации цифровых данных.

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

Линейные фильтры являются неотъемлемой частью любого приемного устройства. С их помощью осуществляется как додетекторная, так и последетекторная обработки сигналов. С помощью линейных фильтров сигналы разделяются в многоканальных системах передачи. Требования к этим фильтрам могут быть весьма различными в зависимости от их назначения. Здесь рассмотрим теорию оптимальной линейной фильтрации.

Пусть сигнал на входе линейного фильтра с импульсной реакцией g(t) представляет сумму переданного сигнала s(t) и помехи n(t)

z(t) = s(t) + n(t). (7.57)

Требуется найти такую функцию g(t)y которая минимизирует средний квадрат ошибки

где s(t) -оценка сигнала на выходе фильтра. Здесь считаем, что время запаздывания сигнала s(t) в фильтре t 0 = 0, а среднее значение берется по ансамблям сигналов S и помех N. Будем полагать, что s(t) и n(t) - стационарные взаимно-некоррелированные процессы с известными энергетическими спектрами G s (f) и G N (f). В такой постановке задача была решена независимо друг от друга академиком А. Н. Колмогоровым (1939 г.) и Н. Винером (1942 г.) и поэтому оптимальный (в указанном смысле) линейный фильтр называют фильтром Колмогорова - Винера. Требование физической реализуемости фильтра, как известно, сводится

к тому, что импульсная реакция фильтра должна удовлетворять условию g(t)=0 для всех t

где область интегрирования γ для физически реализуемого фильтра есть интервал (0, ∞), а для не реализуемых фильтров - (-∞, ∞). Можно доказать, что необходимым и достаточным условием оптимальной линейной фильтрации является условие

(7.60)

для всех τ из γ.

Это означает, что фильтр нужно выбрать так, чтобы ошибка ε (t) = s (t) - S(t) была не коррелирована со входным сигналом Z(i) во все моменты времени в области γ. Если бы имела место корреляция между ошибкой и принимаемым сигналом, то при последующей обработке можно было бы получить лучшую оценку.

Докажем справедливость условия (7.60). Пусть g 1 (t) - импульсная характеристика оптимального фильтра, удовлетворяющего условию (7.60), g 2 (t)-импульсная характеристика любого другого линейного фильтра. Отклики фильтров соответственно обозначим через S 1 (t) и S 2 (t). Тогда

Поскольку функция S (t)- Ŝ 1 (t) = ε(t) удовлетворяет условию(7.60), то

Следовательно,

Очевидно, это выражение будет минимальным, когда Ŝ 2 (t) = Ŝ 1 (t), что и доказывает справедливость условия (7.60). Смысл этого условия состоит в том, что случайный вектор Ṡ должен быть строго ортогональной проекцией Ṡ на линейное подпространство, порождаемое случайным вектором z. Представим условие (7.60) в виде

для всех τ из γ. Отсюда с учетом (7.59)

В том случае, когда сигнал S(t) и помеха n(t) некоррелированы,

(7.61) принимает вид

Это основное интегральное уравнение теории линейной фильтрации называется уравнением Винера - Хопфа. Его решением является искомая функция g(t), минимизирующая средний квадрат ошибки

Не следует путать оптимальные линейные фильтры, определяемые (7.61) или (7.62) с согласованными фильтрами, рассмотренными в § 6.4. Если основное назначение рассматриваемых здесь фильтров состоит в наилучшем воспроизведении неизвестной формы сигнала, то задача согласованных фильтров заключается в формировании максимально возможного пика сигнала известной формы в момент отсчета на фоне шума.

Уравнение (7.62) легко решается для нереализуемых фильтров, т. е. когда γ=(-∞, ∞). Для этого случая, применив преобразование Фурье к обоим частям (7.62), получим в частотной области

G s (f) = k(f). (7.63)

Отсюда коэффициент передачи оптимального линейного фильтра

или в более общем случае, когда учитывается время запаздывания t 0 в фильтре,

Ошибка при этом

Легко заметить, что ошибка

только в том случае, когда G S (f)G N (f) = 0, т. е. когда спектры сигнала и помехи не перекрываются. Во всех других случаях оптимальный фильтр пропускает различные частоты с тем меньшим весом, чем больше отношение G N (f)/G S (f) при данной частоте.

Увеличивая время запаздывания t 0 , можно приблизиться к (7.66) и в случае реализуемого фильтра. Задача существенно усложняется, когда требуется, чтобы оптимальный фильтр k(f) был реализуем без существенного запаздывания. Для получения передаточной функции k(f) реализуемого фильтра используют полученное выше решение для γ=(-∞, ∞). С этой целью не реализуемый фильтр (7.65) раскладывают на несколько фильтров и выделяют из него оптимальную реализуемую часть . В общем случае оптимальным по критерию минимума среднеквадратической ошибки является нелинейный фильтр. Исключение представляет случай, когда сигнал и помеха гауссовские, так как для них оптимальный фильтр всегда линеен.

Результаты оптимальной фильтрации можно существенно улучшить, если применить так называемое предыскажение сигнала с последующей его коррекцией на приеме. Сущность метода предыскажения состоит в том, что на передающей стороне сигнал s(t) пропускается через фильтр с коэффициентом передачи k 1 (f). Полученный таким образом видоизмененный сигнал s"(i) передается по каналу. На приемной стороне включен другой фильтр k 2 (f). Характеристики фильтров k 1 (f) и k 2 {f) выбираются так, чтобы обеспечить минимум среднеквадратической ошибки. Расчеты показывают, что предыскажение дает тем больший выигрыш, чем меньше относительная ширина полосы перекрытия спектров сигнала и помехи. Предыскажение позволяет перераспределить мощность полезного сигнала в полосе частот канала так, чтобы обеспечить лучшие условия согласования источника сигнала с каналом (в общем случае полезно стремиться к тому, чтобы сумма спектральных плотностей мощности сигнала и мощности помехи была постоянной в пределах полосы частот канала). Это означает, что предыскажение можно рассматривать как "линейное кодирование" непрерывного сигнала, позволяющее уменьшить ошибку и улучшить использование пропускной способности канала.

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

Следует отметить, что в результате предыскажений формируется новый сигнал с необходимыми свойствами. Так, например, в радиовещании и многоканальной радиорелейной и спутниковой связи с частотной модуляцией несущей используется предыскажение, близкое к дифференцированию. В этом случае на вход частотного модулятора поступает не первичный сигнал b(t), как это делается при ЧМ без искажений, а его производная db/dt. Поэтому пропорционально b(t) изменяется не мгновенная частота, а мгновенная начальная фаза несущего колебания, т. е. формируется не ЧМ, а ФМ сигнал. Так как спектр шума на выходе демодулятора ФМ сигнала равномерный (7.39), то тем самым в многоканальных системах обеспечивается одинаковая помехоустойчивость во всех частотных каналах, а в случае радиовещания - более качественное воспроизведение речевых и музыкальных передач }