Доклады БГУИР
2010
№6 (52)
УДК 621.385.6
МНОГОСЕКЦИОННЫЕ ЛБВ НА ВОЛНООБРАЗНО ИЗОГНУТЫХ ПРЯМОУГОЛЬНЫХ ВОЛНОВОДАХ СУБМИЛЛИМЕТРОВОГО ДИАПАЗОНА ДЛИН ВОЛН
А.В. АКСЕНЧИК, И.Ф. КИРИНОВИЧ
Белорусский государственный университет информатики и радиоэлектроники, П. Бровки 6, Минск220013, Беларусь
Поступила в редакцию 26 мая 2010
Описаны методы расчета полей пространственного заряда и полей в зазорах волноводов. Приведены результаты оптимизации одно, двух и трех секционных ЛБВ на волнообразно изогнутых прямоугольных волноводах в диапазоне частот 600-3000 ГГц. Показано, что коэффициент усиления многосекционных ЛБВ может достигать 30-60 дБ. Проведен анализ физических процессов и определен механизм режима самовозбуждения. Отмечено, что полоса частот усиления зависит от длины трубок дрейфа. Показана возможность электронной перестройки полосы усиления.
Ключевые слова: лампа бегущей волны (ЛБВ), субмиллиметровый, многосекционная, оптимизация, волновод.
ЛБВ и ЛОВ на волнообразно изогнутых прямоугольных волноводах являются одними из перспективных приборов для работы в субмиллиметровом диапазоне [1-4].
В работе [3] рассчитаны частотные характеристики (ЧХ) для ЛБВ субмиллиметрового диапазона и проведен подробный анализ влияния различных параметров ЛБВ на вид частотных характеристик. В этой работе отмечается, что даже односекционные ЛБВ с коэффициентом усиления по мощности 10-20 дБ склонны к самовозбуждению на определенных частотах. Двух, трех секционные ЛБВ могут иметь коэффициент усиления 30-60 дБ и также могут самовозбуждаться. Для подавления самовозбуждения таких ЛБВ используются методы, предложенные в работе [4].
На рис. 1. приведена схема трехсекционной ЛБВ на волнообразно изогнутом прямоугольном волноводе.
у: СУ-у
Рис. 1. Схема трехсекционной ЛБВ на волнообразно изогнутом прямоугольном волноводе
Как видно из рис. 1, секции по электромагнитной (ЭМ) волне не связаны, а связаны по электронному потоку. Концы волноводов этих секций нагружены на поглощающие нагрузки. В первой секции электронный поток получает предварительную модуляцию по скорости и плотности. На вход второй секции поступает модулированный поток электронов, возбуждающий в этой секции ЭМ-волну, которая получает дальнейшее усиление и одновременно увели-
Введение
г
чивается модуляция электронного потока по скорости и по плотности. В третьей секции модулированный электронный поток возбуждает ЭМ-волну с нарастающей амплитудой и с выходного конца третьей секции усиленная ЭМ-волна мощностью Pout поступает в нагрузку (антенну, волновод и т.д.). Развязка секций по ЭМ-полю позволяет значительно уменьшить возможность самовозбуждение ЛБВ.
В работе [3] сформулирована математическая модель ЛБВ и ЛОВ на волнообразно изогнутых прямоугольных волноводах: описаны уравнения возбуждения волновода, уравнения возбуждения цепочки эквивалентных четырехполюсников, кратко описаны уравнения движения электронов. Отметим, что адекватность математической модели физическим процессам в ЛБВ зависит от многих факторов, в том числе от точности решения дифференциальных уравнений движения электронов. Анализ сравнения результатов расчетов параметров ЛБВ с экспериментальными данными свидетельствует о том, что точность расчета параметров и режима работы ЛБВ определяет точный расчет сил пространственного заряда в электронном пучке и расчет полей для реальной конфигурации зазора. Поэтому в данной работе в дополнение к математической модели, сформулированной в [3], дано описание методов расчета электрических полей пространственного заряда и электрических полей зазора волновода для реальной конфигурации зазоров.
Нами впервые исследована возможность электронной перестройки полосы усиления многосекционных ЛБВ на волнообразно изогнутых прямоугольных волноводах, а также проведен анализ влияния длины трубок дрейфа на полосу усиления многосекционных ЛБВ.
Математическая модель нерегулярной ЛБВ на волнообразно изогнутом волноводе
Описанный прибор (рис. 1) моделируется цепочкой эквивалентных четырехполюсников, где один четырехполюсник моделирует один изогнутый отрезок волновода. В работе [3] сформулирована математическая модель нерегулярной ЛБВ на волнообразно изогнутом прямоугольном волноводе с учетом потерь и диэлектрического заполнения волновода. Здесь рассмотрим уравнения движения электронов.
В математической модели используется система координат, совпадающая с общепринятой системой координат для описания полей в волноводе. Учитывая, что электронный поток движется перпендикулярно и посредине широкой стенки волновода, принимаем за продольную координату, совпадающую с направлением движения электронов, координату У. Все продольные электрические поля будем записывать относительно координаты У. На каждый электрон действуют поля рядом расположенных электронов (поля пространственного заряда) и поля в зазорах волноводов.
Поля пространственного заряда. В данном приборе для моделирования электронного потока применяется метод крупных частиц. Электронный поток цилиндрической формы может быть сплошной или полый. Соответственно и форма крупных частиц выбирается в виде дисков или колец конечной толщины по У (рис. 2).
Рис. 2. Периодизация полей заряженных частиц Принимаем равномерное распределение плотности заряда кольца по объему кольца, тогда радиус ведущего центра кольца г0 определяется так: г = ^(г12 + г22)/ 2 . Поля пространственного заряда, электрические поля крупных частиц можно получить, если выбрать в качестве модели бесконечно тонкое кольцо и применить методику усреднения полей частицы по радиусам точки наблюдения и точки источника, изложенную В.А. Солнцевым в [5]. После со-
ответствующих преобразований продольную составляющую электрического поля Еу крупной частицы можно представить следующим выражением:
Еу =
р 6
2реа
г 302(п _ „гм
_а_е ' аЬ
2 г2, е
х =1
3ГК)
(1)
где 3о(по* -) =
;(г0/а) да
й г г г и '
к2(к —) - — 3\ (к ~)ъ; У , Г — координаты точки наблюдения;
а ы
у, г0 — координаты центра масс источника; а — радиус трубки дрейфа; Р 6 — заряд частицы;
Ь = дД - Ь2 — релятивистская поправка, учитывающая изменение поля Еу при переходе от движущейся системы координат к неподвижной.
Для ускорения времени счета по одномерной модели (исключается экстраполяция электронов вперед и назад) вводятся периодизированные поля пространственного заряда. Учитывая то, что сгустки электронов формируются в трубке дрейфа на расстоянии электронной длины волны I = 10Ь0, периодизация позволяет учесть влияние источников поля пространственного заряда, отстоящих друг от друга на расстояние I .
В данной модели периодизация проводится численно. Процесс ее проведения поясняет рис. 2. Поле в точках от I до 21 (рис. 2) определяется суммированием полей от 4 источников, расположенных на расстоянии I друг от друга. Штриховой линией показано периодизирован-ное поле, рассчитанное с учетом других, периодически расположенных колец справа и слева от конкретного кольца, например расположенного на расстоянии I . Ограничение 4 источниками при расчете периодизированного поля вызвано тем, что полями от источников, расположенных на расстояниях, больших 21 , от точки наблюдения в виду их экспоненциального ослабления, можно пренебречь. Для сокращения времени расчета полей во время интегрирования уравнений движения (чтобы не рассчитывать каждый раз поле по формуле (1 ) при изменении координаты у) составляются таблицы полей пространственного заряда Е р . Таблицы полей Е р составляются для заданного г0, Г, Г и ряда фиксированных значений \у - у | с применением формулы
(1 ). Для использования их при численном решении уравнений движения применяется линейная интерполяция.
Квазистатические электрические поля зазора в волноводе. Расчет поля зазора в волноводе сводится к решению электростатической задачи о распределении потенциала в некотором объеме сложной формы с заданными граничными условиями-решению уравнения Лапласа. В данном случае объем сложной формы представляет собой отрезок волновода с трубками дрейфа, расположенными посредине широкой стенки и образующими между собой внутри волновода некоторый зазор. Учтем, что трубки дрейфа и небольшую часть окружающего их объема волновода можно рассматривать как аксиально-симметричную систему, найдем решение уравнения Лапласа в цилиндрической системе координат. Для этого воспользуемся численными методами и получим решение в узлах некоторой сетки, построенной в этой области (рис. 2). Применяя метод сеток, уравнение Лапласа заменим приближенно уравнением в конечных разностях [6], которое для осесимметричной задачи в цилиндрической системе координат имеет вид
2(г -1/2/Г)
)
и
з -1,;
и
2(Г + 1 / 2 / Г)
К (\
)
и
з + и
+
К. К
и
2
(2)
з,; -1
К К.., + К)
и
з, •+1
-и =о
з, •
где К = г. , - г.; К = у.,, - у.; г.,,,, = — (г + г.,,).
^ г з- 1 з у, ^ 1 1 ± 1/2 2
При г = 0 вместо (5) получаем
а
п
1 -1
2
2
у,-1 у
---и. ,. —— и.. +---и.+1. +
к (Иг + Иг) 1-1-' к нг 1■' Нг (Иг + Иг) 1 г1-1 г1-1 Г1 г1-1 Г1 Г1 Г1 -1 Г1
1 11 (
+-и.. , +-и. .+1--и. . = 0.
к (К + к ) 1'1 к (К + к ) 1'1 к к 1 1
Пу,-1 1 + ПУ, ) к 1 + ПУ, ) /V,
Таким образом, задача сводится к двумерной задаче расчета распределения потенциала в области, ограниченной контуром сложной формы с заданными граничными условиями. С учетом того, что трубки дрейфа могут выступать внутрь волновода, на рис. 3 изображена область зазора (шириной $) в волноводе с такими трубками дрейфа. Так как система симметрична относительно оси У, то достаточно найти решение уравнения Лапласа для положительных г.
Рис. 3. Область зазора в волноводе: а — радиус трубки дрейфа; b — ширина узкой стенки волновода; с=5а; d — длина зазора в волноводе (может совпадать с b),.
Граничные условия задаются следующим образом (рис. 3): на металлической поверхности левой трубки дрейфа 4 и стенке 2 потенциал полагаем равным +j 0 = const, на металлической поверхности правой трубки дрейфа 5 и стенке 3 потенциал полагаем равным - j 0. Учитывая тот факт, что на достаточном расстоянии (2а) от края зазора потенциал внутри трубок дрейфа постоянен, можно ограничить длину трубок дрейфа, поставив металлические поверхности 6 и 7 и принимая потенциал на них равным потенциалу соответствующей трубки. Граничные условия на поверхности 1 задаются следующим образом: принимается линейное изменение потенциала вдоль поверхности 1 от +j 0 до - j 0 пропорционально шагу сетки вдоль оси Y.
Расположение поверхности 1 выбирается из условия малого влияния изменения объема волновода на поле в зазоре и трубках дрейфа.
Для решения разностных уравнений (2), (3) применяем метод последовательной верхней релаксации [6]. Данный метод реализован в виде программы решения уравнения Лапласа, аппроксимированного пятиточечными разностными уравнениями (2), (3). Задаваемая точность решения порядка 10-5 достигалась за 30-60 итераций в зависимости от формы задаваемой области (рис. 3).
Для данной одномерной модели необходима только продольная составляющая электрического поля Ey, которая рассчитывается путем численного дифференцирования полученного распределения потенциала на известной сетке по следующей формуле:
Ey,= - ^ = Uj+1 - U +2j+l)/ & + dЩ (4)
С использованием формулы (4) составляется одномерная таблица распределения поля Ey в зазоре, и для ее использования при численном интегрировании уравнений движения применяется линейная интерполяция.
Уравнения движения электронов. После того как все ЭМ-поля, действующие на электроны, определены запишем уравнение движения электрона без учета торможения излучением (до v0 / c < < 0,99) в одномерном приближении:
£=- <-01К)2 Г * о
где <0 = е / т0; е, т0 — соответственно заряд и масса покоя электрона; V — продольная скорость электрона; с — скорость света в вакууме; v0 — начальная скорость электронного потока.
Напряженность электрического поля действующего на электроны представим как = Е&0 + , где = Re(£meJ (ш+3)) — напряженность продольной составляющей ВЧ поля
в зазоре волновода (4); Е& — напряженность продольной составляющей поля пространственного заряда. Поле пространственного заряда Е& рассчитываем исходя из моделирования электронного потока N6 крупными частицами (расположенными на электронной длине волны / = 10Ъ0) по описанной выше методике. Введем следующие
V - 1/2
безразмерные переменные: Т = у/ Ь,3 = Ш,g0 = |1 - (—)2Ц ,
V = \ и, = - шу / У0,30= —, gi = (ЭД , Л = С1, , (6)
V. V- и С ш
С
- 1/ 2
, = I ; х = е = & = е1 О
а Хк С т т , и т Е ± тпа, р 2 2 ,
2рЪ 0 ё т и 0 wv т 0^0Т
Ток = У1/ Ь Т1к = У 2 / Ь ,
где у — продольная координата электрона, Ь — длина ЛБВ, \к — длина отрезка волновода; г — радиус электронного луча, d=y2гy^ — длина зазора, ш — опорная частота, хк - относительное напряжение на зазоре, полученное с учетом прямых и обратных волн, возбуждаемых электронным потоком в других зазорах, и рассчитанные по специальным алгоритмам [3].
После нормировки уравнения (5) получим релятивистское уравнение движения 7-й крупной частицы:
^ = "4"[- *а Re(Xkej(и'+3оТ+3к>)е°(Т - Ток) + s^] (7)
% = ^-1), (8)
аТ V,
1 т ,1 и, - и,
ер(
где ^ = — с е (—'-—^^(и - и); е°(Т - Т0) — весовая функция пространственного
Р ^е , = 1 Р Р ,
распределения поля зазора - рассчитывалась методом сеток; е (у) — весовая функция распределения напряженности поля пространственного заряда, определялась с использованием пери-одизированных полей крупных частиц (дисков, колец) и подробно описана выше.
Уравнения движения (7), (8) и уравнения возбуждения волновода [3] образуют самосогласованную нелинейную систему уравнений ЛБВ, ЛОВ на волнообразно изогнутом прямоугольном волноводе. После определения наведенных в зазорах волновода токов используется алгоритм возбуждения цепочки четырехполюсников, описанный в [3].
Результаты расчетов многосекционных ЛБВ
Для проведения расчетов многосекционных ЛБВ по описанной математической модели составлена программа. Проведены расчеты односекционной регулярной ЛБВ для работы на длине волны /0 =0,01 см (вариант А1). Ускоряющее напряжение ?У0=12 кВ, ток электронного луча /0=0,15 А. Размеры волновода а=0,006 см, 6=0,0005 см, эквивалентное сопротивление вол-
новода 7^=113,67 Ом, радиус трубки дрейфа Л=0,0002 см, радиус электронного луча г=0,00016 см, длина зазора ё=Ъ. Волновод изготовлен из меди с проводимостью х = 57 Ч106 См/м. В результате оптимизации параметров ЛБВ получен коэффициент усиления по мощности К = 10^(РОМ( / Р1Я )=13,6 дБ, электронный КПД К =0,0006, выходная мощность ЛБВ Рои=0,62 Вт, входная мощность =0,027 Вт. Коэффициент фазы эквивалентного четырехполюсника 2,615 рад. Число зазоров (изогнутых отрезков волновода) N=42.
Для этого варианта ЛБВ характерно устойчивое самовозбуждение [3], на относительных частотах от 1,08 до 1,11 наблюдается "полоса генерации". В результате самовозбуждения на входе прибора появляются электромагнитное (ЭМ) излучение встречной волны с мощностью 18,7 Вт. Мощность ЭМ-волны на выходе прибора равна 4,7 Вт. Для устранения самовозбуждения воспользуемся методиками работы [4].
Уменьшение оптимальной длины трубок дрейфа на 10% (частотная характеристика (ЧХ) сдвигается влево по оси частот) и уменьшение и0 до 11,5 кВ (ЧХ сдвигается вправо по оси частот) устраняет генерацию. ЧХ данного варианта-зависимость коэффициента усиления по мощности Кр от относительной частоты w = // /0 приведена на рис. 4, кривая 1. "Полосы генерации" нет.
Рис. 4. Зависимость Кр от относительной частоты w, I =0,01 см: кривая 1 — односекционная ЛБВ; кривая 2 — двухсекционная ЛБВ, Р,-и=0,0027-10~5 Вт; трехсекционная ЛБВ вариант В1, и0=12 кВ: кривая 3 — входная мощность Р,-и=0,27-10-5 Вт; кривая 4 — Р,-и=0,27-10-7 Вт
Анализ физических процессов, приводящих к самовозбуждению, показал, что имеет место самовозбуждение ЛБВ за счет положительной обратной связи, создаваемой обратной волной в волноводе. Выполняется условие для баланса фаз при самовозбуждении колебаний в ЛБВ [3] ЬРЬц + КЬб = 2рт , где ЬРп — постоянная распространения прямой ЭМ-волны в замедляющей системе; Ь — длина системы; Ь^об и Ьоб — соответственно постоянная распространения обратной ЭМ-волны и длина системы по цепи обратной связи.
Разделение ЛБВ на две или три секции, не связанных по ЭМ-волне, уменьшает склонность прибора к самовозбуждению и позволяет значительно увеличить коэффициент усиления. Проведены расчеты двух- и трехсекционных ЛБВ с большим коэффициентом усиления на длине волны 10 =0,01 см. Параметры луча и размеры трубок дрейфа и волновода такие же, как и в варианте А1: ускоряющее напряжение Цо=12 кВ, ток электронного луча /0=0,15 А, число зазоров в секциях одинаково, по N=42. Для первой секции выбран коэффициент фазы четырехполюсника К =2,615 рад, для второй и третьей секций — коэффициенты фаз К. =1,3 рад. Это
позволяет разнести по частоте "полосы генераций" в каждой секции.
На рис. 4 приведены ЧХ двухсекционной ЛБВ, кривая 2 — входная мощность /'„,=0.0027 ■ 10 ' Вт. На относительной частоте и—0,94 (рис. 4, кривая 2) коэффициент усиления по мощности Кэтого варианта двухсекционной ЛБВ равен 38 дБ, электронный КПД
К =0,0134, выходная мощность ЛБВ ^^=17 Вт при входной мощности Р(и =0,0027 Вт. Полоса частот усиления составляет ® 10%. В полосе относительных частот 1,08-1,1 наблюдалось самовозбуждение прибора. Однако при уменьшении входной мощности до 0.0027 -10' Вт амплитуда сигнала на этих частотах отсутствует, рис. 4, кривая 2. Следовательно, прибор работает в режиме регенеративного усиления.
Кр [сШ]
0,135 0,901 0,967 1,033 1,099 1,165
На рис. 4 (кривые 3, 4) приведены частотные характеристики трехсекционной ЛБВ. На относительной частоте ^=0,934 коэффициент усиления по мощности Кр трехсекционной
ЛБВ (вариант В1, рис. 4, кривая 3) равен 62 дБ, электронный КПД К =0,0071, выходная мощность ЛБВ Рои=4,1 Вт при входной мощности р =0,0000027 Вт. Полоса частот усиления составляет «15%.
В полосе относительных частот 1,08-1,11 наблюдается самовозбуждение прибора, рис. 4, кривая 3. Анализ режима работы ЛБВ на этих частотах показал, что самовозбуждение возникает в 1-й секции. При уменьшении входной мощности до 0,27-10 7 Вт амплитуда сигнала на этих частотах отсутствует (рис. 4, кривая 4) следовательно, прибор работает в режиме регенеративного усиления.
Расчеты показали, что ЛБВ с длинными трубками дрейфа (в трубке дрейфа угол пролета электронов более 10 рад) имеют более узкую полосу усиления, чем с короткими (угол пролета электронов менее 6 рад).
На рис. 5 (кривая 1) приведена ЧХ трехсекционной ЛБВ (вариант В2) с короткими длинами труб дрейфа (угол пролета электронов 4,8 рад), рассчитанной для работы на длине волны 10 =0,05 см. Ускоряющее напряжение и0=12 кВ, ток электронного луча /0=0,1 А. Оптимальная длина трубки дрейфа между смежными зазорами короткая (угол пролета электронов менее 6 рад).
Рис. 5. Зависимость Кр от относительной частоты w: кривая 1 — вариант В2 длина волны 10 =0,05 см, и0=12 кВ, входная мощность Р,-я=0,27-10~5 Вт; кривая 2 — и0=10,4 кВ, Ри=0,27-10~5 Вт; кривая 3 — и0=10,4 кВ, Рш=0,27-10~7 Вт
Для устранения самовозбуждения на частотах 1,23...1,25 согласно методике работы [4], уменьшим оптимальную длину трубок дрейфа на 10% (ЧХ сдвигается влево по оси частот) и уменьшение и0 до 10,4 кВ (ЧХ сдвигается вправо по оси частот) устраняет генерацию. ЧХ этого варианта приведена на рис. 5, кривая 2. На частотах 1,2. 1,23 есть полоса регенеративного усиления, так как уменьшение входной мощности в 100 раз устраняет сигнал на этих частотах, (рис. 5, кривая 3). Полоса частот усиления этого варианта более 15-20%.
На рис. 6 приведена ЧХ варианта двухсекционной ЛБВ с длинными трубками дрейфа (угол пролета электронов 11,8 рад) с ускоряющим напряжением 11,,= 12 кВ (кривая 2). Полоса частот усиления составляет «6%. Рис. 6 демонстрирует влияние изменения ускоряющего напряжения на сдвиг ЧХ этого варианта. Отметим, что при изменении и0 от 11,8 кВ до 12,8 кВ частотная характеристика ЛБВ сдвигается вправо по оси частот в значительных пределах-почти на 10%. Такую ЛБВ можно использовать в селекторах частоты с электронной перестройкой частоты.
Рис. 6. Зависимость ^ от относительной частоты w, двухсекционная ЛБВ, 10 =0,05 см: кривая 1 — (70=11,8 кВ; кривая 2 — ^=12 кВ; кривая 3 — ^=12,4 кВ; кривая 4 — ^=12,8 кВ;
Заключение
Проведена оптимизация и рассчитаны частотные характеристики многосекционных ЛБВ. Получены варианты приборов с большим коэффициентом усиления — 30-62 дБ, выходной мощностью 5-18 Вт в диапазоне частот 600-3000 ГГц. Проведен анализ физических процессов, в частности процесса самовозбуждения. Одним из методов подавления самовозбуждения является разделение прибора с большим усилением на две, три секции и применением секций в таком гибридном приборе с разными значениями коэффициентов фаз четырехполюсников. Показана и эффективность применения других методов устранения самовозбуждения, предложенных в [4]. В многосекционных ЛБВ самовозбуждение чаще возникает в первой секции.
Установлено, что при больших длинах труб дрейфа (угол пролета электронов >10 рад) полоса частот усиления ЛБВ достаточно узкая «6%, и такую ЛБВ можно использовать как эффективный селектор частоты с электронной перестройкой частоты. При меньших длинах труб дрейфа (угол пролета электронов менее 6 рад) полоса частот усиления увеличивается в 3 раза.
MULTISECTION TWT ON WAVY BENT RECTANGULAR WAVEGUIDES OF THE SUBMILLIMETER RANGE OF LENGTHS OF WAVES
A.V. AKSENCHYK, IF. KIRINOVICH
Abstract
Methods of calculation of fields of a space charge and fields in gaps of waveguides are described. Results of optimization one, two and three section TWT on wavy bent rectangular waveguides in a range of frequencies of 600-3000 GHz are resulted. It is shown, that the gain multisection TWT can reach 30-60 dB. The analysis of physical processes in such a TWT is carried out, the mechanism of a mode of self-excitation is defined. It is noticed, that the band of frequencies of gain depends on length of tubes of drift. Possibility of electron tuning of the band gain is shown.
Литература
1. Seong-Tae Han, Jung-I1 Kim, Gun-SikPark // 3rd IEEE Int. Vacuum Electronics Conf. 2002. DoubleTree Hotel Monterey. California. April 23-25. 2002. P. 94-95.
2. AksenchykA. V., Kurayev A.A. // 6th International Kharkov Symposium "Physics and Engineering of Millimeter and Submillimeter Waves (MSMW'07) / Kharkov. Ukraine. June 25-30. 2007. P. 266-268.
3. АксенчикA.B., Кураев А.А., Киринович И.Ф. // Весщ НАН Беларуси Сер. фiз.-тэхн. навук. 2009. № 3. С. 113-124.
4. АксенчикA.В., Кураев А.А. Киринович И.Ф. // Весщ НАН Беларуи. Сер. фiз.-тэхн. навук. 2009. № 4. С. 102-110.
5. СолнцевВ.А. // Лекции по электронике СВЧ. Саратов. 1972. Кн. 1. С. 46-89.
6. Ильин В.П. Численные методы решения задач электрооптики. М., 1974.