Модифицированный метод последовательных конформных отображений наперед заданных многоугольных областей Текст научной статьи по специальности «Математика»

2016 Математика и механика № 1(39)

УДК 517.54

DOI 10.17223/19988621/39/3

В.М. Радыгин, И.С. Полянский

МОДИФИЦИРОВАННЫЙ МЕТОД ПОСЛЕДОВАТЕЛЬНЫХ КОНФОРМНЫХ ОТОБРАЖЕНИЙ НАПЕРЕД ЗАДАННЫХ МНОГОУГОЛЬНЫХ ОБЛАСТЕЙ

Разработан метод решения обратной задачи Кристоффеля - Шварца - конформное отображение наперед заданной многоугольной области на каноническую область - единичный круг. В основе метода лежит применение обобщенной функции Жуковского для решения модельных задач отображения на верхнюю полуплоскость луночки, в общем случае заданной сектором произвольной области с гладкой границей в расширенной плоскости комплексного переменного. Работоспособность решения проверена на конкретных примерах.

Ключевые слова: конформное отображение, многоугольная область, единичный круг, обратная задача Кристоффеля - Шварца.

Большое число работ [1-4 и др.] рассматривают решение задачи конформного отображения канонической области, например полуплоскости или единичного круга, на односвязную область расширенной комплексной плоскости, ограниченную ломаной, - многоугольную область. При этом наибольший практический интерес широкого класса прикладных задач вызывает решение обратной задачи (конформное отображение наперед заданной многоугольной области). Поскольку дополнение методов прямого отображения (интеграл Кристоффеля - Шварца [4]) обратными [5] существенно расширяет эффективное применение конформного отображения при решении различных краевых задач физики и техники [1], например в постановке [6] барицентрического метода [7, 8]. Поэтому целью статьи является разработка метода конформного отображения наперед заданных многоугольных областей.

1. Постановка задачи

Определим последовательность неповторяющихся N точек А1,А2,...,Лы комплексной плоскости С, которые являются вершинами ^-угольной области Ос С с заданными значениями величины внутренних углов а1,а2,...,аN , соответствующих по п-му индексу (п = 1, N) точкам Ап , и таких, что 0 <ап <п,

N а

V—п = N - 2. Геометрический центр области О совпадает с началом координат

п=1 П

комплексной плоскости. Нумерацию точек Ап = гпе1фп определим в порядке положительного обхода области О : ф1 <ф2 <... . Известно, что для заданных точек Ьп , расположенных на единичной окружности О с С и являющихся про-

образами вершин Ап, конформное отображение единичного круга < 1 на область О определяется известным интегралом Кристоффеля - Шварца [3]:

*=f (t )=с j

П(1 -V b„)

.n=1

a J П-1

d т + C1

(1)

где z eQ ; C0 Ф 0 и C1 - комплексные константы.

С учетом использованных соотношений найдем обратное отображение области Q на единичный круг |t| < 1, т. е. функцию t = f_1 (z) для заданного интеграла Кристоффеля - Шварца (1).

2. Построение обратного отображения

Известно [2, 3 и др.], что прямое ш = ( z -i)/(z + i) и обратное z = i (1 + ю)/(1 -ю) преобразования Кэли переводят любую прямую в дугу единичной окружности и обратно. Используя дробно-линейное преобразование

Z = i ( R + z )/( R - z ), (2)

где R - минимальный радиус окружности Q', содержащей в себе все вершины An многоугольника, отобразим многоугольную область z eQ на верхнюю полуплоскость Im Z > 0 с вырезанными круговыми двуугольниками (круговые

луночки) [2]. При этом вершины, расположенные на границе окружности Q', будут отображаться в точки действительной оси an = - sin (фп )/(1 - cos (фп))

полуплоскости. Вершины, расположенные внутри окружности Q', будут отображаться в точки an , отступающие от действительной оси на величину

(R2 -|An|2)Д|An|2 + R2 -2|An\cos(фп)) , а от мнимой оси соответственно на величину -2R|An|sin(фп)j(|An12 + R2 -2|An\cos(фп)) . На рис. 1 приведен пример отображения семиугольника на верхнюю полуплоскость дробно-линейным отображением (2).

Im z

Im Z

Рис. 1. Пример отображения семиугольника на верхнюю полуплоскость дробно-линейным отображением (2)

Для последующего отображения полученной области с вырезанными круговыми двуугольниками (см., например, рис. 1) на верхнюю полуплоскость применим последовательность отображений. Подобный метод с учетом рассмотренного обобщения в [1] принято называть методом последовательных конформных отображений (метод исчерпывания). Для реализации указанного метода рассмотрим решение трех модельных задач (обратных отображений).

3. Модельные задачи обратного отображения

Первой является задача на построение отображения полуплоскости с выброшенным сегментом круга на полуплоскость [3]. В ней указанное отображение производится для ребра en = {An, An+1} исходного многоугольника, заданного парой соседних вершин An и An+1, расположенных на границе окружности Q'. Решение задачи отображения полуплоскости с выброшенным сегментом круга на полуплоскость без вырезов известно [1 - 3] и производится обобщенной функцией Жуковского. Параметры удаляемого n-го сегмента: радиус Rn и положение центра координат круга mn определяются по трем точкам, принадлежащим ребру en . При этом добавление произвольного кругового сегмента в полуплоскость обобщенной функцией Жуковского приведет к изменению остальных круговых двуугольников.

Лемма 1. Применение обобщенной функцией Жуковского w(Z)=S3 (S2 (Sj (Z))),

где S (Z) = (Z-a )/(Z + a); S2 (Z) = ZY; S3 (Z) = ay(1 + Z)/(1 -Z) при yeR + , для

добавления n-го (n = 1, N) кругового сегмента в верхнюю полуплоскость Im Z> 0

c N выброшенными круговыми сегментами приводит к отображению исходной области на верхнюю полуплоскость Im та> 0 c N -1 выброшенными эллиптическими сегментами.

Доказательство. Известно [9], что применение любого дробно-линейного преобразования к окружности приводит к ее отображению на другую окружность, при этом под термином «окружность» понимается не только окружность в собственном смысле слова, но и прямая. С учетом этого следует рассматривать влияние степенного отображения S2 (Z) = ZY на форму N -1 выброшенных сегментов круга из верхней полуплоскости Im Z > 0 при решении задачи отображения n-й (n = 1, N) круговой луночки.

Пусть функция S2 (Z) задана степенным рядом

ад

S2 (Z) = Zy=S CnZn; Cn = an + ibn. (3)

n=0

Тогда, представив степенную функцию равенством ZY = eY 1n(Z) = ^(y In (Z))/n!,

n=0

где Z= r (cos (ф) + i sin (ф)), найдем коэффициенты an и bn :

nj 2 /

an = (y/r)n cos(nаrctg[(фc0sф - Inr sinф)/(1пr cosф + фsinф)])(2 r + ф2) /n!;

nj 2 /

bn =(y/r)n sin (n arctg[(фc0sф- 1nr sinф)/(1nr cos ф + фsin ф)])(2 r + ф2) n!. (4)

С учетом полученных соотношений (5) ряд (4) будет иметь вид

ад

S2 (Z) = S [уnPn (r, ф)( (Фп (г, ф)) + i sin(ф n (r, ф)))];

n=0

Pn (r, ф) = (in2 r +ф2) Yn!; (5)

Фп (r, ф) = n [arctg ((ф cos ф-in r sin ф)/( in r cos ф + ф sin ф)) + ф].

Известно [10], что применение степенного преобразования к кривой второго порядка не приводит к повышении порядка кривой. Пусть на границе произвольного k -го (k, n = 1, N, Re mk < Re mn, k Ф n) кругового сегмента с цент г r\ i г'ф i + 2ln

тром в mk исходной полуплоскости ImZ> 0 задано пять точек hj = j ' (У i = 0,1,2,...;j = 1,5) при условии, что точки hx и h5 расположены на противоположных концах выбранной k-й дуги, т. е. ф2 = ф5 = 0. Тогда коэффициенты уравнения кривой второго порядка в декартовой системе

anx2 + a22 y2 + 2a12 xy + 2a13 x + 2a23 y = 1 (6)

задаются равенствами

% = \H\/\M\; % = 0,5| H2\¡\M\; a22 = \H3\¡\M\;

a13 = 0,51H41/|M|; a23 = 0,51H5\¡|M|, (7)

где M = (mj i )5 5 - матрица, элементы которой определяются j-ми точками кривой для i = 0 с учетом выражения (5) в виде

/ N Л2 í N

MJ,1 =1 Srnp j,n cosф j,n I ; Mi,3 =1 Synp j.n sinф j,n

V n=0 / V n=0

NN

Mj,4 = S Ynp j,n cos ф j,n ; Mj,5 = S Ynp j,n sinф j,n ; (8)

n=0 n=0

N N

Mj,2 = SYn P j ,n cos ф j,n SYn P j ,n sin ф j,n , при N .

n=0 n=0

Hk = (Hj,) - матрица, формируемая матрицей M при замене ее k-го столбца

вектором е =(1 1 1 1 1)Г ; pj n = Pn (rj, фj); фj,n = Фп (rj,фj).

Из (8) следует, что при фиксированном N в матрице M элементы 1-го и 3-го столбца можно представить суммой из N + ^ ^ j слагаемых, элементы 2-го столбца - суммой из N2 слагаемых, элементы 4-го и 5-го столбцов - суммой из N слагаемых. Поскольку в указанном случае элементы j-х столбцов - суть суммы равного числа слагаемых, то определитель матрицы M можно представить в виде

10N+2[ N j+N2-14

суммы из N' = 2 V 2 ^ определителей. В этих определителях элементы

соответствующих строк заменяются отдельными слагаемыми, представленными

произведением у2" или уn и р2 n cos2 ф. n или р2 n sin2 ф ■ n, или

2Ру,и Sin ф j,n C0S ф j,n , или 2Р j1,nl Р;2,n2 Sinф;i,n1C0S фj2,n2 ' или Р j,n C0S ф j,n , или

p j n sin ф j n множителей. При этом все элементы отдельного слагаемого определителя произвольного столбца матрицы содержат общий множитель у2n или уn , который с учетом свойств определителя можно вынести за его знак. Аналогичным образом представляются определители матриц H1,H2,...,H5 в виде суммы из

10N+Í N |+N2 -13 10N + 2| N 1-12 10N+| N |+N2 -13 5N+2| N |+N2-9 5N+ 2| N |+N2-9

2 ^2 J ,2 ^2 J , 2 ^2 J , 2 ^2 J ,2 ^2 J слагаемых соответственно с вынесенными множителями, кратными уn .

Поскольку исходно кривая (6) в области Z является окружностью при у = 1 в выражении (8), то a11 = a22 и a12 = 0 , а инварианты кривой

11 = a11a22a33 + 2a12a13a23 + a22a13 + a11a23; 12 = a11a22 - a12; 13 = a11 + a22 (9)

удовлетворяют условиям 12 > 0; 11 • 13 < 0; = 412 [10].

С учетом указанного вида определителей матриц M и H1,H2,...,H5 при условии M12 = M5 2 = M13 = M5 3 = M15 = M5 5 = 0, задающих коэффициенты (7)

уравнения кривой второго порядка (6) и, как следствие, вид инвариантов (9): условия 12 > 0 и 11 • 13 < 0 для инвариантов кривой (6), останутся неизменными.

Следовательно, кривые при преобразовании S2 (Z) будут эллипсами. ■

Указанное изменение остальных круговых двуугольников определяет постановку второй модельной задачи: отображение полуплоскости с выброшенным сегментом эллипса (эллиптической луночки) на полуплоскость. Решение данной задачи базируется на следующем утверждении.

Теорема 1. Конформное отображение верхней полуплоскости с выброшенным сегментом, образованным пересечением полуплоскости с произвольной одно-связной областью © с гладкой границей Г расширенной плоскости комплексного переменного, на верхнюю полуплоскость при условии За (ф), b (ф) для Уфе[0; п] осуществляется обобщенной функцией Жуковского вида ЦС) = S3 (S2 (S1 (Z))) при

S1 (Z) = [Z-a (ф)][- b (ф)]-1; S2 (0 = Су(ф); S3 (Z) = 0,5у (ф)[а (ф) -b (ф)](1+ Z V(1 -Z)

и нормировке = та'(да) = 1, (10)

где ф = arg(Z); a(ф), b (ф) - точки пересечения действительной оси с соприкасающейся окружностью; у(ф) = л/(л^(ф)); ß^) - угол, образованный соприкасающейся окружностью в точке ее пересечения a (ф) с действительной осью. При этом соприкасающаяся окружность формируется границей выброшенного

из полуплоскости сегмента односвязной области в окрестности точки теГ при arg (т) = ф.

Доказательство. Пусть граница Г области © в верхней полуплоскости определена кривой, заданной уравнением т(ф) = u (ф) + iv (ф) с параметризацией по 0 <ф<п. Для произвольной точки кривой Г определим соприкасающуюся окружность (рис. 2). Данная окружность имеет радиус р , соответствующий радиусу кривизны кривой в окрестности т, и центр, расположенный в z0 = x0 + iy0.

Введем обозначения

ц(С) = ргга=|+ in = (Z - a)/(£- b), (11)

где S = arg[(Z-aV(Z-b)] - угол, под которым из заданной точки Z исходной верхней полуплоскости с выброшенным сегментом при arg (Z) = ф видны точки a и b; р = |(Z - a)/(Z -b)| - отношение расстояний от точки Z до точек a и b. Так как Z = x + iy, то преобразование (12) приведем к виду

in = (l-(x - b )(a - b )/[(x - b )2 + y2 ]) + i (y (a - b )/[(x - b ) + y2]). (12)

Рис. 2. Геометрическое представление задачи отображения полуплоскости с выброшенным сегментом

С учетом (12) параметры р ид могут быть определены через соответствующие выражения:

р2 = 1 -(a - b )(2x - b - a ))[(x - b )2 + y2]; tg d= y (a-b))[ x2 + y2 -x (a + b ) + ab ]. (13)

Предположим, что обобщенная функция ra(Z) Жуковского вида (10) отображает кривую т на действительную ось. Тогда для ц(х) = S2 (S1 (т)) справедливо равенство

Im (ц(т)) = 0. (14)

При подстановке (11) в (14) получим

Im(ц(т)) = ру sinуЗ = 0, или sin(лЗ/(л-р)) = 0. (15)

Условие (15) выполняется при равенстве аргумента под знаком синуса nk для k = 0,1,2,..., следовательно:

3 = л-р. (16)

Подставляя в равенство (16) определение угла $ из (13) и угла Р = arctg ((a -x0 )/| y01) (см. рис. 2), получим тождество

y0y (a - b) = (a - x0)(x2 + y2 - x(a + b) + ab). (17)

Принимая во внимание равенства:

a = x0 + h; b = x0 - h; h = -\]p2 - y^; x - x0 =p cos T; y - y0 =p sin T, (18) приведем (17) преобразованиями

y0y2h = h (x2 + y2 - 2xx0 + x^ - h2);

2y0y = (x - x0 )2 + y2 - h2; (19)

2y0у = P2 cos2 T + у2 - P2 + Уо

к виду (y - Уо )2 =p2sin2 T, (20)

которое соответствует соотношениям (18). Следовательно, исходное предположение и равенство (14) верны.

Поскольку при отображении M-(Z) точка Z = ^ переходит в точку ц(да) = 1, то для выполнения условий нормировки (10) при конформном отображении на верхнюю полуплоскость та применим дополнительное преобразование

ЦС) = k (1+ц(0)(1 -Ц(О)-1 , ч 1 + (Z-a)Y(Z-b)-Y (1 -b/Z/+(1 -a/Z)y

или та(0 = k—^--'— = k±--i-, (21)

™ 1 -(Z-a)Y(Z-b )-Y (1 - b/C)Y-(1 - a/Z/

где k - некоторая комплексная константа. Найдя производную функции ra(Z):

2k íy [a(1 - db¡5Z) -b (1 - da/dZ) -Z (da/dZ - 5b/dZ)]/Z2 дта = |+(1 - a/Z)(1 - b/Z)[ln (1 - a/ Z)- ln (1 - b/Z)]/Z

(22)

Д-У

dZ [(1 -a/Z/-(1 -b/Z)y] (1 -a¡Z)1-Y (1 -b/Z)1

путем разложения в биномиальный ряд по степеням (1 - a/Z/ , (1 - b¡Z/ знаменателя, получим для больших |Z| выражение

дта^ 2k (1 - a/Z/-1 (1 - b¡Z/-1 íy[a (1-Sb/SZ)-b(1 -Sa/dZ)-Z(da/dZ-db/dZ/+ 1 ^^^^ 5Z" y2 (a - b)2 l+Z2 (1 - a/Z)(1 - b/Z/ln (1-a/ Z)- ln (1-^//Т ( )

Применяя в (23) к da/dZ и db¡dZ правило дифференцирования параметрически заданных функций и учитывая, что

d<p/5Z = - Im (Z) IZI 2 = -ri sin <p(cos ф + i sin ф)/г 2 = - sin ф(ш ф-i cos ф)/г ,

при ^ ^го имеем

та'(го) = 2к (у (а - Ь ))-1 или к = 0,5у(а - Ь). (24)

Подставляя теперь константу к в (21), получим искомое отображение (10). ■ С учетом теоремы 1 определим зависимость переменных у(ф), а (ф),Ь (ф) в выражении (10) от угла ф в полуплоскости с эллиптической луночкой, заданной эллипсом с полуосями а0 и Ь0, центром в 20 = х0 + гу0 и углом а наклона к действительной оси. При этом зависимость угла р = агссоБ(|у0|/р) для эллиптической луночки от ф определяется с помощью соотношений

y0 = y-psinТ; p = (sin2a+^cos2a) /a0b0; T = ±2arctg (yjk° +1 + \)jl

>0;

k = aO cos2 ф + bosin2 ф

kH = -

cosф sinф( - b02) - (u - x0 )a0b0 ( cos2 ф + bo sin2 ф -(u - x0 )2)

1Ф/ , • ч f arccos (a,) if |ф- ф| < 0,5п;

t = u + iv = е1ф (a0cos a + b0 sina) + z0; a = <! ; , (25)

[arccos (a 2) otherwise,

a, 2 = (a0 (sinф - tgфcosф)0 ±|b0 (cosф + tgфsinф)_0 - G2 ;

G = (x0 tgф-y0); D = a° (sinф-tgфcosф)2 + b^ (cosф + tgфsinф)2 .

Параметры отображения (10) a (ф), b (ф) находятся по соответствующим тождествам (18) с использованием соотношений (25).

Построение обратного отображения с применением рассматриваемых модельных задач осуществляется с использованием леммы 1 для случая, когда все вершины исходной многоугольной области Q располагаются на границе окружности Q'. В случае же, когда вокруг исходной многоугольной области Q невозможно описать такой окружности, возникает ситуация, приводящая к необходимости отображения полуплоскости с выброшенным сегментом с кусочно-гладкой границей Г (см., например рис. 1, граница между точками a2, a4) на полуплоскость Im Z > 0 . В этом случае отображение выполняется с применением обобщенной функции Жуковского вида (10) для каждого гладкого участка границы Г - отдельного эллипса или круга. При этом добавление в полуплоскость кругового или эллиптического сегмента согласно лемме 1 не означает, что все границы оставшихся луночек будут являться границами окружностей или эллипсов. В такой ситуации возможен случай формирования выброшенного сегмента, образованного пересечением верхней полуплоскости с гиперболической луночкой.

С учетом теоремы 1 определим решение третьей модельной задачи. Она заключается в нахождении отображения полуплоскости с гиперболической луночкой на полуплоскость обобщенной функцией Жуковского вида (10). Граница луночки задана ветвью гиперболы с полуосями a0 и b0, центром в z0 = x0 + iy0 и углом наклона a к действительной оси.

Зависимость угла р = arccos (|y0|/p) для гиперболической луночки от ф и параметров a (ф), b (ф) задается с помощью соотношений

(22 22 \3/2 a0 sh а + b0 ch а)

= у-рsin т, р= , ^ , ^ 2-^ |Ч ; T = ±2arctg

(a0b0 |sh2 а-ch2 а|)

VkT+7+

>0;

L =

al cos2 ф- b° sin2 ф

cos ^ sin ^ ( + b0) - (u - x0) a2b0 (b^ (u - x0 )2 - cos2 ^ + b0 sin2

Jarcsh (aj) if + 0,5n| ^ 9; [arcsh (a 2) otherwise,

(b0 (cos ^ + tg 9 sin ^)G ± |a0 (tg 9 cos sin ^)|Vd + G2 D;

t = u + iv = e^(a0ch a + b0sh a) + z0; а =

-1/2 '

(26)

а1,2 =

G = (x0 tg ф - у0); D = a° (tg ф cos ф - sin ф)2 - b° (cos ф + tg ф sin ф)2.

Параметры кривых второго порядка, соответствующих при обратном отображении ребрам en ={ An, An+1} исходной многоугольной области Q , могут быть

получены по заданным пяти точкам hj = Tje1ф/ (j = 1,5) с учетом выражения (7).

Примеры обратного отображения наперед заданных многоугольных областей представлены на рис. 3.

Im t

т™ -

а21

Re t

Re t

К

Рис. 3. Примеры обратного отображения наперед заданной четырехугольной и семиугольной областей на единичный круг

Заключение

Таким образом, полученные результаты позволяют за конечное число конформных отображений, равных числу вершин N многоугольника, выполнить конформное отображение заданной односвязной многоугольной области на единичный круг. В случае, если сегмент, формирующий луночку в верхней полуплоскости, на границе содержит точки, для которых не существует a (ф), b (ф), то добавление выброшенного сегмента в полуплоскость выполняется последовательностью конформных отображений, удовлетворяющих условиям теоремы 1. При алгоритмической реализации сформированного метода на этапах добавления n -го сегмента в верхнюю полуплоскость, соответствующего n -му ребру en ={An, An+1} исходного многоугольника, рекомендуется предварительно выполнять поворот многоугольной области в расширенной комплексной плоскости на угол у таким образом, чтобы точка, являющаяся серединой ребра en , располагалась на отрицательной части действительной оси. В целом, полученные результаты расширяют применимость обобщенной функции Жуковского, заданной в виде (10).

ЛИТЕРАТУРА

1. Фильчаков П. Ф. Приближенные методы конформных отображений. Справочное руководство. Киев: Наукова думка, 1964. 536 с.

2. ЛаврентьевМ.А., Шабат Б.В. Методы теории функций комплексного переменного. 4-е изд., испр. М.: Наука, 1973. 736 с.

3. Иванов В.И., Попов Ю.В. Конформные отображения и их приложения. М.: Едиториал УРСС, 2002. 374 с.

4. Driscoll T.A., Trefethen L.N. Schwarz - Christoffel mapping. Cambridge, UK: Cambridge University Press, 2002. 132 p.

5. Канторович Л.В. Эффективные методы в теории конформных отображений // Изв. АН СССР. Сер. матем. 1937. Т. 1. Вып. 1. С. 79-90.

6. Архипов Н.С., Полянский И.С., Степанов Д.Е. Барицентрический метод в задачах анализа поля в регулярном волноводе с произвольным поперечным сечением // Антенны. 2015. № 1(212). С. 32-40.

7. Полянский И.С. Барицентрические координаты Пуассона для многомерной аппроксимации скалярного потенциала внутри произвольной области (Часть 1) // Вестник СГТУ. 2015. № 1(78). С. 30-36.

8. Полянский И.С. Барицентрические координаты Пуассона для многомерной аппроксимации скалярного потенциала внутри произвольной области (Часть 2) // Вестник СГТУ. 2015. № 1(78). С. 36-42.

9. СмирновВ.И. Курс высшей математики. Т. 3. Ч. 2. Изд. 9. М.: Наука, 1974. 672 с. 10. Голованов Н.Н. Геометрическое моделирование. М.: ИФМЛ, 2002. 472 с.

Статья поступила 08.10.2015 г.

Radygin V.M., Polyanskii I.S. MODIFIED METHOD OF SUCCESSIVE CONFORMAL MAPPINGS OF POLYGONAL DOMAINS

DOI 10.17223/19988621/39/3

The method for solving the inverse Schwarz - Christoffel problem - conformal mapping of a given polygonal area on the canonical domain, the unit circle - is developed in the paper. The method is based on the use of a sequence of conformal mappings related to the mapping of the polygon onto the upper half with discarded segments by a Cayley linear fractional transformation

followed by sequential addition of the discarded segments to the upper half-plane. Model problems consider conformal mappings of circular, elliptical, and hyperbolic lunes to the upper halfplane. The solution of the presented model problems is based on the generalized Zhukovsky function; the obtained results expand the applicability of the Zhukovsky function as compared to existing methods of its application in the implementation of conformal mappings. The working ability of the solution was tested by specific examples in solving the problem of conformal mapping of quadrilateral and heptagonal domains to the unit circle. Recommendations for the algorithmic implementation of the method are presented.

Keywords: conformal mapping, polygonal domain, unit circle, inverse Schwarz-Christoffel problem.

RADYGIN Vladimir Mikhailovich (Doctor of Physics and Mathematics, Prof., Academy of the Federal Security Service Russia, Orel, Russian Federation) E-mail: [email protected]

POLANSKY Ivan Sergeyevich (Candidate of Technical Sciences, Academy of the Federal Security Service Russia, Orel, Russian Federation) E-mail: [email protected]

REFERENCES

1. Fil'chakov P.F. Priblizhennye metody konformnykh otobrazheniy. Spravochnoe rukovodstvo. Kiev, Naukova dumka Publ., 1964. 536 p. (in Russian)

2. Lavrent'ev M.A., Shabat B.V. Metody teorii funktsiy kompleksnogo peremennogo. Moscow, Nauka Publ., 1973. 736 p. (in Russian)

3. Ivanov V.I., Popov Yu.V. Konformnye otobrazheniya i ikh prilozheniya. Moscow, Editorial URSS Publ., 2002. 374 p. (in Russian)

4. Driscoll T.A., Trefethen L.N. Schwarz - Christoffel mapping. Cambridge, UK, Cambridge University Press, 2002. 132 p.

5. Kantorovich L.V. Effektivnye metody v teorii konformnykh otobrazheniy. Izv. ANSSSR. Ser. matem., 1937, vol. 1, no. 1, pp. 79-90. (in Russian)

6. Arkhipov N.S., Polyanskiy I.S., Stepanov D.E. Baritsentricheskiy metod v zadachakh analiza polya v regulyarnom volnovode s proizvol'nym poperechnym secheniem. Antenny, 2015, no. 1(212), pp. 32-40. (in Russian)

7. Polyanskiy I.S. Baritsentricheskie koordinaty Puassona dlya mnogomernoy approksimatsii skalyarnogo potentsiala vnutri proizvol'noy oblasti (Chast' 1). Vestnik SGTU, 2015, no. 1(78), pp. 30-36. (in Russian)

8. Polyanskiy I.S. Baritsentricheskie koordinaty Puassona dlya mnogomernoy approksimatsii skalyarnogo potentsiala vnutri proizvol'noy oblasti (Chast' 2). Vestnik SGTU, 2015, no. 1(78), pp. 36-42. (in Russian)

9. Smirnov V.I. Kurs vysshey matematiki, vol. 3, part. 2. Moscow, Nauka Publ., 1974. 672 p. (in Russian)

10. Golovanov N.N. Geometricheskoe modelirovanie. Moscow, IFML Publ., 2002. 472 p. (in Russian)