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

книги / Твёрдотельная фотоэлектроника. Физические основы-1

.pdf
Скачиваний:
16
Добавлен:
12.11.2023
Размер:
15.42 Mб
Скачать

5.7. Согласованная фильтрация сигнала на фоне помех [15,37]

Как известно, синтез электронных цепей состоит из двух этапов:

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

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

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

обнаружение фотосигнала;

оценка параметров фотосигнала с наилучшей точностью (дальнометрия, на­ ведение);

разрешение или различение сигналов (волоконно-оптические линии связи, тепловидение) и ряд других.

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

5.7.1. Согласованная фильтрация импульсного сигнала при белом шу­ ме. Пусть на вход линейного четырехполюсника с постоянными параметрами

ипередаточной функцией К(и>) подается аддитивная смесь u(t) сигнала s(t)

ишума иш(t) с нормальным законом распределения. Сигнал задан своей спек­ тральной функцией G(ш). Энергетический спектр шума — равномерный (белый шум) W (и>) = WQ.

Пиковое максимальное значение выходного сигнала четырехполюсника до­ стигается в некоторый момент времени tm (о значении tm смотрите ниже) и, согласно (4.5.14),

ОО

 

: {tm) = ~ J G{u)k(uj)exp{jutm)cLj.

(5.7.1)

Среднеквадратичное значение помехи на выходе фильтра (уравнение

(5.5.8))

 

 

 

 

О’вых —

[ W(u))K2(u)dw =

f K 2(UJ) CLJ.

(5.7.2)

\

27г J

\ 2тг

J

 

—оо

)

—ОО

 

Для последующей оптимизации воспользуемся неравенством Шварца:

ь

2

b

2

b

J F1 (X )F 2 (X ) dx

 

p

Г

 

< и

Рг (х) dx x

Ы х ) dx.

 

 

 

J

a

a

Видно что это неравенство обращается в равенство (и, стало быть левая часть достигает максимально возможного значения) при выполнении условия

F2{x) = AFS(x)t

то есть когда ^Ь(х) пропорциональна функции, комплексно сопряженной с Р\ (ж). Здесь А — произвольный постоянный коэффициент.

Приравняв Р\ (ж) = G(ш) и Р2(х) = К (u))exp(ju)tm), квадрат пикового зна­ чения сигнала зВых(*т) согласно уравнению (5.7.1) может быть представлен в

форме

 

 

 

«вых(*т)=

J G (и>) К (ш) exp (jutm) du)

£

 

 

— ОС

 

 

 

ОО

 

ОО

 

< — [

G2(uj)doj—

! K 2{u)duj. (5.7.3)

 

27г J

2п

J

 

— ОО

 

— ОО

В подынтегральной функции последнего интеграла отсутствует экспонен­ циальный член, так как \exp(ju>tm)\2 = 1. Поэтому

2

. 2

к (w)exp (jutm)

= К{ш) |ехр {jojtm)\2 = К (и)

Из формулы (5.7.3) следует, что пиковое значение выходного сигнала огра­ ничивается величиной

ОООО

( tm ) ^

\

а его отношение к среднеквадратичному значению выходного шума (5.7.2) со­ ставляет _________ _______________

^

ОО

ОО

 

/ G2(<jS)dw\ f K 2 ( u ) d u j

 

5 ВЫХ (tm)

—оо

—ОО

 

G

 

оо

 

 

 

% / К 2(ш)с1ш

 

 

 

—оо

 

 

оо

----

 

TWO j G

4 u ) d » = J ± .

(5.7.4)

Здесь ^ J ^ G 2(ш)ёш — полная энергия & входного сигнала согласно урав­ нению (4.1.33). Отношение сигнала к шуму получается наибольшим и равным

\J&/WQ при

К(ш) = К(ш)ехр\jtpk(w)] = AG* (w)exp( - ju>tm) =

= A G (w )e x p { - /[¥ > e(w) + u>tm]}. (5.7.5)

Из последнего соотношения вытекают два условия для фазовой и ампли­

тудной характеристик согласованного фильтра при белом шуме:

 

<Pk(w) = -

[Vs(w) + wim],

(5.7.6)

K(u)

= AG(u>).

(5.7.7)

Первое условие (5.7.6) определяет компенсацию начальных фаз в спектре сигнала. После прохождения через фильтр фазовая характеристика выходного сигнала (—wtm) обеспечивает образование его пика при t = tm. Очевидно, что это возможно не ранее окончания сигнала.

Второе условие (5.7.7) показывает, что амплитудно-частотная характери­ стика К(ш) оптимального фильтра усердно повторяет спектр сигнала G(u>) (с точностью до произвольной постоянной А). Теперь, когда соотношение (5.7.7) выведено, это замечательное правило кажется очевидным: на тех частотах, где спектральная плотность G(a>o) сигнала больше и усиление должно быть боль­ ше. Зачем же усиливать сигнал на частотах, где G(u) —*• 0, ведь вклад этих сигнальных гармоник мал, мы только усилим лишний шум. Напомним, что рассматривается случай белого шума W(UJ) = Wo = const.

Поясним вывод соотношения (5.7.7) более наглядным методом «от простого к сложному», не опираясь на формулу Шварца.

Пусть спектр шума белый, а симметричный во времени сигнал имеет по­ стоянную спектральную плотность G\ в узкой полосе частот Д /с. Если полоса усилителя еще уже (Д /у < Д /с), то отношение сигнал/шум на его выходе

Nc/ui~ U C A

GiKiAfy

> Щ г

л/Тш

у / Щ Ц ,

 

Здесь К\ обозначает коэффициент усиления. NQ/ш увеличивается с ро­ стом A f y до тех пор, пока полоса усилителя не превысит полосу сигнала: при A f y > Д /с сигнал уже не возрастает с увеличением полосы, продолжает рас­ ти только шум. Следовательно, полосу усиления надо выбирать равной полосе сигнала.

Усложним задачу. Дополним спектр сигнала второй частотой /2 так, что он теперь состоит из двух равных узких полос Д /с на частотах /1 и /2 с раз­ ной спектральной плотностью G\ и G2 соответственно (рис. 5.7.16). Рассуждая аналогично первому случаю, находим, что полосы усиления опять должны со­ ответствовать частотным полосам сигнала. Обозначив коэффициент усиления

на частоте /г через К

2 , получим

 

 

_Gi Ki Af y + G2 K2 &f y __G i Ki Af y

1 + в2&2

NQ/UI =

 

 

(5.7.8)

 

V K 'tA f y + K l A f y

V K ' t A f y

у / Г + Щ ’

/1

ed -k i

ec2’

Д .

Строчными буквами ег и /гг обозначе­ ны относительные значения спектральной плотности и коэффициента усиления. Вто­ рой множитель в соотношении (5.7.8) по­

/казывает относительное изменение iVc/ш при использовании обеих полос сигнала. Продифференцировав его по /гг (выпишем только числитель производной), получим

к2

-■ ~ e 2^/l + к % - (1 + е 2/г2)

ТУ

Приравняв производную нулю, найдем, что максимальное отношение достигается при к2 = в2 .

Эта простая выкладка позволяет сде­ лать два важных вывода:

оптимальный фильтр, при котором от­ ношение сигнал/шум получается макси­ мальным, существует;

 

 

 

 

усиление

на частоте /2

оптимально­

 

 

 

 

го фильтра

пропорционально спектраль­

 

 

 

 

ной плотности сигнала на этой частоте.

_ _,

 

В последний раз усложним задачу.

оптимального

фильтра

при белом

шуме Рассмотрим

сигнал, спектр которого со-

повторяет спектр сигнала [37]

держит не две, а п узких равных полос

 

 

 

 

Д/с на частотах f i , f 2 ,---fn

со спектраль­

ными

плотностями

G i,G 2 -..Gn (рис. 5.7.1в). Обобщив уравнение (5.7.8) на

рассматриваемый случай

 

 

 

 

 

 

Д1.

1 + ег/гг + ез/гз +

+ еп кп

 

 

 

 

N c m ~

^

 

 

и приравняв нулю частные производные d N c / u i / d k n , приходим к соотношениям

/гг = сг,

к з = 63,

к п = еп

или

 

(5.7.9)

К\ = Gi,

K2 = G2,- -

K N = GN -

Выражение (5.7.9) практически уже представляет собой частотную характе­ ристику оптимального фильтра при белом шуме и сигнале произвольной формы (рис. 5.7.15).

Определим импульсную характеристику согласованного фильтра. Из соот­ ношений (4.5.19) и (5.7.5) следует

 

 

 

ОО

 

 

 

 

 

 

 

ОО

 

 

h ( t ) =

— / К ( и ) e x p ( j u t ) d u

 

А_ J G * ( u ) e x p \ j u ( t - t m ) ] d w .

 

 

2 п

J

 

 

 

 

 

 

2 7 Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-ОО

 

 

Учитывая, что G * ( u )

= G

( — и ) ,

перейдя к новой переменной и \ = —и

и вспом­

нив уравнение (4.1.22), получим

 

 

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

h ( t )

=

^ ~

[

G ( u i ) e x p

[ j u (

t m

-

t ) ] d u =

^ - s ( t m - t ) .

(5.7.10)

 

 

 

 

Z7T

J

 

 

 

 

 

 

 

Z7T

 

 

 

 

 

—oo

 

 

 

 

 

 

 

 

 

Построение графика функции s

(

t m

-

t ) , где

 

 

t m выбрано равным длительности

сигнала Тс,

 

 

показано

на

рис.

5.7.2.

Импульсный

 

отклик

 

 

h ( t ) фильтра,

согласованного

с

импульсным

 

 

сигналом

s

( t

) ,

есть

зеркально

отраженная и

 

 

смещенная вправо на время сигнала Тс функ­

 

 

ция s ( t ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сигнал на выходе согласованного фильтра

 

 

(см. (4.5.19))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

$ в ы х ( t )

=

^

J

G

( и ) К

( и ) exp ( j u t ) d

u

=

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ^ J

G ( u ) G * ( u ) e x p ( - j u t m ) e x p ( j u t ) d u .

Рис. 5.7.2. Сигнал S i ( t )

и согласо­

 

 

 

 

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

 

(5.7.11)

ванный с ним импульсный отклик

 

 

 

 

 

 

 

 

 

 

 

 

фильтра h { ( t )

 

При t = t m

из (5.7.11)

 

 

 

 

 

 

 

 

 

ОО

«выхМ = ^ J G 2 ( u ) d u = A $ .

—ОО

С учетом того, что G * ( u ) = G ( - u ) , а также соотношений (4.1.19)

( S ( - u ) e x p ( - j u t m ) = F [ s ( t m - t ) } ) и (4.1.27)

оо

(t) = ^- J G(oo) [ G ( - w )e xp ( - ju tm) exp (jcot) dio =

—ОО

оо

= A J s(y)s[y - ( t - t m)]dy = ABa( t - t m). (5.7.12)

Таким образом, сигнал на выходе согласованного фильтра с точностью до

постоянного коэффициента совпадает с корреляционной функцией

входного

сигнала.

 

 

 

 

 

 

 

Вспомним,

что

закон

распределения шума на

выходе

четырехполосни-

ка остается нормальным. Энергетический спектр выходного

шума WBblx (и>) =

= К 2(u>) Wo. а корреляционная функция согласно (5.7.5)

 

 

 

ОО

 

 

 

 

 

 

Д>ых (т ) = ^

J ^вых (и>)ехр Цшт) (ко =

 

 

 

 

— ОО

 

 

 

 

 

 

 

 

 

ОО

 

оо

 

 

 

= %

/

К 2 (ш) exp (jcor) (ко = А2^

J G2 (со) exp (juT) du>.

 

 

 

—оо

 

—оо

 

 

Учитывая (5.7.11) и (5.7.12), в которых (t tm) обозначим как т,

 

Ввых (т) = А2 —

J

G2 (to) exp [jui (t - tm)} (ко =

 

 

 

 

— ОО

 

 

 

 

 

 

 

 

 

= ^W 0SBUx {tm + т) = A2WQBS (т).

(5.7.13)

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

Дисперсия (средняя мощность шума) на выходе

Ввых (0) = A2W0BS(0) = A2W0S.

(5.7.14)

Таким образом, отношение пикового значения сигнала sBUx (tm) к средне­ квадратичному значению шума составляет А §/у7A2WQS = \/&]WQ, ч т о , конеч­ но, совпадает со значением (5.7.4), полученным из частотного анализа.

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

Пусть относительная форма сигнала задана функцией времени

1

при

0 < t < Тс,

(5.7.15)

s(t) =

при

t < О, t > T c.

О

 

Спектральная плотность такого сигнала

 

 

6 И = — [1 -

exp(—jwTc)].

(5.7.16)

Передаточная функция согласованного фильтра согласно (5.7.5)

 

К (w) = А 1 [1 - exp(ju)Tc)\exp{—jioTc) = А ^ - [1 - exp(-ja)Tc)].

(5.7.17)

(-■?<*>)

 

ju

 

Так как К(ш) отличается от спектра сигнала G(u>) только постоянным коэффициентом, то импульсная характеристика согласованного фильтра h(t), определяемая соотношением (4.5.19), совпадает по форме с сигналом

h(t) = As(T{

ь)

 

(5.7.18)

С

 

при f < 0 ,

 

 

Структурная схема фильтра, согласованного с прямоугольным сигналом, имеет вид, представленный на рис. 5.7.3. Первый множитель 1 /jui реализуется интегрирующим звеном, а второй [1 —ехр(—juiTc)) — устройством вычитания, к которому сигнал попадает без задержки и с задержкой Тс (передаточная функция линии задержки без потерь равна ехр(—ju>Tc)).

Корреляционная функция прямоугольного импульса (5.7.15) имеет вид рав­ нобедренного треугольника с основанием Тс и высотой, равной относительной энергии S импульса, выде-

ляемои

на

единичной

на­

 

грузке

(<£ =

1 Тс = Тс).

В

 

соответствии с (5.7.12) мак­

 

симум

выходного

сигнала

 

АТС достигается при t = ТС

Р и с. 5.7.3. Блок-схема согласованного фильтра для пря­

(к концу действия

входного

моугольного сигнала и белого спектра шума

сигнала). При этом отношение сигнала к среднеквадратичному значению шума согласно (5.7.4) равно \J&/WQ = \/TC/W Q.

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

На практике используют квазиоптимальные фильтры из более простых ком­ понентов, подбирая их параметры по тому же критерию — максимальному от­ ношению сигнала к шуму. Часто в качестве квазиоптимального фильтра ис­ пользуют одно или несколько #С-звеньев.

Рассмотрим прохождение прямоугольного видеосигнала через интегрирую­ щую ЛС-цепочку с передаточной функцией К (и>) = 1/(1 + jujRC).

Максимальное значение входного сигнала в этом случае, достигаемое в момент tm = Тс, составляет

 

 

звых (Тс) =

1 -

exp (- % / R C ),

(5.7.19)

а среднеквадратичное напряжение шума

 

 

 

 

 

[ W0K 2( u ) d u =

 

/

-

du>

W0

^ВЫХ--

\

9

(5.7.20)

2тг

J

2тг

J

1 +

iRC)2

2RC'

\

 

—оо

\

—оо

 

 

 

В результате

 

 

 

 

 

 

 

 

SBHX(Тс) _ 1

ехр( Tc/RC)

Ik

 

2RC

1 —exp

(5.7.21)

&вых

 

\JWQ/2RC

W0

~тГ

 

 

При Tc/R C = 1,28 отношение сигнала к шуму достигает максимального значения, равного 0,9y/Tc/Wo. Таким образом, переход от согласованного с сигналом фильтра к квазиоптимальному однозвенному ЯС-фильтру приводит к уменьшению отношения сигнала к шуму всего на 10% При этом отношение сигнала к шуму оказывается некритичным и к настройке ЯС-цепочки: измене­ ние Tc/R C в 4 раза (от 0,64 до 2,56) приводит к дополнительной потере всего

в8%.

5.7.3.Согласованная фильтрация сигнала при небелом шуме. Пере­ даточная функция к (ш) согласованного фильтра определяется в результате приведения заданного спектра шума к белому.

Рассмотрим

структурную схему условного

фильтра,

показанную

на

рис. 5.7.4.

На этой схеме

К {и)

обозначает

искомую

функцию

 

 

I

1

синтезируемого

фильтра, а

 

 

к\{ш)

И \/к \(ш )

являют­

 

 

 

 

 

 

 

 

ся

передаточными

функция­

 

 

 

 

ми двух условных четырехпо­

 

 

I

1

люсников. Введение этих че-

 

 

при тырехполюсников,

очевидно,

Р и с .

5.7.4. Согласованная фильтрация

сигнала

небелом шуме

 

 

не оказывает влияния на ра­

боту устройства так как их результирующая передаточная функция равна единице.

Зададим модуль функции К\ (и>) в виде

=

< 5 7 '2 2 )

где W (ш) — энергетический спектр шума на входе фильтра, WQ постоянная величина. Теперь на выходе первичного четырехполюсника наблюдается шум с равномерным энергетическим спектром (белый шум)

Wi(bj) = W ( w) Kt (w) = Wo-

Очевидно, что и сигнал на выходе этого четырехполюсника отличается от входного сигнала. Его спектральная плотность

G i H = G ( w ) X i M .

(5.7.23)

Однако из раздела 5.7.1 следует, что для максимизации отношения сигнала к шуму сохранение формы сигнала не требуется.

Так как на выходе первого четырехполюсника шум белый, то для получения на выходе всего фильтра максимального отношения сигнала к шуму передаточ­ ная функция фильтра согласно (5.7.5) должна быть равна

~ \ К (ш) = AG1 (oj)exp(-ju)tm) .

Kiiuj)

Откуда

К(и>) = AG\{UJ) K I (w)exp(—jujtm) = AG* {ш)К{ (u>)exp(-ju>tm) =

= AG*(u)KUu)exp(-ju;tm) = AWQ^~— exp{—j u t m) . (5.7.24)

Физический смысл этого соотношения очевиден. Как и в случае белого шу­ ма, для получения максимального отношения сигнала к шуму фильтр должен обеспечить компенсацию начальных фаз спектра входного сигнала G(u). По­ этому в (5.7.24) входит комплексно-сопряженная функция G* (ш). Однако для усиления тех компонентов спектра сигнала, для которых интенсивность шу­ ма меньше, модуль передаточной функции не только пропорционален модулю G(UJ) (как при белом шуме), но и обратно пропорционален энергетическому спектру шума на входе фильтра.

5.8.Минимально обнаруживаемые оптические сигналы

5.8.1.Ограничение флуктуациями сигнала. Фундаментальное огра­ ничение минимально обнаруживаемой оптической энергии обусловливается квантовой природой излучения: в импульсе оптического сигнала не может быть меньше одного фотона. В идеальном фотоприемнике — счетчике фотонов один фотон преобразуется в один электрон, при этом минимальная площадь

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

р

_ и

he

q

1,24 <7 _ 2

10 - 1 9

© m in

 

^

A ”

A

 

 

 

 

Здесь размерности [£т т] = Джоуль; [А] = микрон; Si — монохроматическая токовая чувствительность идеального квантового фотоприемника (Si = q/hv).

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

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

Случайное отклонение числа фотонов в импульсе от своего среднего зна­

чения N

(и,

соответстенно, флуктуация энергии импульса) характеризует­

ся дисперсией

A N 2 = (N N ) 2 Для пуассоновского процесса, как известно,

A N 2 = N

При уменьшении N все чаще появляются импульсы с малым числом

квантов и вероятность пропуска импульсов возрастает (малый сигнал как бы сам мешает себя обнаружить).

Предположим, что фотоны регистрируются идеальным фотоприемником без потерь, при этом не создается дополнительных шумов. Рассмотрим одиночные импульсы со средним числом фотонов в импульсе N Из статистики Пуассона следует, что вероятность реализации импульса с числом фотонов N составляет

 

е х р ( - « К

v

т

В этом случае естественное правило принятия решения «сигнал есть» — это когда фиксируется хотя бы один электрон. Тогда вероятность пропуска сигнала