Лабы / ЧМ.3 вариант / ЛР4 3вар
.odtМИНИСТЕРСТВО ЦИФРОВОГО РАЗВИТИЯ, СВЯЗИ И МАССОВЫХ
КОММУНИКАЦИЙ РОССИЙСКОЙ ФЕДЕРАЦИИ
Ордена Трудового Красного Знамени федеральное государственное бюджетное образовательное учреждение высшего образования
«Московский технический университет связи и информатики»
Кафедра «Информатика»
Лабораторная работа №4
«Численное интегрирование»
по дисциплине
«Численные методы»
Выполнил: студент гр. БИК2107 Мальцев В.В.
Проверил:
Москва, 2022 г.
4.1. Вопросы, подлежащие изучению
Постановка задачи численного интегрирования.
Методы прямоугольников, трапеций, Симпсона.
Оценка погрешности численного интегрирования. Правило Рунге.
Графическая иллюстрация методов прямоугольников, трапеций и Симпсона.
4.2. Задание
Выбрать индивидуальное задание из табл.3-1 для численного интегрирования:
f(x) – подынтегральную функцию;
a, b– пределы интегрирования;
методы интегрирования для выполнения п.2 – значение в столбце tиm;
начальный шаг интегрирования h0.
При этом значения в столбцах t и m означают: 1 –интегрирование методом средних прямоугольников, 2 – методом трапеций, 3 – методом Симпсона.
В сценарии пакета Scilab создать функцию для вычисления интеграла по 1-му заданному методу, определяя значения (столбец m)из табл. 3-1, с шагом и h0 и h0 / 2 (I и I).
Провести оценку погрешностей полученных результатов по правилу Рунге.
Провести оценку погрешностей полученных результатов по правилу Рунге.
Написать и выполнить программу вычисляется интеграла по 2-му заданному методу (столбец t из табл. 3-1) с точностью 10-4.
Вычислить заданный интеграл с использованием функции intg пакета Scilab.
4.3. Варианты задания
Таблица 3-1
-
№
Подынтегральная функция
a
b
t
m
1
f(x) = 8 e-x sin(-2x)
2
3
1
3
0.25
2
f(x) = e-x sin(2x)
0
2
2
1
0.5
3
f(x) = x3/2 – 2 x sin(x)
3
4
3
2
0.25
4
f(x) = e-xcos(-2x)
2
4
1
3
0.5
5
f(x) = cos(2x) + 2 sin(x)
1
3
2
1
0.5
6
f(x) = 8 sin(2x) – x
0.2
1.2
3
2
0.25
7
f(x) = 5 cos(-2x) e-x
-0.5
0.5
2
3
0.25
8
f(x) = x sin(x + 1) – cos(x – 5)
1
2
1
2
0.25
9
f(x) = 0,25 x3 + cos(x/4)
1
3
1
3
0,5
10
f(x) = sin(2x) – 2 sin(x)
3.5
5
1
3
0.5
11
f(x) = sin(ex) – e-x +1
0
1
2
1
0.25
12
f(x) = 5 x sin(x + 1) + 2 cos(x)
1
2
1
2
0.25
13
f(x) = 5 e-x + 4 x + x3/3
-1
1
1
2
0.5
14
f(x) = -2 sin(4x) ln(-x) + 5
-2.5
-1.5
1
3
0.25
15
f(x) = sin(x – 1) – x cos(x + 3)
-4
-2
3
1
0.5
16
f(x) = 4 sin (x) – x1/2
1
2
2
3
0.25
17
f(x) = 5 sin3(x) + cos3(x)
1
2
2
1
0.25
18
f(x) = cos(2x + 1) ln (2 / x) + 3
1
3
3
2
0.5
19
f(x) = 3 cos(x2) / ln(x + 5)
-1
1
1
3
0.5
20
f(x) = sin(x2) + 1 / (2 – x)
-1.5
0.5
2
1
0.5
21
f(x) = x sin(x) + cos(x) + 5
0
2
1
2
0.5
22
f(x) = – cos(x) – cos(2x) – x + 5
1
3
3
1
0.5
23
f(x) = 1 + sin(4x) / ln(x)
1.5
2.5
1
3
0.25
24
f(x) = (1 + x2)1/2 + e-x
-1
2
2
1
0.75
25
f(x) = sin(x + 1) e2 / x
1
2
3
2
0.25
26
f(x) = 2 (1 + x) e-x – 2 cos(x)
1
4
2
3
0.75
27
f(x) = – 8 sin(– x3) e-x
0.4
1.4
1
3
0.25
28
f(x) = – 10 sin(x3) cos(– x)
-1.4
-0.4
2
1
0.25
29
f(x) = x2cos(x + 3) – 4
3
4
3
1
0.25
30
f(x) = – cos(x – 5) e2x / 3
1
3
1
3
0.5
4.4. Содержание отчета
Индивидуальное задание.
Сценарий пакета Scilab для проведения расчета двумя заданными методами интегралов с шагом h0 и h0/2 (IH0 и IH0/2 ) и значения погрешностей по правилу Рунге.
Программа вычисления интеграла по 2-му заданному методу с точностью 10-4.
Результаты программы
Ход работы
Задания для численного интегрирования:
f(x) = x3/2 – 2*x*sin(x)– подынтегральная функция;
a=3, b=4–пределы интегрирования;
методы интегрирования – средних прямоугольников, трапеций,
начальный шаг интегрирования h0=0.25.
Вычисление интегралов с шагом h0 и h0/2 (IH0 и IH0/2 ) и оценка его погрешности по правилу Рунге
Правило Рунге применяют для вычисления погрешности путём двойного просчёта интеграла с шагами h/2 и h,при этом погрешность вычисляется по формуле R = |Ih - Ih/2| / (2p - 1) .
Полагают, что интеграл вычислен с точностьюЕ, если |R|<=E тогда I = Ih/2 + R , где I – уточненное значение интеграла, p – порядок метода.
Вычислим интеграл по формуле
трапеций и оценим погрешность интегрирования методом двойного просчёта:
1) a = 3, b = 4, h = 0.25, n = (b – a) / h + 1 = 5
xi = [3, 3.25, 3.5, 3.75, 4]
f(xi) = [4.3494324, 6.5622892, 9.003383, 11.5485537, 14.05442]
I0.25=0.25*(4.3494324+14.05442+2*(6.5622892+9.003383+11.5485537)= 18.146411675
2) a = 3, b = 4, h = h/2 = 0.125, n = (b – a) / h + 1 = 9
xi = [3, 3.125, 3.25, 3.375, 3.5, 3.625, 3.75, 3.875, 4 ]
f(xi) = [4.3494324, 5.4205724, 6.5622892, 7.7615041, 9.003383, 10.2715856 11.5485537, 12.8158325, 14.05442]
I0.125=0.125*(4.3494324+14.05442+2*(5.4205724+6.5622892+7.7615041+ 9.003383+10.2715856+11.5485537+12.8158325) = 18.158076050000002
3) R = |I0.125 - I0.25 | / 15 = 0.01166437500000228/15 = 0.000777625000000152
4.5 Вычисление определенных интегралов программой
Метод Симпсона:
from math import * def f(x): return x**1.5 + 2*x*sin(x) a = 3 b = 4 n = 8 I = f(a) + f(b) def simpson(a, I): h = (b - a) / n c=4 for i in range(1, n): x = a + i*h I = I + c*f(x) c = 6 - c I = h*I/3 return I, n+1, h simps = simpson(a, I) print(simps) a = 3 b = 4 I = f(a) + f(b) def simpson1(a, I): global n n = 2*n h = (b - a) / n c = 4 for i in range(1, n): x = a + i*h I = I + c*f(x) c = 6 - c I = h*I/3 return I, n+1, h simps1 = simpson1(a, I) print(simps1) I = simps[0] I1 = simps1[0] abc = 0 def R1(I,I1, e=0.0001): global abc if abs(I - I1)/15 <= e: abc = abs(I - I1) else: pass return abc R1 = R1(I,I1, e) print(float(R1), n//2, n)
И его результат
(4.0579724717112535, 9, 0.125)
(4.057966472243474, 17, 0.0625)
5.999467779460588 e-06, 8, 16
Вывод:
В данной работе я сделал операцию численное интегрирование. На компьютере и на бумаге