Учебное пособие: Численные методы для решения нелинейных уравнений

Методы простой итерации, Зейделя, модифицированный метод Ньютона, метод наискорейшего спуска (см. [1] , [2] , [3] , [4] ) являются методами первого порядка – это значит, что имеет место неравенство , k =1, 2, . . . , где – константа, своя у каждого метода, зависящая от выбора начального приближения , функции fi , i = 1, 2, . . . , n , и их частных производных первого и второго порядков – точнее их оценок в некоторой окрестности искомого решения, которой принадлежит начальное приближение.

Метод Ньютона является методом второго порядка, то есть для него имеет место неравенство , k =1, 2, . . . , где – константа, зависящая от тех же величин, что и константа .

А теперь рассмотрим достаточные условия сходимости метода простой итерации и метода Ньютона.

Сходимость процесса простой итерации зависит от двух условий. Первое условие состоит в том, что какая-нибудь точка должна оказаться близкой к исходному решению . Степень необходимой близости зависит от функций j 1 , j 2 , . . . , j n . Это требование не относится к системам линейных уравнений, для которых сходимость процесса простой итерации зависит только от второго условия.

Второе условие связано с матрицей, составленной из частных производных первого порядка функций j 1 , j 2 , . . . , j n – матрицей Якоби

,

вычисленных в точке .

В случае, когда рассматривается система линейных алгебраических уравнений, матрица M состоит из постоянных чисел – коэффициентов, стоящих при неизвестных в правой части уравнения (3) . В случае нелинейных уравнений элементы матрицы M зависят, вообще говоря, от . Для сходимости процесса простой итерации достаточно, чтобы выполнялось неравенство: для из некоторой окрестности точного решения , которой должно принадлежать начальное приближение .

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

Предположим, что имеется начальное приближение к искомому решению системы (2) , функции непрерывны и имеют непрерывные частные производные до второго порядка в шаре , тогда, если выполнены условия:

1) Матрица Якоби системы (2) на начальном приближении имеет обратную и известна оценка нормы обратной матрицы ,

2) Для всех точек шара выполнено неравенство

при i , j = 1, 2, . . . , n ,

3) Выполнено неравенство

,

где L – постоянная 0 £ L £ 1,

4) Числа b , N , r подчинены условию a = nbNr < 0,4, тогда система уравнений (2) в шаре имеет единственное решение, к которому сходятся последовательные приближения (8) или (7’) , (9’) .

Для других методов условия сходимости имеют сложный вид, и мы отсылаем читателя к специальной литературе [1] , [2] , [3] , [4] .


6. Примерный перечень возможных исследований

1) Сравнение различных методов на экономичность при решении конкретной задачи:

· по числу операций на одной итерации;

· по числу итераций, необходимых для достижения заданной точности;

2) Зависимость числа итераций для достижения заданной точности:

· от выбора вида нормы;

· от выбора критерия окончания итерационного процесса по или по невязке ;

· от выбора начального приближения;

· от погрешности задания коэффициентов в уравнении.


7. Контрольные вопросы

1) Понятие о нелинейных системах уравнений в Rn .

2) Понятие приближенного и точного решения нелинейной системы уравнений.

3) Сущность графического метода отделения решения для системы двух нелинейных уравнений, каковы его преимущества и недостатки?

К-во Просмотров: 355
Бесплатно скачать Учебное пособие: Численные методы для решения нелинейных уравнений