Файл: Амитей Н. Теория и анализ фазированных антенных решеток.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 248

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Методы решения

121

при достаточно большом Р, считая, что при этом получается приблшкенное решение с необходимой точностью. Следуя описанной выше процедуре, выразим хР через x N:

Х р — У] ( — I)9 (bpiv)9 x N i

(68)

9 = О

где L pn = Kj} (К Р — KjY). Таким образом, используя выражение (68), можно уточнить первоначальное приближенное решение.

Заметим, что при получении выражений (67) и (68) необходимо увеличить размерность (т. е. следует расширить ранг K N), так как приближенное ядро K N обычно имеет меньшую размерность, чем ядра К и К Р. Это можно сделать разными способами. Удобно, например, левую верхнюю подматрицу N X N расширенной мат­ рицы приравнять первоначальной матрице K jY, а затем дополнить ее так, чтобы диагональные элементы были равны диагональным элементам К или К Р, а недиагопальные элементы были равны нулю, т. е.

(

для i,

j ^ N ,

 

 

N или j > iV,

KnE;j = КриЬи или К и 8ij для i >

 

 

 

(69)

где K ne .- — элементы расширенной матрицы.

Аналогично

x N

можно расширить до x NE, полагая x NE = {{%}, 0, 0, . . .,

0}.

В результате обращение расширенной матрицы будет не сложнее, чем обращение первоначальной матрицы.

Определим

теперь

последовательность векторов (х)\

i = 1 , 2 , . . . ,

как частичную сумму ряда Неймана (66), т. е.

 

(a;)° =

a;jY,

 

(а;)1 = хЕ LhXn ,

 

 

(70)

(я)8 = 2 ( - 1 rb%xN. 5=0

Каждый вектор можно представить в форме

(a;)r = {(a:)i, (xf2, (х)\, . . .}.

Интересно отметить, что первая коррекция равна нулю для первых N компонент (а;)1 — (а;)0:

(а:)1 — (х)° = {0, 0, ... , 0, Да:Л-+1, AxN+2, •••,},

(71)

где Aa:jY+1 = (а?)дг+1—(х)%+1 и т. д. Поэтому необходимо суммиро­ вать ряд Неймана, начиная со второго члена, чтобы получить уточнение для первых N компонент приближенного решения.


122

Глава 3

Выражение (65) представляет значительный интерес, так как с его помощью можно оценить отклонение приближенного решения от точного только через оператор L N и приближенное решение x N.

Возможны другие способы преобразования уравнения (1) в форму, удобную для применения итерационного процесса. Например, можно представить матрицу К в виде

К = N — Р.

(72)

Тогда уравнение (1) можно записать в виде

Nx = Рх + у.

(73)

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

xL= N^Pxo-^N-hj,

x2 = N -1Px1-\-N-1ij= N -1PN-1Px0-\-N-1P (N^y+y) . (74)

Очевидно, матрицу N надо выбрать так, чтобы легко было найти обратную к ней матрицу. Достаточное условие сходимости процес­ са состоит в том, чтобы наибольшее собственное значение матрицы N~XP было меньше 1. Это условие оказывается также необходимым, если процесс должен сходиться независимо от первоначального приближения £0.

Известно много способов расщепления исходной матрицы. Спе­ циальный случай N;j = К и 8и , P tj = Кц8и K tj, называется итерационным процессом Якоби, или процессом одновременных итераций. Процесс Гаусса — Зайделя, пли процесс последователь­ ных итераций, приводит к выбору N в виде

М

f

K-ijt

 

о,' i < j .

Особый способ построения расширенной матрицы K N, указанный

в соотношении (69), иногда

называется процессом блочных

итераций.

 

 

П Р И Л О Ж Е Н И Е 1

ЗАКОН СОХРАНЕНИЯ ЭНЕРГИИ И МЕТОД МОМЕНТОВ

Чтобы доказать, что решение, полученное методом Ритца — Галеркина, автоматически удовлетворяет закону сохранения энер­ гии, рассмотрим уравнение (22) для электрического поля в апер­


Методы решения

123

туре:

СО

2 ^ ( I V ( r ) =

j

{ 2

y j O j { T ) O j ( r ' ) +

 

 

A

i=1

 

 

 

 

 

+ 2

(!•)

(r') } • E( (r') dr' ■ (22)

 

 

 

7 7 1 = 1

 

Напомним,

что

уравиеипе

выведено для

случая возбуждения

решетки типом волны Ф,- (г).

Для перехода от уравнения (22) к системе лииейиьтх алгебраи­ ческих уравнений применим метод моментов с системой базисных

функций rjj (г), (i =

1, 2,

. .

I) и с системой весовых функций

т|* (г) (i

= 1, 2, . .

I).

В результате получим

 

 

 

di

I — [Lji]

di 1

,

(П.1)

где

 

 

 

 

/

 

 

 

dj = 2y}>{Ф}>,

 

 

(П.2а)

 

 

л*).

 

 

 

 

 

м

 

 

(П.26)

Ljt = 2

Уп ('Hi) ф») (Фп,

ф ) +

2 Y„ (Л*.

'Em) ('Em, Г]*)

7 1 = 1

 

 

 

7 7 1 = 1

 

 

 

и аг — неизвестные коэффициенты разложения. Заметим, что при получении системы (П.1) введено приближенное выражение для ядра путем обрыва двух рядов прп значениях N и М. Это экви­ валентно использованию N типов волн в волноводе и М гармоник Флоке в свободном пространстве при применении метода сшивания гармоник. Таким образом, электрическое поле в апертуре пред­ ставляется в виде

N

 

I

Eg ~ (1+ Д ,0 фу (г)+ 2

и,)Ъ Ж (Г) ^

2 ащ(г) (П.За)

Еа ~ 2 ^т'Ещ (г) ~ 2 аЭ1г (г).

(П.Зб)

771=1

1 = 1

 

Индексы g жа соответственно обозначают поле в волноводе и в сво­ бодном пространстве. Магнитные поля Hg и На определяются выражениями

 

N

 

— z х Нg=ijj' (1 — Ry) Фх (г) — 2 УпЪ'пФп (г),

(П.4а)

М

71=1

 

 

(П.46)

— Z X На= 2 УгНт'Е(г).

 

т=1

 

 


124

Глава 3

Из выражений (И.За) и (П.Зб) можно найти неизвестные модаль­ ные функции с помощью коэффициентов аь:

I

1 + Д ,'=

2

ОН.

Фу')в/.

(П.5а)

 

i=i

 

 

 

 

I

СПг,

Фп) at,

(П.56)

Ь п =

2

 

1=1

 

 

 

 

/

 

 

 

/1т = 2 ( ф , ЧГ^а,.

(П.5в)

 

г=1

 

 

 

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

 

 

 

Jz.E*xHJdr' = jz.EaxH;dr.

 

 

(П.6)

 

 

 

А

 

 

 

 

А

 

 

 

 

 

 

Используя выражения

(П.З) — (П.5),

найдем

 

 

 

2

2

Уго-i (Фг>

110 =

 

 

 

 

 

 

 

 

 

 

 

i= 1

 

I

и

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

2

2

aiaj [2

у*(1li,

Фп) (Фп,

11*)+

 

 

 

 

 

i = l j = l

 

п = 1

м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

2

Y % ( r u ,

4 + )

('F m, + ) ] .

(П.7)

 

 

 

 

 

 

7 7 1 = 1

 

 

 

 

 

 

Из сравнения выражений (П .1)

и (П .7)

вытекает,

что выражение

(П .7)

является квадратичной формой, полученной из выражения

(П .1):

 

 

/ .

\ *

 

 

 

 

/

\

*

 

 

 

 

 

 

 

 

 

 

 

 

{ • • • a-i

■■■} I

di

= { . . . at

...} [L}i\* I

at

I .

(П.8)

Таким образом,

мы показали,

 

что

выбор Ln = (t]|, Lip)

ведет

к сохранению энергии. Соотношение взаимности при таком выборе не удовлетворяется при комплексных значениях {тр}. Однако в ряде задач используются действительные значения {тр} и тогда приближенное решение по методу Ритца — Галеркина удовлетво­ ряет как принципу взаимности, так и закону сохранения энергии.


Методы решения

125

П Р И Л О Ж Е Н И Е 2

ПОСТРОЕНИЕ АНАЛИТИЧЕСКОЙ ФУНКЦИИ В МЕТОДЕ ВЫЧЕТОВ

Функция g (и?) выбирается в следующей форме:

СО

л

 

р И ( и»-Гй [] (ш -Г ')(ш -П д)(й/25я)2Р0^

 

_________ Я=±_________________________

(П.9)

оо

 

{^ + Уо)(ч’ — У'о) П (Tn — w) W nn) ewd/nx

п—1

где р (w) — целая функция, определяемая ниже. При выборе g (w) распределение нулей и полюсов принимается в соответствии с усло­ виями 1 и 2, указанными в разд. 6. Экспоненциальные множители в бесконечных произведениях введены для обеспечения сходимо­ сти. Присутствие р (w) необходимо для выполнения условия 3, т. е. | g (if) | ->- 0 при | if | — оо. Для определения функциональ­ ного вида р (if) полезно заметить, что для больших т

т п

И

2тл

Ут ■

d

 

 

Следовательно, g (w) в окрестности бесконечно удаленной точки отличается на ограниченную функцию от функции

Рм ( » - г » ц ( » - ^ )

( 4 г ) 2*” "к

glW = -------------------------------- - ---------------------------------------------

=

<»+«<”-»» п ( т - “) Ш

7 1 = 1

М » ) ( » - Г 1 ) Г ( ^ ) Г ( ^ ) , | Д ( ^ * ) , , „ ( ^ )

я (“' + ?£) (u> —Yo) г ( ^ - ) sinwd

(П.10)

где Г (if) — гамма-функция. Итак, для определения поведения функции g (if) при больших if необходимо только исследовать изменение функции gx (w) при w -> оо. Если w стремится к беско­ нечности, первый член формулы Стирлинга для гамма-функции имеет вид

Г (z) ~ У 2п/z е2 |п2-2.

Подстановка асимптотического выражения для гамма-функции в выражение (П.10) приводит к следующей асимптотической фор­ муле:

gi М ~ !п 2.