А.В. Матисов Решение нелинейных уравнений

10

10.Искомый корень равен х и функция равна у.

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

2.3.4. Метод хорд

Метод хорд по эффективности и идее можно сравнить с предыдущим. Пусть имеется интервал изоляции корня [x1, x2]. Функция f(x) имеет на [x1, x2] непрерывные производные f '(х) и f ''(x). Точки М1(х1, у1) и М2(x2, у2) (рис. 2.6), где y1=f(x1) и y2=f(x2), соединим хордой М1М2, которая пересечет ось абсцисс в точке N(x3, 0). Уравнение этой хорды

( x-х1 ) / ( х2- х1 ) = ( y-y1 ) / ( у2-у1 ). (2.12)

При у=0 получим

x = x1 –y1 ( х2-х1 ) / ( у2-у1 ) (точка хЗ на рис. 2.6)

(2.13)

Значение х=хЗ является первым приближением к искомому корню. Найдем следующее приближение. Вычислим уЗ=f(хЗ). В зависимости от знака f(x3) новым интервалом изоляции будет либо [x1, x3], либо [x3, x2]. В нашем примере (рис. 2.6) корень находится в интервале [x3, x2]. Как видно из рис. 2.6, второе приближение корня определяется с помощью

хорды М3М2.

Алгоритм метода хорд:

1.Выберем точность вычислений корня Е.

2.Установим границы интервала изоляции х1 и х2.

3.Вычислим y1=f(x1) и y2=f(x2).

4.Вычислим приближение x3=x1-y1 (х2-х1)/(у2-у1).

5.Вычислим |уЗ|=f(х3).

6.Если y3<E, то идти к пункту 9.

7.Если у1 у3<0, то x2=x3 и y1=y3, иначе х1=х3 и y1=y3.

8.Идти к пункту 4.

9.Искомый корень равен х3 и функция равна у3.

11

Рис. 2.6. Графическая иллюстрация метода хорд

Как следует из приведенного алгоритма, вычисления по методу хорд прекращают, когда модуль функции f(x) при очередном приближении станет меньше Е.

2.3.5. Метод Монте-Карло (статистических испытаний)

Метод заключается в использовании случайных чисел для уточнения корня нелинейной функции. Равномерно-распределенные случайные числа обычно генерируются ЭВМ в отрезке [0, 1], причем любое значение из этого отрезка равновероятно. Перевод случайных чисел из интервала [0, 1] на интервал [а, b] производится по формуле

Xi = a + ( b – a ) Ui ,

(2.14)

где Xi – случайное число в интервале [а, b]; Ui – случайное число в интервале [0, 1].

Рассмотрим алгоритм метода. Предполагаем, что решена задача отделения корней и интервал [а, b] (рис. 2.7) является интервалом изоляции корня, а производная f '(x) в интервале [а, b] знак не меняет.

1.Выберем точность вычислений корня Е.

2.Установим границы изоляции корня a и b.

3.Генерируем x случайным образом на интервале [a, b].

4.Вычисляем значение функции F=f(x) в случайной точке x.

12

5.Если |F|≥Е, то идти к пункту 3.

6.Искомый корень равен х и функция равна F.

Рис. 2.7. Иллюстрация к методу Монте-Карло

Поскольку x в интервале [а, b] каждый раз выбирается случайным образом, т.е. любое значение из этого интервала равновероятно, то и достижение результата по точности вычислений корня возможно на любом шаге итерации, а значит, выполнение алгоритма может продолжаться сколь угодно долго.

На практике чаще применяют разновидность метода Монте-Карло, ограничивая выполнение алгоритма числом шагов итераций (например, N = 1 000, 10 000, 1 000 000 и т.д. – в зависимости от конкретных условий), после выполнения которых результатом считают ту точку х, в которой функция F=f(x) имеет меньшее значение по сравнению со всеми другими, полученными в ходе испытаний.

Алгоритм метода Монте-Карло по числу испытаний:

1.Установим число испытаний N.

2.Установим границы интервала изоляции a и b.

3.Генерируем x случайным образом на интервале [a, b].

4.Вычисляем значение функции F=f(x) в случайной точке x.

5.Установим счетчик испытаний I=1.

6.Если I>N, то идти к пункту 10.

7.Генерируем x1 случайным образом на интервале [a, b].

8.Вычисляем значение функции F1=f(x1) в случайной точке x1.

9.Если |F|≥|F1|, то F=F1 и х=х1 (запомним новые значения).

10.I=I+1.

11.Идти к пункту 6.

12.Искомый корень равен х и функция равна F.

13

3. ВАРИАНТЫ ЗАДАНИЙ

Найти с точностью 0,000001 наименьший корень уравнения. Перед применением указанного в варианте метода провести отделение корней. Сделать проверку найденного решения.

Варианты заданий приведены в табл. 3.1.

Таблица 3.1

Варианты заданий

Метод	Вариант								Уравнение
Половинного деления	1	2 ln( x ) − 1 +0,5 =0
								x
	6		x −1 = cos( 0,5 x )
	11		x +1 = 1
						x
	16		E−x =0,5 +							x
	21		x4 + 2 x3 − x −1 = 0
	26		2 sin		2	( x )		−		3 cos	2	( x )		= 0
			3							4
Простых итераций	2	0,2 x3 −0,6 x2 −1,6 x +3,7 =0
	7	3 x − E x =0
	12	2 x −ln( x ) −7 =0
	17	3 x −sin( x ) =7
	22		x3 +12 x − 2 = 0
	27		x −	sin( x )			−1 = 0
					2
Касательных	3		x3 −3 x2 +6 x −5 =0
	8		x6 + x2 −3 =0
	13	3 x2 −cos( 2 x ) −1 =0
	18		x3 +3 x +2 =0
	23	tg( x ) − x = 0
	28		x3 −6 x2 + 20 = 0

		14
		Продолжение табл. 3.1
Метод	Вариант	Уравнение
Хорд	4	x5 − x −0,2 =0
	9	x − sin( x ) =0,25
	14	x3 + x2 −3 =0
	19	x8 −0,4 x3 −1,24 =0
	24	x2 −1,3 ln( x +0,5 ) −2,8 x +1,15 =0
	29	x2 − sin( 5 x ) = 0
Монте-Карло	5	x4 −6 x3 +9 x2 −16 =0
	10	x3 −2 x −5 =0
	15	x =1 +cos( x )
	20	x +ln( 1 + x ) =1,5
	25	1,8 x4 − sin( 10 x ) =0
	30	5 x3 +10 x2 + 5 x −1 = 0

СПИСОК ЛИТЕРАТУРЫ

1.Инженерные расчеты на ЭВМ: Справ. пособие / Под ред. В.А.Троицкого. – Л.: Машиностроение. Ленингр. отд-ние, 1979. – 288 с.

2.Тихонов А.Н. Вводные лекции по прикладной математике / А.Н.Тихонов, Д.П.Костомаров. – М.: Наука. Гл. ред. физико-математ.

лит., 1984. – 192 с.

3.Васильков Ю.В. Компьютерные технологии вычислений в математическом моделировании: Учеб. пособие / Ю.В.Васильков, Н.Н.Василькова. – М.: Финансы и статистика, 1999. – 256 с.

4.Дьяконов В.П. Справочник по алгоритмам и программам на языке Бейсик для персональных ЭВМ. – М.: Наука. Гл. ред. физико-

математ. лит., 1987. – 240 с.

5.Мак-Кракен Д. Численные методы и программирование на Фортране / Д.Мак-Кракен, У.Дорн – 2-е изд., стереотипное. – М.: Мир, 1977. – 584 с.

6.Боглаев Ю.П. Вычислительная математика и программирование: Учеб. пособие для студентов втузов. – М.: Высш. шк., 1990. – 544 с.

15

	ОГЛАВЛЕНИЕ
1.	Цель работы.............................................................................................	1
2.	Теоретические сведения.........................................................................	1
	2.1. Основные положения.......................................................................	1
	2.2. Отделение корней нелинейного уравнения....................................	2
	2.3. Вычисление корней с заданной точностью....................................	4
3.	Варианты заданий.................................................................................	13
Список литературы.....................................................................................		14

16

Александр Вениаминович Матисов Галина Алексеевна Алексеева Виктор Иванович Штефан

РЕШЕНИЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙ

Методические указания к лабораторной работе по курсу "Информатика" для студентов специальности 210200 "Автоматизация технологических процессов и производств (в машиностроении)", 071900 "Информационные системы и технологии", 120100 "Технология машиностроения", 120200 "Металлорежущие станки и инструменты", 120500 "Оборудование и технология сварочного производства"

Редактор З.М.Савина

Подписано в печать 21.11.02. Формат 60х84/16. Бумага офсетная.

Отпечатано на ризографе. Уч.-изд. л. 1,00. Тираж 75 экз. Заказ ГУ КузГТУ 650026, Кемерово, ул. Весенняя, 28.

Типография ГУ КузГТУ 650099, Кемерово, ул. Д.Бедного, 4 А.