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