Файл: Амитей Н. Теория и анализ фазированных антенных решеток.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 10.04.2024
Просмотров: 247
Скачиваний: 2
Методы решения |
117 |
сти (гармоники Флоке). При выборе одного вида функций в каче стве базисных, а другого — в качестве весовых элементы матрицы получаются более сложными, чем при использовании только одно го вида функций. Однако, если стенки волноводов бесконечно тон кие, выбор смешанного базиса может привести к значительному упрощению элементов матрицы, так как апертура волновода теперь занимает почти всю площадь единичной ячейки. Более того, функ циональная форма элементов матрицы такова, что можно получить точное аналитическое решение бесконечной системы линейных уравнений.
Для иллюстрации возьмем решетку, возбуждаемую низшим типом волн ТМ (можно также взять волну ТЕМ) и сканирующую в плоскости Е (случай возбуждения волной ТЕ анализируется аналогично). Интегральное уравнение для электрического поля имеет вид
d/2 °о
j{ 2 i # W ) (W ) +
-c l/2 5=0
oo
+ 2 |
(53) |
Неизвестное электрическое поле можно представить в виде следую щего разложения:
оо |
|
Еу(у) = 2 УпФпЫ- |
(54) |
п=0 |
|
Подставляя выражение (54) в уравнение (53) и вычисляя моменты с функциями {XF| (у), q = 0, + 1, ± 2, . . .}, после ряда алгебраи ческих преобразований получим (см. гл. 4).
00
sin (ф/2) |
2 |
~Тг^ Р ~ Вп для 9= 0>± !> ± 2> • • ■. |
(55) |
|||
V 2 т а - ж |
||||||
Yi (Г'+Vi) |
7 1 = 0 |
1 а |
гп |
|
|
|
|
|
|
|
|
|
|
где неизвестные |
{В п} |
связаны с {v'n) соотношением |
|
|||
|
|
|
sin |
Ф |
для четных значении п, |
|
|
|
|
-д- |
(55а) |
||
Еп— ( — 4)п/2 У (2 — &on)/d |
|
i |
|
|||
cos |
Ф |
для нечетных значении п , |
|
|||
|
|
|
|
|||
где и' = 1 -{- у" и v'n = v'n при п > |
0. Система уравнений |
(55) |
||||
имеет особую форму. |
Коэффициент |
1/(Гд — у'п) при Вп весьма |
прост и зависит только от разности (Г^ — уй). Уравнения такого типа решаются с помощью функционально-теоретического метода, известного под названием метод вычетов. В этом методе вводится интеграл от специальной аналитической функции и затем вычеты
118 Глава 3
подынтегральной функции связываются с неизвестными коэффи
циентами |
Вп. |
|
|
|
|
Рассмотрим контурный интеграл |
|
||
|
|
|
8 И |
|
|
|
CsIш —Г—д dw, |
(56) |
|
где g (w) — аналитическая |
функция, |
обладающая следующими |
||
свойствами: |
|
—у„ и у'т, т — 0, 1, 2, . . |
||
1) g (w) имеет простые полюсы при w = |
||||
2) |
g (w) имеет простые нули при w = |
T'q, д = 0, ± 1, ± 2, . . . ; |
||
3) |
I g (w) I |
0 при | w| |
оо; |
|
4)lim (w + у'0) g (w) = — Y2UI --1П„У/2-.
u’^-v0' |
'о |
В выражении (56) Cs обозначает простой замкнутый контур, окру жающий первые s полюсов. По теореме Коши о вычетах контурный интеграл равен произведению 2nj на сумму вычетов в полюсах. Когда Cs уходит в бесконечность, согласно условиям (3) н (4) кон турный интеграл равен
2 |
|
'•(Yn) |
| |
Y~2jd' sina|)/2 |
— = 0, |
(57) |
|
|
Уп~Га |
Уо |
Уо+ |
Г |
|
||
п —О |
|
|
|
|
|
|
|
где |
|
г (у’п) = |
lim (w —y'n)g(w). |
|
|
||
|
|
|
|
||||
|
|
|
|
w-*y' |
|
|
|
Из сравнения выражений (55) и (57) формально получаем |
(58) |
||||||
|
В п = г |
(у'п) = l i m ( w — |
У п ) g ( w ) . |
||||
|
|
|
|
W-+V |
|
|
|
Таким образом, определение неизвестных коэффициентов Вп свелось к нахождению соответствующей функции g (w). Детали построения такой функции обсуждаются в приложении II.
Пределы применения метода вычетов, очевидно, очень ограни ченны. В работе [30] область его приложения была расширена на за дачи для решеток из параллельных пластин конечной толщины. Это развитие метода называется модифицированным методом вычетов.
Решетка из тонких пластин может быть проанализирована также методом Винера — Хопфа. Применение этого метода
взадачах для сканирующих решеток рассмотрено в гл. 4.
7.ИТЕРАЦИОННЫЙ метод
Наша основная задача состоит в том, чтобы решить оператор ное уравнение (1), т. е. для заданного у надо найти оператор К -1, обратный оператору К. Если обратный оператор найден, решение
Методы решения |
119 |
уравнения (1) можно записать в виде
х = К~1у. |
(59) |
Вопрос о существовании обратного оператора выходит за рамки данной книги. Нахождение обратного оператора является трудной задачей и поэтому для ее решения приходится использовать при ближенные методы (например, метод моментов). Кроме приближе ния, связанного с усечением бесконечной матрицы в методе мо ментов, причиной погрешности могут быть ошибки округления при обращении матрицы на ЭВМ. Однако в настоящее время ошиб ки округления можно свести к минимальным. Из-за ошибок, свя занных с усечением бесконечной матрицы и больших затрат време ни на вычисление обратной матрицы, иногда желательно исполь зовать итерационные методы решения линейных уравнений. Пере пишем уравнение (1) в виде
Кх = (1-\-М)х = у или х — у — Мх, |
(60) |
где I — единичный оператор. Предположим, что мы нашли приб лиженное решение х0. Подставив это решение в правую часть урав нения (60), получим новое приближенное решение хг. Мы надеемся, что новое решение дает лучшее приближение к точному решению, чем предыдущее, и что, повторяя процесс, в конце концов придем к точному решению. Такой процесс называется итерационным. Сходимость процесса зависит от свойств оператора М. Достаточ ные условия сходимости будут рассмотрены ниже. Предположим что Мх в некотором смысле мало. Тогда этим слагаемым можно пренебречь и рассматривать у как приближение нулевого поряд ка, т. е.
— У> |
(61) |
Заметим, что нижний индекс обозначает число выполненных ите раций. Подставляя выражение (61) в урвнение (60), получим
х^ |
у |
— Мх0 = у — Му. |
(62) |
Продолжая процесс п раз, |
найдем |
|
|
*п = {/ _ |
М + М*+ .. . + ( - 1Г М) у . |
(63) |
Таким образом, итерационный процесс порождает решение в виде ряда. Если процесс продолжается бесконечно, получаем беско нечный ряд, называемый рядом Неймана. Сразу же возникает вопрос, сходится ли ряд, а если сходится, то действительно ли к точ ному решению? В работах [1, 2, 4] показано, что если || М || < 1, || М || — норма оператора М, то итерационный процесс сходится абсолютно. Необходимо подчеркнуть, что это условие является достаточным, но не необходимым. При выполнении этого условия ряд не только сходится к точному решению, но его сходимость
120 |
Глава 3 |
ие зависит от начального приближения. Конечно, чем ближе пуле вое приближение к точному решению, тем быстрее сходится про цесс. При использовании итерационных методов выбор начального приближения оказывается важным с точки зрения более быстрого нахождения приемлемого по точности решения. Заметим, что выра жение (63) можно вывести, если переписать уравнение (60) в виде х = (I + М)~1у и затем формально разложить (I + М)~г в бино миальный ряд.
Ряд Неймана можно также использовать для улучшения при ближенного решения. Предположим, что приближенное решение уравнения (1) x N получено решением уравнения
У = K Nx N, |
(64) |
где К Лг обозначает А-мерное вырожденное приближение оператора К[ 1 ]. Здесь предполагается, что приближение для уравнения нахо дится в подпространстве, в котором у содержится полностью, так что у = у N. Это выполняется во многих задачах для ФАР и неод нородностей в волноводах.
Если же у ие полностью содержится в подпространстве, исполь зуемом в уравнении (64), ряд необходимо представить в виде
оо
Х = Z ( — i ) q L N q { K N - i y ) при XN->-KN-llJ.
9= 0
Переписав уравнение (1) в виде *)
y= [KN-\-(K — KN)] (х — л;Л-)-]-Аа;,у
ииспользуя уравнение (64), найдем
Хх — £ = |
(7-f-Ljv) 1 LjsrXx, |
(65) |
где L n = Kjl (К — K N). Разлагая множитель (I + |
L jy)-1 в би |
|
номиальный ряд, запишем ряд Неймана в виде |
|
|
* = 2 |
( ~ i ) qb y N. |
(66) |
9= 0 |
|
|
Для сходимости этого ряда необходимо, чтобы || LjV ||< 1 . Это условие также является достаточным, но не необходимым. Из вы ражения (66) следует, что, начав с приближенного решения, можно получить решение первоначальной задачи. Этот метод особенно полезен для повышения точности уже известного приближенного решения. Заменим уравнение (1) выражением
у = КрХр |
(67) |
!) Операторы К и K N должны иметь одинаковую размерность (в данном случае предполагается, что это условие выполнено). Способ расширения размерности K N и x N будет рассмотрен ниже.