Инвариантные кубатурные формулы одиннадцатой степени, содержащие оператор Лапласа

Тип работы:
Реферат
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

А. К. Пономаренко
ИНВАРИАНТНЫЕ КУБАТУРНЫЕ ФОРМУЛЫ ОДИННАДЦАТОЙ СТЕПЕНИ, СОДЕРЖАЩИЕ ОПЕРАТОР ЛАПЛАСА1
В работе рассматривается вопрос о построении инвариантных кубатурных формул
[1] с 11-свойством для интеграла по пространству Д& quot-. Формулы содержат значение оператора Лалпаса подынтегральной функции в одном из узлов (начале координат). Построение инвариантных формул связано с системами нелинейных алгебраических уравнений для координат узлов и коэффициентов, не всегда имеющими вещественные
решения, поэтому получение формул с вещественными параметрами представляет как теоретический, так и практический интерес.
Для повышения степени точности соответствующей формулы приближенного вычисления интегралов широко используется включение в кубатурную (квадратурную) сумму производных подынтегральной функции. Для интерполяционных кубатурных формул в связи с этим необходимо отметить теоремы существования И. П. Мысовских
[2], на основании которых в работах [3−5] получены содержащие производные кубатурные формулы не выше девятой степени.
В статье [6] для уточнения кубатурных формул применяется аналог теоремы С. Л. Соболева об инвариантных кубатурных формулах [1−2].
В настоящей работе при построении инвариантных относительно группы ОпО (группы всех ортогональных преобразований гипероктаэдра Оп в себя) формул вве-
П 2 П
ден оператор Лапласа X АА) двойственный с формой X хА инвариантной ато) Х) ^ сительно той же группы. (Здесь Оп -выпуклая оболочка 2п точек (±1,0,…, 0), (0, ±1,…, 0),…, (0,
0,…, ±1).) Значение лапласиана подынтегральной функции в начале координат выбрано из соображений простоты получающихся ниже нелинейных систем алгебраических уравнений. Будем рассматривать интеграл
I (/) = / Р (г)? (х) dx, х = (х1,…, Хп), г = Х& quot-=1 X.)2, р (г) -неотрицательная весовая функция такаяДЕП
& gt- 1 ½ ^ = I-*--1 '- Ю
что существуют моменты гкр (г) dr, к = 0, 1,. ., Уо & gt- 0.
о
Получим кубатурные формулы, инвариантные относительно группы ОпО и точные для любого алгебраического многочлена относительно х1,…, Хп не выше одиннадцатой степени.
Первая из рассматриваемых формул имеет вид:
2 2п
1(/) =А (О) + ВД/(О) +? А^а / (а^°))+
И 1
1Работа выполнена при финансовой поддержке РФФИ (грант № 05−01−984). © А. К. Пономаренко, 2008
Здесь использованы следующие обозначения: g л = а=Ь=(1,.. ., 1, 0 ,.. , 0), в = 0,1,., в + 1
= (Р1,№, у/1-: ЕЧ-?%, о,…, 0),
а,], Ь],] = 1, 2,
с, й, г, е — радиусы сфер с центрами в точке 0(0,…, 0), А, В, А], В],] = 1, 2, С, О, Е, Е — коэффициенты, въ в2 -параметры, Д — оператор Лапласа.
В формуле (1) суммирование распространяется на все элементы 0"С-орбит точек, указанных в круглых скобках после знака функции ?
Для нахождения параметров формул используется модификация указанной выше теоремы С. Л. Соболева, согласно которой рассматриваемые формулы должны быть точны для многочленов, а 2к, к = 0,1,…, 5, а& quot-40Т2Гс, к = 0,1, 2,3, а42, о& quot-42о"-2, '-~ек2к,
п
к = 0, 1, 2, ава4, аяа2к, к = 0, 1, аю, где а2 = У х]2, а4 = У х] 2Ак2, …, а2п =
.1 = 1)& lt-к
Х12Х22 … Хп2 -базисные инвариантные формы группы ОпО.
Это требование для формулы (1) приводит к системе 19-и нелинейных алгебраических уравнений
ДО I ___________.,. 10, _^_п /Ю, (° ~ ^)~(д ~ «) /.1. 10
128 Х& lt-'- 12. '-. '- !
у-1 10 | /.1 Ш | '- /•) до.
-Си- & lt-МчГ I -Di. il I
/ т 1и
Ь]_с —
1 л, («- 1)(и-2)(«-а),
250 1 --------24^-----------/& quot- т*||Й!
к 4,5.
(«- 1)(«- 2) (/г — 3)(н — 4)
120/г
(2)
згносительно 20 неизвестных: А, В, Ау1 — 2пА^. В$ 1 — 2п (п-1)В^,-/ - 1. 2, С — 48(. Яр, 0- Ь (. '-^Т). Е- 2пЕ, 1. 2, с, г, (. 1, с, /1|, А.
При записи системы (2) использовались обозначения:
I 2/.: 1, /г -0. 1, и — 1
I 4 _ [11уп I 2 к
2(п 2)
. 2, 'Л, 4. 5.
)Пк]!, — 1Шп | -& gt-к 1-
— 1)(& quot-
1111- (9 — I 2к 177
2)(п + 4) (и + 6)
(п — 1) (/& gt- - 2
, к — 3,4,5,
к 4,5.
'-6(я I 2)(и | 4)'-
11 — 1 (п — 2)(п2 | 1 1 77- 1 1 2)
I 1 п±& gt-2{п | 2) (}) | 4)(п | 0)(п | 8) ¦
… — 1)(& gt-г — 2)(и. — 3)
11гУц — 2к — ] I ?. ! ч г-. -
24(п I 2)(и 4)(п I («)
(п — 1)(«. — 2)(т& gt- - 3)(л — 4)
1 М ^ & quot- 12П-,) | 2)(/- I 4)(п I С)(п I 8)
2Д-Ц/2
//. — ---г-- I (^)-гахша-футтктцтя Эйлера.
\nj2) & quot-
Ступпичатая система (2) ]-к!тас-тся с частичтттлм исполтлчпваттием методики, рачраг Зоташнж и раГхп'-е |7]. Укажем ее решение,акнснщее (^г свиГх1, но|'-о параметра с-
1 20/) ' |й
п — I){п — 2 }(п — - 4) '-
т-17 — гУч
. а — с~
¦!п | ' - ¦е½ '-
у _ («~ 1^?? ~ 2^П ~ 1)1 _ 9 Г, С-П1Т ~ °У2
24 п
г/& quot-1
га*_1 V& quot- ?. -1 -- /?1 & lt-11к------Е-е-и к 2,3,4, 5,
4 2п '- …
II 1'-21& gt- 0 /т -1А: / 1 г
«А- | 1 — '-Шк ¦ и ~ ----Г^& gt--- К — 4, 5,
9
64
1 V-
(«_!)(«-а) к_
ЪтГ-
-----------Е{
} 25
ли
${ 1.1.5 -I | /'-• - I г. С 0. 5(-/ - 2 — 4г),
2
/ любой вещественный корень уравнения
t3 I 212 I (1 I z) t I у -r (l I:.
= 0.
Зл
У =
— з
7i =
(«з -'-las I 3o: s)(a7-r4 — 2сїйГ& quot- | ag)
r2 (& lt-y,'7<-y, ij —
a^ag — ag
Pi
С
О- 4 — 46
O: o — Ла: -|о :2 ад — у? г8 a- - уі& quot-ть '-
,"7 — ї/і7!'-/'-0
ifr5(a7r4 — 2agr2 + ao) ' c, J
j — 1,2, корни квадратного уравнения -и2 + рм + q — 0, «11(11−1 — «12& lt-*13 Q'-13 — «12 «її
где р
tt-v2 — 0: 110:14 в: п-і-у 1
••'--I — - 0: 12
Q-к і 9 — '-l (afc
3Clt
гі*іг
2/c
a12 — ttHQl-l J- 1−2.
fr = 2. 3. 4,5
o: 1s+, — ffl, — h: — + r, r-: — + Fvr2k + Т)^2к + /tlS2fe)
3 1
'-1 — -OMpi — jvT — 4(/i), «2 — 0. 5(-pi I yVi -¦%).
А/ij д 2
Pi — -7-, & lt-Л — -7-, A (j — «1* - Ul?"19, iif! Ao
Ap, — air"20 — & lt-*18ai9t Д91 — - C*18& lt-*20>-
А, А і л A _42
Ao-1 —
Л11
Ад,
4 4/'-'- 2
, Д, ТІ «И Г-: — «і)'-
— & lt-*is), A. 4−2 — -a^afair —
Діл '- Діл
Приведем пример параметров формулы (1) при п — 5 для весовой функции р (г) —
r-h-
А — 5. 188 488 962 699е + 00, В — 2. 578 621 491 614е — 01, «1 — 1. 284 003 969 924е + 01, Лі = 8. 949 365 779 309 г — 11, а-2 — 1. 91 684 565 7915s+ 00, Л-2 — 2. 58 626 830 9151s — 02,
с — 4. 659 402 407 556е 4- 00,
С — -2. 66 247 906 2243t — 06, г — 2. 44 948 974 2783ft 4- 00,
F = 3. 337 598 865 408с- 01, d — 3. 632 864 387 059е4- 00,
Т) — 6. 72 422 899 5560ft — 05,
61 — 1. 6 109 481 9922(к 4−00,
В г — 1. 38 427С90СС25е-01, b2 — 2. 44 849 990 0388ft 4- 00, В-2 -3. 99 952 791 51 G2ft + 00,
t — 1. 79 000 000 0000ft+ 00, Е — 6. 82 668 111 868е — 02.
Я
5. 11 669 285 3074ft — 01.
/3-і 4. 87 843 089 4815ft- 01.
fi. i:
Рассмотрим далее формулу, получаемую из (1) в результате замены последнего слагаемого в
правой части на ЕА 12Сп? (ед (4)).
Решение соответствующей, аналогичной (2), системы, зависящее от параметра е, отличается от
решения системы (2) лишь выражениями
г —
d- - & lt-¦ -
а 4 —
т-1-
«Мее'--* - 25m 1 я
I Cl. I & gt-!••!.
n: t — 4 as | За»
і де обозначено
З
a*--1 — VT)-. +4 ah | і = mk (¦4-і
D-| dlk
¦K
: 2k
6 — - 01((1к — -Eic 64 25
2k
25 3
128'-
— 2, 3. 4. 5, к = 4, 5, к — Л. 4, 5,
Здесь el = mclE.
І Іриведем гужмер параметров последней формулы при п — 6 п случае весопой функции р (г) — с_г:
А — 7. 108 611 484 509с 4- 00, В = 3. 884 616 869 089с — 01, a-i — 7. 678 622 224 934с 4- 00, А: — 0. 848 729 338 124с — 08, «2 — 2. 147 406 607 715с 4- 00, Л-2 = 5. 116 321 188 763с — 02, bi — 1. 640 356 801 648с 4- 00, R-i — 1. 815 163 736 702с — 01, b-& gt- = 2. 543 851 025 693с I 00, Во — -2. 112 305 855 193с-02.
с — 2. 581 981 447 715с + 00,
С = -6. 173 963 472 088с — 02, г — 2. 64 575 131 1065с+ 00, F — 14 51 646 899 478с-02, d — 2. 34 788 722 766с + 00,
В = 5. 96 668 578 653с — 02, с — 3. 110 000 000 000с+ 00,
Е — 1. 117 881 561 393с-03, в = 4. 851 330 019 827с — 01, ЗІ - 3. 135 831 959 235с-01.
(1. 2)
Как показали вычисления, формула (1} имеет вещественные координаты узлов и коэффициенты при п = 5, 7, 8,…, 11, вторая формула — при п = 5,6 (при п = 5 фо]& gt-
'-г'-Ч] - 0 & lt- г & lt- 1.
о, 1 г.
мулы совпадают) дня весовых функции р (г) — гас г. р (г) —
и соответствующих значении е.
Следует заметить, что хотя формулы (1. 1) и (1. 2) имеют по одному отрицательному коэффициенту, вычисление всех моментов до 11-го порядка включительно показало совпадение 14-и первых значащих цифр и полученных точных значений (использовался пакет BC31, переменные, характеризующие промежуточные и окончательные результаты, имели описание double). Это свидетельствует о достаточно хорошей устойчивости формул.
Число узлов формулы (1) равно
2
Ni = -n (n3 + 8n2 — 25n + 22) + 2™ + 2,
второй формулы----
N2 = -n (2n4 — 15n3 + IlOn2 — 225n + 158) + 2. 15
В аналогичных формулах автора (например, [9]), не содержащих оператор Лапласа, число узлов больше на 2(n — 1).
Для сравнения количества узлов полученных и других формул используем теорему о нижней границе числа узлов кубатурной формулы в случае центральной симметрии области интегрирования [2], формулу — сферическое произведение формулы типа Гаусса с шестью узлами для интеграла по радиусу и формулы с 11-свойством для сферы [8] (с минимальным числом узлов, кроме n = 8, 9) и формулу (6. 26) из [2]. 2Число узлов рассматриваемой формулы типа сферического произведения
N3 = -n (2n4 — 15n3 + 70n2 — 105n + 63), 5
формулы (6. 26) и формулы — декартова произведения формул типа Гаусса с 6-ю узлами по каждой из координат x1, Х2, • •
•, xn равно
N4 = 6n.
Нижняя граница числа узлов, полученных формул по теореме Мысовских, имеет вид Nmin = jAn-(& quot--4 + IOn3 + 55n2 + 110n + 184).
Следует отметить, что эта нижняя граница практически не достигается.
Как показали вычисления, отношения А1, А2 при n = 5,6,---, 11 заключены между 2. 32 и 2. 70, А между 9. 96 и
22. 26, а-между 23. 35 и 54 960.9.
Summary
A. K. Ponomarenko. Some invariant cubature formulae of the eleventh degree containing Laplasian operator.
For integral over the n-space with a radial symmetric weight function two cubature formulae of degree eleven invariant on the group of the hyperoctahedron are constructed. Cubature summas formulas contain the value of the Laplasian operator of the integrand function in the point O (O, 0, • • •, 0). The examples of approximate values of parameters of these formulae are given.
Литература
1. Соболев С. Л. О формулах механических кубатур на поверхности сферы // Сибирск. мат. журн. 1962. T. 3, № 5. С. 769−791.
2. Мысовских И. П. Интерполяционные кубатурные формулы. М., 1981. 336 с.
3. Мысовских И. П. Применение ортогональных многочленов к построению кубатурных формул // Журн. вычисл. мат. и мат. физ. 1972. Т. 12, № 2. С. 467−475.
4. Исматуллаев Г. П. О построении кубатурных формул для гипершара и поверхнсти сферы // Вопросы вычислит. и прикл. математики. Ташкент, 1978. Вып. 51. С. 191−203.
5. Пономаренко А. К. О кубатурных формулах девятой степени, содежащих производные // Методы вычислений. Л., 1991. Вып. 16. С. 30−41.
6. Шамсиев С. Ш. Об инвариантных кубатурных формулах, содержащих производные // Численное интегрирование и смежные вопросы. Ташкент: Фан, 1990. С. 77−85.
7. Лебедев В. И. О квадратурах на сфере // Журн. вычисл. мат. и мат. физ. 1976. Т. 16, № 2. С. 293−306.
8. Стоянова С. Б. Кубатурная формула одиннадцатой степени точности для сферы // Методы вычислений. Л., 1980. Вып. 12. С. 38−46.
9. Пономаренко А. К. Инвариантные кубатурные формулы одиннадцатой степени, инвариантные относительно группы гипероктаэдра // Кубатурные формулы и их приложения: VI международный семинар-совещание. Уфа: ИМВЦ УфНЦ РАН, БГПУ, 2001. С. 96−103.
Статья поступила в редакцию 13 сентября 2007 г.
2АВТору неизвестны другие формулы 11-й степени по Кп, отличные от формул типа декартового и сферического произведений и формул автора.

ПоказатьСвернуть
Заполнить форму текущей работой