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

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

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

При п — 4 (m = 2) будем иметь:

/4= 0,04 (1,00000+ 1,61607)+ 2 (1,27125 + + 4(1,12750 + 1,43333) - 0,6160756.

Сравнивая

/8 и 14, замечаем,

что совпадают все де­

сятичные знаки до пятого знака

после запятой. Значит,

в найденном

значении /8 пять

верных десятичных зна­

ков. Следовательно,

можно принять

0 , 4 8

J

f(x )d x & 0,61608.

о

 

По формуле (9.7) находим, что погрешность не превы­ шает 0,000002.

§ 26. ПРИБЛИЖЕННОЕ ВЫЧИСЛЕНИЕ ГЕОМЕТРИЧЕСКИХ ВЕЛИЧИН

102°. Вычисление площади плоской фигуры. Пусть требуется вычислить площадь плоской фигуры (рис. 67). Примем некоторую прямую за ось абсцисс. Обозначим а наименьшую из абсцисс точек этой фигуры, b— наиболь­ шую. Возьмем произвольную прямую, перпендикулярную оси Ох, пересекающую фигуру по отрезку M N. Длину этого отрезка обозначим d. Величина d есть функция

абсциссы х. Обозначим ее:

d = D(x). Из

математического

 

анализа

 

известно,

что

 

площадь фигуры S вы­

 

ражается

интегралом

 

ь

 

 

 

 

 

а

 

 

 

 

 

 

Если

у

 

нас

нет

удоб­

 

ного

 

аналитического

 

выражения

функции

 

D(x),

то вычисляем ин­

 

теграл

 

 

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

 

х В частности,

если D(x)

Рис. 67

вообще

 

не

выражена

210

формулой, то произво­

 

дим измерения и при­

 

меняем формулы

меха­

 

нических

квадратур

по

 

данным измерений.

На

 

практике

это

делается

 

так. Выбираем удобное

 

расположение

оси

аб­

 

сцисс,

устанавливаем

 

крайние

абсциссы

 

а

и

1

Ъ. Делим отрезок

[а;

b]

х

на некоторое

число

п

Рис. 68

равных

частей,

через

 

точки деления проводим прямые, перпендикулярные оси Ох, и измеряем длины сечений d0, du d2, . . . , dn (рис. 68). Расстояние между двумя соседними прямыми обозначим h. По сути дела мы получаем таблицу функции D(x).

Применяя теперь для приближенного вычисления ин­ теграла (9.10), например, формулу трапеций, получим:

(9.11)

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

к

0

1

2

3

4

5

6

7

8

d k

0

120

162

174

186

188

175

124

0

При этом h 50 м. Площадь фигуры вычисляем по формуле (9.11), получаем: S = 56 450 мг.

103°. Вычисление объемов тел. Пусть требуется найти объем тела. Примем некоторую прямую за ось Ох. На­ именьшую и наибольшую абсциссы точек тела обозна­ чим а и Ь. Проведем произвольную плоскость, перпен­ дикулярную оси Ох, пересекающую тело по плоской фигуре, площадь которой 5 (рис. 69). Величина S есть функция абсциссы х. Обозначим эту функцию: S = P(x).

5

211

Из математического анализа известно, что объем тела выражается формулой

V = \p { x )d x .

(9.12)

Если у нас нет удобного аналитического выражения функции Р(х), то интеграл (9.12) вычисляем прибли­ женно. В частности, если Р (х) вообще не выражена формулой, то производим измерения и находим прибли­

женное значение V на основании измерений.

 

 

 

Практически это делается так. Во­

 

ображаем, что дана некоторая ось Ох,

 

мысленно проводим плоскости, пер­

 

пендикулярные этой оси, так, чтобы

 

расстояние

между

двумя

соседними

 

плоскостями было одно и то же; обоз­

 

начим его h. Измеряем площади

S0,

 

Si, . . . , S „

сечений тела этими плоско­

 

стями. Применяя теперь формулу тра­

 

пеций, получаем:

 

 

 

 

V = A

+ «■»■- f

S , +

S 2 + . . .

Sn_j

 

 

 

 

 

 

 

(9.

3)

 

П р и м е р .

Вычислить емкость

со-

Рис. 70

суда

имеющего форму тела

вращения,

212

осевое сечение которого изображено на рисунке 70. Производим измерения диаметров dk круговых сечений плоскостями, перпендикулярными оси. Расстояние между двумя соседними плоскостями — 10 см, диаметры сечений даны в таблице (в дм). По этим значениям диаметров вычислены площади сечений.

k

0

1

2 ,

3

4

5

6

7

8

dk

2,75

3,16

4,38

5,23

6,04

4,86

2,94

1,63

1,87

Sk

5,94

7,84

15,07

21,48

28,65

18,55

6,79

2,08

2,75

Применяя формулу (9.13), получим;

V — 105 дм3.

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

1. Вычислить приведенные ниже интегралы двумя способами: а) по одной из формул приближенного вычис­ ления интегралов; б) по формуле Ньютона—Лейбница. Сравнить результаты.2

2

2

dx

9

^ x* dx\

f‘

S

1

 

 

 

0

 

1

8

Л

 

1

3

 

CY x dx<

\ sin xdx\

I

и

J

 

0

0

 

0

2. Вычислить приведенные ниже интегралы но фор­ муле Симпсона, оценить погрешность результата. Число 2п частичных отрезков указано в скобках.

3

6

а) Г \^4~+~xsdx( 2 н =

16); б) J у = ? ( 2 н = 1 6 ) ;

о

0

 

я

10

6

в) ( Д (2я = 8);

г) J | / t g * d * ( 2r t = 12);

4

0

213

7 л

д ) | ^ ( 2« = 6); е) | ^ Л с ( 2я = 12).

1

2

3. Вычислить приближенное значение числа л, при­

меняя

формулу

Симпсона

к

вычислению

интеграла

1 т т ? < 2" = ,6>-

 

 

 

 

 

о

 

 

 

 

 

 

4.

Вычислить

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

значение In 10,

приме­

няя

формулу

Симпсона

к

вычислению

интеграла

! т$ 1 (2я=|6>-

 

 

 

 

 

о

 

 

 

 

 

 

5.

Вычислить

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

значение

lg 3,

приме­

няя

формулу

Симпсона

к

вычислению

интеграла

з

 

 

 

 

 

 

(2rt= 16)-

У к а з а н и е , lg N = М - In N, где М = lge = 0,43429...

6. Вычислить по формуле Симпсона интеграл

ь

\j f{x)dx, пользуясь следующей таблицей значений функ-

а

ций f(x).

3,00

3,25

3,50

3,75

4,00

4,25

4,50

10,392

11,969

13,668

15,490

17,436

19,505

21,698

 

 

 

 

 

П родолжение

4,75

5,00

 

5,25

5,50

5,75

6,00

24,016

26,458

29,024

31,715

34,530

37,470

Оценить погрешность по правилу удвоения.

214

7. Вычислить

площадь

 

зеркала водоема по дан­

 

ным чертежа (рис. 71),

 

применяя формулу трапе­

 

ций (длины указаны в мет­

 

рах).

 

 

 

8. Для вычисления ра­

 

боты пара в цилиндре па­

 

ровой машины вычисляют

 

площадь

индикаторной

 

диаграммы,

представля­

Рис. 71

ющей собой

графическое

изображение

зависимости

цилиндре и ходом поршня.

между давлением

пара в

На рисунке 72

изображена индикаторная диаграмма

паровой машины.

Ординаты

точек линий ABD и ED,

соответствующие

абсциссам

хк, обозначены соответст­

венно yh й yk (k — 0,

1, 2, . . . ,

10).

 

 

 

 

 

х к

Ч

*1

Ч

*3

X*

Ч

Ч

Ч

Ч

х я

x to

 

 

 

 

 

 

 

 

 

 

 

Ч

6 , 2

2 , 4

1 . 2

0 , 9

0 , 8

0 , 7

0 , 7

0 , 8

0 , 9

1 , 3

2 , 9

 

 

 

 

 

 

 

 

 

 

 

Ук 5 4 , 6

4 8 , 7

31 ,4

2 2 , 6

2 0 , 1

1 8 , 9

1 5 , 9

1 2 , 6

1 0 , 4

6 , 3

3 , 8

Вычислить с помощью формулы Симпсона площадь диа­ граммы, если лг0= 0, х10—90 мм.

215

V

9.

Скорость

v

движения

автомобиля (в км/ч) в зави­

симости от времени выражается графиком

(рис. 73).

Если v =

/(/),

то

путь,

пройденный

телом,

выражается

 

 

 

 

т

 

 

 

 

 

 

 

 

 

интегралом

s =

^ / (/) dt.

 

Пользуясь

графиком

скорости,

найти

 

 

 

*0

путь

за

время от 0 ч до

2 ч.

пройденный

10.

Длина дуги

 

 

 

 

 

k

 

 

I гиперболы У = — от точки с абсцис­

сой а

до

точки

с

абсциссой

&(0 < а < Ь )

выражается

интегралом

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

V~P

 

 

 

 

 

 

 

 

 

 

-f

dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Применяя формулу Симпсона, полагая k = 10,

вычислить

длину

дуги

гиперболы,

 

заключенной между

точками

А (2; 5) и В (10; 1).

 

 

 

с

полуосями а и Ъ выражается

1 1.

 

Длина

эллипса

формулой:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f Л

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

S = С V

1 — е2sin <р dq>, е2= а ~2Ь .

 

Вычислить с помощью формулы Симпсона длину эллип­

са, если

а — 20 м, Ь= 10 м (2п=16).

 

 

21G

12. В некоторых кол­

 

 

 

хозах и совхозах для вы­

 

 

 

числения

объема

стога

с

 

 

 

широким

 

основанием

и

 

 

 

закругленной

вершиной

 

 

 

применяется

приближен­

 

 

 

ная формула:

 

 

 

 

 

 

V=z (25~~8§) С*м3’

 

 

 

 

где Я —перекид,

С—дли-

 

Рис. 74

на окружности основания.

 

 

 

Проведено вычисление объема

по этой формуле при дан­

ных: Я =

9,7 м\ С = 12,5 м. Для

контроля

проведено вы­

числение

по Формуле трапеций:

измерили

высоту стога

Я = 4 м и

длины

окружностей

горизонтальных сечений

на высотах

1,

2,

3 м.

Длины

окружностей оказались

равными

соответственно

12,2;

10,9; 8,3 м. Считая второй

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

13. Объем тела, вырезанного из цилиндра радиуса R цилиндром радиуса г (г < R), если оси цилиндров пере­ секаются под прямым углом, выражается интегралом

Т

V — 8 ^ У (R2х%) 2х2) dx.

о

Вычислить объем такого тела, пользуясь формулой тра­

пеций

(п = 20), полагая R = 50 см, г = 20 см (рис. 74).

14.

Вычислить интеграл

 

ю

 

^ (х — 5)2(10—х) dx,

 

о

пользуясь малой формулой Симпсона. (Заметим, что результат будет точным.)

15. Длина s полуволны синусоиды выражается интег­

ралом

31

s = J Y 1 + cos2х dx.

о

Вычислить s, пользуясь формулой Симпсона (2п = 32).

217

Г л а в а ДОКАЗАТЕЛЬСТВА НЕКОТОРЫХ ТЕОРЕМ

10И ВЫВОДЫ ФОРМУЛ

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

§ 27. ОЦЕНКА ПОГРЕШНОСТИ МЕТОДАМИ ДИФФЕРЕНЦИАЛЬНОГО ИСЧИСЛЕНИЯ

104°. Погрешности значений функции. В главе 2 было показано, что истинную погрешность значения функции при малой погрешности аргумента можно принять рав­ ной дифференциалу функции, т. е. если y = f(x), то Ау « dy, или Ау ж /' (х) dx, или Ау « f (х) Ах (Ах и Ау соответственно погрешности аргумента и функции). На основании этого можно доказать справедливость формул таблицы VI, пользуясь формулами производных, извест­ ными из математического анализа (см. приложение).

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

Ay « (uv)' dx ж (u’v + uv') dx = vur dx + uv' dx.

Ho u 'd x =^du, v’ dx = dv. Итак, Au m vdu-}-и dv. Ho du, dv

можно заменить через Aи, Av. Поэтому

| Amj < j v j • j Am| + 1и j • | Ay |.

Отсюда следует, что можно принять:

А« = М-А«-Н“|-А*-

Это и есть формула 10 таблицы VI.

Подобным образом в случае произведения трех фуню ций получаем формулу 11 таблицы VI.

218

Погрешность частного. Пусть у — ^ .

Отсюда

[ A# | < -l£jl !..4“. I.+ li± L M ..

Стало быть, можно принять

д

_ М ’

+ I И 1■

.

У

v2

(формула 12 таблицы VII).

Погрешность показательной и логарифмической функ­

ции. Пусть

у = ах

или z — ]gx.

Имеем (ах)' = ах Лпа,

(loga x)' = y loga e.

Поэтому

 

« ах In а ■Ах,

A

л; — • lg е- Ах.

Отсюда заключаем

(обозначаем М =

lg е):

Ау ах | In а I ■Ах,

Аг = ~ • Ах

 

(формулы 2 и 3 таблицы VI).

доказываются остальные

Подобными рассуждениями

формулы таблицы VI.

 

 

 

§ 28. О СХОДИМОСТИ МЕТОДОВ ПРИБЛИЖЕННОГО РЕШЕНИЯ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ

105°. Постановка задачи. В главе 4 были рассмот­ рены методы приближенного решения (уточнения кор­ ней) уравнений с одним неизвестным. Были рассмотрены, в частности, метод хорд, метод касательных и метод

итераций.

Сущность применения

этих методов состоит

в следующем.

 

Дано

уравнение

 

 

/(* )= 0.

(10.1)

Будем предполагать, что уравнение имеет корень 1, ко­ торый отделен на отрезке [а; £>]. В результате примене-

219

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