УДК 519.61 Вестник СПбГУ. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 2 МБС 65F22
Определение точек разрыва и величины скачка оригинала по его изображению по Лапласу*
А. В. Лебедева, В. М. Рябов
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Лебедева А. В., Рябов В. М. Определение точек разрыва и величины скачка оригинала по его изображению по Лапласу // Вестник Санкт-Петербургского университета. Математика. Механика. Астрономия. 2024. Т. 11 (69). Вып. 2. С. 316-323. https://doi.org/10.21638/spbu01.2024.205
Применение интегрального преобразования Лапласа для широкого класса задач приводит к более простому уравнению относительно изображения искомого оригинала. На следующем шаге возникает задача обращения, т. е. нахождения оригинала по его изображению. Как правило, осуществить этот шаг аналитически не удается. Возникает задача использования приближенных методов обращения. При этом приближенное решение представляется в виде линейной комбинации образа и его производных в ряде точек комплексной полуплоскости, в которой изображение регулярно. Однако оригинал, в отличие от изображения, может даже иметь точки разрыва. Несомненный интерес представляет задача разработки методов определения возможных точек разрыва оригинала и величины скачка оригинала в этих точках. В предлагаемых методах используются значения производных образа высокого порядка с целью получения удовлетворительной точности приближенных решений. Указаны приемы ускорения сходимости получаемых приближений. Приведены результаты численных экспериментов, иллюстрирующие эффективность предлагаемых методов.
Ключевые слова: интегральное преобразование Лапласа, задача обращения, точки разрыва оригинала, скачок оригинала.
1. Применение интегрального преобразования Лапласа приводит к более простому уравнению относительно изображения искомого оригинала. На следующем шаге возникает задача обращения, т. е. нахождения оригинала /(£) по его изображению Е(р) из уравнения
Теория преобразования Лапласа и аналитические методы его обращения содержатся в классических работах [1, 2].
Формула обращения задается интегралом Римана — Меллина:
где y — абсцисса сходимости интеграла Лапласа (1).
* Статья подготовлена при поддержке гранта Санкт-Петербургского государственного университета «Мероприятие 3» (Pure ID 75207094).
В приведенных примерах вычисления проводились с помощью пакета Maple 15 с использованием оборудования Ресурсного центра. Вычислительный центр СПбГУ (http://www.cc.spbu.ru/). © Санкт-Петербургский государственный университет, 2024
(1)
(2)
К сожалению, формула (2) мало пригодна для вычислений, и потому возникает необходимость в разработке и использовании приближенных методов обращения.
Методам обращения посвящены книги [3, 4]. Выбор метода обращения определяется той априорной информацией об оригинале и его образе, которой мы располагаем. Известно, что образ F(p) как функция комплексной переменной p регулярна в полуплоскости Re (p) > j, так что мы можем при построении метода обращения использовать значения образа либо в окрестности начала координат, либо в окрестности бесконечно удаленной точки, либо на вещественной оси, либо во всей полуплоскости регулярности и т. д.
Частные случаи реализации возможных ситуаций и аналитический вид приближений подробно описаны в книге [4]. Задача обращения некорректна, это, в частности, проявляется в том, что искомый оригинал может быть и разрывным, в то время как приближенное решение во всех методах обращения есть линейная комбинация образа и его производных в некоторых точках, являющихся бесконечно гладкими функциями в области регулярности.
В случае разрыва оригинала первого рода в точке t методы обращения доставляют значение (f (t + 0) + f (t — 0))/2.
Для вычисления величины f (t + 0) — f (t — 0) воспользуемся доказанным в книге [5] следующим утверждением.
Теорема 1 ([5], с. 298). Если интеграл
/• ж
A(p)= e-pt da(t) J о
сходится при некотором p, то
lim (-¡)"i(n) (^=a(i+0)-a(i-0).
2. Пусть оригинал f (t) имеет ограниченную вариацию на отрезке [0,T] при любом T > 0, его изображение равно F(p) и производная f '(t) также является функцией-оригиналом. Из условия теоремы 1 вытекает существование интеграла
/ e-pt df (t) о
и его равенство величине A(p) = pF(p) — f (0). По теореме 1 скачок оригинала в точке t равен
f (t + 0) — f (t — 0) = lim Ф„^), (3)
ф„м = Н)'»©- <4>
Если удается аналитически точно вычислять производные высокого порядка, то эти формулы пригодны для проведения расчетов. Заметим, что число точек разрыва может быть любым.
Пример 1. Рассмотрим функцию
1, 0 < t < 1, — 1, 1 <t < 2, 1, 2 <t < 3, 0, t> 3.
Ее изображение равно
F (p)
1 — 2exp(— p) + 2exp(—2p) — exp(—3p)
и функция (4) легко записывается в явном виде. На рис. 1, а представлен график функции Ф3000 (£), на котором отчетливо видны величины скачков оригинала в точках £ = 1, 2, 3. На рис. 1, б представлен график функции Ф500(£) — Ф490(£). Очевидно, функции Ф500(£) и Ф490(£) должны быть близки друг к другу и иметь вид кривой рис. 1, а, а их разность в точках разрыва близка к нулю и имеет в этих точках почти нулевую производную. Заметим, что точки разрыва можно найти как корни уравнения Ф500(£) — Ф490 (£) = 0, а затем в них вычислить функцию (4).
Рис. 1. Графики функций Ф3000(4) (а) и Ф500(4) - #490(4) (б).
3. Перейдем к случаю приближенного вычисления величины (4) и определению скорости сходимости метода (3).
Введем два оператора Виддера:
pn+1 t ^ Wn(f,t) = (-1 )"P—FM(p)
n!
n+1
\p=n/t
Sr^f,t) = (-ir1^F^(p)[p={n+l)/t
n = 1, 2,....
n = 0,1,...
(5)
Запишем их в виде сингулярных интегралов:
Wn(f,t)= / ivn(x)f(tx)dx, wn{x) = -—e
Jo n!
/»Ж>
Sn(f,t)= Sn(x)f (tx) dx, Sn(x) o
(n + 1)n+1
e-(n+1)xxn
Ядра wn(x), sn(x) удовлетворяют условиям теорем И. П. Романовского и Д. К. Фад-деева [6], откуда сразу следует сходимость, т.е.
lim Wn(f,t) = f(t), lim Sn(f,t) = f(t),
п^ж п^ж
если f (x) интегрируема и непрерывна в точке x = t.
p
~nxx7
Однако нас интересует не только факт сходимости, но и скорость сходимости, т. е. поведение величин
eW(f,t) = Wn(f,t) - f (t), eSn(f,t) = Sn(f,t) - f (t).
Свойства введенных операторов Виддера и методы их приближенного вычисления изучены в книге [4].
Теорема 2. Пусть оригинал f (t) имеет ограниченную вариацию на отрезке [0, T] при любом T > 0, его изображение равно F(p) и производная f '(t) также является функцией-оригиналом. Тогда при t > 0 справедлива формула
V2n lim л/п (Wn(f, t) — Sn-i(f, t)) = f(t+ 0) — f(t — 0),
где Wn(f,t),Sn(f,t) — операторы Виддера (5).
Доказательство. По теореме 1 скачок оригинала в точке t равен
По формуле Лейбница находим
Непосредственной проверкой убеждаемся в справедливости тождества
"п\ , n\tn Л
А(П) (?) = " ~ • (6)
Подставляя эту величину в (4) и заменяя факториал по формуле Стирлинга, придем к утверждению теоремы. Итак, с помощью приближений Шп) мы получаем представление не только о значении оригинала в некоторой точке, но и о наличии разрыва в ней.
Для приближенного вычисления величины — Бп—1(/,1) в книге [4] были
получены приближенные формулы
1 m
Wn(f,t) -Sn-iiLt) « (7(1 -гет(Я)) (rem(j))-n,
m z—' Vi /
j=i
cp(p) =pF(p), em(x) = exp ( — x ) ,
m
(7)
и приведены оценки их погрешности как функции входящих в нее параметров.
Очевидно, при наличии разрыва разность Wn(f,t)-Sn—i(f,t) в соответствии с теоремой 2 убывает, вообще говоря, как 1 /л/п. Значит, вычисление скачка с удовлетворительной точностью на основе формулы (6) может потребовать проведения вычислений для очень больших n, что вызовет большие затруднения при фактической реализации методов Виддера. Поэтому необходим детальный анализ алгоритма вычисления скачков, скорости его сходимости и приемы ее увеличения.
4. В книге [4] получена асимптотическая формула *n(i) = fit + о) - f(t -0) + (/'(t + 0) + f'(t - 0)) + О Q , n ->• 00.
Следовательно, в общем случае скорость сходимости невелика. Она увеличивается, если в точке Ь выполняется условие I' (Ь + 0) + I'(Ь — 0) = 0 (в частности, при I'(Ь + 0) = I'(Ь — 0) = 0, как в предыдущем примере).
Пример 2. Рассмотрим функцию
I (*) =
Ь, 0 < Ь < 1, 0, Ь> 1.
Ее изображение равно
Р (р)
1 — ехр(—р) — р ехр(— р)
р2
В этом случае в точке разрыва I'(Ь + 0) + I'(Ь — 0) = 0, и скорость сходимости уменьшается, что хорошо видно на рис. 2. Входящие в формулу (7) параметры таковы: г = 0.95, т = 251.
Рис. 2. Графики функций ФюоМ (а) и ФюоМ _ #97(4) (б).
В книге [4] предложены методы ускорения сходимости рассмотренных способов определения точек разрыва искомого оригинала и величины скачков в точках разрыва.
Замечание. Для вычисления величины (6) с помощью представления (7) можно прийти иным путем, если для вычисления входящих в формулу (4) производных воспользоваться известной формулой Коши вида
д(п\г)
п!
2^
д(м) ¿гл
{го - г)п+1'
д(г) = А(г),
где С — окружность с центром в точке г и таким радиусом, при котором вся линия интегрирования расположена в области регулярности изображения. Затем для вычисления интеграла применяется составная квадратурная формула прямоугольников. Однако на этом пути не видна скорость сходимости величины (4) к своему предельному значению.
5. Предположим, что с помощью описанной выше техники определены все точки разрыва оригинала в количестве т штук и для каждой из них вычислены пары
с
чисел (Ьк ,вк), к = 1, 2,...,т, где Ь1 < Ь2 < ...,< Ьт — точки разрыва, а вк I(Ьк + 0) — I(Ьк — 0) — величины скачков оригинала в этих точках.
Построим новую последовательность изображений {Рк(р)}т=1 по формулам
= ^(р) - ^6ХР( Р*к\ к = 0,1,..., т — 1, =
р
Изображение Рт(р) соответствует непрерывному оригиналу ^(Ь), получающемуся из искомого оригинала I(Ь) удалением всех разрывов. Функция ^(Ь) находится с помощью известных методов обращения. Затем необходимо к этой функции добавить все скачки вк.
Пример 3. Рассмотрим функцию
I (*)
0, 0 < Ь < 1, Ь, Ь> 1.
Ее изображение равно
Р (р) =
- + Л ) ехр(-р)-р р2
Рис.3. Приближение по Виддеру (а); функция (б); функция
!(4) (в).
Соответствующий ему оригинал, найденный методом Виддера, представлен на рис. 3, а. На нем точка разрыва практически отсутствует.
Имеется одна точка разрыва, и в ней ¿1 = 1,в1 = 1. Следовательно,
р р 2
Соответствующий этому изображению сглаженный оригинал Il(t) изображен на рис. 3, б.
После добавления к нему скачка в точке Ь =1 находим искомый оригинал, представленный на рис. 3, в.
Литература
1. Лаврентьев М.А., Шабат Б. В. Методы теории функций комплексного переменного. Москва, Лань (2002).
2. Крылов В. И., Скобля Н.С. Методы приближенного преобразования Фурье и обращения преобразования Лапласа. Москва, Наука (1974).
3. Cohen A.M. Numerical methods for Laplace transform inversion. New York, Springer (2007).
4. Рябов В.М. Численное обращение преобразования Лапласа. Санкт-Петербург, Изд-во С.-Петерб. ун-та (2013).
5. Widder D. V. The Laplace transform. Princeton (1946).
6. Натансон И. П. Теория функций вещественной переменной. Москва, Наука (1974).
Статья поступила в редакцию 4 сентября 2023 г.;
доработана 11 октября 2023 г.; рекомендована к печати 9 ноября 2023 г.
Контактная информация:
Лебедева Анастасия Владимировна — канд. физ.-мат. наук, доц.; a.v.lebedeva@spbu.ru Рябов Виктор Михайлович — д-р физ.-мат. наук, проф.; v.ryabov@spbu.ru
Determination of breakpoints and the magnitude of the jump of the original according to its Laplace image*
A. V. Lebedeva, V. M. Ryabov
St. Petersburg State University, 7—9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Lebedeva A. V., Ryabov V. M. Determination of breakpoints and the magnitude of the jump of the original according to its Laplace image. Vestnik of Saint Petersburg University. Mathematics. Mechanics. Astronomy, 2024, vol. 11 (69), issue 2, pp. 316-323. https://doi.org/10.21638/spbu01.2024.205 (In Russian)
The application of the integral Laplace transform for a wide class of problems leads to a simpler equation for the image of the desired original. At the next step, the inversion problem arises, i.e., finding the original by its image. As a rule, it is not possible to carry out this step analytically. The problem arises of using approximate inversion methods. In this case, the approximate solution is represented as a linear combination of the image
*This paper was prepared with the support by a grant from St. Petersburg State University "Event 3" (Pure ID 75207094).
Research was carried out using computational resources provided by Resource Center "Computer Center of St. Petersburg State University" (http://www.cc.spbu.ru/en).
and its derivatives at a number of points of the complex half-plane in which the image is regular. However, the original, unlike the image, may even have break points. Of undoubted interest is the task of developing methods for determining the possible break points of the original and the magnitude of the original jump at these points. The proposed methods use the values of high-order image derivatives in order to obtain a satisfactory accuracy of approximate solutions. Methods for accelerating the convergence of the obtained approximations are indicated. The results of numerical experiments illustrating the effectiveness of the proposed methods are presented.
Keywords: Laplace integral transform, inversion problem, original discontinuity points, original jump.
References
1. Lavrent'ev M. A., Shabat B. V. Methods of the theory of functions of a complex variable. Moscow, Lan' Publ. (2002). (In Russian)
2. Krylov V. I., Skoblya N. S. Methods of the approximate Fourier transform and the inversion of the Laplace transform. Moscow Publ. (1974). (In Russian)
3. Cohen A.M. Numerical methods for Laplace transform inversion. New York, Springer (2007).
4. Ryabov V. M. Numerical inversion of the Laplace transform. St. Petersburg, St. Petersburg University Press (2013). (In Russian)
5. Widder D. V. The Laplace transform. Princeton (1946).
6. Natanson I. P. Theory of functions of a real variable. Moscow, Nauka Publ. (1974). (In Russian)
Received: September 4, 2023 Revised: October 11, 2023 Accepted: November 9, 2023
Authors' information:
Anastasia V. Lebedeva — a.v.lebedeva@spbu.ru Victor M. Ryabov — v.ryabov@spbu.ru