Файл: Амитей Н. Теория и анализ фазированных антенных решеток.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.