
- •Лабораторная работа №1 Вычисление погрешностей величин и результатов арифметических действий
- •Ход работы
- •24,5643 - В широком смысле
- •Лабораторная работа №2 Приближенное решение уравнений методом деления пополам и методом простой итерации
- •Ход работы:
- •Лабораторная работа №3 Приближенное решение уравнений методом хорд и касательных
- •Ход работы:
- •Лабораторная работа №4 Решение слау методом Гаусса
- •Ход работы:
- •Лабораторная работа №5 Решение слау методом итераций
- •Ход работы:
- •Лабораторная работа №6 Составление интерполяционных формул Лагранжа и Ньютона
- •Ход работы:
- •Лабораторная работа №7 Вычисление интегралов по формулам Ньютона - Котеса
- •Ход работы:
- •Лабораторная работа № 8 Реализация формул Ньютона-Котеса на эвм
- •Ход работы:
Лабораторная работа №6 Составление интерполяционных формул Лагранжа и Ньютона
Задание 1:
А) Построить правильную таблицу разностей заданной функции f(x) = р(х) + q(x) на участке [0,15; 0,25] с шагом h = 0,01. Значения f(x) вычислить с помощью ЭВМ и округлить до 4 знаков после запятой.
б) С помощью формул Ньютона вычислить значения функции f(x) в двух заданных точках c и d. Провести апостериорную оценку точности вычислений.
g(t) = tg 2x;
p(x) = e2x;
d = 0,238;
c = 0,155;
Задание 2:
С помощью таблицы, содержащей 5 узлов, задана функция у = f(x). Требуется:
а) построить для нее интерполяционный полином Лагранжа (найти его коэффициенты):
б) построить кубический сплайн с непрерывной второй производной;
в) нарисовать графики полученных полинома Лагранжа и сплайна.
Из верхней таблицы выбрать значения f(x) в 5 узлах. Номера узлов определить из нижней таблицы по номеру варианта.
N узла |
x |
y |
1 |
0,2 |
0,6 |
4 |
0,65 |
0,45 |
5 |
1,1 |
0,3 |
8 |
1,9 |
0,6 |
10 |
2,35 |
0,45 |
Ход работы:
Задание 1:
х |
y |
∆y |
∆2y |
0,15 |
2,1315 |
0,0964 |
0,0033 |
0,16 |
2,2279 |
0,0997 |
0,0035 |
0,17 |
2,3276 |
0,1032 |
0,0036 |
0,18 |
2,4309 |
0,1069 |
0,0038 |
0,19 |
2,5377 |
0,1106 |
0,0039 |
0,2 |
2,6483 |
0,1146 |
0,0041 |
0,21 |
2,7629 |
0,1187 |
0,0043 |
0,22 |
2,8817 |
0,1231 |
0,0045 |
0,23 |
|
|
|
3,0047 |
0,1276 |
0,0047 |
|
0,24 |
3,1323 |
0,1323 |
|
0,25 |
3,2646 |
|
|
б)
Если х < x0 => t*(x-x0)/h < 0
h=0,01
x=0,155
x0=0,15
t = (x-x0)/h
t = (0,155-0,15)/0,01=0,5
Pn(x) = y0+t* Δy0+t*(t-1)/2!* Δ2y0
Pn(x) = 2,1315+0,5*0,0964+0,5(0,5-1)/2!* 0,0033=2,1792;
Если х > xn => t*(x-xn)/h > 0
h=0,01
x=0,238
xn=0,25
t = (x-xn)/h
t = (0,238-0,25)/0,01= -1,2;
Pn(x) = yn + t * Δyn-1 + t*(t+1)/2!* Δ2n-2
Pn(x) = 3,2646-1,2*0,1323-1,2(-1,2+1)/2!* 0,0047 = 3,1064;
Задание 2:
x |
0,2 |
0,65 |
1,1 |
1,9 |
2,35 |
f(x) |
0,6 |
0,45 |
0,3 |
0,6 |
0,45 |
L2(x)
= 0,6
+0,45
+
+
0,3
+ 0,6
+
+
0,45
= 13,2x4+18,1х3-6,3х2+8,7х+2,4;