Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Пулькин, С. П. Вычислительная математика пособие для учащихся 9-10 классов по факультативному курсу

.pdf
Скачиваний:
17
Добавлен:
21.10.2023
Размер:
9.64 Mб
Скачать

4) Барометрическая формула (рис.

10).

Здесь h— высота над уровнем

моря,

р0—давление на

уровне моря («нормальное»),

р —давление на данной

высоте.

мы рассмотрели два вида

Погрешность формулы. В § 4

погрешностей: неустранимая погрешность и погрешность округления. Теперь мы встретили еще один вид погреш­ ности; погрешность формулы. Она происходит от того, что некоторую точную формулу мы заменяем другой, неточной, но более удобной для вычисления. Если мы пользуемся приближенной формулой, то полная погрешность вычисли­ тельного процесса определяется, таким образом, в ре­ зультате влияния трех видов погрешностей: неустранимой погрешности, погрешности округления, погрешности фор­

мулы.

____

 

П р и м е р .

Вычислить a — l ^ l + x при

х = 0,186 по

приближенной формуле (2.9),

 

 

5000м

405мм

 

3000м

52Вмм

Применяя формулу, получаем:

К Т Т М « 1,093.

Вычисление по приближенной формуле мы провели точно. Учтя погрешность формулы (см. таблицу V), ви-

дим, что Д < -g-, т. е. Д < 0,005. Значит, в полученном

результате верны только две цифры после запятой. Отбрасывая цифру 3, мы округляем результат, полу­ чаем 1,09. Погрешность этого окончательного результата равна сумме погрешности округления (0,003) и погреш­ ности формулы (0,005). Итак, полная погрешность не превышает 0,008, т. е. меньше 0,01.

Эти расчеты мы провели,

предполагая, что х = 0,186

есть точное число. Если же

указанное значение х есть

приближенное число, то следовало бы еще учесть по­ грешность этого числа.

У праж нения к гла ве 2

1. Дана функция у —

З д;2 _ !_ 7

• Найти значение у при

х = 3,142, сохраняя

в ответе лишь верные цифры. Оце­

нить абсолютную и

относительную погрешности резуль­

тата. Значение х считать приближенным.

2. Составить таблицу значений функции y==-JL^p^L

для значений х: 1,0; 1,2; 1,4; 1,6; 1,8; 2,0. В ответе сохранить лишь верные цифры. Значения х считать точ­ ными, коэффициент 3,284 — приближенным числом.

3. Составить таблицу значений функции

5 ,64- У'х2+ 1,24

8 ,32+ £/х 2+ 1,24

для х, равного 2,0; 2,4; 2,8; 3,2; 3,6. В ответе сохра­ нить лишь верные цифры. Значения х считать точными,

числа

5,64; 8,32; 1,24— приближенными.

4.

Составить таблицу значений функции

У~х

Уy i T T + J / x - i

для значений х: 1; 1,2; 1,4; 1,6; 1,8; 2,0. Значения х считать точными. Значения квадратного и кубического

 

 

Т а б л и ц а VI

Формула погрешностей основных элементарных функций

Функция

Предельная абсолютная

Предельная относитель­

погрешность

ная погрешность

1 ^ II * 8

2.у = а*

3.у = Ig x

4.у = sinx

5.j/= c o sx

6.y = t g x

7.у —ctg x

f# —arcsine

\(/= arccosA;

9 1 (/ = arctgx

\t/ = arcctg x

10.y = u - v

11.y = u - v - w

12. y t/

13. y = /" й

Ду=

| а 1.1 х |а 1• Л*

Ду= a-v-| In а | • ДЛ

Л

м ,

д у = Т Л*

Л4 = 0 ,4 3 4 2 9 ...

Ду= |

cosx|-A x

Ду==1Sin х | • Дд.

А —

1

. д„

У

COSSX

А

д „ = ^ 4 - • д *

у

siti"

х

 

Л

 

-1..-

. А..

у

V I

- V

х

Ау 1 +1х 2 А

бу =

| а |- б х

 

 

Sy=

| In а |-Дх

 

 

х

y“

м

 

 

s

 

 

|l g x |

 

 

Sj,= |ctgx|-A *

 

8j,= |t g x |.A x

 

б

 

 

2

. д

° у

| sin 2х |

 

 

° и

 

 

2

 

*

 

|sin 2х |

 

Л

J/

— —

 

* ..

2

. Л .

 

 

1/1

 

 

 

 

 

у V

1 —х2

 

б

 

А

-6

Д у = |« |- До + М ' Ди 6 у = б ц + б г,

Д = |м у |-Д и+ |к -а ||-Дг,+ 6у = ба + б® + 4- \uv 1■д®

лм - д „ + м - д ®

д - д “

&у— ~2 6U

у2 ^ ' «

62

корней брать с четырьмя верными десятичными знаками. В ответе сохранить лишь верные цифры.

5. Вычислить у =

ное число). В ответе сохранить верные цифры.

8.

Вычислить у = |/ х ф

х ф j / х ф

х ф ]/^х, если

х — 5,64 (точное

число).

 

 

 

7.

Вычислить

г/ = х4ф |/

х,

если

 

а)

х = 1 ,3 8

(точное число),

число).

 

б)

лг = 1,38

(приближенное

 

В случае (а) вычислить с точностью до 0,00005.

8.

Составить расчетный бланк для вычисления зна­

чений функций:

 

 

 

 

9.Подготовить расчетный бланк для составления

таблицы функции с шагом h (шаг таблицы— разность между двумя соседними значениями аргумента) на отрез­ ке а Вычислить два-три значения функции без строгого учета погрешностей:

Глава УСТРОЙСТВО И УПОТРЕБЛЕНИЕ 3 МАТЕМАТИЧЕСКИХ ТАБЛИЦ

§ 9. ОБЩИЕ СВЕДЕНИЯ О ТАБЛИЦАХ

29°. Математические таблицы. В вычислительной прак­ тике широкое применение находят таблицы значений раз­ личных функций. Таблицы являются важным средством вычислений. Существуют различные сборники таблиц и спра­ вочники, в которых помещены таблицы значений разных функций (чаще всего приближенных значений). Наиболее часто применяются в практике вычислений таблицы таких

функций: а) алгебраических: х2, я3, У х , У х (квадратов,

кубов, корней квадратных, корней кубических), (обрат­

ных величин); б) трансцендентных: sin х, cos х, tg х, ctg х, логарифмической, показательной функции.

Простейшие таблицы устроены так. В первом столбце помещается ряд значений аргумента (х0, лч, хг, ...), во втором

столбце — соответствующие

значения

функции

(у0, уи

у.2, ...) (см. табл. VII). В

качестве конкретного

примера

возьмем отрывок из таблицы квадратов,

помещенной в БС

(у нас табл. VIII).

Разность между двумя соседними значениями аргу­ мента в таблице называется, как известно, ш а г о м т а б ­ л и ц ы:

Лл=**+1— хк.

В большинстве таблиц шаг есть величина постоянная; будем его обозначать просто h. В таблице VIII Л=0,01.

Разность между двумя соседними значениями функции называется табличной разностью: 6.yk= y h+ ly k.

Табличные разности обычно выражаются в единицах низшего разряда. В таблице VIII табличные разности за­ писаны в сотых долях, фактически они равны: 0,12; 0,12; 0,13; 0,12.

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

64

Т а б л и ц а V I I

 

Т а б л и ц а V I I I

X

У

X

Х %

Ду

х 9

Уо

6,10

37,21

12

х\

Уi

 

 

 

Хг

У-2

б,п

37,33

12

Х ц - \

У к - \

6,12

37,45

13

6,13

 

12

А>

Ук

37,58

Хц + 1

Ук +1

6,14

37,70

 

 

 

 

30°. Линейная интерполяция. Таблица функции y= f(x) дает непосредственно значения этой функции только для тех значений аргумента, которые содержатся в таблице. Между тем часто требуется определить значения функции для значений аргумента, не содержащихся в таблице, т. е. для промежуточных значений аргумента. Например, пусть требуется найти х 2 при х —6,134 по таблице VIII. Значения х=6,134 нет в таблице, это значение заключено между соседними табличными значениями 6,13 и 6,14.

Задача об отыскании значения функции для промежу­ точных значений аргумента решается при помощи интер­ поляции. Простейшим видом интерполяции является ли­ нейная интерполяция. Принцип линейной интерполяции разберем на только что приведенном примере. Итак, мы знаем из таблицы значения функции х2при таких значениях аргумента:

*в =

6,130

(«начальное значение

аргумента») и

*х =

6,140

(«конечное

значение

аргумента»):

 

6,1ЗО3 =

37,58

(</„ = 37,58)

 

6,140* ~

37,70

(у, =37,70).

Табличная разность Ay0= l/i—р»=12 (в единицах низшего разряда, т. е. в сотых долях).

Мы видим, что когда аргумент увеличился на 10 единиц (в тысячных долях) по сравнению с начальным значе­ нием 6,130, то функция увеличилась на 12 единиц (в сотых долях). На сколько же увеличится функция, когда аргу­ мент увеличится на 4 единицы? Обозначим это искомое приращение функции Ау. При малых изменениях аргу­ мента можно допустить (и в этом заключается основная

3 № T J7 2

65

идея метода линейкой интерполяции), что приращение функции пропорционально приращению аргумента. Значит,

можно записать пропорцию:

Д^__4

12 10

Отсюда Ау—4,8. Округляя, получим Дг/=5 (в сотых долях). Итак, искомое значение функции равно 37,58+0,05=37,63.

С геометрической точки зрения линейная интерполяция - означает, что на отрезке 1л0; x j дугу кривой, которая является графиком функции, заменяем хордой этой дуги, т. е. истинный криволинейный график функции прибли­ женно заменяем прямолинейным (рис. 11), Для прямоли­ нейного графика как раз и характерно то, что прираще­ ние у пропорционально приращению х.

Проведенные рассуждения можно в общем виде вы­ разить посредством формул. Пусть х0 и хх—два соседних табличных значения аргумента, им соответствуют значе­ ния функции у0 и ух, которые можно найти по таблице.

Требуется

найти значение функции у

при промежуточном

значении

х:

х0 < С х< хх;

хх—-х0 есть шаг таблицы й;

ух—йо есть

табличная

разность,

обозначим ее Ау&.

Обозначим: х —х0 = Дх; у у0~ А у . Достаточно найти Ду (эта величина называется поправкой). Так как мы до­ пускаем, что приращение функции пропорционально приращению аргумента, то

Ду Ах

66

Отсюда находим поправку:

 

 

 

 

 

 

 

(3.1)

Получаем:

 

 

 

 

 

 

 

У= Уо + ^~- А«/о-

 

(3.2)

Это и есть формула линейной интерполяции.

таб­

П р и м е р

1.

Найти хя при х = 1,487, пользуясь

лицей

2 из

БС

(«Квадраты, кубы, квадратные и куби­

ческие корни»). Отрывок из этой таблицы—таблица

IX.

 

Т а бл и ц а

I X

 

 

 

X

 

V

*0=1.40,

у0 = 3,242,

 

 

 

 

 

 

 

хх= 1,490,

у, = 3,308.

 

1 .4 8

3 ,2 4 2

 

/1 = 1 0 (тысячных),

 

 

Д#0 = 66 (тысячных),

 

 

 

66

 

 

3 ,3 0 8

Ах — 7 (тысячных).

 

1 .4 9

 

 

Пользуясь формулой (3.1), найдем поправку:

7

Ay = jq • 66 = 46,2 (тысячных).

Отсюда у = 3,242 + 0,046 = 3,288. Итак, 1,487* = 3,288.

П р и м е р

2. Найти

, пользуясь таблицей обрат­

ных величин

(см. БС, табл.

4).

 

Таблица X

X

1

 

X

 

 

 

2 ,1 2 0

0 ,4 7 1 7

2 ,1 3 0

 

— 22

0 ,8 6

9 5

h = 10, Дх = 3, Ау()= —22.

По формуле (3.1) получим:

дУ = - Ш * ( - 2 2 ) « - 7 .

Следовательно, у = 0,4717 — 0,0007 = 0,4710.

О т в е т .

= 0,4710.

3*

67

П р и м е р

3. Найти

sin35°26\ пользуясь

таблицей

«Натуральные

значения

тригонометрических

функций»

(БС, табл. 9).

 

 

 

X

sin х ЛУо

 

 

35°20'

0 ,5 7 8 3

 

24

35°30'

0 ,5 8 0 7

/г — 10 (минут),

Ах: = 6 (минут), Ау0 — 24 (десятитысячных).

— 0

По формуле (3.1) находим Аг/ = ~ • 24 = 14,4 л* 14 (деся­

титысячных). Стало быть,

sin 35°26' = 0,5783 4-0,0014 =

= 0,5797.

 

 

 

 

 

Для облегчения вычислений при интерполировании

употребляются таблицы пропорциональных

частей (РР).

В них приводятся для каждой

табличной

разности про­

изведения этой табличной разности на числа

 

При помощи

таких таблиц

можно быстро находить А

по формуле (3.1).

Таблицы

РР

помещаются обычно на

отдельном листе-вкладыше.

 

 

 

П р и м е р

4. Пользуясь четырехзначной таблицей ку­

бов и таблицей РР,

найти х 3 при х — 1,483; х= 1,488. Зна­

чение аргумента х заключено между 1,48 и 1,49, поэтому можем опять воспользоваться таблицей IX. Табличная разность равна 66. Таблица XI содержит часть таблицы РР

для разности

66.

Т а б л и ц а X I

По таблице IX:

 

при х = 1,480 л3 = 3,242.

66По таблице РР (таблица XI) по­ правка на последний десятичный

1

6 ,6

21 3,2

31 9 ,8

42 6 ,4

53 3 ,0

63 9 ,6

74 6 ,2

8

52,8

;

9

5 9 ,4

I

знак равна:

на 3 — 19,8; на 8 — 52,8.

Поэтому, округляя поправки, по­ лучи м:

1,483* = 3,242 + 0,020 = 3,262,

1,4883 = 3,242 + 0,053 = 3,295.

£8

Условие допустимости линейной интерполяции. Приме­ нить линейную интерполяцию формально всегда можно, но получим ли мы в результате достаточно точное значение функции? Применение линейной интерполяции может удов­ летворить в том случае, если значение функции, полученное посредством интерполяции, будет иметь ту же точность, что и табличное значение функции. Если таблица (или некоторая ее часть) удовлетворяет этому требованию, то говорят, что линейная интерполяция для данной таблицы (или для соответствующей ее части) допустима.

Методами математического анализа доказывается, что

таблица с постоянным шагом допускает

линейную

интер­

поляцию, если

соседние

табличные разности отличаются

друг'от друга не более чем на 4 единицы

низшего

разряда

табличного значения функции.

 

 

В качестве

примера

приведем таблицы:

 

 

Т а б л и ц а А

 

Т а б л и ц а В

 

 

X

tg*

ь у

 

 

78°00'

4,705

210

 

 

78°30'

4,915

 

 

230

 

 

79°00'

5,145

 

 

251

 

 

79°30'

5,396

 

 

 

Таблица А допускает линейную интерполяцию. Таб­ лица В не допускает линейной интерполяции.

Обычно при составлении таблиц стараются удовлетво­ рить требованию допустимости линейной интерполяции.

§ 10= ОБЗОР НЕКОТОРЫХ ТАБЛИЦ

3!°. О таблицах логарифмической функции. Основанием логарифмов может быть любое положительное число а, отличное от 1. Логарифмическая функция с основанием а записывается так: y —\ognx. Широко применяются десятич­ ные логарифмы, т. е. логарифмы по основанию 10. Во всех сборниках таблиц помещают таблицы десятичных лога­ рифмов.

69

Соседние файлы в папке книги из ГПНТБ