Файл: Ермолаева Э.Н. Элементы численного анализа учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 07.07.2024
Просмотров: 129
Скачиваний: 0
Д3 — погрешность округления результата; Л, = 13,8—13,79 = 0,01;
А= 0,051 + 0,0005 + 0,01« 0,062.
Так как |
Л >0,05, то последняя цифра |
результата 8 — со |
|||||||||||
мнительная. |
|
|
|
|
|
|
|
|
|
||||
Продолжаем |
округление |
результата: |
13,8~ 14. Новая по |
||||||||||
грешность |
есть сумма |
прежней |
|
погрешности |
и погрешности |
||||||||
этого |
округления: |
|
|
|
|
|
|
|
|
||||
Д' |
= |
0,062 + |
(14 - |
|
13,8) = |
0,062 + 0,2 ^ |
0,27 ^ 0,3 . |
||||||
Так как |
Д'<0,5, цифра 4 — верная. |
|
|
||||||||||
Результат: 14 ±0,3. |
|
|
|
|
|
|
|
||||||
д) Вычисляем относительную погрешность, используя бо |
|||||||||||||
лее точное значение абсолютной |
погрешности: |
|
|||||||||||
|
|
|
|
|
|
|
|
0 'Л |
|
|
|
|
|
|
|
|
|
|
|
|
|
- |
~ |
0,02 = 2«о • |
|
||
|
|
|
|
|
|
|
|
14 |
|
|
|
|
|
2 . 0,862 - 7,15832: |
|
|
|
|
|
|
|
||||||
а) |
|
Оставляем |
|
без изменения 0,862, так как в нем 3 верные |
|||||||||
значащие цифры, а во втором множителе — пять. |
|||||||||||||
б) |
|
Округляем |
7,1583 до четырех значащих |
цифр: 7,1583^ |
|||||||||
«7,158. |
|
|
|
|
|
|
|
|
|
|
|
||
в) |
|
Производим вычисления |
и в результате тоже оставляем |
||||||||||
3 значащие цифры: 0,8621582 ~44,2. |
|
|
|||||||||||
г) Вычисляем относительную |
погрешность |
|
|||||||||||
|
|
|
|
|
|
|
|
3 = |
Sj + |
о , , |
|
|
|
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
З; |
|
— относительная |
погрешность числа 0,862; |
||||||||||
° 1 |
|
0,0005 |
_ |
|
|
|
|
|
|
|
|||
= |
~0"862~ |
О'ШЪЭ ; |
|
|
|
|
|||||||
5, |
|
— относительная погрешность числа 7.1582, равная удво |
|||||||||||
|
|
|
енной погрешности числа 7,158. |
|
|
||||||||
Абсолютная погрешность числа 7,158 равна сумме погреш |
|||||||||||||
ности |
|
исходного |
|
числа |
7,1583 и погрешности округления: |
||||||||
|
|
|
|
|
|
|
|
0,00005 + 0,0003, |
|
|
|||
таким |
|
образом, |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
8 |
, |
, |
2 |
• Д > ° 0 ° 3 5 |
« 0,000098 ; |
||||
|
|
|
|
|
о ж |
0,00069^0,0007 = 0,07%. |
|
д) Зная о , вычисляем Д:
А= 0,00069-44,2 «0,03.
А<0,05, следовательно, 2 — цифра верная. Результат: 44,2 ±0,03.
|
|
3. 83,112 |
: 2,88 « 83,11 |
: 2,88 ^ |
28,9. |
||
- |
- - х |
0,0005 |
+ 0,002 |
, |
0,005 |
|
|
о - |
о , + 6 3 |
= |
|
то |
-г• - T e g - - 1-0019-^0,19% |
||
|
|
Д = |
0,0019 • 28,9 |
0,055 . |
|||
Так как Д'>0,05, последняя |
цифра результата 9 — сомни |
||||||
тельная. |
|
|
|
|
|
|
|
|
28,9 |
^ 29; |
Д' = 0.055 4- 0,1 |
1С ^ 0,2 . |
А' <0,5, следовательно, в числе 29 все цифры верные. Результат: 29±0,2.
Г л а в а 2
ПР И Б Л И Ж Е Н Н О Е РЕШЕНИЕ АЛГЕБРАИЧЕСКИХ
ИТРАНСЦЕНДЕНТНЫХ УРАВНЕНИЙ
§ 2.1. П О С Т А Н О В К А З А Д А Ч И
Для некоторого уравнения (алгебраического или транс цендентного) существует формула решения, если его корни можно выразить через входящие в уравнение величины при по мощи арифметических операций, извлечения корней, показа тельной, логарифмической, тригонометрических и обратных тригонометрических функций. В этом смысле у квадратного уравнения x2 + px + q = 0 есть формула решения, имеющая вид:
- - - л - * V / J f ^ -
Есть формула и для решения кубического уравнения хг + рх + д = 0. Она имеет вид:
Однако практическое применение этой формулы наталки
вается на ряд трудностей и требует использования |
комплекс |
ных чисел. |
|
Существует формула и для решения уравнений |
четвертой |
степени, но настолько сложная, что приводить ее |
не будем. |
Однако доказано, что в общем случае при п> 5 не существует
формулы, выражающей |
решение алгебраического уравнения |
а0 хп 4 ал |
х"-1 + . . . 4- а„- = О |
при помощи арифметических операций и извлечения корней.
22
Кроме того, в практических инженерных задачах чаще все го уравнение содержит коэффициенты, известные лишь при близительно, и тогда сама задача о точном определении корней уравнения теряет смысл. Да и практически точное значение корня знать не надо, а надо знать его значение с определенной, заранее заданной степенью точности.
Поэтому существуют способы приближенного нахождения корней уравнения и оценки степени их точности. Большая часть этих способов основана на идее последовательных приближе ний.
Будем рассматривать уравнение
/ ( л - ) - 0 . . |
(2.1) |
где функция |
/ |
У /( . V) |
(2.2) |
существует и непрерывна в некотором отрезке.
Корнем (или нулем) уравнения (2.1) называют всякое зна
чение с, обращающее функцию |
(2.2) в нуль, т. е. |
|
|||
|
|
/ ( = ) |
= 0 . |
|
(2.3) |
В случае нахождения |
действительных |
корней |
равенство |
||
\2.3) геометрически |
означает, |
что график |
функции |
(2.2) про |
|
ходит через точку с |
оси |
ОХ. |
|
|
|
Приближенное нахождение действительного корня урав нения (2.1) состоит в том, что надо сначала отделить этот ко рень от других, т. е. найти Такой отрезок [а; Ь], в котором со держится только один этот корень уравнения, и потом уточ нить приближенное значение этого корня, т. е. найти его зна чение с заданной степенью точности.
Графический метод
Для отделения корней можно использовать график функ ции (2.2). Абсциссы точек пересечения этого графика с осью абсцисс являются корнями уравнения (2.1). Построение гра фика даже с малой точностью дает обычно представление о расположении интересующих нас корней и тем самым позво ляет отделить их. Если построение указанного графика вызы вает затруднения, то следует преобразовать исходное уравне ние к виду ?i(jr) = c2 (jc) так, чтобы графики функций \>=ъ,(х) и _у ="-<?2(JC) можно было легче построить. Абсциссы точек пе ресечения этих графиков и будут корнями рассматриваемого уравнения.
Для примера решим кубическое уравнение.
х3— \,75х + 0,75 = 0.
Перепишем уравнение в виде я 3 =1,75%—0,75. Тогда иско мые корни найдутся как абсциссы точек пересечения параболы у = х3 с прямой у=\,7Ьх—0,75. Построив графики параболы и прямой, из рис. 2.1 находим, что уравнение имеет три действи тельных корня: Х\ = —1,5; ^2 = 0,5; х 3 = 1.
Р и с . 2.1.
Несмотря на то, что графические методы решения уравне ний удобны и просты, они применимы лишь для грубого опре деления корней. Особенно неблагоприятным в смысле потери точности является случай, когда линии пересекаются под очень острым углом и практически сливаются по некоторой дуге.
Тогда применяют метод проб.
§ 2.2. М Е Т О Д П Р О Б
Метод проб основан на теореме: если непрерывная функция y=f(x) принимает значения разных знаков на концах отрезка fa; b], т. е. f(a)f(b) <0, то внутри этого отрезка содержится по меньшей мере один корень £ уравнения f(x) — 0.
Если существует f'(x) и сохраняет постоянный знак в [а; Ь], то корень q будет единственным (рис. 2.2 и 2.3).
Процесс отделения корней надо начинать с определения знака функции (2.2) в граничных точках отрезка [а; Ь]. Если полученные знаки противоположны, надо выбрать промежу точную точку Хо и определить знаки функции на концах отрез ков [а; Хо] и [х0; Ь). Корень будет в отрезке, на концах которого
3
и
Рис. 2.4
знаки |
функции |
противоположны |
и т. д. |
(рис. 2.4), |
при этом |
хо, х\,... |
можно |
брать более или |
менее |
произвольно, |
исходя |
из удобства вычислений. О единственности корня будем судить
по тому, меняет |
знак |
f'(x) в последнем рассматриваемом от |
|||||||||
резке или нет. |
|
|
|
|
|
|
|
|
|
|
|
Для отделения корней практически часто бывает достаточ |
|||||||||||
но делить каждый раз промежуток |
пополам |
(метод |
половин |
||||||||
ного деления). |
|
|
|
|
|
|
|
|
|
|
|
Если в одном из промежутков |
|
|
|
|
|
|
|
||||
|
|
|
или |
а |
+ b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
находится корень, то обозначим этот промежуток |
[а\\ Ь{\. Но |
||||||||||
вый промежуток [аі; b\] снова делим |
пополам н находим |
отре |
|||||||||
зок [й2 ; Ь2], в котором |
находится корень и т. д. |
На л-м |
шаге |
||||||||
получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
On |
|
|
|
|
|
|
|
|
Так как последовательность а„ монотонно неубывающая и |
|||||||||||
ограничена числом b, |
a b,t — последовательность |
|
монотонно |
||||||||
невозрастающая |
и ограничена |
числом |
а, то существуют |
пре |
|||||||
делы Игл ап |
и |
iim b„, а из последней |
формулы |
следует, что |
|||||||
эти пределы равны искомому корню |
с. |
|
|
— ап |
|
|
|||||
Очевидно, |
что абсолютная |
погрешность |
; |
удовлет |
|||||||
воряет неравенству |
|
|
|
|
|
|
|
|
|
||
|
|
О < ? - а„ < ~ ~ (Ь -а) . |
|
|
|
|
|||||
Пример 2.1. Методом половинного деления уточнить корень |
|||||||||||
уравнения |
|
|
|
|
|
|
|
|
|
|
|
|
fix) |
2х[) |
х |
— 1 |
О |
|
|
|
|
лежащий на отрезке [0; 1].
Р е ш е н и е . Последовательно вычисляем:
/(0) = - |
1; / ( 1 ) = |
1; |
/ ( 0 , 5 ) |
|
= - 1 , 1 9 ; |
/ ( 0 , 7 5 ) . - - 0,59; |
|
/ (0,875) - |
0.05; |
/ |
(0,8125) = |
- 0,304; |
|
f |
(0,8438) = |
- |
0,135; |
/ |
(0,8594) = |
— 0.Г43. |
Принимаем |
1 |
|
|
|
|
|
|
|
(0,859 + 0,875) =0,867. |
Промежуточных выкладок в этом примере не показали, но видно, что в общем случае для подсчета корня -с достаточно