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

КОММУНИКАЦИЙ РОССИЙСКОЙ ФЕДЕРАЦИИ

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

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

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

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

по теме «Методы решения нелинейных уравнений»

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

Выполнил: студент гр. БИК2107 Мальцев В.В.

Проверил:

Москва, 2022 г.

1.1. Вопросы, подлежащие изучению

1. Постановка задачи численного решения нелинейных уравнений.

2. Этапы численного решения уравнения.

3. Аналитический и графический методы отделения корней.

4. Уточнение корня методами половинного деления, итерации, Ньютона и хорд.

5. Графическая иллюстрация методов половинного деления, итерации, Ньютона и хорд.

6. Условие окончания вычислений при использовании методов половинного деления, итерации, Ньютона и хорд.

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

8. Теорема о сходимости метода Ньютона и оценка погрешности метода.

9. Правило выбора неподвижной точки, начальной точки и условие сходимости метода хорд.

10. Условия окончания вычислений в методах итерации, Ньютона и хорд.

11. Сравнение методов половинного деления, итерации, Ньютона и хорд.

1.2. Общее задание

1. Выбрать индивидуальное задание из табл. 1-1:

• нелинейное уравнение;

• методы решения нелинейного уравнения для выполнения 3-х итераций;

2. Отделить корни заданного уравнения графическим и аналитическим методом с использованием средств пакета Scilab.

3. Для каждого из заданных методов провести исследование функции нелинейного уравнения:

• проверить выполнение условий сходимости вычислительного процесса, в случае расходящегося процесса – сделать необходимые преобразования для обеспечения сходимости;

• выбрать начальное приближение к корню;

• сформулировать условие окончания этапа уточнения корня.

4. С использованием итерационной формуле 1-го заданного методу провести расчет трех итераций с использованием средств пакета Scilab. Результаты расчета свести в табл. 1-2.

5. Оценить погрешность результата после 3-х итераций.

6. Для 2-го заданного метода выполнить решение уравнения с точностью 10^-4, создав программу, реализующую заданный метод. Произвести расчет, а результаты решений свести в табл. 1-2.

7. Найти решение нелинейного уравнения на отделенном отрезке с использованием функции fsolve пакета Scilab.

Варианты заданий

Таблица — 1

Индивидуальное задание

1) Задание для решения уравнения:

3 sin (x^1/2) + x – 3 = 0 – уравнение

Методы, с использованием которого требуется провести 3 итерации по уточнению корня заданного уравнения. (В моем примере рассматривается метод половинного деления, метод Ньютона).

2) Отделение корней

Рисунок 1 — Графическое отделение корней

## Функция для вывода x, f(x),f'(x),f''(x) import math import numpy as np x = np.array([i for i in range(0, 12, 2)]) / 10 f = np.array([3. * math.sin(i ** 0.5) + i - 3. for i in x]) df1 = np.array([3. * math.cos(i ** 0.5)/(2 * i ** 0.5) + 1. for i in x]) df2 = np.array([3*(-1 * math.sin(i**0.5)*1/(2* i ** 0.5) * 2 * i ** 0.5 - 3*math.cos(i ** 0.5) * 2 * 1/(2*i ** 0.5))/(2*i ** 0.5) ** 2 for i in x]) x.shape = (6,1) f.shape = (6,1) df1.shape = (6,1) df2.shape = (6,1) p = np.array([x, f, df1, df2]) p.shape = (4,6) p = p.T print(f'x = __________f(x) = _______ df1(x) =_______ df2(x) =') print(p)

И его результат

x = ___f(x) = ___________df1(x) =______df2(x) =

[[ 0. -3. inf nan]

[ 0.2 -1.50263548 4.0242448 -24.30354168]

[ 0.4 -0.82661865 2.91296867 -8.28199584]

[ 0.6 -0.30171639 2.38401678 -4.33432679]

[ 0.8 0.13955195 2.04977627 -2.69944048]

[ 1. 0.52441295 1.81045346 -1.84678343]]

Таблица 2 — результаты

x	F(x)	F’(x)	F’’(x)
0	-3	+∞	∞
0.2	-1.5026355	4.0242448	-9.1823177
0.4	-0.8266186	2.9129687	-3.4995742
0.6	-0.3017164	2.3840168	-2.0276322
0.8	0.139552	2.0497763	-1.3872202
1	0.524413	1.8104535	-1.03633

Вывод: На концах отрезка (0;1] функция имеет противоположные знаки, а 1-я производная знакопостоянна, следовательно, на этом отрезке уравнение 3 sin (x^1/2) + x – 3 = 0 имеет единственный корень.