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

⇐ ПредыдущаяСтр 5 из 5

Основная идея метода Ньютона состоит в выделении из уравнений системы линейных частей, которые являются главными при малых приращениях аргументов. Это позволяет свести исходную задачу к решению последовательности систем линейных уравнений.

Рассмотрим систему двух нелинейных уравнений с двумя неизвестными вида:

(3.10)

Пусть известно некоторое приближение , корня , . Тогда поправки , можно найти, решая систему:

(3.11)

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

(3.12)

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

Таким образом, решение системы уравнений по методу Ньютона состоит в построении итерационной последовательности:

(3.13)

где , - решения систем линейных уравнений, вида (3.12) на каждом шаге итерации.

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

Пример 3.2. Найти решение системы (3.8) методом Ньютона с точностью .

Решение. Начальные приближения , . Определим частные производные:

;

и, используя (3.12), построим систему линейных уравнений относительно поправок

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

определяем поправки на первом шаге итерации

Далее начальное приближение уточняем по формулам (3.13)

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

определяем поправки на втором шаге итерации

Далее иуточняем по формулам (3.13)

Определяем погрешностьпо формуле :

Таким образом, имеем решение: , .

Программа, реализующая метод Ньютона для указанной задачи, представлена на рис. 3.2. Исходные данные – начальные приближения , , точность и максимальное число итераций (табл. 3.2).

Таблица 3.2.

Исходные данные к программе решения системы

нелинейных уравнений

методом Ньютона

	A	B
	x0	-1
	y0	-0, 7
	e	0, 001
	n
	x	-1, 11149
	y	-0, 72253

Sub program6() x = Cells(1, 2) y = Cells(2, 2) e = Cells(3, 2) n = Cells(4, 2) For k = 1 To n F = 2 * Sin(x + 1) – y - 0.5 G = 10 * Cos(y - 1) – x + 0.4 Fx = 2 * Cos(x + 1) Fy = -1 Gx = -1 Gy = -10 * Sin(y - 1) D = Fx * Gy – Gx * Fy Dx = (G * Fy – F * Gy) / D Dy = (F * Gx – G * Fx) / D xk = x + Dx yk = y + Dy If Abs(xk - x) < e And Abs(yk - y) < e Then Cells(5, 2) = xk Cells(6, 2) = yk End End If x = xk y = yk Next k MsgBox " решение не найдено" End End Sub

Рис. 3.2. Программа, реализующая метод Ньютона на языке VBA.

Пример 3.3. Найти решение системы (3.8) с помощью программы Excel.

Порядок решения.

1) Подключить надстройку «Поиск решения» через Кнопка «Офис»-Параметры Excel-Надстройки-Надстройки Excel-Перейти (рис. 3.3);

2) Ввести в ячейки A1, B1, C1, D1 заголовки столбцов (рис. 3.4а);

3) В ячейку A2 – начальное приближение для :

4) В ячейку B2 – начальное приближение для :

5) В ячейку C2 – формулу =2*SIN(A2+1)-B2-0, 5

6) В ячейку D2 – формулу =10*COS(B2-1)-A2+0, 4

7) Вызвать диалоговое окно «Поиск решения»: Данные-Поиск решения (рис. 3.5)

8) В качестве целевой ячейки указываем результат вычисления левой части одного из уравнений, например, , т.е. ячейку C2

9) Для решения уравнения значение , поэтому выбираем переключатель «значение», а в соответствующее поле вводим 0

10) Установив курсор в поле «Изменяя ячейки», выделяем ячейки незвестных , , т.е. A2: B2

11) Остальные уравнения системы рассматриваются как дополнительные ограничения ( ). Нажимаем кнопку «Добавить», отмечаем мышью ячейку D2и вводим=0

12) Нажимаем кнопку «Выполнить». Если решение найдено, появляется окно сообщения (рис. 3.6). Нажимаем кнопку ОК.

13) В ячейках A2: B2 - решение системы (рис. 3.4б), т.е ,

		а) б)
Рис. 3.3. Подключения надстройки «Поиск решения».		Рис. 3.4. Рабочий лист до и после выполнения поиска решения.


Рис. 3.5. Параметры окна «Поиск решения».

Рис. 3.6. Сообщение о завершении поиска решения.

ЛИТЕРАТУРА

1. Калиткин Н.П. Численные методы. М.: Наука, 1978. - 512 с.

2. Турчак Л.И., Плотников П.В. Основы численных методов: Учебное пособие. М.: ФИЗМАТЛИТ, 2003. – 304 с.

3. Васильев А.Н. Научные вычисления в Microsoft Excel. М.: Издательский дом " Вильямс", 2004. – 512 с.

4. Ларсен У.Р. Инженерные расчеты в Excel. М.: Издательский дом " Вильямс", 2004. – 544 с.

5. Попов В.И. Численные методы расчета мостовых конструкций на ЭВМ. М.: 1981. – 78 с.

6. Методические указания к лабораторным и самостоятельным работам по курсам «Информатика» и «Вычислительная математика». Численные методы. Часть 1. / Сост.: Ф.Г.Ахмадиев, Ф.Г.Габбасов, Р.Ф.Гиззятов, И.В.Маланичев. Казань: КГАСУ, 2011. – 32 с.

7. Методические указания к лабораторным и самостоятельным работам по курсам «Информатика» и «Вычислительная математика». Численные методы. Часть 2. / Сост.: Ф.Г.Ахмадиев, Ф.Г.Габбасов, Р.Ф.Гиззятов, И.В.Маланичев. Казань: КГАСУ, 2011. – 36 с.

⇐ Предыдущая 1 2 3 45