DOI: 10.21122/1029-7448-2017-60-4-334-340 УДК 517.958:519.6
Приближенное решение одной задачи об электрических колебаниях в проводах с помощью полилогарифмов
П. Г. Ласый1), И. Н. Мелешко1
'-Белорусский национальный технический университет (Минск, Республика Беларусь)
© Белорусский национальный технический университет, 2017 Belarusian National Technical University, 2017
Реферат. В статье рассматривается смешанная задача с однородными краевыми условиями для одномерного однородного волнового уравнения. Такая задача может возникнуть, например, при изучении колебаний силы тока и напряжения в проводнике, по которому проходит электрический ток, и линия свободна от искажения. Решение можно найти методом Фурье в виде тригонометрического ряда. Данное представление имеет только теоретический интерес, поскольку для реального вычисления необходимо, во-первых, находить большое число коэффициентов-интегралов, что само по себе - задача нетривиальная и, во-вторых, практически невозможно провести оценку погрешности вычислений. Предлагается альтернативный способ решения этой задачи, основанный на использовании трансцендентных функций - полилогарифмов, которые представляют собой комплексные степенные ряды специального вида. Точное решение задачи выражается через мнимую часть полилогарифма первого порядка на единичной окружности, а приближенное - через действительную часть дилогарифма. Кроме того, если начальные условия в задаче являются элементарными функциями, то и решение также осуществляется через элементарные функции. Найдена простая и вместе с тем эффективная оценка погрешности приближенного решения задачи. Она не зависит от времени и имеет первый порядок точности относительно шага разбиения отрезка числовой оси, на котором рассматривается задача. Указанная оценка является равномерной относительно переменных задачи - как пространственной, так и временной.
Ключевые слова: электрические колебания, решение задачи, трансцендентные функции, оценка погрешности, полилогарифм
Для цитирования: Ласый, П. Г. Приближенное решение одной задачи об электрических колебаниях в проводах с помощью полилогарифмов / П. Г. Ласый, И. Н. Мелешко // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2017. Т. 60. № 4. С. 334-340. DOI: 10.21122/1029-7448-2017-60-4-334-340
Approximate Solution of One Problem on Electrical Oscillations in Wires with the Use of Polylogarithms
P. G. Lasy1), I. N. Meleshko1
1-Belarusian National Technical University (Minsk, Republic of Belarus)
Abstract. The article considers a mixed problem with homogeneous boundary conditions for one-dimensional homogeneous wave equation. Such a problem can arise, for example, when studying
Адрес для переписки Address for correspondence
Ласый Петр Григорьевич Lasy Petr G.
Белорусский национальный технический университет Bekrusian National Technical University
ул. Я. Коласа, 12, 12 Ya. Кд^а str.,
220013, г. Минск, Республика Беларусь 220013, Minsk, Republic of Belarus
Тел.: +375 17 292-82-73 Tel.: +375 17 292-82-73
[email protected] [email protected]
oscillations of current and voltage in the conductor through which electric current flows, while the line is free from distortion. The solution can be found with the use of the Fourier method in the form of trigonometric series. This representation is of purely theoretical interest, because the real calculation should be, first, to find a large number of coefficients of the integrals, which in itself is not a trivial task and, second, it is almost impossible to assess the error of the calculations. An alternative way of solving this problem based on the use of transcendental functions i. e. poly-logarithms that represent complex power series of a special kind. The exact solution of the problem is expressed through the imaginary part of a polylogarithm of the first order on the single circle and the approximate one - via the real part of the dilogarithm. In addition, if the initial conditions in the problem are elementary functions, then the solution is also computed using elementary functions. A simple and effective error estimate of the approximate solution has been found. It does not depend on time and it has the first-order of accuracy regarding the step of a partitioning segment of the numerical axis on which the problem is considered. This valuation is uniform with respect to the variables of the problem - both spatial and temporal.
Keywords: electric oscillations, solution of the problem, transcendental functions, error estimate, polylogarithm
For citation: Lasy P. G., Meleshko I. N. (2017) Approximate Solution of One Problem on Electrical Oscillations in Wires with the Use of Polylogarithms. Energetika. Proc. CIS Higher Educ. Inst. and Power Eng. Assoc. 60 (4), 334-340. DOI: 10.21122/1029-7448-2017-60-4-334-340 (in Russian)
Колебания силы тока и напряжения при прохождении по проводу длиной l > 0 электрического тока удовлетворяют телеграфному уравнению
dxxw = a0dttw + 2b0dtw + c0W, (1)
где w = w(x, t) - неизвестная функция переменных 0 < x < l и t > 0; t > 0, a0, b0, c0 - положительные постоянные [1-5].
Другие приложения уравнения (1) можно найти, например, в [6-8]. После замены в (1) искомой функции по формуле
-b0 t an
w = e 0 u
получим следующее линейное уравнение в частных производных:
dttu = a 2д xxu + b 2u,
где a = ,— ; b =
л/b
0 a0c0
В линии, свободной от искажения, Ь = 0, и приведенное выше выражение принимает вид
д ttu = a 'д **и- (2)
То есть (2) является одномерным однородным волновым уравнением. Рассмотрим смешанную задачу для (2) с заданными начальными условиями:
u(х,0) = f (х); д^(х,0) = F(х); 0 < х < I (3)
a0
и постоянными краевыми условиями, которые для удобства будем считать однородными:
u(0, t) = u(l, t) = 0; t > 0. (4)
Это не является ограничением общности, так как в случае неоднородных условий:
u (0, t) = u0; u (l, t) = u¡; t > 0,
где u0, ut - действительные числа, заменой искомой функции
u¡ u¡0
u = v + —-0 X + un
l 0
можно свести данную задачу к аналогичной относительно новой неизвестной функции v, но уже с однородными краевыми условиями.
Предположим, что функция f(x) дифференцируема на отрезке [0, ¡], а ее производная f'(x) удовлетворяет условию Липшица с константой L1 > 0, т. е. для любых (xi, x2) е [0, l] выполняется неравенство
\f'( Xi) - f'(x—)| < al Xi - X2i.
Функцию ^(x) на отрезке [0, l] будем предполагать кусочно-монотонной и удовлетворяющей условию Липшица с константой L2 > 0. Тогда существует единственно обобщенное решение смешанной задачи (2)-(4) [1, с. 144], которое представляется рядом
u(X,t) = ^(ak coskaat + bk sinkaat)sinkaX, (5)
k=i
где ю = п/l, и при любом натуральном k:
2 l 2 l
ak = — Jf(s)sinkasds; bk =-J^(s)sinkasds. (6)
l „ kna
0
Формула (5) малопригодна для приближенного вычисления решения, так как необходимо находить коэффициенты-интегралы (6), да и не существует общей оценки погрешности. Авторы статьи предлагают более эффективный метод решения этой задачи, основанный на использовании полилогарифмов [9, 10]. В дальнейшем будем использовать мнимую часть полилогарифма первого порядка на единичной окружности
Ы\х) = 1т Ь1(в1х) = £^^, хеЯ, (7)
к=1 к
и там же действительную часть дилогарифма
М2(х) = ЯеЬ2(в,х) = £С°?2кх, хеЯ. (8)
к=1 к
При X е [—п, п] эти периодические функции задаются выражениями:
где sgn(x)=
N\x) = 2(nsgn(x) - x); M2(x) = fc-ML - П2, (9)
X, x Ф 0;
\x\ - функция знака.
0, x = 0
Выполнив в первом из интегралов (6) интегрирование по частям и учитывая, что /(0) = /(I) = 0, получим
2 1
ak =—| f '(s) coskœsds.
kn 0
Используя это выражение для коэффициентов ак и Ьк (6), после подстановки их в (5), несложных преобразований и, принимая во внимание (7), найдем следующее представление для искомой функции:
/( x, t ) = -1J f '(s) + - F (s) (( N1 (œ( x - s + at )) + N1 (œ( x + s - at )) ) + f '(s) - - F (s) (( N1 (œ( x - s - at )) + N 1(ю( x + s + at )) )
(10)
Теперь выполним приближенное вычисление функции и(х, 1). Разобьем отрезок [0, I] на п равных частей точками хк = кк, к = 0, п (где к = 1/п - шаг разбиения) и заменим под знаком интеграла в правой части формулы (10) на каждом из частичных отрезков [хк-1, хк ], к = 1, п функции А'(х) и F(х) ее значениями в средней точке отрезка. В результате, учитывая, что первообразной функции Ы1(х) является, очевидно, функция -М2(х), получим:
u( x, t ) - ¿Ê Î ii-f [xk-1 + f + -F fxk-1 + h]! ( N1 (œ ( x - s + at )) + N1 (œ ( x + s - at ))) +
2n ^ J 1Г l k-1 2 J a
+1 f '[ xk-1 + f ( - - F [ xk-1 + f ( ) ( N1 (œ ( x - s - at )) + N1 (œ ( x + s + at )) ) I ds --
l Ê i(f(xk-1+ f( + aF fxk-1+ ff)I(M 22œ( x - s + at)) - M2(œ(x + s - at))) +
2п2 Ц I к-1 2( а
+(-/"^хк-1 + "кк)-^ [хк-1 + 2))(м2(®(х - ^ - а)) - М 2(ю( х + ^ + а))) Таким образом, приближенным решением задачи (2)-(4) является функция
». (х,')=¿2 х (( /'( хк-.+2 ) +^ ( хк-,+2))«
X (М2 (<(х - 5 + М)) -М2 (<(х + 5 - м))) +
+ (1'(хк-1 + ) - 1Р(хк-1 + 2))(М2 (<х - 5 - «0) -М2 (<(х + 5 + М)))
(11)
хк
хк-1
Найдем оценку допускаемой при вычислении решения по формуле (11) погрешности, принимая во внимание, что
Л
/(х, X) - ип (х, X) = | ((I'(5) - 1'(хк-1 + |) + М(Р(5) - р (хк-1 + 2
х(N (<(х - 5 + М)) + N (<(х + 5 - М))) +
(1 '(■5) - 1'( хк-,+ 2)-1 ( Е ( 5) - Е ( * + 2 ))
х(N (<(х - 5 - М)) + N (<(х + 5 + М))))й5 .
Тогда, учитывая липшициевость функций /'(х) и Е(х), а также тот факт, что ввиду (9)
N Ч х)| <-,
I 'I 2
получим
1 п хк
|и(х, X) - ип (х, X)| < —X |
к=1
1'(5) - 1'( хк-1
Е (5) - Р\хк-1 + 2
< (|N (<(х - 5 + М))| + |N (<(х + 5 - М))|) +
1 '(5 ) -1 '\хк-1 + 2
Е (5) - Е|хк-1
<(|N (<(х - 5 - М))| + N (<(х + 5 + М))|))й5 :
1 п хк (
' 2(*,
к=1 х
хк-1
<
5 хк-1 2
1 и
5 хк-1 2
п п
2 + 2
— + — Ш5 <
хк-1
ХДЧ + 1 Ч
1 Г И Ъ -1 + Ь2 )
2М
к-1
Следовательно, равномерно по x е [0, l] и t > 0
|и(x, t) -ип(x, t)| < l^ +^)h. (12)
2а
Проведенные выше исследования позволяют сформулировать следующее утверждение.
Теорема. При сделанных выше предположениях относительно функций f (x) и F(x) точное решение смешанной задачи (2)-(4) для однородного волнового уравнения представляется через мнимую часть полилогарифма первого порядка на единичной окружности (7) по формуле (10), а приближенное решение находится через действительную часть дилога-рифма (8) по (11). Абсолютная погрешность вычислений оценивается по (12), т. е. она равномерна по x е[0, l], t > 0 и имеет первый порядок
малости относительно шага разбиения h.
Пример. Найти приближенное решение следующей смешанной задачи:
д ttu = 4 д **и;
и(x,0) = sh 8т2га:; ди(x,0) =-4п cos2лx sh 8т2тсг; 0 < x < 1; u(0, t) = u(1, t) = 0; t > 0.
Непосредственной проверкой мы можем убедиться в том, что точным решением этой задачи является функция
u( x, ^ = e-sm4nt С(к2п^т2лх ^4я0.
Вычисление значений приближенного решения данной задачи по формуле (11) в узлах сетки (xk, t]); k, / = 0,10 (где xk = 0,1к; tj = 0,1/') показывает, что уже при п = 50 абсолютная погрешность меньше 0,002. Графики точного и(х, 0 и приближенного и50(х, ^ решений данной задачи в квадрате 0 < х, X < 1 представлены на рис. 1.
Рис. 1. График решений: а - точного u(x, t); b - приближенного u50(x, t) Fig. 1. The graphs of the decisions: a - exact one u(x, t); b - approximate one u50(x, t)
ВЫВОД
С помощью полилогарифмов найдено точное и приближенное представления решений задачи об электрических колебаниях в линии, свободной от искажения. Приведена эффективная равномерная оценка погрешности приближенного решения.
ЛИТЕРАТУРА
1. Кошляков, Н. С. Дифференциальные уравнения математической физики / Н. С. Кошляков, Э. Б. Глинер, М. М. Смирнов. М.: Физматгиз, 1962. 767 с.
2. Араманович, И. Г. Уравнения математической физики / И. Г. Араманович, В. И. Левин. М.: Наука, 1969. 288 с.
3. Смирнов, В. И. Курс высшей математики: в 2 т. / В. И. Смирнов. М.: Наука, 1974. Т. 2. 479 с.
4. Мышкис, А. Д. Лекции по высшей математике / А. Д. Мышкис. СПб.: Лань, 2007. 688 с.
5. Остапенко, В. Телеграфное уравнение. Краевые задачи / В. Остапенко. Саарбрюккен: LAP Lambert Academic Publishing, 2012. 272 с.
6. Новиков, Ю. Н. Электротехника и электроника. Теория цепей и сигналов, методы анализа / Ю. Н. Новиков. СПб.: Питер, 2005. 384 с.
7. Бычков, Ю. А. Основы теории электрических цепей / Ю. А. Бычков, В. М. Золотницкий, Э. П. Чернышев. СПб.: Лань, 2002. 464 с.
8. Дубнищев, Ю. Н. Колебания и волны / Ю. Н. Дубнищев. СПб.: Лань, 2011. 384 с.
9. Пыхтеев, Г. Н. Полилогарифмы, их свойства и методы вычисления / Г. Н. Пыхтеев, И. Н. Мелешко. Минск: Изд-во БГУ, 1976. 68 с.
10. Мелешко, И. Н. Специальные формулы для интегралов типа Коши и их приложения / И. Н. Мелешко. Минск: ВУЗ-ЮНИТИ, 1999. 197 с.
Поступила 28.10.2016 Подписана в печать 07.01.2017 Опубликована онлайн 28.07.2017
REFERENCES
1. Koshlyakov N. S., Gliner E. B., Smirnov M. M. (1962) Differential Equations of Mathematical Physics. Moscow, Fizmatgiz Publ. 767 (in Russian).
2. Aramanovich I. G., Levin V. I. (1969) Equations of Mathematical Physics. Moscow, Nauka Publ. 288 (in Russian).
3. Smirnov V. I. (1974) A Course of Higher Mathematics. Vol. 2. Moscow, Nauka Publ. 479 (in Russian).
4. Myshkis A. D. (2007) Lectures in Higher Mathematics. St.-Petersburg, Lan' Publ. 688 (in Russian).
5. Ostapenko V. (2012) Telegraph Equation. Boundary Value Problems. Saarbrucken, LAP Lambert Academic Publishing. 272 (in Russian).
6. Novikov Yu. N. (2005) Electrical Engineering and Electronics. Circuit and Signals Theory, Methods of Analysis. St.-Petersburg, Piter Publ. 384 (in Russian).
7. Bychkov Yu. A., Zolotnickij V. M., Chernyshev Je. P. (2002) Fundamentals of Electric Circuit Theory. St.-Petersburg, Lan' Publ. 464 (in Russian).
8. Dubnishchev Yu. N. (2011) Oscillations and Waves. St.-Petersburg, Lan' Publ. 384 (in Russian).
9. Pykhteev G. N., Meleshko I. N. (1976) Polylogarithms, their Properties and Calculation Method. Minsk, Belarusian State University Publ. 68 (in Russian).
10. Meleshko I. N. (1999) Special Formula for Integrals of Cauchy-Type Integrals and their Applications. Minsk, VUZ-YuNITI Publ. 197 (in Russian).
Received: 28 October 2016
Accepted: 7 January 2017
Published online: 28 July 2017