- •Методические указания
- •Часть 1
- •1. Численное решение нелинейных уравнений.
- •1.1. Метод деления отрезка пополам.
- •1 Input a, b, e
- •1.2. Метод Ньютона (метод касательных).
- •Input X, e
- •1.3. Метод простой итерации.
- •Input X, e, m
- •2. Методы решения систем линейных
- •2.1. Метод Гаусса.
- •2.2. Метод прогонки.
- •2.3. Метод простой итерации (метод Якоби).
- •2.4. Метод Зейделя.
- •3. Численные методы решения систем
- •3.1. Метод простой итерации (метод Якоби).
- •3.2. Метод Зейделя.
- •Input X,y, m1,m2
- •Input tt
- •3.3. Метод Ньютона.
- •Input X, y
- •Input tt
- •4.1. Приближение функции по методу наименьших
- •I XI xi2 yi xiyi
- •I XI xi2 xi3 xi4 yi xiyi xi2yi
- •4.2. Интерполяционный полином
- •4.3. Интерполяционный полином
- •Литература
I XI xi2 xi3 xi4 yi xiyi xi2yi
1 -2 4 -8 16 6 -12 24
2 -1 1 -1 1 2 -2 2
3 0 0 0 0 -1 0 0
4 1 1 1 1 -2 -2 -2
5 2 4 8 16 -1 -2 -4
0 10 0 34 4 -18 20
Тогда система линейных уравнений (4.8) относительно значений a, b, c примет вид:
34 a+ 0 b+10 c= 20
0 a+10 b+ 0 c= -18 (4.9)
10 a+ 0 b+ 5 c = 4
Решая систему (4.9) получим следующие значения параметров a=0.857; b=-1.8; c=-0.914. Таким образом, искомый полином имеет вид:
y=0.857 x2 - 1.8 x - 0.914.
Таблица 4.3
i xi yi yi (yi - yi)2
1 -2 6 6.114 0.012
2 -1 2 1.743 0.066
3 0 -1 -0.914 0.007
4 1 -2 -1.857 0.020
5 2 -1 -1.086 0.007
По таблице 4.3 можно определить сумму квадратов отклонений экспериментальных данных yi от расчетных yi S=0.112.
Каким образом можно уменьшить погрешность S ? Для этого необходимо увеличить число коэффициентов полинома y. Это приведет к увеличению размерности системы. Таким образом, стараясь улучшить расчеты, т.е. уменьшить погрешность, приходится увеличивать объем вычислений.
4.2. Интерполяционный полином
в форме Лагранжа.
При решении дифференциальных, интегральных уравнений численными методами вместо искомой функции, обычно, определяют ее значения в узлах. На следующем этапе проводят интерполирование функций, т.е. восстановление функции по заданным узлам, замену графически заданной функции аналитической. Интерполяция интересует нас главным образом как один из способов построения многочлена, приближающего функцию на данном отрезке.
Пусть на некотором промежутке [a,b] заданы n различных узлов x1,x2,x3, ..., xn, а также значения некоторой функции y1,y2,y3, ..., yn в этих узлах. Необходимо построить полином P(x), проходящий через заданные точки, т.е.
P(xi)=yi
Этот полином называется интерполяционным полиномом, является единственным полиномом степени не выше n-1, и может быть записан, например, в форме Лагранжа или Ньютона.
Интерполяционный полином Лагранжа имеет следующую формулу:
P(x)= Ln-1(x)= i li(x), (4.10)
(x-x1)...(x-xi-1)(x-xi+1)...(x-xn)
где li(x) =
(xi-x1)...(xi-xi-1)(xi-xi+1)...(xi-xn)
фундаментальные полиномы Лагранжа. Они удовлетворяют равенствам
1, если i = k
lk(xi) = (4.11)
0, если i k ,
и зависят лишь от заданных узлов xi , но не от интерполируемой функции yi .
Пример. Пусть задана таблица 4.4
Таблица 4.4
x i -1 0 1/2 1
yi 0 2 9/8 0
Необходимо построить интерполяционный полином Лагранжа, проходящий через заданные точки
Ln-1(xi)=yi , i=1,2,...,n.
Решение. Запишем фундаментальные полиномы Лагранжа:
(x-x2)(x-x3)(x-x4) (x-0)(x-1/2)(x-1)
l1(x)= = = - (x3-3/2 x2+1/2 x)/3,
(x1-x2)(x1-x3)(x1-x4) (-1-0)(-1-1/2)(-1-1)
(x-x1)(x-x3)(x-x4) (x+1)(x-1/2)(x-1)
l2(x)= = = 2(x3-1/2 x2- x+1/2),
(x2-x1)(x2-x3)(x2-x4) (0+1)(0-1/2)(0-1)
(x-x1)(x-x2)(x-x4) (x+1)(x-0)(x-1)
l3(x)= = = -8(x3- x)/3,
(x3-x1)(x3-x2)(x3-x4) (1/2+1)(1/2-0)(1/2-1)
(x-x1)(x-x2)(x-x3) (x+1)(x-0)(x-1/2)
l4(x)= = = 1(x3+1/2 x2- 1/2x).
(x4-x1)(x4-x2)(x4-x3) (1+1)(1-0)(1-1/2)
Например, для l4(x) можно проверить свойство (5.2).
l4(-1)=0, l4(0)=0, l4(1.2)=0, l4(1)=1.
Подставляя li(x) в полином Лагранжа находим:
L3(x)=y1 l1(x)+y2 l2(x)+y3 l3(x)+y4 l4(x) =
= 0 l1(x)+2 l2(x)+9/8 l3(x)+0 l4(x) = x3 - 2 x2 - x + 2.