Файл: Ермолаева Э.Н. Элементы численного анализа учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 07.07.2024
Просмотров: 133
Скачиваний: 0
|
- |
1 |
> |
|
|
|
Если |
q |
2 |
1 |
0 |
|
|
|
|
I |
с |
- |
x„ і < дгл - х„ , І . |
(2.14) |
Итак, |
формулы |
(2.13) и (2.13') имеют место при |
| 'f' (л') | < |
|||
ц < 1 , а формула |
|
(2.14) — только при |
|
Вывод. Если при ! 'f' (х) | < в окрестности корня ; два
последовательных значения х„ і и хп совпадают между собой с заданной точностью г (например, для этих приближений уста новились т первых десятичных знаков), то.с той же точностью
справедливо равенство с ^ |
хп. |
|
|
|
|
Однако при -g- < |
| ъ' (х) |
| < 1 |
из неравенства |
; х„ — х„.. \ \ < s |
|
может не следовать неравенство |
| с — хп\ < г. |
Поэтому |
урав |
||
нение (2.1) желательно преобразовать в уравнение (2.8) |
так, |
||||
чтобы Ф' (х) | < ~ |
в окрестности корня. |
|
|
||
Так, в примере 2.6 |
|
|
|
|
3 1 (А- ; 2)-
при АГ є [1; 2], поэтому оценка погрешности произведена верно. Если уравнение j(x) = 0 переписать в виде
f(x) |
(b — x) |
, т, е. |
X = X — - ' |
|
f ( b ) - f ( x )
то формула (2.9) примет вид (2.4).
Если же уравнение (2.1) переписать в виде
х = х |
f, <*>- , т. е. ? (х) = ? ї ( д 0 = * - - Д * > - - |
го формула (2.9) превратится в формулу (2.7).
Вывод. Методы хорд и касательных являются частными случаями метода итераций.
Поэтому все выводы настоящего параграфа остаются спра ведливыми для методов хорд и касательных.
|
Так, |
если |
в окрестности |
корня |
\'i'i(x) |
j < |
, то |
в |
смысле |
|
быстроты сходимости выгоднее применять метод хорд. |
||||||||||
|
В примере 2.2 |
|
|
|
|
|
|
|
||
. л _ |
_ |
1 - |
х |
, . |
_ |
(х' + З х - 1 ) ( 1 - х ) _ |
||||
|
|
г 1 - 1 |
|
_ |
З I 2x3 |
- Зх2 |
+ 1 I |
1 |
||
|
x3 |
+ 3x — 4 ' |
• 1 v ; |
|
(x3 + 3x - 4)2 |
|
^ 2 |
|||
і:ри |
x є |
[0; 1] и поэтому заключение о точности найденного кор |
||||||||
ня |
верно. |
|
|
|
|
|
|
|
|
|
|
Если |
же |
в окрестности |
корня |
\ъ'.,(х) |
J *?"2~> т о |
выгоднее |
применять метод касательных. Из выражения для производ ной
_£{х)Г_{х)_
следует, что метод касательных при одном и том же числе ша-
юв п тем точнее, чем меньше |
\f'(x)\ |
в окрестности корня, что |
|||
согласуется с выводами, сделанными |
ранее. |
|
|||
В примере |
2.3 |
|
|
|
|
, , , , . |
2 і х (3я — Зх - 5) |
^ 1 |
|||
1?»(*>1 = |
' |
3 ( * » - Ц ' |
|
||
при х є [0; 3], |
поэтому |
корень с заданной |
точностью найден |
||
верно. |
|
|
|
|
|
В примере 2.4 для уравнения / ( х ) =? хк — а = 0 имеем
|
» / \ |
|
k |
\ |
X. |
ее |
' |
(Х) |
= |
г |
|
|
. — , |
- К ' |
|
k |
|
|
хк |
|
откуда 'f'2 (х) % 0 |
при |
хк—а^О, |
|
т. е. вблизи корня; поэтому |
метод Ньютона — самый лучший для приближенного извлече ния корней, причем абсолютная погрешность вычисляется по формуле (2.14).
Необходимо отметить, что при нахождении корней уравне ния рассмотренными методами считаем коэффициенты уравне ния числами точными. Поэтому корни уравнения могут быть найдены с любой степенью точности в том смысле, что при под становке их в решаемое уравнение (2.1) можно добиться тою,
что левая часть будет практически сколь угодно мало отли чаться от нуля.
§ 2.7. В Ы П О Л Н Е Н И Е Л А Б О Р А Т О Р Н О Й Р А Б О Т Ы
Рассмотрим пример. Найти корень уравнения 3 lgx + х2—6 = 0 с точностью до четвертого десятичного знака методом итера ций.
У
Ри с . 2.16
1.Сначала определим грубо приближенное значение корня ірафическим методом. Для этого данное уравнение запишем в виде 3 lgx = 6—х2 и по нескольким точкам построим графики
кривых yi=3\gx |
и г/2 = 6—х2, учитывая, что из условия следует |
||||
х > 0 (табл. 2.2 |
и рис. 2.16). |
|
|
|
|
|
|
|
|
|
Таблица 2.2 |
X |
0 |
1 |
2 |
у'"П) |
10 |
Уі |
— оо |
0 |
0,9 |
1,5 |
3 |
Уа |
6 |
5 |
2 |
— 4 |
— 94 |
|
|
|
|
|
Абсциссу точки пересечения кривых, найденную грубо цз графика рис. 2.16: Xo~2,3, считаем первым приближением ис комого корня.
2. Приведем данное уравнение к виду |
х=ъ(х) |
|
и проверим, |
||||||
|
|
|
. // |
м |
1 |
|
|
|
|
выполняется |
ли условие | ? (Л:) | < — в окрестности корня: |
||||||||
х = | / |
6 - |
3 lg х ; |
© (х) = | / 6 |
3 Ig ж ; |
|||||
|
|
' |
(х) = |
3 • 0.43 |
|
|
|
|
|
|
|
•2хУ |
Є — |
3\gx |
|
|
|
||
|
|
|
|
|
|
|
|||
?' |
(2, |
3) |
|
3 |
• 0,43 |
|
|
|
|
2 • 2,3 | / 6 - |
3 lg 2,3 |
|
|
||||||
|
|
|
|
|
|||||
3. Уточним приближенное |
|
|
|
|
0 |
~2,3 но фор |
|||
муле (2.9): |
|
|
|
значение |
корня |
л: |
|
< |
|
хп |
= |
у 6 — 3 lg л:„_ |
и = 1, 2, 3 , . |
|
|
Так как требуется вычислить корень с точностью до четвер того десятичного знака, то согласно формуле (2.14) процесс итерации надо закончить, когда у двух последовательных при ближений совпадут четыре десятичных знака.
Все промежуточные вычисления записываем в табл. 2.3 с пятью верными десятичными знаками. Значения lg* и значе ния корня находим с помощью таблиц логарифмов.
|
|
|
|
|
|
|
Таблица 2.3 |
|
|
|
|
6 — |
l g (6 - |
• l g ( 6 - |
|
хп 1 !g хп-\ 3 l S x n 1 - 3 i g * „ _ , |
|
- 3 1 g x „ |
, ) |
||||
|
|
|
|
|
|
||
2,3 |
J 0,36173 1,08519 |
4,0148 |
0,69150 |
0,34575 |
2,2169 |
||
2,2169! |
0,34575 |
1,03725 |
4,9628 |
0,69572 |
0,34786 |
2,2277 |
|
|
|
|
|
|
|||
2.227?! |
0,34786 |
1,04358 |
4,9564 |
0,69517 |
0,34758 |
2,2263 |
|
|
|
|
|
|
|||
2,2263 |
|
0,34758 |
1,04274 |
4,9573 |
0,6952.) |
0,34762 |
2,2265 |
|
|
|
|
|
|
||
2,22:5' |
0,34762 |
1,04286 |
4,9571 |
0,69523 |
0,31762 |
2,2265 |
|
|
|
|
|
|
Итак, лг4 = *5 = 2,2265. Следовательно, с точностью до 0,0001 число 2,2265 является корнем уравнения 3\gx + х2—6 = 0.
К О Н Т Р О Л Ь Н Ы Е В О П Р О С Ы |
|
1. Когда применяются приближенные |
методы решения |
уравнений? |
|
2. В чем суть графического метода? |
( |
3. |
Что такое метод проб, метод половинного деления? |
||||||
4. |
Сделать вывод |
формулы метода |
хорд в случае, если |
||||
/ ( а ) > 0 , / ( 6 ) < 0 , |
а<Ь. |
|
|
|
|
||
5. |
Будет |
ли |
формула метода |
хорд |
совпадать |
с формулой |
|
'2.4), если /"(*) <0,1(a) <0, f(b) |
>0? |
|
|
||||
6. Сделать вывод формулы метода касательных для случая, |
|||||||
если f(a)>0, |
f(b)<0. |
Какой знак должно принимать произве |
|||||
дение f(xo)f"(xo) |
для |
сходимости |
метода, где х0 |
— начальное |
|||
приближение? |
|
|
|
|
|
7.Когда применяют метод хорд, а когда — метод касатель
ных?
8.В чем суть комбинированного метода?
9.В чем суть итерационного процесса? При каких условиях З'Тот процесс сходится?
10.Как можно геометрически истолковать процесс итера ции? Какие случаи здесь могут быть?
11.Как оценивается погрешность метода итерации? Сде лать вывод соответствующих формул.
12. Как оценивается погрешность итерации в случае, если
і ?' |
(х) |
і < 4 - |
? |
|
13. Как уравнение f(x) |
= 0 можно привести к виду х |
— ъ (л)? |
||
Однозначно ли это делается? |
Как |
надо выбирать |
функцию |
(х)для быстрой сходимости метода итерации?
14.Существует ли связь между методами хорд, касатель ных и методом итерации?
15.Каким методом выбирают первое, грубое приближение
корня?
§ 2.8. М Е Т О Д Л О Б А Ч Е В С К О Г О Р Е Ш Е Н И Я |
* |
А Л Г Е Б Р А И Ч Е С К И Х У Р А В Н Е Н И Й |
|
Предыдущая группа методов разыскания приближенных корней алгебраических уравнений предусматривала, прежде всего, отделение корней. Этот процесс достаточно затрудните лен в случае большого числа корней и особенно при наличии комплексных корней. В большинстве же практических вопро сов, связанных с электро- и радиотехникой, нужно знать не только вещественные, но и комплексные корни.
Метод Лобачевского предоставляет такую возможность, будучи достаточно простым в вычислительном отношении и применимым к любому алгебраическому уравнению.
I. Идея метода. Рассмотрим уравнение
/ (х) == Xя + а, ж»-1 + . . . + ап = 0 , |
(2.15) |
имеющее только вещественные корни, причем
46
'•v., > -v, > . . . > ] А „ | . (2.16)
Предположим, |
|
что по |
уравнению |
(2.15) |
удалось |
найти |
|
уравнение |
|
|
|
|
|
|
|
|
/ • ( А ) .-== |
A " -f- Л, х" 1 + . . . 4- Ап |
= 0 |
(2. і7) |
|||
с корнями |
|
|
|
|
|
|
|
|
|
— |
X ' " . — |
А''" |
А " " , |
|
|
где /// — достаточно |
велико. |
|
|
|
|||
Коэффициенты |
и корни |
уравнения |
(2.17) |
связаны |
соотно |
||
шениями |
Виета: |
|
|
|
|
|
|
Л, ; X'" |
-\- X? f |
. . . 4 X™ |
— |
А., |
1 |
|
|
|
|
|
|
|
|
|
|
• + |
4 |
|
|
|
|
|
|
|
|
|
Л2 |
« |
(А'| |
х2)" |
1 + |
|
|
|
|
|
|
|
A'J |
|
|
|
|
|
|
|
|
|
|
(2.18) |
As |
= |
( А , А , . . . А , ) " ' 1 4- |
|
|
|
||
|
|
|
|
+ |
Xn—s |
• • • |
Хп |
|
|
|
|
...+ - |
. |
х. |
|
А„= |
( А , |
А , . |
..Х„Г . |
|
|
|
|
|
Если //і — достаточно велико, то, так как корпи пронумеро |
||||||
ваны |
в порядке убывания их абсолютных |
величин, дроби в |
квадратных скобках значительно меньше единицы, и прибли женно можно положить:
Л, ^ л'" ;
Л, » (xj х,)т ;
(2.19)
As » (л, х , . . . л-,.)'";
Ап =-- (х, х2 . . . хп)>" .