Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. Т. 62, № 5 (2019), с. 413-421 Energetika. Proc. CIS Higher Educ. Inst. and Power Eng. Assoc. V. 62, No 5 (2019), pp. 413-421 413
https://doi.org/10.21122/1029-7448-2019-62-5-413-421 УДК 517.958:519.6
Применение полилогарифмов к приближенному решению неоднородного телеграфного уравнения для линии без искажений
П. Г. Ласый1*, И. Н. Мелешко1*
'^Белорусский национальный технический университет (Минск, Республика Беларусь)
© Белорусский национальный технический университет, 2019 Belarasian National Technical University, 2019
Реферат. В статье рассматривается смешанная задача для хорошо известного в электротехнике и электронике телеграфного уравнения при условии, что линия свободна от искажений. Эта задача сводится к аналогичной для одномерного неоднородного волнового уравнения. Ее решение можно найти как сумму решения смешанной задачи с однородными краевыми условиями для соответствующего однородного волнового уравнения и решения неоднородного волнового уравнения с однородными краевыми и нулевыми начальными условиями. Решения обеих задач можно отыскать методом разделения переменных в виде ряда по тригонометрическим функциям точки линии с коэффициентами, зависящими от времени. Такие решения неудобны для реального применения, поскольку требуют вычисления большого числа интегралов и трудно оценить погрешность их вычислений. Предлагается альтернативный способ решения этой задачи, основанный на использовании специальных функций - полилогарифмов, которые представляют собой комплексные степенные ряды со степенными же коэффициентами, сходящиеся в единичном круге. Точное решение задачи выражается в интегральной форме через мнимую часть полилогарифма первого порядка на единичной окружности, а приближенное - в виде конечной суммы через действительную часть дилогарифма и мнимую часть полилогарифма третьего порядка. Все указанные части полилогарифмов являются периодическими функциями, имеющими полиномиальные выражения соответствующих степеней на отрезке длиной в период. Это позволяет эффективно находить приближенное решение задачи. Также найдена простая и удобная оценка погрешности приближенного решения задачи. Она линейна относительно шага разбиения линии и шага разбиения временного диапазона, на котором рассматривается задача. Оценка является равномерной по длине линии в каждый фиксированный момент времени. Приведен конкретный пример решения задачи разработанным способом, построены графики точного и приближенного решений.
Ключевые слова: телеграфное уравнение, волновое уравнение, смешанная задача, приближенное решение, оценка погрешности, полилогарифм
Для цитирования: Ласый, П. Г. Применение полилогарифмов к приближенному решению неоднородного телеграфного уравнения для линии без искажений / П. Г. Ласый, И. Н. Ме-лешко // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2019. Т. 62. № 5. С. 413-421. https://doi.org/10.21122/1029-7448-2019-62-5-413-421
Адрес для переписки
Ласый Петр Григорьевич Белорусский национальный технический ул. Я. Коласа, 12,
220013, г. Минск, Республика Беларусь Тел.: +375 17 292-82-73 [email protected]
Address for correspondence
Lasy Petr G.
Belarusian National Technical University 12 Ya. Kolasa str., 220013, Minsk, Republic of Belarus Tel.: +375 17 292-82-73 [email protected]
Application of Polylogarithms to the Approximate Solution of the Inhomogeneous Telegraph Equation for the Distortionless Line
P. G. Lasy1), I. N. Meleshko1
'-Belarusian National Technical University (Minsk, Republic of Belarus)
Abstract. The paper deals with a mixed problem for the telegraph equation well-known in electrical engineering and electronics, provided that the line is free from distortion. This problem is reduced to the analogous one for the one-dimensional inhomogeneous wave equation. Its solution can be found as the sum of the solution for a mixed homogeneous boundary value problem for the corresponding homogeneous wave equation and for the solution of a non-homogeneous wave equation with homogeneous boundary data and zero initial conditions. Solutions to both problems can be found by separating the variables in the form of a series of trigonometric functions of the line point with time-dependent coefficients. Such solutions are inconvenient for real application because they require calculation of a large number of integrals, and it is difficult to estimate the miscalculation. An alternative method for solving this problem is proposed, based on the use of special functions, viz. polylogarithms, which are complex power-series with power coefficients converging in a unit circle. The exact solution of the problem is expressed in the integral form via the imaginary part of the first-order polylogarithm on the unit circle, and the approximate one is expressed in the form of a finite sum via the real part of the dilogarithm and the imaginary part of the third-order polylogarithm. All these parts of the polylogarithms are periodic functions that have polynomial expressions of the corresponding powers on the segment of the length equal to the period. This makes it possible to effectively find an approximate solution to the problem. Also, a simple and convenient error estimate of the approximate solution of the problem is found. It is linear with respect to the step of splitting the line and the step of splitting the time range in which the problem is considered. The score is uniform along the length of the line at each fixed point of time. A concrete example of solving the problem according to the proposed mode is presented; graphs of exact and approximate solutions are constructed.
Keywords: telegraph equation, wave equation, mixed problem, approximate solution, error estimation, polylogarithm
For citation: Lasy P. G., Meleshko I. N. (2019) Application of Polylogarithms to the Approximate Solution of the Inhomogeneous Telegraph Equation for the Distortionless Line. Energetika. Proa CIS Higher Educ. Inst. and Power Eng. Assoc. 62 (5), 413-421. https://doi.org/10. 21122/10297448-2019-62-5-413-421 (in Russian)
Введение
Электромагнитное поле, возникающее вокруг проводника при прохождении по нему электрического тока, вызывает колебания напряжения и силы тока.
Пусть ось Ox совпадает с осью прямолинейного проводника длиной l, а один из его концов примем за начало отсчета. Известно [1-6], что величины разности потенциалов и силы тока в любой точке х е [0, l] проводника в любой момент времени t > 0 являются решениями телеграфного уравнения
д xx w = LCd tt w + (RC + GL)d t w + GRw, (1)
где w = w(x, t) - неизвестная функция (сила тока или разность потенциалов); L, C, R, G - величина самоиндукции, емкости, сопротивления и проводимости изоляции соответственно, рассчитанные на единицу длины линии.
Уравнение (1) имеет многочисленные приложения в электротехнике и электронике, с которыми можно ознакомиться, например, в [1, 7-9]. Для того чтобы решение (1) существовало и было единственным, должны быть заданы начальные и краевые условия. Такая задача называется смешанной.
Если рассматриваемая линия не имеет искажений, т. е.
R G
— = — = о, L C
то телеграфное уравнение заменой искомой функции ^ = ехр(-с^) ы приводится к одномерному однородному волновому уравнению
дии = а 2дxxu, (2)
1
где а = ,- - скорость распространения волн напряжения и силы тока,
^JLC
возникающих в линии. Основная часть
Решение смешанной задачи с неоднородными, т. е. ненулевыми, краевыми условиями для (2) сводится к аналогичной задаче с однородными краевыми условиями для неоднородного волнового уравнения. Эту последнюю задачу мы и рассмотрим.
Требуется найти решение уравнения
д ии = а 2д хи + д( х, t) (3)
в прямоугольнике П = {(х^)| х е [0,1], t е [0,Т]}, Т > 0 при заданных начальных
и(х, 0) = /(х), ди(х,0) = Е(х), х е [0,1] (4)
и краевых условиях первого рода
и(0,0 = Ф0(0, и(I,0 = ф(t), t е [0,Т]. (5)
Наложим ограничения на известные функции д(х, t), /(х), Е(х), ф0(t), ф;(0. Функцию д(х, t) будем считать удовлетворяющей условию Липшица в прямоугольнике П, т. е. существует положительная постоянная L > 0 такая, что для любых двух точек данного прямоугольника выполняется неравенство
| д(х1, О - Я(х2, t2)|< L(| х1 - х2 | + | tl - t2 |). Коротко этот факт будем записывать следующим образом:
д(х,t) е Lip(L,П). Далее функция /х) дифференцируема на отрезке [0, /] f'(х) е Ь1р(А,[0,I]); F(х) е Llp(L2,[0,I]).
Функции ф0 (t) и ф1 (t) дважды дифференцируемы на отрезке [0, T]
Ф0(t) е Lip(Lo,[0,T]); ф,'(t) е Lip(L,,[0,T]).
Проведем в уравнении (3) подстановку
u = V + ф(x, t) (6)
с функцией ф( x, t) = ф0 (t) +1-1 x(фl (t) - ф0 (t)).
В результате получим также неоднородное волновое уравнение
д ttV = а 2д xxV + Ч\( x, t), (7)
где qx (x, t) = q(x, t) - dttф(x, t).
Несложно проверить, что qx (x, t) е Lip(L3, П), причем L3 = L + 2L0 + L, + +1-lM, M = max | ф0 (t) - ф'(t) |.
te[0,T ]
Начальные условия принимают вид:
v(x,0) = f (x); дt (x, 0) = F (x) (8)
с функциями f (x) = f (x) - ф(x, 0), f(x) eLip(L1 ,[0, l]), F1 (x) = F(x) -
- д,ф(x, 0) е Lip(L2 ,[0, l]), где L2 = L2 +1| ф0 (0) - ф, (0)|.
Краевые условия, ввиду подстановки (6), становятся однородными:
v(0, t) = 0; v(l, t) = 0; t е [0, T]. (9)
Решение смешанной задачи (7)-(9) можно найти, следуя [2]:
v( x, t) = V (x, t) + v*( x, t), (10)
где V (x, t) - решение смешанных задач (8), (9) для однородного волнового уравнения
д^ = а 2д xxv; (11)
V* ( x, t) - решение неоднородного уравнения (7) при нулевых начальных и краевых условиях.
В [10], используя полилогарифмы [11, 12], точное решение задач (8), (9) для уравнения (11) представлено в виде
1 l 2
v(x, t) = —JX (f(y) + (-1)k+1 a-lFx(y))Pk (x, t, y )dy, (12)
2П 0 к=1
где при к = 1 или к = 2 Pk (x, t, y) = (-1)k+1(N'(ю(at + (-1)k+1 x + (-1)k y)) -
- N1 (ro(at + (-1)k x + (-1)k y))); ю = у; N1(x) - мнимая часть полилогарифма первого порядка на единичной окружности
W • 1
N (x) = ImL (exp(ix)) = X, x е R. (13)
k=1 k
Там же найдено и удобное для вычислений приближенное решение
1 п 2 х уп(х,0 = — ££(Дхк_0,5) + (-1)5+1 а"1 F (х^))(х,t,у) * ,
2п k=1 s=1
(14)
где при фиксированном натуральном n и s = 1 или s = 2 x0 = 0; xk = kh;
h = l /n; xk-05 = xk-1 + 0,5h; k = 1,n; gs (x, t, y) = M2 (<(at + (-1)s+1 x + (-1)s y)) -
- M2 (< at + (-1)s x + (-1)sy)); M2(x) - действительная часть дилогарифма на единичной окружности
M2 (x) = Re L (exp(zx)) = ¿ , x e R.
k=1 k
Абсолютная величина погрешности вычисления решения (12) по формуле (14) оценивается величиной
^К+М^ (15)
2a
То есть погрешность равномерна по x e [0, l], t > 0 и имеет первый порядок малости относительно шага h разбиения отрезка [0, I].
Выразим теперь решение v*(x, t) уравнения (7) через полилогарифмы. В [2] оно записано в виде
ад
v* (x, t) = ^Tk (t )sin k <x, (16)
k=1
где при любом натуральном k
Tk (t) = —[[ q1 (y, s)sin k<a(t - s) • sin k<ydyds (17)
k na U
и интегрирование ведется по прямоугольнику Пlt = { (y, s) | y e [0, l], S e [0,t]}
Подставив выражение (17) в (16), выполнив несложные преобразования и учитывая (13), получим:
V* (x, t) = 2na И q\( y, s) P( x, t, y, s)dyds, (18)
где Р(х,^у,5) = £ (-1)*- 5) + (-1)[(к-1)/2] х + (-1)[к/2]у)); через [•]
к=1
обозначена целая часть действительного числа.
Отыщем приближенное выражение для решения у*( х, t). Набросим на прямоугольник Пы сетку с узлами в точках (х^, tк ), где х^ = кхк*,
к1 = 0,п1, tk = к2т*, к2 = 0,п; к* = 1 / п{, т* = t / п{ - шаги сетки по переменным х и t соответственно. Далее заменим в (18) под знаком интеграла
it
в кажд°м го n¡nt прям°уг°льник°в П^ = {(y,s) | xki- <y < xki, tki- < s < X^ } функцию q1(x,t) ее значением в средней точке (x^-05, t^-05), где
Xk1 -0,5 = Xk1 -1 + 0,5 h , k1 = 1, П1, tk2-0,5 = tk2 -1 + 0,5 , k2 = 1, П .
В результате получим
1 n n
V (x, t) « Vn¡n¡ (x, t) = — XX^ (xk1 -0,5 , tk2-0,5 ) Ц P(x t, У, s)dyds. (19) 2n°k1 = k2 = Ц^
Поскольку первообразными для функций Nx) и M 2( x) являются соответственно функции -M2(x) и N3(x), причем N(x) = ImZ3(exp(/x) =
■A sin kx ,
= > —з—, x e R - мнимая часть полилогарифма третьего порядка на еди-
k=1 k
ничной окружности, то
ff P(x,t,y,s)dyds = Qk,k2(x,t), ara2 1 2
■"•¡¡1-1
'k2-1
где аккг (х,0 = X (-1)к+[к/2]N4.ш(а(г - 5) + (-1)[(к-1)/2] х + (-1)[к/2]у))[
к=1
Следовательно, приближенным значением для функции V* (х, t) служит выражение V* (х, t), вычисляемое по формуле
1 "' п
^ЩЩ ( х t) = --ГТ XX 91 ( Хк1-0,5 , tkг-0,5 )бк1кг ( x, t). (г0)
2па ® к,=1 кг =1
Найдем оценку абсолютной погрешности вычисления решения (18) по формуле (20). Из (18) и (19) следует, что
1 "'"'¡г
V (х I) - v*'"'(х, I) = -— XX Д (9 (У, 5) - 91(хк!-0,5 , (к2-0,5 ))Р(х(, y, 5
гп°к1 =1 кг=1 Пк!кг
Нечетная периодическая функция Nх) на промежутке (0, л] имеет выражение х) = 0,5(п- х) и, значит, | N1(х)|< 0,5п для всех х е R. Следовательно, при всех х, t, у, 5
| Р(х,Г ,у,5)|< 2л.
Учитывая это неравенство и тот факт, что 91(х, t) е Ыр(Х3, П), получаем:
IV* (х, ^ - (х, 0|< -Л- X X Ц 1?1(У, 5) -гп°к1=1 кг =1 Пкк,
Ч\(хк1 -0,5>^2-0,5) II Р(х, I,y,s) | dyds <
< 2"" 11 13(1 У - Хк1 -0,5 1 + 1 • - (к2-0,5 |) 1 2п ^ <
2П°к1 =1 к2 =1 Пк1к2
Т "' "' .. т и
< -22 22 Л +°'5т*)dyds = ТР(^ + т*).
0 к1 =1 к2 =1 Пк!к2 20
Таким образом, в прямоугольнике П имеет место следующая оценка погрешности вычисления решения х,') уравнения (7):
IV* (х,')-у;„' (х, ')|< М (к* + т*). (21)
' 2о
Проведем теперь сборку решения поставленной задачи (3)-(5). Принимая во внимание (6) и (1°), точное решение находится по формуле
и (х,') = V (х,') + V* (х,') + ф( х,'), (22)
где V(х,'), v*(х,') - находятся по (12) и (18), а приближенное - по формуле
иппл (х') = ^(х,') + ^(х') + Ф(х'X (23)
в которой функции уп (х,') и V* (х,') имеют, соответственно, выражения (14) и (2°).
Из (15) и (21) следует общая оценка погрешности решения задачи
u(x,t)-ыппл (x,t) |< 2-((aL +L2)h + L3t(h* + т*)).
Проведенные исследования позволяют сформулировать следующее утверждение.
Теорема. При указанных выше предположениях относительно функций д(х,'), /(х), ^(х), ф0('), ф' (') точное решение смешанной краевой задачи (4), (5) для неоднородного волнового уравнения (3) находится с помощью полилогарифма первого порядка по формуле (22), а приближенное -с помощью дилогарифма и полилогарифма третьего порядка по (23). Абсолютная погрешность вычисления оценивается величиной
((аЦ +ТТ2 ^ +Тъг (И* + т*)).
То есть она линейна по шагам к и к* разбиения отрезка [0,'] и шагу т* разбиения отрезка [0,']. Кроме того, эта оценка является равномерной по х е [0,'] при каждом фиксированном ' е [0, Т].
Следует заметить, что, несмотря на кажущуюся громоздкость формул (14), (2°), (23), они достаточно эффективны, так как использующиеся
в вычислениях функции М2 (х) и N3 (х) являются элементарными, а именно:
М2 (х) = (П-4х 1)2 - П-; N3 (х) = 12 х(п-1 Х |)(2п-1 х I), х е [-п, п]. (24)
Для вычисления значений этих функций в произвольной точке х е R следует, учитывая их периодичность с периодом 2п, заменить в формулах (24) х на а(х), где
х- | х |, х Ф 0,
а( x) = sign( x)
Г x Л
|x|-2п
V _ 2п_ /
sign( x) = -
I 0, x = 0.
Пример. Найти решение неоднородного волнового уравнения
дии = 9дхи + 40л(л(36х2 - 1^п2л(х2 + ^-9^2л(х2 + 0)
при начальных и(х,0) = х2-—х + — + 108т2лх2, д,и(х,0) = 2х-1 + 20лсо82лх2,
2 16
9 1 2 9
х е [0, 1] и краевых и(0,¿) = 9г - Г + — + 10зт2л^, и(1,Г) = 9г + t + — +
16 16
+ 108т2л^ + 1), t > 0 условиях.
Точным решением этой задачи является функция
и(х,0 = х2 + 9Г2 + 2xt -1 х -1 + — + 10эт2л;(х2 +1).
2 16
Вычисление приближенного решения задачи по формуле (23) в квадрате {(х, t)| х е [0,1], t е [0,1]} при п = п1 = п = 100 дает максимальную погрешность, не превышающую 0,026, что является вполне приемлемым для данной сетки с суммой шагов Н +Н* +т* = 0,03. Графики точного и (х, ^ и приближенного и(х.1) = и]пп |пп |пп(х./) решений данной задачи представлены на рис. 1.
а Ь
Рис. 1. Графики решений: a - точного u(x, t); b - приближенного u(x, t) = u100 100 100(x, t)
Fig. 1. The graphics of the solutions: a - of exact one u(x, t); b - of approximate one u(x, t) = u100 100 100(x, t)
ВЫВОД
С использованием полилогарифмов первого-третьего порядков найдены точное и приближенное решения телеграфного уравнения для линии без искажений при произвольных начальных и граничных условиях. Получена равномерная по длине линии оценка погрешности приближенного решения.
ЛИТЕРАТУРА
1. Анго, А. Математика для электро- и радиоинженеров / А. Анго. М.: Наука, 1964. 772 с.
2. Кошляков, Н. С. Дифференциальные уравнения математической физики / Н. С. Кошля-ков, Э. Б. Глинер, М. М. Смирнов. М.: ГИФМЛ, 1962. 767 с.
3. Араманович, И. Г. Уравнения математической физики / И. Г. Араманович, В. И. Левин. М.: Наука, 1969. 288 с.
4. Смирнов, В. И. Курс высшей математики: в 5 т. / В. И. Смирнов. М.: Наука, 1974. Т. 2. 479 с.
5. Мышкис, А. Д. Лекции по высшей математике / А. Д. Мышкис. СПб.: Лань, 2007. 688 с.
6. Остапенко, В. Телеграфное уравнение. Краевые задачи / В. Остапенко. Саарбрюккен: LAP Lambert Academic Publishing, 2012. 272 с.
7. Новиков, Ю. Н. Электротехника и электроника. Теория цепей и сигналов, методы анализа / Ю. Н. Новиков. СПб.: Питер, 2005. 384 с.
8. Бычков, Ю. А. Основы теории электрических цепей / Ю. А. Бычков, В. М. Золотницкий, Э. П. Чернышев. СПб.: Лань, 2002. 464 с.
9. Дубнищев, Ю. Н. Колебания и волны / Ю. Н. Дубнищев. СПб.: Лань, 2011. 384 с.
10. Ласый, П. Г. Приближенное решение одной задачи об электрических колебаниях в проводах с помощью полилогарифмов / П. Г. Ласый, И. Н. Мелешко // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2017. Т. 60, № 4. С. 334-340. https://doi.org/10.21122/1029-7448-2017-60-4-334-340.
11. Пыхтеев, Г. Н. Полилогарифмы, их свойства и методы вычисления / Г. Н. Пыхтеев, И. Н. Мелешко. Минск: Изд-во БГУ, 1976. 68 с.
12. Мелешко, И. Н. Специальные формулы для интегралов типа Коши и их приложения / И. Н. Мелешко. Минск: ВУЗ-ЮНИТИ, 1999. 197 с.
Поступила 26.03.2019 Подписана в печать 04.06.2019 Опубликована онлайн 30.09.2019
REFERENCES
1. Ango А. (1964) Mathematics for Electrical and Radio Engineers. Moscow, Nauka Publ. 772 (in Russian).
2. Koshlyakov N. S., Gliner E. B., Smirnov M. M. (1962) Differential Equations of Mathematical Physics. Moscow, State Publishing House of Physical and Mathematical Literature. 767 (in Russian).
3. Aramanovich I. G., Levin V. I. (1969) Equations of Mathematical Physics. Moscow, Nauka Publ. 288 (in Russian).
4. Smirnov V. I. (1974) Learning Course of Higher Mathematics. Vol. 2. Moscow, Nauka Publ. 479 (in Russian).
5. Myshkis A. D. (2007) Lectures on Higher Mathematics. Saint-Petersburg, Lan' Publ. 688 (in Russian).
6. Ostapenko V. (2012) Telegraph Equation. Boundary Value Problems. Saarbrücken, Lambert Academic Publishing. 272 (in Russian).
7. Novikov Yu. N. (2005) Electrical Engineering and Electronics. Theory of Circuits and Signals, Methods of Analysis. Saint-Petersburg, Piter Publ. 384 (in Russian).
8. Bychkov Yu., A., Zolotnitskii V. M., Chernyshev E. P. (2002) Fundamentals of the Theory of Electrical Circuits. Saint-Petersburg, Lan' Publ. 464 (in Russian).
9. Dubnishchev Yu. N. (2011) Oscillations and Waves. Saint-Petersburg, Lan' Publ. 384 (in Russian).
10. Lasy P. G., Meleshko I. N. (2017) Approximate Solution of One Problem on Electrical Oscillations in Wires with the Use of Polylogarithms. Energetika. Izvestiya Vysshikh Uchebnykh Zavedenii i Energeticheskikh Obedinenii SNG = Energetika. Proceedings of CIS Higher Education Institutions and Power Engineering Associations, 60 (4), 334-340 (in Russian). https://doi.org/10.21122/1029-7448-2017-60-4-334-340.
11. Pykhteev G. N., Meleshko I. N. (1976) Polylogarithms, their Properties and Methods of their Computation. Minsk, BSU Publ. 68 (in Russian).
12. Meleshko I. N. (1999) Special Formulas for Cauchy-Type Integrals and their Applications. Minsk, VUZ-YuNITI Publ. 197 (in Russian).
Recеived: 26 March 2019
Accepted: 4 June 2019
Published online: 30 September 2019