книги из ГПНТБ / Пулькин, С. П. Вычислительная математика пособие для учащихся 9-10 классов по факультативному курсу
.pdfПри п — 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 = 4а С 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