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

Содержание

Слайд 2

Введем обозначения

Введем обозначения

Слайд 3

Тогда исходную систему запишем относительно векторной функции F векторного аргумента Следовательно,

Тогда исходную систему запишем

относительно векторной функции F
векторного аргумента

Следовательно,

исходную задачу можно

рассматривать как
задачу о нулях нелинейного отображения

и

В такой постановке задача является
обобщением задачи о нахождении решения
нелинейного уравнения для случая задачи
большой размерности.

Слайд 4

Однако, переход от к вносит в задачу нахождения нулей свою специфику.

Однако, переход от

к

вносит в задачу нахождения нулей свою
специфику.
Метод Ньютона решения

систем нелинейных уравнений.
Пусть исходная система приведена к виду:

……………………..

Слайд 5

или в компактной форме: Для задачи о неподвижной точке нелинейного отображения запишем

или в компактной форме:
Для задачи о неподвижной точке нелинейного
отображения
запишем

Слайд 6

формальное равенство: где определяет метод простых итераций и Пусть известно -

формальное равенство:

где

определяет метод простых

итераций и

Пусть известно

- е приближение

одного из

изолированных корней

векторного уравнения

Слайд 7

Тогда точный корень уравнения можно представить: где поправка (погрешность корня). Имеем

Тогда точный корень уравнения можно
представить:

где

поправка (погрешность корня). Имеем

Предполагая, что

непрерывно дифференцируема

в некоторой
выпуклой области, содержащей

и

разложим левую часть уравнения

по степеням малого

вектора

ограничиваясь линейными

Слайд 8

членами Под производной следует понимать матрицу Якоби системы функций относительно переменных т.е.

членами

Под производной

следует понимать

матрицу Якоби системы функций

относительно переменных

т.е.

Слайд 9

Тогда Если то Следовательно, метод Ньютона для решения исходной системы состоит

Тогда

Если

то

Следовательно, метод Ньютона для решения
исходной системы состоит в построении
итерационной

последовательности:

k=0,1,2,

Если поправки достаточно малы, счет закончен

Слайд 10

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

Пример. Найти методом Ньютона решение
системы уравнений

исходя из начального приближения

Слайд 11

Полагая имеем

Полагая

имеем

Слайд 12

Подставляя данные, получаем Составим матрицу Якоби:

Подставляя данные, получаем

Составим матрицу Якоби:

Слайд 13

При этом Т.к. то найдем обратную ей матрицу

При этом

Т.к.

то найдем обратную ей матрицу

Слайд 14

Получим первое приближение:

Получим первое приближение:

Слайд 15

а Аналогично находим дальнейшие приближения В результате получим ( )

а

Аналогично находим дальнейшие приближения
В результате получим ( )

Слайд 16

Решение нелинейных систем методами спуска Общий недостаток всех рассмотренных ранее методов

Решение нелинейных систем методами спуска
Общий недостаток всех рассмотренных ранее
методов решения систем

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

Пусть задана система Образуем новую функцию Т.к. эта функция не отрицательна,

Пусть задана система

Образуем новую функцию

Т.к. эта функция не отрицательна, то всегда
найдется

точка

такая, что

Слайд 18

Слайд 19

Т.е., если удается получить точку минимизирующую функцию и если при этом

Т.е., если удается получить точку

минимизирующую функцию

и если

при этом окажется, что

то

точка

истинное решение системы, т.к.

Слайд 20

Последовательность точек приближений к точке минимума функции получают по формуле где

Последовательность точек

приближений к точке

минимума

функции

получают по формуле

где

вектор, определяющий

направление минимизации, а

скалярная

величина, характеризующая

величину шага
минимизации.
Слайд 21

Исходя из геометрического смысла задачи, итерационный метод называется методом спуска. При

Исходя из геометрического смысла задачи,
итерационный метод называется методом
спуска. При выборе направления

спуска
определяющим является направление, в
котором функция убывает наиболее быстро.
Слайд 22

Направление наибольшего возрастания функции в данной точке показывает её градиент в этой точке. Поэтому антиградиент функции

Направление наибольшего возрастания
функции в данной точке показывает её
градиент в этой точке.

Поэтому

антиградиент функции

Слайд 23

Тогда суть градиентного метода Достоинство градиентного метода решения нелинейных систем –

Тогда суть градиентного метода

Достоинство градиентного метода решения
нелинейных систем – глобальная сходимость.
Главный

недостаток – медленная сходимость.
Слайд 24

Пример. Найти максимум функции Методом скорейшего спуска при ограничениях: Функция является

Пример. Найти максимум функции

Методом скорейшего спуска при ограничениях:

Функция

является выпуклой, поэтому

её

локальный максимум совпадает с глобаль-
ным. Оптимум достигается внутри области
решений системы ограничений. Градиент
Слайд 25

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

Пусть исходная точка

Подставляя

координаты

в выражение градиента,

получим

Используя формулу

получим

Используя скалярное произведение векторов

Слайд 26

найдём Отсюда а=0,7449 – величина длины шага. Тогда координаты следующей точки и градиента:

найдём

Отсюда а=0,7449 – величина длины шага.
Тогда координаты следующей точки и
градиента:

Слайд 27

Выполняя в цикле представленные расчёты, процесс итерации заканчиваем при достижении заданной

Выполняя в цикле представленные расчёты,
процесс итерации заканчиваем при
достижении заданной точности отклонения

В

нашем примере итоги расчёта:
Слайд 28

Пример. Найти направление наискорейшего возрастания функции в точке и вычислить значение

Пример. Найти направление наискорейшего
возрастания функции

в точке

и

вычислить значение производной в этом
направлении.
Решение. Координаты

градиента данной
функции:

Итак,

В точке

градиент имеет координаты