Численные методы

Содержание

Слайд 2

Цель лекции Изучить два метода вычисления корней нелинейных уравнений, а именно:

Цель лекции

Изучить два метода вычисления корней нелинейных уравнений, а именно:
Метод хорд

(метод секущих)
Метод касательных(метод Ньютона)
Слайд 3

Прежде чем приступить к изучению новых методов решения нелинейных уравнений, вспомним

Прежде чем приступить к изучению новых методов решения нелинейных уравнений, вспомним

(смотри лекцию «Способы отбора корней нелинейных уравнений») как графическим способом отобрать корни данного уравнения:
f(х) = 0.
Слайд 4

Для того, чтобы найти графически интервалы изоляции действительных корней данного скалярного

Для того, чтобы найти графически интервалы изоляции действительных корней данного скалярного

уравнения
f(х) = 0 необходимо:
1). Представить уравнение f(х) = 0 в виде
f1 (х) = f2 (х).
2). Построить графики функций
у = f1(х) и у = f2 (х),
3). Определить приближенно по графику абсциссу точки пересечения этих графиков х0.
4). Определить промежуток изоляции [а; b], содержащий корень х0.
Слайд 5

Метод хорд (метод секущих) Пусть требуется вычислить действительный корень уравнения f

Метод хорд (метод секущих)
Пусть требуется вычислить действительный корень уравнения
f

(х) = 0,
изолированный на отрезке [а;Ь]. Рассмотрим график функции у = f(х). Пусть f (а) < 0 и f (b) > 0.
Слайд 6

График функции у = f(х).

График функции у = f(х).

Слайд 7

Точки графика A(a; f(a)) и B(b; f(b)) соединим хордой. За приближенное

Точки графика A(a; f(a)) и B(b; f(b)) соединим хордой.
За

приближенное значение искомого корня примем абсциссу х1 точки пересечения хорды АВ с осью Ох.
(смотри следующий рисунок)
Слайд 8

Слайд 9

Это приближенное значение находится по формуле (b – a) f (a)

Это приближенное значение находится по формуле
(b – a) f (a)

x1 = a - --------------------- ,
f (b) –f (a)
где х1е[а; b].
Слайд 10

Пусть, например, f(х1) (b – х1) f (х1) x2 = х1

Пусть, например, f(х1)<0, тогда за новый (более узкий) промежуток изоляции корня

можно принять [x 1; b]. Соединив точки А(х1; f(x1)) и В(b; f(b)), получим в точке пересечения хорды с осью Ох второе приближение х2, которое вычислим по формуле:
(b – х1) f (х1)
x2 = х1 - --------------------- и т. д.
f (b) –f (х1)
Слайд 11

Последовательность чисел а, х1, х2,... стремится к искомому корню уравнения. Если

Последовательность чисел а, х1, х2,... стремится к искомому корню уравнения.
Если

было бы f(х1) > 0, то за новый промежуток изоляции корня можно было бы принять [а; х1] и тогда второе приближение х2 вычисляли бы по формуле
(х1 – a) f (a)
x2 = a - --------------------- ,
f (х1) –f (a)
Слайд 12

Вычисление приближенных значений корней уравнения следует вести до тех пор ,

Вычисление приближенных значений корней уравнения следует вести до тех пор ,

пока не перестанут изменятся те десятичные знаки, которые мы хотим сохранить в ответе ( т. е. пока не будет достигнута заданная степень точности).
Слайд 13

Алгоритм выполнения задачи Методом хорд решить уравнение f(х)=0 с точностью до

Алгоритм выполнения задачи

Методом хорд решить уравнение f(х)=0 с точностью до

ε.
Решение:
Вычислим приближенное значение корня с заданной точность ε. (уточним корень, найденный графически) .
Слайд 14

1). Найдем f(а), для этого в f(х) вместо х подставим а.

1). Найдем f(а), для этого в f(х) вместо х подставим а.

Определим знак f(а).
2). Найдем f(b), для этого в f(х) вместо х подставим b. Определим знак f(b).
3). Найдем первое приближенное значение корня по формуле :
(b – a) f (a)
x1 = a - --------------------- ,
f (b) –f (a)
Слайд 15

4). Найдем f(x1) для этого в f(x) вместо х подставим x1.

4). Найдем f(x1) для этого в f(x) вместо х подставим x1.
5).

Определим знак f(x1) ;
6). Найдем новый (более узкий) промежуток изоляции:
а)Если f(x1) имеет знак противоположный знаку f(а), то за новый промежуток примем [а;x1,].
б)Если f(x1) имеет знак противоположный f(b), то за новый промежуток примем |х1;b].
Слайд 16

7). Найдем второе приближение корня в случае а) по формуле: (b

7). Найдем второе приближение корня в случае а) по формуле:

(b – х1) f (х1)
x2 = х1 - --------------------- ,
f (b) –f (х1)
в случае б) по формуле:
(х1 – a) f (a)
x2 = a - --------------------- ,
f (х1) –f (a)
Слайд 17

8). Найдем f(x2) , для этого в f(х) вместо x подставим

8). Найдем f(x2) , для этого в f(х) вместо x подставим

x2.
9). Определим f(x2) . Сравним его со знаками на концах промежутка изоляции (найденного в п.6).).
Если знак f(x2) противоположен знаку f(x1), то за новый промежуток изоляции примем отрезок [х1;х2]
Если знак f(x2) противоположен знаку f(b), то за новый промежуток изоляции примем [х2; b].
Слайд 18

10). Вычисление приближенных корней уравнения ведем до тех пор, пока не

10). Вычисление приближенных корней уравнения ведем до тех пор, пока не

перестанут изменяться те десятичные знаки, которые мы хотим сохранить в ответе.
11) Результаты вычислений занесем в таблицу:
Слайд 19

Метод касательных (метод Ньютона) Пусть корень уравнения f(х) = 0 изолирован

Метод касательных (метод Ньютона)

Пусть корень уравнения f(х) = 0 изолирован

на отрезке [а: b].
Пусть снова f(а) < 0 и f(b) > 0, причем первая производная на этом отрезке не меняет своего знака.
Тогда в отрезке [а; b ] имеется один корень уравнения f(х) = 0.
Возьмем на отрезке [а;b] такое число х0., при котором f´ (х0) имеет тот же знак, что и f ´´(х0), т.е. f´(х0) f´´(х0 ) > 0 ( в частности, за х0 может быть принят тот из концов отрезка [а;b], в котором соблюдено это условие).
Сохранение знака второй производной на отрезке означает, что кривая либо только выпукла, либо только вогнута на нем.
Слайд 20

Проведем в точке М0(x0; f(x0)) касательную к кривой у = f(х).

Проведем в точке М0(x0; f(x0)) касательную к кривой у = f(х).
За

приближенное значение корня примем абсциссу точки пересечения этой касательной с осью Ох.
Чтобы найти эту абсциссу напишем уравнение касательной в точке Мо. y – f(x0) = f' (x0) ( x – x0), при у=0 и х =х1 получим:
– f(x0) = f' (x0) ( x1 – x0), или x1 – x0 = - f(x0)/ f' (x0) ,
отсюда найдем х1 = х0 - f(x0)/ f' (x0).
Слайд 21

График заданной функции

График заданной функции

Слайд 22

Применив этот прием вторично в точке М1(x1;f(x1)), найдем f (x1) x2

Применив этот прием вторично в точке М1(x1;f(x1)), найдем
f (x1)
x2

= x1 - ----------
f´ (x1)
И т. д f (xn-1)
xn = xn-1 - ----------
f´ (xn-1)
Слайд 23

Полученная таким образом последовательность х0, х1, х2,…имеет своим пределом искомый корень.

Полученная таким образом последовательность х0, х1, х2,…имеет своим пределом искомый корень.

Если х - точный корень уравнения f(х) = 0, изолированный на отрезке [а;b ], а ξ –приближенное значение корня , найденное методом хорд, то оценка погрешности этого приближенного значения такова:
│f (ξ) │2 f´ (x)
│ x - ξ │< ------------------- max ----------- .
2 [ a:b] (f´´(x))^2
Приняв за a и b концы промежутка изоляции, на котором найдено приближенное значение корня.
Слайд 24

Алгоритм выполнения задачи Методом касательных (методом Ньютона) решить уравнение f(х) =

Алгоритм выполнения задачи

Методом касательных (методом Ньютона) решить уравнение f(х) = О

с точностью до ε
Решение:
Вычислим приближенное значение корня с заданной точностью ε. (уточним корень, найденный графически )
1). Найдем f'(х) и f"(х) для данной функции f (х).
2). Возьмем на отрезке изоляции [а;b] такое число х0 , при котором f´(х0) имеет тот же знак, что и вторая производна, т. е. f'(х) ∙ f´´(х)>0
( в частности за х 0 может быть принят тот из концов отрезка [а; b], в котором соблюдено это условие
Слайд 25

3) Найдем f'(х0 ). 4). Найдем первое приближенное значение корня х1

3) Найдем f'(х0 ).
4). Найдем первое приближенное значение корня х1 по

формуле:
x1 = x0 -. f (x0)/ f´ (x0)
5). Найдем значения f/(х1), подставив x1 в f /(х) вместо х0 .
6)Найдем значение f/(х2), затем по формуле найдем второе значение корня:
x2 = x1 - f (x1)/f´ (x1)
Слайд 26

7). Таким образом находим n-ое приближенное значение корня хn по формуле:

7). Таким образом находим n-ое приближенное значение корня хn по формуле:

f (xn-1)
xn = xn-1 - ----------
f´ (xn-1)
Слайд 27

8). Вычисление приближенных значений корней уравнения ведем до тех пор, пока

8). Вычисление приближенных значений корней уравнения ведем до тех пор,

пока не перестанут изменяться те десятичные знаки, которые мы хотим сохранить в ответе (т.е. пока не будет достигнута заданная степень точности │хn-хn-1│<ε
Слайд 28

9). Результаты вычислений занесем в таблицу:

9). Результаты вычислений занесем в таблицу: