Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ЧМ_4

.docx
Скачиваний:
1
Добавлен:
07.06.2024
Размер:
123.13 Кб
Скачать

МИНИСТЕРСТВО ЦИФРОВОГО РАЗВИТИЯ, СВЯЗИ И МАССОВЫХ КОММУНИКАЦИЙ РОССИЙСКОЙ ФЕДЕРАЦИИ

Ордена Трудового Красного Знамени федеральное государственное бюджетное образовательное учреждение высшего образования

«Московский технический университет связи и информатики»

Факультет Сети и системы связи

Кафедра Информатики

Лабораторная работа № 4

по дисциплине «Численные методы»:

  

«Численное интегрирование»

Подготовил:

студент 2-го курса,

17 вариант

Проверил:

Юсков И.О.

Москва

2024

Оглавление

Индивидуальное задание к работе 2

Решение интеграла с помощью С++ 4

Индивидуальное задание к работе

Рисунок 1 Индивидуальное задание к работе

Где в столбцах t и m. 1 – интегрирование методом средних прямоугольников, 2 – метод трапеций; a и b – пределы интегрирования; f(x) = 5sin3(x)+cos3(x) – подынтегральная функция; h0 – начальный шаг интегрирования

Выполнение задания:

Интегрирование методом средних прямоугольников (m = 1)

Вычисление интегралов с шагом h0 = 0.25 и h0/2 = 0.125 (Ih0 и Ih0/2) и оценка его погрешности по правилу Рунге

Правило Рунге применяют для вычисления погрешности путём двойного просчёта интеграла с шагами h/2 и h, при этом погрешность вычисляется по формуле

Полагают, что интеграл вычислен с точностью E, если |R|<E тогда I=Ih/2+R, где I – уточненное значение интеграла, p – порядок метода.

Вычислим интеграл по формуле средних прямоугольников и оценим погрешность интегрирования методом двойного просчёта:

Решение интеграла с помощью программы Scilab:

deff('y=f(x)','y=(5*(sin(x)^3)+(cos(x)^3))');

a=1;b=2;

[s,ir]=intg(a,b,f)

s =

4.4150773

ir =

4.902D-14

Интегрирование методом трапеций (t = 2)

Здесь мы используем метод трапеций, вычисляя интеграл с точность 10-4

n=4

Необходимая точность

Решение интеграла с помощью С++

Вывод: результаты интегрирования подынтегральной функции на отрезке   методами трапеций ( = 12.861; =5.721405; R=-2.379865; I=3.341540), средних прямоугольников ( ; ; ) и методом, заложенным в функции intg() пакета Scilab, совпали в пределах допустимой погрешности - значит, расчёты проведены верно.

Соседние файлы в предмете Численные методы