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

книги / Теория и расчет импульсных и цифровых оптико-электронных систем

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

Аналогичным образом для конуса с

углом при

вершине 2а

и высотой Як, образованного вращением прямой

х = (tg а) г

вокруг оси

oz, найдем

 

 

и

F (z) = 2яр (tg3a)z{(1+

tg2а)1/2

(2.34)

£э (0 = яре2(tg3а)//[2 (1 +

tg2а)1/2],

(2.35)

 

где

0 < t < ть тЛ= 2Я;г/с.

(2.36)

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

мостью г (у) = а + b cos" у, то

в

этом случае

выражение для

& (0 будет

 

 

 

Ум <*>

 

Ь|dxjdт|

 

|дх/дт|

+

 

go (0 —

1+п •II dy. (2.37)

р2+ ?)'12

 

(1 + Р2 + <72) 2

J lt=*

ym(t) И

 

 

 

 

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

В частности, для сферического объекта при п = 1 на основе выражения (2.37) получим

^ ( / ) = w

(Tc _ ° 2 [ а + i b (Tc“ / ) ]-

(2-38)

По известным импульсным характеристикам объектов могут быть рассчитаны их частотные характеристики и ЭПО.

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

1. Сфера.

Кэ (/со) =

- ~ГзС 12 (“ V — sin юте) + /

12 0 — cos WTc) — «от?]);

 

 

 

(2.39)

 

Qs *=* 2ряЯ2/3.

 

(2.40)

Кзо(w) =

3 {4 (соте — sin юте)2+ [ 2 ( 1 — cos штс) — са2т212)1/2.

(2.41)

2. Конус.

 

 

 

2npRl sin a

,. .

 

Кь (M — -------ъ~2 * [(coscoTft —

0)T/t sin m k — 1) +

 

 

CO-T*

 

 

 

+ / (сот* cos o)T/t — sin (0Tft)],

(2.42)

31

где Rk — радиус основания конуса; 2а — угол при вершине;

<2э= Jt/?£psina;

(2.43)

Кэо (и) = [2/(шМ)] [(1 - cos сот*)2+ (сох* - sin <от*)2]1/2.

(2.44)

Формулы (2.40) и (2.43), выражающие известные результаты относительно ЭПО, получены как частные случаи более общих

соотношений.

На рис. 2.5 изображены нормированные амплитудно-частот­ ные характеристики и соответствующие им импульсные харак­ теристики сферы и конуса, при­ чем принято тс = тк = т. При­ веденные кривые характери­

\^

 

 

 

 

 

 

 

 

зуют

свойства

объектов

как

0,8

 

 

 

 

 

 

 

линейных

фильтров.

При оди­

0,7

\

 

 

/ - 4

 

наковой

радиальной

протяжен­

0,6

 

KJ \

 

ности

объектов

их

спектраль­

0,5 -

 

 

 

у

/

 

 

ные функции по своей ширине

Ofi

 

 

 

 

 

 

 

(в общем

случае)

могут

раз­

 

 

 

 

 

 

 

личаться

 

достаточно

суще­

0,3-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ственно.

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— u

По известным

характеристи­

0,1

 

 

 

 

 

 

 

11

з к

5 t|

 

- о

кам gэ (/)

и

/Сэ (/©) объекта

/ 1 f

 

 

0 f

 

:

-

П

и L «i

 

 

 

 

 

 

 

 

 

 

______ !_______ 1______ 1_______ 1_______1_______ i

и заданному зондирующему сиг­

0

0,1 0,г 0,3 ф 0,5 0,6

0,7 0,8 0,9 1,0

2/г

Рис.

2.5.

Нормированные амплитуд­

налу

могут

быть

определены

характеристики

 

отраженного

но-частотные

характеристики

(1

сигнала во временной и частот­

сфера;

2 — конус)

и

соответствую­

щие

им

 

импульсные

 

характеристики

ной областях.

 

 

 

что нахо­

 

 

(3— сфера; 4 — конус)

 

Следует

отметить,

 

 

 

 

 

 

 

 

ждение

временных

характери­

стик отраженного сигнала посредством интеграла Дюамеля

имеет

в нашем случае свои особенности, связанные с тем, что функ­

ции & ( / ) — финитные. Это предопределяет

необходимость раз-

 

 

 

 

 

Таблица 2.1

 

It

ВрсменнбП

интервал

 

 

 

. 0 а ^

Ти

0

t

То

т и

ТИ 5 $

Т0

0

т и

 

 

Т0 г $ ^

Т1{ + Т0

t — Т0

Ти

 

 

 

 

/

 

 

 

0 а $

Т0

0

t

То ^

т и

. Т0 ==5

Ти

< — Т0

t

 

 

т и ^ ^ т и + т 0

< — То

Ти

32

дельного вычисления интеграла для соответствующих временных интервалов

4 (0= 4" j Еэ(т) ё* V — т) dT>

где значения пределов tHи £в зависят от соотношения между дли­ тельностью излучаемого импульса ти и временной протяженностью объекта т0 (табл. 2.1).

§ 2.3. СИСТЕМНЫЕ ОТРАЖАТЕЛЬНЫЕ ХАРАКТЕРИСТИКИ ОБЪЕКТОВ

Рассмотренные выше характеристики Qa, /Сэ (/©) и g:) (t) позво­ ляют определить отраженный сигнал в эталонных условиях,

предусматривающих равномерность поля излучения

и приема

в рабочей зоне, а также наблюдение с больших по

сравнению

с габаритными размерами объекта дистанций. При этом харак­ теристики <3Э, Ка(/со) и ga (t) зависят лишь, от свойств самих объ­ ектов и поэтому могут считаться характеристиками собственнообъектов.

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

Таким образом, можно говорить о системных отражательных характеристиках, зависящих от параметров ОЭС и условий наблюдения.

Пусть излучающее и приемное устройства ОЭС характери­

зуются

диаграммами

направленности

(Р*, ру) и

ЧС (p.v, Р,,),

где рх,

ру — углы,

отсчитываемые от

совпадающих

осей полей

излучения и приема по двум взаимно перпендикулярным направ­ лениям. Положим также, что система координат объекта связана с некоторой его характерной точкой (в качестве которой могут быть выбраны, например, наиболее близкая точка поверхности, геометрический или энергетический центр объекта и т. п.), а опти­ ческая ось ОЭС ориентирована на произвольную точку поверх­ ности объекта [х0, у0]. Тогда при сохранении исходной предпо­ сылки о существенном превышении дистанции над габаритными1 размерами объекта выражение для системной ЭПО объекта может быть записано в виде

QJ = j j

(Р* - Д«. Р» - Д») Ч'з (Р* - Д*. Р» - Д»)г (у) cos’ у ds,

 

(2.45>

2 Лебедько £. Г,

за

где Дх = x0)/L0; Д у = ( у У

о ) / Ь 0 \ s* — облучаемая

область

поверхности объекта.

сигналов от объекта и

эталон-

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

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

Следует заметить, что при заданном расположении объекта различным ориентациям оси ОЭС (характеризуемым углами А*, Ау) будут соответствовать разные значения Q3, что, естественно, осложняет процедуру расчета или экспериментальных иссле­ дований.

С учетом формулы (2.45) системная импульсная характеристика объекта получит вид

gi (0 = JJ (Р* - Д*. P# ~ &у) Ч'г (Р* - Д*.

—Д») X

S*

 

X г (у) cos2 уб (t — 2zlc) ds.

(2.46)

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

(Р) = ехр

 

(2.47)

где ©!, 02 — ширина поля

излучения и приема

на уровне~0,5

от максимального значения.

углом при вер­

1. Найдем значения Q3

и g£ (t) для конуса с

шине 2а при его центральном облучении (совпадение оси опти­ ческого поля ОЭС с осью конуса).

Сучетом осесимметричного характера задачи (Ах = Ау s= 0)

представим выражение

для

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

ОЭС

в виде

 

 

 

(р7е?) (1 + е?/0|)).

 

 

У (Р) = % ф) ъ ((?) =

ехр { - я

 

(2.48)

Учитывая,

что для

произвольной

точки поверхности

конуса

р = z tg a/L0,

запишем

на основании выражений (2.45)

и

(2.48)

Q; = f f e",,z’cos2yds =

 

(2.49)

\ e** f (г) <fe,

 

 

s

 

 

0

 

 

34

где 2м = # ь g2 =

£зх tg2 ot/(Z,oG?)] О +

0?/0i)»

a F (г)

определяется,

согласно выражению (2.34). Тогда

 

 

 

 

я,

 

 

 

 

2яр tg!>ос

ze-9'2* dz =

 

 

(1 -h tg2^ 1'2 /

 

 

о

 

 

 

=

р sin aLo0? (1 + Oi/ef)'1(1 -

е"’^ ) .

(2.50)

Обозначим через Dx = L0Q1 раствор поля излучения в кар­ тинной плоскости объекта. Тогда с учетом Hh = RJ\g а получим

 

,2

Q g ^ p s in a D H i+ e i/e i)'1

nR^

- e x p

 

D2.

( i + e ? / e i )

(2.51)

При 0f/02 < 1, что типично для некоторых ОЭС активного типа, имеем

Ol = рD\ sin a [1 — exp (—siRl/Df)].

(2.52)

При согласованных полях излучения и приема (0Х= 02)

имеем

 

 

Q* =

(1/2) pD?sin a[ 1— exp (—2nRu/Di)]-

(2.53)

Например, при Dx — 2Rh. из последней формулы

получим

Ql =

Rl sin a (l — е"я/2) = l,58p#| sin a.

(2.54)

Сравнивая выражения (2.54) и (2.43), заключаем, что 2. Таким образом, в рассматриваемом случае неравномерность

диаграммы направленности приводит к уменьшению отраженного сигнала в два раза. Аналогичным образом можно показать, что

при Di =

2Rh и 0i =

02 отношение <2э/<2э =

1,44.

Отметим также,

что. при

01 = 02 -*

оо (Dx -> оо) имеем

Ql =

Q3 — Jip/?|sina,

т. е. получаем выражение (2.43) как частный случай более общей формулы.

Системная импульсная характеристика конуса с учетом выра­ жений (2.46) и (2.48) примет вид

 

 

 

 

 

 

(2.55)

где

F* (z) = er^'F (2),

0 < / <• (2Hh/c) =

%k.

 

 

С учетом

соотношения

(2.33) получим

 

 

gt (ё) = ёэ V) exp

 

 

(2.56)

где

g„ (t) определяется соотношением

(2.35).

 

 

Например,

при

8. >

0Х и Dx =

2Rk = 2Нкtg а

 

 

 

gl (t) = g, (/) exp [—я Л 2/(16//?)].

(2.57>

 

2. Случай наклонного облучения неограниченной плоской

поверхности.

 

 

 

 

 

2*

35

Имея в виду узкопольные ОЭС, будем считать, что поверх­ ность облучается плоской волной. Пусть оптическая ось поля совпадает с ог, а поверхность проходит через начало координат под углом у к плоскости хог. Тогда системная ЭПО плоской по­ верхности будет

Qs = COS v -fJ-oo -f-oo Y,(p)4J \§)dxdy,

что при гауссовых диаграммах направленности приводит к вы­ ражению

 

=

cos у \

ехРJ [ ~ Л (Р* +

Ру)

+

(2.58)

 

 

 

•^-ОО—00

 

 

 

 

Поскольку

 

= xlL0i

ру = y/L0,

где L0 — осевая

дальность

до поверхности,

получим

 

 

 

 

 

 

-f-OO

-f-OO

р

*2 + У1

 

 

cos у

QZ=

cos у

J

 

е х|р [ - п

dxdy =

1 +

 

 

(0i/e2)2*

 

 

— 00

— 00

 

 

 

 

(2.59)

В

частности,

при 02

0i, Ql =

LoGfcosy.

 

 

 

Для системной импульсной характеристики плоской поверх­ ности при тех же условиях имеем

-f-oo -f-oo

 

 

 

gZ(t)=*. cosy j

expj

—я

*2 + ya

! ) ] « ( - - j - ) dxdy.

 

1 +

•—oo

—oo

 

Ф 1

(2.60)

 

 

 

 

Для упрощения вычислений положим, что рассматриваемая поверхность проходит через ось оу. Тогда z — х tg у и в соответ­ ствии с этим

 

-f-00

-f-OO

р

 

x* + y*

 

gl (t) = cos у

 

 

exp I

—n

 

 

 

 

 

’ Jdy j

Ф ,2

 

 

—00

—00

 

 

 

 

c cos у

 

яу2

dy X

 

 

2 tg у

 

 

Ф 2

 

 

 

 

 

ЯС2Т2

x

4 t g 2 yL\b\

 

c cos Y^O®I

(2.61)

 

где т4 = 2/.09itg T/{с / 1 + в?/в1).

56

Таким образом, при гауссовых диаграммах направленности импульсная характеристика наклонно облучаемой плоской лам­ бертовой поверхности имеет гауссову форму с длительностью тб ло уровню ^ 0,5 от максимального значения.

В частности, при 0а > 0i и тв = 2Ь$Хtg у/с нормированная импульсная характеристика плоской поверхности имеет вид

glo (0 = (1/^а) ехр [—я т ) ] .

(2.62)

При относительно широком поле излучения аппроксимация облучения плоской волной может оказаться неправомерной. Расчет импульсных характеристик в таких условиях рассмотрен

вработе [49].

§2.4. НЕКОТОРЫЕ

ЭКСПЕРИМЕНТАЛЬНО-РАСЧЕТНЫЕ МЕТОДЫ ОПРЕДЕЛЕНИЯ ОТРАЖЕННЫХ СИГНАЛОВ

Импульсные характеристики объектов сложной конфигурации и сигналы, отраженные от таких объектов, могут быть получены экспериментально-расчетными методами.

Один из таких методов заключается в том, что определяется переходная характеристика объекта hb (£) посредством измерения величины потока, отраженного от модели реального объекта при облучении ее модулированным параллельным пучком света. При этом облучаемая площадь модели постепенно увеличивается посредством перемещения в направлении излучения плотно охва­ тывающего модель объекта раздвижного экрана (с последующим учетом вносимого им вклада в отраженный сигнал). Импульсная характеристика g6 (/) вычисляется посредством дифференциро­ вания полученной переходной характеристики h3 (t).

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

 

 

т

п

 

 

 

 

 

 

g, (/) =

s Sг’ч As?/ red (t -

х„),

(2.63)

 

 

i=1/=1

 

 

 

 

 

где r*j — гц cos у ц — эффективный коэффициент

яркости

эле­

мента поверхности

в направлении локационного наблюдения;

As*/ =

ASi/cos ytf — площадь

проекции

элемента

поверхности

Asu на картинную плоскость;

— 21ц/с — запаздывание сигнала,

отраженного от элемента

поверхности;

— расстояние от

эле­

мента

поверхности до

опорной плоскости,

перпендикулярной

37

к направлению наблюдения и проходящей через некоторую выб­ ранную точку, за которую удобно принять ближайшую к ОЭС точку объекта;

1

при

О < t < At

rect (/) =

при

i <

— стандартный импульс.

О

0, t> A t

При достаточно малом шаге разбиения картинной проекции объекта на элементарные зоны, пренебрегая краевыми нерегу­ лярностями, можно принять As*/ = As* = const.

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

величин зональных яркостей или г//) и суммированием

зональных

сигналов с учетом присущих им значений T tJ. При этом

величины

зональных расстояний могут быть определены

расчетным

путем либо с помощью тех или иных измерительных приспособ­ лений.

Зональный метод может быть распространен и- на определение

системных

импульсных характеристик (/)• При этом должны

быть учтены значения диаграмм направленности излучения

и приема

в центрах элементарных зон при

выбранном спо­

собе ориентации оси

ОЭС на объект. Тогда получим

 

 

m

п

 

(2.64)

 

gl (/) = Е Ег'чAs.VFi ,7Т2 ,7 rect (t -

%„).

 

1=1 /=1

 

 

Таким образом, в процедуру нахождения импульсных харак­

теристик входит определение матриц зональных яркостей

г,*. т ,

времён запаздывания

тп>т и диаграммы направленности

ЧГ„,т

с последующим расчетом gl (?) по соотношению (2.64), что может быть осуществлено на ЭВМ. Характеристика gH(£) является частным случаем характеристики £э (t) при равномерной диаграмме

направленности (Ч ^ =

= 1).

В некоторых случаях экспериментально-расчетным методом удобнее получать непосредственно характеристики отраженных сигналов при заданном виде зондирующих импульсов. При этом для определения характеристик элементарных сигналов (ампли­ туд и запаздываний) помимо рассмотренного выше метода может быть использовано сканирование узким непрерывным лучом по поверхности (модели) объекта (по принципу строчно-кадровой развертки) с измерением величин зональных сигналов и расстоя­ ний до элементов поверхности [33].

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

28

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

определен

в виде

 

а

b

 

 

 

 

 

Xi+T

yi+T

 

 

 

 

Is if (0 =

It if fit —т//)

j

J

lF (x,

y) dxdy,

(2.65)

 

 

 

a

h

 

 

 

 

 

xi~~ 2

‘2

 

 

 

где / tif— величина зонального

сигнала,

полученного

экспери­

ментально;

(х,

у) — значение

диаграммы

направленности

в точке х, у картинной плоскости объекта;

а и b — размеры сто­

рон элементов поверхности.

от объекта, выражается

в виде

Общий

сигнал,

отраженный

 

 

т

п

 

 

 

(2.66)

 

 

/э(0=2 Е /.«»■

 

 

 

 

/=1

 

 

 

 

Расчет отраженного сигнала по формуле (2.66) с учетом выра­ жения (2.65) может быть проведен на ЭВМ методом дискретной свертки.

Рис. 2.6. Нормированные временные функции сигналов, отраженных от крупногабаритного объекта: а — гаус­ сова диаграмма направленности; 6 — прямоугольная диаграмма направленности [точки—импульс вида (2.69), сплошные линии — вида (2.70)]

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

Рассмотрим гауссову диаграмму направленности

Y (X, у) = -^ п ехр [ - 4 (* - * " )2^ (у-

уо)а] ,

(2.67)

где D — диаметр пятна излучения в плоскости

объекта; x0f у0

координаты центра пятна излучения; kx — коэффициент пропор­ циональности.

39

Интегрирование по формуле (2.65) проводится с учетом зави­ симости (2.67), согласно выражению

* (+J т » (+4- У (х, jу) dxdy = -L [erf (

+ ° ■2х°- ) -

аи

* 1 - т ч — Г

_

erf ( —

д ~ 2ж° ) ]

[erf( 2yi + p ~ 2- ^ - ) -

erf ( 2g(^ * ~ 2y°) ] ,

 

 

 

z

 

 

 

 

(2.68)

 

 

 

 

 

 

 

 

где

erf (z)=

-~L=-J er^dt— функция

Крампа.

В

качестве

вре­

менных функций зондирующего импульса возьмем

 

 

/1/ти

при

О с / < т „

 

 

импульс

(2.69)

iv ) — \ л

___

^

v _ Q — прямоугольный

 

О

при

/ > т и, t

 

 

 

 

и

 

 

 

4.9 f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(t) =

fe T“

= а2/е-а *.

 

(2.70)

На рис. 2.6, а и б приведены нормированные временные функ­ ции сигналов, отраженных от крупногабаритного объекта при облучении его оптическими импульсами вида (2.69) и (2.70) дли­ тельностью 25-10"9 с.

§2.5. ОСНОВНЫЕ ВИДЫ ПОМЕХ В ИМПУЛЬСНЫХ' ОЭС

Косновным видам помех в импульсных ОЭС активного типа могут быть отнесены следующие: 1) флуктуацйонные шумы прием­ но-усилительного тракта; 2) естественные (фоновые) лучистые помехи, а также помехи, обусловленные отражением солнечной

радиации от рассматриваемых объектов; 3) помехи, связанные

спрохождением рабочего излучения через среду распространения,

ккоторым относятся диффузные помехи, обусловленные отра­ жением излучения от рассеивающей среды (обратное рассеяние),

и ослабление сигнала вследствие потерь лучистой энергии в среде, а также флуктуации сигнала, вызываемые турбулент­ ностью среды; 4) помехи, возникающие из-за отражения рабочего излучения от подстилающей поверхности; 5) внутренние лучистые помехи, вызванные частичным проникновением излучаемого сиг­ нала непосредственно в приемный тракт вследствие рассеяния излу­ чения на элементах конструкций (в том числе оптических); 6) внут­ ренние электромагнитные наводки, создаваемые излучающим устройством и воздействующие на приемный тракт; 7) микро­ фонный эффект (шумы вибраций); 8) различного рода искусствен­ ные (организованные) помехи.

40