- •Введение
- •Раздел I. МЕТОДЫ ПРИБЛИЖЕННЫХ ВЫЧИСЛЕНИЙ
- •1. Приближенное решение уравнения
- •2. Метод последовательных приближений решения уравнения
- •3. Метод Ньютона (метод касательных) решения уравнения
- •4. Метод секущих (метод хорд) решения уравнения
- •8. Приближенные вычисления значений функций с помощью рядов
- •9. Приближенные вычисления пределов с помощью рядов
- •10. Задачи для самостоятельного решения
- •Раздел II. ПРИБЛИЖЕННОЕ ВЫЧИСЛЕНИЕ ОПРЕДЕЛЕННЫХ ИНТЕГРАЛОВ
- •1. Квадратурные формулы
- •3. Формула трапеций
- •4. Метод парабол (метод Симпсона)
- •5. Вычисление определенных интегралов с помощью рядов
- •6. Задачи для самостоятельного решения
- •Раздел III. ЧИСЛЕННОЕ ИНТЕГРИРОВАНИЕ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •1. Метод Пикара последовательных приближений
- •2. Метод Эйлера
- •3. Улучшенный метод Эйлера
- •4. Метод Эйлера–Коши
- •5. Метод Рунге–Кутта
- •6. Метод Адамса
- •8. Метод конечных разностей решения краевых задач для обыкновенных дифференциальных уравнений
- •10. Метод наименьших квадратов
- •11. Задачи для самостоятельного решения
- •Раздел IV. УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ ГИПЕРБОЛИЧЕСКОГО ТИПА
- •1. Виды уравнений математической физики
- •2. Вывод уравнения колебания струны
- •7. Задачи для самостоятельного решения
- •Раздел V. УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ ПАРАБОЛИЧЕСКОГО ТИПА
- •4. Распространение тепла в неограниченном стержне
- •6. Задачи для самостоятельного решения
- •Раздел VI. УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ ЭЛЛИПТИЧЕСКОГО ТИПА
- •1. Уравнения эллиптического типа и краевые задачи для них
- •3. Уравнение Лапласа в цилиндрических координатах
- •5. Примеры решения задач
- •6. Задачи для самостоятельного решения
- •Раздел VII. ЧИСЛЕННОЕ ИНТЕГРИРОВАНИЕ УРАВНЕНИЙ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ
- •1. Метод сеток численного решения дифференциальных уравнений с частными производными
- •3. Метод сеток для уравнения гиперболического типа
- •5. Метод прогонки для уравнения теплопроводности
- •6. Решение уравнения движения грунта
- •7. Задачи для самостоятельного решения
- •Раздел VIII. ВАРИАЦИОННЫЕ МЕТОДЫ
- •1. Некоторые сведения из функционального анализа
- •2. Теоретические основы метода Ритца
- •5. Задачи для самостоятельного решения
- •Раздел IX. КОРРЕЛЯЦИЯ
- •1. Понятие корреляции
- •4. Корреляция
- •5. Задачи для самостоятельного решения
- •Библиографический список
Вычислим теперь исходный интеграл по формуле Ньютона– Лейбница:
sin
0
3x
2
|
1 |
|
|
|
2 |
|
|
|
d x |
|
|
cos |
|
|
|
|||||
|
2 |
|
|
|
3 |
|
3x |
|
1 |
|
|
|
|
|
||||
|
|
|
|||
|
|
|
x |
|
|
|
|
||||
2 |
|
2 |
|
|
0 |
|
|
|
|
СибАДИ |
|||||||||||||||||||||
|
|
2 |
|
3 |
|
|
|
|
|
2 |
|
3 |
0 |
|
0 |
|
|
|
2 |
|
|
|
|
|
cos |
|
|
|
|
|
|
|
cos |
|
|
|
|
|
|
|
|
|
2,237463 . |
|
|
3 |
|
2 2 |
|
|
|
3 |
|
|
2 |
|
2 |
|
2 3 |
|
Результаты совпадают с точностью до десятитысячных. Точность метода Симпсона (парабол) выше точности метода
прямоугольн ков трапеций для заданного n (это видно из оценки абсолютной погрешности), так что его использование предпочтительнее.
Следует помн ть о влиянии вычислительной погрешности на результат при больш х n, что может отдалить приближенное значение от точного.
5. Вычисление определенных интегралов с помощью рядов
При вычислении определенных интегралов иногда удобно разложить функции в степенные ряды.
Рассмотрим примеры приближенного вычисления определенных интегралов.
Примеры.
|
|
|
|
|
|
|
|
|
0,3 |
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
1. Вычислим интеграл |
|
|
|
|
|
с точностью до 0,0001. |
||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||
3 1 x2 |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
0,2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Решение. Разложим подынтегральную функцию в степенной |
||||||||||||||||||||||||||||||
ряд, пользуясь биномиальным разложением: |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
1 |
|
4 |
|
|
|
|
|
1 |
|
4 |
|
|
7 |
|
|
|
|
||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
1 |
|
|
|
3 |
|
3 |
|
2 |
|
|
|
3 |
|
3 |
|
|
3 |
|
|
3 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
1 x |
3 |
1 |
|
|
x |
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|||
3 |
|
2 |
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
1 n 1 4 7 3n 2 xn (разложение верно при |
1 x 1). |
|||||||||||||||||||||||||||||
|
|
|
3n n! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
41
Теперь подставим x2 вместо x , получим
1 x2 |
1 |
|
|
1 n 1 4 |
7 |
3n 2 x2n |
|
||
3 |
1 |
, при 1 x 1. |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
n 1 |
|
|
|
3n n! |
|
СибАДИ |
|||||||||||||||||||||||||||||||
|
|
Воспользуемся возможностью почленного интегрирования сте- |
|||||||||||||||||||||||||||||
пенного ряда: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
0,3 |
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,2 3 1 x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
0,3 |
|
|
|
1 |
|
|
2 |
|
2 |
|
|
4 |
|
|
14 |
|
|
6 |
|
|
n 1 4 7 3n 2 |
|
2n |
|
|
||||||
|
1 |
|
|
|
x |
|
|
|
|
x |
|
|
|
|
|
|
|
x |
|
|
1 |
x |
|
dx |
|
||||||
0,2 |
|
|
|
3 |
|
|
|
|
9 |
|
|
|
|
|
27 |
|
|
|
|
|
3n n! |
|
|
|
|
||||||
|
|
1 |
|
|
|
3 |
|
|
2 |
|
|
|
5 |
|
|
|
2 |
|
|
|
7 |
|
|
|
0,3 |
|
|
|
|
||
x |
|
|
x |
|
|
|
|
x |
|
|
|
|
|
|
x |
|
|
|
0,1 0,00211 0,00009 0,09798 |
|
|||||||||||
|
|
9 |
|
|
|
|
|
|
45 |
|
|
|
|
81 |
|
|
|
|
|
|
0,2 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,0980.
После почленного интегрирования получили знакочередующийся ряд, удовлетворяющий условиям теоремы Лейбница. Восполь-
зуемся оценкой остатка сходящегося |
знакочередующегося ряда |
|||
|
Rn |
|
an 1 . Поскольку четвертый член |
ряда по модулю оказался |
|
|
меньше заданной точности 0,0001, то для вычислений достаточно взять сумму первых трех слагаемых.
Итак,
0,3 |
|
dx |
0,0980. |
|
|
|
|||
3 |
1 x2 |
|||
0,2 |
|
Отметим, что если после почленного интегрирования получится знакоположительный ряд, то для определения необходимого для вычислений числа слагаемых проводят оценку остатка ряда (обычно
42
оценивают с помощью геометрической прогрессии).
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2. Вычислим интеграл e x3 dx с точностью до 0,0001. |
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Решение. Разложим подынтегральную функцию в ряд, исполь- |
|||||||||||||||||||||||||||||||||||||||
зуя разложение в ряд Маклорена функции ex , |
заменяя x на x3 . По- |
||||||||||||||||||||||||||||||||||||||
лучим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
СибАДИ |
|||||||||||||||||||||||||||||||||||||||
1 |
x |
3 |
|
|
|
|
1 |
|
|
|
|
x3 |
|
x6 |
|
x9 |
|
x12 |
|
|
x15 |
|
|
x18 |
|
|
|
|
|||||||||||
e |
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
1 |
1! |
2! |
3! |
4! |
5! |
|
6! |
dx |
||||||||||||||||||||||||||||
0 |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
4 |
|
|
|
|
|
7 |
|
|
|
10 |
|
13 |
|
|
16 |
|
|
19 |
|
|
1 |
|
|||||||||||
x |
x |
|
|
|
|
x |
|
|
|
x |
|
|
|
|
x |
|
|
|
x |
|
|
x |
|
|
|
||||||||||||||
|
|
|
|
|
4 1! |
|
7 2! |
10 3! |
13 |
4! |
16 5! |
19 6! |
|
|
|
|
|||||||||||||||||||||||
|
1 |
|
|
1 |
|
|
|
1 |
|
|
1 |
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
0 |
|||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
4 |
|
14 |
|
|
60 |
312 |
1920 |
|
|
13680 |
|
|
|
|
|
|
|
|
|||||||||||||||||||
1 0,25 0,07142 |
0,01666 |
0,00320 |
0,00052 |
0,00007 |
|||||||||||||||||||||||||||||||||||
0,80744 0,8074 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Поскольку мы использовали знакочередующийся ряд, то для достижения точности вычислений необходимо было найти слагаемое меньшее точности по модулю, отбросить его и остальные члены ряда.
Так как 1 0,00007 0,0001 , то можно отбросить это слагаемое
13680
все остальные. Возникающая при этом погрешность будет меньше, чем 0,00007, то есть меньше 0,0001. Все вычисления мы проводили с пятью знаками после запятой, результат округлили до четырех знаков после запятой.
Итак,
1 e x3 dx 0,8074
.
0
43
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 sin x |
|
|
|
|
|
|
|
||||||
3. Вычислим |
|
|
теперь |
|
|
интеграл |
|
|
|
|
|
|
|
|
dx |
|
с |
точностью |
до |
|||||||||||||||||
|
|
|
|
|
|
|
|
x |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|||
0,5 10 3 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Решение. |
Заметим, |
что неопределенный интеграл |
sin x |
dx |
не |
|||||||||||||||||||||||||||||||
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
СибАДИ |
||||||||||||||||||||||||||||||||||||
вычисляется в конечном виде («неберущийся» интеграл). |
|
|||||||||||||||||||||||||||||||||||
Разлож м sin x в ряд и поделим почленно на x . Получим ряд |
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
sin x |
1 |
x2 |
|
|
x4 |
|
x6 |
|
, |
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
3! |
|
|
|
5! |
|
|
7! |
|
|
|
|
|
|
|
|
|
|
||||
сходящ йся при лю ом значении x . |
нтегрируя, получим |
|
||||||||||||||||||||||||||||||||||
2 sin x |
|
|
|
|
|
|
|
|
|
x3 |
|
|
|
x5 |
|
|
|
|
x7 |
|
|
2 |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
dx x |
|
|
|
|
|
|
5 5! |
|
|
|
|
|
|||||||||||||||||||||
0 |
|
x |
|
|
|
|
|
|
|
|
3 |
3! |
7 |
7! |
|
|
0 |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
3 |
|
|
|
|
|
|
|
5 |
|
|
|
7 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
2 |
|
|
|
|
|
600 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
18 |
|
|
35280 |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
1,5708 0,2153 0,0159 0,0007 1,3707 |
1,371. |
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Первый отброшенный член |
|
много меньше, |
чем 0,5 10 3 . |
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 9! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Итак, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 sin x |
dx 1,371 . |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6. Задачи для самостоятельного решения
1. Написать формулу прямоугольников приближенного вычисления определенного интеграла, соответствующего рис. 8.
44
у
|
y2 |
y1 |
y3 |
Си0 бАДИ значением интеграла.
x0 |
x1 |
x2 |
|
x3 |
x |
|
|
|
|
Рис. 8 |
|
2. Выч сл ть |
нтеграл |
0,2 sin x |
dx : |
|
|
|
x |
|
|||
|
|
0 |
|
|
|
а) по формуле прямоугольников; |
|
||||
б) по формуле трапеций; |
|
|
|||
в) по формуле Симпсона. |
|
|
|||
|
|
2 |
x sin x 2 dx: |
|
|
3. Вычислить интеграл |
|
|
|||
|
|
0 |
|
|
|
а) по формуле прямоугольников; |
|
||||
б) по формуле трапеций; |
|
|
|||
в) по формуле Симпсона. |
|
|
|||
4. Применяя формулу прямоугольников (n=12), приближенно вычис- |
|||||
|
2 |
|
|
|
|
лить интеграл |
x sin x dx |
результат вычислений сравнить с точным |
5. С помощью формулы трапеций вычислить интегралы и оценить их погрешности:
1 |
1 |
|
|
а) |
dx (n=8); |
||
1 x |
|||
0 |
|
45
|
1 |
1 |
|
|
|
|
||
б) |
|
|
dx |
(n=12); |
||||
|
|
|
||||||
1 x3 |
||||||||
|
0 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
2 |
|
1 |
1 |
sin 2 x dx (n=6). |
|||
в) |
|
|
||||||
4 |
||||||||
|
0 |
|
|
|
|
|||
СибАДИ |
||||||||
6. помощью формулы Симпсона вычислить интегралы: |
||||||||
|
9 |
|
|
|
|
|
|
|
а) |
|
|
x dx (n=4); |
|||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 sin x |
|
|
|
|
|||
б) |
|
|
x |
dx |
(n=10); |
|||
|
0 |
|
|
|
|
|
|
|
|
|
|
3 cos x dx (n=6); |
|||||
в) |
|
|
||||||
|
0 |
|
|
|
|
|
|
|
|
1 |
|
x |
|
|
|
||
г) |
|
|
ln 1 x |
dx (n=6). |
||||
|
0 |
|
|
|
|
|
|
7. Вычислить определенные интегралы с точностью , раскладывая подынтегральную функцию в ряд:
|
0,2 sin x |
|
|
|
|
, 0,0001 ; |
||||
а) |
|
|
|
|
dx |
|||||
|
x |
|||||||||
|
0 |
|
|
|
|
|
|
|
||
) |
0,5 |
1 cos x |
dx , 0,0001 ; |
|||||||
|
|
x2 |
|
|
|
|
||||
|
0 |
|
|
|
|
|
|
|
||
|
0,1 ex 1 |
dx, 0,001 ; |
||||||||
в) |
|
|
|
|
|
|
||||
|
x |
|
|
|||||||
|
0 |
|
|
|
|
|
|
|
||
|
0,1 ln 1 x |
dx , 0,001 ; |
||||||||
г) |
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
||||
|
0 |
|
|
|
|
|
|
|
||
|
0,5 |
|
x ln 1 x2 dx, 0,001 ; |
|||||||
д) |
|
|
||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
е) |
e x 2 dx , |
0,001 . |
||||||||
|
0 |
|
|
|
|
|
|
|
|
|
46
Ответы:
4.6,2832 ;
5.а) 069315;
б) 0,83566; в) 1, 4675.
СибАДИ6. а) 17,333; б) 5,4024;
в) 1,37039; г) 0,2288.
7. а) 0,1996;
б) 0,2483; в) 0,102; г) 0,098; д) 0,015; е) 0,747.
Вопросы задан я для самопроверк [1, 2, 3, 4, 5, 6, 9, 10]
1. Как составляется интегральная сумма?
2. Дайте определение определенного интеграла.
3. В чем заключается метод Эйлера приближенного вычисления определенного интеграла?
4. В чем заключается метод парабол приближенного вычисления определенного интеграла?
5. В чем заключается метод трапеций приближенного вычисления определенного интеграла?
6. В чем заключается метод приближенного вычисления определенных интегралов с помощью рядов?
7. Какова точность различных методов приближенного вычисления определенного интеграла?
47