Ak = 2(1 - x2 k )-1 (P’n (xk ))-2
Квадратурная формула Чебышева существует при Ak = l/n, l = b - а и xk Î [a, b ] лишь для n = 1,..., 7, 9; в ней m = n - 1. Применение равных весов минимизирует вероятностную ошибку, если значения f (x ) содержат независимые случайные ошибки с одинаковой дисперсией.
При вычислении интегралов от функций с периодом l наиболее употребительны квадратурные формулы типа Гаусса:
.
Существуют квадратурные формулы для вычисления интегралов вида
где р (х ) — фиксированная, т. н. весовая функция. Её подбирают так, чтобы для всех f Î W функции f (x ) хорошо приближалась линейными комбинациями функций wq (x ).
Для приближённого вычисления неопределённых интегралов их представляют как определённые интегралы с переменным верхним пределом и далее применяют указанные выше формулы.
Таблицы узлов и весов, а также оценки погрешности квадратурных формул приводятся в специальных справочниках.
Квадратурные формулы вычисления кратных интегралов иногда называются кубатурными формулами. Кратные интегралы можно вычислять как повторные интегралы, применяя описанные квадратурные формулы. Т. к. при увеличении кратности существенно возрастает количество узлов, то для вычисления кратных интегралов разработан ряд специальных формул.
Вычисление интегралов на ЭВМ обычно осуществляется с помощью стандартных программ. В случае однократных интегралов наиболее употребительны стандартные программы с автоматическим выбором шага.
Лит.: Крылов В. И., Приближенное вычисление интегралов, 2 изд., М., 1967; Бахвалов Н. С., Численные методы, М., 1973; Никольский С. М., Квадратурные формулы, М., 1958; Березин И. С., Жидков Н. П., Методы вычисления, 3 изд., ч. 1, М., 1966; Соболев С. Л., Введение в теорию кубатурных формул, М., 1974; Коробов Н. М., Теоретикочисловые методы в приближенном анализе, М., 1963.
В. И. Лебедев.
Приближённое решение
Приближённое реше'ние дифференциальных уравнений, получение аналитических выражений (формул) или численных значений, приближающих с той или иной степенью точности искомое частное решение дифференциального уравнения.
П. р. дифференциальных уравнений в виде аналитического выражения может быть найдено методом рядов (степенных, тригонометрических и др.), методом малого параметра, последовательных приближений методом , Ритца и Галёркина методами , Чаплыгина методом . Каждый из этих методов определяет один или несколько бесконечных процессов, с помощью которых при выполнении определённых условий можно получить точное решение задачи. Для получения П. р. останавливаются на некотором шаге процесса.
Если решение ищется в виде бесконечного ряда, то за П. р. принимают конечный отрезок ряда. Например, пусть требуется найти решение дифференциального уравнения y' = f (x, у ), удовлетворяющее начальным условиям у (х0 ) = y0 , причём известно, что f (x, у ) — аналитическая функция х, у в некоторой окрестности точки (х0 , y0 ). Тогда решение можно искать в виде степенного ряда:
y (x ) - y (x0 ) = .
Коэффициенты Ak ряда могут быть найдены либо по формулам:
A1 = y’0 = f (x0 , y0 );
либо с помощью неопределенных коэффициентов метода . Метод рядов позволяет находить решение лишь при малых значениях величины х — х0 .
Часто (например, при изучении периодических движений в небесной механике и теории колебаний) встречается случай, когда уравнение состоит из членов двоякого вида: главных и второстепенных, причём второстепенные члены характеризуются наличием в них малых постоянных множителей. Обычно после отбрасывания второстепенных членов получается уравнение, допускающее точное решение. Тогда решение основного уравнения можно искать в виде ряда, первым членом которого является решение уравнения без второстепенных членов, а остальные члены ряда расположены по степеням малых постоянных величин, входящих во второстепенные члены (малых параметров). При этом уравнения для коэффициентов при степенях малых параметров линейны, что облегчает их решение. В роли малого параметра иногда выступают начальные значения (например, при изучении колебаний около положения равновесия). Метод малого параметра был использован при решении задачи о возмущённом движении в небесной механике Л. Эйлером и П. Лапласом . Теоретическое обоснование этого метода дали А. М. Ляпунов и А. Пуанкаре .
К численным методам относятся методы, позволяющие находить П. р. при некоторых значениях аргумента (т. е. получать таблицу приближённых значений искомого решения), пользуясь известными значениями решения в одной или нескольких точках. Такими методами являются, например, метод Эйлера, метод Рунге и целый ряд разностных методов.
Поясним эти методы на примере уравнения
y’ ’ = f (x, у )
с начальным условием у (х0 ) = y0 . Пусть точное решение этого уравнения представлено в некоторой окрестности точки х0 в виде ряда по степеням h = х — х0 Основной характеристикой точности формул П. р. дифференциальных уравнений является требование, чтобы первые k членов разложения в ряд по степеням h П. р. совпадали с первыми k членами разложения в ряд по степеням h точного решения.
Основная идея метода Эйлера заключается в применении метода рядов для вычисления приближённых значений решения у (х ) в точках x1 , x2 ,..., xn некоторого фиксированного отрезка [х0 , b ] Так, для того чтобы вычислить у (х1 ), где х1 = х0 + h, h = (b — x0 )/n, представляют у (х1 ) в виде конечного числа членов ряда по степеням h = х1 — х0 . Например, ограничиваясь первыми двумя членами ряда, получают для вычисления у (xk ) формулы:
,
Это т. н. метод ломаных Эйлера (на каждом отрезке [xk, xk+1 ] интегральная кривая заменяется прямолинейным отрезком — звеном ломаной Эйлера). Погрешность метода пропорциональна h2 .
В методе Рунге вместо того, чтобы отыскивать производные, находят такую комбинацию значений f (x, у ) в некоторых точках, которая даёт с определённой точностью несколько первых членов степенного ряда для точного решения уравнения. Например, правая часть формулы Рунге:
,
где
;
;
;
дает первые пять членов степенного ряда с точностью до величин порядка h5 .
В разностных формулах П. р. удаётся несколько раз использовать уже вычисленные значения правой части. Решение ищется в виде линейной комбинации у (xi ), hi и разностей Di hj , где
hj = hf (xj , yj ); Dhj = hj+1 - hj ;
Di hj = Di-1 hj+1 - Di-1 hj .
Примером разностной формулы П. р. является экстраполяционная формула Адамса. Так, формула Адамса, учитывающая «разности» 3-го порядка:
даёт решение у (х ) в точке xk с точностью до величин порядка h4 .
Для уравнений 2-го порядка можно получить формулы численного интегрирования путём двукратного применения
