Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

МАОС без 1 главы

.pdf
Скачиваний:
47
Добавлен:
22.03.2016
Размер:
1.43 Mб
Скачать

ние равно нулю, и в этом случае ковариационная функция совпадает с корреляционной:

K( ) = R( ).

Таким образом, согласованный фильтр в принципе выполняет ту же операцию, что и корреляционный приемник; в этом смысле они эквивалентны. Вопрос о применении корреляционного приемника или согласованного фильтра в каждом конкретном случае решается в зависимости от простоты технической реализации. Следует отметить, что при одном и том же входном сигнале s(t, ) характер сигнальной и шумовой функций на выходе корреляционного приемника и согласованного фильтра различен. Отличие в вычислении корреляционного интеграла корреляционным приемником и согласованным фильтром заключается в том, что коррелятор определяет единственную точку ковариационной функции сигнала, а именно максимальную точку, а согласованный фильтр вычисляет ковариационную функцию сигнала полностью.

Дисперсия шума на выходе согласованного фильтра может быть определена через выходной шум, равный на основании (1.12) и (1.13б):

 

 

 

 

 

 

 

 

 

 

kN

 

 

 

 

 

 

 

 

 

 

 

 

 

 

nв(t)

 

0

q0n(t).

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда корреляционная функция шума

 

 

 

 

 

R(t

t

2

) M n

t

n t

2

 

 

 

 

 

 

 

 

1

 

 

 

в

1

в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 2M

 

s(t0 t1 u1,λ)n(u1)du1

s(t0 t2 u2,λ)n(u2 )du2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 2

 

s(t0

t1 u1,λ)s(t0 t2 u2 ,λ)n(u1 )n(u2 )du1du2

 

 

.

Но n(u1)n(u2 ) есть корреляционная функция белого шума:

n(u1)n(u2 ) N20 (u1 u2 ) .

Тогда в силу фильтрующего свойства -функции интеграл отличен от нуля лишь при u1 = u2 = u и двойной интеграл переходит в однократный:

 

 

k 2 N0

 

R(t1

t2 )

s(t0 t1 u,λ)s(t0 t2 u,λ)du ,

 

 

2

 

 

 

 

или, вводя замену t0 – t1 + u ,

13

 

 

2

N

 

 

 

 

R(t1

t2 )

k

0

 

s( ,λ)s(τ t1

t2 ,λ)dτ.

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.15)

Из формулы (1.15) видно, что корреляционная функция выходного шума имеет вид корреляционной функции входного сигнала. Тогда дисперсия выходного шума

Dв= R(t

t )

k 2N 0Es

.

(1.16)

 

1

1

2

 

 

 

 

 

 

Из формул (1.14) и (1.16) находим отношение наибольшего значения выходного сигнала к среднеквадратическому значению выходного шума:

q

s

в max

(t

0

)

 

2E

s

,

 

 

 

 

 

 

D

 

 

 

N

 

 

 

 

 

 

 

 

0

 

 

 

в

 

 

 

 

 

где q – отношение сигнал/помеха по амплитуде.

1. 5. Физически возможные фильтры. Квазиоптимальные фильтры

(1.17)

При практическом построении оптимальных и согласованных линейных фильтров кроме найденных соотношений надо также учитывать условия физической возможности и практической реализуемости фильтров. Условие физической возможности фильтра записывается в виде [6]:

h(t) = 0 при t 0;

 

log K jω

 

 

2

dω .

 

1 ω

 

 

 

Если сигнал s(t), с которым должен быть согласован фильтр, начинается в момент времени 0 и полностью прекращается при t 0+ и, то первое из условий выполняется при t0 0+ и. Только при этом условии будет использована вся энергия сигнала для формирования сигнального пика на выходе фильтра в момент t0. Увеличение t0 сверх 0 + и, не влияя на значение пика, сдвигает его в сторону большего запаздывания, что обычно нежелательно. Поэтому следует брать t0 = 0+ и, т. е. момент наблюдения должен совпадать с окончанием входного сигнала. Иногда для аппроксимации реальных импульсных сигналов используют бесконечно длинные импульсы (гауссовский, экспоненциальный и т. д.). Тогда приходится искусственно выбирать конечное значение длительности аппроксимирующего сигнала, содержащей основную долю энергии реального сигнала.

14

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

Обозначим через отношение значения сигнал / помеха на выходе произвольного линейного фильтра к значению сигнал / помеха на выходе согласованного фильтра. Используя выражение (1.6) и заменяя Sn( ) на N0 / 2 (для белого шума), получаем

 

 

 

 

 

j t

 

 

 

 

 

 

Ф j K j e

0

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

1/ 2

 

 

 

 

 

 

 

 

 

 

Ф j

2

K j

2

 

 

 

d

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В таблице приведены максимальные значения max для различных форм полезных радиоимпульсных сигналов и разных видов частотных

Радиоимпульс

 

Вид АЧХ фильтра

max

a

Прямоугольный

Идеально прямоугольная

0,91

1,37

Гауссовский

Идеально прямоугольная

0,94

0,72

Прямоугольный

Гауссовская

0,94

0,72

Гауссовский

То же

1,00

0,63

Прямоугольный

Одиночный резонансный контур

0,90

0,40

То же

2

несвязанных резонансных контура

0,93

0,61

- “ -

5

несвязанных резонансных контуров

0,94

0,67

характеристик реализуемых фильтров при наилучших значениях их полос пропускания. При этом полоса выбирается из условия и = a, f – ширина полосы пропускания на уровне 0,5 по мощности, и – эффективная длительность импульса.

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

При проектировании квазиоптимальных фильтров задаются структурой фильтра исходя из конструктивных соображений, а полосу его пропускания на уровне 0,707 от максимума определяют, максимизируя величину ρ при

15

изменении полосы пропускания. Квазиоптимальные фильтры для радиоимпульсных и вообще для высокочастотных сигналов выполняются на базе колебательных контуров или активных полосовых фильтров. Число контуров обычно задается из конструктивных соображений. Полосу пропускания оптимизируют, изменяя добротность колебательной системы.

На практике часто приходится работать с сигналами, имеющими случайную амплитуду и фазу. Как следует из (1.8) и (1.9), форма частотной характеристики не зависит от амплитуды. Поэтому для сигнала со случайной амплитудой можно использовать тот же фильтр, что и для сигнала с детерминированной амплитудой. Фазочастотная характеристика фильтра зависит от фазы сигнала. Однако при непрерывном случайном изменении фазы сигнала мы в подавляющем большинстве не имеем возможности перестраивать фильтр. Поэтому случайная фаза сигнала при проектировании фильтра принимается равной своему среднему значению, что несколько снижает отношение сигнал/помеха на выходе.

1.6. Синтез оптимальных фильтров

Рассмотрим различные способы синтеза оптимальных фильтров. Фильтры для выделения сигнала на фоне коррелированного шума строятся обычно на основе спектрального метода, т. е. при использовании для комплексной частотной характеристики фильтра выражения (1.8).

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

16

1.6.1. Синтез согласованного фильтра для прямоугольного видеоимпульса

Рассмотрим временной и спектральный методы синтеза фильтра на примере прямоугольного видеоимпульса:

A

при

0 t и

 

;

x(t)

0

при

t < 0, t

 

.

 

и

 

 

 

 

Последовательность действий при синтезе временным методом иллюстрируется рис. 1.2.

x1(t)

(t)

t

x2(t)

1

t

x3(t)

1

x4(t)

τи

t

 

 

 

τи

t

 

 

-1

x5(t)

1

τи

t

 

Рис. 2.2

Известно, что единичная ступенька (перепад), или функция Хевисайда Y(t), (x2(t)) есть интеграл от -функции (x1(t)):

17

Y(t) x2 (t) (t)dt 1. 0

После задержки единичной ступеньки на длительность импульса и (x3(t)), ее инвертирования (x4(t)) и вычитания из x2(t) получим заданный прямоугольный импульс x5(t), амплитуду которого можно изменять, меняя коэффициент передачи устройства. Отсюда следует, что искомый фильтр (рис. 1.3, а) состоит из интегратора 1, линии задержки на и 2, инвертора 3, сумматора 4 и усилителя 5. Инвертор и сумматор могут быть заменены вычитающим устройством 6 (рис. 1.3, б). Работа обоих вариантов фильтра идентична.

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

3

-1

 

 

 

 

 

x5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

x4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

1

 

 

 

 

 

 

 

6

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

x5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и x3

б

Рис. 1.3

Рассмотрим синтез фильтра спектральным методом. Комплексный спектр прямоугольного видеоимпульса

 

 

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jωt

 

 

 

и

 

 

 

jωt

 

 

 

 

A

 

 

jωτ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(jω)

 

x(t)e

 

 

dt

Аe

 

 

 

 

dt

 

1 e

 

и

..

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

jω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для согласованного фильтра

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jωt

 

 

 

 

A

 

 

 

 

 

 

jωτ

jωt

 

 

 

 

 

 

K0(jω) k *(jω)e

 

 

0 k

 

 

 

 

 

 

1

e

 

 

 

и e

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

jω

 

 

 

 

e

 

 

 

 

 

 

 

 

.

 

 

A

jωτ

 

 

jωt

 

 

 

 

 

A

 

jωτ t

 

 

jωt

 

k

 

 

 

и

 

1 e

 

 

 

 

0

k

 

 

 

 

и 0

 

e

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jω

 

 

 

 

 

 

 

 

 

 

 

 

 

jω

 

 

 

 

 

 

 

 

 

Полагая и = t0 , окончательно получим

18

K

0

( j )

kA

(1

e

j

и ).

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим члены, входящие в выражение (1.18). Оператор

(1.18)

1

, как из-

j

 

вестно, представляет собой оператор идеального интегрирования гармониче-

ского сигнала; kA – коэффициент передачи линейного устройства; e j и – задержку на время и. Видно, что структурная схема такого фильтра соответствует рис. 1.3.

Механизм работы согласованного фильтра (см. рис. 1.3) можно выяснить, рассматривая прохождение через него импульса сигнала и шума

(рис. 1.4).

x1(t)

x2(t)

 

 

 

 

t

 

u

 

 

 

 

 

x3(t)

 

t

 

 

t

x4(t)

 

 

 

 

u

2 u

t

Рис. 1.4

Входной сигнал x1(t) с помощью интегрирующего устройства накапливается в течение времени и = t0 до своего пикового значения (x2 на рис. 1.4). Задержка на и (x3) и вычитание прекращают накопление сигнала, который

19

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

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

Для прямоугольного видеоимпульса ковариационная функция

K

 

и

A dt

 

s(t)s(t )dt

 

 

 

2

 

 

 

 

 

 

 

и

2

dt

A

 

 

A

(

 

2

 

 

 

 

и

)

,

или, в силу симметрии ковариационной функции:

 

2

 

 

 

 

K ( ) A

 

1

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

и

 

 

 

 

 

  

.

Заменяя в последнем выражении на t t0 и полагая и = t0 , на основании формулы (1.14) получаем

2

 

 

t

и

 

 

 

 

 

sв (t) kK(t t0 ) kA

 

 

 

и 1

и

 

.

 

 

 

 

 

График этой функции совпадает с приведенным на рис. 1.4. Видно, что длительность полезного сигнала на выходе фильтра удваивается.

1.6.2. Синтез оптимального фильтра для приема прямоугольного импульса на фоне коррелированного шума

Пусть шум на входе фильтра имеет спектральную плотность Sn( ), отличную от равномерной:

Sn(ω)

2ag 2

 

,

2

g

2

 

ω

 

 

где 2a – спектральная плотность шума при = 0; g–- постоянная, характеризующая ширину энергетического спектра.

В соответствии с выражением (1.8) тогда можно получить формулу для комплексной частотной характеристики оптимального фильтра для приема прямоугольного видеоимпульса на фоне коррелированного шума:

 

 

*(jω)

 

jωt

 

 

 

 

A

e

jωτ

1

2

g

2

jωt

 

 

K(jω) c

e

 

c

 

ω

e

 

 

 

 

 

 

 

 

0

 

 

и

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

(ω)

 

 

 

 

 

 

jω

 

 

 

 

2ag

2

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cA

 

 

(1 e

jωτ

)

1

g

2

 

2

.

 

 

 

 

 

 

 

 

 

 

 

 

и

 

 

jω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ag

2

 

 

 

 

 

 

jω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.19)

20

Оператор j соответствует оператору идеального дифференцирования. Структурная схема оптимального фильтра, построенного в соответствии с формулой (1.19), изображена на рис. 1.5:

x(t)

1

 

cA

 

 

2

 

 

 

 

 

 

 

 

2ag2

 

 

 

 

 

 

 

 

 

 

8

7

3

5

g2

d dt

6

4

и

d dt

Рис. 1.5

На рисунке 6,7 – устройства дифференцирования; 8 – вычитатель. Назначение остальных блоков ясно из предыдущего.

Выражение (1.19) может быть преобразовано к виду

 

cA

 

1 e

j

 

 

g

2

 

 

 

K ( j )

 

и

 

 

j

 

.

 

 

 

 

 

j

 

 

2ag

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

1.6.3. Синтез фильтров, согласованных с радиоимпульсом

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

21

z

a

(t)

 

 

Его комплексный спектр

   

Ae

j

0

t

,

0 t

и

;

 

 

 

 

 

 

 

 

 

0

вне интервала .

 

 

 

 

и

 

 

 

 

 

 

 

τ

 

 

 

(jω) za(t)e jωt dt A e j

ω0t e j ωt dt

 

 

 

 

0

 

 

 

 

τ

 

 

 

 

j(ω ω)τ

и

j(ω

ω)t

 

e

A e

dt A

 

0

и 1

0

 

 

 

j(ω ω)

 

0

0

 

Тогда комплексная частотная характеристика согласованного фильтра при t0 = и:

K (jω) k *(j ω)e

j ωt

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jωτ e

j(ω ω )τ

1

 

e

j ω τ

e

j ωτ

 

 

kAe

 

 

0

и

kA

0 и

 

и

.

 

 

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j(ω ω )

 

 

j(ω ω )

 

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

Если 0 и = 2n , где

число периодов, то

e

jω τ

0

и

 

 

n–- целое

e

jω τ

 

0

и

 

 

число, т. е. импульс содержит целое

1

и

 

1 e

jωτ

K (jω) kA

и

 

 

 

0

j(ω ω )

 

 

0

(1.20)

Множитель j соответствует передаточной функции высокодобротного колебательного контура (при сопротивлении потерь, равном нулю). Сравнивая выражения (1.18) и (1.20), можно заметить, что при переходе от приема прямоугольного видеоимпульса к приему радиоимпульса с прямоугольной огибающей необходимо идеальный интегратор заменить на высокодобротный колебательный контур. Структурная схема согласованного фильтра для радиоимпульса с прямоугольной огибающей и формы входного и выходного сигналов показаны на рис. 1.6. Из него видно, что согласованный фильтр для радиоимпульса искажает входной сигнал, затягивает его, причем форма огибающей совпадает с формой сигнала на выходе фильтра, согласованного с прямоугольным видеоимпульсом.

Найдем выражение для выходного сигнала такого согласованного фильтра. При этом для упрощения вычислений используем соотношение [7]:

 

1

 

 

 

 

K(τ)

Re

za(t) z*a(t τ)dt

,

 

 

2

 

 

 

 

 

 

 

 

где za(t) – аналитический сигнал; Re - реальная часть. Подставляя сюда значение za(t), приведенное выше, получим:

22