УДК 517.98
DOI 10.21685/2072-3040-2020-3-10
И. В. Бойков, П. В. Айкашев
ПРИМЕНЕНИЕ НЕПРЕРЫВНОГО ОПЕРАТОРНОГО МЕТОДА К РЕШЕНИЮ УРАВНЕНИЙ ПОКЛИНГТОНА И ГАЛЛЕНА ДЛЯ ТОНКИХ ПРОВОЛОЧНЫХ АНТЕНН
Аннотация.
Актуальность и цели. Одной из центральных задач в электронике сверхвысоких частот является построение миниатюрных антенн, обладающих высокими характеристиками. Основными уравнениями, используемыми при моделировании проволочных антенн различной конфигурации, являются уравнения Поклингтона, Галлена, сингулярные и гиперсингулярные интегральные уравнения. При численном решении уравнений Поклингтона и Галлена в основном используются методы моментов и Галеркина. Так как уравнения Поклинг-тона и Галлена относятся к классу некорректных задач, то при реализации методов моментов и Галеркина возникают дополнительные трудности, связанные с неустойчивостью вычислительных схем. В данной работе для решения уравнений Поклингтона и Галлена предлагается применить непрерывный метод решения операторных уравнений, обладающий эффектом регуляризации. Этот эффект обусловлен тем, что непрерывный метод решения операторных уравнений построен на основе Ляпуновской теории устойчивости решений дифференциальных уравнений.
Материалы и методы. Исследуются приближенные методы решения уравнений Поклингтона и Галлена. Метод построения вычислительных схем заключается в следующем. Исходные уравнения аппроксимируются системой линейных алгебраических уравнений, построенных по технологии метода сплайн-коллокации. Система линейных алгебраических уравнений решается непрерывным операторным методом. В качестве достоинств предложенного метода следует отметить: его устойчивость при возмущении ядер уравнений и правых частей и возможность при построении вычислительной схемы учесть граничные условия на концах вибратора.
Результаты. Построены новые устойчивые численные методы решения уравнений Поклингтона и Галлена. Эффективность предложенных методов продемонстрирована решением модельных примеров.
Выводы. Предложенный в работе метод построения и обоснования сходимости вычислительных схем может быть распространен на уравнения, моделирующие различные модификации антенн и являющиеся аналогами уравнений типа Поклингтона и Галлена.
Ключевые слова: уравнения Поклингтона и Галлена, метод коллокаций, вибраторы.
I. V. Boykov, P. V. Aykashev
APPLICATION OF THE CONTINUOUS OPERATOR METHOD TO THE SOLUTION OF THE POCKLINGTON AND GALLEN EQUATIONS FOR THIN WIRE ANTENNAS
© Бойков И. В., Айкашев П. В., 2020. Данная статья доступна по условиям всемирной лицензии Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), которая дает разрешение на неограниченное использование, копирование на любые носители при условии указания авторства, источника и ссылки на лицензию Creative Commons, а также изменений, если таковые имеют место.
Abstract.
Background. One of the central tasks in microwave electronics is the construction of miniature antennas with high performance. The main equations used in modeling wire antennas of various configurations are the Pocklington, Gallen equations, singular and hypersingular integral equations. In the numerical solution of the Pocklington and Gallen equations, the methods of moments and Galerkin are mainly used. Since the Pocklington and Gallen equations belong to the class of ill-posed problems, when implementing the methods of moments and Galerkin, additional difficulties arise due to the instability of computational schemes. In this paper, to solve the Pocklington and Gallen equations, it is proposed to apply a continuous method for solving operator equations, which has the effect of regularization. This effect is due to the fact that the continuous method for solving operator equations is based on the Lyapunov theory of stability of solutions of differential equations.
Materials and methods. The article investigates approximate methods for solving the Pocklington and Gallen equations. The method for constructing computational schemes is as follows. The initial equations are approximated by a system of linear algebraic equations constructed using the spline collocation method. The system of linear algebraic equations is solved by a continuous operator method. The advantages of the proposed method should be noted: 1) its stability under perturbation of the kernels of equations and right-hand sides; 2) the construction of a computational scheme allows taking into account the boundary conditions at the ends of the vibrator.
Results. New stable numerical methods for solving the Pocklington and Gallen equations are constructed. The effectiveness of the proposed methods is demonstrated by solving model examples.
Conclusions. The proposed method for constructing and justifying the convergence of computational schemes can be extended to equations that model various antenna modifications and are analogous to equations of the Pocklington and Gallen
type.
Keywords: Pocklington and Gallen equations, collocation method, vibrators.
Введение
Более ста лет одним из основных математических аппаратов, применяемых при моделировании электромагнитных процессов в электрическом вибраторе, являются интегродифференциальные уравнения Поклингтона и Хар-рингтона и интегральное уравнение Галлена [1]. В последнее время к исследованию этих процессов привлекаются также сингулярные интегральные уравнения [2, 3] и гиперсингулярные интегральные уравнения [4].
В качестве численных методов для решения уравнений Поклингтона и Галлена привлекается, как правило, метод моментов и его модификации [5].
В отличие от известных работ, в которых упомянутые уравнения решаются методом моментов, в данной работе используется непрерывный метод решения операторных уравнений. Преимущества этого метода заключаются в его устойчивости к возмущениям ядер уравнений и правых частей. Уравнения Поклингтона и Галлена являются интегральными уравнениями Фредгольма первого рода и их решения являются некорректными задачами. Поэтому целесообразно решать упомянутые уравнения непрерывным методом решения операторных уравнений, который обладает регуляризующими свойствами.
Статья посвящена приближенному решению уравнений Поклингтона и Галлена, моделирующих тонкие изогнутые проволочные антенны.
В работе предложены новые численные методы и дано их обоснование. Исследуются вопросы быстроты сходимости этих методов при решении уравнений Поклингтона и Галлена.
Статья построена следующим образом. В первом разделе представлен непрерывный метод решения операторных уравнений. Во втором разделе даны постановка задачи и краткий обзор литературы. Третий раздел посвящен приближенному решению уравнения Поклингтона. Четвертый раздел посвящен приближенному решению уравнения Галлена. В пятом разделе приведены решения модельных задач.
1. Непрерывный метод решения операторных уравнений
Для решения уравнений Поклингтона и Галлена в статье используется непрерывный метод решения нелинейных операторных уравнений [6]. Приведем здесь его краткое описание.
Пусть X — банахово пространство; К — оператор, действующий из X в X; B(a,r) = {x,aе X:|| x — a ||< r}; S(a,r) = {x,ae X:|| x — a ||= r}; Л(К) -логарифмическая норма линейного оператора К, определяемая [7] выражением Л(К) = limh^o(|| I + hK || —1) / h, где символ h 0 означает, что h стремится к нулю, убывая.
Для матриц в часто используемых пространствах логарифмические нормы известны.
Пусть дана вещественная матрица A = {a^}, i,j = 1,2,...,n, в n-мерном пространстве Rn векторов x = (xj,..., xn) с нормой
n n
H x |ll= 2 1 xk W x ^2 = [2 1 xk |2]1/2, || x | 3= max 1 xk |. k=1 k=1 1<k<n
Логарифмическая норма матрицы A равна [8]:
f n >
Л1(A) = max a^ + 2 I aj
j I i=1,i*j
Л2(A) = Xmax ((A + AT )/2),
f
Л3 (A) = max
aii + 2 1 aij
}=1, &
Здесь Xтах((А + Ат)/2)- наибольшее собственное значение матрицы (А + Ат) / 2.
Рассматривается нелинейное операторное уравнение
А( х) - / = 0, (1)
где А: X ^ X - нелинейный оператор, отображающий банахово пространство X само в себя.
Известия высших учебных заведений. Поволжский регион Поставим уравнению (1) в соответствие задачу Коши:
^ = А(х(г))-(2)
х(0) = хо. (3)
Справедливо следующее утверждение.
*
Теорема 1 [6]. Пусть задача (2)-(3) имеет решение х , и на любой
*
дифференцируемой кривой g(г), расположенной в шаре В(х , г), выполняются следующие условия:
г
1) при любом г (г > 0 ) выполняется неравенство |л( А (g (т)))т < 0;
о
1 Г '
2) справедливо неравенство lim - |Л(А (g(т)))^т < -аg, аg > 0.
г^те г
о
Тогда решение задачи Коши (2)-(3) при г, стремящемся к беско-
*
нечности, сходится к решению х уравнения (1).
2. Уравнения Поклингтона и Галлена
Интегральное уравнение электрического вибратора описывает связь между электрическим полем и распределением тока по проводнику. Оно получено из уравнений Максвелла и имеет следующий вид [9]:
г те 2 2 е-Л1г-г'1
\г) = -^ Г/(г')[к2 + У2]-^--М. (4)
к •>, 4п | г - г |
V
Et/
Здесь Е1 (г)(г = (Х1, Х2, Х3)) - вектор электрического поля в точке наблюдения г , I(г') (г ' = (ху, х^, Х3')), - плотность текущего распределения
/ 2 2
тока (искомая величина) в точке источника г , [к + У ], - оператор Гельмгольца, функция Грина которого равна
( 3 >
3
e~ik\r-г'\ / | г _ /
Г _ r'\=TIXj _ XJ
V
]=1
к и п- соответственно волновое число и внутренний импеданс в свободном пространстве: к = 0)^Дое(),П = VДо / е0 , где ю - угловая частота, До и £о -магнитная и диэлектрическая проницаемость соответственно.
Интегральное уравнение Поклингтона получено из уравнения (4) и граничного условия для тонкой проволоки произвольной формы. Граничное условие для совершенного проводника определяется из условия, что полное электрическое поле исчезает на его поверхности. Тогда можно считать, что
точка г лежит на поверхности проводника и Е1 (г5) = Е1 (г5) + Е5 (г5) = 0, где
Е1 (г*) и Е* (г*) - наложенное и рассеянное электрическое поле соответственно. Для высоких частот благодаря скин-эффекту ток может быть расположен только на поверхности провода. Следовательно, наложенное электрическое поле имеет вид
■_ —г£\г—г'\
ГГ7 „ чг/21т72п е
E(Г = ^) = -J Я^(Г ')(r ' = ^ )[k2 +V2]
s' ф
4п|r-r
—ad ф ds ,
(5)
где переменная 5 представляет длину дуги по проводу; ф — азимутальным угол вокруг поперечного сечения провода. Текущее распределение тока трансверсально постоянно, учитывая, что радиус проволоки а очень мал по сравнению с длиной волны.
В результате получена [10] следующая модификация уравнения Поклингтона:
E*
(s) = --L j>s (s') юе J
k 2( s • s') +
dsds'
e-ik\r-r\
4п \ r - r \
(6)
где Е* - тангенциальное наложенное электрическое поле. Учитывая скин-эффект в тонкой проволоке, в уравнении (6) электрическое поле можно выразить в виде криволитейного интеграла по длине дуги .
Уравнение (6) можно использовать при любой топологии тонкой проволочной антенны.
Из уравнения (6) следует, что тангенциальное наложенное электрическое поле на поверхности совершенного проводника выражается гиперсингулярным интегральным уравнением.
В уравнении (5) оператор V действует по переменной 5 . Таким образом, уравнение (5) эквивалентно следующему:
~к2 + V2 ] Я7* (Гх/ = Г*)/-1*'Г Г|, adф ж', (7)
E (r = rs ) = ■
in
4п\r-r
s ф
Классическое уравнение Поклингтона прямолинейного вибратора имеет вид [1, с. 222]:
,2 > I
+ — IJ Iz (z )G (a, z - z' )dz' = Az (a, z),
(8)
'-l
где Лг (а, г) - продольная составляющая векторного электродинамического потенциала для электрического тока 12 (г); 21 - длина вибратора; а - радиус вибратора,
1
G(a, z - z') = ■
—exp{-ikRa(z - z,)},
4nRa (z - z') Ra (z - z') = \j(z - z')2 + a2, k = ю/eja /<
£ и ц - диэлектрическая и магнитная проницаемость среды, окружающая вибратор; ю - угловая частота; с - скорость света.
Это уравнение соответствует модели тонкого вибратора, работающего на малых частотах. В этом случае ток распределяется по всему поперечному сечению провода и моделируется бесконечно тонкой линией, проходящей через центр цилиндра (провода) и точку разрыва.
Следует отметить, что классическое уравнение Поклингтона описывается уравнением Фредгольма первого рода с гладким ядром [1, 11] и, следовательно, является некорректной задачей.
Приближенному решению уравнения (8) посвящено много работ, в которых в большинстве своем применяются методы Галеркина и моментов [5, 10, 12].
В работе [12] метод моментов применяется для решения уравнения Поклингтона линейной антенны:
l
i1 (z')
-l
£ + * 2
-ikRa ( z-z')
z = iroeEz (z), (9)
4nRa (z - z')
где Ra (z - z') = yja2 + (z - z')2 .
Замечание. Это уравнение может быть получено из (7) коммутацией операторов дифференцирования и интегрирования.
Отметим некоторые модификации уравнения (7). Так как электромагнитное поле моделируется бесконечно тонкой линией, проходящей по центру вибратора и точку разрыва, в уравнении (7) при очень малом радиусе провода естественно перейти к пределу при а ^ 0. В результате получаем уравнение [4]:
k2 +-
2 ^ 1 l —ik\z-z'\
z 2
ilz (z')--rrdz' = 4 (z), (10)
4n Ji | z - z \
которое является гиперсингулярным интегральным уравнением.
Таким образом, гиперсингулярные интегральные уравнения естественно возникают как при моделировании антенн, работающих на высоких частотах, так и при моделировании очень тонких проволочных антенн.
Приближенные методы решения уравнения (10) исследовались в работе [4], где было предложено несколько методов. В частности, был предложен следующий метод. Так как на торцах вибратора ток 12 (I) = 12 (-1) = 0, то
уравнение (10) было продолжено на числовую ось и к нему было применено преобразование Фурье. Полученное в спектральной области уравнение решалось итерационным методом.
В работе [2] показано, что физически корректной является математическая модель произвольного вибратора, представленная следующим сингулярным интегральным уравнением первого рода:
1 Г1 Ь^! & + Г1 4 (/)К (2, / = 4/Ю£0£ЕсГ (2 ),
где
K (z, z' ) = -¡- feih(z'-z W (h )dh, Ag(h) = Jo (-Wh2 - Y2H)2) (-/Wh2 -Y2 ]) + g(h),
(2)
Jo и H0 - функции Бесселя первого рода и Ханкеля второго рода соответственно.
Приближенные методы решения гиперсингулярных интегральных уравнений первого рода на различных классах функций исследовались в работах [13-17].
Наряду с уравнением Поклингтона при исследовании вибраторов часто используется уравнение Галлена
l . 2
ГIz (z')G(a, z - z')dz = A cos pz + B sin pz -—sin p | z (11)
J w
-l
где А и В - произвольные константы, определяемые из граничных условий обращения тока в нуль на концах вибратора; w = в / (юе).
При стремлении a к нулю уравнение (11) трансформируется в гиперсингулярное интегральное уравнение
1 /2nv .
w
ГIz (z')G(z - z')dz = A cos pz + B sin pz -—sin p | z (12)
w
-l
где G(z - z') = ——1——exp{-/k | z - z |}. 4п | z - z |
3. Приближенное решение уравнения Поклингтона
При исследовании излучения прямого вибратора исходным может служить уравнение Поклингтона
( d Л l ~/fRa (z - z')
[ ЬM*=>
где Ect - стороннее электрическое поле; Rz -z') = ^(z -z )2 + a2, 2 2
Y = k ец, k = ю / c - волновое число; е, ц - диэлектрическая и магнитная проницаемость пространства, окружающего вибратор, a - толщина вибратора.
Первая вычислительная схема
Для простоты описания вычислительной схемы введем обозначение
f (z) = -/roeoeEct (z) и рассмотрим уравнение
Л i -/fRa (z - z')
+ Y2 I Г Iz (z )-dz = f (z), -1 < z <l. (13)
JJ-l zV UnRa (z - z') W
v dz2
Известия высших учебных заведений. Поволжский регион Введем обозначение
( г—г')
U(z) = \jz (Z)
e
4nRa (z - z')
Подставив U(z) в уравнение (13), приходим к дифференциальному уравнению
' U(z) = f(z). (14)
V dz2
+ Y2
Присоединим к уравнению (14) граничные условия
U (—l ) = abU (l ) = a2, (15)
где a1 и - произвольные вещественные числа, которые в дальнейшем подлежат определению.
Решение уравнения (14) имеет вид
U (z) = F (z) + A cos (yz) + B sin (yz), (16)
где F (z) — частное его решение.
Подставив в уравнение (16) z = —l, z = l и воспользовавшись граничными условиями (15), приходим к системе уравнений
a1 = F (— l) + A cos (yl) — B sin (yl),
a2 = F (l) + A cos (yl) + B sin (yl).
Для простоты выкладок предположим, что рассматриваемый диполь симметричный и Iz (— l) = Iz (l), U (— l) = U (l). Тогда
U (z) = F (z)+ a1 — a2 — F H) — F (l) cos (yz) + a2 — a1 — F () + F (—l) sin(yz) W W 2cos(yl) 2sin(yl) U '
в предположении, что cos (yl) ^ 0, sin (yl0 .
Для определенности положим в уравнении (14) f (z) = 4iro£o£sin(roz), т.е. полагаем Ect(z) = —4sinroz . Тогда
, . 4iroeo£sin (rot)
F (z) =-2-2-.
у — ro
Таким образом, задача свелась к решению уравнения
e~iyR ' ^ 4iroe0esin(roz)
b')
/ч e
-dz =-^-^—1 +
4nR y2 —ro2
a1 — a2 — 4iro£o£sin (rol) + 4iro£o£sin (ro(—l))
+--—-—--——— cos (yz) +
2cos(yl) v '
a2 - a1 - 4iroe0£sin (ю4) + 4iroe0£sin (ю(-l))
+--—--- sin (yz)
2sin (yl) U '
с неизвестными коэффициентами ai и a^ .
Замечание. Отметим, что полученное уравнение совпадает с уравнением Галлена, полученным в [18] из физических соображений.
Так как в работе рассматривается симметричный вибратор, то окончательно имеем
U' 0')
е 4Ra (z z') ^/_4/'юе0е sin (coz) +
4nRa (z - z') y2 _
c
+ 1 ■ 2. cos (yz )+ 2 ■ 1. sin (yz). (17)
2cos(Yl) V ' 2sin(Yl) V '
Вычислительную схему решения уравнения (17) построим, пользуясь технологией метода сплайн-коллокации со сплайнами нулевого порядка.
lk
Введем узлы Xk = -l +—, k = 0,1,...,2N.
Приближенное решение уравнения (17) будем искать в виде функции
2 N -1
IzN (z )= X akVk (z ). (18)
k=0
где
Vk к >={0; [, k,k=0-1....2 n -1-
Ak = [xk,xk+1), k = 0,1,...,2N - 2,Д2N-1 = [x2N-1,x2N].
Так как на торцах вибратора Iz (±l) = 0, то полагаем a0 = a2n-1 = 0. Коэффициенты {ak} , k = 1,...,2N - 2, ищутся из следующих систем уравнений.
Первый вариант. Вычисляя интеграл в левой части уравнения (17) по квадратурной формуле прямоугольников, имеем
l 2^2 e~iYRa(kj ^esin (rax-) a1 - a2 , ,
— X ak-=---^—— + ——T^cos( yx,) +
N k=1 k 4nRa(k,j) y2-ю2 2cos(yI) V]>
+ a2 -(a1;)sin(yx- ), j = 0,1,...,2N-1. (19)
2sin(Yl) v '
Здесь Ra (k, j) = ^(( _ xj )
2 2
+ a , a - радиус стержня антенны.
Таким образом, получаем систему из 2N уравнений с 2N неизвестными.
Эту систему более удобно представить в следующем виде, вынеся в левую часть неизвестные:
, 2N-2 -iyRa (k,j)
l vn e a a1 - a2 / \ — > ak---i——cos[yx,- )-
N k=i k 4nRa(k,j) 2cos(yl)
-sin(yx- ^), j = 0,1,...,2N-1. (20)
2sin(yl) у2-ю2
Система (20) имеет единственное решение, если существует такой вектор v = (vo, •••, V2n-i), vk =±1, k = 0,1, •••,2N -1, что матрица, стоящая в левой части системы
, 2N-2 -iyR„ (k,j) v ( l V a e a ai - a2
j(N ak 4nRa(k, j) 2cos(12/)C0S(Yrj)
, 47w£q£ sin (oar.-)
2 • ( ;^in (j )) = VJ-2 2 , j = 0,1,...,2N -1,
2sin(yl) v 7 y2 -w2
a2 - ai 2 1 -sin
имеет в некотором банаховом пространстве отрицательную логарифмическую норму.
Это следует из теоремы Лозинского [19].
Погрешность решения уравнения (20) зависит от класса функций, к которому принадлежат искомые токи. Если они принадлежат классу функций, имеющих ограниченные производные первого порядка, то можно показать
(по аналогии с работой [20]), что погрешность имеет порядок где N -
число узлов коллокации.
Второй вариант. Подставляя функцию (18) в уравнение (17), имеем
2^2 4/юе0е(у) а _а2 , V
£ ак = у2 _М2 + 2^*)«»^ ) +
a2 - ai + 2 , \ sin 2sin (y/ )
in (yx- ), j = 0,1,...,2N -1, (21)
здесь
Xk+1 e~iYRa(z'-x;) r-
Pjk = J -(-) , Ra (( - xj )) ( - xj ) + a , j = Q,1,...,2N -1.
xk 4nRa (z - Xj )
Как и в первом варианте, систему уравнений удобно представить в виде
2 N-2 - -
У pjk a k ———г^т cos (yx. )—°2 , sin (yx. ) = k=1 j k 2cos(y/) 1 j 2sin(Yl) 1 j
4/w£o£ sin ((OXj )
Y2-(
j = 0,1, ...,2N -1.
(22)
При реализации вычислительных схем (20), (22) применяется непрерывный метод решения операторных уравнений. Так как он одинаково реализуется для обоих уравнений, ограничимся рассмотрением уравнения (22). Поставим в соответствие уравнению (22) задачу Коши
dai (а)
= vo
d а
ai(a) - а2(а)
2cos(y/)
2 N- 2 4/юе0е sin ((ox0)
X P0k ak(а)--2-2--
k=1 Y
cos (yxo )- а2(а) ,ai(a) sin (yxO )
2sin(y/)
(23)
d а у (а) d а
= V;
2 N-2 4(£o£ sin ((xj) а1(а) - а2(а)
X pjkak(а) Y^2 2cos(y/)
, j = 1,...,2N - 2;
k=1 Y2
cos
((xj)-
а2(а) - a1(a)sin
2sin(Yl)
(Yxj)
da2 (а)
: V2 N-1
2N 2 4/юе0еsin(wx2N-1)
d а
а1(а) - а2(а)
X P2N-1,k ak (а)-k=1
Y
(
cos(YX2N-1)- а2(а) /а1;(а) sin(Yx2N-1)
2cos(yI ) 2sin(Yl)
a1 (0) = 0, a2 (0) = 0, a j (0) = 0, j = 1,...,2 N - 2.
(24)
Константы = ±1, ] = 0,1,...,2 N — 1, подбираются таким образом, чтобы логарифмическая норма матрицы, стоящей в левой части системы уравнений (22), была отрицательной в метрике соответствующего банахова пространства.
Тогда согласно теореме 1 решение задачи Коши (23), (24) сходится в метрике выбранного банахова пространства к решению системы уравнений (22).
Замечание. В качестве начальных условий в задаче Коши (23), (24) может быть взят любой набор констант.
Замечание. Для решения задачи Коши (23), (24) может быть использован любой численный метод решения обыкновенных дифференциальных уравнений. В приведенных ниже примерах задача Коши (23), (24) решалась методом Эйлера.
Вторая вычислительная схема
Приведем еще одну вычислительную схему для решения уравнения Поклингтона.
Будем исходить из уравнения
± fЧ )
- + k2
e-iYRa (Z-Z )
Ra (Z - z')
—dz' = iweEZ (z),
где
a (Z - z') = Va2 + (z - z')2 .
я
I
Проведя несложные вычисления, приходим к уравнению 1 Г1
где
(
h(z,z )=-
( z - z' )
Ra (Z - z')
k2 +-
(z - z' )2
a2 + (z - z ')2
3/2
1
a 2 + (z - z')"
+
3(z - z' )2
(a2 + (z - z' )2 )
+
+ik
(z - z')2
2 (z - z' )2
a2 + (z - zf (a2 + (z - z ')2 ) + (a2 + (z - z')2 )
(25)
Вычислительные схемы для решения уравнения (25) строятся по аналогии с вычислительными схемами (20) и (22). Здесь они не приводятся.
4. Приближенное решение уравнения Галлена
При выводе уравнения Галлена тонкого вибратора полагают, что а ^ I, а ^ X. Здесь а - радиус поперечного сечения вибратора, 21 - его длина, X - длина волны.
При этом делаются следующие предположения [18]:
1. Поверхностные электрические токи 1г (z) вместе с магнитными эквивалентными токами (z) заменяются бесконечно тонкой нитью продольного электрического тока ^^) = 2па1-, и полагают I- (-1) = I- (I). Продольный электрический ток в области зазора считается непрерывным. Торцевые токи игнорируются.
2. Касательная составляющая вектора напряженности электрического поля Ez ^), создаваемая током ^ (z) на боковой поверхности цилиндра радиуса а , обращается в нуль всюду, кроме окрестности зазора шириной Ь :
0.1 Щ >
Е') = ^ \ , Ь
гшЕ (г), Щ < —,
где Ь - ширина зазора.
Замечание. Для узких зазоров функция Ег (г),—2 < г < 2, обычно
полагается постоянной.
С учетом этих предположений получено (подробности см. в [18]) уравнение Галлена:
J-2Iz (/)K(z -z')dz' =
2ra'V
= C1cos kz + C2 sin kz--—- sin k | z |, - /1 < z < 12, (26)
где
—ik'J ( z-z ')2+a2
K (z-z)=e
V( z - z')
2 2 2 + a
Правая часть уравнения (26) может рассматриваться, как суперпозиция
2кУ
трех бегущих волн, одна из которых имеет амплитуду -, где V - ЭДС
W
генератора, W = k / (юе).
В случае симметричного вибратора имеем
J* lz (z')K(z - z)dz = C3cos kz - sin k | z |, - / < z < /. (27)
Замечание. Уравнение Галлена может быть получено из уравнения Поклингтона тождественными преобразованиями. Это продемонстрировано выше при выводе вычислительной схемы для уравнения Поклингтона.
Построение численных методов для уравнения Галлена подобно аналогичному построению для уравнения Поклингтона, поэтому опускается.
5. Вычислительный эксперимент
В этом разделе приведены результаты численных экспериментов решения уравнения Галлена с 5 -функцией в правой части. Таким образом, ставится задача решения непрерывным методом операторных уравнений следующего уравнения:
l -iYRa (z - z')
J!z (zl-hr,--^z' = 5(z). (28)
- 4nRa (z - z)
При этом 5 -функция аппроксимируется 5 -образными функциями. Physical and mathematical sciences. Mathematics 139
Пример 1.
l
f h (Z)
-iyRa (z - z')
-dz' = f (z),
-l
4nRa (z - z')
l = 10, a = 0,5, y = 0,0001366732941, 61 = 1,52 = 0,8, h = 10, f0, -1 < z <-51,
-51 < z <-62, -62 < z 62 < z <61,
(29)
f (z ) =
h (61 + z)
61 -62
h,
h (61 - z)
61 -62
0, 61 < z < I.
График функции, стоящей в правой части уравнения (29) приведен на
рис. 1.
Рис. 1. Правая часть уравнения (29)
Уравнение (29) решалось при числе узлов коллокации N = 50. Количество итераций - 5000.
Результаты вычислений представлены на рис. 2 и 3. Пример 2.
г e~iYR
f Iz (Z' )e- dZ = f (Z ),
- zV ' 4kR w
(30)
l = 10, a = 0,5, Y = 0,0001366732941, 61 = 0,2, 62 = 0,1,
f (z ) =
0, - / < z <-51,
h (51 + z)
51 -5 2
h (51 - z)
-5i < z <-52
■ 52 < z < 52
52 < z <5
51 -52 0, 51 < z < /.
Рис. 2. Действительная часть волны возбужденного тока
Рис. 3. Мнимая часть волны возбужденного тока
Уравнение (30) решалось при числе узлов коллокации N = 50. Количество итераций - 5000.
Результаты вычислений представлены на рис. 4 и 5.
Рис. 4. Действительная часть волны возбужденного тока
Рис. 5. Мнимая часть волны возбужденного тока
На рис. 6, 7 представлена зависимость точности решения уравнений Галлена от числа итераций при реализации непрерывного операторного метода вычислительной схемой Эйлера.
0.7
0.6
1000 2000 3000 4000 J000 6000 7000 8000 9000 10000 Число итераций
Рис. 6. Погрешность при числе узлов коллокации N = 50
1.0
0.9 0.8 $ 0.7
1000 2000 3000 4000 i000 6000 7000 8000 9000 10000 Числ о ит ераций
Рис. 7. Погрешность при числе узлов коллокации N = 100 Библиографический список
1. Сазонов, Д. М. Антенны и устройства СВЧ / Д. М. Сазонов. - Москва : Высшая школа, 1988. - 434 с.
2. Неганов, В. А. Метод сведения уравнения Поклингтона для электрического вибратора к сингулярному интегральному уравнению / В. А. Неганов, И. В. Матвеев, С. В. Медведев // Письма в журнал теоретической физики. - 2000. - № 12. -С. 86-94.
3. Неганов, В. А. Сингулярные интегральные уравнения как метод физической регуляризации некорректных электродинамических задач радиотехники и связи /
B. А. Неганов // Успехи современной радиотехники. - 2005. - № 12. - С. 16-24.
4. Бойков, И. В. Применение гиперсингулярных интегральных уравнений к численному моделированию электрического вибратора / И. В. Бойков, Д. В. Тарасов // Известия высших учебных заведений. Поволжский регион. Технические науки. - 2008. - № 8. - С. 94-106.
5. Вычислительные методы в электродинамике / под ред. Р. Митры. - Москва : Мир. -1977. - 488 с.
6. Бойков И. В. Об одном непрерывном методе решения нелинейных операторных уравнений / И. В. Бойков // Дифференциальные уравнения. - 2012. - Т. 48, № 9. - С. 1308-1314.
7. Далецкий, Ю. Л. Устойчивость решений дифференциальных уравнений в банаховом пространстве / Ю. Л. Далецкий, М. Г. Крейн. - Москва : Наука, 1970. -534 с.
8. Деккер, К. Устойчивость методов Рунге-Кутты для жестких нелинейных дифференциальных уравнений / К. Деккер, Я. Вервер. - Москва : Мир, 1988. - 334 с.
9. Burke, G. J. Method of Moment / G. J. Burke, A. J. Poggio // Numerical Electromagnetics Code (NEC). - Lawrence Livermore Laboratory, 1981.
10. Barrera-Figueroa, V. Simplification of Poklington's integral equation for arbitrary bent thin wires / V. Barrera-Figueroa, J. Sosa-Pedroza and J. Lopez-Bonilla // WIT Transactions on Modelling and Simulation. - 2005. - Vol. 39. - P. 563-574.
11. Pocklington, H. C. Electrical Oscillations in Wires / H. C. Pocklington. - London, England, Cambridge Phil. Soc. Proc., 1897. - Р. 324-332.
12. Rawle, W. D. The Method of Moments: A Numerical Technique for Wire Antenna Desing / W. D. Rawle // Higt Fraquency Electronics. - 2006. - Vol. 2. - P. 43-47.
13. Вайникко, Г. М. Численные методы в гиперсингулярных интегральных уравнениях и их приложения / Г. М. Вайникко, И. К. Лифанов, Л. Н. Полтавский. -Москва : Янус-К, 2001. - 508 с.
14. Бойков, И. В. Об одном приближенном методе решения линейных гиперсингулярных интегральных уравнений на разомкнутых контурах интегрирования / И. В. Бойков, А. И. Бойкова // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2016. - № 2. - С. 68-84.
15. Бойков, И. В. Приближенные методы решения гиперсингулярных интегральных уравнений первого рода с особенностями второго порядка на классах функ-
2 —1/2
ций с весами (1 -1 ) / И. В. Бойков, А. И. Бойкова // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2017. -№ 2. - С. 79-90.
16. Бойков, И. В. Приближенное решение гиперсингулярных интегральных уравнений первого рода с особенностями второго порядка на классе функций с весом
±1 /2
((1 -1 )/(1 +1 )) / И. В. Бойков, А. И. Бойкова // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. 2019. - № 3. -
C. 76-92.
17. Бойков, И. В. Аналитические и численные методы решения гиперсингулярных интегральных уравнений / И. В. Бойков // Динамические системы. - 2019. -Т. 9 (37), № 3. - С. 244-272.
18. Войтович, Н. И. УКВ вибраторные антенны / Н. И. Войтович, А. В. Ершов, А. Н. Соколов. - Челябинск : Изд-во ЮУрГУ, 2002. - 85 с.
19. Лозинский, С. М. Оценка погрешности численного интегрирования обыкновенных дифференциальных уравнений. I / С. М. Лозинский // Известия вузов. Математика. - 1958. - № 5. - С. 52-90.
20. Boykov, I. V. New iterative method for solving linear and nonlinear hypersingular integral equations / I. V. Boykov, V. A. Roudnev, A. I. Boykova, O. A. Baulina // Applied Numerical Mathematics. - 2018. - Vol. 127. - P. 280-305.
References
1. Sazonov D. M. Antenny i ustroystva SVCh [Microwave antennas and devices]. Moscow: Vysshaya shkola, 1988, 434 p. [In Russian]
2. Neganov V. A., Matveev I. V., Medvedev S. V. Pis'ma v zhurnal teoreticheskoy fiziki [Letters to the Journal of Theoretical Physics]. 2000, no. 12, pp. 86-94. [In Russian]
3. Neganov V. A. Uspekhi sovremennoy radiotekhniki [Advances in modern radio engineering]. 2005, no. 12, pp. 16-24. [In Russian]
4. Boykov I. V., Tarasov D. V. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Tekhnicheskie nauki [University proceedings. Volga region. Engineering sciences]. 2008, no. 8, pp. 94-106. [In Russian]
5. Vychislitel'nye metody v elektrodinamike [Computational methods in electrodynamics]. Ed. by R. Mitry. Moscow: Mir, 1977, 488 p. [In Russian]
6. Boykov I. V. Differentsial'nye uravneniya [Differential Equations]. 2012, vol. 48, no. 9, pp. 1308-1314. [In Russian]
7. Daletskiy Yu. L., Kreyn M. G. Ustoychivost' resheniy differentsial'nykh uravneniy v banakhovom prostranstve [Stability of solutions of differential equations in a Banach space]. Moscow: Nauka, 1970, 534 p. [In Russian]
8. Dekker K., Verver Ya. Ustoychivost' metodov Runge-Kutty dlya zhestkikh nelineynykh differentsial'nykh uravneniy [Stability of Runge-Kutta methods for stiff nonlinear differential equations]. Moscow: Mir, 1988, 334 p. [In Russian]
9. Burke G. J., Poggio A. J. Numerical Electro-magnetics Code (NEC). Lawrence Livermore Laboratory, 1981.
10. Barrera-Figueroa V., Sosa-Pedroza J. and Lopez-Bonilla J. WIT Transactions on Modelling and Simulation. 2005, vol. 39, pp. 563-574.
11. Pocklington H. C. Electrical Oscillations in Wires. London, England, Cambridge Phil. Soc. Proc., 1897, pp. 324-332.
12. Rawle W. D. Higt Fraquency Electronics. 2006, vol. 2, pp. 43-47.
13. Vaynikko G. M., Lifanov I. K., Poltavskiy L. N. Chislennye metody v gipersingulyarnykh integral'nykh uravneniyakh i ikh prilozheniya [Numerical methods in hypersingular integral equations and their applications]. Moscow: Yanus-K, 2001, 508 p. [In Russian]
14. Boykov I. V., Boykova A. I. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2016, no. 2, pp. 68-84. [In Russian]
15. Boykov I. V., Boykova A. I. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2017, no. 2, pp. 79-90. [In Russian]
16. Boykov I. V., Boykova A. I. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2019, no. 3, pp. 76-92. [In Russian]
17. Boykov I. V. Dinamicheskie sistemy [Dynamical systems]. 2019, vol. 9 (37), no. 3, pp. 244-272. [In Russian]
18. Voytovich N. I., Ershov A. V., Sokolov A. N. UKVvibratornye antenny [VHF vibrator antennas]. Chelyabinsk: Izd-vo YuUrGU, 2002, 85 p. [In Russian]
19. Lozinskiy S. M. Izvestiya vuzov. Matematika [University proceedings. Mathematics]. 1958, no. 5, pp. 52-90. [In Russian]
20. Boykov I. V., Roudnev V. A., Boykova A. I., Baulina O. A. Applied Numerical Mathematics. 2018, vol. 127, pp. 280-305.
Бойков Илья Владимирович
доктор физико-математических наук, профессор, заведующий кафедрой высшей и прикладной математики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
Е-таП: boikov@pnzgu.ru
Айкашев Павел Владимирович
аспирант, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
Е-таП: math@pnzgu.ru
Образец цитирования:
Бойков, И. В. Применение непрерывного операторного метода к решению уравнений Поклингтона и Галлена для тонких проволочных антенн / И. В. Бойков, П. В. Айкашев // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2020. - № 3 (55). -С. 127-146. - DOI 10.21685/2072-3040-2020-3-10.
Boykov Il'ya Vladimirovich Doctor of physical and mathematical sciences, professor, head of the sub-department of higher and applied mathematics, Penza State University (40 Krasnaya street, Penza, Russia)
Aykashev Pavel Vladimirovich Postgraduate student, Penza State University (40 Krasnaya street, Penza, Russia)