
- •Иркутского государственного технического университета
- •2202 “Автоматизированные системы обработки информации и управления ”
- •Раздел 1 Погрешности вычислений
- •Причины возникновения погрешностей
- •Лабораторная работа 1
- •Вспомогательные материалы
- •5.1 Разработка программы в MathCad
- •Раздел 2 Численные методы решения систем линейных алгебраических уравнений
- •Решение систем линейных алгебраических уравнений методом Гаусса
- •Решение систем линейных алгебраических уравнений методом итераций
- •Лабораторная работа 2 Решение систем линейных алгебраических уравнений методом Гаусса
- •Вспомогательные материалы
- •Пример программы в MathCad
- •Лабораторная работа 3 Решение систем линейных алгебраических уравнений методом итераций
- •Раздел 3 Решение нелинейных уравнений
- •Итерационный метод решения нелинейных уравнений
- •Лабораторная работа 4 Решение нелинейных уравнений
- •Раздел 4
- •1. Первая интерполяционная формула Ньютона
- •3. Оценки погрешностей интерполяционных формул Ньютона
- •4. Интерполяционная формула Лагранжа
- •6 Обратное интерполирование
- •Блок-схема программы построения кубического сплайна и построения полинома Лагранжа представлены на рис. 4.3 и рис. 4.4
- •Лабораторная работа 5 Методы интерполирования функций
- •Раздел 5
- •Квадратурные формулы интерполяционного типа (формулы Ньютона - Котеса)
- •Формула прямоугольников
- •Формула трапеций
- •Формула Симпсона
- •Оценка погрешности квадратурных формул
- •Квадратурные формулы Гаусса
- •Блок-схемы решения задачи рассмотренными методами приведены на рис. 5.4 и 5.5.
- •Лабораторная работа 6 Методы численного интегрирования
- •Раздел 6 Численное решение обыкновенных дифференциальных уравнений
Квадратурные формулы Гаусса
Поскольку формулы Ньютона - Котеса являются интерполяционными, очевидно, что они не могут успешно использоваться для получения формул высокой точности по причине неустойчивости интерполяционного процесса для многочленов высокого порядка. По этой причине обычно используются полиномы степени от нуля до трех (соответственно, формулы прямоугольников со средней точкой, трапеций, Симпсона). Вычисление с их помощью интегралов от функций, обладающих высокой степенью гладкости, например, близким к полиномам высокой степени, представляется нерациональным. В выражение для погрешности этих формул входят первая, вторая или четвертая производные. Погрешность определяется низким порядком производной при высокой степени гладкости интегрируемой функции. Этих недостатков лишены квадратуры Гаусса.
Формулировка задачи построения квадратурных формул, поставленная Гауссом, такова.
Для заданного количества точек, а именно, для (N + 1) точки, найти такое расположение узлов и такие веса ci, чтобы квадратурная формула
была точной для полиномов как можно более высокой степени
Пусть отрезок интегрирования [a,b] непрерывной функции f(x) разбит на n равных частей точками (h – шаг разбиения . Обозначим через S(x) сплайн-функцию, аппроксимирующую подынтегральную функцию f(x).
Пусть
на каждой части разбиения [xi-1
,xi]
(i=1,2,…,n)
расположено m узлов
[xi1,…,xim]
, в которых подынтегральная функция
f(x) на каждой
i-й части аппроксимируется
многочленом
степени p, х Є [xi-1,…,xi],
(i=1,2,…,n).
При этом на многочлен
накладываются 2 ограничения:
а) значения многочлена и подынтегральной функции равны в узлах интерполяции: Si(xij)=f(xij) (i=1,…,n; j=1,…,m);
б) определенный интеграл от функции на отрезке [xi-1, xi] выражается через значения подынтегральной функции f(xij) в узлах в виде их линейной комбинации
(4)
Квадратурные формулы Гаусса для выбранной степени p сплайна будут определены, если из условий а) и б) удастся найти m неизвестных коэффициентов ci и координаты m узлов xij ( j=1,…,m).
Задача решается одновременно для всех n частей разбиения отрезка [a,b] , если выразить х Є [xi-1,…,xi], (i=1,2,…,n) через переменную t Є [-1;1]:
;
Положим
;
;
Тогда
(i=1,…,n;
j=1,…,m)
и соотношение (4) перепишем в виде
(5)
Выведем
квадратурную формулу Гаусса с тремя
узлами (m=3). Для этого
необходимо определить шесть величин:
.
Функция
-многочлен степени p, общий
вид которого
(6)
Подставив
соотношение (6) и (5) и учитывая, что
(j=1,2,3,…,p), получим тождество относительно
коэффициентов ak (k=0,1,2,…,p):
Шесть неизвестных будут определены однозначно из системы шести уравнений. В общем случае степень p аппроксимирующего многочлена всегда является нечетным числом и связана с числом узлов m соотношением p=2m-1. В частности, для трех узлов имеем многочлен пятой степени.
Множители при ak в левой части тождества вычисляется так:
Приравняем
подобные выражения в левых и правых
частях тождества при одинаковых
коэффициентах
.
Получим:
(7)
Эту систему позволяет упростить следующее свойство ее решения: неизвестные ti (i=1,2,…,m) системы 2m уравнений вида
где
k=0,1,…,2m-1
являются нулями многочлена Лежандра
;
нули принадлежат интервалу (-1;1) и
расположены симметрично относительно
середины интервала. В нашем случае m=3
и
.
Находим нули многочлена из уравнения
5t3-3t=0.
Подставив корни уравнения
,
t2=0,
в (7), получим систему трех линейных
уравнений относительно переменных
:
,
,
Теперь подставим найденные значения в соотношение (4):
,
Где
,
,
.
Итак, квадратурная формула Гаусса с тремя узлами записывается в виде
Если подынтегральная функция имеет непрерывную производную шестого порядка, то для оценки погрешности формулы Гаусса с тремя узлами можно использовать неравенство:
При вычислении интеграла до достижения заданной точности ε методом двойного пересчета условие окончания вычислений имеет вид
,где
k=2m (m
– число узлов в квадратурной формуле
Гаусса). Полагают, что
с точностью ε. Формулы Гаусса обеспечивают
высокую точность уже при небольшом
количестве узлов.