Алгебраические и трансцендентные уравнения

Август 3, 2022

Главная
Математика
Алгебраические и трансцендентные уравнения

Содержание

2. 2.1 Методы локализации корней 2.1.1 Аналитический метод Теоретической основой алгоритма отделения корней служит теорема Коши о
3. Пример 2.1. Отделить корни уравнения sin 5x + x2 – 1 = 0. Решение. Построим таблицу
4. Теорема 2.2. Если непрерывная функция f(x) монотонна на отрезке [a, b] и на его концах принимает
5. 2.1.2 Графический метод Для отделения корней уравнения зачастую бывает удобно применять графический метод. График функции y
6. рис 2.2 Из рис. 2.2 видно, что графики пересекаются в четырех точках, и данное уравнение имеет
7. 1) Если функция f(x) на концах отрезка [a, b] принимает значения разных знаков, то делим отрезок
8. рис 2.3 Рис. 2.3 показывает, что корень находится в интервале (0,44; 0,45), так как функция меняет
9. Повторив всю процедуру еще раз, заменим в диапазоне A2:A12 интервал (0,44; 0,45) на интервал (0,441; 0,442)
10. Табл. 2.2 Изложенный метод можно охарактеризовать как метод деления отрезка на 10 частей. Метод применим в
11. 2.2 Методы уточнения корней После того как найден интервал, содержащий корень, применяют методы последовательных приближений, или
12. Метод половинного деления (другие названия: метод бисекций, метод дихотомии) для решения уравнения f(x) = 0 заключается
13. Пример 2.4. Уточнить методом бисекций с точностью до 0,01 корень уравнения (x – 1)3 = 0,
14. Таблица 2.4 Приведенный алгоритм учитывает возможный случай «попадания в корень», т.е. равенство f(x) нулю на очередном
15. Создадим в программе Excel пользовательские функции f(x) и bisect(a, b, eps) для решения уравнения методом половинного
16. 2.2.2 Метод итераций Метод простых итераций для уравнения f(x) = 0 заключается в следующем: 1) Исходное
17. На рис. 2.4 показан процесс получения очередного приближения по методу итераций. Последовательность приближений сходится к корню
18. рис. 2.5 Пример 2.5. Найти корень уравнения sin 5x + x2 – 1 = 0 с
19. Все условия теоремы выполнены, мы можем применить метод итераций. Выполним вычисления в программе Excel: 1) Вводим
20. По этому критерию приближенным значением с точностью ε = 0,001 является x5 = 0,244707 ≈ 0,245.
21. 2.2.3 Метод хорд Метод хорд заключается в замене кривой y = f(x) отрезком прямой, проходящей через
22. Алгоритм метода хорд: 1) Пусть k = 0; 2) Вычислим следующий номер итерации: k = k
23. Решение. Проведем расчеты в программе Excel: 1) В ячейки A1:H1 запишем заголовки столбцов как в табл.
24. Решение в программе Mathcad: Как видим, результаты расчетов согласуются с предыдущими ответами.
25. 2.2.4 Метод Ньютона (метод касательных) Пусть найдено приближенное значение корня уравнения f(x) = 0, обозначим его
26. В точке пересечения касательной с осью OX переменная y = 0. Приравнивая y нулю, выразим x
27. Теорема 2.4. Пусть на отрезке выполняются условия: 1) функция и ее производные и непрерывны; 2) производные
28. , следовательно, функция возрастающая, но тогда последовательность является монотонной. Рис. 2.8. Достаточные условия сходимости метода Ньютона
29. Пример 2.7. Уточнить до 0,000001 методом Ньютона корень уравнения sin 5x + x2 – 1 =
30. 2.2.5 Метод секущих Метод секущих может быть получен из метода Ньютона при замене производной приближенным выражением
31. 2.3 Системы нелинейных уравнений Система n нелинейных уравнений с n неизвестными имеет вид fk(x1, x2, …,
32. Рис. 2.9. Графическое решение системы двух уравнений Пример 2.8. Решить графически систему двух уравнений Решение. Сначала
33. 3) Изменяя, т.е. уменьшая отрезок [a, b], уточняем решение (xs, ys). Решение в программе Excel. Так
34. Щелкнем правой кнопкой мыши по оси X диаграммы и выберем «формат оси», в появившемся окне выберем
35. 5) На диаграмме увидим, что графики пересекаются между значениями x = 0,6 и x = 0,7.
36. 2.3.1 Метод итераций Приведем систему (2.8) к виду, удобному для итераций xk = φk(x1, x2, …,
37. 2.3.2 Метод Ньютона Строгие формулировки теорем об условиях сходимости метода Ньютона достаточно громоздки, на практике часто
38. (2.14) Пример 2.10. Решить методом Ньютона систему нелинейных уравнений из примера 2.7 Решение. Найдем определитель матрицы
39. Теперь для данной системы метод Ньютона можно записать в виде итерационных формул: В следующей таблице приведены
40. Контрольные вопросы. 1. Как определяется понятие «корень уравнения»? 2. Какие уравнения называются алгебраическими? 3. Какие уравнения
42. Скачать презентацию

Слайд 2

2.1 Методы локализации корней 2.1.1 Аналитический метод
Теоретической основой алгоритма отделения корней

служит теорема Коши о промежуточных значениях непрерывной функции:
Теорема 2.1. Если функция f(x) непрерывна на отрезке [a, b] и f(a) = A, f(b) = B, то для любой точки C, лежащей между A и B на этом отрезке существует точка , что f(ξ) = C.
Следствие. Если функция f(x) непрерывна на отрезке [a, b] и на его концах принимает значения разных знаков, то на этом отрезке существует хотя бы один корень уравнения f(x) = 0.
Пусть область определения и непрерывности функции является конечным отрезком [a, b]. Разделим отрезок на n частей:
Вычисляя последовательно значения функции в точках a0, a1, … an, находим такие отрезки [ak, ak+1], для которых выполняется условие
(2.1)
т.е. или Эти отрезки и содержат хотя бы по одному корню.

Слайд 3

Пример 2.1. Отделить корни уравнения sin 5x + x2 – 1

= 0.
Решение. Построим таблицу значений функции y = sin 5x + x2 – 1 на отрезке [–4; 4] с шагом изменения аргумента h = 1, пользуясь калькулятором или электронными таблицами (табл. 2.1).
Табл. 2.1
Табл. 2.1 показывает, что данное уравнение имеет корни в интервалах (–1; 0) и (1; 2), так как функция меняет знак в этих промежутках. Пока мы не можем утверждать, что в найденных интервалах содержится ровно по одному корню и, что в других интервалах корней нет. Чтобы уточнить информацию о числе корней можно построить таблицу значений функции с меньшим шагом, например h = 0,1.

Слайд 4

Теорема 2.2. Если непрерывная функция f(x) монотонна на отрезке [a, b]

и на его концах принимает значения разных знаков, то на этом отрезке существует единственный корень уравнения f(x) = 0.
Если функция f(x) дифференцируема и её производная сохраняет знак на отрезке [a, b], то f(x) монотонна на этом отрезке.
Если производная легко вычисляется и нетрудно определить её корни, то для отделения корней уравнения f(x) = 0 можно применить следующий алгоритм:
1) Найти критические точки, т.е. точки, в которых производная равна
нулю или не существует, и определить интервалы знакопостоянства
производной (на этих интервалах функция f(x) может иметь только по
одному корню);
2) Составить таблицу знаков функции f(x), приравнивая переменную x
критическим и граничным значениям, или близким к ним;
3) Определить отрезки, на концах которых функция принимает значения
разных знаков.
Пример 2.2. Отделить корни уравнения sin x + x – 1 = 0.
Решение. Найдем производную функции f(x) = sin x + x – 1 и её корни:
Функция f(x) = sin x + x – 1 монотонна на отрезках [–π + 2πk, π + 2πk]. Очевидно, что лишь отрезок [– π, π] содержит корень и он единственный, так как при x> π f(x)>0; при x< –π f(x)<0.

Слайд 5

2.1.2 Графический метод
Для отделения корней уравнения зачастую бывает удобно применять графический

метод. График функции y = f(x) с учетом ее свойств дает много информации для определения числа корней и отрезков их расположения для уравнения f(x) = 0.
Рис 2.1
Из графика функции y = sin 5x + x2 – 1 на рис.2.1 видно, что на отрезке [0; 0,5] есть два корня уравнения sin 5x + x2 – 1 = 0, а на отрезках [–1; –0,5] и [1; 1,5] по одному корню. Чтобы убедиться в том, что больше корней нет, преобразуем уравнение к виду sin 5x = 1 – x2 и построим графики двух функций f1(x) = sin 5x и f2(x) = 1 – x2. Корням соответствуют абсциссы точек пересечения этих графиков.

Слайд 6

рис 2.2
Из рис. 2.2 видно, что графики пересекаются в четырех точках,

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

Слайд 7

1) Если функция f(x) на концах отрезка [a, b] принимает значения

разных знаков, то делим отрезок на 10 равных частей и находим ту часть, которая содержит корень (таким способом мы можем уменьшить длину отрезка, содержащего корень, в 10 раз).
2) Повторим действия предыдущего пункта для полученного отрезка.
Этот процесс можно продолжать до тех пор, пока длина отрезка не станет меньше заданной погрешности.
Пример 2.3. Вычислить графически с точностью до 0,0001 корень уравнения sin 5x + x2 – 1 = 0, принадлежащий интервалу (0,4; 0,5).
Решение. Построим график функции y = sin 5x + x2 – 1 на отрезке [0,4; 0,5] с шагом h = 0,01 (делим отрезок на 10 частей) в программе Excel:
1) В диапазоне A2:A12 введем значения переменной x. Для этого в ячейке A2 запишем 0,40, в ячейке A3 — значение 0,41. После этого выделим диапазон A2:A3 и с помощью маркера заполнения присвоим значения остальным ячейкам до ячейки A12.
2) В ячейку B2 введем формулу =SIN(5*A2)+A2^2–1 и скопируем B2 с помощью маркера заполнения в остальные ячейки до ячейки B12.
3) Выделим диапазон A2:B12 и с помощью мастера диаграмм (тип диаграммы “Точечная”!) построим график функции.
Лист Excel отображен на рис. 2.3.

Слайд 8

рис 2.3
Рис. 2.3 показывает, что корень находится в интервале (0,44; 0,45),

так как функция меняет знак в точках 0,44 и 0,45.
Заменим значения переменной x на том же листе в диапазоне A2:A12, то есть вместо интервала (0,4; 0,5) подставим интервал (0,44; 0,45) с шагом h = 0,001. Для этого в ячейке A2 запишем 0,440, а в ячейке A3 — значение 0,441. Затем выделим диапазон A2:A3 и с помощью маркера заполнения присвоим значения остальным ячейкам до ячейки A12. Формулы в ячейках B2:B12 не трогаем! В результате этого получим новую таблицу значений функции, из которой получаем уточненный интервал (0,441; 0,442).

Слайд 9

Повторив всю процедуру еще раз, заменим в диапазоне A2:A12 интервал (0,44;

0,45) на интервал (0,441; 0,442) с шагом h = 0,0001. Искомый корень содержится в интервале (0,4410; 0,4411). Длина этого интервала равна 0,0001 и любое число из этого интервала можно принять за приближенное значение корня с погрешностью 0,0001. Выберем середину отрезка, т.е. положим x ≈ 0,44105.
В таблице 2.2 приведены все три этапа уточнения корня. Здесь мы не приводим соответствующие графики, так как для отделения корня достаточно рассмотреть таблицу значений функции и найти последовательные значения переменной x, в которых функция изменяет знак.
Аналогично можно уточнить значения других корней данного уравнения.
Для этого достаточно на том же расчетном листе вместо отрезка [0,4; 0,5] рассмотреть любой из оставшихся трех отрезков [– 0,8; – 0,7], [0,2; 0,3], [1,1; 1,2].
Обратите внимание на то, что здесь рассматриваются начальные отрезки длиной 0,1 для того, чтобы после каждого уточнения мы получили уточненную верную десятичную цифру приближенного значения корня.

Слайд 10

Табл. 2.2
Изложенный метод можно охарактеризовать как метод деления отрезка на 10

частей. Метод применим в случае, когда левая часть уравнения f(x) = 0 задана аналитическим выражением через известные функции, непрерывна на данном отрезке, и на концах его принимает значения разных знаков. Этот метод особенно удобен для применения в электронных таблицах.

Слайд 11

2.2 Методы уточнения корней
После того как найден интервал, содержащий корень,

применяют методы последовательных приближений, или итерационные методы вычисления корня с заданной точностью.
2.2.1. Метод половинного деления
рис. Метод половинного деления

Слайд 12

Метод половинного деления (другие названия: метод бисекций, метод дихотомии) для решения

уравнения f(x) = 0 заключается в следующем. Пусть известно, что функция непрерывна и принимает на концах отрезка [a, b] значения разных знаков, тогда корень содержится в интервале (a, b). Разделим интервал на две половины и дальше будем рассматривать ту половину, на концах которой функция принимает значения разных знаков. Этот новый отрезок снова делим на две равные части и выбираем из них ту, которая содержит корень. Этот процесс продолжается до тех пор, пока длина очередного отрезка не станет меньше требуемой величины погрешности.
Более строгое изложение алгоритма метода половинного деления:
1) Вычислим x = (a + b)/2; вычислим f(x);
2) Если f(x) = 0, то переходим к пункту 5;
3) Если f(x)∙ f(a) < 0, то b = x, иначе a = x;
4) Если |b – a| > ε, переходим к пункту 1;
5) Выводим значение x;
6) Конец.

Слайд 13

Пример 2.4. Уточнить методом бисекций с точностью до 0,01 корень уравнения

(x – 1)3 = 0, принадлежащий отрезку [0,95; 1,1].
Решение в программе Excel:
1) В ячейках A1:F4 введем обозначения, начальные значения и формулы, как показано в таблице 2.3.
2) Каждую формулу скопируем в нижние ячейки маркером заполнения до десятой строки, т.е. B4 — до B10, C4 — до C10, D3 — до D10, E4 — до E10, F3 — до F10.
Таблица 2.3
Результаты расчетов приведены в табл. 2.4. В столбце F проверяем значения длины интервала b – a. Если значение меньше чем 0,01, то в данной строке найдено приближенное значение корня с заданной погрешностью. Потребовалось 5 итераций для достижения требуемой точности. Приближенное значение корня с точностью до 0,01 после округления до трех знаков равно 1,0015625 ≈ 1,00.

Слайд 14

Таблица 2.4
Приведенный алгоритм учитывает возможный случай «попадания в корень», т.е. равенство

f(x) нулю на очередном этапе. Если в примере 2.3 взять отрезок [0,9; 1,1], то на первом же шаге попадаем в корень x = 1. Действительно, запишем в ячейке B3 значение 0,9. Тогда таблица результатов примет вид 2.5 (приведены только 2 итерации).
Таблица 2.5

Слайд 15

Создадим в программе Excel пользовательские функции f(x) и bisect(a, b, eps)

для решения уравнения методом половинного деления, пользуясь встроенным языком Visual Basic. Их описания приведены ниже:
Function f(x)
f = (x - 1) ^ 3
End Function
Function bisect(a, b, eps)
1 x = (a + b) / 2
If f(x) = 0 Then GoTo 5
If f(x) * f(a) < 0 Then
b = x
Else
a = x
End If
If Abs(a - b) > eps Then GoTo 1
5 bisect = x
End Function
Функция f(x) определяет левую часть уравнения, а функция bisect(a, b, eps) вычисляет методом половинного деления корень уравнения f(x) = 0. Обратим внимание на то, что в функции bisect(a, b, eps) используется обращение к функции f(x).

Слайд 16

2.2.2 Метод итераций
Метод простых итераций для уравнения f(x) = 0

заключается в следующем:
1) Исходное уравнение преобразуют к виду, удобному для итераций:
x = φ(x); (2.2)
2) Выбирают начальное приближение x0 и вычисляют последующие приближения по итерационной формуле
xk = φ(xk-1), k = 1, 2, … (2.3)
Если существует предел итерационной последовательности , то он является корнем уравнения f(x) = 0, то есть f(ξ) = 0.
Рис. 2.4. Сходящийся процесс итераций

Слайд 17

На рис. 2.4 показан процесс получения очередного приближения по методу итераций.

Последовательность приближений сходится к корню ξ.
Теоретические основы для применения метода итераций дает следующая теорема
Теорема 2.3. Пусть выполняются условия:
корень уравнения x = φ(x) принадлежит отрезку [a, b];
все значения функции φ(x) принадлежат отрезку [a, b], т.е. a ≤ φ(x) ≤ b;
существует такое положительное число q < 1, что производная во всех точках отрезка [a, b] удовлетворяет неравенству .
Тогда:
Итерационная последовательность сходится при любом .
Предел итерационной последовательности является корнем уравнения т.е. если , то ξ = φ(ξ).
Справедливо неравенство, характеризующее скорость сходимости итерационной последовательности
(2.4)
Как мы видим, эта теорема ставит довольно жесткие условия, которые необходимо проверить перед применением метода итераций. Если производная функции φ(x) по модулю больше единицы, то процесс итераций расходится (рис.2.5).

Слайд 18

рис. 2.5
Пример 2.5. Найти корень уравнения sin 5x + x2 –

1 = 0 с точностью до ε = 0,001 методом итераций, используя найденный выше отрезок [0,2; 0,3].
Решение. Преобразуем уравнение f(x) = 0 к виду x = φ(x):
sin 5x + x2 – 1 = 0, => sin 5x = 1 – x2, => x = arcsin(1 – x2)/5
φ(x) = arcsin(1 – x2)/5
Проверим условия теоремы. Так как функция φ(x) монотонна на отрезке [0,2; 0,3], то нетрудно показать, что верны следующие неравенства:

Слайд 19

Все условия теоремы выполнены, мы можем применить метод итераций. Выполним вычисления

в программе Excel:
1) Вводим в ячейки обозначения и формулы, как показано в табл.2.6;
2) Ячейку B3 с помощью маркера заполнения копируем вниз до ячейки B6; аналогично копируем ячейку C2 до ячейки C6, ячейку D2 — до ячейки D6; Выделим диапазон A2:A3 и с помощью маркера заполнения копируем вниз, до ячейки A6.
Таблица 2.6
Результаты расчетов приведены в табл. 2.7.
Таблица 2.7
В качестве условия сходимости итерационных методов часто используется неравенство
|xk – x k–1| ≤ ε. (2.5)

Слайд 20

По этому критерию приближенным значением с точностью ε = 0,001 является

x5 = 0,244707 ≈ 0,245.
Создадим в программе Excel функции для решения уравнения методом итераций.
Приведем текст программы-функции iter для решения уравнения методом итераций в программе Mathcad и результат вычисления корня:
Параметры программы iter(φ, x0, ε):
φ — имя функции в правой части уравнения x = φ(x);
x0 — начальное приближение;
ε — точность приближения соответствующая формуле (2.4.3).
Результат расчета 0,244627588 ≈ 0,2446 с начальным значением 0,2 и точностью 0,0001 содержит больше верных знаков, чем корень, полученный в программе Excel с меньшей точностью 0,001.

Слайд 21

2.2.3 Метод хорд
Метод хорд заключается в замене кривой y = f(x)

отрезком прямой, проходящей через точки (a, f(a)) и (b, f(b)) (см. рис. 2.6). Абсцисса точки пересечения прямой с осью OX принимается за очередное приближение.
Чтобы получить расчетную формулу метода хорд, запишем уравнение прямой, проходящей через точки (a, f(a)) и (b, f(b)) и, приравнивая y нулю, найдем x:
рис.2.6. Метод хорд

Слайд 22

Алгоритм метода хорд:
1) Пусть k = 0;
2) Вычислим следующий номер итерации:

k = k + 1; Найдем очередное k-ое приближение по формуле: xk = a – f(a)(b – a)/(f(b) – f(a)); Вычислим f(xk).
3) Если f(xk) = 0 (корень найден), то переходим к 5).
Если f(xk)f(b)>0, b = xk, иначе — переменной a = xk.
4) Если |xk – xk-1| > ε, то переходим к шагу 2);
5) Выводим значение корня xk.
6) Конец.
Замечание. Действия третьего пункта аналогичны действиям метода половинного деления. Однако в методе хорд на каждом шаге может сдвигаться один и тот же конец отрезка (правый или левый), если график функции в окрестности корня выпуклый вверх (рис. 2.6, a)) или вогнутый вниз (рис. 2.6, b)). Поэтому в критерии сходимости используется разность соседних приближений.
Пример 2.6. Применим метод хорд к уравнению sin 5x + x2 – 1 = 0 и отрезку [0,2; 0,3] для определения корня с точностью до ε = 0,001.

Слайд 23

Решение. Проведем расчеты в программе Excel:
1) В ячейки A1:H1 запишем заголовки

столбцов как в табл. 2.6;
2) В ячейку B3 запишем формулу =ЕСЛИ(C2*E2<0;B2;D2) и затем ячейку B3 протянем маркером заполнения до ячейки B10;
3) В ячейку C2 запишем формулу =SIN(5*B2)+B2^2-1 и затем ячейку C2 протянем маркером заполнения до ячейки C10;
4) В ячейку D2 запишем формулу =B2-C2*(F2-B2)/(G2-C2) и затем ячейку D2 протянем маркером заполнения до ячейки D10;
5) В ячейку E2 запишем формулу =SIN(5*D2)+D2^2-1 и затем ячейку E2 протянем маркером заполнения до ячейки E10;
6) В ячейку F3 запишем формулу =ЕСЛИ(C2*E2<0;D2;F2) и затем ячейку F3 протянем маркером заполнения до ячейки F10;
7) В ячейку G2 запишем формулу =SIN(5*F2)+F2^2-1 и затем ячейку G2 протянем маркером заполнения до ячейки G10;
8) В ячейку H2 запишем формулу =ABS(F2-B2) и затем ячейку H2 протянем маркером заполнения до ячейки H10;
В таблице 2.8 приведены результаты. Необходимая точность достигается на шаге k = 4.
Таблица 2.8

Слайд 24

Решение в программе Mathcad:
Как видим, результаты расчетов согласуются с предыдущими ответами.

Слайд 25

2.2.4 Метод Ньютона (метод касательных)
Пусть найдено приближенное значение корня уравнения

f(x) = 0, обозначим его xn. Расчетная формула метода Ньютона для определения очередного приближения xn+1 может быть получена двумя способами.
Первый способ выражает геометрический смысл метода Ньютона и состоит в том, что вместо точки пересечения графика функции y = f(x) с осью OX, мы ищем точку пересечения с осью OX касательной, проведенной к графику функции в точке (xn, f(xn)) как показано на рис. 2.6. Уравнение касательной имеет вид .
Рис. 2.7. Метод Ньютона (касательных)

Слайд 26

В точке пересечения касательной с осью OX переменная y = 0.

Приравнивая y нулю, выразим x и получим формулу метода касательных:
(2.6)
Второй способ. Разложим функцию f(x) в ряд Тейлора в окрестности точки x = xn:
Ограничимся линейными относительно (x – xn) слагаемыми, приравняем нулю f(x) и, выразив из полученного уравнения неизвестное x и обозначив его через xn+1, мы получим формулу (2.6).
Приведем достаточные условия сходимости метода Ньютона.

Слайд 27

Теорема 2.4. Пусть на отрезке выполняются условия:
1) функция и ее производные

и непрерывны;
2) производные и отличны от нуля и сохраняют определенные постоянные знаки;
3) (функция меняет знак на отрезке).
Тогда существует отрезок , содержащий искомый корень уравнения на котором итерационная последовательность сходится. Если в качестве нулевого приближения выбрать ту граничную точку , в которой знак функции совпадает со знаком второй производной, т.е. , то итерационная последовательность сходится монотонно (рис.2.8).
Доказательство. Так как непрерывна, меняет знак и монотонна на , то
— интервал изоляции корня. Обозначим искомый корень через .
Рассмотрим функцию и найдем ее производную
Итак, непрерывна на , обращается в нуль в точке , так как в этой точке обращается в нуль . Следовательно, существует такой отрезок
( ),что . Если возьмем ту часть отрезка, где , то

Слайд 28

, следовательно, функция возрастающая, но тогда последовательность является монотонной.
Рис.

2.8. Достаточные условия сходимости метода Ньютона
Замечание. Отметим, что метод хорд как раз идет с противоположной стороны, и оба этих метода т.о. могут друг друга дополнять, а возможен и комбинированный метод хорд-касательных.

Слайд 29

Пример 2.7. Уточнить до 0,000001 методом Ньютона корень уравнения sin 5x

+ x2 – 1 = 0. За начальное значение принять x0 = – 0,7.
Решение. Найдем производную .
В программе Excel введем расчетные формулы:
1) Введем формулы и обозначения в ячейках диапазона A1:D3 и скопируем вниз маркером заполнения ячейки с формулами: B3 — до B5, C2 — до C5, D2 — до D5;
Таблица 2.9
Результаты расчетов приведены в таблице 2.10. Получено значение корня – 0,726631609 ≈ – 0,726632 с погрешностью 0,000001.
Таблица 2.10

Слайд 30

2.2.5 Метод секущих
Метод секущих может быть получен из метода Ньютона при

замене производной приближенным выражением — разностной формулой:
(2.7)
В формуле (2.7) используются два предыдущих приближения xn и xn–1. Поэтому при заданном начальном значении x0 необходимо вычислить следующее приближение x1 каким-нибудь методом, например, методом Ньютона с приближенной заменой производной по формуле
Алгоритм метода секущих:
1) Заданы начальное значение x0 и погрешность ε. Вычислим x1 = x0 – f(x0)ε/( f(x0 + ε) – f(x0));
2) Для n = 1, 2, … пока выполняется условие |xn – xn–1| > ε вычисляем xn+1 по формуле (2.7):

Слайд 31

2.3 Системы нелинейных уравнений
Система n нелинейных уравнений с n неизвестными

имеет вид
fk(x1, x2, …, xn) = 0, 1 ≤ k ≤ n. (2.8)
Систему двух нелинейных уравнений
(2.9)
можно решить приближенно графическим способом. Для этого достаточно преобразовать систему к виду
(2.10)
построить графики функций y = y1(x), y = y2(x) и найти координаты точек пересечения графиков (рис. 2.9). При использовании электронных таблиц или математических пакетов решение можно уточнить графически, сужая отрезок [a, b] около корня xs.

Слайд 32

Рис. 2.9. Графическое решение системы двух уравнений
Пример 2.8. Решить графически систему

двух уравнений
Решение. Сначала приведем алгоритм определения решения системы двух уравнений графическим методом:
1) Преобразуем систему к виду .
2) Построим графики функций и y = cosx, подбирая отрезок [a, b] изменения переменной x так, чтобы графики пересекались.

Слайд 33

3) Изменяя, т.е. уменьшая отрезок [a, b], уточняем решение (xs, ys).

Решение в программе Excel. Так как область определения функции задается условием x ≥ 0, выберем для построения графиков отрезок [0; 1] с шагом изменения 0,1. Если графики не будут пересекаться, то вместо отрезка [0; 1] возьмем отрезок [1; 2] и т.д.
1) В ячейки A2, A3 запишем соответственно 0 и 0,1; выделим диапазон A2:A3 и маркером заполнения протянем вниз до ячейки A12.
2) В ячейку B2 запишем формулу = корень(A2); выделим B2 и маркером заполнения протянем вниз до ячейки B12.
3) В ячейку C2 запишем формулу =cos (A2); выделим C2 и маркером заполнения протянем вниз до ячейки C12.
4) Выделим диапазон A2:C12 и построим диаграмму «Точечная». Графики, как видим, пересекаются. Проведем настройку диаграммы.
Щелкнем правой кнопкой мыши по диаграмме и выберем «параметры диаграммы», вкладку «Легенда» и снимем флажок с параметра «показать легенду».
Щелкнем правой кнопкой мыши по диаграмме и выберем «параметры диаграммы», вкладку «линии сетки», отметим «промежуточные линии» оси X и «промежуточные линии» оси Y.

Слайд 34

Щелкнем правой кнопкой мыши по оси X диаграммы и выберем «формат

оси», в появившемся окне выберем вкладку «Шкала» и введем «минимальное значение — 0», «максимальное значение — 1», «цена основных делений — 0,1», «цена промежуточных делений — 0,1».
Аналогично, для оси Y диаграммы выберем «цена основных делений — 0,1», «цена промежуточных делений — 0,1». Полученная диаграмма приведена на рис.2.10.
Рис.2.10. Графическое решение системы двух уравнений

Слайд 35

5) На диаграмме увидим, что графики пересекаются между значениями x =

0,6 и x = 0,7. Заменим отрезок [0; 1] на отрезок [0,6; 0,7], для чего введем в ячейки A2, A3 числа 0,60 и 0,61, выделим диапазон A2:A3 и маркером заполнения протянем вниз до ячейки A12. Графики изменятся. Щелкнем правой кнопкой мыши по оси X диаграммы и выберем «формат оси», в появившемся окне выберем вкладку «Шкала» и введем «минимальное значение — 0,6», «максимальное значение — 0,7», «цена основных делений — 0,01», «цена промежуточных делений — 0,01». И для оси Y диаграммы внесем изменения: «цена основных делений — 0,01», «цена промежуточных делений — 0,01». Теперь мы увидим, что графики пересекаются между значениями x = 0,64 и x = 0,65.
6) Аналогичными действиями заменим отрезок [0,6; 0,7] на новый отрезок [0,64; 0,65] с шагом изменения 0,01. Получим x ≈ 0,641; y ≈ 0,801.
Процесс уточнения можно продолжать и дальше. Погрешность полученного решения составляет приблизительно 0,001 для обеих неизвестных.
Для уточнения решения (xs, ys) можно также применить метод итераций или метод Ньютона, которые рассматриваются ниже.

Слайд 36

2.3.1 Метод итераций
Приведем систему (2.8) к виду, удобному для итераций
xk =

φk(x1, x2, …, xn), 1 ≤ k ≤ n. (2.11)
Выберем начальное приближение к корню (x10, x20, …, xn0) и последующие приближения вычислим по формулам
xks +1 = φk(x1s, x2s, …, xns), 1 ≤ k ≤ n, s = 0, 1, 2, … (2.12)
Приведем без доказательства достаточные условия сходимости метода итераций. Обозначим точное решение системы (2.8) . Назовем ε-окрестностью точки множество точек x = (x1, x2, …, xn), удовлетворяющих условиям
Теорема 2.5. Пусть в некоторой ε-окрестности точного решения частные производные существуют и удовлетворяют одному из трех неравенств
(2.13)
где . Если начальное приближение принадлежит ε-окрестности точного решения, то метод простой итерации (2.12) сходится к точному решению.

Слайд 37

2.3.2 Метод Ньютона
Строгие формулировки теорем об условиях сходимости метода Ньютона достаточно

громоздки, на практике часто ограничиваются следующим рассуждением.
Пусть для системы нелинейных уравнений
fk(x1, x2, …, xn) = 0, 1 ≤ k ≤ n,
в некоторой ε-окрестности точного решения не равен нулю определитель
матрицы частных производных (матрицы Якоби):
Тогда существует начальное приближение , принадлежащее ε-окрестности точного решения (достаточно близкое к точному решению), что метод Ньютона сходится к точному решению.

Слайд 38

(2.14)
Пример 2.10. Решить методом Ньютона систему нелинейных уравнений из примера 2.7
Решение.

Найдем определитель матрицы Якоби:
Очевидно, что в некоторой окрестности точки (0,641; 0,801) определитель матрицы Якоби не равен нулю. Найдем матрицу, обратную к матрице Якоби:

Слайд 39

Теперь для данной системы метод Ньютона можно записать в виде итерационных

формул:
В следующей таблице приведены результаты расчетов по этим формулам с начальным приближением (0,5; 0,5):
Третий шаг итераций дает результаты, совпадающие до трех цифр с решением примера 2.7, а пятый и шестой шаги дают значения, совпадающие друг с другом точно. Это говорит о том, что достигнута максимальная точность. Эти результаты объясняются высокой скоростью сходимости метода Ньютона.

Слайд 40

Контрольные вопросы.
1. Как определяется понятие «корень уравнения»?
2. Какие уравнения называются алгебраическими?
3.

Какие уравнения называются трансцендентными?
4. В чем заключается процедура отделения корней уравнения?
5. В чем суть аналитического метода отделения корней?
6. В чем суть графического метода отделения корней?
7. На какой теореме математического анализа основан алгоритм отделения корней уравнения?
8. Перечислите методы уточнения корней?
9. В чем суть метода половинного деления?
10. В чем суть метода простых итераций?
11. В чем суть метода Ньютона?
12. В чем суть метода секущих?
13. В чем суть метода хорд?
14. Как оценивается погрешность метода половинного деления?
15. Как оценивается погрешность метода итераций?

Алгебраические и трансцендентные уравнения

Содержание

2.1 Методы локализации корней 2.1.1 Аналитический метод Теоретической основой алгоритма отделения корней

Пример 2.1. Отделить корни уравнения sin 5x + x2 – 1

Теорема 2.2. Если непрерывная функция f(x) монотонна на отрезке [a, b]

2.1.2 Графический метод Для отделения корней уравнения зачастую бывает удобно применять графический

рис 2.2Из рис. 2.2 видно, что графики пересекаются в четырех точках,

1) Если функция f(x) на концах отрезка [a, b] принимает значения

рис 2.3Рис. 2.3 показывает, что корень находится в интервале (0,44; 0,45),

Повторив всю процедуру еще раз, заменим в диапазоне A2:A12 интервал (0,44;

Табл. 2.2Изложенный метод можно охарактеризовать как метод деления отрезка на 10

2.2 Методы уточнения корней После того как найден интервал, содержащий корень,

Метод половинного деления (другие названия: метод бисекций, метод дихотомии) для решения

Пример 2.4. Уточнить методом бисекций с точностью до 0,01 корень уравнения

Таблица 2.4Приведенный алгоритм учитывает возможный случай «попадания в корень», т.е. равенство

Создадим в программе Excel пользовательские функции f(x) и bisect(a, b, eps)

2.2.2 Метод итераций Метод простых итераций для уравнения f(x) = 0

На рис. 2.4 показан процесс получения очередного приближения по методу итераций.

рис. 2.5Пример 2.5. Найти корень уравнения sin 5x + x2 –

Все условия теоремы выполнены, мы можем применить метод итераций. Выполним вычисления

По этому критерию приближенным значением с точностью ε = 0,001 является

2.2.3 Метод хорд Метод хорд заключается в замене кривой y = f(x)

Алгоритм метода хорд:1) Пусть k = 0;2) Вычислим следующий номер итерации:

Решение. Проведем расчеты в программе Excel:1) В ячейки A1:H1 запишем заголовки

Решение в программе Mathcad:Как видим, результаты расчетов согласуются с предыдущими ответами.

2.2.4 Метод Ньютона (метод касательных) Пусть найдено приближенное значение корня уравнения