Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Болтушкин Л.С., группа 712-2, лабораторная 1.docx
Скачиваний:
6
Добавлен:
04.10.2024
Размер:
132.57 Кб
Скачать

Министерство науки и высшего образования Российской Федерации

Федеральное государственное бюджетное образовательное учреждение высшего образования

ТОМСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ СИСТЕМ УПРАВЛЕНИЯ И РАДИОЭЛЕКТРОННИКИ (ТУСУР)

Кафедра комплексной информационной безопасности электронно-вычислительных систем (КИБЭВС)

ЧИСЛЕННОЕ РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ (СЛАУ). ЧИСЛЕННОЕ РЕШЕНИЕ НЕЛИНЕЙНОГО УРАВНЕНИЯ (НУ)

Отчет по лабораторной работе №1

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

Студент гр. 712-2 ___________ Л.С. Болтушкин ___________

Руководитель Старший преподаватель

_______ __________ Е.С. Катаева

__________

Томск 2023

Введение

Целью лабораторной работы является освоение методов решения систем линейных алгебраических уравнений и решения нелинейных уравнений с одной переменной.

Ход работы

1 Решение СЛАУ

1.1 Пример системы уравнений

Перед решением своего варианта программы нужно протестировать пример.

Для тестирования использовать систему уравнений, изображенная на рисунке 1.1.

Рисунок 1.1 – Пример системы

Из каждого уравнения были выражены x1, x2, x3:

x1= (-2 * x2 – 3 * x3 + 105) / 100;

x2 = (-x1 - 3 * x3 +104) / 100;

x3 = (-x1 - 2 * x2 + 103) / 100.

На рисунке 1.2 представлена работа программы для тестировачной системы уравнений.

Рисунок 1.2 – Работа программы тестировочной системы

1.2 Индивидуальная система уравнений

На рисунке 1.3 представлена индивидуальная система уравнений.

Рисунок 1.3 – Система уравнений

Из каждого уравнения были выражены х1, х2, х3:

x1= (-6,09 * x2 – 3,44 * x3 – 8,39) / -10,25;

x2 = (-9,64 * x1 – 9,39 * x3 -1,1) / -20,01;

x3 = (-0,6 * x1 + 0,94 * x2 + 4,87) / -6,06.

На рисунке 1.4 представлена работа программы для индивидуальной системы уравнений.

Рисунок 1.4 – Работа программы индивидуальной системы

Таблица 1 – Результаты решения тестировочной системы и индивидуальной системы.

Решение

Начальная точка

Число итераций

Тестировочная система

0,99; 0,99; 1.

0;

0;

0.

3

Индивидуальная система

0,55;

-0,02;

-0,74.

-8,39;

-1,1; 4,87

8

2 Решение нелинейного уравнения

2.1 Отделение корня

Рассмотрим функцию F(х) = 0,33x2 + 0,98x + 0,68 и ее производную

F’(x) = 0,66x + 0,98. Найдем точки перегиба, то есть точки, где производная функции F’(x) = 0. Точка перегиба x = -1,48.

Для x < -1,48 значение производной отрицательно, а для x > -1,48 значение производной положительно. Значит, функция монотонно возрастает в интервале (-∞, -1,48) и монотонно убывает в интервале (-1,48, +∞). Таким образом, отрезок [-∞, -1,48] является отрезком, на котором функция монотонна и имеет ровно один корень. График функции изображен на рисунке 2.1.

Рисунок 2.1 – График индивидуальной функции

2.2 Нахождение корня методом Ньютона

В данном методе сначала приводят касательную к функции в отрезке, выбранном пользователем, затем находится точка пересечения касательной с осью ОХ и эта точка будет первым приближением для искомого корня. Затем каждому последующему Х присваивается разница между Х на предыдущей итерации и делением функции от предыдущего Х на первую производную. Условие остановки определяется значением функции от Х на текущей итерации, если оно меньше Е, которое равно 0,01.

На рисунке 2.2 представлено решение методом Ньютона тестировочного уравнения, а на рисунке 2.3 решение тем же методом, только для индивидуальной функции.

Рисунок 2.2 – Решение тестировочного уравнения

Рисунок 2.3 – Решение индивидуального уравнения

2.3 Нахождение корня методом простых итераций

В данном методе Х присваивается сумма Х, на предыдущей итерации, и произведения функции Х, на предыдущей итерации и L, которое ровно -2 деленому на первую производную от функции, Х на прошлой итерации. Программа останавливается, тогда, когда модуль разности Х будет меньше заданной точности.

На рисунке 2.4 представлено решение тестировочного уравнения с помощью метода простых итераций, на рисунке 2.5 решение тем же методом, только для индивидуальной функции.

Рисунок 2.4 – Решение тестировочного уравнения

Рисунок 2.5 – Решение индивидуального уравнения

Таблица 2 – Результаты по тестировочному уравнению

Результаты Метод

Начальная точка

Число итераций, за которые найден корень

Найденный корень

Метод Ньютона

0,5

3

1,00

Метод простых итераций

0,5

5

1,00

Таблица 3 – Результаты по индивидуальному уравнению

Результаты Метод

Начальная точка

Число итераций, за которые найден корень

Найденный корень

Метод Ньютона

0

5

-1,1049

Метод простых итераций

0

33

-1,1113