Файл: Ермолаева Э.Н. Элементы численного анализа учеб. пособие.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)

[)

х

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.

Промежуточных выкладок в этом примере не показали, но видно, что в общем случае для подсчета корня -с достаточно