CC BY
40
DOI: 10.21122/2227-1031 -2017-16-4-343-347 УДК 519.6
К приближенному интегрированию сильно осциллирующих функций
Докт. физ.-мат. наук, проф. И. Н. Мелешко1), Д. А. Нифонтова4, В. В. Сорокин4
^Белорусский национальный технический университет (Минск, Республика Беларусь)
© Белорусский национальный технический университет, 2017 Бе1аг^ап National Technical University, 2017
Реферат. Построены и исследованы простейшие приближенные формулы для численного интегрирования функций, содержащих осциллирующие множители специального вида с параметром. Общие квадратурные формулы в этом случае могут быть использованы только при достаточно малых значениях параметра. Следовательно, чтобы получить формулы численного интегрирования, пригодные при изменении параметра в широких границах, необходимо заранее учитывать наличие сильно осциллирующих множителей. Это можно сделать, принимая, например, такие множители за весовые функции. Кроме того, поскольку параметр способен принимать значения, которые заранее предвидеть не всегда можно, приближенные формулы для вычисления таких интегралов необходимо строить так, чтобы они содержали этот параметр в буквенном виде и были пригодны для вычисления при любых, в частности при больших, значениях параметра. Вычислительные правила, обладающие такими свойствами, обычно получают путем разбиения промежутка интегрирования на элементарные с последующим приближением плотности интеграла на каждом элементарном промежутке многочленами первой, второй и третьей степеней, принимая при этом осциллирующие множители за весовые функции. В статье рассмотрен тот вариант, когда плотность интегралов на каждом элементарном промежутке аппроксимируется многочленом нулевой степени - константой, равной значению плотности в середине этого промежутка. Попутно сконструирована одна приближенная формула для вычисления несобственного интеграла по бесконечному промежутку от функции, содержащей осциллирующий множитель специального вида. При этом предполагали, что плотность несобственного интеграла достаточно быстро стремится к нулю, когда модуль аргумента неограниченно возрастает. Другими словами, она считается пренебрежимо малой вне некоторого конечного отрезка. Получены равномерные по параметру оценки погрешностей приближенных формул, позволяющие вычислять интегралы с заданной точностью.
Ключевые слова: интегралы, приближенное интегрирование, приближенные формулы, осциллирующая функция
Для цитирования: Мелешко, И. Н. К приближенному интегрированию сильно осциллирующих функций / И. Н. Ме-лешко, Д. А. Нифонтова, В. В. Сорокин // Наука и техника. 2017. Т. 16, № 4. С. 343-347. DOI: 10.21122/2227-10312017-16-4-343-347
Approximate Integration of Highly Oscillating Functions
I. N. Melashko1), D. A. Nifontova1), U. U. Sorokin1)
1)Belarusian National Technical University ^irak, Republic of Be^rus)
Abstract. Elementary approximate formulae for numerical integration of functions containing oscillating factors of a special form with a parameter have been proposed in the paper. In this case general quadrature formulae can be used only at sufficiently small values of the parameter. Therefore, it is necessary to consider in advance presence of strongly oscillating factors in order to obtain formulae for numerical integration which are suitable in the case when the parameter is changing within wide limits. This can be done by taking into account such factors as weighting functions. Moreover, since the parameter can take values which cannot always be predicted in advance, approximate formulae for calculation of such integrals should be constructed in such a way that they contain this parameter in a letter format and they are suitable for calculation at any and particularly large values of the parameter. Computational rules with such properties are generally obtained by dividing
Адрес для переписки
Мелешко Иван Николаевич Белорусский национальный технический ул. Я. Коласа, 12/2,
220013, г. Минск, Республика Беларусь Тел.: +375 17 292-82-73 kafvm2@bntu.by
Address for correspondence
Melashko Ivan N.
Belarusian National Technical University 12/2, Ya. Kolasa str., 220013, Minsk, Republic of Belarus Tel.: +375 17 292-82-73 kafvm2@bntu.by
Наука
итехника. Т. 16, № 4 (2017)
an interval of integration into elementary while making successive approximation of the integral density at each elementary interval with polynomials of the first, second and third degrees and taking the oscillating factors as weighting functions. The paper considers the variant when density of the integrals at each elementary interval is approximated by a polynomial of zero degree that is a constant which is equal to the value of density in the middle of the interval. At the same time one approximate formula for calculation of an improper integral with infinite interval of the function with oscillating factor of a special type has been constructed in the paper. In this case it has been assumed that density of the improper integral rather quickly goes to zero when an argument module is increasing indefinitely. In other words it is considered as small to negligible outside some finite interval. Uniforms in parameter used for evaluation of errors in approximate formulae have been obtained in the paper and they make it possible to calculate integrals with the required accuracy.
Keywords: integrals, approximate integration, approximate formulas, oscillating function
For citation: Melashko I. N., Nifontova D. A., Sorokin U. U. (2017) Approximate Integration of Highly Oscillating Functions. Science and Technique. 16 (4), 343-347. DOI: 10.21122/2227-1031-2017-16-4-343-347 (in Russian)
Введение
При решении разных задач прикладной математики, механики, физики и техники приходится вычислять интегралы вида:
J с (X) = J ( f; X)= J f (t)cos Xtdt;
-T T
Js (X) = Js (f; X) = J f (t)sin Xtdt;
-T T
Je (*)= Je (f; X)= J f (tК
-T
J(X) = J(f; X)=J f (t)e':
jxt,
ш,
(1) (2)
(3)
(4)
где X - положительный числовой параметр [1-5].
Формулы для приближенного вычисления подобных интегралов строятся обычно путем разбиения промежутка интегрирования на элементарные с последующим приближением функции ДО на каждом элементарном промежутке многочленами первой, второй и третьей степеней, принимая при этом осциллирующие множители за весовые функции [6-10]. Авторами исследован тот вариант, когда функция ДО на элементарном промежутке аппроксимируется многочленом нулевой степени -константой, равной значению ДО в середине этого промежутка.
1. Зададим на промежутке [- Т, Т] систему
2Т
точек ^ = кк; к = -и,...,-1,0,1,..., и; к =-
2п +1
и аппроксимацию функции ДО на этом промежутке определим формулой
f (t) - f (t) =£e* (t) f (tk),
(5)
k=-n
в которой ek (t) = 1, если t e
h 2 '
tk--, tk +-
и ek (t) = 0, когда t £
tk 2' tk + 2
Нетрудно убедиться, что если ДО непрерывна на отрезке [- Т, Т] , то
/(х) -/ (х )| <ю( /; к ), х е [-Т, Т ], (6)
где ю(f; h) = max If(t")- f (t')| - модуль
v ' t\t"e[-T, T]' 1
непрерывности функции f(t).
Если же ДО - непрерывно дифференцируемая функция на этом отрезке, то с помощью формулы Тейлора легко установить, что
f (t) - f (t)
< M h' t e[-T' T ]' (7)
где Ml = max f'(t).
te[-T, T ]' 1
2. Подставив в (3) для интеграла Je (f; Я)
вместо функции ДО ее приближение формулой (5), получим приближенное уравнение для этого интеграла
J,
(f; X)«Je(f; X) = Je( X)=^Ake)(X)f(tkX (8)
k=-n
2
в котором коэффициенты А1(в) (X) = | в^Ж,
к
к — к 2
к = - и, ...,-1, 0,1, ..., п .
Вычислив интегралы в правой части последнего равенства и выполнив несложные преобразования, находим единое представление коэффициентов
Наука
итехника. Т. 16, № 4 (2017)
h
U, +
A<e)(X) = -sin -^eiXt>
X 2n +1
(9)
Коэффициенты из (9) подставим в правую часть (8). В результате приближенная формула для интеграла (3) запишется в виде
Je (f; X)« 2 £
( . XT eiXtk Л sin-
k=-n
2n +1 X
f(tk). (10)
Оценим погрешность приближенной формулы (10).
Теорема 1. Если функция ДО непрерывна на отрезке [- Т, Т], то для всех X имеет место оценка погрешности
\1е (X)- Зе (X) < 2Т«(/; И). (11)
Если же fit) непрерывно дифференцируемая функция, то
\1е (X)- 1е (Х)|< ТМИ (12) Доказательство. Очевидно, что \Зе (X)- 1е (Х)| = |Зе (/; X)- Зе (/; х)| =
= Я / (')- / (')] е ^
<
(13)
< ¡f (t)-f (t М Iе
i X t
= 2T шах I/(t)-f (t)|.
te[-T; T]l 4 7 4 /l
Из соотношений (13) и неравенств (6), (7) вытекают оценки (11), (12).
3. Заменим функцию ft) в интегралах (1), (2) по формуле (5). В результате получим следующие приближенные выражения для инте-грал°в Jc (f; X), Js (f; X) :
J с (/; X)- Jc (/; X)=Jc (X) = £ 4c) (X) / (t*); (14)
k--n
Js (/; X) - Js (f; X) = Js (X) = £ 4s) (X)/(t*). (15)
k=-n
Коэффициенты 4c)(X) и Al*'(X) соответственно равны:
Наука
итехника. Т. 16, № 4 (2017)
2 2 Akc)(X) = í cos Xtdt; A(s)(X) = í sin Xtdt,
k = - n, ...,-1, 0,1, ..., n.
Выделив действительную и мнимую части в формуле (9), находим следующие представления для этих коэффициентов:
2 XT
Akc) (X) =т-sin -—:cos Xtk; X 2n +1
2 XT Aks)(X) = -sin — sin Xtk, X 2n +1
k = - n, ...,-1, 0,1, ..., n.
Подставим полученные выражения А^ (X)
и Ак-'(X) в правые части (14) и (15) соответственно. В результате приближенные формулы для интегралов (1) и (2) запишутся в виде:
tí f i\ о V ( • XT cos Xt,
X
k 14tk); (16)
Js (f; X)« 2 £(si
k=- n
XT sin Xt,
sin-
2n +1 X
(tk). (17)
Погрешности приближенных формул (16), (17) можно сверху оценить теми же величинами, что и погрешность приближенного уравнения (10).
Теорема 2. Пусть функция fit) непрерывна на отрезке [- T, T ]. Тогда для всех X справедливы неравенства:
\JC (X)- Jc (X)|< 2T«(f; h); \JS (X)- J (X)| < 2T«(f; h).
(18)
Если же fit) является непрерывно дифференцируемой функцией на этом отрезке, то
\J.
(X) - Jc (X)| < TMA I Js (X) - J s (X)| < TM 1h . (19)
Доказательство. Запишем очевидные соотношения:
Jc (X)- Jc (X)| = Jc (f; X)- Jc (f; X)|
J [ f (t)- f (t)] cos Xtdt
< max Jf (t)-f (t)| J |cosXt|dt <
<2T max f (t)-f (t) ;
te[-T; T]l 4 7 4
Js (X)- Js (X)| = Js (f; X)- Js (f; X)|
j
J[ f (t)- f (t)] sin X tdt
(20)
< max If (t)- f (t)| Jin Xt\dt <
te[-T; T]l V ' I J 1
< 2T max I f (t)-f (t)|.
te[-T; T]l 4 7 4 /l
Оценивая в (20) max f (t)- f (t)
te[-T; T]
с помо-
щью (6), (7), получаем неравенства (18), (19).
4. Предположим, что функция /Х) в интеграле (4) достаточно быстро стремится к нулю при |х| ^ <». Тогда, считая ее пренебрежимо малой вне некоторого конечного отрезка [-Т, Т], в пределах заданной точности справедливо приближенное равенство
Т
3(X) = | /(1)в1хйХ = 3е (X). (21)
-Т
Подставив в (21) приближенное представление интеграла (3) формулой (10), получаем приближенное уравнение для интеграла (4)
J (X)« Je (X) = 2 £
k=-n
XT
iXtk \
sin-
2n +1 X
f(tk). (22)
При оценке погрешности приближенной формулы (22) следует учитывать погрешность приближенного уравнения (21). Будем предполагать, что при больших значениях Х выполня-
ется неравенство
C
f (t) =
t
0 <C <<х> 5 > 0.
(23)
Тогда IJ(X)- Je (X)| =
T
J f (t) eiXtdt + J f (t) eiXtdt
< 2J f (t) dt < 2C J-
2C
(24)
1+5
5T
8 '
Далее
| J (X) - Je (X)| = IJ (X) - Je (X) + Je (X) - Je (X)| < < |J (X)-Je (X) + I Je (X)-Je (X)|.
Если функция /Х) имеет непрерывную производную на промежутке [- Т, Т ], то с учетом неравенства (12) теоремы 1 приходим к следующей оценке погрешности приближенной формулы (22):
|J(X)- Je (X)|<
2C
5T
5 . TMih.
Примечание. Очевидно, что приближенную формулу (22) и оценку ее погрешности можно использовать для вычисления интеграла (4) только в том случае, когда функция f(t) достаточно быстро убывает с возрастанием |t| (5 достаточно велико).
5. Пример. Коэффициенты тригонометрического ряда Фурье функции f (t) = et определяются формулами:
1 п t (-1)" 2shn
— I e cos mtdt = n J
( fet
1
— I el sin mtdt = n J
n(1 + m2)
(-1)m+: 2mshn n (1 + m2)
Jc
V v
—, m n
J J
J ( et
V v
' m
J
Результаты численного эксперимента по точным формулам и с помощью приближенных уравнений (16), (17) при п = 100 приведены в табл. 1
Таблица 1
m Jc (m) Jc (m) Js (m) J s (m)
1 -3,6760800 -3,6756300 3,6760800 3,6762300
10 0,0727936 0,0721969 -0,7279360 -0,7279660
100 0,0007351 0,0007262 -0,0735142 -0,0735126
Примечание. Предположим, что значения функции fx) в узлах квадратурной формулы (8) вычисляются приближенно, т. е. вместо f (tk)
имеем f (tk), так что погрешность вычисления
Наука
итехника. Т. 16, № 4 (2017)
равна гк = f (tk)-f (tk). Пусть |et|<s, к = -и,..., -1, 0,1,... n. Тогда для погрешности вычисления квадратурной суммы в (8) получается неравенство
£ л, (e)(x) f (tk )-£ a, (e)(x) f {tk)
<s£| A, (X)|,
из которого следует, что точная верхняя граница погрешности вычисления квадратурной
и
суммы в (8) пропорциональна Ак(в)(Х).
-и
Оценим эту сумму числом, не зависящим от X:
a A,(e)wi=£
2
J e'Xtdt
n .2
■a J
iXt
- n h
n r
a J dt = 2T.
-n h
Следовательно, при всех X справедлива оценка
£ Ak (e)(X)f (tk )-£ Ak (e)(X) f (tk)
<s2T.
Очевидно, что погрешности в вычислениях квадратурных сумм (14), (15) оцениваются такой же величиной. Это означает, что при больших и погрешности в вычислениях упомянутых квадратурных сумм имеют тот же порядок, что и погрешность в вычислении функции /Х). В таких случаях говорят, что квадратурные формулы численно устойчивы.
ВЫВОДЫ
1. Для численного интегрирования функций, содержащих сильно осциллирующие множители специального вида с параметром, построены простейшие приближенные формулы.
2. Попутно сконструирована одна приближенная формула для вычисления несобственного интеграла по бесконечному промежутку от функции с осциллирующим множителем специального вида в случае, когда плотность несобственного интеграла достаточно быстро стремится к нулю при неограниченном возрастании модуля аргумента.
3. Получены равномерные по параметру оценки погрешностей приближенных формул,
■■ Наука
итехника. Т. 16, № 4 (2017)
позволяющие вычислять интегралы с заданной точностью. Приведен численный анализ квадратурных сумм на устойчивость.
ЛИТЕРАТУРА
1. Канторович, Л. В. Приближенные методы высшего анализа / Л. В. Канторович, В. И. Крылов. М.; Л.: Физматгиз, 1962. 708 с.
2. Лаврентьев, М. А. Методы теории функций комплексного переменного / М. А. Лаврентьев, Б. В. Шабат. М.: Наука, 1973. 736 с.
3. Фихтенгольц, Г. М. Курс дифференциального и интегрального исчисления / Г. М. Фихтенгольц. М.: Наука, 1970. Т. 3. 656 с.
4. Боголюбов, А. Н. Задачи по математической физике /
A. Н. Боголюбов, В. В. Кравцов. М.: Изд-во МГУ, 1988. 349.
5. Ланцош, К. Практические методы прикладного анализа / К. Ланцош. М.: Физматгиз, 1956. 524 с.
6. Крылов, В. И. Приближенное вычисление интегралов /
B. И. Крылов. М.: Наука, 1967. 500 с.
7. Крылов, В. И. Методы приближенного преобразования Фурье и обращения преобразования Лапласа / В. И. Крылов, Н. С. Скобля. М.: Наука, 1974. 224 с.
8. Завьялов, Н. С. Методы сплайн-функций / Н. С. Завьялов, Б. И. Квасов, В. П. Мирошниченко. М.: Наука, 1980. 352 с.
9. Бахвалов, Н. С. Численные методы / Н. С. Бахвалов, Н. П. Жидков, Г. М. Кобельков. М.: Наука, 1987. 598 с.
10. Калиткин, Н. Н. Численные методы / Н. Н. Калиткин. М.: Наука, 1978. 512 с.
Поступила 16.03.2016 Подписана в печать 27.05.2016 Опубликована онлайн 28.07.2017
REFERENCES
1. Kantarovich L. V., Krylov V. I. (1962) Approximate Methods for Highest Analysis. Moscow-Leningrad, Phiz-matgiz Publ. 708 (in Russian).
2. Lavrentiev M. A., Shabat B. V. (1973) Methods for Theory of Functions of a Complex Variable. Moscow, Nauka Publ. 736 (in Russian).
3. Fikhtengolts G. M. (1970) Calculus course. Vol. 3. Moscow, Nauka Publ. 656 (in Russian).
4. Bogolyubov A. N., Kravtsov V. V. (1988) Exercises in Mathematical Physics. Moscow, Moscow State University. 349 (in Russian).
5. Lantsosh K. (1956) Practical Methods for Application Analysis. Moscow, Phizmatgiz Publ. 524 (in Russian).
6. Krylov V. I. (1967) Integral Approximation. Moscow, Nauka Publ. 500 (in Russian).
7. Krylov V. I., Skoblya N. S. (1974) Methods of Approximate Fourier Transformation and Inverse Laplace Transform. Moscow, Nauka Publ. 224 (in Russian).
8. Zavyalov N. S,. Kvasov B. I., Miroshnichenko V. P. (1980) Methods of Spline-Functions. Moscow, Nauka Publ. 352 (in Russian).
9. Bakhvalov N. S., Zhidkov N. P., Kobelkov G. M. (1987) Numerical Technique. Moscow, Nauka Publ. 598 (in Russian).
10. Kalitkin N. N. (1978) Numerical Technique. Moscow, Nauka Publ. 512 (in Russian).
Received: 16.03.2016 Accepted: 27.05.2016 Published online: 28.07.2017
h
t,, +
k
2
