Документ: Лабораторные (заочн), страницы 1-5

Тема 2. Лабораторная работа Методы решения нелинейных уравнений

2.1.Вопросы, подлежащие изучению

1.Постановка задачи численного решения нелинейных уравнений.

2.Этапы численного решения уравнения.

3.Аналитический и графический методы отделения корней.

4.Уточнение корня методами половинного деления, итерации, Ньютона и хорд.

5.Графическая иллюстрация методов половинного деления, итерации, Ньютона и хорд.

6.Условие окончания вычислений при использовании методов половинного деления, итерации, Ньютона и хорд.

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

8.Теорема о сходимости метода Ньютона и оценка погрешности метода.

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

10. Условия окончания вычислений в методах итерации, Ньютона и хорд.

11.Сравнение методов половинного деления, итерации, Ньютона и хорд.

12.Алгоритмы и программы решения нелинейных уравнений на языке программирования.

13.Решение нелинейных уравнений средствами математических пакетов

2.2.Задание

1.Выбрать индивидуальное задание (нелинейное уравнение) из таблицы.2-1.

2.Отделить корни уравнения.

3.Уточнить корень уравнения 4-мя численными методами, для чего провести исследование нелинейного уравнения для каждого метода:

проверить выполнение условий сходимости вычислительного процесса, в случае расходящегося процесса – сделать необходимые преобразования для обеспечения сходимости;

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

сформулировать условия окончания этапа уточнения корня.

4.Провести «ручной расчет» трех итераций.

5.Оценить погрешность результата «ручного расчета».

6.Решить нелинейное уравнение средствами математического пакета.

Тема1.2.Методы решениянелинейных уравнений (Лабораторныйпрактикум)

Страница 1

2.3. Варианты задания

										Таблица 2-1
№	Уравнение			№	Уравнение
1	tg(0,36x +0,4) = x2			19	ex + ln(x) – x = 0
2	x + ln(4x) – 1 = 0			20	1–x+sin(x)–ln(1+x) = 0
3	ex – 4 e-x – 1 = 0			21			– cos(1–x) = 0
						1 x
4	x ex – 2 = 0			22	sin(x2)+cos(x2)–10x = 0
5	4 (x2 + 1) ln(x) – 1 = 0			23	x2 – ln(1 + x) – 3 = 0
6	2 – x – sin(x / 4) = 0			24	cos(x / 2) ln(x – 1) = 0
7	x2 + ln(x) – 2 = 0			25	cos(x/5)					– x = 0
									1 x
8	cos(x)–(x + 2)1/2 + 1 = 0			26	3x – e-x = 0
9	4 (1 + x1/2) ln(x) – 1 = 0			27	4(1+
								x) ln(x)–10 = 0
10	5 ln(x) –			28	(x-3)2 lg(x-2) = -2
			x = 0
11	ex + x3 – 2 = 0			29	x – 1 / (3 + sin(3.6x)) = 0
12	3 sin (			30	0.25x3 + cos(x / 4) = 0
		x) + x – 3 = 0
13	0.1x2 – x ln(x) = 0			31	2 – x = ln x
14	cos(1 + 0.2x2) – x = 0			32	x2 + 4sinx = 0
15	3 x – 4 ln(x) – 5 = 0			33	0,5x -3 = (x+2)2
16	sin(1 – 0.2x2) – x = 0			34	1 + lgx = 0,5
17	ex – e-x – 2 = 0			35	2 lgx – x/2+ 1 = 0
18	x – sin(1 / x) = 0			36	x – sinx = 0,25

2.4.Содержание отчета

1.Индивидуальноезадание (уравнение, метод решения).

2.Результаты этапа отделения корней (интервалы изоляции корня уравнения).

3.Результаты исследования задания для «ручного расчета»:

условие сходимости вычислительного процесса;

начальное приближение;

условие окончания процесса уточнения.

4.Результаты «ручного расчета», представленные по форме табл. 2.-2а для метода половинного деления или по форме табл. 2.-2б для остальных методов.

						Таблица 2-2а
к	a	b	f(a)	f(b)	(a+b)/2	f( (a+b)/2)	b-a
1
2
3
4

		Таблица 2-2б
к	x	f(x)

1

2

3

4

5. Оценкипогрешностей результатов «ручного расчета».

6. Результаты решения задачи, полученные с помощью математического пакета. Страница 2

Тема1.2.Методы решениянелинейных уравнений (Лабораторныйпрактикум)

2.5.Пример выполнения задания

1.Задание для решения нелинейных уравнений:

уравнение f(x) 1 3x cosx 0;

методы решения нелинейных уравнений для ручного расчета – половинного деления, итерации, Ньютона и хорд;

2.Отделение корней с использованием MathCad

Отделение корней производим графическим методом (а) с обязательным подтверждением результата аналитически (б)

а)

f (x) 1 3x cos(x)

							10
							10
									5






f(x)		2						5			0			2		4



						10


						15						x
						15

f1(x) d			f (x)							sin(x)					3
	dx
f2(x)		d	2		f (x) cos(x)
		d
				2
		dx
x 0 0.1 1
x				f(x)									f1(x)				f2(x)

0							2							-3			-1

0.1							1.695							-3.1			-0.995

0.2							1.38							-3.199			-0.98

0.3							1.055							-3.296			-0.955

0.4							0.721							-3.389			-0.921

0.5							0.378							-3.479			-0.878

0.6							0.025							-3.565			-0.825

0.7				-0.335										-3.644			-0.765
0.8				-0.703										-3.717			-0.697

0.9				-1.078										-3.783			-0.622

1							-1.46							-3.841			-0.54

б) На отрезке [0; 1] функция f(x) меняет знак, т.е. существует, по крайней мере, один корень. Поскольку знак первой производной f1(x)= -sin(x) - 3 < 0 на выбранном отрезке остается постоянным, то можно сказать, что функция на этом отрезке монотонна. Знакопостоянство второй производной f2(x)= -cos(x)<0 на выбранном отрезке является необходимым условием применения метода Ньютона и метода хорд. Следовательно, уравнение 1-3х+cos(x)=0 имеет единственный корень на отрезке [0;1].

Тема1.2.Методы решениянелинейных уравнений (Лабораторныйпрактикум)

Страница 3

3.Уточнение корня с использованием MathCad

Метод половинного деления

1) Исследование задания

Метод половинного деления сходится, если на выбранном отрезке отделен один

корень. Так как на отрезке [0;1] функция f(x) 1 3x cosx 0 меняет знак (f(0) f(1) 0) и монотонна (f’(x)<0), то условие сходимости выполняется.

Выберем за начальное приближение середину отрезка	x	0	a b	=0.5.
		0	2

Условие окончания процесса уточнения корня. Для оценки погрешности метода половинного деления справедливо условие |bn – an|<ε , т.е. длина отрезка, полученного на n-ом шаге должна быть меньше заданной точности - b2na

2)Результаты «ручного расчета» трех итераций

f(x) 1 3x cos(x)

1 итерация

a

b

a 0

b 1

x0

x0 0.5

f(x0)=0.377

2

f(a)>0,

f(x0)>0 и f(b)<0

следовательно, a1 x0

b1 b

2 итерация

a1 b1

a1 0.5

b1

1 x1

x1 0.75

f(x1)=-0.518

2

f(a)>0,

f(x1)<0 и f(b)<0

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

a2 a1

b2 x1

3 итерация

a2 0.5

b2

0.75

x2

a2 b2

x2 0.625

f(x2)=-0.064

2

f(a)>0,

f(x2)<0 и f(b)<0

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

a3 a2

b3 x2

и т.д.

Результаты вычислений представлены в форме табл. 2-2а.

n

a

b

f(a)

f(b)

(a+b)/2

f( (a+b)/2)

b-a

1

0

1

2

-1.459

0.5

0.377

0.5

2

0.5

1

0.377

-1.459

0.75

-0.518

0.25

3

0.5

0.75

0.377

-0.518

0.625

-0.064

0.125

4

0.5

0.625

0.377

-0.064

0.5625

0.158

0.0625

После трех итераций приближение к корню – середина отрезка x3=0.5625.

3) Погрешность численного решения нелинейных уравнений

Оценим погрешность результата после трех итераций

b3 a3

0.125.

Страница 4

Тема1.2.Методы решениянелинейных уравнений (Лабораторныйпрактикум)

Метод итераций

1) Исследование задания для «ручного расчета»

Приведем

уравнение

f(x)=0

к виду x (x). Тогда

рекуррентная формула

xn 1 (xn),

n 0,1,....

Для сходимости процесса итерации необходимо, чтобы

'(x)

1,

при x [a;b]. Если

'(x)

1, то сходимость не обеспечена.

Приведем

уравнение f(x) 1 3x cosx 0 к виду x =

(cos(x)+1)/3 и проведем

исследование.

(x)

1 cos(x)

' (x) d (x)

3

dx

x 0 0.1 1

x

' (x)

0	0	0	0.2	0.4	0.6	0.8	1
0.1	-0.033	0	0.2	0.4	0.6	0.8	1
0.1	-0.033
0.2	-0.066	0.1
0.3	-0.099	' (x)
0.3	-0.099	0.2
0.4	-0.13	0.2
0.4	-0.13
0.5	-0.16	0.3
0.6	-0.188	0.3
0.6	-0.188
0.7	-0.215				x
0.8	-0.239	'(x)<1 для всех значений аргумента
0.9	-0.261	'(x)<1 для всех значений аргумента
0.9	-0.261	на отрезке [0,1]
		на отрезке [0,1]

1-0.28

Вприведенном примере условие сходимости выполняется и можно использовать

итерирующую функцию	(x)	1 cosx	в рекуррентной формуле для уточнения корня
		3

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

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

Построим функцию (x) λf(x) x, где параметр λ может быть определен по правилу:

если f (x) 0,то

1

λ 0;

r

1,

если f (x) 0,

то 0 λ

где r max(

f (a)

;

f (b))

.

r

x [a;b]

Приведем пример выбора параметра λ и итерирующей функции. Для заданного уравнения f(x) 1 3x cosx 0исследовано, что

f `(x)<0, тогда 0 1.
r
f `(0)= -3, f `(1)=-3.841,
r =max{ \|-3\|, \|-3.841\| }=3.841, тогда 0<λ<0.26.
Полагаем λ=0.25. Тогда рекуррентная формула xn+1 = φ(xn),
где φ(x)= 0.25(1 - 3x + cos x) + x = 0.25(1+x+cos x).	Страница 5
Тема1.2.Методы решениянелинейных уравнений (Лабораторныйпрактикум)	Страница 5

Материал: Лабораторные (заочн)

Смотрите также: