ДОКЛАДЫ АКАДЕМИИ НАУК РЕСПУБЛИКИ ТАДЖИКИСТАН _______________________________________2013, том 56, №7___________________________________
МАТЕМАТИКА
УДК 517.5
Г.А.Юсупов, А.А.Шабозова
ТОЧНЫЕ ОЦЕНКИ ПРИБЛИЖЁННОГО ИНТЕГРИРОВАНИЯ КРИВОЛИНЕЙНЫХ ИНТЕГРАЛОВ ПЕРВОГО РОДА НА НЕКОТОРЫХ
КЛАССАХ ФУНКЦИЙ И КРИВЫХ
Таджикский национальный университет
(Представлено академиком АН Республики Таджикистан М.Ш.Шабозовым 14.05.2012 г.)
Для приближённого вычисления криволинейного интеграла первого рода построена усложнённая квадратурная формула Симпсона и для классов дифференцируемых функций, задаваемых модулями непрерывности, найдена оценка её погрешности.
Ключевые слова: криволинейные интегралы первого рода - усложнённая квадратурная формула Симпсона - модуль непрерывности.
Рассмотрим задачу приближённого интегрирования криволинейного интеграла первого рода в форме линейной комбинации конечного числа значений подынтегральной функции
N
J f (M )ds = £ ptf (Mt) + Rn (f ,Г), (1)
Г k=0
где Г - некоторая плоская простая кривая, f (M) = f (x, y) — определённая на кривой Г функция, Mk єГ - узлы, pk - коэффициенты, R^(f,Г) - погрешность квадратурной формулы на функции
f (M).
Через Жд (L) обозначим класс плоских спрямляемых кривых Г, у которых длина равна L , кривизна кусочно-непрерывна и все они расположены в области Q = {(x, y) : x2 + у2 < L2}.
Известно [1], что параметрические уравнения кривой Г є^ (L), отнесённой к длине дуги s как параметру, в прямоугольной системе координат Oxy имеют вид
x = x(s), у = y(s), 0 < s < L. (2)
Обозначим через sfc є [0,1], k = 0, N значения длины дуги s кривой Г, которые соответст-
вуют точкам Mk єГ, и перепишем формулу (1) в виде
L n
J f (x(s), у( s)) ds = £ pkf (x(sk ), y( sk )) + Rn (f, Г),
0 k=0
где x = x(s), у = y(s) - параметрические уравнения кривой Г, представленные в виде (2).
Адрес для корреспонденции: Юсупов Гулзорхон Амиршоевич. 734025, Республика Таджикистан, г.Душанбе, пр. Рудаки, 17, Таджикский национальный университет. E-mail: [email protected]
В этой статье будем рассматривать усложнённые квадратурные формулы [2], которые строятся следующим образом. Отрезок [0, Ь] делят на N равных частей точками
5^ = кк, к = Ь / N, к = 0, N и на каждом интервале [5^ 15 5^ ], к = 0, N применяют заранее квадратурную формулу с узлами , < /, < /2 < • • • < //н < Л'^ и коэффициентами рі (/ = 1,2, • • •, да). В результате получим усложнённую квадратурную формулу
Ь
| / (х(5) у(5)) ж = Е ь (5к-^ sk ■ /)+^ (/),
к=1
где
Ь (5к-і, ^ : /) = Цр,/(х&X У(Іг )) •
В этой статье мы дадим оценку погрешности усложнённой квадратурной формулы Симпсона, имеющей вид:
- 2к - и 1 (2к -1 || +
|/ (х( 5) у(5) ) ^ = ~^{/ (х(0) у(0))+4Е / [х [ ^кщ1 Ь ], у
6^ 4 2N Г V 2 N
N^1 Г /^к^ /^кМ I +2Е / X [ кЬ |, У [ кЬ]| + / (т, У( Ь))} + с/, Г) (3)
к=1
V V
V
N
для классов дифференцируемых на отрезке [0, Ь] сложных функций / (х(5), у(5)) = ¥(5) одного
переменного, задаваемых модулями непрерывности. С этой целью вводим оператор набла 'V ", полагая
V/ = V/ (х(5), У(5))-.= 3/- -Iх + ^ -Iу = ^ (5),
ох 15 |У 15
I2/ Г |2 ~д2/ <Лх dy д2/ Г dy |2 .
V2/ = V(V/) = ^-- — I + 2^------------------- + ^г' — | := ^ (5),
дх V d5) дхду d5 d5 ду ^
и для k = 3,4,... положим
d5
V/ = у(ък-1/). (4)
Через := Шд>>Н<а\0,Ь\ к е N обозначим множество функций /{М) = /(х,у), у
которых почти всюду в области Q существуют частные производные
дк Г
дх 3ду3
0
г=1
и для любых двух точек 5 ,5 е [0, Ь] функция (4) удовлетворяет условию
(5)
где со(д') - заданный модуль непрерывности, то есть непрерывная неубывающая функция, в нуле равная нулю.
Очевидно, что в обозначении сложной функции ¥(5) = / (x(s), у(5)) производные
р(к) = ук/ (х(^), у(^)) , а условие (5)
означает, что
¥№)(5 ) - ¥№)(5 ) <о(|5 - 5 |), Уз , 5 е[0, Ь],
то есть функция ¥(^) е Ж Н°[0,Ь] (см.,например, [2], стр. 19). Записав для произвольной функции / еЩ2' Н как функции одного переменного / (^(^), у(^)) формулу Тейлора с остаточным членом
в интегральной форме Коши, принимая во внимание, что формула Симпсона точна для всех многочленов третьей степени, остаток квадратурной формулы (3) запишем в виде
^ (/,Г) = \ (2^ 2)) 22(5)й2 = 2 | ¥ 'ОЖ
д2(5) =
(5 -кк)^5 -кк -0 кк < 5 < кк + к / 2,к^Ь /N,
(5 - (к + 1)к) Г 5 - (к + 1)к + к\ (к + 1)к - к / 2 < 5 < (к + 1)к.
Так как
д2 (кк + 5) = д^), 0 <5</г, к = 1,2,---,Ы-I,
то для любого / е ЖQ2)Н0
имеем:
|RN(/,Г)|< N вир
/ еЖ(2)Нш
к/2
| У2/ (x(s), У(5)) #2(№
= N вир
д>еН0}
к/2
(6)
Следуя схеме рассуждения работы [4], положим
<
0
9з(5) = | 9і(і)Ж =152 Г5 -И\ 0 < 5 < И
3 ^ 2) 2
и заметим, что д3 (5) строго убывает на [0, к / 3) и строго возрастает на (к / 3, к / 2], причём д3 (к / 2) = 0 . Введём в рассмотрение функцию р((), которая определяется равенством:
03(я) = д3((р(э)), (0 < 5 < к / 3, к / 3 < р(^) < к / 2).
Используя лемму 5.1 из [3], будем иметь
Ь/(2 N ) Ь/)
8ир
д>єН®
| д2(5)Ф¥5 < | 192(5)\®(р(5)-5)
(7)
Из (6) и (7) следует, что для любой функции / є Ж2)Н®
Ь/(3N)
\ ^(/,Г) \< N І \ 9і(5) \®(р(5) - 5)^
где б)(5) - заданный модуль непрерывности.
Если ®(5) выпуклый вверх модуль непрерывности, то в классе Н®[0, Ь] существует возрастающая на [0,И / 2] функция <р0(5) , которая реализует в (7) точную верхнюю грань. Она имеет вид:
%(5) = <
И/2
(р(ї) - Ї) dt, 0 < 5 < к / 3,
5
5
(t- р_1(0)dt, к / 3 < 5 < к / 2,
И/2
где р 1 (х) - функция, обратная к р(х) .
С помощью равенств
(р0(кк + = ср^), 0<5</гу, к = \,2,---,Ы-I
продолжим (t) на весь отрезок [0, Ь]. Очевидно, что (р0 є Н®[0, Ь]. Пусть /й (х(5), у(5)) = ^0(5) - функция, у которой
(5) = ^ (x(5), у(5)) = 0 <5 < Ь
Тогда, используя схему рассуждений работы [4], нетрудно подсчитать, что
0
h/2
h/3
( h
RN (f0 ,Г) = N j f0(s)92(s)ds = N j s I- - s]®(p(s) - s) ds
0 0 V 3
27 N
hr js(1 - s)a
P\ — \-— ds. 1 3N J 3N I
Если о(^) - выпуклый вверх модуль непрерывности, то в силу леммы 5.2 из [3] имеем:
й/2 й/2 | V2/ (x(s), у(5)) д2(^ = | ^0)° '(s)ds,
sup
gW^ k) H
где щ(5) - убывающая перестановка функции | д3 (5) | . В нашем случае у/($) =| У С5) | (0 < 5 < к / 2) и для любой / е Ж2)Н°
|RN (/,Г)|< N
—2
j V/ (2(s), y(s)) g2(s)ds < N j 12(s) I ® (s)ds
-2
L1 2
1
I
sL
2N
а {—J ds. 2 N J
(8)
Для введённой выше функции / (х(5), у(^)) здесь имеет место знак равенства. Если учесть (8), то получим окончательно для выпуклого вверх модуля непрерывности ) :
2 1
sup | Rn(f, Г) |= —L—- j <j 1 - cos ■wPh" 648N3 J 1
л/э ж
arccos—s------
2 6
Y
'(2LN N ds
Оценку сверху для (/,Г), удобную для практического использования, можно получить
таким образом. В силу монотонности у (х) из (8) имеем
А/2
J | y (s') | а ' (s)ds <| y(0) | J а ' (s)ds,
(9)
где а определяется из равенства:
h/2
а | У(0) |= J | y(s) | ds.
9L L3
Легко подсчитать, что а =------------,| y(0) |= max | о,(s) |=------------ . Тогда для любой f g W2)H' из
32N 0<s<L/(2 N ) 3 162N3' e
неравенства (9) следует, что
0
0
|Rn(f, Г)| < —L— • а (—].
1 N 1 162 N \ 32 N J
Поступило 16.05.2012 г.
ЛИТЕРАТУРА
1. Фиников С.П. Курс дифференциальной геометрии. - М.: Наука, 1952.
2. Никольский С.М. Квадратурные формулы. - М.: Наука, 1979, 254 с.
3. Корнейчук Н.П. - Изв. АН СССР, сер. матем., 1971, т.35, №1, с. 93-124.
4. Бусарова Т.Н. - В сб.: Вопросы теории приближений функций и её приложений. - ИМ АН УССР, Киев, 1976.
Г.А.Юсупов, А.А.Шабозова
БАХ,ОИ АНИЦИ ИНТЕГРОНИИ ТАЦРИБИИ ИНТЕГРАЛ^ОИ КАЧ,ХАТТАИ НАВЪИ ЯКУМ БАРОИ БАЪЗЕ СИНФИ ФУНКСИЯ^О ВА ХАТХ,ОИ КА^
Донишго^и миллии Тоцикистон
Барои хдсобкунии такрибии интегралх,ои качхаттаи навъи якум формулаи квадратурии мураккабкардашудаи Симпсон сохта шуда, барои синфи функсиях,ои дифференсиронидашуда, ки ба воситаи модули бефосилагй дода мешаванд, хатогии формулах,о бах,о дода шудааст. Калима^ои калиди: интегралуои кацхаттаи навъи якум - формулаи квадратурии мураккабкардашудаи Симпсон - модули бефосилагй.
G.A.Yusupov, A.A.Shabozova EXACT ESTIMATES OF APPROXIMATE INTEGRATION OF CURVILINEAR INTEGRALS OF THE FIRST KIND ON SOME CLASSES FUNCTIONS AND CURVING LINES
Tajik National University For approximate computation of the curvilinear integrals of the first kind of complicated quadrature Simpson’s formula is constructed and for classes differentiable functions, which assigned by the modulus of continuity, in the error estimate is found.
Key words: curvilinear integrals of the first kind - the complicated quadrature Simpson’s formula - modulus of continuity.