ISSN 0868-5886 НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2019, том 29, № 4, c. 96-109
МАТЕМАТИЧЕСКИЕ МЕТОДЫ И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ
УДК537.213, 537.612, 517.5
© А. С. Бердников, Н. К. Краснова, К. В. Соловьёв,
А. Г. Кузьмин, С. В. Масюкевич, Ю. А. Титов, |Ю. К. Голиков
2019
ГИПЕРГЕОМЕТРИЧЕСКИЙ БАЗИС ДЛЯ ТРЕХМЕРНЫХ ГАРМОНИЧЕСКИХ ФУНКЦИЙ, ОДНОРОДНЫХ ПО ЭЙЛЕРУ С НЕЦЕЛОЧИСЛЕННЫМИ СТЕПЕНЯМИ ОДНОРОДНОСТИ
Электрические и магнитные поля, однородные по Эйлеру, являются удобным инструментом при синтезе электронно- и ионно-оптических систем со специальными свойствами. Известно, что скалярные потенциалы подобных полей представляют собой трехмерные скалярные гармонические функции, однородные по Эйлеру, с заданной степенью однородности. Вопрос о параметризации трехмерных однородных гармонических функций с целочисленными степенями однородности исчерпывающим образом решается с помощью формул Донкина для однородных гармонических функций со степенями однородности 0 и -1, теоремы о дифференцировании трехмерных однородных гармонических функций и формулы Томсона для однородных гармонических функций. Однако число аналитических формул, которые можно использовать для описания трехмерных скалярных гармонических функций при нецелочисленных степенях однородности, к сожалению, в настоящий момент не слишком велико, а вопрос об исчерпывающем описании таких функций весьма далек от своего окончательного решения. В то же время использование трехмерных однородных гармонических потенциалов с нецелочисленными степенями однородности значительно расширяет инструментарий разработчиков электронно- и ионно-оптических систем. Целью данной работы является построение гипергеометрического базиса, составленного из базовых однородных гармонических функций с нецелочисленными степенями однородности, с помощью которого любую трехмерную однородную гармоническую функцию, не имеющую сингулярных точек за исключением луча х = у = 0, г < 0, можно представить в виде бесконечного ряда наподобие ряда Фурье. Полученный результат, по-видимому, отчасти решает проблему исчерпывающего описания трехмерных скалярных однородных гармонических функций при нецелочисленных степенях однородности.
Кл. сл.: электрические поля; гармонические функции; функции, однородные по Эйлеру; принцип подобия траекторий в оптике заряженных частиц; формула Донкина; аналитические решения уравнения Лапласа
ПОСТАНОВКА ЗАДАЧИ
Введение в проблему*
Полезным инструментом при синтезе электронно- и ионно-оптических систем специального вида являются электрические и магнитные поля, однородные по Эйлеру [2-7]. Для полей, однородных по Эйлеру, напряженность электрического поля Е и/или индукция магнитного поля В должны удовлетворять не только уравнениям Максвелла для электромагнитного поля, но и тождеству однородности:
VA> 0: E (Ax,Ay,Az) = Am_1E (x, y, z), VA> 0: B (Ax,Ay,Az) = Am_1B (x, y, z),
(1)
Этот подраздел вставлен авторами статьи для полноты изложения, хотя он почти полностью воспроизводит раздел "Базовые определения" статьи [1] этого же авторского коллектива. (Ред.)
где т — степень однородности поля (не обязательно целочисленная). Траектории движения заряженных частиц в электростатических и магни-тостатических полях, однородных по Эйлеру, подчиняются принципу подобия траекторий Ю.К. Голикова [8, 9]. Отсюда следуют уникальные оптические свойства устройств, управляющих движением заряженных частиц, если эти устройства используют однородные по Эйлеру электрические и магнитные поля [2-6, 10-18].
Однородные по Эйлеру электрические или магнитные поля, подчиняющиеся тождеству (1), характеризуются скалярным электрическим или магнитным потенциалом в виде гармонической функции, которая будет однородной (точнее, положительно однородной) по Эйлеру в смысле, который придается этому термину в классическом математическом анализе [19, 20]:
VA > 0: U (Ax,Ay,Az) = AmU (x,y, z):
(2)
где т — степень однородности функции, совпадающая со степенью однородности электрического или магнитного поля, как она определена с помощью тождеств (1). Единственным исключением из данного общего правила являются поля с нулевой степенью однородности — для таких полей скалярный потенциал и в самом общем случае имеет вид:
и (х, у, z) = и0 (х, у, z) + А 1п (г + д/х2+у^+22), (3)
где ио является однородной гармонической функцией нулевой степени, а А — произвольная константа. Очевидно, что при А ф 0 выражение (3) больше не является однородной по Эйлеру функцией, хотя градиент функции (3) и подчиняется условиям однородности (1) для электрического либо магнитного поля. Вопрос об однородности скалярных и/или векторных потенциалов для полей, однородных по Эйлеру, подробно исследуется в [21].
Исследования Ю.К. Голикова [4, 22, 23] наглядно демонстрируют, каким удобным и эффективным инструментом при синтезе новых электронно- и ионно-оптических систем являются аналитические методы и, в частности, аналитические выражения для потенциалов электрических и магнитных полей. В случае однородных электрических и магнитных полей при получении аналитических выражений для их скалярных потенциалов требуется предпринимать дополнительные и целенаправленные усилия, поскольку таких функций в некотором смысле "много меньше", чем обычных гармонических функций.
Действительно, трехмерная гармоническая функция общего вида и (х, у, г) полностью определяется двумя произвольными функциями двух переменных — например, значением и(0)(х, у) = = и(х,у,0) функции и(х,у,г) вдоль плоскости г = 0 и значением и(0,п) ( х, у ) = ди ( х, у,0 )/дг нормальной производной функции и (х, у, г )
вдоль плоскости г = 0 . Это следует из того факта, что решение такой задачи Коши для трехмерного уравнения Лапласа однозначным образом находится (по крайней мере, в окрестности плоскости г = 0) с помощью ряда Шерцера по переменной г:
и ( х, у, г) = и(0) ( х, у ) + ги(0п > (х, у ) -
(и«0)(х,у) + иу; (х.
/(0)
(^ у ))-
3!
г4 + —
4!
(ЦТ)(х, у) + иуу,п)( х, у)) +
(и(0) (х,у) + 2и(0) (х, у) + и(0) (х,у)) +
у хххх V ' ) ххуу V ) уууу V
+ ^ (и(0'п) ( х, у) + 2и(0'п) ( х, у) + и(0'п) ( х, у ))
-! у хххх V ) ххуу V ) уууу V
(4)
где коэффициенты перед степенями переменной г выражаются однозначным образом через функции и(0), и^ и их частные производные, поскольку функция и обязана удовлетворять уравнению Лапласа.
Для получения трехмерной однородной гармонической функции степени т необходимо и достаточно, чтобы функции и (0)( х, у ) = и (х, у,0 ) и и(0,п)(х, у) = ди (х, у,0)/дг были однородными по Эйлеру со степенями однородности т и (т -1) соответственно (достаточность следует из формулы (4) и того факта, что частные производные однородных функций и (0)( х, у ) и и (0,п)( х, у ) сами являются однородными функциями соответствующей степени [19, 20]). Однако однородная по Эйлеру функция двух переменных задается с помощью произвольной функции одного переменного, как это следует из ее представления в виде /(х,у) = хкg(у/х) [19, 20]. Тем самым трехмерные однородные функции однозначным образом определяются через две произвольные функции всего лишь одного вещественного переменного в отличие от трехмерных гармонических функций общего вида, для однозначного определения которых в соответствии с формулой (4) требуется задать две произвольные функции двух переменных.
Текущее состояние
Для параметризации однородных гармонических функций 0-й степени и степени -1 используются явные формулы Донкина [2, 3, 12, 13, 2428]:
V (х,у,г) = НГ-^-1
^ г + г г + г )
У_1 ( х, у, г ) = - Н
У
г + г г + г
(5)
(6)
(здесь и далее для упрощения формул используется подстановка г = д/х2 + у2 + г2 ). Эти формулы устанавливают взаимно-однозначное соответствие
3
2
х
2
между решениями Н (р, q) двумерного уравнения Лапласа Нрр + = 0 и трехмерными однородными гармоническими функциями со степенями однородности 0 и -1 (где для упрощения внешнего вида формул нижние индексы обозначают частные производные по соответствующим переменным). Переход от формулы (5) к формуле (6) и обратно может быть выполнен с помощью формулы Томсона для однородных гармонических функций [29-42]:
Кт-1 (х, у, г) = г-2т-1 Vm (х, у, г) , (7)
где Vm (х, у, г) и V-m-1 (х, у, г) — это трехмерные
гармонические функции, однородные по Эйлеру с соответствующими степенями однородности.
Задача вычисления однородных гармонических функций с целочисленными степенями однородности может считаться полностью решенной. С помощью дифференциальных операторов Том-сона—Донкина первого порядка [39-42] либо с помощью теоремы о дифференцировании однородных гармонических функций [24, 27, 57] в комбинации с формулой Томсона для трехмерных гармонических функций [29-42] имеется возможность получать из алгебраических формул Донкина (5), (6) алгебраически-дифференциальные формулы для произвольной целочисленной степени однородности. Эти формулы позволяют выразить трехмерные однородные гармонические функции общего вида через решения двумерного уравнения Лапласа [24-27], причем гарантируется, что ни одна трехмерная однородная функция не будет пропущена. Публикации [25, 26] являются, по всей видимости, первыми, в которых исчерпывающим образом решен вопрос о перечислении всех трехмерных однородных гармонических функций с целочисленными порядками однородности. Некоторым недостатком указанной процедуры является то, что для целочисленных степеней однородности, отличных от 0 и -1 , несовпадающие друг с другом решения двумерного уравнения Лапласа могут порождать одинаковые трехмерные однородные функции. При этом чем больше модуль степени однородности, тем выше коэффициент повторяемости для трехмерных однородных функций, генерируемых с помощью соответствующих алгебраически-дифференциальных формул.
К сожалению, решение задачи полного и конструктивного перечисления трехмерных однородных функций с нецелочисленными степенями однородности на настоящий момент неизвестно. Некоторые полезные аналитические выражения для гармонических функций, однородных по Эйлеру, которые можно использовать в качестве скаляр-
ных потенциалов, приводятся в [24, 43-49]. Ограниченную практическую применимость имеют интегральные выражения для однородных гармонических функций с нецелочисленными степенями однородности, которые приводятся в [45, 48, 49]. Однако, если ограничиваться лишь этими частичными аналитическими выражениями, то при оптимизации новых электронно- и ионно-оптических систем есть большой риск пропустить действительно оптимальное решение.
Формулировка задачи
Данная работа посвящена исследованию одного из возможных подходов к решению задачи о полном перечислении однородных гармонических функций с нецелочисленными степенями однородности. А именно, для разложения в ряд общего вида трехмерных однородных гармонических функций с произвольной степенью однородности будут сконструированы базисы из эталонных однородных гармонических функций, являющихся аналогом тригонометрического ряда Фурье. Как известно, с помощью ряда Фурье любую непрерывную периодическую функцию (возможно, обладающую конечным числом точек разрыва первого рода) можно однозначным образом разложить в ряд по тригонометрическим функциям. Аналогичным образом любую трехмерную однородную гармоническую функцию общего вида, у которой нет дополнительных сингулярных точек кроме, возможно, начала координат и/или лучей x = y = 0, z<0 либо x = y = 0, z> 0, можно разложить в сходящийся ряд по рассматриваемым базисным функциям. Хотя такой подход и нельзя рассматривать как окончательное решение проблемы полного и конструктивного перечисления трехмерных однородных гармонических функций (поскольку при разложении в ряд достаточно простые алгебраические выражения способны превращаться в трудно обозримый бесконечный ряд с весьма сложной структурой), подобного рода базисы могут оказаться полезным инструментом при исследовании трехмерных однородных гармонических функций.
СВЯЗЬ МЕЖДУ ТРЕХМЕРНЫМИ ОДНОРОДНЫМИ ГАРМОНИЧЕСКИМИ
ФУНКЦИЯМИ И ДВУМЕРНЫМИ ЭЛЛИПТИЧЕСКИМИ УРАВНЕНИЯМИ
Для дальнейших выкладок потребуется взаимно-однозначный переход от трехмерных однородных гармонических функций к двумерным функциям, которые удовлетворяют некоторым вспомогательным двумерным эллиптическим уравнениям. Данная операция представляет определенный
интерес и сама по себе, поэтому рассмотрим этот момент подробнее.
Для трехмерных гармонических функций степени 0 и степени -1 имеются формулы Донкина (5), (6), которые устанавливают взаимнооднозначное соответствие между решениями Н (р, q) двумерного уравнения Лапласа Нрр + Н = 0 и трехмерными однородными гармоническими функциями со степенями однородности 0 и -1. По аналогии любую однородную функцию и (х, у, z) степени т с помощью взаимно-однозначной замены переменных Донкина
х
р= q =
z + г Л = г ,
z + г У
х=
2 рЛ
1 + р2 + q2'
У =
2qЛ
z =
1 + р2 + q2'
Л(1 - р2 - q2)
1 + р2 + q 2
(8)
можно представить в виде
и ( х, у, z ) = гтН 1
^ z + г z + г )
(9)
(это представление является слегка модифицированной формой универсального представления f (х1,х2,...хп) = х1mh(х2/х^...хп/х1) для однородных функций степени т [20, 21]). С помощью подстановки выражения (9) в трехмерное уравнение Лапласа ихх + иуу + и^ = 0 легко можно установить, что для того, чтобы функция (9) была гармонической (удовлетворяла трехмерному уравнению Лапласа), необходимо и достаточно, чтобы функция Н (р, q) удовлетворяла двумерному эллиптическому дифференциальному уравнению в частных производных определенного вида:
д2 Н ( р, q ) д2 Н (р, q )
+-
др дq
4т ( т +1)
(1+р2 + q2 )2
Н ( р, q ) = 0.
(10)
Подстановка (9), (10) является взаимнооднозначной. Действительно:
а) если и (х, у, z) является однородной гармонической функцией степени т , то существует такая функция Н (р, q) , с помощью которой функция и может быть представлена в виде (9), при-
чем функция Н (р, q) обязана будет подчиняться уравнению (10);
б) если функция Н (р, q) подчиняется уравнению (10), то функция и (х, у, z) , вычисленная в соответствии с правилом (9), будет однородной гармонической функцией степени т ;
в) разным функциям и (х, у, z) соответствуют разные функции Н (р, q), а при выборе разных функций Н (р, q) получаются разные функции и (х, у, z) .
Например, если требуется проанализировать некоторое соотношение между однородной гармонической функцией и (х, у, z) степени т и однородной гармонической функцией V (х, у, z) степени ], то подстановка
х у
и ( х, у, z ) = гтН ,1 ^ z + г z + г )
V ( х, у, z ) = г- ,1 ^ z + г z + г )
позволяет без потери общности свести задачу к анализу эквивалентного соотношения между функциями двух переменных Н (р, q) и — (р, q), подчиняющимися уравнениям
<2Н (p, q) + д2Н ( А q) + 4т (т +Н ^ q) = 0,
др2
V (1 + р2 + q2 )2
д2-(p,q) +д—М + 4](]+1) -() 0
др2 дq2 (1 + р2 + ч 2 )2 ^ •
Подстановка (9) не является единственно возможной. В равной степени можно использовать подстановки
и (х, у, * )=( ,+г ун \, )
и (х, у, г ) = гтН I-, у|,
Г г
и ( х, у, z ) = zmH\-, у|,
z z
(11) (12) (13)
которые, как и формула (9), являются взаимнооднозначным способом параметризации трехмерных однородных функций общего вида с помощью
произвольных вспомогательных двумерных функций Н (р, q) . Для того, чтобы функция вида (11) была однородной гармонической функцией, двумерная функция Н (р, q) должна удовлетворять эллиптическому уравнению
Л 2 2\ Гд2 Н д2 Н^ (1 + Р + q ^ + дд^
- 4mq — + 4т2 Н = 0. дq
, дН
- 4тр--
др
(14)
Для того чтобы функция вида (12) была однородной гармонической функцией, двумерная функция Н (р, q) должна удовлетворять эллиптическому уравнению
(1 - Р2 - 2 pq
д 2Н / 2ч д 2Н
2 -гч —+(1 - д )—г др дрдд дд
„ дН _ дН / - 2 р--2q-+ т ( т +1) Н = 0.
др дq
(15)
Наконец, для того чтобы функция вида (13) была однородной гармонической функцией, двумерная функция Н (р, q) должна удовлетворять эллиптическому уравнению
д2 н д2 Н (1 + р2 )^гг + 2 рд— + (1 +
др
дpдq
('+ч1)
д2 Н
дч2
ности с центром в начале координат. При этом полупространство г < 0 отображается на внешнюю часть р2 + д2 > 1 этой же окружности, а плоскость г = 0 отображается на саму окружность р2 + д2 = 1.
Перейдем в уравнении (10) от декартовых переменных (р,д) к полярным координатам (и,р) :
I и=4 р2+ч2,
[р = а.гс^ ( д/ р )
р = и cos р, д = и sin р.
(17)
После такой замены переменных решения Нт (р, Ч) уравнения (10), полученного для трехмерных однородных гармонических функций степени т, преобразуются в решения Фт (и, р) эквивалентного линейного дифференциального уравнения второго порядка в частных производных:
д 2Ф„
1 д 2Ф т 1 дФ т 4т (т +1)
т + т + т + \
ди2 и2 др2 и др (1 + и2 )2
-Ф т = 0. (18)
дН дН - 2 (т-1) р---2( т-1) д— + т (т-1) Н = 0. (16)
Приведенные здесь способы параметризации однородных гармонических функций являются, по-видимому, самыми очевидными, но далеко не единственными. Выбор того, какую именно подстановку следует использовать для решения конкретной задачи, определяется удобством манипулирования тем или иным двумерным эллиптическим уравнением в частных производных применительно к рассматриваемой задаче, а также эстетическими предпочтениями исследователя.
КОНСТРУИРОВАНИЕ ГИПЕРГЕОМЕТРИЧЕСКОГО БАЗИСА ДЛЯ КРУГА
Рассмотрим параметризацию трехмерных однородных гармонических функций в форме (9), (10). С помощью замены переменных (8) трехмерное полупространство г > 0 отображается на внутреннюю часть р2 + д2 < 1 единичной окруж-
Предположим, что функция Нт (р, д) не имеет особых точек внутри единичного круга р2 + д2 < 1. Тогда функцию Ф: (и,р), которая является периодической по переменной р с периодом 2п , можно разложить в сходящийся ряд Фурье
Ф: (и,р) = СсФт,0 (и) +
+ ЕФт,^ (и)(Ск C0S кр + ''{к крр, (19)
к=1,да
где ск и sk — произвольные константы, а функции Фтк (и) в соответствии с условием (18) обязаны будут удовлетворять обыкновенному дифференциальному уравнению
ф: ,к (и)+-и Фт,к (и)+ и
+
4: (: +1)
Т1
: -
+ и2 )2
к
и
Ф:,к (и) = 0.
(20)
У уравнения (20) имеется единственное (с точностью до множителя) решение, не имеющее сингулярности в точке и = 0 . Это решение можно записать как
Ф.
(и) = ик (1 + и2):+1 х
х 2Е (1 + т, 1 + т + к;1 + к; - ц2), (21)
где 2^ (а,ь;с;z) — гипергеометрическая функция [50-52].
Гипергеометрическая функция, обозначаемая
как 2а,ь;с; z), где а,Ь,с — параметры, а z —
переменная, вообще говоря, комплексная, удовлетворяет гипергеометрическому уравнению Эйлера
z (1 - z ) z ) + [с-( а + Ь +1) z ] Е'(z )-
-аЬЕ ( z ) = 0. (22)
Если с Ф 0, -1, -2,..., то гипергеометрическая функция разлагается в сходящийся при < 1 гипергеометрический ряд
, Ц (a,»;*z)=I+± f+
v ' c I! c(c +1)
2!
- +... =
=i + У
П
(a + 5 )(b + 5 )
A j-I (I + 5 )(c + 5 )
r(a + j)r(b + j) zj
j =i,m
Г( c) у _
Г( a )Г( b ) ^ Г( c + j )
j!
а при условии сходимости рассматриваемого ниже интеграла (т.е. при Re(c) > Re(Ь) > 0) может быть представлена как
, a, b; c; z ) =
Г( c )
Г(Ь )Г( c - b )
х]У-1 (1 -т)с-ь-1 (1 -тz )-а dт. (23)
0
Формула (23) обеспечивает при Re(c) > > Re(b) > 0 корректное аналитическое продолжение гипергеометрического ряда на всю комплексную плоскость с разрезом на действительной оси от 1 до и, в частности, позволяет вычислять такую гипергеометрическую функцию конструктивным образом при любых вещественных значениях z, удовлетворяющих условию -да < z < 1.
Для гипергеометрических функций с параметрами, используемыми в формуле (21), интеграл (23) расходится, если т < -1 либо т > 0 . Однако с учетом того, что с помощью линейных дифференциальных операторов Томсона—Донкина первого порядка [39-42] можно неограниченно понижать либо повышать степень однородности т в соответствии с правилами т ^ т +1 и т ^ т -1, а при т = 0 и т = -1 однородные гармонические функции исчерпывающим образом могут быть
описаны с помощью алгебраических формул Дон-кина (5) и (6), то возможность конструктивно вычислять функции (2I) лишь при -I < m < 0 представляется вполне достаточной. Для упрощения вычисления решений (2I) и/или для вычисления решений (2I) при степенях однородности m, отличных от интервала значений -I < m < 0 , обеспечивающего сходимость для интеграла (23), можно использовать разнообразные линейные и квадратичные гипергеометрические тождества [50-52], позволяющие достаточно гибко изменять аргументы гипергеометрических функций.
Функции Фm,0 (м) , Фт,к (м)C0sФщк (м) sin заданные с помощью формул (2I), образуют базис для решений Фт (м,ф) уравнения (I8), не имеющих сингулярных точек внутри круга p2 + q2 <I. После обратной замены переменных (I7) и обратной "донкиновской" подстановки (8) эти функции образуют базис из однородных гармонических функций.
Однородная гармоническая функция, порождаемая функцией Фm0 (м), соответствует трехмерному осесимметричному решению. Однородные гармонические функции, порождаемые функциями Фmk (м)cosкф , Фтк (м)sinкф, можно интерпретировать как трехмерные мультипольные гармоники, характеризующие однородные гармонические функции соответствующей степени.
С помощью полученного базиса любая однородная гармоническая функция степени m, не имеющая дополнительных особых точек на полупространстве z > 0, за исключением неизбежной особой точки в начале координат и составленного из особых точек луча x = y = 0, z < 0 (результат использования "донкиновских" аргументов в формуле (9)), может быть представлена в виде сходящегося ряда. Особыми точками считаются точки с сингулярным поведением функции, а также точки, в окрестности которых функция не может быть разложена в сходящийся ряд Тейлора, т.е. не является аналитической.
КОНСТРУИРОВАНИЕ ГИПЕРГЕОМЕТРИЧЕСКОГО БАЗИСА ДЛЯ КОЛЬЦА
У уравнения (20) существует еще одно решение ^тк (ц) , линейно независимое от приведенного ранее решения Фтк (ц), не имеющего сингулярности в точке ц = 0, которое демонстрирует сингулярное поведение при ц ^ 0 . С учетом того, что вронскиан
2
z
Ж (и) = Ф:,к (и)^ :,к (и) - Ф :,к (и)^:,к (и) =
:,к (и))
2 А
dи
Ф:,к (и) ^:,к (и)
составленный из решений линейного дифференциального уравнения (20), удовлетворяет уравнению
Ж'( и) + — Ж (и) = 0, и
это решение может быть представлено в соответствии с формулой Лиувилля—Остроградского в интегральной форме как
и
^:,к (и) = Ф:,к (и)|
1
ds
и0
Ф2:,к (S )
s) s
^ :,к (и) = О
и
| } {-к - :, - :} {-к,0} {}
(24)
где
о
• •,а„}, {а_,,-••,ар}
К-
{¿1,...,ь}, {ьы,..„ьч}
1 ПГ(Ь, - s)ПГ(1
к Г Л
- а] + s
)
2ж\
ПГ(1 -Ь + s)ПГ(аз -s)
1=п+1
-г' ds
1=1+1
это — О -функция Мейера ([51, п. 5.3]). (Здесь контур интегрирования L выбирается как разрез вдоль мнимой оси от -да до , возможно, с захватом некоторого участка действительной оси,
при котором все полюсы функций Г(Ь - '), где 1 < 1 < I, остаются справа от разреза, а все полюсы функций Г(1 - а^ + ') , где 1 < 1 < п , остаются слева от разреза.) При к ф 0 решение, линейно независимое ¥:к (и) с сингулярностью при и = 0, может быть также выражено с помощью гипергеометрической функции
^:,к (и) =
=и
(1 + и2 ):+1 (1 + :,1 + : - к;1 - к; -и2), (25)
конструктивное вычисление которой, впрочем,
сопряжено с определенными техническими трудностями (интегральная формула (23) для вычисления гипергеометрической функции, которая используется в формуле (25), не может быть использована, поскольку при таких параметрах в случае к > 1 интеграл (23) расходится).
Однако можно поступить проще. Непосредственной проверкой можно убедиться, что если функция Ф(и) удовлетворяет уравнению (20), то этому же уравнению обязана будет удовлетворять также и функция Ф(1/и). (Отсюда, в частности,
следует, что если функция Н (р, д) удовлетворяет уравнению (10), то и функция
Н (р/ (р2 + д2), д/ (р2 + д2)),
С помощью специальных функций парное решение ¥тк (и) выражается как
получаемая после замены переменных, имеющей смысл инверсии плоскости ( р, д) относительно
единичного круга с центром в начале координат, также удовлетворяет уравнению (10).) Поэтому в качестве дополняющего базиса, позволяющего разлагать решения уравнения (18) в кольце 0 < £ < и < 1, можно использовать функции
¥:,к (и) =
= и
1+ии
г
2 ^
1 + :,1 + т + к;1 + к; —1-и2
\
. (26)
В итоге функции
Ф:,0 (и) , Ф:,к (и) C0S кр , Ф:,к (и) sin кр , ¥:,0 (и) , ¥:,к (и) C0S кр , ¥:,к (и) sin кр ,
заданные с помощью формул (21) и (24), (25) либо (26), образуют базис для тех решений Ф: (и,р) уравнения (18), которые не имеют сингулярных точек внутри кольца 0 < £ < р2 + д2 < 1:
Ф: (и,р)= С0Ф:,0 (и) +
+
Е Ф :,к (и)( Ск C0S кр+ 'к Sin кр)
1 +
1 = к ,с
+ :,0 (и) +
+
Е ¥:,к (и)(dk C0S кр + ек Sin кр). (27)
]=к
После обратной замены переменных (17) и (8) эти функции образуют расширенный базис для однородных гармонических функций, а формула (27) переходит в разложение трехмерной однородной гармонической функции степени :. Для трехмерных однородных гармонических функций,
к
г
которые можно представить с помощью ряда, получаемого из разложения (27) при обратном переходе к декартовым координатам (х, у, z), кроме
особой точки в начале координат и состоящего из особых/сингулярных точек луча х = у = 0, z < 0, имеется также состоящий из сингулярных точек луч х = у = 0, z > 0, являющийся прообразом точки
р = q = 0 на плоскости (р, q) . Однако, никаких
других сингулярных точек у этих трехмерных однородных гармонических функций быть не может, поэтому как ряд (19), так и ряд (27) не покрывают все возможные трехмерные однородные гармонические функции степени т.
ВСПОМОГАТЕЛЬНЫЕ ТОЖДЕСТВА, ПОЛУЧАЕМЫЕ ДЛЯ ГИПЕРГЕОМЕТРИЧЕСКОГО БАЗИСА
В публикациях [39-42] приводятся дифференциальные аналоги первого порядка формулы Том-сона для трехмерных гармонических функций. Эти формулы преобразуют трехмерные гармонические функции в новые трехмерные гармонические функции, причем форма трансформирующих выражений сознательно выбрана таким образом, чтобы при подстановке в них однородных по Эйлеру функций на выходе получались новые однородные по Эйлеру функции, возможно, с другой степенью однородности. Указанный список формул выглядит следующим образом:
V(т)( х, у, г) = и(т>(х, у, г) ,
Кт) |
V<т-1)(х, у, г) = ихт)( х, у, г) , VК1)(х, у, г) = иут)( х, у, г) , Кт-1) /V ,, = тт(т)!
,= и{т)(
| = 77(т)!
V(т-1)(х, у, г) = иХт)( х, у, г) ,
(28)
(29)
(30)
(31)
V(т+1)( х, у, г ) = хи (т)( х, у, г ) +
(т) |
+ (х2 - у2 - г2 )ихт)(х,у, г) +
+ 2 хуиут)( х, у, г ) + 2хги[т)( х, у, г ) ,(32)
(т) /
V (т+1)( х, у, г )= уи (т)( х, у, г) +
+ 2хуихт) (х, у, г) + 1-х2 + у2 - г2 )и (т)( х, у, г)
+ (-х2 + у2 - г2)иут)(х,у,г) + 2уги|т)(х,у,г), (33)
Vт+1 (х, у, г) = ги(т) (х, у, г) + 2x2^ (х, у, г) + .у,т)( х, у, г )+(-х2 - у2 + г2 )и[т) (х, у, г
V(т)( х, у, г ) = = хи[т) (х,у, г) + уиут) (х,у, г) + ги^ (х,у,г), (35)
Vт) (х,у, г) = уи{Г] (х,у,г) - хи{т (х,у,г), (36)
¥{т) (х, у, г) = ги[т) (х, у, г)- хи[т) (х, у, г) , (37)
^т) (х,у,г) = ги{т (х,у,г) - уи[т (х,у,г), (38)
V( т-1) ^^ г^-¿гит) (x,y, г),
V т) (х, у, г) = и(:] (х, у, г)
(т),
V -т) (х, у, г ) = ^ иут) (х, у, г ) ^ ^т) (x, У, г ^-¿т и(Г] (x, У, г )
V-т-2) (х, у, г > = -¿3 ит) (х, у, г ) +
(39)
(40)
(41)
(42)
( х,у г ) = и (х,у,г
2 2 2 х - у - г и(т)
+ ^ ихт) (х, у, г) +
+ -ёги(;] (х,у,+ (х,у,г), (43)
Vт2) (х, у, г ) =
_ у Т-Кт),
и(т) (х, у, г) + ихт) (х, у, г) +
-х2 + у2 - 2:2 (т) ( \ 2 у2 т т ( т) / \
+—г^т+3— иу (x,У, г)+Т^т+Тиг} (x,У, г), (44)
Vт2) (х, у, г ) =
г ^М (ху2)+ 2x2 ит)|
ит) (х,у,г)+-5т+3ихт) (х,у,г) +
+ Ту* и? (х,у, г)+-х2 -¿+ ^ игт) (х,у, г), (45) V -т-1) (х, у, г ) = -тг ихт) (х, у, г ) +
^-^т+г иут) (x, y, 2)+-mTигm) (х^ (46)
(т),
г
Д-т-1)
V(-т-1) (х, у, г ) =
+ ^^ ^ y, г) +(-х2 - у2 + г2 ^ (x, y, г), (34) и(-) (х, у, г)--и« (х, у, г)
2т+1 х V / 2т+1 у V ' ' /
г
г
V
(x y,z ) =
U(m)!
( X y,z )--
U(m)
( x, y, z ),
V (-m-1)( x, y, z ) =
t(m)
uym)( ^ y, z )--¿r UZm)( ^ y, z )
Ám) ,
(48)
(49)
7 2 2 2 « х + у + г , а верхнии индекс показывает, какая степень однородности должна быть у соответствующей однородной гармонической функции.
Как входная, так и выходная функции могут быть представлены с помощью подстановки Дон-кина (11), (12), а соответствующие им двумерные функции Н (р, q) и J (р, q) записаны в полярных
координатах (17) и представлены в виде ряда Фурье (19), где коэффициенты ряда определяются единственно возможным образом. Сравнение между собой коэффициентов Фурье справа и слева от знака равенства в тождествах (28)-(49) позволяет установить полезные тождества для радиальных функций Фт к (с точностью до множителя, выбираемого достаточно произвольным образом, поскольку и сами функции Фтк определены
с точностью до множителя). Однако практическая полезность таких соотношений сомнительна, тем более что все они фактически копируют соотношения, получаемые для решений (26) с помощью классических соотношений для гипергеометрических функций, которые можно найти, в частности, в [50-52].
Е I (ck cosk* + sk sink*) /k (l + /л2)
:=0,« I ^ ' \ /
1
xjY-
0
jTm+k (1 -Y)-m-1 (l +/ )-(1+m) d Y
0 i k=0,c
XY 1 - Y
= (1 +
ИНТЕГРАЛЬНАЯ ФОРМУЛА ДЛЯ ОДНОРОДНЫХ ГАРМОНИЧЕСКИХ ФУНКЦИИ СО СТЕПЕНЯМИ ОДНОРОДНОСТИ -1 < т < 0
Интегральная формула (23) для гипергеометрических функций после ее подстановки в решение (21) для функций Фтк и далее в разложение Фурье (19) дает формулу
фт (/*) = Е ф m,k (/)(ck cos k* + sk sin k*) =
m+1
+ Y2 ) x
= (1 + Л2 )m+ j| Е (ck cos k* + sk sin k*)
(1 -y)-m-1 (1 + Y//2 )-(1+m) dY =
(1 + Л2 )m+1 j ÍE лkYk (^ + ^ )]x
0 I k=0,« )
X Ym (1 - Y )-m1 (1 + Y/Л )-(1+m) dY, (50)
которую можно без каких-либо опасений использовать при -1 < m < 0. В этой формуле комплексные коэффициенты ck и sk являются не совсем произвольными, поскольку выбираются так, чтобы результирующее выражение было вещественным. Однако их можно считать совершенно произвольными, если в формуле (50) произвольные вещественные константы ck и sk заменить на произвольные комплексные константы, но при этом для получившегося комплексного выражения использовать лишь его вещественную часть.
Выражение Е /"y" (Cke+Vcv + ske~ík'p) в форму-
k=0,«
ле (50) можно упростить. Действительно,
Е С/ = Kc (t) = Kcr (u,v) + iKci (u,v),
k=0,«
E Sktk = Ks (t) = Ksr (u, v) + iKsi (u,v)
k=0,«
это произвольные аналитические функции комплексного переменного t = u + iv, разложенные в ряд Тейлора в окрестности точки t = 0 . Поэтому
Е YkYk (cke+^ + Ske-^) =
k=0,«
Е ck (/Ye+ Е s\ (/Ye-*)" =
=0,« k=0,«
= Kc (/Ye+* + K (/Ye-*) = = Kcr (Yy cos *, YY sin *) + iKci (/Y cos *, Y Y sin *) +
+ Ksr (/Y COs *, -/Y sin *) + iKsi (YY cos *, -/Y sin *) = = K(/Y cos *, YY sin *),
(51)
где K (u, v) — это произвольная вещественная функция, удовлетворяющая двумерному уравнению Лапласа. Действительно, функции Kci и Ksi в формуле (51) должны сократиться, чтобы получить на выходе вещественное значение. Из усло-
k=0,
вия
Kci (u,v) = -K. (u,-v)
z
x
r
r
k=0,«
и соотношений Коши—Римана [53-56], справедливых для пар функций Ксг (и, V), КС1 (и, V)
и Кг (и,V), Кя (и,V) , следует цепочка равенств
К (и, V) _ 8Ка (и,v)_ 8Кя (и, -V) _ дКг (и, -V)
ди
8v
dv
ди
8Ка (и, v)_ 8Ka (и, v)_8Кя (и, -v)_ 8Ks (и, -v)
8v
ди
ди
8v
откуда следует условие Кг (и,-V) = Ксг (и, V) + С,
где С — это вещественная константа. Следовательно, с точностью до аддитивной вещественной константы функции Кс ^) и К (t) должны быть
комплексно сопряженными: К (t) = Кс^ ) + С . Если ввести в рассмотрение новую функцию К (и, V) = 2КСг (и, V) + С, которая обязана будет удовлетворять двумерному уравнению Лапласа вследствие соотношений Коши—Римана для вещественной и мнимой частей аналитической функции комплексного переменного [53-56], то получим конечное равенство в цепочке (51). Поскольку для любой функции К(и,V) , удовлетворяющей двумерному уравнению Лапласа, можно восстановить подходящую аналитическую функцию К^) комплексного переменного t = и + IV, вещественной частью которой является заданная функция К (и, V) , то никакие дополнительные ограничения на функцию К(и,V) кроме ее гармоничности не накладываются.
Переходя от полярных координат (¡и,ф) к декартовым координатам (р, q) в соответствии с формулами (17), получаем интегральное выражение
Hm(p, q ) = (1 + p2 + q2 )m+1 x i
x J K(pr, qr)
-dr (52)
(1 -т)т+1(1 + т( Р2 + q 2 ))"
для решения эллиптического уравнения (10), когда это решение не имеет особых точек внутри круга р2 + q2 < 1. В этом выражении функция К (и, V)
должна быть произвольной двумерной гармонической функцией, не имеющей особых точек внутри круга и2 + V2 < 1.
В частном случае в качестве функций К (и, V) можно использовать гармонические однородные полиномы степени k, которые вычисляются по фор-
муле К(и, v) = Re (и + iv)к либо по формуле К(и,v) = Im (и + iv)к . Тогда формула (52) с точностью до замены переменных (p,q(и,ф) и константы-множителя будет порождать функции Фт,к (u)cosкф и Фm,k (¡)sinкф, принадлежащие
гипергеометрическому базису.
Прямая подстановка выражения (52) в уравнение (10) показывает, что формула (52) обеспечивает решение для уравнения (10) даже тогда, когда двумерная гармоническая функция К (и, v) имеет
особые точки внутри круга и2 + v2 < 1. Однако вопрос о том, в какой степени формула (52) дает все решения для уравнения (10), требует аккуратного исследования и будет являться предметом отдельной публикации.
ЗАКЛЮЧЕНИЕ
Полученные аналитические выражения для базисных функций (21) и (26), которые содержат гипергеометрические функции, после умножения их на тригонометрические функции cos(m^) и sin(m^), обратной замены полярных координат (и, ф) на переменные (p, q) с помощью формул (17), возвращения от вспомогательных переменных (p, q) к трехмерным декартовым координатам (x, y, z) с помощью подстановки Донкина (8) и окончательного восстановления трехмерных однородных гармонических функций U(x, y, z) степени т в соответствии с формулой (9) позволяют сконструировать для однородных гармонических функций степени т аналог тригонометрического базиса Фурье. С помощью этого базиса любая однородная гармоническая функция степени т, не имеющая "лишних" сингулярных точек, может быть разложена в сходящийся ряд. Использование для этой цели лишь тех базисных функций, которые порождаются формулами (21), позволяет получать однородные гармонические функции степени т, которые не имеют сингулярных точек вдоль луча x = 0, y = 0, z > 0, но при этом неизбежным образом имеют сингулярные точки вдоль луча x = 0, У = 0, z < 0.
Точно также из формулы (52) после подстановки в формулу (9) с учетом замены переменных (8) получается интегральная формула общего вида для трехмерных гармонических функций, однородных по Эйлеру со степенью однородности т, удовлетворяющих условию -1 < т < 0. С помощью процедуры дифференцирования [24, 27, 57] и формулы Томсона (7) либо напрямую с помощью дифференциальных операторов Томсона— Донкина [39-42] из этой формулы можно полу-
т
г
чить общие формулы и для других степеней однородности. Следует, однако, обратить внимание, что в интегральной формуле (52) при дифференцировании под знаком интеграла не надо забывать о предварительном выделении и аналитическом интегрировании сингулярностей интегрального ядра.
Авторы публично заявляют, что у них нет конфликта интересов, в том числе и финансовых. Данная работа частично выполнена в рамках НИР 0074-20190009, входящей в состав гос. задания № 075-00780-1902 Министерства науки и высшего образования Российской Федерации для ИАП РАН. Для проведения и проверки аналитических выкладок использовалась программа Wolfram Mathematica версии 11 (http://wolfram. com/mathematica/).
СПИСОК ЛИТЕРАТУРЫ
1. Бердников А.С., Краснова Н.К., Соловьёв К.В., Кузьмин А.Г., Масюкевич С.В., Титов Ю.А. Скрещенные гармонические потенциалы, однородные по Эйлеру // Научное приборостроение. 2019. Т. 29, № 4. С. 84-95.
2. Габдуллин П.Г., Голиков Ю.К., Краснова Н.К., Давыдов С.Н. Применение формулы Донкина в теории энергоанализаторов. I // Журнал технической физики. 2000. Т. 70, № 2. C. 91-94.
3. Габдуллин П.Г., Голиков Ю.К., Краснова Н.К., Давыдов С.Н. Применение формулы Донкина в теории энергоанализаторов. II // Журнал технической физики. 2000. Т. 70, № 3. С. 44-47.
4. Голиков Ю.К., Краснова Н.К. Теория синтеза электростатических энергоанализаторов. СПб.: Изд-во Политехнического университета, 2010. 409 с.
5. Голиков Ю.К., Краснова Н.К. Электрические поля, однородные по Эйлеру, для электронной спектрографии // Журнал технической физики. 2011. Т. 81, № 2. С. 9-15.
6. Краснова Н.К. Теория и синтез диспергирующих и фокусирующих электронно-оптических сред. Дис. ... д-ра физ.-мат. наук. СПб., 2013. 259 с.
7. Голиков Ю.К., Краснова Н.К. Аналитические структуры электрических обобщенно-однородных спектрографических сред // Научное приборостроение. 2014. Т. 24, №1. С. 50-58.
URL: http://iairas.ru/mag/2014/abst1 .php#abst6
8. Голиков Ю.К., Краснова Н.К. Обобщенный принцип подобия и его применение в электронной спектрографии // Прикладная физика. 2007. № 2. С. 5-11.
9. Аверин И.А., Бердников А.С., Галль Н.Р. Принцип подобия траекторий при движении заряженных частиц с разными массами в однородных по Эйлеру электрических и магнитных полях // Письма в Журнал технической физики. 2017. Т. 43. № 3. С. 39-43. DOI: 10.1134/S106378501702002X
10. Бердников А.С., Галль Л.Н., Антонов А.С., Соловьев К.В. Синтез краевых магнитных полей для статических масс-анализаторов спектрографического типа // Масс-спектрометрия. 2018. Т. 15, № 1. С. 26-43.
11. Голиков Ю.К., Бердников А.С., Антонов А.С., Красно-
ва Н.К., Соловьев К.В. Синтез электродных конфигураций, сохраняющих для краевых электрических полей свойство однородности по Эйлеру // Журнал технической физики. 2018. Т. 88, № 4. С. 609-613. DOI: 10.21883/JTF.2018.04.45732.2483
12. Голиков Ю.К., Бердников А.С., Антонов А.С., Краснова Н.К., Соловьев К.В. Применение формулы Донкина в теории электростатических призм // Журнал технической физики. 2018. Т. 88, № 11. С. 1711-1719. DOI: 10.21883/JTF.2018.11.46635.2498
13. Голиков Ю.К., Бердников А.С., Антонов А.С., Краснова Н.К., Соловьев К.В. Применение формулы Донкина в теории отражающих и поворотных устройств // Журнал технической физики. 2019. Т. 89, № 12. С. 1946-1963. DOI: 10.21883/JTF.2019.12.48496.201-18
14. Бердников А.С., Аверин И.А. О невозможности двойной фокусировки в комбинированных электрических и магнитных полях, однородных по Эйлеру // Масс-спектрометрия. 2016. Т. 13, № 1. С. 67-70.
15. Бердников А.С., Аверин И.А., Голиков Ю.К. Статические масс-спектрографы нового типа, использующие электрические и магнитные поля, однородные по Эйлеру. I. Общий принцип и однокаскадные схемы // Масс-спектрометрия. 2015. Т. 12, № 4. С. 272-281
16. Бердников А.С., Аверин И.А., Голиков Ю.К. Статические масс-спектрографы нового типа, использующие электрические и магнитные поля, однородные по Эйлеру. II. Условия двойной фокусировки высокого порядка у двухкаскадной схемы // Масс-спектрометрия. 2016. Т. 13, № 1. С. 11-20.
17. Бердников А.С., Аверин И.А. Новый подход к разработке ионно-оптических схем статических масс-спектрографов на основе неоднородных полей, однородных по Эйлеру // Успехи прикладной физики.
2016. Т. 4, № 1. С. 89-95.
18. Аверин И.А., Бердников А.С. Краевые поля бессеточных электронных спектрографов с однородными по Эйлеру электростатическими полями // Успехи прикладной физики. 2016. Т. 4, № 1. С. 5-8.
19. Фихтенгольц Г.М. Курс дифференциального и интегрального исчисления. Т. 1. М.: ФИЗМАТЛИТ, 2001. 616 с.
20. Смирнов В.И. Курс высшей математики. Т. 1. М.: Наука, 1974. 480 с.
21. Бердников А.С., Аверин И.А., Краснова Н.К., Соловьев К.В. Об однородности скалярных и векторных потенциалов электрических и магнитных полей, однородных по Эйлеру // Успехи прикладной физики.
2017. Т. 5, № 1. С. 10-27.
22. Голиков Ю.К., Уткин К.Г., Чепарухин В.В. Расчет элементов электростатических электронно-оптических систем. Учебное пособие. Л.: Изд-во ЛПИ, 1984. 79 с.
23. Голиков Ю.К., Соловьев К.В. Электростатические ионные ловушки. СПб.: Изд-во Политехнического университета, 2008. 152 с.
24. Бердников А.С., Аверин И.А., Краснова Н.К., Соловьев К.В. Общие формулы для трехмерных электрических и магнитных потенциалов, однородных по Эйлеру с целочисленным порядком однородности // Науч-
ное приборостроение. 2016. Т. 26, № 4. С. 13-30. URL: http://iairas.ru/mag/2016/abst4.php#abst2
25. Donkin W.F. On the Equation of Laplace's Functions &c. // Philosophical Transactions of the Royal Society of London. 1857. Vol. 147. P. 43-57.
DOI: 10.1098/rstl. 1857.0005
26. Donkin W.F. On the Equation of Laplace's Functions &c. // Proceedings of the Royal Society of London. 1856-1857. Vol. 8. P. 307-310.
27. Гобсон Е.В. Теория сферических и эллипсоидальных функций. М.: Изд-во иностранной литературы, 1952. 476 с.
28. Уиттекер Э.Т., Ватсон Дж. Курс современного анализа. Ч. 2: Трансцендентные функции. М.: ГИФМЛ, 1963. 516 с.
29. Thomson W. Extraits de deux Lettres adressées à M. Liouville // Journal de mathématiques pures et appliquées. 1847. T. XII. P. 256-264.
30. Томсон У. (лорд Кельвин), Тэт П.Г. Трактат по натуральной философии. Ч. I. Москва-Ижевск: НИЦ "Регулярная и хаотическая динамика", 2010. 572 с.
31. Томсон У. (лорд Кельвин), Тэт П.Г. Трактат по натуральной философии. Ч. II. Москва-Ижевск: НИЦ "Регулярная и хаотическая динамика", 2011. 560 с.
32. Сретенский Л.Н. Теория ньютоновского потенциала. М.-Л.: ОГИЗ-ГИТТЛ, 1946. 318 с.
33. Смирнов В.И. Курс высшей математики. Т. 4. Ч. 2. М.: Наука, 1981. 297 с.
34. Владимиров В.С. Уравнения математической физики. М.: Наука, 1981. 512 с.
35. Kellogg O.D. Foundations of Potential Theory. SpringerVerlag, Berlin, Heidelberg, New-York, 1967. 386 p.
36. Уэрмлер Дж. Теория потенциала. М.: Изд-во "Мир", 1980. 134 с.
37. Helms L.L. Potential Theory. Second Edition. Springer, London, Heidelberg, New-York, Dordrecht, 2014. 485 p.
38. Голиков Ю.К. Аналитические способы описания гармонических функций // Вестник Актюбинского регионального государственного университета им. К. Жу-банова. Физико-математические науки. 2016. № 2. С. 165-181.
39. Бердников А.С., Галль Л.Н., Галль Н.Р., Соловьев К.В. Обобщение формулы Томсона для гармонических функций общего вида // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2019. Т. 12, № 2. С. 32-48. DOI: 10.18721/JPM.12203
40. Бердников А.С., Галль Л.Н., Галль Н.Р., Соловьев К.В. Обобщение формулы Томсона для гармонических однородных функций // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2019. Т. 12, № 2. С. 49-62. DOI: 10.18721/JPM.12204
41. Бердников А.С., Галль Л.Н., Галль Н.Р., Соловьев К.В. Базисные дифференциальные операторы Донкина для однородных гармонических функций // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2019. Т. 12, № 3. С. 26-44. DOI: 10.18721/JPM.12303
42. Бердников А.С., Галль Л.Н., Галль Н.Р., Соловьев К.В. Дифференциальные операторы Донкина для однородных гармонических функций // Научно-технические ведомости СПбГПУ. Физико-математические науки. 2019. Т. 12, № 3. С. 45-62. DOI: 10.18721/JPM.12304
43. Бердников А.С., Аверин И.А., Краснова Н.К., Соловьев К.В. Простейшие аналитические электрические и магнитные потенциалы, однородные по Эйлеру // Вестник Актюбинского регионального государственного университета им. К. Жубанова. 2016. № 2. С. 17-32.
44. Бердников А.С., Аверин И.А., Краснова Н.К., Соловьев К.В. Трехмерные электрические и магнитные потенциалы, однородные по Эйлеру // Вестник Актю-бинского регионального государственного университета им. К. Жубанова. 2016. № 2. С. 147-165.
45. Бердников А.С., Аверин И.А., Краснова Н.К., Соловьев К.В. Интегральные формулы для трехмерных электрических и магнитных потенциалов, однородных по Эйлеру с нецелочисленными порядками однородности // Научное приборостроение. 2016. Т. 26, № 4. С. 31-42. URL: http://iairas.ru/mag/2016/abst4.php#abst3
46. Бердников А.С., Аверин И.А., Краснова Н.К., Соловьев К.В. Квазиполиномиальные трехмерные электрические и магнитные потенциалы, однородные по Эйлеру // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета. Физико-математические науки. 2017. Т. 10, № 1. С. 71-80. DOI: 10.18721/JPM.10107
47. Краснова Н.К., Бердников А.С., Соловьев К.В., Аверин И.А. О квазиполиномиальных трехмерных потенциалах электрических и магнитных полей // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета. Физико-математические науки. 2017. Т. 10, № 1. С. 81-92. DOI: 10.18721/JPM.10108
48. Бердников А.С., Краснова Н.К., Соловьев К.В. Анализ интегральной формулы Уиттекера общего вида для электрических и магнитных потенциалов, однородных по Эйлеру // Научное приборостроение. 2017. Т. 27, № 4. С. 63-71.
URL: http://iairas.ru/mag/2017/abst4.php#abst8
49. Бердников А.С., Краснова Н.К., Соловьев К.В. Интегральная формула Уиттекера для электрических и магнитных потенциалов с нулевым порядком однородности и ее следствия // Научное приборостроение. 2017. Т. 27, № 4. С. 72-89.
URL: http://iairas.ru/mag/2017/abst4.php#abst9
50. Абрамовиц М., Стиган И. (ред.). Справочник по специальным функциям с формулами, графиками и математическими таблицами. М.: Наука, 1979. 832 с.
51. Бейтмен Г., Эрдейи А. Высшие трансцендентные функции. Т. 1. Изд. 2-е. М.: Наука, 1973. 296 с.
52. Aomoto K., Kita M. Theory of Hypergeometric Functions (Springer Monographs in Mathematics Series, Vol. 305). Springer, 2011. 317 p.
ISSN 0868-5886
NAUCHNOE PRIBOROSTROENIE, 2019, Vol. 29, No. 4, pp. 96-109
53. Лаврентьев М.А., Шабат Б.В. Методы теории функций комплексного переменного. М.: Наука, 1965. 716 с.
54. Маркушевич А.И. Теория аналитических функций. М.: Наука, 1968. Т. 1. 486 с., т. 2. 624 с.
55. Евграфов М.А. Аналитические функции. 3-е изд. пе-рераб., доп. М.: Наука, 1991. 448 с.
56. Гурвиц А., Курант P. Теория функций. М.: Наука, 1968. 646 с.
57. Бердников А.С., Краснова Н.К., Соловьев К.В. Теорема о дифференцировании трехмерных электрических и магнитных потенциалов, однородных по Эйлеру // Научное приборостроение. 2017. Т. 27, № 3. С. 107119. URL: http://iairas.ru/mag/2017/abst3.php#abst13
Институт аналитического приборостроения РАН, Санкт-Петербург (Бердников А.С., Соловьёв К.В., Кузьмин А.Г., Масюкевич С.В., Титов Ю.А., Голиков Ю.К.)
Санкт-Петербургский политехнический университет Петра Великого (Краснова Н.К., Соловьёв К.В., Голиков Ю.К.\)
Контакты: Бердников Александр Сергеевич, asberd@yandex.ru
Материал поступил в редакцию 23.09.2019
HYPERGEOMETRIC BASIS FOR THREE-DIMENSIONAL HARMONIC FUNCTIONS WHICH ARE HOMOGENEOUS IN EULER TERMS WITH A NON-INTEGER POWER OF HOMOGENEITY
1 2 12 1 A. S. Berdnikov , N. K. Krasnova , K. V. Solovyev , A. G. Kuzmin ,
S. V. Masyukevich1, Yu. A. Titov1, Yu.K. Golikpy12
1 Institute for Analytical Instrumentation of RAS, Saint Petersburg, Russia 2Peter The Great Saint Petersburg Polytechnic University, Russia
Electric and magnetic fields which are homogeneous in Euler terms, are a convenient tool for the synthesis of electron and ion optical systems with special properties. It is known that the scalar potentials of such fields are three-dimensional scalar harmonic functions which are homogeneous in Euler terms with a given power of homogeneity. The problem of the exhaustive parameterization of three-dimensional homogeneous harmonic functions with integer powers of homogeneity is already solved using the Donkin formulas for homogeneous harmonic functions with powers of homogeneity 0 and -1, the theorem on the differentiation of three-dimensional homogeneous harmonic functions and Thomson's formula for homogeneous harmonic functions. However, the number of analytical formulas that can be used to describe three-dimensional scalar harmonic functions with non-integer powers of homogeneity is, unfortunately, not very large at the moment, and the problem of an exhaustive description of such functions is very far from its final solution. At the same time, the usage of three-dimensional homogeneous harmonic potentials with non-integer powers of homogeneity significantly expands the toolkit of developers of electron and ion optical systems. The goal of this work is to construct a hypergeometric basis composed of basic homogeneous harmonic functions with non-integer powers of homogeneity, with the help of which any three-dimensional homogeneous harmonic function without singular points exept the line x = y = 0, z < 0 can be represented as an infinite series like the Fourier series. It seems that this result partially solves the problem of an exhaustive description of three-dimensional scalar homogeneous harmonic functions with non-integer powers of homogeneity.
Keywords: electric fields; harmonic functions; functions homogeneous in Euler' terms; similarity principle for charged particle trajectories; Donkin formula; analytical solutions of Laplace equation
potentials which are homogeneous in Euler' terms]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2019, vol. 29, no. 4, pp. 84-95 (In Russ.). 2. Gabdullin P.G., Golikov Yu.K., Krasnova N.K., Davy-
REFERENCES
1. Berdnikov A.S., Krasnova N.K., Solovyev K.V., Kuzmin A.G., Masyukevich S.V., Titov Yu.A. [Crossed harmonic
dov S.N. [The use of Donkin's formula in the theory of energy analyzers. I]. Zhurnal tekhnicheskoj fiziki [Journal of Applied Physics], 2000, vol. 70, no. 2, pp. 91-94. (In Russ.).
3. Gabdullin P.G., Golikov Yu.K., Krasnova N.K., Davydov S.N. [The use of Donkin's formula in the theory of energy analyzers. II]. Zhurnal tekhnicheskoj fiziki [Journal of Applied Physics], 2000, vol. 70, no. 3, pp. 44-47. (In Russ.).
4. Golikov Yu.K., Krasnova N.K. Teoriya sinteza ehlek-trostaticheskih ehnergoanalizatorov [Theory of synthesis of electrostatic power analyzers]. Saint Petersburg, Polytechnical university Publ., 2010. 409 p. (In Russ.).
5. Golikov Yu.K., Krasnova N.K. [Application of electric fields uniform in the Euler sense in electron spectrogra-phy]. Zhurnal tekhnicheskoj fiziki [Journal of Applied Physics], 2011, vol. 81, no. 2, pp. 9-15. (In Russ.).
6. Krasnova N.K. Teoriya i sintez dispergiruyushchih i fokusi-ruyushchih ehlektronno-opticheskih sred. Diss. dokt. fiz.-mat. nauk [The theory and synthesis of the dispersing and focusing electron-optical environments. Dr. phys. and math. sci. diss.]. Saint Petersburg, 2013. 259 p. (In Russ.).
7. Golikov Yu.K., Krasnova N.K. [Analytical structures of electric spectrographs the fields of which are expressed in a uniform generalized form]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2014, vol. 24, no. 1, pp. 50-58. (In Russ.).
URL: http://iairas.ru/en/mag/2014/abst1.php#abst6
8. Golikov Yu.K., Krasnova N.K. [The generalized principle of similarity and its application in electronic spec-trography]. Prikladnaya fizika [Applied physics], 2007, no. 2, pp. 5-11. (In Russ.).
9. Averin I.A., Berdnikov A.S., Gall N.R. [The principle of similarity of trajectories for the motion of charged particles with different masses in electric and magnetic fields that are homogeneous in Euler terms]. Pis'ma v ZhTF [Letters in ZhTF], 2017, vol. 43, no. 3, pp. 39-43.
DOI: 10.1134/S106378501702002X (In Russ.).
10. Berdnikov A.S., Gall L.N., Antonov A.S., Soloviev K.V. [Synthesis of Fringing Magnetic Fields for Static Mass Analyzers of the Spectrographic Type]. Mass-spektrometriya [Mass-spectrometry], 2018, vol. 15, no. 1, pp. 26-43. (In Russ.).
11. Golikov Yu.K., Berdnikov A.S., Antonov A.S., Krasnova N.K., Solov'ev K.V. [Synthesis of Electrode Configurations that Conserve Fringing Electric Field Homogeneity in Euler Terms]. Zhurnal tekhnicheskoj fiziki [Journal of Applied Physics], 2018, vol. 88, no. 4, pp. 609-613. (In Russ.). DOI: 10.21883/JTF.2018.04.45732.2483
12. Golikov Yu.K., Berdnikov A.S., Antonov A.S., Krasnova N.K., Solov'ev K.V. [Application of the Donkin Formula in the Theory of Electrostatic Prisms]. Zhurnal tekhnicheskoj fiziki [Journal of Applied Physics], 2018, vol. 88, no. 11, pp. 1711-1719. (In Russ.). DOI: 10.21883/JTF.2018.11.46635.2498
13. Golikov Yu.K., Berdnikov A.S., Antonov A.S., Krasnova N.K., Solovyev K.V. [Applications of Donkin formula in the theory of reflecting and turning devices]. Zhurnal tekhnicheskoj fiziki [Journal of Applied Physics], 2019, vol. 89, no. 12, pp. 1946-1963. (In Russ.). DOI:
10.21883/JTF.2019.12.48496.201-18
14. Berdnikov A.S., Averin I.A. [About impossibility of double focusing in the combined electric and magnetic fields uniform in Euler]. Mass-spektrometriya [Mass-spectrometry], 2016, vol. 13, no. 1, pp. 67-70. (In Russ.).
15. Berdnikov A.S., Averin I.A., Golikov Yu.K. [The static mass spectrogaphs of new type using the electric and magnetic fields uniform in Euler. I]. Mass-spektrometriya [Mass-spectrometry], 2015, vol. 12, no. 4, pp. 272-281. (In Russ.).
16. Berdnikov A.S., Averin I.A., Golikov Yu.K. [The static mass spectrogaphs of new type using the electric and magnetic fields uniform in Euler. II. Conditions of highorder double focusing for two-cascade schemes]. Mass-spektrometriya [Mass-spectrometry], 2016, vol. 13, no. 1, pp. 11-20. (In Russ.).
17. Berdnikov A.S., Averin I.A. [New approach to development of ion-optical schemes of static mass spectrographs on the basis of the non-uniform magnetic fields uniform in Euler]. Uspekhi prikladnojfiziki [Achievements of applied physics], 2016, vol. 4, no. 1, pp. 89-95. (In Russ.).
18. Averin I.A., Berdnikov A.S. [Fringing fields of gridless electron spectrographs with electrostatic fields homogeneous in Euler terms]. Uspekhi prikladnoj fiziki [Achievements of applied physics], 2016, vol. 4, no. 1, pp. 5-8 (In Russ.).
19. Fihtengolc G.M. Kurs differencial'nogo i integral'nogo ischisleniya [Course of differential and integral calculus], vol. 1. Moscow, Fizmatlit Publ., 2001. 616 p. (In Russ.).
20. Smirnov V.I. Kurs vysshej matematiki [Course of the higher mathematics], vol. 1. Moscow, Nauka Publ., 1974. 480 p. (In Russ.).
21. Berdnikov A.S., Averin I.A., Krasnova N.K., Solovyev K.V. [On homogeneity of scalar and vector potentials of electric and magnetic fields, homogeneous in Euler terms]. Uspekhi prikladnoj fiziki [Achievements of applied physics], 2017, vol. 5, no. 1, pp. 10-27. (In Russ.).
22. Golikov Yu.K., Utkin K.G., Cheparuhin V.V. Raschet elementov elektrostaticheskich elektronno-opticheskich system. Uchebnoe posobie [Calculation of elements of electrostatic electron-optical systems. Education book]. Leningrad, LPI Publ., 1984. 79 p. (In Russ.).
23. Golikov Yu.K., Solovyov K.V. Elektrostaticheskie ionnye lovushki [Electrostatic ion traps]. Saint Petersburg, Saint-Petersburg Polytechnic University Publ., 2008. 152 p. (In Russ.).
24. Berdnikov A.S., Averin I.A., Krasnova N.K., Solovyev K.V. [Universal expressions for 3d electric and magnetic potentials which are uniform in Euler terms]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2016, vol. 26, no. 4, pp. 13-30. DOI: 10.18358/np-26-4-i1330 (In Russ.).
25. Donkin W.F. On the Equation of Laplace's Functions &c. Philosophical Transactions of the Royal Society of London, 1857, vol. 147, pp. 43-57.
DOI: 10.1098/rstl. 1857.0005
26. Donkin W.F. On the Equation of Laplace's Functions &c. Proceedings of the Royal Society of London, 1856-1857, vol. 8, pp. 307-310. DOI: 10.1098/rspl.1856.0075
27. Gobson E.V. Teoriya sfericheskih i ellipsoidal'nyh funkcij [Theory of spherical and ellipsoidal functions]. Moscow, Inostrannaya literatura Publ., 1952. 476 p. (In Russ.).
28. Whittaker E.T., Watson G.N. A course of modern analy-
sis. Cambridge, University Press, 1920. 608 p. (Russ ed.: Uitteker E.T., Vatson G. Kurs sovremennogo analiza. Chast 2. Transcendentnye funkcii. Moscow, GIFML Publ., 1963. 516 p.). (In Russ.).
29. Thomson W. Extraits de deux Lettres adressées à M. Liouville. Journal de mathématiques pures et appliquées, 1847, tome XII, pp. 256-264.
30. Thomson W, (later Lord Kelvin), Tait P.G. Treatise on Natural Philosophy. Oxford University Press, 1867. (Russ. Ed.: Tomson U. (Lord Kelvin), Tet P.G. Traktat po natural'noj filosofii. Part. I. Moscow-Izhevsk, NIC "Re-gulyarnaya i haoticheskaya dinamika", 2010. 572 p.). (In Russ.).
31. Thomson W, (later Lord Kelvin), Tait P.G. Treatise on Natural Philosophy. Oxford University Press, 1867. (Russ. Ed.: Tomson U. (Lord Kelvin), Tet P.G. Traktat po natural'noj filosofii. Part. II. Moscow-Izhevsk, NIC "Re-gulyarnaya i haoticheskaya dinamika", 2011. 560 p.). (In Russ.).
32. Sretenskiy L.N. Teoriya nyutonovskogo potenciala [Theory of Newtonian potential]. Moscow-Leningrad, OGIZ-GITTL Publ., 1946. 318 p. (In Russ.).
33. Smirnov V.I. Kurs vysshej matematiki [Course of the higher mathematics], vol. 4, part 2. Moscow, Nauka Publ., 1981. 297 p. (In Russ.).
34. Vladimirov V.S. Uravneniya matematicheskoj fiziki [Equations of mathematical physics]. Moscow, Nauka Publ., 1981. 512 p. (In Russ.).
35. Kellogg O.D. Foundations of Potential Theory. SpringerVerlag, Berlin, Heidelberg, New-York, 1967. 386 p.
36. Wermer J. Potential theory. Second edition. Ser. "Lecture Notes in Mathematics", vol. 408. Springer-Verlag, Berlin, Heidelberg, New-York, 1981. 166 p. (Russ. Ed.: Uermler G. Teoriya potenciala. Moscow, Mir Publ., 1980. 134 p.). (In Russ.).
37. Helms L.L. Potential Theory. Second Edition. Springer, London, Heidelberg, New-York, Dordrecht, 2014. 485 p.
38. Golikov Yu.K. [Analytical ways of describing harmonic functions]. Vestnik Aktyubinskogo regional'nogo gosu-darstvennogo universiteta im. K. Zhubanova. Fiziko-matematicheskie nauki [Bulletin of the Aktyubinsk regional state university of K. Zhubanov. Physical and mathematical sciences], 2016, no. 2, pp. 165-181. (In Russ.).
39. Berdnikov A.S., Gall L.N., Gall N.R., Solovyev K.V. [Generalization of the Thomson formula for general harmonic functions]. Nauchno-tekhnicheskie vedomosti SPbGPU. Fiziko-matematicheskie nauki [Scientific and technical sheets SPbGPU. Physical and mathematical sciences], 2019, vol. 12, no. 2, pp. 32-48. DOI: 10.18721/JPM.12203 (In Russ.).
40. Berdnikov A.S., Gall L.N., Gall N.R., Solovyev K.V. [Generalization of the Thomson formula for homogeneous harmonic functions]. Nauchno-tekhnicheskie vedomosti SPbGPU. Fiziko-matematicheskie nauki [Scientific and technical sheets SPbGPU. Physical and mathematical sciences], 2019, vol. 12, no. 2, pp. 49-62. DOI: 10.18721/JPM.12204 (In Russ.).
41. Berdnikov A.S., Gall L.N., Gall N.R., Solovyev K.V. [Basic Donkin's differential operators for homogeneous harmonic functions]. Nauchno-tekhnicheskie vedomosti
SPbGPU. Fiziko-matematicheskie nauki [Scientific and technical sheets SPbGPU. Physical and mathematical sciences], 2019, vol. 12, no. 3, pp. 26-44. (In Russ.). DOI: 10.18721/JPM.12303
42. Berdnikov A.S., Gall L.N., Gall N.R., Solovyev K.V. [Donkin's differential operators for homogeneous harmonic functions]. Nauchno-tekhnicheskie vedomosti SPbGPU. Fiziko-matematicheskie nauki [Scientific and technical sheets SPbGPU. Physical and mathematical sciences], 2019, vol. 12, no. 3, pp. 45-62. (In Russ.). DOI: 10.18721/JPM.12304
43. Berdnikov A.S., Averin I.A., Krasnova N.K., Solovyev K.V. [The elementary analytical electric and magnetic potentials uniform in Euler]. Vestnik Aktyubinskogo regional'nogo gosudarstvennogo universiteta im. K. Zhubanova. Fiziko-matematicheskie nauki [Bulletin of the Aktyubinsk regional state university of K. Zhubanov. Physical and mathematical sciences], 2016, no. 2, pp. 17-32. (In Russ.).
44. Berdnikov A.S., Averin I.A., Krasnova N.K., Solovyev K.V. [Three-dimensional electric and magnetic potentials uniform in Euler]. Vestnik Aktyubinskogo regional'nogo go-sudarstvennogo universiteta im. K. Zhubanova. Fiziko-matematicheskie nauki [Bulletin of the Aktyubinsk regional state university of K. Zhubanov. Physical and mathematical sciences], 2016, no. 2, pp. 147-165. (In Russ.).
45. Berdnikov A.S., Averin I.A., Krasnova N.K., Solovyev K.V. [Integral expressions for 3D electric and magnetic potentials which are uniform in Euler terms and have non-integer orders of uniformity]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2016, vol. 26, no. 4, pp. 31-42.
DOI: 10.18358/np-26-4-i3142 (In Russ.).
46. Berdnikov A.S., Averin I.A., Krasnova N.K., Solovyev K.V. [Quasi-Polynomial 3D Electric and Magnetic Potentials Homogeneous in Euler Terms]. Nauchno-tekhnicheskie vedomosti SPbGPU. Fiziko-matematicheskie nauki [Scientific and technical sheets SPbGPU. Physical and mathematical sciences], 2017, vol. 10, no. 1, pp. 71-80. (In Russ.). DOI: 10.18721/JPM.10107
47. Krasnova N.K., Berdnikov A.S., Solovyev K.V., Averin I.A. [On the Quasi-Polynomial 3D Potentials of Electric and Magnetic Fields]. Nauchno-tekhnicheskie vedomosti SPbGPU. Fiziko-matematicheskie nauki [Scientific and technical sheets SPbGPU. Physical and mathematical sciences], 2017, vol. 10, no. 1, pp. 81-92. (In Russ.). DOI: 10.18721/JPM.10108
48. Berdnikov A.S., Krasnova N.K., Solovyev K.V. [Analysis of the general whittaker' formula for 3D electric and magnetic potentials homogeneous in Euler terms]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2017, vol. 27, no. № 4, pp. 63-71. DOI: 10.18358/np-27-4-i6371 (In Russ.).
49. Berdnikov A.S., Krasnova N.K., Solovyev K.V. [Whittaker' formula for 3D electric and magnetic potentials with a zero order of homogeneity and its consequencies]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2017, vol. 27, no. 4, pp. 72-89. DOI: 10.18358/np-27-4-i7289 (In Russ.).
50. Abramowitz M., Stegun I. A. (eds.) Handbook on Mathematical Functions With Formulas, Graphs, and Mathematical Tables. Applied Mathematics Series, vol. 55, 10th
ed. USA National Bureau of Standards, 1972. 1046 p. (Russ. Ed.: Abramovic M., Stigan I. (red.). Spravochnik po specialnym funkciyam s formulami, grafikami i mate-maticheskimi tablicami. Moscow, Nauka Publ., 1979. 832 p.). (In Russ.).
51. Bateman H., Erdélyi A. Higher Transcendental Functions, vol. I. New York, Toronto, London, McGraw-Hill Book Company, 1953. 302 p. (Russ. Ed.: Beytmen G., Erdeji A. Vysshie transcendentnye funkcii. T. 1. Izd. 2. Moscow, Nauka Publ., 1973. 296 p.). (In Russ.).
52. Aomoto K., Kita M. Theory of Hypergeometric Functions. Springer Monographs in Mathematics Series, vol. 305. Springer, 2011. 317 p. DOI: 10.1007/978-4-431-53938-4
53. Lavrentev M.A., Shabat B.V. Metody teorii funkcij kom-pleksnogo peremennogo [Methods of complex variable function theory]. Moscow, Nauka Publ., 1965. 716 p. (In
Contacts: Berdnikov Aleksandr Sergeevich, asberd@yandex. ru
Russ.).
54. Markushevich A.I. Teoriya analiticheskih funkcij. Т. 1, 2. [Theory of analytical functions. V. 1, 2]. Moscow, Nauka Publ., 1968. 486 p., 624 p. (In Russ.).
55. Evgrafov M.A. Analiticheskie funkcii. 3-e izd. pererab., dop [Analytical functions. 3rd ed.]. Moscow, Nauka Publ., 1991. 448 p. (In Russ.).
56. Gurvic A., Kurant P. Teoriya funkcij [Theory of functions]. Moscow, Nauka Publ., 1968. 646 p. (In Russ.).
57. Berdnikov A.S., Krasnova N.K., Solovyev K.V. [Theorem on integration and differentiation of 3D electric and magnetic potentials which are homogeneous in Euler terms]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2017, vol. 27, no. 3, pp. 107-119. DOI: 10.18358/np-27-3-i107119 (In Russ.).
Article received by the editorial office on 23.09.2019