Об оптимальных квадратурных формулах для приближённого вычисления криволинейных интегралов первого рода на некоторых классах функций и кривых Текст научной статьи по специальности «Математика»

ДОКЛАДЫ АКАДЕМИИ НАУК РЕСПУБЛИКИ ТАДЖИКИСТАН ________________________________________2012, том 55, №12______________________________________

МАТЕМАТИКА

УДК 517.5

Л.Г.Файзмамадова ОБ ОПТИМАЛЬНЫХ КВАДРАТУРНЫХ ФОРМУЛАХ ДЛЯ ПРИБЛИЖЁННОГО ВЫЧИСЛЕНИЯ КРИВОЛИНЕЙНЫХ ИНТЕГРАЛОВ ПЕРВОГО РОДА НА НЕКОТОРЫХ КЛАССАХ ФУНКЦИЙ И КРИВЫХ

Горно-металлургический институт Таджикистана

(Представлено академиком АН Республики Таджикистан М.Ш.Шабозовым 30.10.2012 г.)

В работе рассматривается задача отыскания наилучших квадратурных формул для приближённого вычисления криволинейных интегралов первого рода для некоторых классов функций и классов кривых.

Ключевые слова: определённый интеграл - погрешность - верхняя грань погрешности - формула трапеций.

1.Пусть функция f (M) = f (X, у) определена и интегрируема вдоль кривой R2 и

J(f ;Г) :=J f (M)ds = J f (x, y)ds. (1)

г г

Пусть на кривой Г установлено положительное направление так, что положение точки

M = M (X, у) еГ может быть определено длиною дуги s = AM, отсчитываемой от начальной точки

A . Тогда кривая Г параметрически выразится уравнениями

x = x(s), у = y(s), 0 < s < L,

а функция f (M) , заданная в точках кривой, сведётся к сложной функции f (x(s), y(s)) от переменной s . В этом случае криволинейный интеграл (1) запишется в виде следующего определённого интеграла

L

J (f; Г) = j f (x(s), y(s))ds. (2)

0

Всякая квадратурная формула вида

N

J (f; Г) * Ln (f; Г; P, S) := £ p,f (x(s) y(s)) (3)

k=1

Адрес для корреспонденции: Файзмамадова Лолазор Гадомамадовна. 735730, Республика Таджикистан, Согдийская область, г. Чкаловск, ул. Московская, 6, Горно-металлургический институт Таджикистана. E-mail: [email protected]

для приближённого вычисления интеграла (2) задаётся векторами коэффициентов P = {рк }^=1 и S = {^ }^=1, где р, p2,pN - произвольные действительные числа, 0 < £1 < ^ <... < < L .

Через (L) обозначим класс плоских спрямляемых кривых Г с непрерывной кривизной,

расположенных в области Q = {(x, у) : x2 + у2 ; <L2} и у которых длина не более L, а через »£’№6) := »„(2) L, (К;0),1 < p < » - класс функций {/ (х^), y(s))}, для которых

/(х(0), у(0)) = 0, /'х(х(0), у(0)) = 0, /’ (х(0), у(0)) = 0 , имеющих почти всюду в области е ча-

д 2f

дх 2-1 ду1

стные производные 2_г.—-, 1 = 0,1,2 и удовлетворяющих условия

( L 1Р

v2 f (x(sX у (s))|| = I j I v2 f (x(sX y(s)) |pds < к, 1 < p < да

Lp [0,L]

Р V 0 У

II'v2/(x(s), у(s)^ = supvrai {V2/(x(s), j(s))| : (x, у) е q} < K, p = да. Здесь, как обычно, оператор "V " определяется равенством

V = I- + !-■■dУ. а V2-. f (x(s), у^)) := V( Vf (x(s), у(*))) =

dx ds ду ds

д2f (dx Y 2 д2f dx 2 д2f (2 '-2 '

дx

ds

1 ds) дxдy ds ds ду2

Для каждой функции f е Wq2J(K; Q) и каждой кривой Ге^е (L) остаток квадратурной формулы (3) имеет вполне определённое значение

| Rn (f; Г; P, S) |=| J (f; Г) - Ln (f; Г; P, S) |.

За величину, характеризующую точность квадратурной формулы для всех функций f е W^ (K; Q) , заданных на кривой Ге Ng (L), примем число [1,2]

Rn (W0(2)(K; ,;Г; P, S) = sup {| Rn (/;Г; P, S) |: f е W0(2)(K;Q)}.

Наибольшую погрешность квадратурной формулы (3) всего класса функций Wq2J(K;Q) на классе кривых (L) обозначим

Rn (w0, (K; Q);q (L); P, S) = sup {Rn (W0(2) (2; Q); Г; 2, S) : Г е Nq (L)}.

Для получения квадратурной формулы, которую можно было бы считать оптимальной для всех функций f є Wq2 (K; Q) и всех кривых Ге Шд (L), полагаем, что соотношение (3) является точным

N

для функций f (x(s), y(s)) = c = const, что приводит к выполнению равенства J ds = S рк = L .

Г к=1

Задача состоит в отыскании величины [3]

£n(W0Q);,(L)) = inf,(,,;Q);,(L);P,S): (,S) є A}, (4)

где A - множество всех векторов (P, S), для которых квадратурная формула (3) имеет смысл. Если существуют P0 = {p°}N=1 и S0 = {s°}N=1 такие, для которых

£» (W0; Q); , >))=Rn (W0; Q); «q ( l); p” > S”),

то квадратурная формула с этим вектором коэффициентов и вектором узлов называется наилучшей.

Используя определение класса Wq2J(K;Q) для произвольной функции, принадлежащей этому классу, получим представление

L

f (x(s)> y(s)) = J (s -tX V2 f (x(t)> y(t ))dt> (5)

0

где (s — t)+ = max(s — t, 0). Подставляя функцию (5) в квадратурную формулу (3), для погрешности получим формулу

L

Rn(f ;Г) := Rn(f; Г; P, S) = Jv2f (x(s), y(s))Ф2(s)ds, (6)

0

где

1 N

Ф2 (S) = “ (L — S)2 — S Pk (Sk — S)+ .

2 к=1

Из равенства (6) с учётом неравенства Гльдера имеем:

L

R(f ;Г; P, S)| < J| Vіf (x(s),y(s)) | • І Ф2(s) ds <

0

L

V1P f L V1 q

<

J! V2f (x(s), y(s))|Pds •( J^2(s)i«ds

ч0

'L V/q

<

<

K J^2(s)|'ds

, р 1 + q 1 = 1, 1 < q < да. (7)

Непосредственными вычислениями легко проверить, что для кривой Г* ^ (Ь), заданной пара-

метрическими уравнениями

/ N ^ / N 4

= ^2’ >,(5) ^^2’ 0 ~ ^ “ Ь’

и определенной на кривой Г*, функция

Ф)ґ і

У (*)ґ і

о V о

о V о

где

ч-1/р

Ф) = К /|Ф2(^)| Ш ■ |ф 2(м)| Ч—^ф 2 (и),

1 N ^

Ф2(и) = ~ (^ —'12и) — 2 рк (я. — '^2м)+, ф2( —^) = ф2(я),

2 ¿=1 >/2

в неравенстве (7) достигается равенство. В самом деле, поскольку

(Ь Т1/Р

V2/о(х(^)’у^)) = КI || Ф2(4> № ■ I ф2(^) Г_1 41ф2(5)’

то

\1/Ч

Ду (/о; Г; Р, X) = | V/х(4), у(4)) ^(4)* = XI || Ф2 (4) \’Ж

О V 0 ,

Таким образом, правая часть неравенства (7) является точной верхней гранью погрешности

квадратурной формулы (3) на классах функций (К;0) и кривых (Ь):

Лу «р>(К; Є); N (і); Р, 5) = КI }| «ЗД |ЧЛ

\1/Ч

(8)

(р + ч = 1,1 ^ Ч ^да)-

Полагая “к = рк /і, &к = ^ /і, к = 1,2,...,N , перепишем функцию Ф2(я) в следующем виде

N

= і2

211—7 1—2“. к* — т

к=1

:= і2 / і).

(9)

Учитывая (9), из правой части равенства (8) получаем

/ 1

Лу0С(К;0; Пд(і);Р, 5) = Кіі+1|| ф**(*) | Ш

\1/Ч

(р + ч = 1, 1 ^ Ч ^да).

Отсюда, согласно равенству (4), имеем

Є,

N

«Ц (К; 0; N (і)) = КІг'1 іпґ [{ |ф*(я) |Ш

“к ,кк^

1/Ч

В работе [4], в частности, доказано, что величина

\1/Ч

принимает минимум для значений

2 N + 7^4(1)

1 + ^^2,ч (1)

, к = 1,2,..., N — 1,

“у =

2N

2к

2N + 7Л2,ч (1) ’

к = 1,2,..., N,

(10)

(11)

где Я (?) - многочлен вида Ї2 + й + Ь, наименее уклоняющийся от нуля в метрике

АД—1,1], 1 ^ ^.

При этом

В„ = 2—

2 N + ^27 (1) ^

Таким образом, мы приходим к следующему результату.

Теорема. Среди всех квадратурных формул вида (3) для приближённого вычисления криволинейного интеграла первого рода (2) на классе функций Щ(2)(К;б) (1 ^ р < *) и классе кривых (і) наилучшей является формула

. N

I / (М Ш = 2 р*/(М *) + Я,, (/;Г),

Г к=1

2

*

=

*

=

2

где рк=акЬ,М;= М(х(а*Ь), у(акЬ)),ак и а* определены равенствами (10) и (11), х = х(4), у = у(4) - параметрические уравнения кривой Г, Ь - её длина. Для минимальной оценки погрешности формулы (12) на указанных классах функций и кривых имеет место оценка

^2+1 КТ

£„ (Жо‘;>(К; 0; N (Ь)) = ------——2, (р~' + = 1, 1 < д <»),

2

Отметим некоторые частные случаи доказанной теоремы.

Следствие 1. Пусть р = 1(ч = ю). Тогда Я оо(/) есть многочлен Чебышева

1 2 1 1 А (ґ) = ~ сов(2агссо8 Ґ) = Ґ — —, Я2 ^(1) = —, коэффициенты и узлы имеют вид:

р*= ^ , * = 1,2,..., N—1; 2^/2N + 1

* (^ + 1)і

р,у= •

При этом

* 2уі2N + 1’

* _*г 2^2кі 1 1 _ дг

",: ^ 2^ + 1 ’ * 1’2’'..’N'

КЇ2

£» «?(К ;Є); N (і)) =

(2,/2N +1)2'

Следствие 2. Пусть р = 2(ч = 2). Тогда Я 2(ґ ) есть многочлен Лежандра

2 1 2

/2 (ґ) = х ——, Я 2 (1) = —, а коэффициенты и узлы имеют вид:

* 2>/зі , 1 0 ЛГ 1

р* =2/3^^, к=1,2,...,N—1;

* _ (^3+-4/2)і _

^ = 2^/зN + >/2 ’

* * 2>/3к

"к''=акі = 2/зNГW2 , * =1 2,..., N,

и тогда

3Кі5/2

^ (ЮК; ®; % (і» = +^г.

Следствие 3. Пусть р = да(ч = 1). В этом случае Яг(ї) есть многочлен Чебышва второго

рода Q2(t) =

1

3

= Ґ---, Я1 (1) = — и коэффициенты и узлы имеют вид:

При этом

р.=

4і

4 N + >/з ’

к = 1,2,..., N — 1;

* (2 + л/3)і

р = 4ЛГ + ;

* * Т 4кі

^ :=а*і=тт;—к

4 N + Ы 3

, к = 1,2,..., N.

Ъ (Ок; Q); N (і)) =

(4N + л/з)2

ЛИТЕРАТУРА

1. Вакарчук С.Б. - Укр. матем. журн., 1986, т.38, №5, с.643-645 с.

2. Шабозов М.Ш., Мирпоччоев Ф.М. - ДАН РТ, 2010, т.53, №6, с.415-419.

3. Никольский С.М. Квадратурные формулы. - М.: Наука, 1988, 256 с.

4. Аксень М.Б., Турецкий А.Х. - ДАН СССР, 1966, т.166, №5, с.1019-1021.

Поступило 01.11.2012 г.

Л.Г.Файзмамадова

ДАР БОРАИ ФОРМУЛАМИ КВАДРАТУРИИ ОПТИМАЛИИ ТАЦРИБИ Х,ИСОБ НАМУДАНИ ИНТЕГРАЛ^ОИ КА^ХАТТАИ НАВЪИ ЯКУМ БАРОИ БАЪЗЕ СИНФ^ОИ ФУНКСИЯ^О ВА ХАТХ,ОИ КА^

Донишкадаи ку^й-металлургии Тоцикистон

Дар макола масъалаи кофтукови формулами квадратурии бехтарини такрибан х,исоб намудани интегралх,ои качхаттаи навъи якум барои баъзе синфи функсиях,о ва хатх,ои кач х,ал карда шудааст.

Калима^ои калидй: интеграли муайян - хатоги - уудуди болоии хатоги - формулаи трапесияуо.

L.G.F ayzmamadova

ON OPTIMAL QUADRATURE FORMULAS FOR APPROXIMATION INTEGRATION OF CURVILINEAR INTEGRALS OF FIRST KIND FOR SOME CLASSES OF FUNCTION AND CURVES

The Institute of Mining and Smelting of Tajikistan In this paper is considered a problem of finding the best quadrature formula for approximate calculation of curvilinear integrals of first kind for some classes of functions and curves.

Key words: definite integral - error - upper bound error - the trapezium formula.