Математическая энциклопедия - кубатурная формула

формула для приближенного вычисления кратных интегралов вида

Интегрирование выполняется по -множеству в евклидовом пространстве К. ф. наз. приближенное равенство

Подинтегральная функция записана в виде произведения двух функций: первая р(х).считается фиксированной для данной К. ф. и наз. весовой функцией, вторая f(x) принадлежит достаточно широкому классу функций, напр., непрерывных и таких, что интеграл I(f) существует. Сумма в правой части (1) наз. кубатурной суммой, точки наз. узлами К. ф., а числа ее коэффициентами. Обычно хотя это требование не является обязательным. Нахождение приближенного значения I(f) с помощью формулы (1) сводится к вычислению кубатурной суммы. При n=1 формула (1) и сумма в правой ее части наз. квадратурными (см. Квадратурная формула).

Пусть

мультииндекс, здесь неотрицательные целые числа;

одночлен степени от ппеременных;

число одночленов степени не выше т от п переменных; одночлены, занумерованные так, что одночлены меньшей степени имеют меньший номер, а одночлены одной и той же степени нумеруются в любом порядке, напр. лексикографическом. При указанной нумерации а среди содержатся все одночлены степени не выше т. Пусть многочлен степени т. Множество точек в комплексном пространстве С ⁿ, удовлетворяющих уравнению наз. алгебраич. гиперповерхностью порядка т.

Один из способов получения К. ф. основан на алгебраич. интерполировании. Точки выберем так, чтобы они не лежали на алгебраич. гиперповерхности порядка т или, что равносильно, чтобы матрица Вандермонда

была неособенной. Интерполяционный многочлен функции f(x).по ее значениям в в форме Лагранжа имеет вид

где многочлен влияния j-го узла: ( символ Кронекера). Умножение приближенного равенства на р(х).и интегрирование по W приводит к К. ф. (1), в к-рой N=m и

Существование интегралов (2) равносильно существованию моментов весовой функции Здесь и далее предполагается, что требуемые моменты р(х).существуют. К. ф. (1), узлы к-рой не лежат на алгебраич. гиперповерхности порядка т, их число N=m и коэффициенты определяются равенствами (2), наз. интерполяционной К. ф. Формула (1) обладает m-с в о й с т в о м, если она обращается в точное равенство, когда f(х) - любой многочлен степени не выше т; интерполяционная К. ф. обладает m-свойством. Для того чтобы К. ф. (1), обладающая m-свойством и с числом узлов была интерполяционной, необходимо и достаточно, чтобы ранг матрицы

был равен N. При n=1 это условие выполнено, так что квадратурная формула, обладающая m-свойством и с числом узлов является интерполяционной. Фактич. построение интерполяционной К. ф. сводится к выбору узлов и вычислению коэффициентов. Коэффициенты Cj можно находить из линейной алгебраич. системы

которую получают, записывая, что К. ф. (1) (при N=m) точна для одночленов степени не выше т. Матрица системы совпадает с V.

Пусть требуется построить К. ф. (1), к-рая обладает m-свойством, а число ее узлов меньше m. За счет выбора коэффициентов это сделать невозможно, поэтому неизвестными в (1) считаются не только коэффициенты, но и узлы. Таким образом, имеется N(n+1). неизвестных. Так как К. ф. должна обладать m-свойством, то это дает m уравнений

Естественно потребовать, чтобы число неизвестных совпадало с числом уравнений: N(n+1)=m.. Это равенство ориентировочно определяет число узлов искомой К. ф. Если N=m/(n+1) не целое, то полагаем означает целую часть m/(n+1) Требуемая К. ф. с указанным числом узлов не всегда существует. Если К. ф. существует, то она имеет в (n+1) раз меньше узлов, чем интерполяционная К. ф. Однако в этом случае узлы и коэффициенты определяются из нелинейной системы уравнений (3). Метод неопределенных параметров построения К. ф. состоит в том, что К. ф. ищется в таком виде, к-рый приводит к упрощению системы (3). Это возможно в тех случаях, когда W и р(х).обладают симметрией. Расположение узлов согласуется с симметрией W и р(х), при этом симметричным узлам сопоставляются одинаковые коэффициенты. Упрощение системы (3) связано с риском: система (3) может иметь решение, а упрощенная не имеет.

Пример. Пусть Требуется построить К. ф. с 7-свойством; n=2, m=M, (2, 7)=36 и число узлов равно 12. Узлы расположим следующим образом. Первую группу узлов образуют точки пересечения окружности с центром в начале координат и радиусом ас координатными осями. Вторая группа узлов состоит из точек пересечения окружности с тем же центром и радиусом bс прямыми Третья группа узлов образуется как и вторая, при этом радиус окружности обозначается с. Коэффициенты, отвечающие узлам одной и той же группы, считаются одинаковыми и равными А, В, С для узлов первой, второй и третьей групп соответственно. Такой выбор узлов и коэффициентов обеспечивает точность К. ф. для одночленов у к-рых хоть одно из чисел iи jнечетно. Чтобы К. ф. имела 7-свойство, достаточно потребовать, чтобы она была точна для

Это приводит к нелинейной системе шести уравнений относительно шести неизвестных а, Ъ, с, А, В, С. Решая ее, получаем К. ф., узлы к-рой принадлежат К% и коэффициенты положительны.

Пусть G конечная подгруппа группы ортогональных преобразований О(п).пространства оставляющих неподвижным начало координат. МножествоW и функция р(х).наз. инвариантными относительно G, если для любого Совокупность точек вида ga, где а - фиксированная точка и g пробегает все элементы группы G, наз. орбитой, содержащей а. К. ф. (1) наз. инвариантной относительно G, если W и р(х).инвариантны относительно G и совокупность ее узлов представляет собой объединение орбит, при этом узлам одной и той же орбиты сопоставляются одинаковые коэффициенты. Инвариантными относительно G множествами являются все пространство шар и сфера с центром в начале координат, а если G группа преобразований правильного многогранника Uв себя, то инвариантен и U. Таким образом, инвариантные К. ф. можно рассматривать, когда в качестве W берутся шар, сфера, куб и любой правильный многогранник, а в качестве р (г) любая инвариантная относительно G функция, напр., р(r), где Теорема 1. Чтобы инвариантная относительно G К. ф. обладала m-свойством, необходимо и достаточно, чтобы она была точна для тех многочленов степени не выше т, к-рые инвариантны относительно G (см. [5]). Метод неопределенных параметров можно определить как метод построения инвариантных К. ф., обладающих m-свойством. В приведенном выше примере в качестве G можно взять группу симметрии квадрата. Теорема 1 имеет существенное значение при построении инвариантных К. ф.

Для простых областей интегрирования, таких, как куб, симплекс, шар, сфера, и для веса р(х)=1 можно построить К. ф. кратным применением квадратурных

формул. Напр., когда куб, то с помощью квадратурной формулы Гаусса с kузлами i/ и коэффициентами А;может быть получена К. ф.