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

книги из ГПНТБ / Зелкин Е.Г. Линзовые антенны

.pdf
Скачиваний:
97
Добавлен:
24.10.2023
Размер:
9.22 Mб
Скачать

но противоположную точку, а линза в виде полуцилин­

дра или полусферы

преобразует

расходящийся пучок

в

параллельный.

 

 

 

 

 

Чтобы определить траектории лучей в линзе, вернем­

ся

к уравнению

(4.9),

подставив

в него только

что най­

денное

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

n(R):

 

 

 

 

в - « = f

 

*П + * ) < « ^

( 4 . 2 4 )

 

 

 

'

 

 

 

Введем

обозначение:

 

 

 

 

 

 

 

m =

ft/2«i

(4.25)

и учтем, что /?л=Т,

тогда (4.24)

запишется в виде:

 

 

«

f

 

 

(4.26)

Этот интеграл нетрудно свести к табличному,

обозначив

через

 

 

 

 

 

 

 

 

z, = (Я8 1) m/Л ] Л - 4 т 2 .

(4.27)

В новых переменных уравнение (4.26) можно записать так:

Z

 

 

 

 

Ъ - * = [

*

Ь

= Т .

(4.28)

о

У

1

~ г '

 

Интеграл справа определяет

arcsinzi, а

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

sin (6—л) =Zi

 

или

 

 

 

 

3 — 1)/п//?]Л

— 4m* = sin9.

(4.29)

Удобно ввести декартову систему координат: x=i/?cos8J

/y=J?sin0. При этом

(4.29)

преобразуется к

виду:

 

 

х- -+- у2

+ 2у ] Л -

4m2/2m - 1 = 0 .

 

(4.30)

Это

уравнение

окружностей, центр которых

смещен

вдоль оси Y на величину

 

 

 

 

 

уа = -У

1 - 4 m 2 / 2 m = - c t g e , .

 

(4.31)

Радиусы

окружностей

равны:

 

 

.

; i

? o = n

i /

/ x = cscei.

..

(4.32)

110

Здесь '0i — по-прежнему угол выхода

луча

из

фокуса

линзы.

Траектории

лучей, определяемые

уравнением

(4.30), могут быть записаны в другой

форме через пара­

метр 01.-

 

 

 

 

 

 

х22+2у ctg0i—1=0.

 

 

(4.33)

Задавая

различные значения или 0ь можно найти тра­

ектории всех лучей в линзе.

 

 

 

Формулы (4.22),

(4.30) — (4.33) полностью

решают

поставленную задачу. Действительно, уравнение (4.30) является уравнением окружностей, симметричных отно­ сительно оси Y, а значит, все лучи, вышедшие из точки А, пересекают ось Y ортогонально. Если отбросить пра­ вую часть линзы и иметь правее Y среду с постоянным коэффициентом преломления, то лучи не изменят своего направления и будут идти параллельным пучком. Край­ ние лучи, выходящие из облучателя под углами ±л/2, согласно уравнениям (4.31) и (4.32), имеют вид полу­ окружностей единичного радиуса с центром в начале координат.

Итак, линза Максвелла в виде половины кругового цилиндра преобразует гомоцентрический пучок лучей в параллельный; облучатель такой линзы должен быть линейным.

Линза Максвелла со сферической симметрией может

быть получена заменой

R на Rlt где Rl =

~]/rx2-\-y2-\-z" ,

т. е.

 

 

я (/?,) =

я (0)/(1+/?;).

(4.34)

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

Коэффициент преломления в такой линзе, в частно­

сти, можно определить по формуле

 

 

(4.35)

а траектории лучей

 

х2+y2+2yR ctg 0i—R2 = 0.

(4.36)

111

4.1.3. АМПЛИТУДНОЕ РАСПРЕДЕЛЕНИЕ В РАСКРЫВЕ

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

p = tg

(8,/2).

(4.37)

Рассмотрим два луча, выходящие

из облучателя

под углами 0i и

Si + dQj. Они выйдут из раскрыва в точках р и p + dp соответственно. Обозначим через P{&i) диаграмму направленности облучателя лин­ зы, а через Q(p) —распределение энергии в ее раскрыве. На основа­ нии закона сохранения энергии имеем следующее соотношение:

/>(e,)</e,=Q(p)dP.

(4.38)

Из уравнения (4.37) найдем:

 

det /dp=cos2 (9i/2),

 

следовательно:

(4.39)

Q(p)=P(G,) cos2 (0,/2).

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

0.75

1

1

1

1

1

0.5

1

1

1

1

1

 

 

 

 

 

О

0.2

OA 06

0,8 j>

Рис. 4.2. Перераспределение мощности, вызываемое линзой Максвелла.

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

T=4n(p)/[l+;i(p)]2 ,

таким образом, реальное распределение энергии в раскрыве будет

 

п

, ,

р , в ,

4/; (р) cos2(91/2)

(4.40)

 

Q '

(Р) =

я ( 9 - )

[ l - M ( p ) ] 2

Обозначим второй сомножитель в формуле (4.40) через В(р),

тогда

Qi(p)=P(8,)5(р).

 

 

 

На

рис. 4.2 приведены значения В(р), рассчитанные для 8] от 0

до ±я/2 . Умножив В(р)

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

облучателя

линзы, получим распределение мощности в раскрыве.

 

112

Распределение амплитуды равно корню квадратному из распре­ деления мощности (4.40). Диаграмма направленности линзы Макс­ велла может быть рассчитана теперь обычным способом.

4.2.Линза Люнеберга

В1944 г. Люнеберг предложил линзу из неоднородно­

го диэлектрика с радиальной симметрией, получившую

вдальнейшем название линзы Люнеберга.

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

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

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

определена

им в предположении,

что ti(R) — непрерыв­

ная монотонно убывающая функция в интервале

0 ^ # г ^ .

< ; 1 и что коэффициент преломления на краю

линзы

(при R = l)

равен

единице. Последнее предположение

является важным

для хорошего

согласования

линзы

с окружающим пространством: в этом случае почти от­ сутствует отражение лучей на краю линзы.

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

8—342

ИЗ

линзы, т. е. #1=1, а другой в бесконечности (#2 = ° ° ) - В этом случае линза преобразует сферически расходя­ щуюся волну в плоскую.

4.2.1. КОЭФФИЦИЕНТ ПРЕЛОМЛЕНИЯ

 

Рассмотрим,

как и раньше, цилиндрическую

линзу

и ход лучей в

какой-нибудь плоскости, например

z=0.

Уравнение эйконала и его решение определяются попрежнему выражениями (4.1) и (4.6). Из (4.6) следует, что уравнение луча, проходящего через заданную точку Ro, 0о, будет иметь вид:

 

R

 

 

 

8 = e0 =tA j f l W / f l l A r ^ - f t 3 .

 

(4.41)

 

Ro

 

 

 

Определим пределы

изменения

постоянной

Клеро

/г. Для

центрального луча

01 = 0, /г = 0.

Крайнему

лучу,

выходя­

щему из фокуса под углом 01 = л/2, соответствует n(R) =

= 1 и, следовательно, А = 1, т. е. связка лучей

определя­

ется параметром h, изменяющимся в пределах

O ^ / i ^ l .

Из (4.11) следует, что в точке С, для которой <а=л/2, луч претерпевает полное отражение, т. е. до этой точки луч распространяется из менее плотной среды в более плотную, а после нее, наоборот, из среды более плотной в менее плотную. Следовательно, в этой точке коэффи­ циент преломления на пути данного луча имеет макси­ мальное значение и, согласно (4.11), точка С находится на наикратчайшем расстоянии от центра линзы, т. е.

Rc=Rmin-

Точку С(9', Rmin) назовем

точкой поворота.

В этой точке меняется знак dQjdR,

а

также знак перед

интегралом в выражении (4.41).

 

 

Уравнение луча, проходящего через точку поворота,

согласно

(4.41),

можно

записать в

виде

 

 

R

 

 

 

 

6 =

6'=t J

dR/RVn-R'-h2.

 

(4.42)

Поскольку показатель преломления стоит под интегра­ лом в комбинации с R, то удобнее ввести новую пере­ менную р:

p(R)=n(R)R. (4.43)

В выражении (4.42) закреплен пока один предел. Второй предел установим из условия прохождения лучей через 114

фокусы с координатами

(Ri, я) и Л2 с координатами

(R2, я) . Для этого положим

в выражении (4.42) 0 = я и

О соответственно:

 

 

 

R,

 

* = 6' + h

Г

dRIRYf-lf-

R,mm

 

0 = 6' -h

J

dR/RVf-h2.

R

'.

 

 

mm

 

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

через точку

поворота. Исключая 0',

получаем

 

 

 

Rmi.n

 

 

« = А

- h [

™ .

(4.44)

Разобъем каждый интеграл на два так, чтобы отдельно

интегрировать

по участкам

пути внутри и вне линзы:

 

1

«*

 

 

 

 

Ri

1

 

 

 

+ Л

f

+ * f

/ *

(4-45)

Первый и четвертый интегралы являются табличными.

Действительно, вне линзы

в свободном пространстве

n=l, a p = R, следовательно,

первый и четвертый инте­

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

i

 

 

 

 

 

 

 

 

h f r

d R

- h 2

= — arcsin /г -4- arcsin

 

(4.46)

)

RVR2

 

 

 

Ri

к

'

Ri

 

 

 

 

 

 

 

 

Ri

 

 

 

 

 

 

 

 

h Г , d R

h2

=

acrsin-J -

+ arcsinh.

 

(4.47)

J

R Vr.2

 

#2

 

4

y

1

 

и

(4.47),

можно

переписать

 

(4.45)

Учитывая (4.46)

 

в виде:

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

" 1 - * 7 Т = Р - т

 

 

( 4 - 4 8 )

Rmin

115

 

где

f (h) = arccos h + -тр ^arcsin-^- + a r csin-^-^.

Задача сведена теперь к решению интегрального урав­ нения относительно переменной р. Можно привести это интегральное уравнение к уравнению типа Абеля, перей­

дя к интегрированию по р. Пусть R есть функция

р, т. е.

# = # ( р ) , a dR=R'(p)dp.

При R=l n(R) = l и р

также

равно единице. При R=Rmin,

согласно

(4.11),

р = 1г, так

как в точке поворота

угол,

образуемый

лучом

и

радиу­

сом R, равен я/2. Уравнение (4.48) можно записать те­ перь следующим образом:

h

Г

 

R'<?}^_=-h

f _ * E _ = / ( A ) ,

(4.49)

 

ft

 

 

 

Л

 

 

 

 

где через

Q(p)

обозначена величина

 

 

 

 

 

 

 

Q(p)= — lntf(p) .

 

(4.50)

Выражение

(4.49)

есть интегральное уравнение

типа

Абеля. Для его решения предварительно умножим

пра­

вую и левую

части

на dhj У 2 —р2

и проинтегрируем

от

р до 1:

 

 

 

 

 

 

 

 

 

 

'

'

 

hdhdQ

_. Г

f (h) dh

, 4

5

1 )

 

И V р2 — Л 2 / Л 2 - ?2

J / Л 2 — р2 *

 

 

 

 

р Л

 

 

р

 

 

 

 

Меняя порядок интегрирования в левой части и приме­

няя

формулу Дирихле,

получаем:

 

 

1 Л

 

р г (

 

dCl.

С Г

hdhdQ.

Г

hdh

J J

K F ^ F * К A2 — p2

J

J yt* -

№ VW

 

p 1

 

1

*-p

 

 

Внутренний интеграл равен я/2, и потому правая часть выражения (4.51) равна n[Q(p)—Й(1)]/2. Подставляя вычисленный интеграл в (4.51), имеем:

- : i o W - o ( i ) i - j j ^

116

Подставив f(h), можно записать

0 ( р ) - 0 ( 1 ) = 4 - ( : т = М +

" JP К Л 2 f-

1

 

 

arcsin (A//?,) + arcsin (h:R2)

\ ^

^

V 2 p2

1

 

p

Первый интеграл легко вычислить. Действительно, чи­ слитель подынтегрального выражения совпадает с пра­ вой частью выражения (4.46), если положить в нем Ri =

— h:

 

 

 

 

 

 

i

 

 

d f *

 

- = arccos/z.

 

(4.53)

 

 

 

 

 

 

J

 

 

 

 

 

 

 

 

 

 

R VR* -

h-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя

(4.53) в первый

интеграл и меняя порядок

интегрирования,

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 arccos hdh

= —1пр.

 

 

(4.54)

 

 

 

 

 

 

 

 

 

Vh? — р2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р

 

 

 

 

 

 

 

 

 

 

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

функциям.

Следуя

Люнебергу,

обозначим

через

q{p,

а)

интеграл*

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

,

 

1

(* arcsin (1/а)

 

 

 

 

гг.,

 

 

 

 

 

Ч(Р, a)=—\-y=±J-dt.

 

(4.55)

 

 

 

 

 

 

 

 

 

 

р

 

 

 

 

 

 

 

 

 

Учитывая

(4.50)

и

(4.54),

можно

 

 

(4.52)

переписать

в виде:

 

 

 

 

 

1п(р//?)= 9 (р,

Я,)+<7(р, Я а ) .

 

(4.56)

Но p = nR,

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

In n = q(p, Ri) + q(p,

 

R2),

 

 

 

отсюда

получаем

-

>

 

и /? =

ре

*

 

-

 

* «"*»> . (4.57)

n =

e

q ( p

-

R

'

) + 4

]

f p

л о

 

 

 

 

{f R

 

 

 

-

 

 

 

(

я

Уравнение (4.57) в параметрической форме определяет закон изменения коэффициента преломления как функ­ цию радиуса R.

Рассмотрим некоторые частные случаи.

* Таблица значений функции <?(р, а) для а от 1 до 2 и р от 0 до 1 приведены в приложении к книге.

117

1. Пусть R2=oo.

Из

(4.55) следует,

что q(p,

оо)=0,

следовательно:

 

 

 

 

 

 

 

 

 

n = tq{?'Rt\

R^pe-"^^

 

.

 

(4.58)

В этой линзе пучок лучей, выходящих из точки R[,

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

пучок.

 

2. Пусть

Ri=\,

# 2 = 00,

тогда

 

 

 

 

 

д = е 9 ( р '

 

^ = Р е - " р ' " .

 

 

Вычислим функцию

 

 

I

 

 

 

 

 

 

,

 

I f

arcs

 

 

 

 

 

9 ( р '

 

1

(*

arcsin/z .dh.

 

 

 

 

l ) = - ) v w . •p-

 

 

 

Согласно

(4.54)

 

 

p

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2?(p,

l) = I n P +

J F = ^ = = l n ( l + V T = p T ) ,

следовательно

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- i - m (i + VT=F)

 

+ "J/ 1 -

p3 ,

(4.59)

/г = е 2

 

 

= \f1

Отсюда

 

 

 

 

 

 

 

 

 

 

/г +

/?= =

2,

 

л =

У 2 -

R\

 

(4.60)

При расчете линз предполагалось, что их радиус ра­ вен единице. В случае, если линза имеет радиус, равный Ь, в формулах (4.48) — (4.60) вместо R необходимо по­ ставить величину Rib.

Итак, формула (4.59) выражает закон изменения ко­ эффициента преломления в обычной линзе Люнеберга. Она получила наибольшее распространение, так как ока­ залась наиболее удобной для практического„применения.

Линза Люнеберга выгодно отличается от линзы Макс­ велла тем, что допускает качание диаграммы направлен­ ности без искажения ее формы в широком секторе углов вплоть до 360°. Это обеспечивается тем, что линза имеет центральную симметрию и форму целой сферы или ци­ линдра. Плоский фазовый фронт образуется вдоль неко­ торой плоскости, находящейся вне тела линзы. Благода­ ря этому при перемещении облучателя линзы по сфере US

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

/

1^Y

 

/

 

 

/

 

\ 11 X

 

' \ ^

^

V

7

\ \

Р

Рис. 4.4. Качание диаграммы направленности в линзе Люнеберга.

Отметим, что линза Максвелла является также ча­ стным случаем рассмотренной выше линзы с центральной симметрией. Ей соответствуют значения Ri=\ и Rz=\. Согласно (4.57), коэффициент преломления в данном случае выражается следующим образом:

а так

как

2<?(р, 1) =

1п(1

— р2),

как было

показано

выше,

то

коэффициент преломления в

линзе с

R,= l и

1

равен:

 

 

 

 

Подставляя вместо

p=nR,

легко показать, что

 

 

 

 

/г = 2 / ( 1 + ^ .

 

(4.61)

Выражение (4.61) совпадает с выражением для коэффи­ циента преломления в линзе Максвелла (4.22) при пред­ положении, что коэффициент преломления на краю лин­ зы равен единице.

119

Соседние файлы в папке книги из ГПНТБ