Файл: Амитей Н. Теория и анализ фазированных антенных решеток.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 10.04.2024
Просмотров: 245
Скачиваний: 2
Методы решения |
101 |
коэффициентов связи между гармониками. Напомним, что в урав нениях для электрического поля в элементы матрицы входили модальные проводимости. Важно отметить, что модальные сопро тивления внешнего пространства присутствуют только в диаго нальных элементах-матрицы, а модальные сопротивления внутрен ней области — во всех элементах.
Обращая матрицу системы (37), найдем амплитуды простран ственных гармоник (модальные токн) тангенциального магнитного поля. Амплитуды типов волн в волноводе можно определить из условия непрерывности [выраженпе, аналогичное формуле (28)]. В результате получпм
м
(38)
Р = 1
После вычисления амплитуд гармоник можно найтп коэффициенты отражения и передачи.
Элементы матрицы выражений (24а) и (37а) содержат бесконеч ные суммы, получающиеся из произведений модальных проводи мостей (или сопротивлений) н коэффициентов связи между гармо никами. Можно показать, что эти суммы являются сходящимися рядами по крайней мере для решеток из прямоугольных и круглых волноводов, для которых известны выражения для 4§]q в явном виде. Однако суммирование рядов с точностью до пяти и более значащих цифр требует учета многих членов. Как мы увидим ниже, это эквивалентно учету большего числа гармоник в одной из областей. Кроме того, для большинства ФАР размер апертуры элемента сравним с размером единичной ячейки. Это позволяет ожидать, что для описания поля в плоскости решетки требуется примерно одно и то же число гармоник независимо от того, из какой области они берутся. На этом основании суммы в элементах матрицы можно оборвать на номере, сравнимом с порядком матри цы. Например, в выражении (24а) суммирование можно прекра тить при т — М, где М ~ N. Таким же образом в выражении (37а) сумма вычисляется до / = N, где N ~ М. Очевидно, при этом в решение вводится новый источник приближения. Однако это приближение незначительно ухудшает точность решения.
Если раскрыв волновода и единичная ячейка имеют существен но различные площади, число используемых гармоник в каждой из областей сильно зависит от геометрических факторов. Интерес но отметить, что при использовании конечного числа гармоник, например N типов волн в волноводе и М пространственных гармо ник во внешней области, решения, полученные из уравнения для электрического поля и полученные из уравнения для магнитного поля, совпадают. В частном случае решения системы (24) с матрич-
102 |
Глава 3 |
iimra элементами (24а), суммируемыми до М, и решения систе мы (37) с матричными элементами (37а), суммируемыми до N, одинаковы.
4.2. Метод сшивания гармоник
Системы уравнений (24) и (37) выведены в два этапа. Сначала граничная задача строго формулировалась в виде интегральных уравнений для тангенциальных составляющих электрического плп магнитного поля. Затем с помощью метода моментов интег ральные уравнения сводились к системам линейных алгебраи ческих уравнений. Такие же системы уравнений можно получить 'с помощью приема, называемого методом сшивания гармоник [10, 11]. По этому методу находят разложения поля на гармоники в каждой области, а из условий непрерывности на границе между областями получают два отдельных соотношения между коэффи циентами разложения: одно соотношение для непрерывности тан генциальных составляющих электрического поля, а другое — для непрерывности тангенциальных составляющих магнитного поля. Этн соотношения можно записать в виде бесконечных систем
уравнений. Обычно считается, что |
в апертуре А возбуждается |
||
конечное число типов волн. В результате получается |
конечная |
||
система уравнений. |
|
|
|
Рассмотрим бесконечную решетку, для которой Ч*",,, (г) являются |
|||
векторными модальными |
функциями внутренней |
области, |
|
а Ф7- (г) — гармониками |
внешнего |
пространства (гармониками |
Флоке). Еслп решетка возбуждается типом волны Ф,- (г), танген
циальные |
составляющие |
электрического п |
магнитного |
полей |
||
в апертуре |
определяются |
выражениями |
|
|
||
|
|
|
лг |
|
м |
|
E t( г) ~ (1 + Ry) Фу (г) + I |
^ щ Ф } (г) ~ |
У Fm'Fm (г); |
|
|||
|
|
|
1— 1 |
7П =1 |
|
|
|
— z X Н/ (г) ~ уу (1 — Rj.) Фу (г) — |
|
||||
|
N |
Р ^ . ф . ( г) ~ |
М |
|
(39) |
|
|
- У |
у YmF;/Fm(r). |
||||
|
j = |
i |
|
771=1 |
|
|
Знак приближенного равенства используется в этих выражениях, потому что учитывается конечное число гармоник. Уравнение для магнитного поля справедливо только в пределах апертуры волно вода, в то время как уравнение для электрического поля верно на всей единичной ячейке, причем Е* = 0 в области С — А. Умножая первое уравнение на Фр (г), интегрируя в пределах С п учитывая ортонормироваиность гармоник, найдем
П ~ (1 + Д ,- 0 « /.р+ 2 Щрг'у Р = 1 , 2, .. . , М, |
(40а) |
i = i
Методы решения |
103 |
где ’Sjp — коэффициенты связи между гармониками, определяемые по формуле (246). Точно так же, умножая второе уравнение (39) на Фп (г) и интегрируя в пределах А, получим
м
|
— R j ' ) S j ' n — (1 — 6j-n) l/nV'n~ |
2 |
|
YmV'm^nm, |
|
|||||||||
|
|
|
|
|
|
|
|
m=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n = |
l, |
2, ..., |
W. |
(406) |
||
Рассматривая v'n (n — 1, |
. . ., |
N) при |
= |
R y |
|
и 7 и ( и |
= |
1, . . . |
||||||
. . ., M) как |
неизвестные величины, |
запишем |
|
выражения (40а) |
||||||||||
и (406) в матричной форме: |
|
|
|
|
|
|
|
|
|
(41) |
||||
где |
|
|
|
[А] а = |
х, |
|
|
|
|
|
|
|||
|
|
N |
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
ср* |
• • ■ ^N 1 |
— 1 |
0 |
|
. . |
|
0 |
|
|
|||
“ « f l |
® 2 1 |
|
|
|
|
|||||||||
м- |
|
ср* |
• • ■ ^ N 2 |
|
|
|
|
|
|
|
|
|
|
|
Щ г ® 2 2 |
0 |
- |
1 |
|
|
|
|
0 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
ср* |
|
Ср* |
б |
|
|
|
- |
1 |
|
|
|
|
[А] |
Щ м &2М • • |
|
|
|
|
|
|
(41а) |
||||||
У1 |
|
|
. 0 |
Г & и Г 2 в 12 • • • У а г « Ш |
|
|||||||||
|
о |
. |
|
|
||||||||||
N- |
0 |
Уг |
|
|
Y ^ 21 У 2 « 2 2 |
|
|
|
|
|
|
|||
_ б |
|
|
Уя Y & N 1 |
— |
|
(pin |
Y m %n m |
|
|
|||||
|
|
|
v'i |
|
|
|
-- |
|
— |
|
|
|||
|
|
|
|
|
|
|
W j ' l |
|
|
|
|
|||
|
|
|
V2 |
|
|
|
|
|
Cp* |
|
|
|
|
|
|
|
|
|
|
|
|
— V j ' 2 |
|
|
|
||||
|
|
|
V j ’- l |
|
|
|
— Щм |
|
|
|
||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
R; |
|
II |
X = |
|
0 |
|
|
|
|
(416) |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
v’n |
|
|
|
У1 |
|
|
|
|
|
||
|
|
|
VI |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V'2 |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
_ |
Vm _ |
|
|
|
|
6 |
|
|
|
|
|
|
Это система из М + |
N уравнений с М + |
N неизвестными. Инте |
ресным свойством матрицы является ее разделение на четыре под матрицы, из которых две оказываются диагональными. Это свойство наводит на мысль, что выражения (40а) и (406) можно записать в следующей форме:
V' = [i?]v' + b,
d - [ y ] v ' = [S][Y]r.
104 Глава 3
Получим связанную систему уравнений, в которой неизвестными являются
|
|
|
|
|
|
v l’ |
|
|
|
у; - |
|
|
|
V2 |
|
У ' = |
|
У 2 |
и |
v' = |
|
v’j'-l |
(43a) |
|
|
У м . |
|
|
|
Rу |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
_ |
v 'n |
_ |
Матрицы [Л] и [iS] имеют |
порядки |
М х N и N X М соответ |
|||||
ственно и задаются формулами |
|
|
|
|
|||
|
|
|
*© 21 |
• • • |
® iV 1 |
|
|
[ i ? ] = |
^ 1 2 |
^ 2 2 |
• |
|
|
(436) |
|
|
|
|
|
|
|||
и |
|
_ « * M « 2 М • • • « jY M |
|
||||
|
|
|
|
|
|
|
|
|
|
“ « 1 1 |
« 1 2 |
• • • « 1 M |
|
||
[ S ] |
= |
|
« 2 2 |
■ • |
. |
|
(43в) |
|
|
|
|
|
|
||
|
|
_ « * ! « 1 V 2 ■ ■ • « ivivt „ |
|
||||
Матрицы проводимостей являются диагональными |
|||||||
|
|
~ y 1 |
0 |
. . . |
0 |
" |
|
|
|
0 |
г/2 |
|
|
|
|
[ * /] = |
|
|
|
|
|
||
|
|
. 0 |
|
U n |
_ |
|
|
|
|
|
0 |
. . . |
0 |
■ |
|
m |
|
0 |
r 2 |
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
_ o |
0 |
. • • |
У м . |
|
|
Векторы Ъ и d имеют |
вид |
|
|
|
|
|
|
|
|
|
|
|
“ 0 “ |
|
|
|
« Ь |
|
|
|
0 |
|
|
b = |
« Ь |
|
|
|
|
||
и |
(1 = |
|
УУ |
|
|||
|
«Ум |
|
|
|
0 |
|
|
|
|
|
|
|
|
0
Методы решения |
105 |
Подставляя V' из выражения (42а) в выражение (426), исключим У', чтобы получить уравнение только для и'. После небольших преобразований находим
т + [5] [У] [/?]) V' = d - [Y] [5] b. |
(44) |
Ранг матрицы ([г/] + [£] [У] Ш]) равен N. Число неизвестных компонент вектора У' также равно N. Запишем q-e уравнение в форме
N N М
2 |
{г/Ал + (2 |
Ym^qm'Spm)} Vp =ijj'8qj' —2 Y mcSqm%fm. (44а) |
p = i |
7?i=i |
m = i |
Полагая v'j> = 1 + Rj'i уравнение (44a) можно записать в виде
N |
М |
|
2 |
{ у А р + ( 2 Y nlVqmWpm) ) v'v = 2yj.6qr |
(446) |
р = 1 |
m = l |
|
Это уравнение идентично уравнению (24), выведенному для электрического поля с помощью метода моментов при использова нии в качестве базиса волноводных мод. После того как найдены о', из выражения (42а) можно определить V'. Преобразуя выражение (42а), можно показать, что оно совпадает с выражением (29).
Исключая г/ из выражений (42а) и (426), можно записать урав нение только для V':
№ + [Л] [z] [5]) Г = |
[R ] [2 ] d + Ь, |
(45а) |
|||
где |
"1/У, |
0 |
. |
0 |
|
|
|
||||
[2] = [ГГ 1 |
0 |
1/У2 |
|
|
|
|
|
|
|
|
|
|
.0 |
0 |
. |
1/У |
|
|
[2/Г1 п |
l' = [Y] V'. |
|
Выражение (45а) является системой из М уравнений для М неиз вестных Г . Уравнение с номером q имеет вид
МN
2 { Z A * + |
2 г „ ^ л?) / ; = 2 ^ , |
(456) |
р = 1 |
771=1 |
|
Сравнивая уравнения (456) и (37), можно заметить, что эти два уравнения отличаются только множителем zp в правой части. Это вызвано тем, что при выводе уравнения (37) падающая волна имела единичный модальный ток для волны Фр, а при выводе уравнения (456) падающая волна имела единичное модальное напряжение для той же волны. Отношение амплитуд падающих волн равно модальному импедансу zp. Так как система линейна, это отношение оказывается просто множителем пропорциональ ности между двумя решениями.